# Infinitely Many Elliptic Solutions to a Simple Equation and Applications

Authors: Long Wei; Yang Wang
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

Based on auxiliary equation method and Bäcklund transformations, we present an idea to find
infinitely many Weierstrass and Jacobi elliptic function solutions to some nonlinear problems.
First, we give some nonlinear iterated formulae of solutions and some elliptic function solutions
to a simple auxiliary equation, which results in infinitely many Weierstrass and Jacobi elliptic
function solutions of the simple equation. Then applying auxiliary equation method to some
nonlinear problems and combining the results with exact solutions of the auxiliary equation,
we obtain infinitely many elliptic function solutions to the corresponding nonlinear problems.
The employed approach is powerful and can be also applied to solve other nonlinear differential
equations.


---

## Body

## 1. Introduction

Partial differential equations (PDEs) describe various nonlinear phenomena in natural and applied sciences such as fluid dynamics, plasma physics, solid state physics, optical fibers, acoustics, mechanics, biology, and mathematical finance. Their solution spaces are infinite dimensional and contain diverse solution structures. It is of significant importance to solve nonlinear PDEs from both theoretical and practical points of view. In some decades, many authors are devoted to investigate the exact solutions to nonlinear PDEs. With the help of exact solutions, when they exist, the mechanism of complicated physical phenomena and dynamical processes modeled by these nonlinear PDEs can be better understood. They can also help to analyze the stability of these solutions and to check numerical analysis for these nonlinear PDEs.In recent years, reducing PDEs into ordinary differential equations (ODEs) has been proved to be a successful idea to generate exact solutions of nonlinear wave equations. Many approaches to exact solutions in the literature follow such an idea, which contain tanh-function method [1–4], sech-function method [5, 6], homogeneous balance method [7, 8], extended tanh-function method [9], Jacobi elliptic function method [10], exp-function method [11], F-expansion method [12], and transformed rational function method [13]. The tanh-function method and the G′/G expansion methods are special cases of general Frobenius’ idea [14] or the transformed rational function method [13]. Recently, the multiple exp-function was also used to present three waves in (3 + 1) dimension [15]. There is also a generalized theory of the Bell polynomials [16], which describe generalized bilinear differential equations [17]. Among those, auxiliary equation methods are widely used to find exact solutions (especially elliptic solutions) of PDEs by some authors (see [1, 3, 18–24], etc.). But in many literatures, the authors constructed exact solutions of PDEs making use of finite solutions of the auxiliary equations employed and have not obtained infinitely many solutions, especially Weierstrass and Jacobi elliptic solutions.In this work, we will seek for nonlinear iterated formulae of solutions to a simple auxiliary equation. Simultaneously, in order to obtain infinitely many new elliptic solutions of some nonlinear PDEs, we will also give some Weierstrass and Jacobi elliptic function solutions of the simple auxiliary equation. Based on these facts, we can construct infinitely many Weierstrass and Jacobi elliptic function solutions of destination equations. The obtained results are new and introduce useful analysis for equations of identical nonlinearities.This paper is organized as follows: in Section2, we provide some nonlinear iterated formulae of solutions and some new elliptic solutions to a simple equation. In Section 3, we give some examples to demonstrate the application of our approach. The remainder of this paper is a discussion.


## 2. Preliminaries

In this section, we try to find some nonlinear iterated formulae of solutions and infinitely many Weierstrass and Jacobi elliptic function solutions to the auxiliary equation(1)(dη(ξ)dξ)2=aη+bη2+cη3.
We will solve some NPDEs by the exact solutions of the auxiliary equation; consequently we can obtain infinitely many Weierstrass and Jacobi elliptic function solutions and periodic and soliton solutions of destination equations.


### 2.1. Nonlinear Iterated Formulae of Solutions to (1)

