# Feng’s First Integral Method Applied to the ZKBBM and the Generalized Fisher Space-Time Fractional Equations

Authors: Huitzilin Yépez-Martínez; Ivan  O. Sosa; Juan  M. Reyes
Journal: Journal of Applied Mathematics
 (2015)

---

## Abstract

The fractional derivatives in the sense of the modified Riemann-Liouville derivative and Feng’s first integral method are employed to obtain the exact solutions of the nonlinear space-time fractional ZKBBM equation and the nonlinear space-time fractional generalized Fisher equation. The power of this manageable method is presented by applying it to the above equations. Our approach provides first integrals in polynomial form with high accuracy. Exact analytical solutions are obtained through establishing first integrals. The present method is efficient and reliable, and it can be used as an alternative to establish new solutions of different types of fractional differential equations applied in mathematical physics.


---

## Body

## 1. Introduction

Fractional differential equations have been the focus of many studies due to their frequent appearance in various applications in physics, biology, engineering, signal processing, systems identification, control theory, finance, and fractional dynamics [1–3]. Recently, a large amount of literature has been provided to construct the solutions of fractional ordinary differential equations, fractional partial differential equations (PDEs), and integral equations of physical interest. Several powerful methods have been proposed to obtain approximate and exact solutions of fractional differential equations, such as the Adomian decomposition method [4, 5], the variational iteration method [6–8], the homotopy analysis method [9–12], the homotopy perturbation method [13–15], the Lagrange characteristic method [16], and the fractional subequation method [17].In [18], Jumarie proposed a modified Riemann-Liouville derivative. With this kind of fractional derivative and some useful formulas, we can convert fractional differential equations into integer-order differential equations by variable transformation.Feng [19] has introduced a reliable and effective method called Feng’s first integral method to look for travelling wave solutions of nonlinear PDEs. The basic idea of this method is to construct a first integral with polynomial coefficients of an explicit form to an equivalent autonomous planar system by using division theorem. This method in comparison with other methods has many advantages; it avoids a great deal of complicated and tedious calculation and provides exact and explicit travelling solutions with high accuracy.Feng’s first integral method [20–26] can be used to construct the exact solutions for some time fractional differential equations.Among the nonlinear PDEs there are two important examples of fundamental interest in mathematical-physical models: the nonlinear ZKBBM equation that describes the wave water phenomena in the large wavelength limit and one special case of nonlinear reaction-diffusion equation better known as the generalized Fisher equation. Let us review some general aspects about these two types of nonlinear PDEs.The well-known Zakharov-Kuznetsov-Benjamin-Bona-Mahony equation (ZKBBM equation) also referred to as the regularized long wave equation (see [27, 28])
(1)Utx,t+Uxx,t-2aUx,tUxx,t-bUtxxx,t=0,t>0,
was originally derived for water waves under the assumption of small wave amplitude and large wavelength. The ZKBBM equation describes approximately the unidirectional propagation of long waves in certain nonlinear dispersive systems and has been introduced in order to give a satisfactory description for the small wave amplitude Ux,t in the large wavelength regime in many other physical systems. It has been used to account adequately for observable phenomena such as the interaction of solitary waves and dissipationless undular shocks; see [27].The details of the derivation of the ZKBBM equation differ, in different examples, but in general the essentials are as follows. The first property of the systems in question is that dispersive effects on infinitesimal waves vanish in the limit as wavelength becomes infinite, and the limiting phase speed is a constant (c0>0). Thus waves of extreme length, along the +x propagation direction, are described by
(2)Utx,t+c0Uxx,t=0,t>0,
whose general solution is an arbitrary differentiable function of x-c0t, where Ut=∂U/∂t and Ux=∂U/∂x. In certain applications (2) already has some validity as an approximation for real waves of sufficiently small amplitude and large wavelength but is not a valid approximation for very large times, during which nonlinear and frequency dispersive effects are accumulated to a significant level. These effects require the introduction of at least a first order approximation model to extend the range of validity of the model (1).For the linearized equation (1) all values of Ux,t are propagated with the same velocity (dx/dt)=c0. When the nonlinear effects are taken into account, the characteristic velocity becomes dependent linearly on Ux,t for the waves of finite but small amplitude; therefore in the first order approximation we obtain
(3)1c0dxdtu=const.=1+βUx,t,
where β is a constant. After the substitution of (3) into (2) the following equation is obtained:
(4)Utx,t+c0Uxx,t+c0βUx,tUxx,t=0,t>0.
Additionally the dispersive effects as suffered by waves of finite but large wavelength in the model (2) are ignored. In the first order approximation the dispersive effects transform the model (2) into the following equation:
(5)Utx,t+c0Ux,t+αUtxx,tx=0,t>0.
Although for a self-consistent theory the two kinds of small effects have to be considered simultaneously, the outcome of such theories can generally be anticipated by considering the two effects separately. That is, in a first approximation, the above corrections are added separately into (2), and higher order terms appear negligible in this approximation.When the physical parameters and the scaling factors presented in a particular example are appropriately absorbed into the definition of the dependent variableUx,t and the independent variables x and t, which are, respectively, proportional to distance in the physical system and to time, the ZKBBM equation is given by: (6)Utx,t+Uxx,t-2aUx,tUxx,t-bUtxxx,t=0,t>0.
Equation (6) can be extended in order to consider fractional order derivatives obtaining the ZKBBM space-time fractional equation [29]:
(7)DtαUx,t=-DxαUx,t+2aUx,tDxαUx,t+bDtαDx2αUx,tt>0,0<α≤1,
where DxαUx,t and DtαUx,t denotes the fractional derivative of order α with respect to the independent variables x and t.The second nonlinear PDE mentioned before, the generalized Fisher equation [30–32], can be obtained as a special case of the reaction-diffusion equation. The reaction-diffusion equation plays a fundamental role in a great number of various models of heat and reaction-diffusion processes, mathematical biology, chemistry, genetics, and many other fields.We start with the reaction-diffusion equation with power nonlinearity (see, for example, [33]):
(8)Utx,t-Uxxx,t=fUx,t,t>0;
here fUx,t is a sufficiently smooth polynomial function satisfying the relations f0=f1=0.One particular expression for the functionfUx,t in (8) which satisfies these conditions is [33]
(9)Utx,t-Uxxx,t=k-k+1Unx,t+λ1Ux,t+λ2Un+1/2x,t+λ3U3-n/2x,t+λ4U2-nx,t,t>0,
where λ1,…,λ4 are constants that satisfy the condition imposed over fUx,t and n is an arbitrary constant. Equation (9) defines a special class of the nonlinear heat equation, which includes many important applications in some of the areas mentioned before. If we consider the special case in (9) where λ2=λ3=λ4=0, n=6, and λ1=k(k+1)=1, we obtain the generalized Fisher equation:
(10)Utx,t-Uxxx,t=Ux,t1-U6x,t,t>0,
where in the original mathematical biology model U(x,t) denotes the population density and Ux,t1-U6x,t represents the population supply due to births and deaths.The above equation has been introduced as a model for the propagation of a mutant gene in mathematical biology [30–32] and can be extended to fractional order derivatives obtaining the generalized Fisher space-time fractional equation:
(11)DtαUx,t=Dx2αUx,t+Ux,t1-U6x,t,t>0,0<α≤1.
The present work investigates the applicability and effectiveness of Feng’s first integral method to obtain new exact analytical solutions for the nonlinear space-time fractional Zakharov-Kuznetsov-Benjamin-Bona-Mahony (ZKBBM) partial differential equation and the nonlinear space-time fractional generalized Fisher equation.We will show that, for the analytical solutions of the ZKBBM fractional partial differential equation, one of the advantages of Feng’s first integral method is the evaluation of the constants involved in the analytical solutions that have not been evaluated within the subequation method [29]. Furthermore for the analytical solutions of the space-time fractional generalized Fisher equation, by applying Feng’s first integral method, we will obtain three new analytical solutions that have not been obtained in previous works [30, 31, 33].


