Abstract and Applied Analysis

Volume 2013 (2013), Article ID 413529, 12 pages

http://dx.doi.org/10.1155/2013/413529

Research Article

A Modified Generalized Laguerre Spectral Method for Fractional Differential Equations on the Half Line

D. Baleanu,1,2,3 A. H. Bhrawy,4,5 and T. M. Taha⁵

¹Department of Mathematics and Computer Sciences, Faculty of Arts and Sciences, Cankaya University, Eskisehir Yolu 29.Km, 06810 Ankara, Turkey

²Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, Jeddah 21589, Saudi Arabia

³Institute of Space Sciences, Magurele-Bucharest, Romania

⁴Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia

⁵Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef 62511, Egypt

Received 21 May 2013; Accepted 7 July 2013

Academic Editor: Soheil Salahshour

Copyright © 2013 D. Baleanu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper deals with modified generalized Laguerre spectral tau and collocation methods for solving linear and nonlinear multiterm fractional differential equations (FDEs) on the half line. A new formula expressing the Caputo fractional derivatives of modified generalized Laguerre polynomials of any degree and for any fractional order in terms of the modified generalized Laguerre polynomials themselves is derived. An efficient direct solver technique is proposed for solving the linear multiterm FDEs with constant coefficients on the half line using a modified generalized Laguerre tau method. The spatial approximation with its Caputo fractional derivatives is based on modified generalized Laguerre polynomials with , , and , and is the polynomial degree. We implement and develop the modified generalized Laguerre collocation method based on the modified generalized Laguerre-Gauss points which is used as collocation nodes for solving nonlinear multiterm FDEs on the half line.

1. Introduction

Due to the fact that fractional calculus has extensive applications in chemistry, engineering, neuron modeling, and biological sciences, the use of fractional calculus has attracted great interest from the mathematical science research community (e.g., see [1–5]).

As is well known, spectral methods use orthogonal polynomials as the basis functions and so usually provide accurate numerical results, see [6–9] and the references therein. In particular, we could employ tau and pseudospectral modified generalized Laguerre approximations and interpolation to solve multiterm fractional differential equations on semi-infinite interval. However, some authors developed the modified generalized Laguerre spectral method for the half line for ordinary, partial, and delay differential equations, see [8, 10–12]. The Laguerre polynomials of fractional orders and its generalization can be expressed by the Bell polynomials, see [13–15].

Because of the difficulty in obtaining exact solutions to many fractional differential equations, there is a large number of papers dealing with the numerical solution of fractional differential equation [16–21]. The usual spectral methods based on Legendre, Chebyshev, and Jacobi polynomials are only available for bounded domains for the numerical solution of FDEs [22, 23]. In [22, 24, 25], the authors discussed the operational matrices of Caputo fractional derivatives of shifted Chebyshev and Jacobi polynomials with spectral techniques to provide numerical approximations to the solution of linear and nonlinear multiterm FDEs on finite intervals. Ahmadian et al. [26] applied the Jacobi operational matrix to present an efficient algorithm for solving a class of fuzzy FDEs. Li et al. [27] employed spectral approximations to compute the fractional integral and the Caputo derivative; they also presented pseudospectral approximation for a class of FDEs. Yüzbaşı [28] proposed Bessel collocation method for the approximate solution of the Bagley-Torvik equation. Atabakzadeh et al. [29] extended the application of the shifted Chebyshev operational matrix for obtaining an analytical approximation solution for linear and nonlinear multiorder FDE. Moreover, the authors of [30] extended the application of the Legendre operational matrix for treating the fractional order fuzzy differential equations. Recently, Bhrawy et al. [31, 32] proposed the operational matrices of Riemann-Liouville fractional integration of Chebyshev and modified generalized Laguerre polynomials which are employed together with spectral tau method for solving linear FDEs on finite and semi-infinite intervals, respectively.

