IJDE International Journal of Differential Equations 1687-9651 1687-9643 Hindawi Publishing Corporation 472030 10.1155/2012/472030 472030 Research Article Solving the Fractional Rosenau-Hyman Equation via Variational Iteration Method and Homotopy Perturbation Method Yulita Molliq R. 1 Noorani M. S. M. 2 Momani Shaher 1 Department of Mathematics Faculty of Mathematics and Natural Science Universitas Negeri Medan UNIMED 20221 Medan, Sumatera Utara Indonesia unimed.ac.id 2 School of Mathematical Sciences Faculty of Science and Technology Universiti Kebangsaan Malaysia UKM 43600 Bangi, Selangor Malaysia ukm.my 2012 6 12 2012 2012 31 05 2012 08 11 2012 2012 Copyright © 2012 R. Yulita Molliq and M. S. M. Noorani. 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.

In this study, fractional Rosenau-Hynam equations is considered. We implement relatively new analytical techniques, the variational iteration method and the homotopy perturbation method, for solving this equation. The fractional derivatives are described in the Caputo sense. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for fractional Rosenau-Hynam equations. In these schemes, the solution takes the form of a convergent series with easily computable components. The present methods perform extremely well in terms of efficiency and simplicity.

1. Introduction

Recent advances of fractional differential equations are stimulated by new examples of applications in fluid mechanics, viscoelasticity, mathematical biology, electrochemistry, and physics. For example, the nonlinear oscillation of earthquake can be modeled with fractional derivatives , and the fluid-dynamic traffic model with fractional derivatives  can eliminate the deficiency arising from the assumption of continuum traffic flow. Based on experimental data fractional partial differential equations for seepage flow in porous media are suggested in , and differential equations with fractional order have recently proved to be valuable tools to the modeling of many physical phenomena . Fractional partial differential equations also have studied and successfully solved such as the space-time fractional diffusion-wave equation [5–7], the fractional advection-dispersion equation [8, 9], the fractional KdV equation , and the linear inhomogeneous fractional partial differential equations .

Most nonlinear differential equations are usually arising from mathematical modeling of many physical systems. In most cases, it is very difficult to achieve analytic solutions of these equations. Perturbation techniques are widely used in science and engineering to handle nonlinear problems and do great contribution to help us understand many nonlinear phenomena. However, perturbation techniques are based on the existence of small/large parameter. Therefore, these techniques are not valid for strongly nonlinear problems.

The homotopy perturbation method (HPM) is the new approach for finding the approximate analytical solution of linear and nonlinear problems. The method was first proposed by He [12, 13] and was successfully applied to solve nonlinear wave equation by He . The convergence of Homotopy perturbation series to the exact solution is considered in . Similarly, applying the variational iteration method, created by He [3, 18, 19], consists in constructing the appropriate correction functional connected with the considered equation. The correction functional contains a Lagrange multiplier, the determination of which leads to a recurrence formula. Convergence of the VIM method is discussed by Tatari and Dehghan in . Both of the methods examined have found application in determining the approximate solutions of different technical problems . Adaptation of the VIM method for solving fractional heat-wave-like equation and fractional Zakharov-Kuznetsov equation were discussed by Yulita and colleagues in [22, 23]. Whereas, Chun  obtained the numerical solution of heat conduction problem by VIM. Recently, the application of the VIM method for solving kuramoto and Sivashinsky equations was presented by Porshokouhi and Ghanbari in . For the application of HPM, this method used for solving fractional vibration equation  and partial differential equations of fractional order in finite domains .

In the present paper, VIM and HPM will be applied for solving fractional Rosenau-Hay nam equation which written as (1.1)Dtαu=uDxxx(u)+uDx(u)+3Dx(u)Dxx(u),t>0, subject to the initial condition (1.2)u(x,0)=-83ccos2(x4), where u=u(x,t), α is a parameter describing the order of the fractional derivative (0<α1), t is the time, and x is the spatial coordinate. Fractional RH equation when α=1.0 has appeared in the study of the formation of patterns in liquid drops .

2. Basic Definitions

Fractional calculus unifies and generalizes the notions of integer-order differentiation and n-fold integration [4, 29]. We give some basic definitions and properties of fractional calculus theory which will be used in this paper:

Definition 2.1.

A real function f(x),x>0 is said to be in the space 𝒞μ,  μ  R if there exists a real number p(>μ), such that f(x)=xpf1(x), where f1(x)𝒞[0,), and it is said to be in the space 𝒞μm if and only if f(m)𝒞m,  mN.

The Riemann-Liouville fractional integral operator is defined as follows.

Definition 2.2.

The Riemann-Liouville fractional integral operator of order α0, of a function fCμ,  μ-1, is defined as (2.1)Jαf(x)=1Γ(α)0x(x-t)α-1f(t)dt,α>0,x>0,J0f(x)=f(x)

In this paper only real and positive values of α will be considered.

Properties of the operator Jα can be found in  and we mention only the following: for fCμ, μ-1,  α,β0, and γ-1:

JαJβf(x)=Jα+βf(x),

JαJβf(x)=JβJαf(x),

Jαxγ=(Γ(γ+1)/Γ(α+γ+1))xα+γ.

The Reimann-Liouville derivative has certain disadvantages when trying to model real-world phenomena with FDEs. Therefore, we will introduce a modified fractional differential operator D*α proposed by Caputo in his work on the theory of viscoelasticity .

Definition 2.3.

The fractional derivative of f(x) in Caputo sense is defined as (2.2)D*αf(x)=Jm-αD*mf(x)=1Γ(m)0x(x-s)m-α-1f(m)(s)ds,for  m-1<αm,m𝒩,x>0,fC-1m In addition, we also need the following property.

Lemma 2.4.

If m-1<αm,  m𝒩 and fCμm,  μ-1, then (2.3)D*αJαf(x)=f(x),JαD*αf(x)=f(x)-i=0m-1f(i)(0+)xii!,x>0.The Caputo differential derivative is considered here because the initial and boundary conditions can be included in the formulation of the problems . The fractional derivative is taken in the Caputo sense as follows.

Definition 2.5.