In order to find infinitely many elliptic solutions to (1), we assume that abc≠0.By the Bäcklund transformations and a modified truncation approach, we can obtain the following nonlinear iterated formulae of solutions; that is, ifw=w(ξ) is a solution to (1), then so is η=η(ξ):
(2)η1=acw;(3)η2±=±2a+(b∓b2-4ac)w∓b+b2-4ac∓2cw;(4)η3±=±∓ab+ab2-4ac+(∓b2+bb2-4ac±2ac)wc[2a+(b∓b2-4ac)w];(5)η4±=±([∓2a+(b2-4ac∓b)w]×[(2ac-b2)(b2-4ac±b)±2abc])×(2c([(2ac-b2)(b2-4ac±b)±2abc]w-ab2-4ac(b±b2-4ac)))-1;(6)η5±=-ab±ab2-4ac-4acw-c(b±b2-4ac)w22c(a+bw+cw2);(7)ηi±=A0,i+A1,iw+A2,iw2B0,i+B1,iw+B2,iw2,(i=6,7,8,9),
where
(8)A0,6=∓R6±,A1,6=1∓baR6±,A2,6=-ab(b2-3ac)∓(2a2c2-4ab2c+b4)R6±a2c(2a∓bR6±),B0,6=1,B1,6=ab∓(b2-2ac)R6±a(2a±bR6±),B2,6=a(2a2c2+b4-4ab2c)∓b(5a2c2-5ab2c+b4)R6±a2[a(b2-2ac)∓b(b2-3ac)R6±],R6±=±b+b2-4ac2c;A0,7=±(c+Ma)R7±,A1,7=1,A2,7=aM[a(b2-2ac)±bc(3ac-b2)R7±]ab(3ac-b2)±c(2a2c2-4ab2c+b4)R7±,B0,7=-1-acM,B1,7=-2a2M-ac±bc(aM+c)R7±ac2R7±,B2,7=M,R7±=±b+b2-4ca2c2,Misanarbitraryconstant;A0,8=M,A1,8=∓R8±,A2,8=1∓baR8±,B0,8=c(2a∓bR8±)M-ab±(b2-2ac)R8±,B1,8=1,B2,8=c(2a∓bR8±)ab2-4ac,R8±=±b+b2-4ac2c,Misanarbitraryconstant;A0,9=M,A1,9=∓NR9±,A2,9=-a(b-2cN)+±c(bN-2a)R9±b2-2ac-bcN±c(2cN-b)R9±,B0,9=c(2a∓bR9±)M-ab±(b2-2ac)R9±,B1,9=N,B2,9=1,R9±=±b+b2-4ac2c,M,Naretwoarbitraryconstants;(9)η10±=±a(b∓2acw+w′)ac(b∓2acw-w′);(10)η11±=-A1+A2w±A3Lw′c(b+2b2-3ac)(4a2c-ab2±3Lw′),
with
(11)A1=18a2b2c-3ab4-24a3c2+3ab(4ac-b2)b2-3ac,A2=2ac(4ac-b2)(b+2b2-3ac),A3=6ac-b2+bb2-3ac,L=13(b2-4ac)[9abc-2b3+2(3ac-b2)b2-3ac];(12)η12±=-A1+A2w±A3Lw′c(-b+2b2-3ac)(4a2c-ab2±3Lw′),
with
(13)A1=-18a2b2c+3ab4+24a3c2-3ab(4ac-b2)b2-3ac,A2=2ac(4ac-b2)(-b+2b2-3ac),A3=-6ac+b2+bb2-3ac,L=13(b2-4ac)[9abc-2b3-2(3ac-b2)b2-3ac];(14)η13±=A1+A2w±A3w′B1±B2w′,
with
(15)A1=-6a[(-b+b2-3ac)(4b4-27ab2c+36a2c2)+6abc(b2-6ac)(-b+b2-3ac)],A2=27ML2,A3=9c(-b+b2-3ac)ML,B1=8(-b+b2-3ac)(27a2bc2-15ab3c+2b5)+6ac(4b4-27ab2c+36a2c2),B2=27ML,L=2cb(2b2-9ac)-2(b2-3ac)3/2,M=c2[(-b+b2-3ac)(6ac-2b2)-3abc];(16)η14±=A1+A2w±A3w′B1±B2w′,
with
(17)A1=6a[(b+b2-3ac)(4b4-27ab2c+36a2c2)-6abc(b2-6ac)(b+b2-3ac)],A2=27ML2,A3=-9c(b+b2-3ac)ML,B1=-8(b+b2-3ac)(27a2bc2-15ab3c+2b5)+6ac(4b4-27ab2c+36a2c2),B2=27ML,L=2cb(2b2-9ac)+2(b2-3ac)3/2,M=-c2[(b+b2-3ac)(6ac-2b2)+3abc];(18)η15=ac(a+bw+cw2+b-2acw′)c(-a-bw-cw2+b-2acw′);(19)η16=ac(a+bw+cw2-b+2acw′)c(a+bw+cw2+b+2acw′).Above, we have given finitely many nonlinear iterated formulae of solutions to (1). In fact, it is possible to obtain infinitely many ones by substituting one into another repeatedly. Here, it should be pointed out that the previous iterative procedure can be done finite times for some formulae but infinite times for some other ones. For example, substituting (3) into (2) we can obtain a new nonlinear iterated formula:
(20)η2-1±=±a[∓b+b2-4ac∓2cw]c[2a+(b∓b2-4ac)w],
but repeating the procedure again will result in (3) itself; substituting η2- into η2-1- results in a/cw; substituting η2+ into η2-1- results in w and so on. But for solutions η10± and so forth, the iterative procedure are infinite For instance, substituting η2± or η5± into η10± and iterating the results by η10± can generate infinitely many iterated formulae of solutions.


### 2.2. Infinitely Many Weierstrass Elliptic Function Solutions to (1)