The fundamental aim of this paper is to extend the application of spectral tau and collocation methods based on modified generalized Laguerre polynomials to solve the linear and nonlinear multiterm fractional initial value problems on the half line. We propose the modified generalized Laguerre spectral tau (MGLT) approximation for obtaining a direct solution technique to solve linear multiterm FDE on the half line. This method requires the Caputo fractional derivatives of any fractional order of the modified generalized Laguerre polynomials of any degree which is already stated and proved.

The modified generalized Laguerre spectral collocation (MGLC) approximation, which is more reliable, is employed to obtain approximate solution of nonlinear multiterm FDE on the half line with leading fractional differential operator of order and initial conditions. To be more precise, in such approximation, the nonlinear FDE is collocated only at zeros of modified generalized Laguerre polynomials. Therefore, the problem reduces to a system of nonlinear algebraic equations which may be solved by any standard technique to find the spectral solution . To the best of the our knowledge, the proposed algorithms have not been introduced for the numerical solution of multiterm FDEs on the half line.

The plan of the paper is as follows. In the next section, we introduce basic properties of modified generalized Laguerre polynomials. In Section 3, the Caputo fractional derivative of the modified generalized Laguerre polynomials is proved. In Section 4, we develop modified generalized Laguerre tau and collocation spectral methods for solving multiterm FDEs. In Section 5, several numerical examples are implemented. Also, a conclusion is given in Section 6.

2. Properties of Modified Generalized Laguerre Polynomials

We recall below some relevant properties of the modified generalized Laguerre polynomials (see [8, 33, 34]).

Now, let and be a weight function on in the usual sense. Define

equipped with the following inner product and norm:

Next, let be the modified generalized Laguerre polynomial of degree for and , that is defined by

According to (2.3)-(2.4) of [11] for and , we have

where and .

The set of modified generalized Laguerre polynomials is the -orthogonal system, namely,

where is the Kronecker function and .

The modified generalized Laguerre polynomials of degree on the interval , are given by

where .

The special value

will be of important use later.

A function may be expressed in terms of modified generalized Laguerre polynomials as

In particular applications, the modified generalized Laguerre polynomials up to degree are considered. Then, we have

We now turn the modified generalized Laguerre-Gausstype quadratures, including modified generalized Laguerre-Gauss and Gauss-Radau interpolations [11, 34, 35]. Let and , , be the zeros of and , respectively. They are arranged in ascending order. Denote , , , the corresponding Christoffel numbers such that

For the modified generalized Laguerre-Gauss weights, we have

For the Gauss-Radau weights, we have

Note that the earlier two types of quadratures have close relations,

3. The Fractional Derivatives of

In this section, we prove the following theorem for expressing explicitly the fractional order derivatives of the modified generalized Laguerre polynomials in terms of the modified generalized Laguerre polynomials themselves.

Lemma 1. Let be a modified generalized Laguerre polynomial; then

Theorem 2. The Caputo fractional derivative of order for modified generalized Laguerre polynomials is given by

where

Proof. With the aid of (7), the -order derivative of modified generalized Laguerre polynomials is

The approximation of by terms of modified generalized Laguerre series yields

where is given by

Thanks to (18)–(20), we can write

where , and

4. Application of the Fractional Derivative of Modified Generalized Laguerre Polynomials

In this section, we consider spectral tau and collocation methods based on the fractional derivative of modified generalized Laguerre polynomials to solve numerically the linear and nonlinear multiterm FDEs on the half line.

4.1. Tau Method for Linear Multiterm FDEs

We apply the modified generalized Laguerre tau (MGLT) method for tackling the linear multiterm FDEs on the half line

subject to initial conditions

where are constant coefficients, , . While is the -order fractional derivative of , and is a source function.

The modified generalized Laguerre polynomials form a complete orthogonal system in weighted space; . Hence, we define

then the standard modified generalized Laguerre tau approximation to (23) is to find such that

Let us denote

then (26) can be written as

Let us also denote

where

then by virtue of the orthogonality relation of modified generalized Laguerre polynomials (6) and making use of (8) and (16) and after some manipulation, we may deduce that the elements of , , , and are given by

Therefore, (28) can be written in the following matrix form:

