Research Article  Open Access
Generalized Taylor Series Method for Solving Nonlinear Fractional Differential Equations with Modified RiemannLiouville Derivative
Abstract
We propose an efficient analytic method for solving nonlinear differential equations of fractional order. The fractional derivative is defined in the sense of modified RiemannLiouville derivative. A new technique for calculating the generalized Taylor series coefficients (also known as “generalized differential transforms,” GDTs) of nonlinear functions and a new approach of the generalized Taylor series method (GTSM) are presented. This new method offers a simple algorithm for computing GDTs of nonlinear functions and avoids massive computational work that usually arises in the standard method. Several illustrative examples are demonstrated to show effectiveness of the proposed method.
1. Introduction
Fractional differential equations are generalizations of classical differential equations of integer order and have recently proved to be valuable tools in the modeling of many phenomena in various fields of science and engineering. Apart from diverse areas of mathematics, fractional differential equations arise in rheology, viscoelasticity, chemical physics, electrical networks, fluid flows, control, and dynamical processes in selfsimilar and porous structures. There has appeared lots of work in which fractional derivatives are used for a better description of considered material properties; mathematical modelling based on enhanced rheological models naturally leads to differential equations of fractional order and to the necessity of the formulation of initial conditions to such equations. Several numerical methods for solving fractional differential equations have been introduced lately. The authors in [1] presented the predictorcorrector approach based on the AdamBashforthMoulton type numerical method that has been successful in obtaining stable approximations for solving many fractional differential equations. Some of the semianalytic methods such as the Adomian decomposition method (ADM) [2, 3], homotopy analysis method (HAM) [4–6], homotopy perturbation method (HPM) [7, 8], variational iteration method (VIM) [9, 10], and generalized differential transform method (GDTM) [11–13] have been introduced to provide analytic or numeric approximations.
In this paper we focus on the generalized Taylor series method (GTSM), which is based on the generalized Taylor series. The fractional derivative is defined in the sense of the modified RiemannLiouville derivative [14]. From the given fractional differential equation, the GTSM provides a simple recurrence relation of the generalized Taylor series coefficients of the solution. We obtain recurrence relations of complex nonlinear functions such as the exponential, logarithmic, and trigonometric functions. The paper is organized as follows. Section 2 introduces some preliminary results from the fractional calculus that we will use. The basic idea and some properties of GTSM are presented in Section 3. We present the recurrence relations of complex nonlinear functions in Section 4. In Section 5, numerical results of several examples are demonstrated by using new recurrence relations. Finally, we give a conclusion in Section 6.
2. Preliminary Results
In [14], Jumarie proposed a definition for fractional derivative which is known as the modified RiemannLiouville derivative in the literature. Since then, many authors have investigated various applications of the modified RiemannLiouville derivative (e.g., see [15–17]) including various fractional calculus formulae, the fractional variational iteration method, and the fractional subequation method for solving fractional partial differential equations. The definition and some of the key properties of modified RiemannLiouville derivative arewhich do not hold for classical RiemannLiouville and Caputo derivatives. Particularly chain rule (also known as “Faà di Bruno’s formula” in fractional calculus) plays the key role in our method.
3. Generalized Taylor Series Method
The Taylor series method (TSM) (also known as differential transform method, DTM) that is based on the Taylor series has been successful in achieving accurate approximate solutions for the linear and nonlinear problems. In TSM, all coefficients of Taylor series of the solution can be determined by solving the recurrence equations induced from the given differential equation. The authors in [18] developed the fractional differential transform method (FDTM) which is based on the classical TSM and generalized with fractional derivative. The generalized Taylor formula is introduced in [11] and has been used to develop GDTM in [12]. The authors in [19] proposed a new algorithm for calculating the differential transforms of several nonlinear functions based on the chain rule. The purpose of this paper is to obtain efficient algorithms to calculate generalized Taylor series coefficients of complex nonlinear functions.
For the fractional differential operator , , in the sense of modified RiemannLiouville definition, let us define ( times).
Theorem 1 (generalized Taylor’s formula [11]). Suppose that for , where ; then one has for all with .
For an analytic function , let us define the generalized differential transform (GDT) of the th derivative as follows:where , , and the generalized differential inverse transform of is defined as follows:Roughly speaking, coefficient of the th term in the generalized Taylor series of the function is called the GDT of the th derivative of . In case of , the GDT reduces to the classical differential transform. Some of the fundamental properties of GDT are listed below.
Theorem 2. Suppose that , , and are the generalized differential transforms of the functions , , and , respectively. Then the following properties are satisfied.(1)If , then .(2)If , where is a constant, then .(3)If , then .(4)If , then .
The proof of these results for Caputo derivative can be found in [12] and clearly holds for modified RiemannLiouville derivative also.
4. Calculating the Coefficients of Generalized Taylor Series for Complex Nonlinear Functions
In this section we will introduce an efficient algorithm to calculate generalized Taylor series coefficients (GDTs) of several complex nonlinear functions. Since the chain rule holds for the modified RiemannLiouville derivative, the algorithms for calculating the GDTs of typical nonlinear functions given in [19] can be totally adopted to nonlinear fractional differential equations with modified RiemannLiouville derivative.
4.1. GDTs for the Exponential and Logarithmic Nonlinearity
Let us consider , where is a constant. From the definition of transform,Now taking fractional derivative of with respect to variable , we haveWith the application of Theorem 2 to (6) to obtain the generalized Taylor series coefficients, we have Now replacing with givesfor . Combining (5) and (8), we obtain the recursive relationship for calculating GDTs of :
Now consider , where and are constants and is a continuous function satisfying for all . By definition of transform, we haveBy taking fractional derivative of with respect to variable , we haveor equivalentlyWith the application of the generalized differential transform to (12) we have, from Theorem 2,Replacing with givesSubstitute in (14) to getFor , (14) can be rewritten asBy combining (10), (15), and (16), we obtain the recursive relationship for calculating GDTM of :
4.2. GDTMs for the Trigonometric and Hyperbolic Nonlinearity
Now consider and , where and are constants. By definition of GDT,Further, taking fractional derivative of (18), we haveSimilarly, applying GDT with the aid of Theorem 2, we haveReplacing with givesCombine (18) and (21) to obtain recurrence relation:For the functions and , we immediately obtain the recursive relationships:
4.3. Error Analysis
From Theorem 1, if for , where , then we have for all . Furthermore, there is a value with so that the error term has the form So, if is GDT of the th derivative of the analytic function , we have and the error term has the form for some . Absolute value of the error done with approximation can be calculated by maximizing the right hand side of following inequality:for . The accuracy of increases when we choose large and decreases as the value of moves away from the center . Hence, we must choose large enough so that the error does not exceed a specified bound.
5. Numerical Illustrations
Example 1. Consider the following nonlinear fractional differential equation:for and , subject to initial condition . Applying GTSM to (30) with the aid of Theorem 2, we have the following recurrence relation:where is the GDT of and is given by (9) asfor . From the initial condition , we have and .
Solving the recurrence relations (31) and (32) together gives all GDTs. Table 1 shows the values of and for , , and . We set the number of iteration which gives the approximation and the absolute error bound . We see that the speed of convergence of is sufficient to obtain accurate approximations even in 5 iterations. In Table 2, we summarize these approximations of solution for and .
Approximations of the solution for different order of fractional derivatives are given in Figure 1. We set iteration number and summarize approximate solutions for .