Let℘(ξ)=℘(ξ;g2,g3) be the Weierstrass elliptic function satisfying
(21)℘′2(ξ)=4℘3(ξ)-g2℘(ξ)-g3,
where g2, g3 are real parameters and called invariants [25], and ℘′(ξ) denotes the first-order derivative. It follows from (21) that the second-order derivative is
(22)℘′′(ξ)=6℘2(ξ)-g2.
From the last previous two equations, we will give some Weierstrass elliptic function solutions to (1), among which there are many new solutions different from those in [26] and some latter references. In fact, in [26] the author considered the following auxiliary equation:
(23)(dθdξ)2=d0+d2θ2+d4θ4,
but this equation can be converted into (1) by the transformation η(ξ)=θ2(ξ), with a=4d0, b=4d2, and c=4d4.From (21) and (22), we have the following Weierstrass elliptic function solutions to (1).(i) Wheng2=(b2-3ac)/12, g3=b(9ac-2b2)/432, then
(24)η=3a12℘(ξ)-b,12℘(ξ)-b3c;
or(25)η=L1℘(ξ)-(b(5ac-b2)L1+a(6ac-b2))/12c(2a+bL1)℘(ξ)-(2(3ac-b2)L1-ab)/12(2a+bL1),a(℘(ξ)-(2(3ac-b2)L1-ab)/12(2a+bL1))c(L1℘(ξ)-(b(5ac-b2)L1+a(6ac-b2))/12c(2a+bL1)),where L1=(-b±b2-4ac)/2c.(ii) Wheng2=-11ac/64-(b2+5bL2)/48, g3=5b3/3456-abc/512+7(b2-3ac)L2/1152, then
(26)η=3a(12℘(ξ)+b+2L2)(12℘(ξ)+L2)2,(12℘(ξ)+L2)23c(12℘(ξ)+b+2L2),
where L2=(-5b±3b2-4ac)/8.(iii) Wheng2=b2/192+ac/16 and g3=b(36ac-b2)/13824, then
(27)η=576℘′2(ξ)c(24℘(ξ)+b)2,a(24℘(ξ)+b)2576℘′2(ξ),(28)b2-36ac-96b℘(ξ)+2304℘2(ξ)24c(24℘(ξ)+b),24a(24℘(ξ)+b)b2-36ac-96b℘(ξ)+2304℘2(ξ),(29)(L3℘(ξ)-A-1)2L3(℘(ξ)+B-1)2,aL3(℘(ξ)+B-1)2c(L3℘(ξ)-A-1)2
with A-1=(4ab+(12ac-b2)L3)/48(b+2cL3), B-1=(4ab-(12ac-5b2)L3)/48(2a+bL3), and L3=(-b±b2-4ac)/2c;
(30)η=L3(℘(ξ)-A-2)2(℘(ξ)-B-2)2,a(℘(ξ)-B-2)2cL3(℘(ξ)-A-2)2
with A-2=(a(12ac-b2)+b(8ac-b2)L3)/48[(2ac-b2)L3-ab], B-2=(ab(5b2-16ac)+(5b4+12a2c2-21ab2c)L3)/48[a(2ac-b2)+b(3ac-b2)L3], and L3=(-b±b2-4ac)/2c; or
(31)η=(A-3℘(ξ)+B-3+E℘′(ξ))2(A-4℘(ξ)+B-4-D℘′(ξ))2,a(A-4℘(ξ)+B-4-D℘′(ξ))2c(A-3℘(ξ)+B-3+E℘′(ξ))2
with A-3=D(2aD2+bE2)L-, B-3=(D/48)(12acE2-b2E2+4abD2)L-, A-4=E(bD2+2cE2)L-, and B-4=(E/48)(12acD2-b2D2+4bcE2)L-, where L-=±(1/2cE4+bE2D2+aD4), and D,E are two arbitrary constants.Note that in any pair of solutions, the latter can be obtained from the former through formula (2). Indeed, combining solutions (24)–(31) with the nonlinear iterated formulae (2)–(20) and their infinitely many derivatives can yield infinitely many Weierstrass elliptic function solutions of (1).


### 2.3. Infinitely Many Jacobi Elliptic Function Solutions to (1)

Assume thatK, C, and D are arbitrary constants, and m∈(0,1) is the elliptic modulus; then (1) has the following Jacobi elliptic solutions.(I)a=4K,  b=-4(1+m2),  c=4m2/K,
(32)η=Kcd(ξ,m)2,Ksn(ξ,m)2,Km2dc(ξ,m)2,Km2ns(ξ,m)2.
Indeed, the former two solutions and the latter two ones can transform each other by formula (2). In what follows, we omit the solutions that can be generated by formula (2).(II)a=4K(1-m2),  b=-4(1-2m2),  c=-4m2/K,
(33)η=Kcn(ξ,m)2,K(1-m2)sd(ξ,m)2.(III)a=-4K(1-m2),  b=4(2-m2),  c=-4/K,
(34)η=Kdn(ξ,m)2,-Kcs(ξ,m)2.(IV)a=K(m2-1),  b=2(1+m2),  c=(m2-1)/K,
(35)η±=K[nd(ξ,m)±msd(ξ,m)]2,-K[nc(ξ,m)±(ξ,m)]2,K1-m2[dn(ξ,m)±mcn(ξ,m)]2,Km2-1[cs(ξ,m)±ds(ξ,m)]2.(V)a=((C2-D2)/K)(m2-1),  b=2(1+m2),  c=-(K/(C2-D2))(m2-1),
(36)η±=[Dcn(ξ,m)+Cdn(ξ,m)]2K[1±1+C2(1-m2)/(D2-C2)sn(ξ,m)]2.(VI)a=K,  b=2(-2+m2),  c=m4/K,
(37)η±=Kcn(ξ,m)2[1-m2±dn(ξ,m)]2,Ksn(ξ,m)2[1±dn(ξ,m)]2,Km2[cd(ξ,m)±i1-m2sd(ξ,m)]2,-Km2[cn(ξ,m)±isn(ξ,m)]2.(VII)a=Km2/(C2-D2),  b=2(-2+m2),  c=(C2-D2)m2/K,
(38)η±=-(K[C+Ddn(ξ,m)]2)×([cn(ξ,m)±D2m2D2-C2-1sn(ξ,m)]2(C2-D2)2×[cn(ξ,m)±D2m2D2-C2-1sn(ξ,m)]2)-1.(VIII)a=K,  b=2(1-2m2),  c=1/K,
(39)η±=K[ns(ξ,m)±cs(ξ,m)]2,K[msn(ξ,m)±idn(ξ,m)]2,K[dc(ξ,m)±1-m2(ξ,m)]2,K[mcd(ξ,m)±i1-m2nd(ξ,m)]2.(IX)a=Km2/(C2-D2),  b=2(1-2m2),  c=(C2-D2)m2/K,
(40)η±=-K[dn(ξ,m)±D2/(D2-C2)-m2sn(ξ,m)]2[C+Dcn(ξ,m)]2.(X)a=K(1-m)2,  b=2(1+6m+m2),  c=(1-m)2/K,
(41)η±=-K[1±msn(ξ,m)]2[1∓msn(ξ,m)]2.(XI)a±=±16Km,  b±=-4(1±6m+m2),  c±=(1±m)2/K,
(42)η±=K[1∓msn(ξ,m)2]2cn(ξ,m)2dn(ξ,m)2,K[ns(ξ,m)±msn(ξ,m)]2(1-m)2.(XII)a±=±16K1-m2,  b±=4(2-m2±61-m2),  c±=4(2-m2±21-m2)/K,
(43)η±=-K[cn(ξ,m)2±1-m2sn(ξ,m)2]2dn(ξ,m)2,-K[dn(ξ,m)±1-m2nd(ξ,m)]22-m2±21-m2,Km4sn(ξ,m)2×[(1-m2∓1-m2)nc(ξ,m)+m2cn(ξ,m)]2,Km4sn(ξ,m)2cn(ξ,m)2[dn(ξ,m)2∓1-m2]2.Combining solutions (32)–(43) with the nonlinear iterated formulae (2)–(20) and their derivatives can give infinitely many Jacobi elliptic function solutions of (1), and these solutions can be degenerated to solutions with hyperbolic or trigonometric functions as long as letting m go to 1 or 0.


