# New Solutions for (1+1)-Dimensional and (2+1)-Dimensional Ito Equations

Authors: A. H. Bhrawy; M. Sh. Alhuthali; M. A. Abdelkawy
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

Using the extended F-expansion method based on computerized symbolic computation technique, we find several new solutions of (1+1)-dimensional and (2+1)-dimensional Ito equations. These solutions contain hyperbolic and triangular solutions. It is shown that the power of the extended F-expansion method is its ease of use to determine shock or solitary type of solutions. In addition, as an illustrative sample, the properties for the extended F-expansion solutions of the Ito equations are shown with some figures.


---

## Body

## 1. Introduction

The nonlinear wave phenomena can be observed in various scientific fields, such as plasma physics, optical fibers, fluid dynamics, and chemical physics. The nonlinear wave phenomena can be obtained in solutions of nonlinear evolution equations (NEEs). The study of NLEEs appear everywhere in applied mathematics and theoretical physics including engineering sciences and biological sciences. These NLEEEs play a key role in describing key scientific phenomena. For example, the nonlinear Schrödinger’s equation describes the dynamics of propagation of solitons through optical fibers. The Korteweg-de Vries equation models the shallow water wave dynamics near ocean shore and beaches. Additionally, the Schrödinger-Hirota equation describes the dispersive soliton propagation through optical fibers. These are just a few examples in the whole wide world of NLEEs and their applications, (see, for instance, [1–4]). While the above mentioned NLEEs are scalar NLEEs, there is a large number of NLEEs that are coupled. Some of them are two-coupled NLEEs such as the Gear-Grimshaw equation [2], while there are several others that are three-coupled NLEEs. An example of a three-coupled NLEE is the Wu-Zhang equation [4]. These coupled NLEEs are also studied in various areas of theoretical physics as well.The exact solutions of these NEEs play an important role in the understanding of nonlinear phenomena. In the past decades, many methods were developed for finding exact solutions of NEEs such as the inverse scattering method [5, 6], improved projective Riccati equations method [7, 8], Cole-Hopf transformation method [9], exp-function method [10–16], bifurcation theory method [17], (G′/G)-expansion method [18, 19], homotopy perturbation method [20], tanh function method [20–24], and Jacobi and Weierstrass elliptic function method [25, 26]. Although Porubov et al. [27–29] have obtained some exact periodic solutions to some nonlinear wave equations, they use the Weierstrass elliptic function and involve complicated deducing. A Jacobi elliptic function (JEF) expansion method, which is straightforward and effective, was proposed for constructing periodic wave solutions for some nonlinear evolution equations. The essential idea of this method is similar to the tanh method by replacing the tanh function with some JEFs such as sn, cn, and dn. For example, the Jacobi periodic solution in terms of sn may be obtained by applying the sn-function expansion. Many similarl repetitious calculations have to be done to search for the Jacobi doubly periodic wave solutions in terms of cn and dn [30].Recently, F-expansion method [31–34] was proposed to obtain periodic wave solutions of NLEEs, which can be thought of as a concentration of JEF expansion since F here stands for every function of JEFs. The objectives of this work are twofold. First, we seek to extend others works to establish new exact solutions of distinct physical structures for the nonlinear equations (1.1) and (1.2). The extended F-expansion (EFE) method will be used to achieve the first goal. The second goal is to show that the power of the EFE method is its ease of use to determine shock or solitary type of solutions. In this paper, we study two well-known PDEs, namely, generalized (1+1)-dimensional and generalized (2+1)-dimensional Ito equations. Many studies are concerning the (1+1)-dimensional Ito equation and the (2+1)-dimensional Ito equation [35–42].The history of the KdV equation started with experiments by John Scott Russell in 1834, followed by theoretical investigations by Lord Rayleigh and Joseph Boussinesq around 1870, and, finally, Korteweg and de Vries in 1895 [39]. The KdV equation was not studied much after this until Zabusky and Kruskal (1965) [40] discovered numerically that its solutions seemed to decompose at large times into a collection of “solitons”: well-separated solitary waves. Ito [41, 42] obtained the well-known generalized (1+1)-dimensional and generalized (2+1)-dimensional Ito equations by generalization of the bilinear KdV equation as
(1.1)utt+uxxxt+3(uxut+uuxt)+3uxx∫∞x‍utdx′=0,(1.2)utt+uxxxt+3(2uxut+uuxt)+3uxx∫∞x‍utdx′+αuyt+βuxt=0.
Also Sawada-Kotera-Ito (SK-Ito) seventh-order equation is the special case of the generalized seventh-order KdV equation as
(1.3)ut+252u2ux+63u3x+378uuxu2x+126u2u3x+63u2xu3x+42uxu4x+21uu5x+u7x=0.
SK-Ito equation is characterized by the presence of three dispersive terms ux,  u3x, and u7x, respectively. SK-Ito seventh-order equation is completely integrable and admits of conservation laws [43]. Moreover, the Ito-type coupled KdV (ItcKdV) equation [44], written in the following form:
(1.4)ut+αuux+βvvx+γuxxx=0;vt+β(uv)x=0,
if we take the special values α=-6,  β=-2, and γ=-1. Equation (1.4) describes the interaction process of two internal long waves which has infinitely many conserved quantities [45, 46].In this paper, we extend the EFE method with symbolic computation to (1.1) and (1.2) for constructing their interesting Jacobi doubly periodic wave solutions. It is shown that soliton solutions and triangular periodic solutions can be established as the limits of Jacobi doubly periodic wave solutions. In addition, the algorithm that we use here is also a computerized method, in which we are generating an algebraic system.


