Solving a Nonlinear Multi-Order Fractional Differential Equation Using Legendre PseudoSpectral Method

In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multiorder fractional differential equations (M-FDEs). The fractional derivative is described in the Caputo sense. The study is conducted through illustrative example to demonstrate the validity and applicability of the presented method. The results reveal that the proposed method is very effective and simple. Moreover, only a small number of shifted Legendre polynomials are needed to obtain a satisfactory result.


Introduction
Many phenomena in engineering physics, chemistry, and other sciences can be described very successfully by models that use mathematical tools of fractional calculus, i.e. the theory of derivatives and integrals of non-integer order [1][2][3].For example, they have been successfully used in modeling frequency dependent damping behavior of many viscoelastic materials.There are numerous research which demonstrate the applications of fractional derivatives in the areas of electrochemical processes, dielectric polarization, colored noise, and chaos.
The numerical solution of differential equations of integer order has been a hot topic in numerical and computational mathematics for a long time.The solution of fractional differential equations has been recently studied by numerous authors.However, the state of the art is far less advanced for general fractional order differential equations.Moreover, to the best of the authors knowledge, very few algorithms for the numerical solution of Multi-order Fractional Differential Equations (M-FDEs) have been suggested [4][5][6], particularly algorithms for analytical solutions and approximate solutions of nonlinear M-FDEs.
As we know, the fractional derivatives are global dependence problems (they are definite by the integral in [0, t]), from this point, the global methods-spectral methods maybe are more suit to solve the FDEs.As we know the standard spectral methods possess infinite order of accuracy for the equations with well regularity solutions, while fail to many complicated problems with singular solutions.So, it is attracted considerable interest that how to enlarge the adaptability of spectral methods, and construct certain simple approximation schemes without loss of any accuracy to even more complicated problems.It is well known that Legendre polynomials are well known family of orthogonal polynomials on the interval [−1, 1].They are widely used because of their good properties in the approximation of functions [7][8][9].For multi-order fractional differential equation, a operational matrix is studied [4], Chebyshev wavelets is considered in [10], adomian decomposition is considered in [11], a variational iteration method is considered in [12], predictor-corrector method is studied in [13], operator splitting method is considered in [14], and Adams method is researched in [15].In this paper, we introduce the Legendre pseudo-spectral method to solve multi-order arbitrary differential equations, which include the linear and nonlinear differential equations.
The outline of this paper is as follows.In Section 2, we review the basic definitions and the properties of the fractional calculus.In Section 3, the approximation of fractional derivative by Legendre Polynomials is obtained.In Section 4, we present the application of the Legendre pseudo-spectral method to multi-order fractional differential equation.Some numerical examples are provided in Section 5. Also a conclusion is given in Section 6.

Fractional Calculus
In this section, we first review the basic definitions and the operational properties of fractional integral and deri-vative for the purpose of acquainting with sufficient fractional calculus theory.Many definitions and studies of fractional calculus have been proposed in the past two centuries.These definitions include Riemann-Liouville, Reize, Caputo and Grnwald-Letnikov fractional operators.The two most commonly used definitions are the Riemann-Liouville operator and the Caputo operator.We give some definitions and properties of the fractional calculus.
c D  denotes Caputo fractional derivative of order  is defined where I  denotes the Riemann-Liouville fractional integral of order order  defined as As we all know, there are some different definitions of fractional operator except the Caputo fractional derivative.From a theoretical point of view the most natural approach is the Riemann-Liouville definition defined as The relationship between the Caputo definition and the Riemann-Liouville definition can be given by the following formula Similar to integer-order differentiation, Caputo's fractional differentiation is a linear operation: For the Caputo and Riemann-Liouville's derivative we have: Amongst a variety of definition for fractional order derivatives, Caputo fractional derivative has been used because it allows physically interpretable initial conditions.

