Abstract
We investigate the numerical solution of linear fractional integro-differential equations by least squares method with aid of shifted Chebyshev polynomial. Some numerical examples are presented to illustrate the theoretical results.
1. Introduction
Many problems can be modeled by fractional Integro-differential equations from various sciences and engineering applications. Furthermore most problems cannot be solved analytically, and hence finding good approximate solutions, using numerical methods, will be very helpful.
Recently, several numerical methods to solve fractional differential equations (FDEs) and fractional Integro-differential equations (FIDEs) have been given. The authors in [1, 2] applied collocation method for solving the following: nonlinear fractional Langevin equation involving two fractional orders in different intervals and fractional Fredholm Integro-differential equations. Chebyshev polynomials method is introduced in [3–5] for solving multiterm fractional orders differential equations and nonlinear Volterra and Fredholm Integro-differential equations of fractional order. The authors in [6] applied variational iteration method for solving fractional Integro-differential equations with the nonlocal boundary conditions. Adomian decomposition method is introduced in [7, 8] for solving fractional diffusion equation and fractional Integro-differential equations. References [9, 10] used homotopy perturbation method for solving nonlinear Fredholm Integro-differential equations of fractional order and system of linear Fredholm fractional Integro-differential equations. Taylor series method is introduced in [11] for solving linear integrofractional differential equations of Volterra type. The authors in [12, 13] give an application of nonlinear fractional differential equations and their approximations and existence and uniqueness theorem for fractional differential equations with integral boundary conditions.
In this paper least squares method with aid of shifted Chebyshev polynomial is applied to solving fractional Integro-differential equations. Least squares method has been studied in [14–18].
In this paper, we are concerned with the numerical solution of the following linear fractional Integro-differential equation: with the following supplementary conditions: where indicates the th Caputo fractional derivative of ; , are given functions, and are real variables varying in the interval , and is the unknown function to be determined.
2. Basic Definitions of Fractional Derivatives
In this section some basic definitions and properties of fractional calculus theory which are necessary for the formulation of the problem are given.
Definition 1. A real function , , is said to be in the space , , if there exists a real number such that , where .
Definition 2. A function , is said to be in the space , , if .
Definition 3. The left sided Riemann-Liouville fractional integral operator of order of a function , , is defined as [19]
Definition 4. Let , . Then the Caputo fractional derivative of is defined as [20–22] Hence, we have the following properties: where denoted the smallest integer greater than or equal to and .
3. Solution of Linear Fractional Integro-Differential Equation
In this section the least squares method with aid of shifted Chebyshev polynomial is applied to study the numerical solution of the fractional Integro-differential (1).
This method is based on approximating the unknown function as where is shifted Chebyshev polynomial of the first kind which is defined in terms of the Chebyshev polynomial by the following relation [23]: and the following recurrence formulae: with initial conditions , , are constants.
Substituting (7) into (1) we obtain Hence the residual equation is defined as
Let where is the positive weight function defined on the interval . In this work we take for simplicity. Thus
So, finding the values of , , which minimize is equivalent to finding the best approximation for the solution of the fractional Integro-differential equation (1).
The minimum value of is obtained by setting Applying (15) to (14) we obtain
By evaluating the above equation for we can obtain a system of linear equations with unknown coefficients ’s. This system can be formed by using matrices form as follows: where
By solving the above system we obtain the values of the unknown coefficients and the approximate solution of (1).
4. Numerical Examples
In this section, some numerical examples of linear fractional Integro-differential equations are presented to illustrate the above results. All results are obtained by using Maple 15.
Example 1. Consider the following fractional Integro-differential equation:
subject to with the exact solution .
Applying the least squares method with aid of shifted Chebyshev polynomial of the first kind , at , to the fractional Integro-differential equation (19) we obtain a system of (6) linear equations with (6) unknown coefficients , . This system can be transformed into a matrix equation and by solving this matrix equation we obtain the inverse which is given in Figure 1 and we obtain the values of the coefficients. Substituting the values of the coefficients into (7) we obtain the approximate solution which is the same as the exact solution and the results are shown in Figure 2.
Example 2. Consider the following fractional Integro-differential equation:
subject to , where
with the exact solution .
Similarly as in Example 1 applying the least squares method with aid of shifted Chebyshev polynomial of the first kind , at , to the fractional Integro-differential equation (20) the numerical results are shown in Figures 3 and 4 and we obtain the approximate solution which is the same as the exact solution.
Example 3. Consider the following fractional Integro-differential equation:
subject to with the exact solution .
Similarly as in Examples 1 and 2 applying the least squares method with aid of shifted Chebyshev polynomial of the first kind , at , to the fractional Integro-differential equation (22) the numerical results are shown in Figures 5 and 6 and we obtain the approximate solution which is the same as the exact solution.
5. Conclusion
In this paper we study the numerical solution of three examples by using least squares method with aid of shifted Chebyshev polynomial which derives a good approximation. We show that this method is effective and has high convergency rate.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.