## 2.1. Nonlinear Iterated Formulae of Solutions to (1)

In order to find infinitely many elliptic solutions to (1), we assume that abc≠0.By the Bäcklund transformations and a modified truncation approach, we can obtain the following nonlinear iterated formulae of solutions; that is, ifw=w(ξ) is a solution to (1), then so is η=η(ξ):
(2)η1=acw;(3)η2±=±2a+(b∓b2-4ac)w∓b+b2-4ac∓2cw;(4)η3±=±∓ab+ab2-4ac+(∓b2+bb2-4ac±2ac)wc[2a+(b∓b2-4ac)w];(5)η4±=±([∓2a+(b2-4ac∓b)w]×[(2ac-b2)(b2-4ac±b)±2abc])×(2c([(2ac-b2)(b2-4ac±b)±2abc]w-ab2-4ac(b±b2-4ac)))-1;(6)η5±=-ab±ab2-4ac-4acw-c(b±b2-4ac)w22c(a+bw+cw2);(7)ηi±=A0,i+A1,iw+A2,iw2B0,i+B1,iw+B2,iw2,(i=6,7,8,9),
where
(8)A0,6=∓R6±,A1,6=1∓baR6±,A2,6=-ab(b2-3ac)∓(2a2c2-4ab2c+b4)R6±a2c(2a∓bR6±),B0,6=1,B1,6=ab∓(b2-2ac)R6±a(2a±bR6±),B2,6=a(2a2c2+b4-4ab2c)∓b(5a2c2-5ab2c+b4)R6±a2[a(b2-2ac)∓b(b2-3ac)R6±],R6±=±b+b2-4ac2c;A0,7=±(c+Ma)R7±,A1,7=1,A2,7=aM[a(b2-2ac)±bc(3ac-b2)R7±]ab(3ac-b2)±c(2a2c2-4ab2c+b4)R7±,B0,7=-1-acM,B1,7=-2a2M-ac±bc(aM+c)R7±ac2R7±,B2,7=M,R7±=±b+b2-4ca2c2,Misanarbitraryconstant;A0,8=M,A1,8=∓R8±,A2,8=1∓baR8±,B0,8=c(2a∓bR8±)M-ab±(b2-2ac)R8±,B1,8=1,B2,8=c(2a∓bR8±)ab2-4ac,R8±=±b+b2-4ac2c,Misanarbitraryconstant;A0,9=M,A1,9=∓NR9±,A2,9=-a(b-2cN)+±c(bN-2a)R9±b2-2ac-bcN±c(2cN-b)R9±,B0,9=c(2a∓bR9±)M-ab±(b2-2ac)R9±,B1,9=N,B2,9=1,R9±=±b+b2-4ac2c,M,Naretwoarbitraryconstants;(9)η10±=±a(b∓2acw+w′)ac(b∓2acw-w′);(10)η11±=-A1+A2w±A3Lw′c(b+2b2-3ac)(4a2c-ab2±3Lw′),
with
(11)A1=18a2b2c-3ab4-24a3c2+3ab(4ac-b2)b2-3ac,A2=2ac(4ac-b2)(b+2b2-3ac),A3=6ac-b2+bb2-3ac,L=13(b2-4ac)[9abc-2b3+2(3ac-b2)b2-3ac];(12)η12±=-A1+A2w±A3Lw′c(-b+2b2-3ac)(4a2c-ab2±3Lw′),
with
(13)A1=-18a2b2c+3ab4+24a3c2-3ab(4ac-b2)b2-3ac,A2=2ac(4ac-b2)(-b+2b2-3ac),A3=-6ac+b2+bb2-3ac,L=13(b2-4ac)[9abc-2b3-2(3ac-b2)b2-3ac];(14)η13±=A1+A2w±A3w′B1±B2w′,
with
(15)A1=-6a[(-b+b2-3ac)(4b4-27ab2c+36a2c2)+6abc(b2-6ac)(-b+b2-3ac)],A2=27ML2,A3=9c(-b+b2-3ac)ML,B1=8(-b+b2-3ac)(27a2bc2-15ab3c+2b5)+6ac(4b4-27ab2c+36a2c2),B2=27ML,L=2cb(2b2-9ac)-2(b2-3ac)3/2,M=c2[(-b+b2-3ac)(6ac-2b2)-3abc];(16)η14±=A1+A2w±A3w′B1±B2w′,
with
(17)A1=6a[(b+b2-3ac)(4b4-27ab2c+36a2c2)-6abc(b2-6ac)(b+b2-3ac)],A2=27ML2,A3=-9c(b+b2-3ac)ML,B1=-8(b+b2-3ac)(27a2bc2-15ab3c+2b5)+6ac(4b4-27ab2c+36a2c2),B2=27ML,L=2cb(2b2-9ac)+2(b2-3ac)3/2,M=-c2[(b+b2-3ac)(6ac-2b2)+3abc];(18)η15=ac(a+bw+cw2+b-2acw′)c(-a-bw-cw2+b-2acw′);(19)η16=ac(a+bw+cw2-b+2acw′)c(a+bw+cw2+b+2acw′).Above, we have given finitely many nonlinear iterated formulae of solutions to (1). In fact, it is possible to obtain infinitely many ones by substituting one into another repeatedly. Here, it should be pointed out that the previous iterative procedure can be done finite times for some formulae but infinite times for some other ones. For example, substituting (3) into (2) we can obtain a new nonlinear iterated formula:
(20)η2-1±=±a[∓b+b2-4ac∓2cw]c[2a+(b∓b2-4ac)w],
but repeating the procedure again will result in (3) itself; substituting η2- into η2-1- results in a/cw; substituting η2+ into η2-1- results in w and so on. But for solutions η10± and so forth, the iterative procedure are infinite For instance, substituting η2± or η5± into η10± and iterating the results by η10± can generate infinitely many iterated formulae of solutions.