## 2. Extended F-Expansion Method

In this section, we introduce a simple description of the EFE method, for a given partial differential equation as(2.1)G(u,ux,uy,uz,uxy,…)=0.
We like to know whether travelling waves (or stationary waves) are solutions of (2.1). The first step is to unite the independent variables x,  y, and t into one particular variable through the new variable as
(2.2)ζ=x+y-νt,u(x,y,t)=U(ζ),
where ν is wave speed, and reduce (2.1) to an ordinary differential equation (ODE) as
(2.3)G(U,U′,U′′,U′′′,…)=0.
Our main goal is to derive exact or at least approximate solutions, if possible, for this ODE. For this purpose, let us simply use U as the expansion in the form
(2.4)u(x,y,t)=U(ζ)=∑i=0N‍aiFi+∑i=1N‍a-iF-i,
where
(2.5)F′=A+BF2+CF4,
the highest degree of (dpU/dζp), is taken as
(2.6)O(dpUdζp)=N+p,p=1,2,3,…,(2.7)O(UqdpUdζp)=(q+1)N+p,q=0,1,2,…,p=1,2,3,…,
where A,  B, and C are constants, and N in (2.3) is a positive integer that can be determined by balancing the nonlinear term(s) and the highest order derivatives. Normally N is a positive integer, so that an analytic solution in closed form may be obtained. Substituting (2.1)–(2.5) into (2.3) and comparing the coefficients of each power of F(ζ) in both sides, we will get an overdetermined system of nonlinear algebraic equations with respect to ν,a0,a1,…. We will solve the over-determined system of nonlinear algebraic equations by use of Mathematica. The relations between values of A,  B,  C, and corresponding JEF solution F(ζ) of (2.4) are given in Table 1. Substituting the values of A,  B,  C, and the corresponding JEF solution F(ζ) chosen from Table 1 into the general form of solution, then an ideal periodic wave solution expressed by JEF can be obtained.Table 1
Relation between values of (A,B,C) and corresponding F.
A
B
C
F
(
ζ
)
1
-
1
-
m
2
m
2
sn
(
ζ
) or cd(ζ)=cn(ζ)/dn(ζ)
1
-
m
2
2
m
2
-
1
-
m
2
cn
(
ζ
)
m
2
-
1
2
-
m
2
−1
dn
(
ζ
)
m
2
-
1
-
m
2
1
ns
(
ζ
)
(
1
/
sn
(
ζ
)
) or dc(ζ) = dn(ζ)/cn(ζ)
-
m
2
2
m
2
-
1
1
-
m
2
nc
(
ζ
)
=
1
/
cn
(
ζ
)
−1
2
-
m
2
m
2
-
1
nd
(
ζ
)
=
1
/
dn
(
ζ
)
1
2
-
m
2
1
-
m
2
sc
(
ζ
)
=
sn
(
ζ
)
/
cn
(
ζ
)
1
2
m
2
-
1
-
m
2
(
-
1
-
m
2
)
sd
(
ζ
)
=
sn
(
ζ
)
/
dn
(
ζ
)
1
-
m
2
2
-
m
2
1
cs
(
ζ
)
=
cn
(
ζ
)
/
sn
(
ζ
)
-
m
2
(
1
-
m
2
)
2
m
2
-
1
1
ds
(
ζ
)
=
dn
(
ζ
)
/
sn
(
ζ
)
1
/
4
(
1
-
2
m
2
)
/
2
1
/
4
ns
(
ζ
)
+
cs
(
ζ
)
(
1
-
m
2
)
/
4
(
1
+
m
2
)
/
2
(
1
-
m
2
)
/
2
nc
(
ζ
)
+
sc
(
ζ
)
1
/
4
(
m
2
-
2
)
/
2
m
2
/
4
ns
(
ζ
)
+
ds
(
ζ
)
m
2
/
4
(
m
2
-
2
)
/
2
m
2
/
4
sn
(
ζ
)
+
i
cs
(
ζ
)sn
(
ζ
), cn(ζ), and dn(ζ) are the JE sine function, JE cosine function, and the JEF of the third kind, respectively. And
(2.8)cn2(ζ)=1-sn2(ζ),dn2(ζ)=1-m2sn2(ζ),
with the modulus m(0<m<1).Whenm→1. the Jacobi functions degenerate to the hyperbolic functions, that is,
(2.9)snζ→tanhζ,cnζ→sechζ,dnζ→sechζ,
when m→0, the Jacobi functions degenerate to the triangular functions, that is,
(2.10)snζ→sinζ,cnζ→cosζ,dn→1.


