Variational iteration method for fractional calculus-a universal approach by Laplace transform

*Correspondence: wuguocheng2002@yahoo.com.cn; dumitru@cankaya.edu.tr 1College of Mathematics and Information Science, Neijiang Normal University, Neijiang, 641112, China 3Department of Mathematics and Computer Sciences, Cankaya University, Balgat, Ankara 06530, Turkey Full list of author information is available at the end of the article Abstract A novel modification of the variational iteration method (VIM) is proposed by means of the Laplace transform. Then the method is successfully extended to fractional differential equations. Several linear fractional differential equations are analytically solved as examples and the methodology is demonstrated. MSC: 39A08; 65K10; 34A12


Introduction
The Lagrange multiplier technique [] was widely used to solve a number of nonlinear problems which arise in mathematical physics and other related areas, and it was developed into a powerful analytical method, i.e., the variational iteration method [, ] for solving differential equations. The method has been applied to initial boundary value problems [-], fractal initial value problems [, ], q-difference equations [] and fuzzy equations [-], etc.
Generally, in applications of VIM to initial value problems of differential equations, one usually follows the following three steps: (a) establishing the correction functional; (b) identifying the Lagrange multipliers; (c) determining the initial iteration. The step (b) is very crucial. Applications of the method to fractional differential equations (FDEs) mainly and directly used the Lagrange multipliers in ordinary differential equations (ODEs) which resulted in poor convergences. This point of view needs some explanations will elucidate the target of the suggested improvement, among them: () When the Riemann-Liouville (RL) integral emerges in the constructed correctional functional, the integration by parts is difficult to apply; () To avoid this problem, the RL integral is replaced by an integer one which allows the integration by parts. This is a very strong simplification but it affects the next steps of the application of the method; () Therefore, the Lagrange multiplier is determined by a simplification not reasonably explained in the literature, so far. To overcome this drawback, the present article conceives a method how the Lagrange multiplier has to be defined from Laplace transform. The technique can be readily and http://www.advancesindifferenceequations.com/content/2013/1/18 universally extended to solve both differential equations and FDEs with initial value conditions.

Basics of the variation iteration method
In order to illustrate the basic idea of the technique, consider the following general nonlinear system: where u = u(t), R is a linear operator, N is a nonlinear operator and g(t) is a given continuous function and d m u/dt m is the term of the highest-order derivative. The basic character of the method is to construct the following correction functional for Eq. (): where λ(t, τ ) is called the general Lagrange multiplier [-] and u n is the nth order approximate solution.
According to the VIM's rules [, , ], readers may note that the integration by parts plays an important role in the derivation of the Lagrange multipliers. But in fractional calculus, generally, the following integration by parts cannot hold (as it is mentioned in point () of the preceding section): where v = v(t), C  D α t and  I α t are the notations of the Caputo derivative and the RL integration, respectively. That's why the VIM was not so successful as other analytical methods such as the Adomian decomposition method (ADM) [-] and the homotopy perturbation method (HPM) [-] in fractional calculus. For this reason, we consider the following reconstruction of the method using the Laplace transform.

New identification of the Lagrange multipliers
Let us revisit the original idea of the Lagrange multipliers in the case of an algebraic equation. Firstly, an iteration formula for finding the solution of the algebraic equation f (x) =  can be constructed as The optimality condition for the extreme δx n+ δx n =  leads to where δ is the classical variational operator. From () and (), for a given initial value x  , we can find the approximate solution x n+ by the iterative scheme for () This algorithm is well known as the Newton-Raphson method and has quadratic convergence. Now, we extend this idea to finding the unknown Lagrange multiplier. The main step is to first take the Laplace transform to Eq. (). Then the linear part is transformed into an algebraic equation as follows: The iteration formula of () can be used to suggest the main iterative scheme involving the Lagrange multiplier as as restricted terms, one can derive a Lagrange multiplier as With Eq. () and the inverse-Laplace transform L - , the iteration formula () can be explicitly given as where the initial iteration u  (t) can be determined by Eq. () also explained why the initial iteration in the classical VIM is determined by the Taylor series. This modified VIM here transfers the problem into the partial differential equation in the Laplace s-domain and removes the differentiation with respect to time. This idea has been used in other analytical methods such as the Laplace ADM [, ] and the Laplace HPM [], respectively.

Illustrative examples
We now consider the applications of the modified VIM to both ODEs and FDEs. http://www.advancesindifferenceequations.com/content/2013/1/18

Ordinary differential equations
Example  Consider the following simple line differential equation: which has the exact solution u(t) = u  e -t .
We can obtain the successive approximate solutions as . . .
For n → ∞, u n (t) tends to the exact solution u  e -t .

Example  The logistic differential equation []
has the exact solution u(t) = e t +e t . By the present VIM, we have the following solutions: . . . There can be various choices of u  (t) and λ(s) which affect the speed of the convergence. We note that the integration by parts is not used and the calculation of the Lagrange multiplier here is much simpler. Furthermore, the VIM can be easily extended to FDEs and this is the main purpose of our work.

Fractional differential equations
In the early application of VIM [] to FDEs, the term C  D α t u is considered as a restricted variation, i.e., and the variational iteration formula is given as and g(τ , u n ) is a nonlinear term. But for the following FDEs, the above popular applications of the VIM were not successful: Now, we consider the application of the modified VIM.
The following Laplace transform [-] of the term C  D α t u holds: Taking the above Laplace transform to both sides of (), the iteration formula of Eq. () can be constructed as As a result, after the identification of a Lagrange multiplier λ = - s α , one can derive Let us apply the above VIM to solve FDEs of Caputo type.
After taking the Laplace transform to both sides of Eq. (), we get the following iteration formula: (   ) Setting L[u n (t)] as a restricted variation, λ(s) can be identified as The approximate solution of Eq. () can be given as which reads . . . u n (t) rapidly tends to the exact solution of Eq. () for n → ∞.
Example  Consider the fourth example, the time-fractional diffusion equation The and λ(s) = - s α is used as earlier. http://www.advancesindifferenceequations.com/content/2013/1/18 As a result, the successive approximation can be obtained as follows: The exact solution can be given in a compact form The method's efficiency for a nonlinear differential equation with variable coefficients is illustrated in []. For other applications of a new modified VIM to ODEs and FDEs, readers are also referred to [-].

Remarks
(a) The conceived modification of the VIM is a universal approach to both ODEs and FDEs. As a result, it becomes possible to design a 'universal' software package in future work. (b) Now one can consider implementing other linearized techniques, i.e., the Adomian series and the homotopy series to handle the nonlinear terms and improve the accuracy of the approximate solutions. (c) This modified VIM can also be used to solve the FDEs of RL type.

Conclusions
A new approach is proposed to identify the Lagrange multipliers of the VIM and a concept of the Laplace-Lagrange multipliers is proposed from the Laplace transform. Especially for the FDEs, to the best of our knowledge, there is no effective method to identify the Lagrange multipliers. With the approach given in this paper, we can easily derive Lagrange multipliers without tedious calculation and new variational iteration formulae can be derived. Some FDEs with the Caputo derivatives are illustrated. The results show the modified method's efficiency compared with other versions of the VIM in fractional calculus. http://www.advancesindifferenceequations.com/content/2013/1/18