4.2. Collocation Method for Nonlinear Multiterm FDEs

It is known that the collocation method approximates differential equations in physical space, so it is easy to apply to various problems, including variable coefficient and nonlinear differential equations (see, for instance, [36]). In this section, an algorithm for solving fractional differential equation is proposed based on modified generalized Laguerre spectral collocation (MGLC) method in a half line. Here, we study the following multiterm FDE:

with initial conditions

where the fractional order derivatives are taken to be of the Caputo type, , , and is a nonlinear function.

We now derive an efficient algorithm for solving nonlinear multiterm FDEs (33)-(34). We expand the numerical approximation in terms of modified generalized Laguerre polynomials

By virtue of (16), the fractional derivatives can be expressed explicitly in terms of the expansion coefficients .

The modified generalized Laguerre spectral collocation method for solving approximately (33), (34) is to find such that

is collocated at , . In other words, we have to collocate (36) at the modified generalized Laguerre roots , which immediately yields

with (33) written in the form

Finally, the problem is reduced to solve a system of nonlinear algebraic equations, (37)-(38), in the unknown expansion coefficients , which can be solved by Newton’s iteration method.

5. Numerical Results

We report in this section some numerical results obtained with the proposed algorithms in the previous section. Several test examples are implemented to ensure the effectiveness of the proposed methods for fractional differential equations in the half line.

Example 3. This example concerns a class of Bagley-Torvik equation given by the initial value problem

The solution is .

This equation arises in the treatment of the motion of a large thin plate immersed in a Newtonian fluid. For the numerical solution of the Bagley-Torvik equation, see [37–41]. The solution of this problem is obtained by applying the modified generalized Laguerre tau method and the right-hand side is treated by modified generalized Laguerre-Gauss (G) and modified generalized Laguerre-Gauss-Radau (GR) quadratures with general parameters and . The maximum absolute errors obtained by MGLT method for various choices of , , and are given in Figures 1 and 2. Moreover, the approximate solutions obtained by the proposed method for , , and two choices of are shown in Figure 3 to make it easier to compare with the analytic solution.

Figure 1: Maximum absolute error by using MGLT (Gauss) with the various choices of .

Figure 2: Maximum absolute error by using MGLT (Radau) with the various choices of .

Figure 3: Comparing the exact and approximate solutions at , 40, with and .

Example 4. Consider the equation

where

and the exact solution is given by (see Table 1).

Table 1: Maximum absolute error, with various choices of , for Example 4 in which .

MGLT(G)
MGLT(GR)

MGLT(G)
MGLT(GR)

MGLT(G)
MGLT(GR)

10

30

50

70

10

30

50

70

10

30

50

70

Example 5. We next consider the following:

where

and the exact solution is given by .

The solution of this problem is obtained by applying the technique described in Section 4. The maximum absolute error for and various choices of , , , and are shown in Tables 2 and 3; moreover, the right-hand side is treated by modified generalized Laguerre Gauss and modified generalized Laguerre-Gauss-Radau quadratures.

Table 2: Maximum absolute error for , with various choices of , for Example 5 in which .

MGLT(G)
MGLT(GR)

MGLT(G)
MGLT(GR)

MGLT(G)
MGLT(GR)

10

30

50

70

10

30

50

70

10

30

50

70

Table 3: Maximum absolute error for , with various choices of , for Example 5 in which .

MGLT(G)
MGLT(GR)

MGLT(G)
MGLT(GR)

MGLT(G)
MGLT(GR)

10

30

50

70

10

30

50

70

10

30

50

70

Example 6. We consider the equation

subject to

The analytic solution of this problem is , we study two different cases of , , , , , , , and .

(I)Regarding , , , , , , and , Figures 4 and 5 display the maximum absolute errors using MGLT method with the treatment of right-hand side by modified generalized Laguerre-Gauss and Gauss-Radau, respectively.(II)Regarding , , , , , , and , Figures 6 and 7 display the maximum absolute errors using MGLT method with the treatment of right-hand side by modified generalized Laguerre-Gauss and Gauss-Radau, respectively.