## 2.2. Infinitely Many Weierstrass Elliptic Function Solutions to (1)

Let℘(ξ)=℘(ξ;g2,g3) be the Weierstrass elliptic function satisfying
(21)℘′2(ξ)=4℘3(ξ)-g2℘(ξ)-g3,
where g2, g3 are real parameters and called invariants [25], and ℘′(ξ) denotes the first-order derivative. It follows from (21) that the second-order derivative is
(22)℘′′(ξ)=6℘2(ξ)-g2.
From the last previous two equations, we will give some Weierstrass elliptic function solutions to (1), among which there are many new solutions different from those in [26] and some latter references. In fact, in [26] the author considered the following auxiliary equation:
(23)(dθdξ)2=d0+d2θ2+d4θ4,
but this equation can be converted into (1) by the transformation η(ξ)=θ2(ξ), with a=4d0, b=4d2, and c=4d4.From (21) and (22), we have the following Weierstrass elliptic function solutions to (1).(i) Wheng2=(b2-3ac)/12, g3=b(9ac-2b2)/432, then
(24)η=3a12℘(ξ)-b,12℘(ξ)-b3c;
or(25)η=L1℘(ξ)-(b(5ac-b2)L1+a(6ac-b2))/12c(2a+bL1)℘(ξ)-(2(3ac-b2)L1-ab)/12(2a+bL1),a(℘(ξ)-(2(3ac-b2)L1-ab)/12(2a+bL1))c(L1℘(ξ)-(b(5ac-b2)L1+a(6ac-b2))/12c(2a+bL1)),where L1=(-b±b2-4ac)/2c.(ii) Wheng2=-11ac/64-(b2+5bL2)/48, g3=5b3/3456-abc/512+7(b2-3ac)L2/1152, then
(26)η=3a(12℘(ξ)+b+2L2)(12℘(ξ)+L2)2,(12℘(ξ)+L2)23c(12℘(ξ)+b+2L2),
where L2=(-5b±3b2-4ac)/8.(iii) Wheng2=b2/192+ac/16 and g3=b(36ac-b2)/13824, then
(27)η=576℘′2(ξ)c(24℘(ξ)+b)2,a(24℘(ξ)+b)2576℘′2(ξ),(28)b2-36ac-96b℘(ξ)+2304℘2(ξ)24c(24℘(ξ)+b),24a(24℘(ξ)+b)b2-36ac-96b℘(ξ)+2304℘2(ξ),(29)(L3℘(ξ)-A-1)2L3(℘(ξ)+B-1)2,aL3(℘(ξ)+B-1)2c(L3℘(ξ)-A-1)2
with A-1=(4ab+(12ac-b2)L3)/48(b+2cL3), B-1=(4ab-(12ac-5b2)L3)/48(2a+bL3), and L3=(-b±b2-4ac)/2c;
(30)η=L3(℘(ξ)-A-2)2(℘(ξ)-B-2)2,a(℘(ξ)-B-2)2cL3(℘(ξ)-A-2)2
with A-2=(a(12ac-b2)+b(8ac-b2)L3)/48[(2ac-b2)L3-ab], B-2=(ab(5b2-16ac)+(5b4+12a2c2-21ab2c)L3)/48[a(2ac-b2)+b(3ac-b2)L3], and L3=(-b±b2-4ac)/2c; or
(31)η=(A-3℘(ξ)+B-3+E℘′(ξ))2(A-4℘(ξ)+B-4-D℘′(ξ))2,a(A-4℘(ξ)+B-4-D℘′(ξ))2c(A-3℘(ξ)+B-3+E℘′(ξ))2
with A-3=D(2aD2+bE2)L-, B-3=(D/48)(12acE2-b2E2+4abD2)L-, A-4=E(bD2+2cE2)L-, and B-4=(E/48)(12acD2-b2D2+4bcE2)L-, where L-=±(1/2cE4+bE2D2+aD4), and D,E are two arbitrary constants.Note that in any pair of solutions, the latter can be obtained from the former through formula (2). Indeed, combining solutions (24)–(31) with the nonlinear iterated formulae (2)–(20) and their infinitely many derivatives can yield infinitely many Weierstrass elliptic function solutions of (1).


## 2.3. Infinitely Many Jacobi Elliptic Function Solutions to (1)

Assume thatK, C, and D are arbitrary constants, and m∈(0,1) is the elliptic modulus; then (1) has the following Jacobi elliptic solutions.(I)a=4K,  b=-4(1+m2),  c=4m2/K,
(32)η=Kcd(ξ,m)2,Ksn(ξ,m)2,Km2dc(ξ,m)2,Km2ns(ξ,m)2.
Indeed, the former two solutions and the latter two ones can transform each other by formula (2). In what follows, we omit the solutions that can be generated by formula (2).(II)a=4K(1-m2),  b=-4(1-2m2),  c=-4m2/K,
(33)η=Kcn(ξ,m)2,K(1-m2)sd(ξ,m)2.(III)a=-4K(1-m2),  b=4(2-m2),  c=-4/K,
(34)η=Kdn(ξ,m)2,-Kcs(ξ,m)2.(IV)a=K(m2-1),  b=2(1+m2),  c=(m2-1)/K,
(35)η±=K[nd(ξ,m)±msd(ξ,m)]2,-K[nc(ξ,m)±(ξ,m)]2,K1-m2[dn(ξ,m)±mcn(ξ,m)]2,Km2-1[cs(ξ,m)±ds(ξ,m)]2.(V)a=((C2-D2)/K)(m2-1),  b=2(1+m2),  c=-(K/(C2-D2))(m2-1),
(36)η±=[Dcn(ξ,m)+Cdn(ξ,m)]2K[1±1+C2(1-m2)/(D2-C2)sn(ξ,m)]2.(VI)a=K,  b=2(-2+m2),  c=m4/K,
(37)η±=Kcn(ξ,m)2[1-m2±dn(ξ,m)]2,Ksn(ξ,m)2[1±dn(ξ,m)]2,Km2[cd(ξ,m)±i1-m2sd(ξ,m)]2,-Km2[cn(ξ,m)±isn(ξ,m)]2.(VII)a=Km2/(C2-D2),  b=2(-2+m2),  c=(C2-D2)m2/K,
(38)η±=-(K[C+Ddn(ξ,m)]2)×([cn(ξ,m)±D2m2D2-C2-1sn(ξ,m)]2(C2-D2)2×[cn(ξ,m)±D2m2D2-C2-1sn(ξ,m)]2)-1.(VIII)a=K,  b=2(1-2m2),  c=1/K,
(39)η±=K[ns(ξ,m)±cs(ξ,m)]2,K[msn(ξ,m)±idn(ξ,m)]2,K[dc(ξ,m)±1-m2(ξ,m)]2,K[mcd(ξ,m)±i1-m2nd(ξ,m)]2.(IX)a=Km2/(C2-D2),  b=2(1-2m2),  c=(C2-D2)m2/K,
(40)η±=-K[dn(ξ,m)±D2/(D2-C2)-m2sn(ξ,m)]2[C+Dcn(ξ,m)]2.(X)a=K(1-m)2,  b=2(1+6m+m2),  c=(1-m)2/K,
(41)η±=-K[1±msn(ξ,m)]2[1∓msn(ξ,m)]2.(XI)a±=±16Km,  b±=-4(1±6m+m2),  c±=(1±m)2/K,
(42)η±=K[1∓msn(ξ,m)2]2cn(ξ,m)2dn(ξ,m)2,K[ns(ξ,m)±msn(ξ,m)]2(1-m)2.(XII)a±=±16K1-m2,  b±=4(2-m2±61-m2),  c±=4(2-m2±21-m2)/K,
(43)η±=-K[cn(ξ,m)2±1-m2sn(ξ,m)2]2dn(ξ,m)2,-K[dn(ξ,m)±1-m2nd(ξ,m)]22-m2±21-m2,Km4sn(ξ,m)2×[(1-m2∓1-m2)nc(ξ,m)+m2cn(ξ,m)]2,Km4sn(ξ,m)2cn(ξ,m)2[dn(ξ,m)2∓1-m2]2.Combining solutions (32)–(43) with the nonlinear iterated formulae (2)–(20) and their derivatives can give infinitely many Jacobi elliptic function solutions of (1), and these solutions can be degenerated to solutions with hyperbolic or trigonometric functions as long as letting m go to 1 or 0.


## 3. Infinitely Many Elliptic Solutions of Some Nonlinear PDEs

Now we apply auxiliary equation methods to some nonlinear PDEs and obtain their infinitely many elliptic function solutions.Example 1.
First we consider the Kadomtsev-Petviashvili (KP) equation(44)(ut-6uux+uxxx)x+3uyy=0,
which was introduced in 1970 [27] and is of both physical and mathematical interests. By using the wave variable ξ=αx+βy+γt and integrating twice, the KP equation (44) can be converted to the ODE
(45)α4u′′+(αγ+3β2)u-3α2u2=0,
if we let the integration constants be zeros.
Assume that the solutions of (45) have the following form:
(46)u(ξ)=a0+a1η(ξ)+a2η2(ξ)b0+b1η(ξ),
where η(ξ) satisfies (1), ai (i=0,1,2) and bj (j=0,1) are constants to be determinated. Substituting (46) into (45) and equating corresponding coefficients of η(ξ) to zero yield an algebraic system of ai, bj, α, β, and γ. Solving this system with the aid of Maple, we obtain the following solutions of (45):
(47)uI(ξ)=α2(b±b2-3ac6+cη(ξ)2),uII(ξ)=α2(b±b2-3ac6+a2η(ξ)),uIII(ξ)=cα2[a+2bη(ξ)+3cη2(ξ)]2[b∓b2-3ac+3cη(ξ)],
where η(ξ) satisfies (1), with γ=(-3β2±α4b2-3ac)/α, and α, β, a, b, and c are arbitrary constants.
(48)uIV(ξ)=α2(b±b2+12ac6+a2η(ξ)+cη(ξ)2),
where η(ξ) satisfies (1), with γ=(-3β2±α4b2+12ac)/α, α, β, a, b, c are arbitrary constants.
Now substituting the solutionsη’s of (1) into (47) and (48), we obtain immediately infinitely many Weierstrass and Jacobi elliptic solutions of (44), and these elliptic solutions can be degenerated to solutions with hyperbolic or trigonometric functions as long as letting m go to 1 or 0.
For instance, from (48), (6), (9), and (24) (or (35)) we can obtain infinitely many Weierstrass or Jacobi elliptic solutions of (44) as follows:
(49)un(x,y,t)=α2(b±b2+12ac6+a2ηn+cηn2),(n=1,2,…),ηn+1±=±a(b∓2acηn+ηn′)ac(b∓2acηn-ηn′),η1±=±a(b∓2acw1+w1′)ac(b∓2acw1-w1′),w1=-ab±ab2-4ac-4acw0-c(b±b2-4ac)w022c(a+bw0+cw02).(i)
For Weierstrass elliptic function solutions,(50)w0=3a12℘(ξ;g2,g3)-b,
with
(51)g2=b2-3ac12,g3=b(9ac-2b2)432,
where a, b, and c are arbitrary constants.(ii)
For Jacobi elliptic function solutions,(52)w0=-K[nc(ξ,m)±sc(ξ,m)]2,
with
(53)a=K(m2-1),b=2(1+m2),c=(m2-1)K,γ=-3β2±α4(2+2m2)2+12(m2-1)2α,ξ=αx+βy+γt,
where α, β, and K are arbitrary constants and m∈(0,1) is the elliptic modulus. If m tends to 1 or 0, we can also obtain multiple soliton or periodic solutions.Example 2.
Next we investigate infinitely many elliptic function solutions to a negative order KdV equation(54)(uxxu)t+2uux=0,
which actually comes from the negative KdV hierarchy and can be transformed to the Camassa-Holm equation through a gauge transform (see [28]). The negative order KdV equation (54) will be transformed to the ODE
(55)αγu′′+u3+λu=0
upon using wave variable ξ=αx+γt and integrating once, where λ is a integration constant to be determinated. Similar to the procedure of Example 1, one can get solutions of (55) as follows,
(56)u(ξ)=±αγ[-b+2(-1)j-1ac][ac-(-1)j-1cη(ξ)]2[ac+(-1)j-1cη(ξ)],(j=1,2),
where η(ξ) satisfies (1), with λ=αγb/2, α, γ, a, b, c are arbitrary constants.
Similarly, substituting the solutionsη’s of (1) into (56), we obtain infinitely many Weierstrass and Jacobi elliptic solutions of (54), among which including hyperbolic and trigonometric functions solutions as long as letting m go to 1 or 0.
For instance, from (56), (3), (9) and (27) (or (34)) we can obtain infinitely many Weierstrass or Jacobi elliptic solutions of (54) as follows
(57)un(x,t)=±αγ[-b+2(-1)j-1ac][ac-(-1)j-1cηn]2[ac+(-1)j-1cηn],ηn+1=±a(b∓2acηn+ηn′)ac(b∓2acηn-ηn′),(n=0,1,2,…,j=1,2),η1±=±a(b∓2acw1+w1′)ac(b∓2acw1-w1′),w1±=±2a+(b∓b2-4ac)w0∓b+b2-4ac∓2cw0,(i)
for Weierstrass elliptic function solutions,(58)w0=576℘′2(ξ;g2,g3)c(24℘(ξ;g2,g3)+b)2,
with
(59)g2=b2192+ac16,g3=b(36ac-b2)13824,
where a, b and c are arbitrary constants;(ii)
for Jacobi elliptic function solutions,(60)w0=-Kcs(ξ,m)2,
with
(61)a=4K(-1+m2),b=4(2-m2),c=-4K,λ=αγ(8-4m2)2,
where α, γ, K are arbitrary constants and m∈(0,1) is the elliptic modulus. If m tends to 1 or 0, we can also obtain multiple soliton or periodic solutions.Example 3.
We consider the two-dimensional Davey-Stewartson equation(62)iut+uxx-uyy-2|u|2u-2uv=0,vxx+vyy+2(|u|2)xx=0.
This equation is completely integrable and often used to describe the long time evolution of a two-dimensional wave packet [29, 30]. In recent years, various methods such as inverse scattering, Darboux transformation, variable separation, and Bäcklund transformation have been used to solve the equation, respectively. In this work, we try to deal with (62) by new idea and give it infinitely many Weierstrass and Jacobi elliptic-like solutions.
For our purpose, we introduce the following transformations:(63)u(x,t)=U(ξ)eiθ,v(x,t)=V(ξ),θ=px+qy+rt,ξ=kx+dy-2(pk-dq)t.
Then system (62) is reduced to the following system of ordinary differential equations:
(64)(q2-p2-r)U+(k2-d2)U′′-2U3-2UV=0,(k2+d2)V′′+2k2(U2)′′=0.
Integrating the second equation in (64) and neglecting the constants of integration we find that
(65)V=-2k2k2+d2U2.
Substituting (65) into the first equation of the system leads to
(66)U′′+q2-p2-rk2-d2U+2k2+d2U3=0,
where prime denotes differentiation with respect to ξ. Further multiplying (66) by U′(ξ) and integrating it, we obtain that
(67)U′2+q2-p2-rk2-d2U2+1k2+d2U4-H=0,
which is regrouped as follows:
(68)U′2=-1k2+d2U4+p2+r-q2k2-d2U2+H,
where H is an arbitrary constant of integration.
In order to write (68) in a simple form that we want, let us introduce the following transformation:
(69)η(ξ)=U2(ξ).
Then (68) can be written as (1) with
(70)a=4H,b=4(p2+r-q2)k2-d2,c=-4k2+d2.
Now from the solutionsη’s of (1) and taking into account (63), (65), and (69), we obtain infinitely many Weierstrass and Jacobi elliptic-like solutions to the two-dimensional Davey-Stewartson equation (62) as follows:
(71)un(x,t)=Un(ξ)eiθ,vn(x,t)=-2k2k2+d2Un2(ξ),(n=1,2,…),Un2(ξ)=ηn(ξ),
where θ=px+qy+rt, ξ=kx+dy-2(pk-dq)t; taking account of (70), p, q, r, k, d, and H satisfy the relations such that (1) has an elliptic solution.
For instance, from (6), (9), and (24) (or (35)) we can obtain infinitely many Weierstrass or Jacobi elliptic-like solutions to the two-dimensional Davey-Stewartson equation (62) as follows:
(72)un(x,t)=Un(ξ)eiθ,vn(x,t)=-2k2k2+d2Un2(ξ),(n=0,1,2,…),Un2(ξ)=ηn(ξ),ηn+1=-a(b+2acηn+ηn′)ac(b+2acηn-ηn′),η0=(-ab±ab2-4ac-4acw0-c(b±b2-4ac)w02)×(2c(a+bw0+cw02))-1.(i)
For Weierstrass elliptic function solutions,(73)w0=-k2+d23(6℘(ξ;g2,g3)-p2+r-q2k2-d2),
with
(74)θ=px+qy+rt,ξ=kx+dy-2(pk-dq)t,g2=4Hk2+d2+4(p2+r-q2)23(k2-d2)2,g3=4(q2-p2-r)27(k2-d2)((p2+r-q2)2(k2-d2)2+9Hk2+d2),a=4H,b=4(p2+r-q2)k2-d2,c=-4k2+d2,
where p, q, r, k, d, and H are arbitrary real constants.(ii)
For Jacobi elliptic function solutions,(75)w0=14(k2+d2)(1-m2)(nd(ξ,m)±msd(ξ,m))2,
with
(76)θ=px+qy+rt,ξ=kx+dy-2(pk-dq)t,a=-14(1-m2)2(k2+d2),b=2(1+m2),c=-4k2+d2,r=12(1+m2)(k2-d2)+q2-p2,
where d, k, p, and q are arbitrary constants and m∈(0,1) is the elliptic modulus. If m tends to 1 or 0, we can also obtain multiple soliton or periodic solutions.


## 4. Discussion

In this paper, we have obtained some nonlinear iterated formulae of solutions to the simple equation (1) and its infinitely many Weierstrass and Jacobi elliptic function solutions. Based on this fact, we can obtain infinitely many Weierstrass and Jacobi elliptic function solutions of some nonlinear equations. Our results suggest that the proposed method is reliable and effective to find infinitely many elliptic solutions to some nonlinear problems. The method is concise, and it can also be applied to other kinds of nonlinear equations in mathematical physics.



---

*Source: 582532-2013-07-31.xml*