## 2. The Modified Riemann-Liouville Derivative and Feng’s First Integral Method

In this section we present the main ideas of Feng’s first integral method. This method considers the Jumarie modified Riemann-Liouville fractional derivative of orderα; we first give some definitions and properties of the modified Riemann-Liouville derivative which are used further in this paper.Assume thatf:R→R, x→f(x) denotes a continuous (but not necessarily differentiable) function. The Jumarie modified Riemann-Liouville derivative of order α is defined by the expression [18]
(12)Dxαfx=1Γ1-α∫0xx-ξ-α-1·fξ-f0dξ,α<0,1Γ1-αddx∫0xx-ξ-α·fξ-f0dξ0<α<1,fα-nxnn≤α≤n+1,n≥1.Some properties of the fractional modified Riemann-Liouville derivative are(13)Dxαxγ=Γγ+1Γγ+1-αxγ-α,Dxαfxgx=gxDxαfx+fxDxαgx,Dxαfgx=fg′gxDxαgx=Dgαfgxg′xα.Now in order to introduce Feng’s first integral method [19], let us consider the space-time fractional differential equation with independent variables x1,x2,…,xm,t and some function u:
(14)Fu,Dtαu,Dx1αu,Dx2αu,Dx3αu,…,Dt2αu,FDx12αu,Dx22αu,Dx32αu,…=0.
Using the variable transformation
(15)ux1,x2,…,xm,t=U(ξ),ξ=l1x1α+l2x2α+⋯+lmxmα+λtαΓ(1+α),
where li and λ are constants to be determined later, the fractional differential equation (14) is reduced to a nonlinear ordinary differential equation
(16)H=HUξ,U′ξ,U′′ξ,…,
where U′ξ=dUξ/dξ.We assume that (16) has a solution in the form
(17)Uξ=Xξ,
and we introduce a new independent variable Y(ξ)=X′ξ, which leads to the following new system of equations:
(18)X′(ξ)=Y(ξ),Y′(ξ)=GX(ξ),Y(ξ).
Now, let us introduce the central idea of Feng’s first integral method. By using the division theorem for two variables in the complex domain C which is based on the Hilbert-Nullstellensatz theorem [34], we can obtain one first integral to (18) which can reduce (16) to a first order integrable ordinary differential equation. An exact solution to (14) is then obtained by solving this equation directly.Division Theorem. Suppose that P(x,y) and Q(x,y) are polynomials in C[x,y], and P(x,y) is irreducible in C[x,y]. If Q(x,y) vanishes at all zero points of P(x,y), then there exists a polynomial H(x,y) in C[x,y] such that
(19)Q(x,y)=P(x,y)H(x,y).


## 3. Applications

In this section, we illustrate the applicability of Feng’s first integral method to solve the nonlinear fractional partial differential equations (7) and (11).Example 1.
As a first example we consider the space-time fractional ZKBBM equation (7):
(20)DtαUx,t=-DxαUx,t+2aUx,tDxαUx,t+bDtαDx2αUx,t,t>0,0<α≤1,
where a and b are arbitrary constants. By considering the travelling wave transformation
(21)Ux,t=U(ξ),with:ξ=lxα+λtαΓ(1+α),
where l and λ are constants, substituting (21) into (20), we can reduce (20) into an ordinary differential equation (ODE):
(22)λdUξdξ=-ldUξdξ+2alUξdUξdξ+bλl2ddξd2Uξdξ2,
or
(23)d3Uξdξ3=λ+lbλl2dUξdξ-2abλlUξdUξdξ.
Integrating (23) and γ being the integration constant, we have
(24)d2Uξdξ2=λ+lbλl2Uξ-abλlU2ξ+γ,
and using (17) and (18), (24) is equivalent to the two-dimensional autonomous system
(25)dXξdξ=Yξ,dYξdξ=λ+lbλl2Xξ-abλlX2ξ+γ.
According to Feng’s first integral method, we suppose that X(ξ) and Y(ξ) are nontrivial solutions of (25), and Q(X,Y) is an irreducible polynomial in the complex domain C such that
(26)QXξ,Yξ=∑i=0maiXξYiξ=0,
where ai(X)(i=0,1,…,m) are polynomials of X and am(X)≠0. Due to the division theorem, there exists a polynomial g(X)+h(X)Y in the complex domain C[x,y] such that
(27)dQdξ=∂Q∂XdXdξ+∂Q∂YdYdξ=gX+hXYX∑i=0maiXYi=0.
Suppose that m=1; by equating the coefficients of Yi(i=2,1,0) on both sides of (27), we have(28a)a˙1X=hXa1X,(28b)a˙0X=gXa1X+hXa0X,(28c)a0XgX=a1Xλ+lbλl2Xξ-abλlX2ξ+γ=a1XεXξ-βX2ξ+γ,where
(29)ε=λ+lbλl2,β=abλl;
since ai(X)(i=0,1) are polynomials, then from (28a) we deduce that a1(X) is constant and h(X)=0. For simplicity, we take a1(X)=1. Balancing the degrees of g(X) and a0(X), we conclude that (30a)gX=gm0Xm0+gm1Xm1,(30b)a0X=gm0Xm0+1m0+1+gm1Xm1+1m1+1+B0,where B0 is an arbitrary constant. Substituting a0(X), a1(X), g(X), and h(X) into (28c) and setting all the coefficients of powers of X to be zero, we obtain a system of nonlinear algebraic equations and by solving it, we obtain
(31)m0=-12,m1=12,B0=0,γ=gm02m0+1=2gm02,
where
(32)gm1=-32β1/2,gm0=38ε-23β1/2orgm1=--32β1/2,gm0=-38ε-23β1/2.
From the above conditions, four solutions can be obtained for the space-time fractional ZKBBM partial differential equation.
Case  1. When gm0=(3/8)ε-2/3β1/2, gm1=-(3/2)β1/2, ε<0, and β<0,
(33)a0X=2gm0X1/2+23gm1X3/2.
Now from the condition QX,Y=0, we obtain
(34)0=a0X+a1XYξ
and therefore taking into account that a1(X)=1 and from (33) it follows that
(35)Yξ=-2gm0X1/2+23gm1X3/2.
Combining this first integral Y(ξ) with the two-dimensional autonomous system of (18), the second order differential equation (24) can be reduced to the following first order differential equation:
(36)dXξdξ=-2gm0X1/2+23gm1X3/2.
We consider the following transformation:
(37)Xξ=Z2ξ,
introduced by Wang [35]. With this transformation (36) can be written as
(38)dZξdξ=-gm0+13gm1Z2
or
(39)dZξdξ=ε332β1/2-β61/2Z2,
where ε=-ε and β=-β.
By solving (39), the following solutions are obtained:
(40)Zξ=3εβ2tanhε/8ξβ,
or
(41)Zξ=3εβ2coth⁡ε/8ξβ,
if we consider the case l>0 and λ>0, and taking into account (29) and (37), the solutions (40) and (41) reduce to
(42)Xξ=6σbλlatanh2-σξ,
or
(43)Xξ=6σbλlacoth⁡2-σξ,
with σ=-ε/8=(l+λ)/8bl2λ, where b<0 and a>0.
Case  2. When gm0=(3/8)ε-2/3β1/2, gm1=-(3/2)β1/2, ε>0, and β<0,
(44)a0X=2gm0X1/2+23gm1X3/2.
Taking into account that a1(X)=1 and from (44) it follows that
(45)Yξ=-2gm0X1/2+23gm1X3/2.
Combining this first integral Y(ξ) with the two-dimensional autonomous system given by (18), the second order differential equation (24) can be reduced to the following first order differential equation:
(46)dXξdξ=-2gm0X1/2+23gm1X3/2.
By using the transformation defined in (37), (46) can be written as
(47)dZξdξ=-ε332β1/2-β61/2Z2,
where β=-β.
By solving (47), the following solutions are obtained:
(48)Zξ=-3εβ2tanε/8ξβ,
or
(49)Zξ=-3εβ2cotε/8ξβ,
if we consider the case l>0 and λ>0, and taking into account (29) and (37), the solutions (48) and (49) reduce to
(50)Xξ=-6σbλlatan2σξ,
or
(51)Xξ=-6σbλlacot2σξ,
with σ=ε/8=(l+λ)/8bl2λ, where b>0 and a<0.
Case  3. When gm0=-(3/8)ε-2/3β1/2, gm1=--(3/2)β1/2, ε<0, and β<0,
(52)a0X=2gm0X1/2+23gm1X3/2.
Taking into account that a1(X)=1 and from (52) it follows that
(53)Yξ=-2gm0X1/2+23gm1X3/2.
Combining this first integral Y(ξ) with the two-dimensional autonomous system given by (18), the second order differential equation (24) can be reduced to the following first order differential equation:
(54)dXξdξ=-2gm0X1/2+23gm1X3/2.
By using the transformation defined in (37), (54) can be written as
(55)dZξdξ=-ε332β1/2+β61/2Z2,
where ε=-ε and β=-β.
By solving (55), the following solutions are obtained:
(56)Zξ=-3εβ2tanhε/8ξβ,
or
(57)Zξ=-3εβ2coth⁡ε/8ξβ,
if we consider the case l>0 and λ>0, and taking into account (29) and (37), the solutions (56) and (57) reduce to
(58)Xξ=6σbλlatanh2-σξ,
or
(59)Xξ=6σbλlacoth⁡2-σξ,
with σ=-ε/8=(l+λ)/8bl2λ, where b<0 and a>0.
Case  4. When gm0=-(3/8)ε-2/3β1/2, gm1=--(3/2)β1/2, ε>0, and β<0,
(60)a0X=2gm0X1/2+23gm1X3/2.
Taking into account that a1(X)=1 and from (60) it follows that
(61)Yξ=-2gm0X1/2+23gm1X3/2.
Combining this first integral Y(ξ) with the two-dimensional autonomous system given by (18), the second order differential equation (24) can be reduced to the following first order differential equation:
(62)dXξdξ=-2gm0X1/2+23gm1X3/2.
By using the transformation defined in (37), (62) can be written as
(63)dZξdξ=+ε332β1/2+β61/2Z2,
where β=-β.
By solving (63), the following solutions are obtained:
(64)Zξ=3εβ2tanε/8ξβ,
or
(65)Zξ=3εβ2cotε/8ξβ,
if we consider the case l>0 and λ>0, and taking into account (29) and (37), the solutions (64) and (65) reduce to
(66)Xξ=-6σbλlatan2σξ,
or
(67)Xξ=-6σbλlacot2σξ,
with σ=ε/8=(l+λ)/8bl2λ, where b>0 and a<0.
For the nonlinear space-time fractional ZKBBM equation, we have successfully recovered the previously known solutions, (42), (43), (48), (49), (58), (59), (66), and (67), that have been found in [29]; furthermore, Feng’s first integral method has evaluated the constant σ which has not been evaluated in this previous work.Example 2.
Next we consider the nonlinear fractional generalized Fisher equation (11):
(68)DtαUx,t=Dx2αUx,t+Ux,t1-U6x,t,t>0,0<α≤1.
By considering the following transformation,
(69)U3x,t=yx,t,
introduced by Wang [35], (68) can be written as
(70)yx,tDtαyx,t=yx,tDx2αyx,t-23Dxαyx,t2+3y2x,t1-y2x,t,t>0,0<α≤1.
By considering the travelling wave transformation,
(71)yx,t=y(ξ),with:ξ=lxα+λtαΓ(1+α),
where l and λ are constants, (70) can be reduced to the following nonlinear ordinary differential equation:
(72)λyξdyξdξ=l2yξd2yξdξ2-23l2dyξdξ2+3y2ξ1-y2ξ,0<α≤1.
We assume that (72) has a solution in the form yξ=Xξ, and using (18), (72) is equivalent to the two-dimensional autonomous system:
(73)dXξdξ=Yξ,l2XξdYξdξ=λXξYξ+23l2Yξ2-3X2ξ1-X2ξ,0<α≤1.
According to Feng’s first integral method, we suppose that Xξ and Yξ are nontrivial solutions of (73), and QX,Y is an irreducible polynomial in the complex domain C such that
(74)QXξ,Yξ=∑i=0maiXξYiξ=0,
where ai(X)(i=0,1,…,m) are polynomials of X and am(X)≠0. Due to the division theorem, there exists a polynomial g(X)+h(X)Y in the complex domain CX,Y such that
(75)dQdξ=∂Q∂XdXdξ+∂Q∂YdYdξ=gX+hXYX∑i=0maiXYi=0.
Suppose that m=1; by equating the coefficients of Yi(i=2,1,0) on both sides of (75), we have(76a)a˙1XXξ=a1XhXXξ-23,(76b)a˙0X+a1Xλl2=a1XgX+a0XhX,(76c)a0XgXXξl2=-3a1X2ξ1-X2ξ.Since ai(X)(i=0,1) are polynomials, then from (76a) we deduce that a1(X) is constant and h(X)=2/3X. For simplicity, we take a1(X)=1. Balancing the degrees of g(X) and a0(X), we conclude that
(77)gX=A0+B0X,a0X=A0+B0X+C0X22.
Substituting a0(X), a1(X), g(X), and h(X) into (76c) and setting all the coefficients of powers of X to be zero, we obtain a system of nonlinear algebraic equations and by solving it, we obtain
(78)A0=0,B0=3A0-λl2,C0=32B0,λ=54A0l2,A02=4l2,B02=4l2;
therefore we obtain the following set of solutions:
(79)A0=±2l,B0=±2l.
From the above conditions, four solutions can be obtained.
Case  1. When A0=2/l and B0=2/l,
(80)a0X=-32lX+32lX2.
Now from the condition QX,Y=0 we obtain
(81)0=a0X+a1XYξ,
and therefore taking into account that a1X=1 and from (81) it follows that
(82)Yξ=32lX-32lX2.
If we consider the definition given in (18) one finally obtains
(83)dXξdξ=32lX-32lX2.
By solving the above equation the solution to (72) is obtained:
(84)Xξ=11+c1e-3ξ/2l;
if we consider the case c1=1, the solution (84) reduces to
(85)Xξ=121+tanh3ξ4l,
with λ=5l/2. Taking into account (69), the solution to the Fisher equation (68) is given by
(86)Ux,t=121+tanh3ξ4l1/3,
with
(87)ξ=lxα+λtαΓ1+α.
Case  2. When A0=-2/l and B0=2/l,
(88)a0X=32lX+32lX2.
Now from the condition QX,Y=0 we obtain
(89)0=a0X+a1XYξ,
and therefore taking into account that a1X=1 and from (88), it follows that
(90)Yξ=-32lX+32lX2.
If we consider the definition given in (18), one finally obtains
(91)dXξdξ=-32lX+32lX2.
By solving the above equation, the solution to (72) is obtained:
(92)Xξ=1c1e3ξ/2l-1;
if we consider the case c1=1, the solution (92) reduces to
(93)Xξ=-121-coth⁡3ξ4l,
with λ=-5l/2. Taking into account (69), the solution to the Fisher equation (68) is given by
(94)Ux,t=-121-coth⁡3ξ4l1/3,
with
(95)ξ=lxα+λtαΓ1+α.
Case  3. When A0=2/l and B0=-2/l,
(96)a0X=-32lX+32lX2.
Now from the condition QX,Y=0, we obtain
(97)0=a0X+a1XYξ,
and therefore taking into account that a1X=1 and from (96), it follows that
(98)Yξ=32lX+32lX2.
If we consider the definition given in (18), one finally obtains
(99)dXξdξ=32lX+32lX2.
By solving the above equation, the solution to (72) is obtained:
(100)Xξ=1c1e-3ξ/2l-1;
if we consider the case c1=1, the solution (100) reduces to
(101)Xξ=-121+coth⁡3ξ4l,
with λ=5l/2. Taking into account (69), the solution to the Fisher equation (68) is given by
(102)Ux,t=-121+coth⁡3ξ4l1/3,
with
(103)ξ=lxα+λtαΓ(1+α).
Case  4. When A0=-2/l and B0=-2/l,
(104)a0X=32lX-32lX2.
Now from the condition QX,Y=0, we obtain
(105)0=a0X+a1XYξ,
and therefore taking into account that a1X=1 and from (104), it follows that
(106)Yξ=-32lX-32lX2.
If we consider the definition given in (18), one finally obtains
(107)dXξdξ=-32lX-32lX2.
By solving the above equation, the solution to (72) is obtained:
(108)Xξ=1c1e3ξ/2l+1;
if we consider the case c1=1, the solution (108) reduces to
(109)Xξ=121-tanh3ξ4l,
with λ=-5l/2. Taking into account (69), the solution to the Fisher equation (68) is given by
(110)Ux,t=121-tanh3ξ4l1/3,
with
(111)ξ=lxα+λtαΓ(1+α).
For the nonlinear space-time generalized fractional Fisher equation, we have successfully recovered the previously known solution (110) that has been found in [31], but to the best of our knowledge the solutions (86), (94), and (102) have not been obtained previously in the literature.


## 4. Conclusions

Feng’s first integral method was applied successfully to obtain new exact analytical solutions of the nonlinear space-time fractional ZKBBM equation and the nonlinear space-time fractional generalized Fisher equation. These solutions can be very useful as a starting point of comparison when some approximate methods are applied to these nonlinear space-time fractional equations. The performance of Feng’s first integral method is reliable and effective to obtain new solutions. This method has more advantages: it is direct and concise. Thus, the proposed method can be extended to solve many systems of nonlinear fractional partial differential equations in mathematical and physical sciences.



---

*Source: 191545-2015-03-02.xml*