Example 2. Consider the nonlinear fractional differential equationfor and , with the initial data . Applying GTSM to (33), we have the following recurrence relation:where is the GDT of and is given by (17) asFrom the initial data , we have and . By using this in (34), we have And by using this in (35) we obtain So we haveSolving the recurrence relations (34) and (38) together gives all GDTs. Table 3 shows the values of and for , , and . From Table 4, we see that the speed of convergence of is sufficient to obtain accurate approximations with 10 iterations. We set the number of iteration which gives the approximation and the error bound . We summarize the approximations of solution for and in Figure 2 and Table 5.



Example 3. Consider the fractional differential equationfor and , with the initial data . Applying GTSM to (39), we have the following recurrence relation:where is the GDT of and is given by (22) asand is the GDT of which is given by (23) asfor .
From the initial data , we have , , and . Solving the recurrence relations (40), (41), and (42) together gives all GDTs. Table 5 shows the values of and for , , and . Table 6 seems a bit complicated and does not give hint about speed of convergence of approximations; however, we see from Table 7 that the speed of convergence of is sufficient to obtain accurate approximations even with 5 steps of iteration. We set the number of iteration which gives the approximation and the error bound . Approximations of the solution for different order of fractional derivatives are given in Table 8 and Figure 3. We set iteration number and summarize approximate solutions for .