Figure 4: Maximum absolute error (case  I) by using MGLT (Gauss) with the various choices of .

Figure 5: Maximum absolute error (case  I) by using MGLT (Radau) with the various choices of .

Figure 6: Maximum absolute error (case  II) by using MGLT (Gauss) with various choices of .

Figure 7: Maximum absolute error (case  II) by using MGLT (Radau) with various choices of .

Example 7. Consider the Bagley-Torvik equation with boundary conditions

where

and the exact solution is given by .

Bagley-Torvik equation involving fractional derivative of order arises in the modeling of the motion of a rigid plate in a Newtonian fluid and a gas in a fluid. We solve this equation with two-point boundary conditions where , , and by using modified generalized Laguerre-Gauss and Gauss-Radau quadratures with general parameters and . The maximum absolute errors obtained by MGLT method for various choices of , , and are given in Table 4. Moreover, the numerical and exact solutions obtained by the proposed method for , , and two choices of are shown in Figures 8-9 to make it easier to compare with the analytic solution.

Table 4: Maximum absolute error for , with various choices of , for Example 5 in which .

MGLT(G)
MGLT(GR)

MGLT(G)
MGLT(GR)

8

16

32

64

8

16

32

64

Figure 8: Comparing the exact and approximate solutions by using MGLT (Gauss) at , with and .

Figure 9: Comparing the exact and approximate solutions by using MGLT (Radau) at , with and .

Example 8. We next consider the following nonlinear initial value problem:

where

and the exact solution is given by .

The solution of this problem is obtained by applying the modified generalized Laguerre collocation method. The approximate solution obtained by the proposed method for , , , and two choices of is shown in Figure 10 to make it easier to compare with the analytic solution. Moreover, the absolute error for , , , and is given in Figure 11.

Figure 10: Comparing the exact and approximate solutions at , 10, 15, and 20, where , , and .

Figure 11: The absolute error for , , , and at .

Example 9. Consider the following nonlinear initial value problem:

whose exact solution is given by .

Comparison between the curves of exact solutions and the approximate solutions at and of proposed problem subject to for the four different fractional orders , 1.4, 1.6, and 1.8 in case of and are shown in Figures 12 and 13, respectively. The maximum absolute errors at for various choices of , , and in the interval are shown in Table 5. Moreover, Figures 14 and 15 display a comparison between the curves of exact solutions and the approximate solutions at and of proposed problem subject to for the four different fractional orders , 0.4, 0.6, and 0.8 in case of and , respectively. From all the Figures 12–15, it can be seen that the numerical solutions are in complete agreement with the exact solutions for all values of . Also, from the numerical results implemented in this example, the classical Laguerre polynomial (, ), which is used most frequently in practice, is not the best one, especially when we are approximating the solution of fractional differential equations.

Table 5: Maximum absolute error for , with various choices of , for Example 9 in which .

MGLC

MGLC

3

6

9

12

3

6

9

12

Figure 12: Comparing the exact and approximate solutions at , , , and , 1.4, 1.6, and 1.8.

Figure 13: Comparing the exact and approximate solutions at , , , and , 1.4, 1.6, and 1.8.

Figure 14: Comparing the exact and approximate solutions at , , , and , 0.4, 0.6, and 0.8.

Figure 15: Comparing the exact and approximate solutions at , , , and , 0.4, 0.6, and 0.8.

6. Conclusions

In this paper, we have proposed two efficient spectral methods based on modified generalized Laguerre polynomials for tackling linear and nonlinear FDEs on the half line. In these methods, the problem is reduced to the solution of a system of algebraic equations in the expansion coefficient of the solution. Numerical examples were given to demonstrate the validity and applicability of the method. The computational results show that the proposed methods can be effectively used in numerical solution of time-dependent fractional partial differential equation and other problems on the half line.

References

D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo, Fractional Calculus Models and Numerical Methods, Complexity, Nonlinearity and Chaos, World Scientific, Singapore, 2012.

K. Miller and B. Ross, An Introduction to the Fractional Calaulus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993.