## 3. Generalized (1+1)-Dimensional Ito Equation

We first consider the generalized (1+1)-dimensional Ito equation (1.1) as follows:
(3.1)utt+uxxxt+3(uxut+uuxt)+3uxx∫∞x‍utdx′=0,
if we use the transformation u=vx, it carries (3.1) into
(3.2)vxtt+vxxxxt+3(vxxvxt+vxvxxt)+3vxxxvt=0,
if we use ζ=x-νt transforms (3.2) into the ODE, we have
(3.3)-νV′′′+V(v)-3ν(V′′V′′+V′V′′′)-3νV′′′V′=0,
where by integrating once we obtain, upon setting the constant of integration to zero,
(3.4)V′′′+3(V′)2-νV′=0,
if we use the transformation W=V′, then (3.4) can be written as follows:
(3.5)W′′+3W2-νW=0.
Balancing the term W′′ with the term W2 we obtain N=2 then
(3.6)W(ζ)=a0+a1ψ+a-1ψ-1+a2ψ2+a-2ψ-2,ψ′=A+Bψ2+Cψ4.Substituting (3.6) into (3.5) and comparing the coefficients of each power of ψ in both sides, we will get an over-determined system of nonlinear algebraic equations with respect to ν,ai, i=0,1,-1,-2,2. Solving the over-determined system of nonlinear algebraic equations by use of Mathematica, we obtain three groups of constants(1)(3.7)a-1=a-1=0,a0=-2B3,a2=-2C,a-2=-2A,ν=-2(B2+12AC)B,(2)(3.8)a-1=a2=a-1=0,a0=-2B3,a-2=-2A,ν=-2(B2-3AC)B,(3)(3.9)a-1=a-2=a-1=0,a0=-2B3,a2=-2C,ν=-2(B2-3AC)B.
The solutions of (3.1) are
(3.10)u1=-2(1+m2)3-2m2sn2(x-2((1+m2)2+12m2)1+m2t)-2ns2(x-2((1+m2)2+12m2)1+m2t),(3.11)u2=-2(1+m2)3-2m2cd2(x-2((1+m2)2+12m2)1+m2t)-2dc2(x-2((1+m2)2+12m2)1+m2t),(3.12)u3=-2(2m2-1)3+2m2cn2(x+2(12m2(m2-1)+(2m2-1)2)2m2-1t)-2(1-m2)nc2(x+2(12m2(m2-1)+(2m2-1)2)2m2-1t),(3.13)u4=-2(2-m2)3+2dn2(x-2((2-m2)2-12(-1+m2))2-m2t)-2(m2-1)nd2(x-2((2-m2)2-12(-1+m2))2-m2t),(3.14)u5=-2(1+m2)3-2ns2(x-2((1+m2)2-3m2)1+m2t),(3.15)u6=-2(2-m2)3-2(1-m2)sc2(x+2(12(1-m2)+(2-m2)2)2-m2t)-2cs2(x+2(12(1-m2)+(2-m2)2)2-m2t),(3.16)u7=2(2m2-1)3+2m2(1+m2)sd2(x+2(-12m2(1-m2)+(-1+2m2)2)2m2-1t)-2ds2(x+2(-12m2(1-m2)+(-1+2m2)2)2m2-1t),(3.17)u8=2(1-2m2)6-12(ns(x+2(0.75+(0.5-m2)2)0.5-m2t)+cs(x+2(0.75+(0.5-m2)2)0.5-m2t))2-12(ns(x+2(0.75+(0.5-m2)2)0.5-m2t)+cs(x+2(0.75+(0.5-m2)2)0.5-m2t))-2,(3.18)u9=-1+m23+1-m22×(nc(x+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5t)+sc(x+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5t))2+1-m22(nc(x+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5t)+sc(x+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5t))-2,(3.19)u10=-m2-23-m22(ns(x+2(0.75m2+(-1+0.5m2)2)0.5m2-1t)+ds(x+2(0.75m2+(-1+0.5m2)2)0.5m2-1t))2-12(ns(x+2(0.75m2+(-1+0.5m2)2)0.5m2-1t)+ds(x+2(0.75m2+(-1+0.5m2)2)0.5m2-1t))-2,(3.20)u11=m2-23-m22(sn(x+2(0.75m4+(-1+0.5m2)2)0.5m2-1t)+ics(x+2(0.75m4+(-1+0.5m2)2)0.5m2-1t))2+m22(sn(x+2(0.75m4+(-1+0.5m2)2)0.5m2-1t)+ics(x+2(0.75m4+(-1+0.5m2)2)0.5m2-1t))-2,(3.21)u12=-2(1+m2)3-2dc2(x-2((1+m2)2-3m2)1+m2t),(3.22)u13=-2(2m2-1)3-2(1-m2)nc2(x+2(-3m2(m2-1)+(2m2-1)2)2m2-1t),(3.23)u14=-2(2-m2)3-2(m2-1)nd2(x-2((2-m2)2+3(-1+m2))2-m2t),(3.24)u15=-2(2-m2)3-2cs2(x+2((2-m2)2-3(1-m2))2-m2t),(3.25)u16=2(2m2-1)3-2ds2(x+2(3m2(1-m2)+(-1+2m2)2)2m2-1t),(3.26)u17=2(1-2m2)6-12(ns(x+2((0.5-m2)2-3/16)0.5-m2t)+cs(x+2((0.5-m2)2-(3/16))0.5-m2t))-2,(3.27)u18=-1+m23+1-m22×(nc(x+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5t)+sc(x+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5t))-2,(3.28)u19=-m2-23-12(ns(x+2((0.5m2-1)2-3/16)0.5m2-1t)+ds(x+2((0.5m2-1)2-3/16)0.5m2-1t))-2,(3.29)u20=m2-23+m22(sn(x+2((0.5m2-1)2-m4/16)0.5m2-1t)+ics(x+2((0.5m2-1)2-m4/16)0.5m2-1t))-2,(3.30)u21=-2(1+m2)3-2m2sn2(x-2((1+m2)2-3m2)1+m2t),(3.31)u22=-2(1+m2)3-2m2cd2(x-2((1+m2)2-3m2)1+m2t),(3.32)u23=-2(2m2-1)3+2m2cn2(x+2((2m2-1)2-3m2(m2-1))2m2-1t),(3.33)u24=-2(2-m2)3+2dn2(x-2((2-m2)2+3(-1+m2))2-m2t),(3.34)u25=-2(2-m2)3-2(1-m2)sc2(x+2((2-m2)2-3(1-m2))2-m2t),(3.35)u26=2(2m2-1)3+2m2(1+m2)sd2(x+2(3m2(1-m2)+(-1+2m2)2)2m2-1t),(3.36)u27=2(1-2m2)6-12(ns(x+2((0.5-m2)2-3/16)0.5-m2t)+cs(x+2((0.5-m2)2-3/16)0.5-m2t))2,(3.37)u28=-1+m23+1-m22×(nc(x+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5t)+sc(x+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5t))2,(3.38)u29=-m2-23-m22(ns(x+2((0.5m2-1)2-m2/16)0.5m2-1t)+ds(x+2((0.5m2-1)2-m2/16)0.5m2-1t))2,(3.39)u30=m2-23-m22(sn(x+2((0.5m2-1)2-m4/16)0.5m2-1t)+ics(x+2((0.5m2-1)2-(m4/16))0.5m2-1t))2.


