Abstract
This paper presents the nonlinear systems of Volterra-type fractional integro-differential equation solutions through a Chebyshev pseudospectral method. The proposed method is based on the Caputo fractional derivative. The results that we get show the accuracy and reliability of the present method. Different nonlinear systems have been solved; the solutions that we get are compared with other methods and the exact solution. Also, from the presented figures, it is easy to conclude that the CPM error converges quickly as compared to other methods. Comparing the exact solution and other techniques reveals that the Chebyshev pseudospectral method has a higher degree of accuracy and converges quickly towards the exact solution. Moreover, it is easy to implement the suggested method for solving fractional-order linear and nonlinear physical problems related to science and engineering.
1. Introduction
Fractional calculus has a long history as classical calculus. The concept of fractional calculus arouse when Leibnitz used a proper representation for the th derivative in his publications. L’hopital raises a question on the particular notation on what happens if “” is a noninteger. It was the beginning of fractional calculus [1]. Recently, mathematicians focused on fractional calculus due to its numerous applications in every field of science: viscoelastic materials [2], economics [3], continuum and statistical mechanics [4], dynamics of interfaces between soft nanoparticles and rough substrates [5], solid mechanics [6], and much more [7–14].
Mathematical formulations solve many problems of nature with the help of converting the physical phenomena to the equation form. Differential equations (DEs) are among those that play the main role in modeling various phenomena. However, some problems are complex and cannot be handled with the help of a differential equation. In this regard, the researchers utilized fractional differential equations (FDEs) that model the phenomenon more accurately than differential equations having order integers. Nowadays, FDEs got the importance of real-world modeling problems: such as electrode-electrolyte polarization [15], electrochemistry of corrosion [16], circuit systems [17], optics and signal processing [18], heat conduction [19], diffusion wave [20], control theory of dynamical systems [21], fluid flow [22], probability and statistics [23, 24], and so on (see [25–28]).
The role of fractional integral and integrodifferential equations is found in every field of engineering and science. When a physical phenomenon is modeled under the differential equation, it finally gives a differential equation, an integral equation, or an integrodifferential equation. Some applications of these types of equations are nanohydrodynamics [29], glass-forming process [30], wind ripple in the desert [31], and drop-wise condensation [32]. The analytical solution of integral and integrodifferential equations does not exist in most cases. Even if it exists in certain cases, it is hard to find. Different numerical methods have been developed for finding an approximate solution of integral and integro-differential equations. The most common among these methods are the Chebyshev polynomials [33], Haar wavelet [34], triangular function method [35], collocation method [36], Legendre wavelet operational method [37], Taylor series expansion method [38], homotopy perturbation method [39], reproducing kernel Hilbert space method [40], Adomian decomposition method [41], Euler wavelet method [42], variational iteration method [43], spectral collocation method [44], least square method [45], homotopy analysis method [46], and differential transform method [47].
We apply Chebyshev pseudospectral method (CPM) to solve nonlinear Volterra integro-differential equation systems in the present work. CPM is a powerful technique for solving linear and nonlinear problems. The obtained results show the higher convergence rate of the present technique. The solution that we get shows that CPM has good agreement with the exact solution. Error analysis reveals the efficiency of the proposed technique that CPM has greater accuracy than other methods.
2. Definitions and Preliminary Concept
This unit shows the preliminary concept and some essential definitions taken from fractional calculus and used in our present research work.
2.1. Definition
The definition for fractional derivative by Caputo of order is showed by the following mathematical expression [48]: for , , , .
2.2. Definition
The fractional derivatives by Jin-Hunan He are described as [48] where does not approach zero.
2.3. Definition
Xiao-Jun explains derivatives having fractional order as [48] where
2.4. Definition
The integral operator by Riemann-Liouville for order is [48]
The Caputo derivative operator and Riemann-Liouville integral operator have the following properties
3. Chebyshev Pseudospectral Method (CPM)
The Chebyshev polynomials are defined in the interval and can be described by the following recurrence formula: where
To apply the Chebyshev polynomials in the interval, we define the Chebyshev shifted polynomials which are defined in the manner of Chebyshev polynomials by relation
And the recurrence formula is as follows: where
A function , in terms of Chebyshev shifted polynomials described as
The Chebyshev shifted polynomials first () terms are considered as
For finding the system of equations, we have whereas
I solved the resultant system using maple software, which provide CPM solution for the given problem.
4. Numerical Representation
4.1. Problem
Consider the nonlinear FIDE system having B.Cs , having , as the exact solution at .
The exact solution and numerical results obtained by means of the proposed method are shown in Table 1. The absolute error comparison of our method and those obtained from OTM is given in Table 2. The behavior of the exact solution and approximate solution (our method) of this example when is presented in Figures 1 and 2 whereas the error comparison of CPM and OTM can be observed in Figures 3 and 4. The graphical representation for different fractional order of is seen in Figures 5 and 6 which confirm that the solution converge to the exact solution as the value of converges from the fractional order to the integer order.
4.2. Problem
Consider the FIDE system with B.Cs , with the exact solution at , where
In Table 3, we give the numerical values of the exact solution and CPM solution for . The absolute errors obtained by the present method are compared with SCM in Table 4. We compare the actual and estimated solution in Figures 7 and 8 which tells us that both the solutions are quite close to each other. Also, Figures 9 and 10 display the error comparison of CPM and SCM which verify that our method is in good agreement with the exact solution.
4.3. Problem
Consider the nonlinear FIDE system having B.Cs , , , and having as the exact solution.
To solve this example, we implement the method suggested in Section 4 for with . The exact solution and estimated solution by CPM are presented in Table 5. The absolute error of our method and those obtained from OTM are given in Table 6. In Figures 11 and 12, it is clear that the numerical solution of the proposed method is in good contact with the exact solution. In order to illustrate the effectiveness of CPM, the error comparison with OTM is shown in Figures 13 and 14. Also, in Figures 15 and 16, we can obtain that as the estimated solution approach to the exact solutions.
5. Conclusion
In this work, we implemented the Chebyshev pseudospectral method for solving nonlinear fractional integral and integro-differential equation systems. The proposed technique reduces this type of systems to the solution of the system of linear and nonlinear algebraic equations. Special attention is given to study the convergence of the proposed method. The results that we get by implementing the suggested technique are in excellent agreement with the exact solution and show more accuracy than the solution obtained using other methods. Also, from the presented figures, it is easy to conclude that the CPM error converges quickly as compared to other methods. The computation work in this article is done using Maple.
Data Availability
The numerical data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no competing interests.
Authors’ Contributions
All authors jointly worked on the results, and they read and approved the final manuscript.
Acknowledgments
This study is supported by the Central Department of Mathematics, Tribhuvan University, Kirtipur, Kathmandu, Nepal.