M. Ortigueira, “Introduction to fraction linear systems. Part 1: continuous-time case,” in IEE Proceedings Vision, Image, Signal Processing, vol. 147, pp. 6–70, 2000.

K. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204, Elsevier, Amsterdam, The Netherlands, 2006.

C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods Fundamentals in Single Domains, Springer, Berlin, Germany, 2006.

E. H. Doha, A. H. Bhrawy, D. Baleanu, and S. S. Ezz-Eldien, “On shifted Jacobi spectral approximations for solving fractional differential equations,” Applied Mathematics and Computation, vol. 219, no. 15, pp. 8042–8056, 2013.

J.-P. Yan and B.-Y. Guo, “A collocation method for initial value problems of second-order ODEs by using Laguerre functions,” Numerical Mathematics. Theory, Methods and Applications, vol. 4, no. 2, pp. 283–295, 2011.

E. H. Doha and A. H. Bhrawy, “An efficient direct solver for multidimensional elliptic Robin boundary value problems using a Legendre spectral-Galerkin method,” Computers & Mathematics with Applications, vol. 64, no. 4, pp. 558–571, 2012.

D. Funaro, “Estimates of Laguerre spectral projectors in Sobolev spaces,” in Orthogonal Polynomials and Their Applications, C. Brezinski, L. Gori, and A. Ronveaux, Eds., pp. 263–266, Scientific Publishing, Singapore, 1991.

B.-Y. Guo, L.-L. Wang, and Z.-Q. Wang, “Generalized Laguerre interpolation and pseudospectral method for unbounded domains,” SIAM Journal on Numerical Analysis, vol. 43, no. 6, pp. 2567–2589, 2006.

M. Gülsu, B. Gürbüz, Y. Öztürk, and M. Sezer, “Laguerre polynomial approach for solving linear delay difference equations,” Applied Mathematics and Computation, vol. 217, no. 15, pp. 6765–6776, 2011.

R. Schimming and S. Z. Rida, “Noncommutative Bell polynomials,” International Journal of Algebra and Computation, vol. 6, no. 5, pp. 635–644, 1996.

A. M. A. El-Sayed and S. Z. Rida, “Bell polynomials of arbitrary (fractional) orders,” Applied Mathematics and Computation, vol. 106, no. 1, pp. 51–62, 1999.

R. A. Van Gorder, “Computation of certain infinite series of the form for arbitrary real-valued ,” Applied Mathematics and Computation, vol. 215, no. 3, pp. 1209–1216, 2009.

B. Shahriyar, F. Ismail, S. Aghabeigi, A. Ahmadian, and S. Salahshour, “An eigenvalueeigenvector method for solving a system of fractional differential equations with uncertainty,” Mathematical Problems in Engineering, vol. 2013, Article ID 579761, 11 pages, 2013.

A. H. Bhrawy, M. M. Tharwat, and Y. Yildirim, “A new formula for fractional integrals of Chebyshev polynomials: application for solving multi-term fractional differential equations,” Applied Mathematical Modelling, vol. 37, no. 6, pp. 4245–4252, 2012.

Y.-L. Jiang and X.-L. Ding, “Waveform relaxation methods for fractional differential equations with the Caputo derivatives,” Journal of Computational and Applied Mathematics, vol. 238, pp. 51–67, 2013.

S. Kazem, S. Abbasbandy, and S. Kumar, “Fractional-order Legendre functions for solving fractional-order differential equations,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 37, no. 7, pp. 5498–5510, 2013.

J. Song, F. Yin, X. Cao, and F. Lu, “Fractional variational iteration method versus Adomian's decomposition method in some fractional partial differential equations,” Journal of Applied Mathematics, vol. 2013, Article ID 392567, 10 pages, 2013.

X. Zhang, L. Liu, Y. Wu, and Y. Lu, “The iterative solutions of nonlinear fractional differential equations,” Applied Mathematics and Computation, vol. 219, no. 9, pp. 4680–4691, 2013.