Example 4. Consider the nonlinear fractional differential equationfor and , with the initial data . Applying GTSM to (43), we have the following recurrence relations:where and are the GDTs of and , which are given byfor . From the initial data , we have , , , , and . Solving the recurrence relations (44), (45) together gives all GDTs.
Table 9 gives the values of , , , , and for , , and . We see from Table 10 that the speed of convergence of is sufficient to obtain accurate approximations with 10 steps of iteration. Approximations of the solution for different order of fractional derivatives are given in Figure 4. We set iteration number and summarize approximate solutions for . The numerical values are also given in Table 11.



6. Conclusion
In this study we proposed a generalized Taylor series method (GTSM) based on the generalized Taylor formula for solving nonlinear differential equations of fractional order. We consider the fractional differential equations with modified RiemannLiouville derivative. In GTSM, it is key to obtain recursive relationships from the given differential equation. However, the recurrence relations for the complex nonlinear functions for fractional differential equations with modified RiemannLiouville derivative have not been derived before. The GTSM with new recurrence relations given in this study offers a simple algorithm to compute the generalized Taylor series coefficients. Thus, it is worthwhile to mention that the GTSM is a straightforward, promising, and powerful method for solving nonlinear fractional differential equations.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
References
 K. Diethelm, N. J. Ford, and A. D. Freed, “A predictorcorrector approach for the numerical solution of fractional differential equations,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 3–22, 2002. View at: Publisher Site  Google Scholar
 S. S. Ray and R. K. Bera, “An approximate solution of a nonlinear fractional differential equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 167, no. 1, pp. 561–571, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Hu, Y. Luo, and Z. Lu, “Analytical solution of the linear fractional differential equation by Adomian decomposition method,” Journal of Computational and Applied Mathematics, vol. 215, no. 1, pp. 220–229, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 I. Hashim, O. Abdulaziz, and S. Momani, “Homotopy analysis method for fractional IVPs,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 674–684, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 A. K. Alomari, M. S. M. Noorani, R. Nazar, and C. P. Li, “Homotopy analysis method for solving fractional lorenz system,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 7, pp. 1864–1872, 2010. View at: Publisher Site  Google Scholar
 M. Zurigat, S. Momani, Z. Odibat, and A. Alawneh, “The homotopy analysis method for handling systems of fractional differential equations,” Applied Mathematical Modelling, vol. 34, no. 1, pp. 24–35, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 O. Abdulaziz, I. Hashim, and S. Momani, “Application of homotopyperturbation method to fractional IVPs,” Journal of Computational and Applied Mathematics, vol. 216, no. 2, pp. 574–584, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 O. Abdulaziz, I. Hashim, and S. Momani, “Solving systems of fractional differential equations by homotopyperturbation method,” Physics Letters A, vol. 372, no. 4, pp. 451–459, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 W. Guocheng, “A fractional variational iteration method for solving fractional nonlinear differential equations,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 2186–2190, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 S. Das, “Analytical solution of a fractional diffusion equation by variational iteration method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 483–487, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 Z. M. Odibat and N. T. Shawagfeh, “Generalized Taylor's formula,” Applied Mathematics and Computation, vol. 186, no. 1, pp. 286–293, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Odibat, S. Momani, and V. S. Erturk, “Generalized differential transform method: application to differential equations of fractional order,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 467–477, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 V. S. Erturk and S. Momani, “Solving systems of fractional differential equations using differential transform method,” Journal of Computational and Applied Mathematics, vol. 215, no. 1, pp. 142–151, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 G. Jumarie, “Modified RiemannLiouville derivative and fractional Taylor series of nondifferentiable functions further results,” Computers & Mathematics with Applications, vol. 51, no. 910, pp. 1367–1376, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 G. Jumarie, “Table of some basic fractional calculus formulae derived from a modified Riemann–Liouville derivative for nondifferentiable functions,” Applied Mathematics Letters, vol. 22, no. 3, pp. 378–385, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 N. Faraz, Y. Khan, H. Jafari, A. Yildirim, and M. Madani, “Fractional variational iteration method via modified RiemannLiouville derivative,” Journal of King Saud University—Science, vol. 23, no. 4, pp. 413–417, 2011. View at: Publisher Site  Google Scholar
 B. Lu, “Backlund transformation of fractional Riccati equation and its applications to nonlinear fractional partial differential equations,” Physics Letters A, vol. 376, no. 2829, pp. 2045–2048, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 A. Arikoglu and I. Ozkol, “Solution of fractional differential equations by using differential transform method,” Chaos, Solitons and Fractals, vol. 34, no. 5, pp. 1473–1481, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S.H. Chang and I.L. Chang, “A new algorithm for calculating onedimensional differential transform of nonlinear functions,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 799–805, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2015 Süleyman Öğrekçi. 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.