For m to be the smallest integer that exceeds α, the Caputo fractional derivative operator of order α>0 is defined as (2.4)Dtαu(x,t)={1Γ(m-α)0t(t-s)m-α-1(mu(x,s)/sm)ds,for  m-1<αm,mu(x,t)tm,for  α=mN.For more information on the mathematical properties of fractional derivatives and integrals, one can consult [4, 29].

3. Basic Idea of He’s Variational Iteration Method

To clarify the basic ideas of VIM, we consider the following differential equation: (3.1)D*tα(u)=f(u,ux,uxx)+g(x,t),m-1<α<m, where u=u(x,t),  D*tα=α/tα is the Caputo fractional derivative of order α, mN, f is a nonlinear function, and g is the source function. According to VIM, we can write down a correction functional as follows: (3.2)un+1(x,t)=un(x,t)+0tλ(s)[msm(un)-f(u~n,u~n,x,u~n,xx)-g(x,s)]ds, where λ(s) is a general Lagrangian multiplier which can be optimally identified via the variational theory . u~n is considered as a restricted variation , that is, δu~n=0 and the subscript n indicates the nth approximation. We have (3.3)δun+1=δun+δ0tλ(s)[msm(un)-f(u~n,u~n,x,u~n,xx)-g(x,s)]ds, where u~n is considered as restricted variations, that is, δu~n=0. For m=1, we have (3.4)δun+1=δun+δ0tλ(s)[sun]ds,δun+1=[1-λ(t)]δun+δλ(t)tun+01δλ′′(s)unds. Thus, we obtain the following stationary conditions: (3.5)δun(t):[1-λ(t)]|s=t=0,δun(t):λ(t)|s=t=0,δun(s):λ′′(t)=0. Solving this system of equations yields (3.6)λ(s)=-1. Furthermore, substituting (3.6) to (3.2), the iteration formula of VIM can be written as follows: (3.7)un+1=un-0t[msm(un)-f(un,un,x,un,xx)-g(x,s)]ds. In this case, we begin with the initial approximation: (3.8)u0(x,t)=h(x). The correction functional (3.8) will give several approximations, and therefore the exact solution is obtained as (3.9)u(x,t)=limnun(x,t).

4. Basic Idea of Homotopy Perturbation Method

Consider the following nonlinear differential equation: (4.1)A(u)-f(r)=0, with boundary conditions (4.2)B(u,un)=0,rΓ, where A is a general differential operator, B is a boundary operator, f(r) is a known analytic function, and Γ is the boundary of the domain Ω.

The operator A can, generally speaking, be divided into two parts, L and N, where L is linear and N is nonlinear, therefore (4.1) can be written as (4.3)L(u)+N(u)-f(r)=0. By using homotopy technique, one can construct a homotopy v(r,p):Ω×[0,1] which satisfies (4.4)H(v,p)=(1-p)[L(v)-L(u0)]+p[A(v)-f(r)]=0, or (4.5)H(v,p)=L(v)-L(u0)+pL(u0)+p[N(v)-f(r)]=0, where p[0,1] is an embedding parameter, and u0 is the initial approximation of (4.1) which satisfies the boundary conditions. Clearly, we have (4.6)H(v,0)=L(v)-L(u0)=0, or (4.7)H(v,1)=A(v)-f(r)=0, the changing process of p from zero to unity is just that of v(r,p) changing from u0(r) to u(r). This is called deformation, and also, L(v)-L(u0) and A(v)-f(r) are called homotopic in topology. If the embedding parameter p(0p1) is considered as a small parameter, applying the classical perturbation technique, we can assume that the solution of (4.3) and (4.4) can be given as a power series in p, that is, (4.8)v=v0+pv1+p2v2+, and setting p=1 results in the approximate solution of (4.1) as (4.9)u=limp1v=v0+v1+v2+.

5. Application of HPM and VIM Methods

We consider the application of VIM to fractional Rosenau-Hynam equations which is rewritten as follows: (5.1)Dtαu=uDxxx(u)+uDx(u)+3Dx(u)Dxx(u), where u=u(x,t) with the initial conditions of (5.2)u(x,0)=-83ccos2(x4).

5.1. VIM Implement for Fractional Rosenau-Hynam Equation

According to the formula (3.8), the iteration formula for (4.9) is given by (5.3)un+1=un-0t[Dsα(un)-unDxxx(un)-unDx(un)-3Dx(un)Dxx(un)]ds. The iteration starts with an initial approximation which is initial condition in (5.2). Furthermore, using the iteration formula in (5.3), we can directly obtain other components as (5.4)u1(x,t)=-23c2sin(x2)t-43c[1+cos(x2)],u2(x,t)=-43c2sin(x2)+[16c2+43c](1+cos(x2))-16c3t2+23c2sin(x2)t(2-α)Γ(3-α),u3(x,t)=76c3t3cos(x2)-9636ccos2(x4)-4836c2tsin(x2)-23cos(x2)t3-αΓ(4-α)+(6-2α)c2sin(x2)t2-αΓ(4-α)-23c2sin(x2)t3-2αΓ(4-2α), and so on.

The exact solution of this equation is given by  (5.5)u(x,t)=-83ccos2(14(x-ct)),|x-ct|2π, where c is arbitrary constant .

6. Analysis of the Homotopy Perturbation Method (HPM)

Now applying the classical perturbation technique for solving (5.1) with initial condition in (5.2). To solve (5.1) by the homotopy perturbation method, we construct the following homotopy: (6.1)(1-p)Dtαu(x,t)+p[Dtαu(x,t)-(u(x,t)Dxxxu(x,t)+u(x,t)Dxu(x,t)(1-p)Dtαu(x,t)+pDtαu(x,t)-+3Dxu(x,t)Dxxu(x,t))Dtα]=0, or (6.2)Dtαu(x,t)=p[u(x,t)Dxxxu(x,t)+u(x,t)Dxu(x,t)+3Dxu(x,t)Dxxu(x,t)], where p[0,1] is an embedding parameter. If p=0, then (6.2) becomes a linear equation, (6.3)Dtαu(x,t)=0. And when p=1, then (6.2) turns out to be (5.1). Assume the solution of (6.2) to be in the form (6.4)u(x,t)=v0+pv1+p2v2+p3v3+. Substituting (6.4) into (6.2) and equating the terms with identical powers of p, we obtain the following set of linear differential equations: (6.5)p0:Dtαu0=0,(6.6)p1:Dtαu1=v0v0x+3v0x2v0x2+v03v0x3,(6.7)p2:Dtαu2=v1v0x+v0v1x+3v1x2v0x2+3v0x2v1x2+v13v0x3+v03v1x3, and so on. Equations (6.6) and (6.7) can be solved by applying the operator Jα, which is the inverse of the operator Dα and then by simple computation, Thus, the solution reads as follows (6.8)v1(x,t)=-2c2tαsin(x/4)3Γ(1+α),v2(x,t)=-c33πt2αcos(x2)Γ(-2α)sin1πα,v3(x,t)=-c4πt3α csc παsin(x/2)6Γ(1-α)Γ(α)Γ(1+3α),v4(x,t)=-c5πt4αcos(x/2) csc απ12Γ(1-α)Γ(α)Γ(1+4α), and so on.

In this manner, the rest of the components of the homotopy perturbation series can be obtained. Finally, we approximate the analytical solutions of u(x,t) by the truncated series (6.9)u(x,t)=limNu5(x,t), where u5(x,t)=n=04vk(x,t).

7. Numerical Results and Discussion

Tables 1 and 2 show the approximate solutions for (5.1) obtained for different values of α using the decomposition method and the variational iteration method in different values of c, that is,   c=1 and c=0.5, respectively. Tables 3 and 4 show the absolute error of (5.1) when α=1.0 in different value of c, that is,   c=1.0 and c=0.5, respectively. Figures 1 and 2 show the approximate solutions for (5.1) in different values of c using the fifth iterates of VIM, the fifth terms of HPM when α=1, and exact solution, respectively. From Tables 3 and 4 show that the approximate solution using the VIM is more accurate than the approximate solution obtained using the HPM. It is to be noted that only the fifth iterates of the variational iteration solution and only fifth terms of the homotopy perturbation series were used in evaluating.

The approximate solution of fractional Rosenau-Hynam equation using fifth iterates of VIM and fifth terms of HPM in different values of α, that is, α=0.7,0.8, and 0.9 and α=1.0 when c=0.5.

x    t VIM HPM VIM HPM VIM HPM VIM HPM Exact
α = 0.7 α = 0.8 α = 0.9 α = 1.0
π / 4 0.2 −1.3027037 −1.3020263 −1.2995666 −1.2994063 −1.2968978 −1.2968865 −1.1139217 −1.2945674 −1.2945674
0.6 −1.3163161 −1.3160219 −1.3159022 −1.3158706 −1.3150924 −1.3150940 −1.3137953 −1.3137952 −1.3137952
1.0 −1.3194666 −1.3204884 −1.3228476 −1.3231121 −1.3251352 −1.3251648 −1.3265574 −1.3265572 −1.3265572

π / 2 0.2 −1.1781765 −1.1775095 −1.1715349 −1.1713677 −1.1658775 −1.1658583 1.1610425 −1.1610425 −1.1610425
0.6 −1.2168471 −1.2161259 −1.2124978 −1.2123777 −1.2080227 −1.2079404 −1.2032236 −1.2032236 −1.2032236
1.0 −1.2404292 −1.2410937 −1.2415934 −1.2417625 −1.2413894 −1.2412100 −1.2400439 −1.2400437 −1.2400437

3 π / 4 0.2 −0.4677769 −0.9752216 −0.4579952 −0.9664930 −0.4496547 −1.1658655 −0.4426464 −0.2194117 −0.4426464
0.6 −0.5359648 −1.0325797 −0.5255989 −1.0258054 −0.5158600 −1.0185228 −0.5064510 −0.2710013 −0.5064510
1.0 −0.5899819 −1.0742478 −0.5848618 −1.0728597 −0.5788060 −1.0701462 −0.5718564 −0.3265443 −0.5718565

π 0.2 −0.7263173 −0.7259589 −0.7160797 −0.7159725 −0.7073528 −0.7073345 −0.6999861 −0.6999861 −0.6999861
0.6 −0.7945239 −0.7933266 −0.7847888 −0.7845575 −0.7754812 −0.7753489 −0.7662921 −0.7662921 −0.7662921
1.0 −0.8456349 −0.8453513 −0.8421984 −0.8421178 −0.8374790 −0.8371566 −0.8316027 −0.8316026 −0.8316026

The approximate solution of fractional Rosenau-Hynam using fifth iterate of VIM and fifth terms of HPM in different values of α, that is, α=0.7,0.8, and 0.9 and α=1.0 when c=1.0.

x    t VIM HPM VIM HPM VIM HPM VIM HPM Exact
α = 0.7 α = 0.8 α = 0.9 α = 1.0
π / 4 0.2 −2.6345127 −2.6293263 −2.6251671 −2.6239643 −2.6172579 −2.6171710 −2.6099581 −2.6099581 −2.6099581
0.6 −2.6345745 −2.6342711 −2.6474590 −2.6477297 −2.6563851 −2.6450822 −2.6609433 −2.6609420 −2.6609420
1.0 −2.5722042 −2.5854045 −2.6080979 −2.6122085 −2.6365766 −2.6372049 −2.6590258 −2.6589984 −2.6589984

π / 2 0.2 −2.4277172 −2.4226797 −2.4036916 −2.4024670 −2.3832715 −2.4054513 −2.3655561 −2.3655561 −2.3655561
0.6 −2.5350335 −2.5308256 −2.5303337 −2.5295597 −2.5232013 −2.5409052 −2.5126533 −2.5126523 −2.5126523
1.0 −2.5563768 −2.5641211 −2.5831444 −2.5854712 −2.6015950 −2.6030469 −2.6127547 −2.6127328 −2.6127328

3 π / 4 0.2 −0.5752642 −2.0501900 −0.5416828 −2.0182039 −0.5130518 −1.9893536 −0.4893583 −0.4893583 −1.9640076
0.6 −0.8341872 −2.5308256 −0.7891721 −2.2292752 −0.7486279 −2.1055837 −0.7112517 −0.7112526 −2.1848220
1.0 −1.0486958 −2.3554614 −1.0195238 −2.3681073 −0.9896491 −2.3733581 −0.9579284 −0.9579473 −2.3716903

π 0.2 −1.5711438 −2.4226797 −1.5304096 −1.5296754 −1.4957251 −1.5301673 −1.4664446 −1.4664446 −1.4664446
0.6 −1.8301141 −2.5308256 −1.7949603 −1.7925916 −1.7612479 −1.7971174 −1.7273603 −1.7273603 −1.7273603
1.0 −1.9968865 −2.5641211 −1.9951864 −1.9932084 −1.9871083 −2.0042754 −1.9725695 −1.9725674 −1.9725674

Errors of the approximate solution of fractional Rosenau-Hynam equation using using fifth term of HPM and fifth iterate of VIM when α=1.0 and c=0.5.

t    x = π / 2 x = π x = 3 π / 2
VIM HPM VIM HPM VIM HPM
0.1 1.0000 E - 11 1.0000 E - 11 5.0000 E - 10 1.0000 E - 11 2.0000 E - 10 1.0000 E - 11
0.2 1.0000 E - 11 1.2173 E - 09 5.0000 E - 10 1.7360 E - 09 3.0000 E - 10 1.2378 E - 09
0.3 1.0000 E - 11 9.2044 E - 09 5.0000 E - 10 1.3182 E - 08 3.0000 E - 10 9.4375 E - 09
0.4 1.0000 E - 10 3.8620 E - 08 1.0000 E - 10 5.5542 E - 08 9.0000 E - 10 3.9929 E - 08
0.5 3.0000 E - 10 1.1734 E - 07 4.0000 E - 10 1.6948 E - 07 1.9000 E - 09 1.2233 E - 07
0.6 1.0000 E - 09 2.9070 E - 07 7.0000 E - 10 4.2165 E - 07 6.1000 E - 09 3.0561 E - 07
0.7 2.6000 E - 09 6.2550 E - 07 1.2000 E - 09 9.1117 E - 07 1.8300 E - 08 6.6309 E - 07
0.8 5.4000 E - 09 1.2140 E - 06 2.1000 E - 09 1.7761 E - 06 3.9300 E - 08 1.2978 E - 06
0.9 1.1110 E - 07 2.1777 E - 06 4.0000 E - 09 3.1998 E - 06 8.2200 E - 08 2.3474 E - 06
1.0 2.1110 E - 07 3.6709 E - 06 8.6000 E - 09 5.4173 E - 06 1.5230 E - 07 3.9903 E - 06

Errors of the approximate solution of fractional Rosenau-Hynam equation using fifth term of HPM and fifth iterate of VIM when α=1.0 and c=0.5.

t    x = π / 2 x = π x = 3 π / 2
VIM HPM VIM HPM VIM ADM
0.1 1.0000 E - 11 2.4346 E - 09 5.0000 E - 10 3.4720 E - 09 2.0000 E - 10 2.4756 E - 09
0.2 1.0000 E - 11 7.7240 E - 08 5.0000 E - 10 1.1109 E - 07 3.0000 E - 10 7.9858 E - 08
0.3 1.0000 E - 11 5.8140 E - 07 5.0000 E - 10 8.4330 E - 07 3.0000 E - 10 6.1121 E - 07
0.4 1.0000 E - 10 2.4280 E - 06 1.0000 E - 10 3.5522 E - 06 9.0000 E - 10 2.5955 E - 06
0.5 3.0000 E - 10 7.3419 E - 06 4.0000 E - 10 1.0835 E - 05 1.9000 E - 09 7.9805 E - 06
0.6 1.0000 E - 09 1.8100 E - 05 7.0000 E - 10 2.6942 E - 05 6.1000 E - 09 2.0004 E - 05
0.7 2.6000 E - 09 3.8743 E - 05 1.2000 E - 09 5.8188 E - 05 1.8300 E - 08 4.3547 E - 05
0.8 5.4000 E - 09 7.4800 E - 05 2.1000 E - 09 1.1335 E - 04 3.9300 E - 08 8.5496 E - 05
0.9 1.1110 E - 07 1.3340 E - 04 4.0000 E - 09 2.0405 E - 04 8.2200 E - 08 1.5512 E - 04
1.0 2.1110 E - 07 2.2370 E - 04 8.6000 E - 09 3.4516 E - 04 1.5230 E - 07 2.6444 E - 04

The approximate solution of fractional Rosenau-Hynam equation when α=1.0 and c=0.5 using such method: (a) HPM, (b) VIM and (c) Exact.

The approximate solution of fractional Rosenau-Hynam equation when α=1.0 and c=1.0 using such method: (a) HPM, (b) VIM, and (c) Exact.

8. Conclusions

The fundamental goal of this work has been to construct an approximate solution of nonlinear partial differential equations of fractional order. For computations and plots, the Mathematica and Maple packages were used. The goal has been achieved by using the variational iteration method (VIM) and the homotopy perturbation method (HPM). The methods were used in a direct way without using linearization or restrictive assumptions. There are four important points that were gotten. First, the VIM and the HPM provide the solutions in terms of convergent series with easily computable components. Second, the approximate solution in (5.1) using the VIM converges faster than the approximate solution using the HPM. Third, the variational iteration method handles nonlinear equations without any need for the so-called He’s polynomials. Finally, the recent appearance of fractional differential equations as models in some fields of applied mathematics makes it necessary to investigate methods of solution for such equations (analytical and numerical) and we hope that this work is a step in this direction.

Acknowledgment

The financial support received from UKM Grant ERGS/1/2011/STG/UKM/01/13 is gratefully acknowledged.

He J. H. Semi-inverse method of establishing generalized variational principles for fluid mechanics with emphasis on turbomachinery aerodynamics International Journal of Turbo and Jet Engines 1997 14 1 23 28 2-s2.0-0031441505 He J. H. Some applications of nonlinear fractional differential equations and their approximations Bulletin of Science, Technology & Society 1999 15 2 86 90 He J. H. Approximate analytical solution for seepage flow with fractional derivatives in porous media Computer Methods in Applied Mechanics and Engineering 1998 167 1-2 57 68 10.1016/S0045-7825(98)00108-X 1665221 ZBL0942.76077 Podlubny I. Fractional Differential Equations 1999 198 San Diego, Calif, USA Academic Press Mathematics in Science and Engineering 1658022 Al-Khaled K. Momani S. An approximate solution for a fractional diffusion-wave equation using the decomposition method Applied Mathematics and Computation 2005 165 2 473 483 10.1016/j.amc.2004.06.026 2137805 ZBL1071.65135 Mainardi F. Luchko Y. Pagnini G. The fundamental solution of the space-time fractional diffusion equation Fractional Calculus & Applied Analysis 2001 4 2 153 192 1829592 ZBL1054.35156 Hanyga A. Multidimensional solutions of time-fractional diffusion-wave equations Proceedings of the Royal Society of London, Series A 2002 458 2020 933 957 10.1098/rspa.2001.0904 1898095 ZBL1153.35347 Huang F. Liu F. The time fractional diffusion equation and the advection-dispersion equation The Australian & New Zealand Industrial and Applied Mathematics Journal 2005 46 3 317 330 10.1017/S1446181100008282 2124926 ZBL1072.35218 Huang F. Liu F. The fundamental solution of the space-time fractional advection-dispersion equation Journal of Applied Mathematics & Computing 2005 18 1-2 339 350 10.1007/BF02936577 2137713 ZBL1086.35003 Momani S. An explicit and numerical solutions of the fractional KdV equation Mathematics and Computers in Simulation 2005 70 2 110 118 10.1016/j.matcom.2005.05.001 2163618 ZBL1119.65394 Debnath L. Bhatta D. D. Solutions to few linear fractional inhomogeneous partial differential equations in fluid mechanics Fractional Calculus & Applied Analysis 2004 7 1 21 36 2077398 ZBL1076.35096 He J. H. Homotopy perturbation technique Computer Methods in Applied Mechanics and Engineering 1999 178 3-4 257 262 10.1016/S0045-7825(99)00018-3 1711041 ZBL0956.70017 He J. H. A coupling method of a homotopy technique and a perturbation technique for non-linear problems International Journal of Non-Linear Mechanics 2000 35 1 37 43 10.1016/S0020-7462(98)00085-7 1723761 He J. H. Periodic solutions and bifurcations of delay-differential equations Physics Letters A 2005 347 4–6 228 230 10.1016/j.physleta.2005.08.014 2176126 ZBL1195.34116 He J. H. Application of homotopy perturbation method to nonlinear wave equations Chaos, Solitons and Fractals 2005 26 3 695 700 2-s2.0-18844426016 10.1016/j.chaos.2005.03.006 ZBL1072.35502 He J. H. Limit cycle and bifurcation of nonlinear problems Chaos, Solitons and Fractals 2005 26 3 827 833 2-s2.0-18844391045 10.1016/j.chaos.2005.03.007 ZBL1093.34520 He J. H. Non-perturbative methods for strongly nonlinear problems [Dissertation] 2006 Berlin, Germany de-Verlag im Internet GmbH He J. H. A new approach to nonlinear partial differential equations Communications in Nonlinear Science and Numerical Simulation 1997 2 4 230 235 2-s2.0-0007140364 10.1016/S1007-5704(97)90007-1 ZBL0923.35046 He J. H. Variational iteration method for delay differential equations Communications in Nonlinear Science and Numerical Simulation 1997 2 4 235 236 2-s2.0-0006990947 10.1016/S1007-5704(97)90008-3 ZBL0924.34063 Tatari M. Dehghan M. On the convergence of He's variational iteration method Journal of Computational and Applied Mathematics 2007 207 1 121 128 10.1016/j.cam.2006.07.017 2332953 ZBL1120.65112 Odibat Z. Momani S. Numerical methods for nonlinear partial differential equations of fractional order Applied Mathematical Modelling 2008 32 1 28 39 2-s2.0-34548384362 10.1016/j.apm.2006.10.025 ZBL1133.65116 Yulita Molliq R. Noorani M. S. M. Hashim I. Variational iteration method for fractional heat- and wave-like equations Nonlinear Analysis: Real World Applications 2009 10 3 1854 1869 10.1016/j.nonrwa.2008.02.026 2502991 ZBL1172.35302 Yulita Molliq R. Noorani M. S. M. Hashim I. Ahmad R. R. Approximate solutions of fractional Zakharov-Kuznetsov equations by VIM Journal of Computational and Applied Mathematics 2009 233 2 103 108 10.1016/j.cam.2009.03.010 2568507 ZBL1173.65066 Chun C. Variational iteration method for a reliable treatment of heat equations with ill-defined initial data International Journal of Nonlinear Sciences and Numerical Simulation 2008 9 4 435 440 2-s2.0-51749092612 Porshokouhi M. G. Ghanbari B. Application of He's variational iteration method for solution of the family of Kuramoto-Sivashinsky equations Journal of King Saud University - Science 2011 23 4 407 411 2-s2.0-77956109792 10.1016/j.jksus.2010.07.023 Mohyud-Din S. T. Yildirim A. An algorithm for solving the fractional vibration equation Computational Mathematics and Modeling 2012 23 2 228 237 10.1007/s10598-012-9133-2 El-Sayed A. M. A. Elsaid A. El-Kalla I. L. Hammad D. A homotopy perturbation technique for solving partial differential equations of fractional order in finite domains Applied Mathematics and Computation 2012 218 17 8329 8340 10.1016/j.amc.2012.01.057 2921326 ZBL1245.65141 Rosenau P. Hyman J. M. Compactons: solitons with finite wavelength Physical Review Letters 1993 70 5 564 567 2-s2.0-12044253199 10.1103/PhysRevLett.70.564 ZBL0952.35502 Gorenflo R. Mainardi F. Carpinteri A. Mainardi F. Fractional calculus: integral and differential equations of fractional order Fractals and Fractional Calculus 1997 New York, NY, USA Springer Caputo M. Linear models of dissipation whose Q is almost frequency independent II Geophysical Journal of the Royal Astronomical Society 1967 13 5 529 539 10.1111/j.1365-246X.1967.tb02303.x Inokuti M. Sekine H. Mura T. Nemat Nasser S. General use of the Lagrange multiplier in nonlinear mathematical physics Variational Method in the Mechanics of Solids 1978 New York, NY, USA Pergamon Press 156 162 Finlayson B. A. The Method of Weighted Residuals and Variational Principles 1972 New York, NY, USA Academic Press Clarkson P. A. Mansfield E. L. Priestley T. J. Symmetries of a class of nonlinear third-order partial differential equations Mathematical and Computer Modelling 1997 25 8-9 195 212 10.1016/S0895-7177(97)00069-1 1465776 ZBL0879.35005