### 3.1. Soliton Solutions

Some solitary wave solutions can be obtained, if the modulusm approaches to 1 in (3.10)–(3.39) as follows:
(3.40)u31=-43-2tanh2(x-16t)-2coth2(x-16t),u32=-23+2sech2(x+2t),u33=-23+2sech2(x-2t),u34=-43-2coth2(x-t),u35=-23-2csch2(x+2t),u36=23+4sinh2(x+2t)-2csch2(x+2t),u37=-13-12(coth(x-4t)+csch(x-4t))2-12(coth(x-4t)+csch(x-4t))-2,u38=-13-12(tanh(x-4t)+icsch(x-4t))2+12(tanh(x-4t)+icsch(x-4t))-2,u39=-13-12(coth(x-14t)+csch(x-14t))-2,u40=-13+12(tanh(x-14t)+icsch(x-14t))-2,u41=-43-2tanh2(x-t),u42=23+4sinh2(x+2t),u43=-13-12(coth(x-14t)+csch(x-14t))2,u44=13-12(coth(x+2t)+csch(x+2t))2,u45=-13-12(tanh(x-34t)+icsch(x-34t))2.


### 3.2. Triangular Periodic Solutions

Some trigonometric function solutions can be obtained, if the modulusm approaches to zero in (3.10)–(3.39) as follwos:
(3.41)u46=-23-2csc2(x-2t),u47=-23-2sec2(x-2t),u48=23-2sec2(x-2t),u49=-43-2tan2(x+16t)-2cot2(x+16t),u50=13-12(csc(x+4t)+cot(x+4t))2-12(csc(x+4t)+cot(x+4t))-2,u51=-13+12(sec(x+7t)+tan(x+7t))2+12(sec(x+7t)+tan(x+7t))-2,u52=23-sin2(x-2t),u53=-43-2cot2(x+t),u54=13-12(csc(x+14t)+cot(x+14t))-2,u55=-13+12(sec(x-12t)+tan(x-12t))-2,u56=23-12sin2(x+38t),u57=-43-2tan2(x+t),u58=13-12(csc(x+14t)+cot(x+14t))2,u59=-13+12(sec(x-12t)+tan(x-12t))2.The modulus of solitary wave solutionsu1,  u2,  u21, and u23 is displayed in Figures 1, 2, 3, and 4, respectively, with values of parameters listed in their captions.The modulus of solitary wave solutionu1 (3.10) where m=0.5.
(a)
(b)The modulus of solitary wave solutionu2 (3.11) where m=0.5.
(a)
(b)The modulus of solitary wave solutionu21 (3.30) where m=0.5.
(a)
(b)The modulus of solitary wave solutionu23 (3.32) where m=0.5.
(a)
(b)