Evaluation of the Fractional Derivative Using Legendre Polynomials
The well known Legendre polynomials are defined on the interval   1, 1  and can be determined with the aid of the following recurrence formula: where . In order to use these polynomials on the interval  0, t L we define the so called shifted Legendre polynomials by introducing the change of variable   Then   i P t can be obtained as follows: where . The analytic form of the shifted Legendre polynomials   i P t of degree given by: i and .The orthogonality condition is: The function   y t , square integrable in   0, L , may be expressed in terms of shifted Legendre polynomials as: where the coefficients are given by: In practice, only the first   1 m  -terms shifted Legendre polynomials are considered.Then we have: In the following theorem we introduce an approximate formula of the fractional derivative of   y t .
Theorem 1 Let   y t be approximated by shifted Legendre polynomials as (10) and also suppose 0 Proof.Since the Caputo's fractional differentiation is a linear operation we have: Employing Equations (3) in Equation ( 6) we have: , by using Equations (3) in Equation ( 6) we get:

Solution of Multi-Order Fractional Differential Equations
Consider the multi-order fractional differential equations of type given in (13): with initial conditions: where denotes the Caputo fractional derivative of order a.It should be noted that F can be nonlinear in general.To solve problems ( 13) and ( 14 From ( 13), (15) and Theorem 1, we have: we now collocate (16) at 1 m        points p t as: For suitable collocation points t p , we use 1 m Also, by substituting (15) in the initial conditions ( 14) we can obtain .

     equations as follows:
ations of the boundary conditions (18), generate   1 m  equations which can be solved using Newton's iterative method for the unknown .Consequently , 0,1, , y t can be calculated.

Numerical Examples
In this section, we will use the Legendre collocation method to solve nonlinear fractional (arbitrary) order differential equation.These examples are considered because closed form solutions are available for them, or they have also been solved using other numerical schemes.This allows one to compare the results obtained using this scheme with the analytical solution or the solutions obtained using other schemes.

Example 1. Consider the following initial value problem
The exact solution of (20) is .1.006 0.0042 0.001 .
My result is in good agreement with exact solution.This demonstrates the importance of my numerical scheme in solving nonlinear multi-order fractional differential equations.
Example 3. Following Odibat and Momani [6], we consider fractional Riccati equation subject to the initial state , which is studied by Odibat [6] by using the modified homotopy perturbation method and Li [10] by using the Chebyshev wavelet operational matrices method.Here we use the Legendre pseudo-spectral method to solve it.
  0 0 y  This is a nonlinear system of algebraic equations.The numerical solution, for m = 8, is shown in Figure 1.The exact solution of this problem, when 1 and we can observe that, as t ,     1 2 y t   .From Figure 1 we can see the numerical solution is very good agreement with the exact solution when 1 , the numerical solution is very good agreement with the result in Ref. [10].Therefore, we hold that the solution for 0.5   and 0.75   is also credible.Numerical results with comparison to Ref. [6] and [10] are given in Table 1 on the interval   0, 1 .The difference between our result or in Ref. [10] and the result in Ref. [6] is obvious.Because only the fourth-order term of the homotopy perturbation solution were used in evaluating the approximate solutions in Ref. [6].We hold that our results are better for 0.5   and 0.75

 
. Compared to the Chebyshev wavelet operational matrices method in Ref. [10] with 192 degrees of freedom, our method reached the same accuracy used 9 degrees of freedom only.

Conclusion
Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations.For that reason we need a reliable and efficient technique for the solution of  fractional differential equations.This paper deals with the approximate solution of a class of multi-order fractional differential equations.The fractional derivatives are described in Caputo sense.Our main aim is to evaluation of fractional derivative using Legendre polynomials and implementing it to solve the nonlinear multiorder fractional differential equations.The main characteristic behind the approach using this technique is that it reduces such problems to those of solving a system of algebraic equations thus greatly simplifying the problem.Illustrative example is included to demonstrate the validity and applicability of the presented technique.The comparison certifies that our method gives good results.

Figure 1 .
Figure 1.The behavior of the exact and approximate solution of example with  m 8 of Example 3.
) we approximate y t as: