In this paper, the exact solutions of space-time fractional telegraph equations are given in terms of Mittage-Leffler functions via a combination of Laplace transform and variational iteration method. New techniques are used to overcome the difficulties arising in identifying the general Lagrange multiplier. As a special case, the obtained solutions reduce to the solutions of standard telegraph equations of the integer orders.
1. Introduction
Fractional differential equations are widely used in many branches of sciences. Many phenomena in engineering, physics, chemistry and other sciences can be described successfully using fractional calculus. Nonlinear oscillation of earthquake, acoustics, electromagnetism, electrochemistry, diffusion processes and signal processing can be modeled by fractional equations [1–5]. The telegraph equations have a wide variety of application in physics and engineering. The applications arise, for example, in the propagation of electrical signals and optimization of guided communication systems [6–9]. It is recently shown by Arbab [10] that a quaternionic momentum eigenvalue produces a telegraph equation. This equation is found to describe the propagation of a quantum particle.
The fractional telegraph equations have been investigated by many authors in recent years. Garg et al. [9] derived a solution of space-time fractional telegraph function in a bounded domain by the method of generalized differential transform and obtained the solution in terms of Mittage-Leffler functions. Chen et al. [11] obtained the solution of nonhomogenous time-fractional telegraph equation with nonhomogenous boundary conditions, namely, Dirichlet, Neumann, and Robin boundary conditions using the method of separation of variables. The solutions are given in the form of the multivariate Mittage-Leffler functions. Ansari [12] derived a formal solution of the time-fractional telegraph equation by applying a fractional exponential operator. Huang [13] considered the time-fractional telegraph equation for the Cauchy problem and signaling problem. He solved the problem by the combined Fourier-Laplace transforms. Also, Huang derived the solution for the bounded problem in a bounded-space domain by means of Sine-Laplace transforms methods. Das et al. [14] used a homotopy analysis method in approximating an analytical solution for the time-fractional telegraph equation and different particular cases have been derived. Jiang and Lin [15] obtained the solution in a series form for the time-fractional telegraph equation with Robin boundary value conditions using the reproducing kernel theorem. Khan et al. [16] used a method based on perturbation theory and Laplace transformation for solving space-time fractional telegraph equations. They considered fractional Taylor series and fractional initial conditions in deriving the solution. Sevimlican [6] considered a one-dimensional space fractional telegraph equations by the variation iteration method; he found the general Lagrange multiplier to be λ=(ξ-x). But, as mentioned by He [17] the exact identification of the general Lagrange multiplier is impossible for most problems and an approximate identification is always followed. He [17], approximated the Lagrange multiplier (as λ=-(ξ-x)) for a one-dimensional space-fractional telegraph equations.
Recently, a method that combined the Laplace transform and variational iteration method (LVIM) has been introduced by many authors in solving various types of problems. Abassy et al. [18] used a combination of variational iteration method, Laplace transform, and Pade’ technique in obtaining solution to nonlinear equations in compact form. Hammouch and Mekkaoui [19] approximated the solutions of a homogenous Smoluchowski coagulation equation by Laplace variational iteration method. Arife and Yildirim [20] developed Laplace variational iteration method (LVIM) for solving eighth-order equations.
In this paper, the authors extend Laplace variational iteration method (LVIM) and apply it to space-time one-dimensional fractional telegraph equations in a half-space domain (signaling problem). This approach enables us to overcome the difficulties that arise in finding the general Lagrange multiplier.
In Section 2, we provide some preliminaries. Section 3 introduces the concept of variational iteration method, while Section 4 illustrates the construction of Laplace variational technique. In Section 5 the authors provide numerical examples. The conclusions are given in Section 6.
2. PreliminariesDefinition 1.
The Caputo fractional derivative of order α>0 of a function f(x), x>0 is defined by [5, 21]
(1)D0xαf(x)={1Γ(n-α)∫0x(x-t)n-α-1f(n)(t)dt,n-1<α≤n∈ℕdndxnf(x),α=n∈ℕ,
where D0xα is called the Caputo derivative operator.
Note 1.
From Definition 1, the following result is obtained:
(2)D0tαtβ={Γ(β+1)Γ(β-α+1)tβ-α,n-1<α≤n,β>n-1,β∈ℛ0,n-1<α≤n,β≤n-1.
Definition 2.
The Laplace transform of fractional order derivative, is defined by [1–3, 21, 22]
(3)ℒ[D0xαf(x)]=sαℒ[f(x)]-∑k=0n-1sα-k-1[f(k)(x)]x=0,n-1<α≤n,n∈ℕ.
Definition 3.
The Mittag-Leffler function with two parameters is defined by [21, 23, 24]
(4)Eα,β(z)=∑n=0∞znΓ(αn+β),α,β,z∈ℂ,Re(α)>0,Re(β)>0.
Note 2.
It follows Definition 3 that
(5)(1)E2,1(x2)=cosh(x),(2)E2,2(x2)=sinh(x)x,(3)E2,3(x2)=1x2[-1+cosh(x)].
Definition 4.
The generalized Mittag-Leffler function is defined by [23, 24]
(6)Eα,βγ(z)=∑k=0∞(γ)kk!Γ(kα+β),α,β,γ∈ℂ,Re(α)>0,Re(β)>0,
where (γ)k=γ(γ+1),…,(γ+k-1)=Γ(γ+k)/Γ(γ). For γ=1, Eα,βγ(z) reduces to Mittage-Leffler function (4).
3. Variational Iteration Method
He [25] developed the variational iteration method (VIM) that is widely used to evaluate either exact or approximate solutions of linear and nonlinear problems [17, 26–28]. The variational iteration method gives the solution in a rapidly infinite convergent series. To illustrate the concept of VIM, we consider the following general nonlinear equation with prescribed auxiliary conditions:
(7)Lu(x,t)+Nu(x,t)=f(x,t),
where u is the unknown function, L and N are linear and nonlinear operators, respectively, and f is the source term. The correction functional for (7) is given as follows:
(8)un+1(x,t)=un(x,t)+∫0xλ[Lun(ξ,t)+Nuˇn(ξ,t)-f(ξ,t)]dξ,
where λ is a general Lagrange multiplier that can be identified optimally via the variation theory. The subscript n indicates the nth approximation and uˇn is considered as a restricted variation δuˇn=0.
4. Laplace Variational Iteration Method (LVIM)
Consider the following general multiterms fractional telegraph equation:(9)∂αu∂xα(x,t)=a1∂βu∂tβ(x,t)+a2∂γu∂tγ(x,t)+a3u(x,t)+f(x,t),
where 1<α, β≤2, 0<γ≤1, x,t≥0, u(0,t)=h(t), ux(0,t)=g(t), and a1, a2, a3 are constants.
The new approach of the Laplace variational iteration technique is based on the following steps.
Step 1.
Removing the fractional derivative of order α with respect to x for unknown function u(x,t) by using Laplace and inverse Laplace transforms.
Step 2.
Differentiating the results obtained in Step 1 with respect to x, then we get the value of the general Lagrange multiplier, for the correction functional (iterative formula) to equal one. The concept of the technique is illustrated in the following context.
By applying Laplace transform with respect to x, on both sides of (9), we get
(10)sαu-(s,t)-sα-1u(0,t)-sα-2ux(0,t)=ℒ[a1∂βu(x,t)∂tβ+a2∂γu(x,t)∂tγ+a3u(x,t)+f(x,t)]∴u-(s,t)=1sh(t)+1s2g(t)+1sα[ℒf(x,t)]+1sαℒ×[a1∂βu(x,t)∂tβ+a2∂γu(x,t)∂tγ+a3u(x,t)].
By taking the inverse Laplace transform to (10), we have
(11)u(x,t)=h(t)+xg(t)+ℒ-1[1sαℒf(x,t)]+ℒ-1×[1sαℒ[a1∂βu(x,t)∂tβ+a2∂γu(x,t)∂tγ+a3u(x,t)]].
Now the fractional derivative of order α with respect to x is removed, and the dependent variable u(x,t) in the left hand side of (11), became free of derivatives. Next step, we differentiate (10) with respect to x to get
(12)∂u(x,t)∂x=g(t)+∂∂x[ℒ-1[1sαℒf(x,t)]]+∂∂x×[ℒ-1[1sαℒ[a1∂βu(x,t)∂tβ+a2∂γu(x,t)∂tγ∂βun(x,t)∂tβ+a3u(x,t)]]].
The above step has been taken to enable us to construct the correction functional for (11) to be
(13)un+1=un(x,t)+∫0xλ{∂un(ξ,t)∂ξ-g(t)-∂∂ξ[ℒ-1[1sαℒf(x,t)]]{[[[a1∂βun(ξ,t)∂tβ]]]}-∂∂ξ[ℒ-1[1sαℒ[a1∂βun(ξ,t)∂tβ+a2∂γun(ξ,t)∂tγ∂βun(ξ,t)∂tβ+a3un(ξ,t)]]]}dξ.
The general Lagrange multiplier for (13) can be identified optimally via variation theory to get
(14)1+λ|ξ=x=0,λ`|ξ=x=0.
From (14), we obtain
(15)λ=-1.
Substituting λ=-1 into (13), we get the iterative formula for n=0,1,2,…, as follows:
(16)un+1=un(x,t)-∫0x{∂un(ξ,t)∂ξ-g(t)-∂∂ξ[ℒ-1[1sαℒf(x,t)]]{[[[a1∂βun(ξ,t)∂tβ]]]}-∂∂ξ[ℒ-1[1sαℒ[a1∂βun(ξ,t)∂tβ+a2∂γun(ξ,t)∂tγa1∂βun(ξ,t)∂tβ+a3un(ξ,t)]]]}dξ.
Start with the initial iteration
(17)u0(x,t)=u(0,t)+xux(0,t)=h(t)+xg(t).
The exact solution is given as a limit of the successive approximations un(x,t), n=0,1,2,…; in other words, u(x,t)=limn→∞un(x,t).
5. Numerical ExamplesExample 1.
Consider the following space-fractional homogenous telegraph equation:
(18)∂αu(x,t)∂xα=∂2u(x,t)∂t2+∂u(x,t)∂t+u(x,t),x,t≥0,1<α≤2,u(0,t)=e-t,ux(0,t)=e-t.
Solution 1.
Applying the Laplace transform with respect to x on both sides of (18), we get
(19)sαu-(s,t)-sα-1u(0,t)-sα-2ux(0,t)=ℒ[∂2u(x,t)∂t2+∂u(x,t)∂t+u(x,t)],u-(s,t)=1se-t+1s2e-t+1sαℒ[∂2u(x,t)∂t2+∂u(x,t)∂t+u(x,t)].
The inverse Laplace transform of (19) yields
(20)u(x,t)=e-t+xe-t+ℒ-1[1sαℒ[∂2u(x,t)∂t2+∂u(x,t)∂t+u(x,t)]].
Differentiating (20) with respect to x, we have
(21)∂u(x,t)∂x=e-t+∂∂xℒ-1[1sαℒ[∂2u(x,t)∂t2+∂u(x,t)∂t+u(x,t)]].
The correction functional for (21) with λ=-1 is given by
(22)un+1(x,t)=un(x,t)-∫0x[∂un(ξ,t)∂ξ-e-t-∂∂ξℒ-1×[1sαℒ[∂2un(ξ,t)∂t2+∂un(ξ,t)∂t[∂2un(ξ,t)∂t2+∂un(ξ,t)∂t+un(ξ,t)]]]dξ.
The initial iteration u0(x,t)=u(0,t)+xux(0,t)=e-t+xe-t; then, we have
(23)u1(x,t)=e-t+xe-t+xαΓ(α+1)e-t+xα+1Γ(α+2)e-t,u2(x,t)=e-t+xe-t+xαΓ(α+1)e-t+xα+1Γ(α+2)e-t+x2αΓ(2α+1)e-t+x2α+1Γ(2α+2)e-t,u3(x,t)=e-t+xe-t+xαΓ(α+1)e-t+xα+1Γ(α+2)e-t+x2αΓ(2α+1)e-t+x2α+1Γ(2α+2)e-t+x3αΓ(3α+1)e-t+x3α+1Γ(3α+2)e-t,⋮
Then the general term in successive approximation is given by
(24)un(x,t)=e-t∑k=0nxkα[1Γ(kα+1)+xΓ(kα+2)].
The solution in a closed form is given by
(25)u(x,t)=limn→∞un(x,t)=e-t[Eα,1(xα)+xEα,2(xα)].
Letting α=2, then
(26)u(x,t)=e-t[E2,1(x2)+xE2,2(x2)]=e-t[cosh(x)+xsinh(x)x]=ex-t.
The solution is the same as that obtained by Wazwaz [26]. The solution surface of this example is graphically presented in Figure 1 for various fractional orders of α.
The surface plot of u(x,t) solution of Example 1: (a) α=1.25, (b) α=1.5, (c) α=1.75, and (d) α=2.
Example 2.
Consider the following space-time fractional nonhomogenous telegraph equation:
(27)∂αu(x,t)∂xα=∂2u(x,t)∂t2+∂γu(x,t)∂tγ+u(x,t)-x2-t+1,1<α≤2,0<γ≤1,x,t≥0,u(0,t)=t,ux(0,t)=0.
Solution 2.
Taking Laplace transform with respect to x to (27), we get
(28)sαu-(s,t)-sα-1u(0,t)-sα-2ux(0,t)=ℒ[∂2u(x,t)∂t2+∂γu(x,t)∂tγ+u(x,t)]-2s3-ts+1s,u-(s,t)=1st-2sα+3-tsα+1+1sα+1+1sαℒ[∂2u(x,t)∂t2+∂γu(x,t)∂tγ+u(x,t)].
The inverse Laplace transform of (28) is given by
(29)u(x,t)=t-2Γ(α+3)xα+2-tΓ(α+1)xα+1Γ(α+1)xα+ℒ-1[1sαℒ[∂2u(x,t)∂t2+∂γu(x,t)∂tγ+u(x,t)]].
Differentiation of (29) with respect to x yields
(30)∂u(x,t)∂x=-2Γ(α+2)xα+1-tΓ(α)xα-1+1Γ(α)xα-1+∂∂x[ℒ-1[1sαℒ[∂2u(x,t)∂t2+∂γu(x,t)∂tγ+u(x,t)]]].
The correction functional for (30) with λ=-1 is given by
(31)un+1=un-∫0x[∂un(ξ,t)∂ξ+2Γ(α+2)ξα+1+tΓ(α)ξα-1-1Γ(α)ξα-1-∂∂ξ[ℒ-1[1sαℒ[∂2un(ξ,t)∂t2+∂γun(ξ,t)∂tγ[∂2un(ξ,t)∂t2]+un(ξ,t)]]]]dξ,u0(x,t)=u(0,t)+xux(0,t)=t,u1(x,t)=t-2xα+2Γ(α+3)+xαΓ(α+1)[1+t1-γΓ(2-γ)],u2(x,t)=t-2xα+2Γ(α+3)+xαΓ(α+1)[1+t1-γΓ(2-γ)]-2x2α+2Γ(2α+3)+x2αΓ(2α+1)[1+t1-γΓ(2-γ)],u3(x,t)=t-2xα+2Γ(α+3)+xαΓ(α+1)[1+t1-γΓ(2-γ)]-2x2α+2Γ(2α+3)+x2αΓ(2α+1)[1+t1-γΓ(2-γ)]-2x3α+2Γ(3α+3)+x3αΓ(3α+1)[1+t1-γΓ(2-γ)],⋮un(x,t)=t-2x2∑k=1n(xα)kΓ(kα+3)+[1+t1-γΓ(2-γ)]∑k=1n(xα)kΓ(kα+1).u(x,t)=limn→∞un(x,t)=t+x2-[1+t1-γΓ(2-γ)]-2x2∑n=0∞(xα)kΓ(kα+3)+[1+t1-γΓ(2-γ)]∑n=0∞(xα)kΓ(kα+1)=t+x2[1-2Eα,3(xα)]+[1+t1-γΓ(2-γ)][Eα,1(xα)-1].
For α=2 and γ=1, then
(32)u(x,t)=t+x2-2-2x2E2,3(x2)+2E2,1(x2)=t+x2.
The solution surface of this example is graphically presented in Figure 2 for fixed γ and various fractional orders of α.
The surface plot of u(x,t) solution of Example 2 for fixed γ=0.5: (a) α=1.25, (b) α=1.5, (c) α=1.75, and (d) α=2.
Example 3.
Consider the following space-time fractional nonhomogenous telegraph equation:
(33)∂αu(x,t)∂xα=∂βu(x,t)∂tβ+∂γu(x,t)∂tγ+u(x,t)-x2-t2-2t,1<α,β≤2,23<γ≤1,x,t≥0,u(0,t)=t2,ux(0,t)=0,2<β+γ≤3.
Solution 3.
Taking Laplace transform with respect to x to (33), we get
(34)sαu-(s,t)-sα-1u(0,t)-sα-2ux(0,t)=ℒ[∂βu(x,t)∂tβ+∂γu(x,t)∂tγ+u(x,t)]-2s3-t2s-2ts,u-(s,t)=1st2+1sαℒ[∂βu(x,t)∂tβ+∂γu(x,t)∂tγ+u(x,t)]-2sα+3-t2sα+1-2tsα+1.
The inverse Laplace transform of (34) is given by
(35)u(x,t)=t2-2xα+2Γ(α+3)-t2xαΓ(α+1)-2txαΓ(α+1)+ℒ-1[1sαℒ[∂βu(x,t)∂tβ+∂γu(x,t)∂tγ+u(x,t)]].
Differentiation of (35) with respect to x yields
(36)∂u(x,t)∂x=-2xα+1Γ(α+2)-t2xα-1Γ(α)-2txα-1Γ(α)+∂∂x[ℒ-1[1sαℒ[∂βu(x,t)∂tβ+∂γu(x,t)∂tγ+u(x,t)]]].
The correction functional for (36) with λ=-1 is given by
(37)un+1=un-∫0x[∂un(ξ,t)∂ξ+2ξα+1Γ(α+2)+t2ξα-1Γ(α)+2tξα-1Γ(α)-∂∂ξ[ℒ-1[1sαℒ[∂βun(ξ,t)∂tβ+∂γun(ξ,t)∂tγ[∂βun(ξ,t)∂tβ+un(ξ,t)]]]]dξ,u0(x,t)=u(0,t)+xux(0,t)=t2,u1(x,t)=t2-2xα+2Γ(α+3)+xαΓ(α+1)[2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t],u2(x,t)=t2-2xα+2Γ(α+3)+xαΓ(α+1)[2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t]-2x2α+2Γ(2α+3)+x2αΓ(2α+1)×[2t2-2γΓ(3-2γ)-2t1-γΓ(2-γ)+2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t],u3(x,t)=t2-2xα+2Γ(α+3)+xαΓ(α+1)[2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t]-2x2α+2Γ(2α+3)+x2αΓ(2α+1)×[2t2-2γΓ(3-2γ)-2t1-γΓ(2-γ)+2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t]-2x3α+2Γ(3α+3)+x3αΓ(3α+1)×[4t2-2γΓ(3-2γ)-4t1-γΓ(2-γ)+2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t],⋮un(x,t)=t2-2x2∑k=1n(xα)kΓ(kα+3)+[2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t]∑k=1n(xα)kΓ(kα+1)+[2t2-2γΓ(3-2γ)-2t1-γΓ(2-γ)]×∑k=2n(k-1)(xα)kΓ(kα+1),n≥2,u(x,t)=limn→∞un(x,t)=t2+x2-[2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t]-2x2∑n=0∞(xα)kΓ(kα+3)+[2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t]∑n=0∞(xα)kΓ(kα+1)+[2t2-2γΓ(3-2γ)-2t1-γΓ(2-γ)]×∑k=2∞(k-1)(xα)kΓ(kα+1),n≥2=t2+x2[1-2Eα,3(xα)]+[2t2-βΓ(3-β)+2t2-γΓ(3-γ)-2t][Eα,1(xα)-1]+[2t2-2γΓ(3-2γ)-2t1-γΓ(2-γ)]×[1+Eα,12(xα)-2Eα,1(xα)].
For α=β=2 and γ=1, we get the standard equation with the solution
(38)u(x,t)=t2+x2-2-2x2E2,3(x2)+2E2,1(x2)+0=t2+x2.
The solution surface of this example is graphically presented in Figure 3 for fixed γ and β and various fractional orders of α.
The surface plot of u(x,t) solution of Example 3 for fixed β=1.25 and γ=0.75: (a) α=1.25, (b) α=1.5, (c) α=1.75, and (d) α=2.
6. Conclusion
In this paper, a combined form of Laplace transform and variational iteration method is presented to handle space-time fractional telegraph equations in a half-space domain. The space and time derivatives are considered in the Caputo sense. Certain techniques are used to overcome the complexity of identifying the general Lagrange multiplier. The solutions are obtained in series form that rapidly converges in a closed exact formula with simply computable terms. The calculations are simple and straightforward. The method was tested on three examples on different situations. The technique is powerful, reliable, and efficient. This technique can be extended to solve various linear and nonlinear fractional problems in applied science.
SabatierJ.AgrawalO. P.MachadoJ. A. T.2007New York, NY, USASpringer10.1007/978-1-4020-6042-7MR2432163DebnathL.Recent applications of fractional calculus to science and engineering2003543413344210.1155/S0161171203301486MR2025566ZBL1036.26004BhattiM.Fractional Schrödinger wave equation and fractional uncertainty principle2007219943950MR2371326ZBL1146.35082GanjiZ. Z.GanjiD. D.RostamiyanY.Solitary wave solutions for a time-fraction generalized Hirota-Satsuma coupled KdV equation by an analytical technique20093373107311310.1016/j.apm.2008.10.034MR2502734ZBL1205.35251GutirrezR. E.RosarioJ. M.MachadoJ. T.Fractional order calculus: basic concepts and engineering applications201020101937585810.1155/2010/375858SevimlicanA.An approximation to solution of space and time fractional telegraph equations by He's variational iteration method201020101029063110.1155/2010/290631MR2610512ZBL1191.65137CamposD.MéndezV.Different microscopic interpretations of the reaction-telegrapher equation20094271307500310.1088/1751-8113/42/7/075003MR2525453ZBL1162.82012GargM.SharmaA.Solution of space-time fractional telegraph equation by Adomian decomposition method20112117MR2827220GargM.ManoharP.KallaS. L.Generalized differential transform method to space-time fractional telegraph equation201120119548982MR2843508ZBL1234.3529910.1155/2011/548982ArbabA. I.Derivation of Dirac, Klein-Gordon, Schrodinger, diffusion and quantum heat transport equations from the universal wave equation20109244000110.1209/0295-5075/92/40001ChenJ.LiuF.AnhV.Analytical solution for the time-fractional telegraph equation by the method of separating variables200833821364137710.1016/j.jmaa.2007.06.023MR2386503ZBL1138.35373AnsariA.Fractional exponential operators and time-fractional telegraph equation20122012article 12510.1186/1687-2770-2012-125MR3016677HuangF.Analytical solution for the time-fractional telegraph equation20092009989015810.1155/2009/890158MR2565356ZBL1190.35224DasS.VishalK.GuptaP. K.YildirimA.An approximate analytical solution of time-fractional telegraph equation2011217187405741110.1016/j.amc.2011.02.030MR2784584ZBL1216.65135JiangW.LinY.Representation of exact solution for the time-fractional telegraph equation in the reproducing kernel space20111693639364510.1016/j.cnsns.2010.12.019MR2787810ZBL1223.35112KhanY.DiblikJ.FarazN.SmardaZ.An efficient new perturbation Laplace method for space-time fractional telegraph equations20122012article 20410.1186/1687-1847-2012-204HeJ. H.An aproximation to solution of space and time fractional telegraph equations by the variational iteration method201220122394212MR2983802ZBL1264.6517210.1155/2012/394212AbassyT. A.El-TawilM. A.El-ZoheiryH.Exact solutions of some nonlinear partial differential equations using the variational iteration method linked with Laplace transforms and the Padé technique2007547-894095410.1016/j.camwa.2006.12.067MR2395632ZBL1141.65382HammouchZ.MekkaouiT.A Laplace-variational iteration method for solving the homogeneous Smoluchowski coagulation equation2012618879886MR2891394ZBL1252.65218ArifeA. S.YildirimA.New modified variational iteration transform method (MVITM) for solving eighth-order boundary value problems in one step2011131021862190PodlubnyI.1999198New York, NY, USAAcademic PressMathematics in Science and EngineeringMR1658022DebnathL.BhattaD.20072ndTaylor & Francis Group,LLCMR2253985HauboldH. J.MathaiA. M.SaxenaR. K.Mittag-Leffler functions and their applications201120115129862810.1155/2011/298628MR2800586ZBL1218.33021Berberan-SantosM. N.Properties of the Mittag-Leffler relaxation function200538462963510.1007/s10910-005-6909-zMR2210207ZBL1101.33015HeJ. H.Variational iteration method—a kind of nonlinear analytical technique: some examples199934469970810.1016/S0020-7462(98)00048-1WazwazA. M.2009Heidelberg, GermanyHigher Education Press, Beijing, China; SpringerNonlinear Physical Science10.1007/978-3-642-00251-9MR2548288HeJ. H.WuX. H.Variational iteration method: new development and applications2007547-888189410.1016/j.camwa.2006.12.083MR2395625ZBL1141.65372KhuriS. A.SayfyA.A Laplace variational iteration strategy for the solution of differential equations201225122298230510.1016/j.aml.2012.06.020MR2967833ZBL1252.65128