## 3.1. Soliton Solutions

Some solitary wave solutions can be obtained, if the modulusm approaches to 1 in (3.10)–(3.39) as follows:
(3.40)u31=-43-2tanh2(x-16t)-2coth2(x-16t),u32=-23+2sech2(x+2t),u33=-23+2sech2(x-2t),u34=-43-2coth2(x-t),u35=-23-2csch2(x+2t),u36=23+4sinh2(x+2t)-2csch2(x+2t),u37=-13-12(coth(x-4t)+csch(x-4t))2-12(coth(x-4t)+csch(x-4t))-2,u38=-13-12(tanh(x-4t)+icsch(x-4t))2+12(tanh(x-4t)+icsch(x-4t))-2,u39=-13-12(coth(x-14t)+csch(x-14t))-2,u40=-13+12(tanh(x-14t)+icsch(x-14t))-2,u41=-43-2tanh2(x-t),u42=23+4sinh2(x+2t),u43=-13-12(coth(x-14t)+csch(x-14t))2,u44=13-12(coth(x+2t)+csch(x+2t))2,u45=-13-12(tanh(x-34t)+icsch(x-34t))2.


## 3.2. Triangular Periodic Solutions

Some trigonometric function solutions can be obtained, if the modulusm approaches to zero in (3.10)–(3.39) as follwos:
(3.41)u46=-23-2csc2(x-2t),u47=-23-2sec2(x-2t),u48=23-2sec2(x-2t),u49=-43-2tan2(x+16t)-2cot2(x+16t),u50=13-12(csc(x+4t)+cot(x+4t))2-12(csc(x+4t)+cot(x+4t))-2,u51=-13+12(sec(x+7t)+tan(x+7t))2+12(sec(x+7t)+tan(x+7t))-2,u52=23-sin2(x-2t),u53=-43-2cot2(x+t),u54=13-12(csc(x+14t)+cot(x+14t))-2,u55=-13+12(sec(x-12t)+tan(x-12t))-2,u56=23-12sin2(x+38t),u57=-43-2tan2(x+t),u58=13-12(csc(x+14t)+cot(x+14t))2,u59=-13+12(sec(x-12t)+tan(x-12t))2.The modulus of solitary wave solutionsu1,  u2,  u21, and u23 is displayed in Figures 1, 2, 3, and 4, respectively, with values of parameters listed in their captions.The modulus of solitary wave solutionu1 (3.10) where m=0.5.
(a)
(b)The modulus of solitary wave solutionu2 (3.11) where m=0.5.
(a)
(b)The modulus of solitary wave solutionu21 (3.30) where m=0.5.
(a)
(b)The modulus of solitary wave solutionu23 (3.32) where m=0.5.
(a)
(b)


## 4. Generalized (2+1)-Dimensional Ito Equation

In this section we consider the generalized (2+1)-dimensional Ito equation (1.2) as follows:
(4.1)utt+uxxxt+3(2uxut+uuxt)+3uxx∫∞x‍utdx′+αuyt+βuxt=0,
if we use the transformation u=vx, it carries (4.1) into
(4.2)vxtt+vxxxxt+3(2vxxvxt+vxvxxt)+3vxxxvt+αvxyt+βvxxt=0,
if we use ζ=x+y-νt carries (4.2) into the ODE, we have
(4.3)(ν-α-β)V′′′-V(v)-3((V′)2)′′=0,
where by integrating twice we obtain, upon setting the constant of integration to zero,
(4.4)(ν-α-β)V′-V′′′-3(V′)2=0,
if we use the transformation W=V′, it carries (4.4) into
(4.5)(ν-α-β)W-W′′-3W2=0.
Balancing the term W′′ with the term W2, we obtain N=2, then
(4.6)W(ζ)=a0+a1ψ+a-1ψ-1+a2ψ2+a-2ψ-2,ψ′=A+Bψ2+Cψ4.Proceeding as in the previous case, we obtain(1)(4.7)a-1=a-1=0,a0=-2B3,a2=-2C,a-2=-2A,ν=α+β-2(B2+12AC)B,(2)(4.8)a-1=a2=a-1=0,a0=-2B3,a-2=-2A,ν=α+β-2(B2-3AC)B,(3)(4.9)a-1=a-2=a-1=0,a0=-2B3,a2=-2C,ν=α+β-2(B2-3AC)B.
The solutions of (4.1) are
(4.10)u1=-2(1+m2)3-2m2sn2(x+(α+β-2((1+m2)2+12m2)1+m2)t)-2ns2(x+y+(α+β-2((1+m2)2+12m2)1+m2)t),u2=-2(1+m2)3-2m2cd2(x+y+(α+β-2((1+m2)2+12m2)1+m2)t)-2dc2(x+y+(α+β-2((1+m2)2+12m2)1+m2)t),u3=-2(2m2-1)3+2m2cn2(x+y+(α+β+2(12m2(m2-1)+(2m2-1)2)2m2-1)t)-2(1-m2)nc2(x+y+(α+β+2(12m2(m2-1)+(2m2-1)2)2m2-1)t),u4=-2(2-m2)3+2dn2(x+y+(α+β-2((2-m2)2-12(-1+m2))2-m2)t)-2(m2-1)nd2(x+y+(α+β-2((2-m2)2-12(-1+m2))2-m2)t),u5=-2(1+m2)3-2ns2(x+y+(α+β-2((1+m2)2-3m2)1+m2)t),u6=-2(2-m2)3-2(1-m2)sc2(x+y+(α+β+2(12(1-m2)+(2-m2)2)2-m2)t)-2cs2(x+y+(α+β+2(12(1-m2)+(2-m2)2)2-m2)t),u7=2(2m2-1)3+2m2(1+m2)sd2(x+y+(α+β+2(-12m2(1-m2)+(-1+2m2)2)2m2-1)t)-2ds2(x+y+(α+β+2(-12m2(1-m2)+(-1+2m2)2)2m2-1)t),u8=2(1-2m2)6-12(ns(x+y+(α+β+2(0.75+(0.5-m2)2)0.5-m2)t)+cs(x+y+(α+β+2(0.75+(0.5-m2)2)0.5-m2)t))2-12(ns(x+y+(α+β+2(0.75+(0.5-m2)2)0.5-m2)t)+cs(x+y+(α+β+2(0.75+(0.5-m2)2)0.5-m2)t))-2,u9=-1+m23+1-m22×(nc(x+y+(α+β+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5)t)+sc(x+y+(α+β+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5)t))2+1-m22×(nc(x+y+(α+β+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5)t)+sc(x+y+(α+β+2(12(0.5-0.5m2)(0.25-0.25m2)+(0.5+0.5m2)2)0.5m2+0.5)t))-2,u10=-m2-23-m22(ns(x+y+(α+β+2(0.75m2+(-1+0.5m2)2)0.5m2-1)t)+ds(x+y+(α+β+2(0.75m2+(-1+0.5m2)2)0.5m2-1)t))2-12(ns(x+y+(α+β+2(0.75m2+(-1+0.5m2)2)0.5m2-1)t)+ds(x+y+(α+β+2(0.75m2+(-1+0.5m2)2)0.5m2-1)t))-2,u11=m2-23-m22(sn(x+y+(α+β+2(0.75m4+(-1+0.5m2)2)0.5m2-1)t)+ics(x+y+(α+β+2(0.75m4+(-1+0.5m2)2)0.5m2-1)t))2+m22(sn(x+y+(α+β+2(0.75m4+(-1+0.5m2)2)0.5m2-1)t)+ics(x+y+(α+β+2(0.75m4+(-1+0.5m2)2)0.5m2-1)t))-2,u12=-2(1+m2)3-2dc2(x+y+(α+β-2((1+m2)2-3m2)1+m2)t),u13=-2(2m2-1)3-2(1-m2)nc2(x+y+(α+β+2(-3m2(m2-1)+(2m2-1)2)2m2-1)t),u14=-2(2-m2)3-2(m2-1)nd2(x+y+(α+β-2((2-m2)2+3(-1+m2))2-m2)t),u15=-2(2-m2)3-2cs2(x+y+(α+β+2((2-m2)2-3(1-m2))2-m2)t),u16=2(2m2-1)3-2ds2(x+y+(α+β+2(3m2(1-m2)+(-1+2m2)2)2m2-1)t),u17=2(1-2m2)6-12(ns(x+y+(α+β+2((0.5-m2)2-3/16)0.5-m2)t)+cs(x+y+(α+β+2((0.5-m2)2-3/16)0.5-m2)t))-2,u18=-1+m23+1-m22×(nc(x+y+(α+β+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5)t)+sc(x+y+(α+β+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5)t))-2,u19=-m2-23-12(ns(x+y+(α+β+2((0.5m2-1)2-3/16)0.5m2-1)t)+ds(x+y+(α+β+2((0.5m2-1)2-3/16)0.5m2-1)t))-2,u20=m2-23+m22×(sn(x+y+(α+β+2((0.5m2-1)2-m4/16)0.5m2-1)t)+ics(x+y+(α+β+2((0.5m2-1)2-m4/16)0.5m2-1)t))-2,u21=-2(1+m2)3-2m2sn2(x+y+(α+β-2((1+m2)2-3m2)1+m2)t),u22=-2(1+m2)3-2m2cd2(x+y+(α+β-2((1+m2)2-3m2)1+m2)t),u23=-2(2m2-1)3+2m2cn2(x+y+(α+β+2((2m2-1)2-3m2(m2-1))2m2-1)t),u24=-2(2-m2)3+2dn2(x+y+(α+β-2((2-m2)2+3(-1+m2))2-m2)t),u25=-2(2-m2)3-2(1-m2)sc2(x+y+(α+β+2((2-m2)2-3(1-m2))2-m2)t),u26=2(2m2-1)3+2m2(1+m2)sd2(x+y+(α+β+2(3m2(1-m2)+(-1+2m2)2)2m2-1)t),u27=2(1-2m2)6-12(ns(x+y+(α+β+2((0.5-m2)2-3/16)0.5-m2)t)+cs(x+y+(α+β+2((0.5-m2)2-3/16)0.5-m2)t))2,u28=-1+m23+1-m22×(nc(x+y+(α+β+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5)t)+sc(x+y+(α+β+2((0.5+0.5m2)2-3(0.5-0.5m2)(0.25-0.25m2))0.5m2+0.5)t))2,u29=-m2-23-m22×(ns(x+y+(α+β+2((0.5m2-1)2-m2/16)0.5m2-1)t)+ds(x+y+(α+β+2((0.5m2-1)2-m2/16)0.5m2-1)t))2,u30=m2-23-m22×(sn(x+y+(α+β+2((0.5m2-1)2-m4/16)0.5m2-1)t)+ics(x+y+(α+β+2((0.5m2-1)2-m4/16)0.5m2-1)t))2.


### 4.1. Soliton Solutions

Some solitary wave solutions can be obtained, if the modulusm approaches to 1 in (4.10) as follows:
(4.11)u31=-43-2tanh2(x+y+(α+β-16)t)-2coth2(x+y+(α+β-16)t),u32=-23+2sech2(x+y+(α+β+2)t),u33=-23+2sech2(x+y+(α+β-2)t),u34=-43-2coth2(x+y+(α+β-1)t),u35=-23-2csch2(x+y+(α+β+2)t),u36=23+4sinh2(x+y+(α+β+2)t)-2csch2(x+y+(α+β+2)t),u37=-13-12(coth(x+y+(α+β-4)t)+csch(x+y+(α+β-4)t))2-12(coth(x+y+(α+β-4)t)+csch(x+y+(α+β-4)t))-2,u38=-13-12(tanh(x+y+(α+β-4)t)+icsch(x+y+(α+β-4)t))2+12(tanh(x+y+(α+β-4)t)+icsch(x+y+(α+β-4)t))-2,u39=-13-12(coth(x+y+(α+β-14)t)+csch(x+y+(α+β-14)t))-2,u40=-13+12(tanh(x+y+(α+β-14)t)+icsch(x+y+(α+β-14)t))-2,u41=-43-2tanh2(x+y+(α+β-1)t),u42=23+4sinh2(x+y+(α+β+2)t),u43=-13-12(coth(x+y+(α+β-14)t)+csch(x+y+(α+β-14)t))2,u44=13-12(coth(x+y+(α+β+2)t)+csch(x+y+(α+β+2)t))2,u45=-13-12(tanh(x+y+(α+β-34)t)+icsch(x+y+(α+β-34)t))2.


### 4.2. Triangular Periodic Solutions

Some trigonometric function solutions can be obtained, if the modulusm approaches to zero in (4.10) as follows:
(4.12)u46=-23-2csc2(x+y+(α+β-2)t),u47=-23-2sec2(x+y+(α+β-2)t),u48=23-2sec2(x+y+(α+β-2)t),u49=-43-2tan2(x+y+(α+β+16)t)-2cot2(x+y+(α+β+16)t),u50=13-12(csc(x+y+(α+β+4)t)+cot(x+y+(α+β+4)t))2-12(csc(x+y+(α+β+4)t)+cot(x+y+(α+β+4)t))-2,u51=-13+12(sec(x+y+(α+β+7)t)+tan(x+y+(α+β+7)t))2+12(sec(x+y+(α+β+7)t)+tan(x+y+(α+β+7)t))-2,u52=23-sin2(x+y+(α+β-2)t),u53=-43-2cot2(x+y+(α+β+1)t),u54=13-12(csc(x+y+(α+β+14)t)+cot(x+y+(α+β+14)t))-2,u55=-13+12(sec(x+y+(α+β-12)t)+tan(x+y+(α+β-12)t))-2,u56=23-12sin2(x+y+(α+β+38)t),u57=-43-2tan2(x+y+(α+β+1)t),u58=13-12(csc(x+y+(α+β+14)t)+cot(x+y+(α+β+14)t))2,u59=-13+12(sec(x+y+(α+β-12)t)+tan(x+y+(α+β-12)t))2.If we takem→1, in the two Sections 3 and 4, we obtain the solutions degenerated by the hyperbolic extended hyperbolic functions methods (tanh, coth, sinh, sech,…, etc.) (see, for example [42]). Moreover, when m→0, the solutions obtained by triangular and extended triangular functions methods (tan, sine, cosine, sec,…, etc.) are found as disused in Sections 3.1, 3.2, 4.1, and 4.2.


## 4.1. Soliton Solutions

Some solitary wave solutions can be obtained, if the modulusm approaches to 1 in (4.10) as follows:
(4.11)u31=-43-2tanh2(x+y+(α+β-16)t)-2coth2(x+y+(α+β-16)t),u32=-23+2sech2(x+y+(α+β+2)t),u33=-23+2sech2(x+y+(α+β-2)t),u34=-43-2coth2(x+y+(α+β-1)t),u35=-23-2csch2(x+y+(α+β+2)t),u36=23+4sinh2(x+y+(α+β+2)t)-2csch2(x+y+(α+β+2)t),u37=-13-12(coth(x+y+(α+β-4)t)+csch(x+y+(α+β-4)t))2-12(coth(x+y+(α+β-4)t)+csch(x+y+(α+β-4)t))-2,u38=-13-12(tanh(x+y+(α+β-4)t)+icsch(x+y+(α+β-4)t))2+12(tanh(x+y+(α+β-4)t)+icsch(x+y+(α+β-4)t))-2,u39=-13-12(coth(x+y+(α+β-14)t)+csch(x+y+(α+β-14)t))-2,u40=-13+12(tanh(x+y+(α+β-14)t)+icsch(x+y+(α+β-14)t))-2,u41=-43-2tanh2(x+y+(α+β-1)t),u42=23+4sinh2(x+y+(α+β+2)t),u43=-13-12(coth(x+y+(α+β-14)t)+csch(x+y+(α+β-14)t))2,u44=13-12(coth(x+y+(α+β+2)t)+csch(x+y+(α+β+2)t))2,u45=-13-12(tanh(x+y+(α+β-34)t)+icsch(x+y+(α+β-34)t))2.


## 4.2. Triangular Periodic Solutions

Some trigonometric function solutions can be obtained, if the modulusm approaches to zero in (4.10) as follows:
(4.12)u46=-23-2csc2(x+y+(α+β-2)t),u47=-23-2sec2(x+y+(α+β-2)t),u48=23-2sec2(x+y+(α+β-2)t),u49=-43-2tan2(x+y+(α+β+16)t)-2cot2(x+y+(α+β+16)t),u50=13-12(csc(x+y+(α+β+4)t)+cot(x+y+(α+β+4)t))2-12(csc(x+y+(α+β+4)t)+cot(x+y+(α+β+4)t))-2,u51=-13+12(sec(x+y+(α+β+7)t)+tan(x+y+(α+β+7)t))2+12(sec(x+y+(α+β+7)t)+tan(x+y+(α+β+7)t))-2,u52=23-sin2(x+y+(α+β-2)t),u53=-43-2cot2(x+y+(α+β+1)t),u54=13-12(csc(x+y+(α+β+14)t)+cot(x+y+(α+β+14)t))-2,u55=-13+12(sec(x+y+(α+β-12)t)+tan(x+y+(α+β-12)t))-2,u56=23-12sin2(x+y+(α+β+38)t),u57=-43-2tan2(x+y+(α+β+1)t),u58=13-12(csc(x+y+(α+β+14)t)+cot(x+y+(α+β+14)t))2,u59=-13+12(sec(x+y+(α+β-12)t)+tan(x+y+(α+β-12)t))2.If we takem→1, in the two Sections 3 and 4, we obtain the solutions degenerated by the hyperbolic extended hyperbolic functions methods (tanh, coth, sinh, sech,…, etc.) (see, for example [42]). Moreover, when m→0, the solutions obtained by triangular and extended triangular functions methods (tan, sine, cosine, sec,…, etc.) are found as disused in Sections 3.1, 3.2, 4.1, and 4.2.


## 5. Conclusion

By introducing appropriate transformations and using extended F-expansion method, we have been able to obtain, in a unified way with the aid of symbolic computation system-mathematica, a series of solutions including single and the combined Jacobi elliptic function. Also, extended F-expansion method showed that soliton solutions and triangular periodic solutions can be established as the limits of Jacobi doubly periodic wave solutions. Whenm→1, the Jacobi functions degenerate to the hyperbolic functions and give the solutions by the extended hyperbolic functions methods. When m→0, the Jacobi functions degenerate to the triangular functions and give the solutions by extended triangular functions methods. In fact, the disadvantage of extended F-expansion method is the existence of complex solutions which are listed here just as solutions.



---

*Source: 537930-2012-11-29.xml*