A. H. Bhrawy, A. S. Alofi, and S. S. Ezz-Eldien, “A quadrature tau method for fractional differential equations with variable coefficients,” Applied Mathematics Letters, vol. 24, no. 12, pp. 2146–2152, 2011.

E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “Efficient Chebyshev spectral methods for solving multi-term fractional orders differential equations,” Applied Mathematical Modelling, vol. 35, no. 12, pp. 5662–5672, 2011.

E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order,” Computers & Mathematics with Applications, vol. 62, no. 5, pp. 2364–2373, 2011.

E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A new Jacobi operational matrix: an application for solving fractional differential equations,” Applied Mathematical Modelling, vol. 36, no. 10, pp. 4931–4943, 2012.

A. Ahmadian, M. Suleiman, S. Salahshour, and D. Baleanu, “A Jacobi operational matrix for solving a fuzzy linear fractional differential equation,” Advances in Difference Equations, 2013, article 104.

C. Li, F. Zeng, and F. Liu, “Spectral approximations to the fractional integral and derivative,” Fractional Calculus and Applied Analysis, vol. 15, no. 3, pp. 383–406, 2012.

S. Yüzbaşı, “Numerical solution of the Bagley-Torvik equation by the Bessel collocation method,” Mathematical Methods in the Applied Sciences, vol. 36, no. 3, pp. 300–312, 2013.

M. H. Atabakzadeh, M. H. Akrami, and G. H. Erjaee, “Chebyshev operational matrix method for solving multi-order fractional ordinary differential equations,” Applied Mathematical Modelling, 2013.

A. Ahmadian, M. Suleiman, and S. Salahshour, “An operational matrix based on Legendre polynomials for solving fuzzy fractional-order differential equations,” Abstract and Applied Analysis, vol. 2013, Article ID 505903, 29 pages, 2013.

A. H. Bhrawy and A. S. Alofi, “The operational matrix of fractional integration for shifted Chebyshev polynomials,” Applied Mathematics Letters, vol. 26, no. 1, pp. 25–31, 2013.

A. H. Bhrawy, M. Alghamdi, and T. M. Taha, “A new modified generalized Laguerre operational matrix of fractional integration for solving fractional differential equations on the half line,” Advances in Difference Equations, vol. 2012, p. 179, 2012.

B.-Y. Guo and X.-Y. Zhang, “A new generalized Laguerre spectral approximation and its applications,” Journal of Computational and Applied Mathematics, vol. 181, no. 2, pp. 342–363, 2005.

G. Szegő, Orthogonal Polynomials, vol. 23, American Mathematical Society. Colloquium Publications, 1985.

B.-Y. Guo and C. Zhang, “The spectral method for high order problems with proper simulations of asymptotic behaviors at infinity,” Journal of Computational and Applied Mathematics, vol. 237, no. 1, pp. 269–294, 2013.

A. H. Bhrawy and A. S. Alofi, “A Jacobi-Gauss collocation method for solving nonlinear Lane-Emden type equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 1, pp. 62–70, 2012.

K. Diethelm and N. J. Ford, “Numerical solution of the Bagley-Torvik equation,” BIT. Numerical Mathematics, vol. 42, no. 3, pp. 490–507, 2002.

S. S. Ray and R. K. Bera, “Analytical solution of the Bagley Torvik equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 168, no. 1, pp. 398–410, 2005.

Y. Çenesiz, Y. Keskin, and A. Kurnaz, “The solution of the Bagley-Torvik equation with the generalized Taylor collocation method,” Journal of the Franklin Institute, vol. 347, no. 2, pp. 452–466, 2010.

M. A. Z. Raja, J. A. Khan, and I. M. Qureshi, “Solution of fractional order system of Bagley-Torvik equation using evolutionary computational intelligence,” Mathematical Problems in Engineering, vol. 2011, Article ID 675075, 18 pages, 2011.

S. S. Ray, “On Haar wavelet operational matrix of general order and its application for the numerical solution of fractional Bagley Torvik equation,” Applied Mathematics and Computation, vol. 218, no. 9, pp. 5239–5248, 2012.