# Exp-Function Method for a Generalized MKdV Equation

Authors: Yuzhen Chai; Tingting Jia; Huiqin Hao; Jianwen Zhang
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

Under investigation in this paper is a generalized MKdV equation, which describes the propagation of shallow water in fluid mechanics. In this paper, we have derived the exact solutions for the generalized MKdV equation including the bright soliton, dark soliton, two-peak bright soliton, two-peak dark soliton, shock soliton and periodic wave solution via Exp-function method. By figures and symbolic computations, we have discussed the propagation characteristics of those solitons under different values of those coefficients in the generalized MKdV equation. The method constructing soliton solutions in this paper may be useful for the investigations on the other nonlinear mathematical physics model and the conclusions of this paper can give theory support for the study of dynamic features of models in the shallow water.


---

## Body

## 1. Introduction

Nonlinear science studies nonlinear phenomena of the world, which is cross-disciplinary [1–4]. In fact, any science, whether natural science or social science, has its own nonlinear phenomena and problems. To study these phenomena and problems, many branches of nonlinear science have been promoted to be built and developed [5–7]. Obviously, most of the phenomena are not linear in the natural sciences and engineering practice; thus many problems cannot be researched and solved by the linear methods, which makes the study of nonlinear science very significant. Actually, nonlinear science, which mainly consists of soliton, chaos, and fractal [3, 4, 7], is not the simple superposition and comprehensive of these nonlinear branches but a comprehensive subject to study the various communist rules in nonlinear phenomena. With the rapid development of nonlinear science, the study of exact solutions of nonlinear evolution equations has attracted much attention of many mathematicians and physicists.In 1834, the solitary wave phenomenon was observed by Huang et al. [8] and Russell [9]. Later, people named the isolation water peak, which kept moving with constant shape and speed on the surface of the water, as a solitary wave [8]. The discovery of a solitary wave turned people into a new field in the study of the waves of the convection. In 1895, Holland’s Professor Kortewrg and his disciples Vries derived the famous KdV equation from the research of shallow water wave motion [10]. The analysis to the KdV equation has improved to recover inverse scattering method, on which people expanded the new research directions of algebra and geometry [8]:
(1)
u
t
±
6
u
u
x
+
u
x
x
x
=
0
.
MKdV equation plays an important role in describing the plasmas and phonon in anharmonic lattice. In the nineteen fifties, physicist Fermi, Pasta, and Ulam made a famous FPU experiment, connecting the 64 particles by nonlinear spring, thereby forming a nonlinear vibrating string. Although the FPU experiment did not gain the solutions of solitary wave, it will expand the study of the solitary wave to the field outside the mechanics [11]. Later, Fermi et al. studied the nonlinear vibration problem of FPU model and obtained the solitary wave solutions [11], which is the right answer for the question of FPU. It is the first time to find solitary wave solutions in the field outside mechanics after the solution was found in the KdV equation, which give rise to the scientists’ interest in researching the solitary wave phenomenon. In 1965, Zabusky and Kruskal studied solitary waves in plasma [12]; they found that the waveform does not change nature before and after collision of nonlinear solitary waves, which is similar to particle collisions, so Zabusky and Kruskal named the solitary wave with the impact properties of collisions as soliton [12]. Soliton concept is an important milepost on the history of the development of the soliton theory. The next few decades, the soliton theory had a rapid development and penetrated into many areas, such as fluid mechanics, nonlinear optical fiber communication, plasma physics, fluid physics, chemistry, life science, and marine science [13–15].Recently many new approaches to nonlinear wave equations have been proposed, such as Tanh-function method [16–18], F-expansion method [19], Jacobian elliptic function method [20], Darboux transformation method [21–26], adomian method [27–29], variational approach [30], and homotopy perturbation method [31]. All methods mentioned above have their limitations in their applications. We will apply Exp-function method to a generalized MKdV equation to gain exact solutions:
(2)
u
t
+
α
u
2
u
x
+
β
u
x
x
x
=
0
,
where α and β are real parameters. When u
(
x
,
t
)
=
c (c is a constant) as t
→
±
∞, solitons can be obtained.This paper will be organized as follows. In Section2, the basic idea of Exp-function method is introduced. In Section 3, carrying on calculating and illustrating by the mathematical software MATHEMATIC, we will solve the solitary wave solutions of (2) based on the Exp-function method. In Section 4, we will obtain shock soliton, bright soliton, two-peak bright soliton, dark soliton, two-peak dark soliton, and periodic wave solutions and analyze the dynamic features of soliton solutions by using some figures. Finally, our conclusions will be addressed in Section 5.


## 2. Basic Idea of Exp-Function Method

In order to illustrate the basic idea of the suggested method, we consider firstly the following general partial differential equation:(3)
F
(
u
,
u
t
,
u
x
,
u
x
t
,
u
x
x
,
…
)
=
0
.
We aim at its exact solutions, so we introduce a complex variable, ξ, defined as
(4)
ξ
=
k
x
+
ω
t
,
where ω and k are constants unknown to be further determined. Therefore we can convert (3) into an ordinary differential equation with respect to ξ:
(5)
F
(
u
,
u
ξ
,
u
ξ
ξ
,
…
)
=
0
.
Very simple and straightforward, the Exp-function method is based on the assumption that traveling wave solutions can be expressed in the following form:
(6)
u
(
ξ
)
=
∑
n
=
-
d
c
a
n
exp
⁡
(
n
ξ
)
∑
m
=
-
q
p
b
m
exp
⁡
(
m
ξ
)
,
where c, d, p, and q are positive integers which are unknown to be further determined and a
n and b
m are unknown constants.We suppose thatu
(
ξ
) which is the solution of (5) can be expressed as
(7)
u
(
ξ
)
=
a
c
exp
⁡
(
c
ξ
)
+
⋯
+
a
-
d
exp
⁡
(
-
d
ξ
)
b
p
exp
⁡
(
p
ξ
)
+
⋯
+
b
-
q
exp
⁡
(
-
q
ξ
)
.To determine the values ofc and p, we balance the linear term of highest order in (5) with the highest order nonlinear term. Similarly we balance the lowest orders to confirm d and q. For simplicity, we set some particular values for c, p, d, and q and then change the left side of (5) into the polynomial of exp
⁡
(
n
ξ
). Equating the coefficients of exp
⁡
(
n
ξ
) to be zero results in a set of algebraic equations. Then solving the algebraic system with symbolic computation system, we can gain the solution. Substituting it into (3), we have the general form of the exact solution expressed as the form of exp
⁡
(
n
ξ
). Taking the form of the solutions into consideration, we can have more extensive solutions via Exp-function method.


## 3. Solving Generalized MKdV Equation via Exp-Function Method

Now we consider the follwing equation [32, 33]:
(8)
u
t
+
u
2
u
x
+
u
x
x
x
=
0
.
This equation is called modified KdV equation, which arises in the process of understanding the role of nonlinear dispersion and in the formation of structures like liquid drops. The KdV equation is one of the most familiar models for solitons and the foundation which studies other equations. Developing upon this foundation, the isolated theories are treated as the milestone of mathematics physical method. During researching the isolated theories, the remarkable application should be laser shooting practice and fiber-optic communication.For researching the variety of the solutions for the modified KdV equation, we carry on the MKdV equation to expand. Thus, we have(9)
u
t
+
α
u
2
u
x
+
β
u
x
x
x
=
0
,
where α and β are real parameters. In the following, we consider (9).Introduce a complex variable,η, defined as
(10)
u
=
u
(
η
)
,
η
=
k
x
+
ω
t
.
By (10), (9) becomes
(11)
ω
u
′
+
α
k
u
2
u
′
+
β
k
3
u
′′′
=
0
,
which denotes the differential with respect to η.The Exp-function method is based on the assumption that traveling wave solutions can be expressed in the following form:(12)
u
(
η
)
=
∑
n
=
-
c
d
a
n
exp
⁡
(
n
η
)
∑
m
=
-
p
q
b
m
exp
⁡
(
m
η
)
,
where c, d, p, and q are positive integers unknown to be further determined and a
n and b
m are unknown constants.We suppose that the solution of (11) can be expressed as
(13)
u
(
η
)
=
a
c
exp
⁡
(
c
η
)
+
⋯
+
a
-
d
exp
⁡
(
-
d
η
)
b
p
exp
⁡
(
p
η
)
+
⋯
+
b
-
q
exp
⁡
(
-
q
η
)
.
To determine values of c and p, we balance the linear term of highest order in (11) with the highest order nonlinear term. By simple calculation, we have
(14)
u
′′′
=
c
1
exp
⁡
[
(
7
p
+
c
)
η
]
+
⋯
c
2
exp
⁡
[
8
p
η
]
+
⋯
,
u
2
u
′
=
c
3
exp
⁡
[
(
p
+
3
c
)
η
]
+
⋯
c
4
exp
⁡
[
4
p
η
]
+
⋯
=
c
3
exp
⁡
[
(
5
p
+
3
c
)
η
]
+
⋯
c
4
exp
⁡
[
8
p
η
]
+
⋯
,
where c
i
(
i
=
1,2
,
3,4
) are coefficients only for simplicity.Balancing highest order of exponential function in (14), we have
(15)
5
p
+
c
=
7
p
+
c
,
which leads to the result
(16)
p
=
c
.
Similarly to make sure of the values of d and q, we balance the linear term of lowest order in (11):
(17)
u
′′′
=
⋯
+
d
1
exp
⁡
[
-
(
7
q
+
d
)
η
]
⋯
+
d
2
exp
⁡
[
-
8
q
η
]
,
u
2
u
′
=
⋯
+
d
3
exp
⁡
[
-
(
q
+
3
d
)
η
]
⋯
+
d
4
exp
⁡
[
-
8
q
η
]
=
⋯
+
d
3
exp
⁡
[
-
(
5
q
+
3
d
)
η
]
⋯
+
d
4
exp
⁡
[
-
8
q
η
]
,
where d
i
(
i
=
1,2
,
3,4
) are determined to be coefficients only for simplicity.Balancing lowest order of exponential function in (17), we have
(18)
-
(
7
q
+
d
)
=
-
(
5
q
+
3
d
)
,
which leads to the result
(19)
q
=
d
.
For simplicity, we set p
=
c
=
1 and q
=
d
=
1, so (13) reduces to
(20)
u
(
η
)
=
a
1
exp
⁡
(
η
)
+
a
0
+
a
-
1
exp
⁡
(
-
η
)
exp
⁡
(
η
)
+
b
0
+
b
-
1
exp
⁡
(
-
η
)
.
Substituting (20) into (11), and by the help of Mathematica, we have
(21)
1
A
[
C
3
exp
⁡
(
3
η
)
+
C
2
exp
⁡
(
2
η
)
hh
+
C
1
exp
⁡
(
η
)
+
C
0
+
C
-
1
exp
⁡
(
-
η
)
hh
+
C
-
2
exp
⁡
(
-
2
η
)
+
C
-
3
exp
⁡
(
-
3
η
)
]
=
0
,
where
(22)
A
=
(
b
-
1
+
e
2
η
+
b
0
e
η
)
4
,
C
3
=
(
h
+
k
α
a
1
2
)
(
-
a
0
+
a
1
b
0
)
,
C
2
=
-
2
k
α
a
0
2
a
1
-
2
a
-
1
(
3
k
3
β
+
h
+
k
α
a
1
2
)
+
2
a
0
(
3
k
3
β
-
h
+
k
α
a
1
2
)
b
0
+
2
a
1
[
(
3
k
3
β
+
h
+
k
α
a
1
2
)
b
-
1
+
(
-
3
k
3
β
+
h
)
b
0
2
]
,
C
1
=
a
0
[
-
k
α
(
a
0
2
+
6
a
-
1
a
1
)
hhhh
+
(
24
k
3
β
-
h
+
5
k
α
a
1
2
)
b
-
1
]
-
h
b
0
2
(
a
0
-
a
1
b
0
)
+
{
a
-
1
(
-
5
h
+
k
α
a
1
2
)
h
h
+
a
1
[
k
α
a
0
2
+
6
(
-
4
k
3
β
+
h
)
b
-
1
]
}
b
0
,
C
0
=
-
4
(
a
-
1
-
a
1
b
-
1
)
×
[
k
α
(
a
0
2
+
a
-
1
a
1
)
-
(
9
k
3
β
-
h
)
b
-
1
+
h
b
0
2
]
,
C
-
1
=
a
-
1
{
a
0
2
6
k
α
a
0
a
1
b
-
1
hhhhh
+
[
-
k
α
a
0
2
+
6
(
4
k
3
β
-
h
)
b
-
1
-
h
b
0
2
]
b
0
}
-
k
α
a
-
1
2
(
5
a
0
+
a
1
b
0
)
+
b
-
1
{
k
α
a
0
3
+
5
h
a
1
b
-
1
b
0
h
h
h
h
h
h
h
+
a
0
[
(
-
24
k
3
β
+
h
)
b
-
1
+
h
b
0
2
]
}
,
C
-
2
=
2
b
-
1
2
(
a
-
1
-
a
1
b
-
1
)
(
3
k
3
β
+
h
)
+
2
b
0
b
-
1
(
a
-
1
b
0
-
a
0
b
-
1
)
(
-
3
k
3
β
+
h
)
-
2
k
α
a
-
1
×
[
a
-
1
(
a
-
1
+
a
0
b
0
)
-
b
-
1
(
a
0
2
+
a
1
a
-
1
)
]
,
C
-
3
=
-
(
k
α
a
-
1
2
+
h
b
-
1
2
)
(
-
a
0
b
-
1
+
a
-
1
b
0
)
,
with h
=
k
3
β
+
ω.Equating the coefficients ofexp
⁡
(
n
η
) to be zero, we have
(23)
C
3
=
0
,
C
2
=
0
,
C
1
=
0
,
C
0
=
0
,
C
-
3
=
0
,
C
-
2
=
0
,
C
-
1
=
0
.
Solving (23), simultaneously, we obtain
(24)
ω
=
-
k
3
β
-
k
α
a
1
2
,
a
-
1
=
3
k
2
β
b
0
2
+
2
α
a
1
2
b
0
2
8
α
a
1
,
a
0
=
3
k
2
β
b
0
+
α
a
1
2
b
0
α
a
1
,
b
-
1
=
3
k
2
β
b
0
2
+
2
α
a
1
2
b
0
2
8
α
a
1
2
.
Substituting (24) into (20) results in a compact-like solution, which reads
(25)
u
(
x
,
t
)
=
a
1
(
1
+
24
y
k
2
β
b
0
3
k
2
β
b
0
2
+
2
α
a
1
2
(
2
y
+
b
0
)
2
)
,
where a
1, b
0, α, β, and k are real parameters and y
=
e
k
x
+
t
(
-
k
3
β
-
k
α
a
1
2
). Besides, the obtained solution equation (25) is a generalized soliton solution of (9).In case,k is an imaginary parameter, the obtained soliton solution can be converted into the periodic solution or compact-like solution. We write
(26)
k
=
i
K
,
where K is a real parameter.Use the transformations(27)
exp
⁡
(
k
x
+
k
2
+
6
a
1
+
1
k
t
)
=
cos
⁡
(
K
x
+
-
K
2
+
6
a
1
+
1
K
t
)
+
i
sin
(
K
x
+
-
K
2
+
6
a
1
+
1
K
t
)
,
exp
⁡
(
-
k
x
-
k
2
+
6
a
1
+
1
k
t
)
=
cos
⁡
(
K
x
+
-
K
2
+
6
a
1
+
1
K
t
)
-
i
sin
(
K
x
+
-
K
2
+
6
a
1
+
1
K
t
)
.
Equation (25) becomes
(28)
u
(
x
,
t
)
=
a
1
(
1
+
(
×
(
-
3
K
2
β
b
0
2
+
2
α
a
1
2
(
2
(
cos
⁡
(
y
1
)
+
i
sin
(
y
1
)
)
+
b
0
)
2
)
-
1
)
-
24
K
2
β
(
cos
⁡
(
y
1
)
+
i
sin
(
y
1
)
)
b
0
h
h
h
h
h
h
h
h
h
×
(
(
2
(
cos
⁡
(
y
1
)
+
i
sin
(
y
1
)
)
+
b
0
)
2
-
3
K
2
β
b
0
2
+
2
α
a
1
2
h
h
h
h
h
h
h
h
h
h
h
h
×
(
2
(
cos
⁡
(
y
1
)
+
i
sin
(
y
1
)
)
+
b
0
)
2
)
-
1
)
)
,
where a
1, b
0, α, β, and K are real parameters and y
1
=
K
x
-
(
α
K
a
1
2
-
β
K
3
)
t. Simplifying the above formula, we have
(29)
u
(
x
,
t
)
=
a
1
+
-
3
K
2
b
0
β
/
α
a
1
b
0
+
(
1
+
ρ
)
cos
⁡
(
y
1
)
+
i
(
1
-
ρ
)
sin
(
y
1
)
,
where a
1, b
0, α, β, and K are real parameters and ρ
=
b
0
2
(
-
3
β
K
2
+
2
α
a
1
2
)
/
8
α
a
1
2.If we search for a periodic solution or compact-like solution, the imaginary part in the denominator of (29) must be zero, which requires that
(30)
1
-
ρ
=
1
-
b
0
2
(
-
3
β
K
2
+
2
α
a
1
2
)
8
α
a
1
2
=
0
.
Solving b
0 from (30), we obtain
(31)
b
0
=
±
8
α
a
1
2
-
3
β
K
2
+
2
α
a
1
2
.
Substituting (31) into (29) results in a compact-like solution, which reads
(32)
u
(
x
,
t
)
=
a
1
-
6
2
K
2
β
g
2
2
α
a
1
g
±
cos
⁡
(
y
1
)
(
-
α
a
1
+
g
2
(
3
K
2
β
-
2
α
a
1
2
)
)
,
where a
1 and K are free parameters, and y
1
=
K
x
-
(
α
K
a
1
2
-
β
K
3
)
t, g
=
α
a
1
2
/
(
-
3
K
2
β
+
2
α
a
1
2
), and it requires that α
a
1
2
/
(
-
3
K
2
β
+
2
α
a
1
2
)
>
0.


## 4. Discussing the Forms of the Solutions

In the above, we gain two cases of solution equations (25) and (32). In order to get the richness of solutions of (9), we have these two cases discussed in the following.


### 4.1. The Solutions from (25)

Through discussions, we obtain the following results.


#### 4.1.1. Whenα
=
0 and β
=
1

whenα
=
0 and β
=
1, (9) becomes u
t
+
u
x
x
x
=
0. Then we can obtain shock solitons. The solution formula is shown as follows:
(33)
u
(
x
,
t
)
=
a
1
(
1
+
8
e
-
k
3
t
+
k
x
b
0
)
=
a
1
+
8
a
1
cosh
⁡
(
k
3
t
-
k
x
)
b
0
-
8
a
1
sinh
⁡
(
k
3
t
-
k
x
)
b
0
,
where a
1, b
0, and k are real parameters.Case 1 (whena
1
=
1, b
0
=
2, k
=
3).
Consider(34)
u
1
(
x
,
t
)
=
1
+
4
e
-
27
t
+
3
x
=
1
+
4
cosh
⁡
(
27
t
-
3
x
)
-
4
sinh
⁡
(
27
t
-
3
x
)
.
Description is as shown in Figure 1(a).Evolution of the solutions of (25). (a) Shock soliton: parameters are α
=
0, β
=
1, a
1
=
1, b
0
=
2, and k
=
3. (b) A shock soliton in collision with a dark soliton: parameters are α
=
0, β
=
1, a
1
=
3, b
0
=
-
2, and k
=
4.
(a)
(b)Case 2 (whena
1
=
3, b
0
=
-
2, and k
=
4).
Consider(35)
u
2
(
x
,
t
)
=
3
(
1
-
4
e
-
64
t
+
4
x
)
=
3
-
12
cosh
⁡
(
64
t
-
4
x
)
+
12
sinh
⁡
(
64
t
-
4
x
)
.
Description is as shown in Figure 1(b).


#### 4.1.2. Whenα
=
1 and β
=
1

whenα
=
1 and β
=
1, (9) becomes u
t
+
u
2
u
x
+
u
x
x
x
=
0. Then we can obtain one-peak bright solitons. The solution formula is shown as follows:
(36)
u
(
x
,
t
)
=
a
1
(
1
+
24
e
y
k
2
b
0
3
k
2
b
0
2
+
2
a
1
2
(
2
e
y
+
b
0
)
2
)
=
a
1
(
1
+
24
k
2
(
cosh
⁡
(
y
)
+
sinh
⁡
(
y
)
)
b
0
3
k
2
b
0
2
+
2
a
1
2
(
2
cosh
⁡
(
y
)
+
2
sinh
⁡
(
y
)
+
b
0
)
2
)
,
where a
1, b
0, and k are real parameters and y
=
k
x
+
t
(
-
k
3
-
k
a
1
2
).Case 3 (whena
1
=
1, b
0
=
2, and k
=
2).
Consider(37)
u
3
(
x
,
t
)
=
1
+
192
e
-
10
t
+
2
x
48
+
2
(
2
+
2
e
-
10
t
+
2
x
)
2
=
1
+
192
(
cosh
⁡
(
10
t
-
2
x
)
-
sinh
⁡
(
10
t
-
2
x
)
)
48
+
2
(
2
+
2
cosh
⁡
(
10
t
-
2
x
)
-
2
sinh
⁡
(
10
t
-
2
x
)
)
2
.
Description is as shown in Figure 2(a).Evolution of the solutions of (25). (a) Bright soliton: parameters are α
=
1, β
=
1, a
1
=
1, b
0
=
2, and k
=
2. (b) A bright soliton in collision with two dark solitons: parameters are α
=
1, β
=
1, a
1
=
1, b
0
=
-
1, and k
=
-
3.
(a)
(b)Case 4 (whena
1
=
1, b
0
=
-
1, and k
=
-
3).
Consider(38)
u
4
(
x
,
t
)
=
1
-
216
e
30
t
-
3
x
27
+
2
(
-
1
+
2
e
30
t
-
3
x
)
2
=
1
-
216
(
cosh
⁡
(
30
t
-
3
x
)
+
sinh
⁡
(
30
t
-
3
x
)
)
27
+
2
(
-
1
+
2
cosh
⁡
(
30
t
-
3
x
)
+
2
sinh
⁡
(
30
t
-
3
x
)
)
2
.
Description is as shown in Figure 2(b).


#### 4.1.3. Whenα
=
1 and β
=
-
1

whenα
=
1 and β
=
-
1, (9) becomes u
t
+
u
2
u
x
-
u
x
x
x
=
0. Then we can obtain bright solitons and dark solitons. The solution formula is shown as follows:
(39)
u
(
x
,
t
)
=
a
1
(
1
-
24
e
y
k
2
b
0
-
3
k
2
b
0
2
+
2
a
1
2
(
2
e
y
+
b
0
)
2
)
=
a
1
(
1
-
24
k
2
(
cosh
⁡
(
y
)
+
sinh
⁡
(
y
)
)
b
0
-
3
k
2
b
0
2
+
2
a
1
2
(
2
cosh
⁡
(
y
)
+
2
sinh
⁡
(
y
)
+
b
0
)
2
)
,
where a
1, b
0, and k are real parameters and y
=
k
x
+
t
(
k
3
-
k
a
1
2
).Case 5 (whena
1
=
-
5, b
0
=
-
2, and k
=
3).
Consider(40)
u
5
(
x
,
t
)
=
-
5
(
1
+
432
e
-
48
t
+
3
x
-
108
+
50
(
-
2
+
2
e
-
48
t
+
3
x
)
2
)
=
-
5
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
1
+
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
432
(
cosh
⁡
(
48
t
-
3
x
)
+
sinh
⁡
(
48
t
-
3
x
)
)
h
h
h
h
h
h
h
h
h
h
×
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
-
108
+
50
(
(
48
t
-
3
x
)
-
2
+
2
cosh
⁡
(
48
t
-
3
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 3(a).Evolution of the solutions of (25). (a) Two-peak bright soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
-
2, and k
=
3. (b) A bright soliton in collision with a dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
3, b
0
=
-
1, and k
=
4.
(a)
(b)Case 6 (whena
1
=
3, b
0
=
-
1, and k
=
4).
Consider(41)
u
7
(
x
,
t
)
=
3
(
1
+
384
e
28
t
+
4
x
-
48
+
18
(
-
1
+
2
e
28
t
+
4
x
)
2
)
=
3
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
)
1
+
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
384
(
cosh
⁡
(
28
t
+
4
x
)
+
sinh
⁡
(
28
t
+
4
x
)
)
h
h
h
h
h
h
h
h
×
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
-
48
+
18
(
(
28
t
+
4
x
)
-
1
+
2
cosh
⁡
(
28
t
+
4
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 3(b).Case 7 (whena
1
=
-
5, b
0
=
2, and k
=
-
3).
Consider(42)
u
6
(
x
,
t
)
=
-
5
(
1
-
432
e
48
t
-
3
x
-
108
+
50
(
2
+
2
e
48
t
-
3
x
)
2
)
=
-
5
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
1
-
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
432
(
cosh
⁡
(
48
t
-
3
x
)
+
sinh
⁡
(
48
t
-
3
x
)
)
h
h
h
h
h
h
h
h
h
×
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
-
108
-
50
(
(
48
t
-
3
x
)
2
+
2
cosh
⁡
(
48
t
-
3
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 4(a).Evolution of the solutions of (25). (a) Dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
2, and k
=
-
3. (b) Two-peak dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
3, and k
=
4.
(a)
(b)Case 8 (whena
1
=
-
5, b
0
=
3, and k
=
4).
Consider(43)
u
8
(
x
,
t
)
=
-
5
(
1
-
1152
e
-
36
t
+
4
x
-
432
+
50
(
3
+
2
e
-
36
t
+
4
x
)
2
)
=
-
5
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
)
1
-
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
1152
(
cosh
⁡
(
36
t
-
4
x
)
-
sinh
⁡
(
36
t
-
4
x
)
)
h
h
h
h
h
h
h
h
h
×
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
-
432
+
50
(
(
36
t
-
4
x
)
3
+
2
cosh
⁡
(
36
t
-
4
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 4(b).


### 4.2. The Solutions from (32)

Based on the form of (32), if we set α
=
0 and β
=
0 at the same time, the solution should be zero. Therefore, it is not significative to study and analyse. So, in the following, we consider α and β under the circumstance that α
=
0 and β
=
0 are not zero in the meantime.Whenα
=
1 and β
=
1, (9) becomes u
t
+
u
2
u
x
+
u
x
x
x
=
0. Then we can obtain periodic wave solutions. The solution formula is shown as follows:
(44)
u
(
x
,
t
)
=
a
1
-
6
2
K
2
g
2
2
a
1
g
±
cos
⁡
(
y
)
(
-
a
1
+
g
2
(
3
K
2
-
2
a
1
2
)
)
,
where y
=
K
x
-
(
K
a
1
2
-
K
3
)
t and g
=
a
1
2
/
(
-
3
K
2
+
2
a
1
2
).To illustrate the analysis, we obtain the following results.Case 9 (whena
1
=
5, K
=
-
2).
Consider(45)
u
9
(
x
,
t
)
=
5
-
24
19
(
10
/
19
+
6
cos
⁡
(
2
(
21
t
-
x
)
)
)
.
Description is as shown in Figure 5.Figure 5
Evolution of the solutions of (32). Periodic wave solution: parameters are α
=
1, β
=
1, a
1
=
5, and K
=
-
2.


## 4.1. The Solutions from (25)

Through discussions, we obtain the following results.


### 4.1.1. Whenα
=
0 and β
=
1

whenα
=
0 and β
=
1, (9) becomes u
t
+
u
x
x
x
=
0. Then we can obtain shock solitons. The solution formula is shown as follows:
(33)
u
(
x
,
t
)
=
a
1
(
1
+
8
e
-
k
3
t
+
k
x
b
0
)
=
a
1
+
8
a
1
cosh
⁡
(
k
3
t
-
k
x
)
b
0
-
8
a
1
sinh
⁡
(
k
3
t
-
k
x
)
b
0
,
where a
1, b
0, and k are real parameters.Case 1 (whena
1
=
1, b
0
=
2, k
=
3).
Consider(34)
u
1
(
x
,
t
)
=
1
+
4
e
-
27
t
+
3
x
=
1
+
4
cosh
⁡
(
27
t
-
3
x
)
-
4
sinh
⁡
(
27
t
-
3
x
)
.
Description is as shown in Figure 1(a).Evolution of the solutions of (25). (a) Shock soliton: parameters are α
=
0, β
=
1, a
1
=
1, b
0
=
2, and k
=
3. (b) A shock soliton in collision with a dark soliton: parameters are α
=
0, β
=
1, a
1
=
3, b
0
=
-
2, and k
=
4.
(a)
(b)Case 2 (whena
1
=
3, b
0
=
-
2, and k
=
4).
Consider(35)
u
2
(
x
,
t
)
=
3
(
1
-
4
e
-
64
t
+
4
x
)
=
3
-
12
cosh
⁡
(
64
t
-
4
x
)
+
12
sinh
⁡
(
64
t
-
4
x
)
.
Description is as shown in Figure 1(b).


### 4.1.2. Whenα
=
1 and β
=
1

whenα
=
1 and β
=
1, (9) becomes u
t
+
u
2
u
x
+
u
x
x
x
=
0. Then we can obtain one-peak bright solitons. The solution formula is shown as follows:
(36)
u
(
x
,
t
)
=
a
1
(
1
+
24
e
y
k
2
b
0
3
k
2
b
0
2
+
2
a
1
2
(
2
e
y
+
b
0
)
2
)
=
a
1
(
1
+
24
k
2
(
cosh
⁡
(
y
)
+
sinh
⁡
(
y
)
)
b
0
3
k
2
b
0
2
+
2
a
1
2
(
2
cosh
⁡
(
y
)
+
2
sinh
⁡
(
y
)
+
b
0
)
2
)
,
where a
1, b
0, and k are real parameters and y
=
k
x
+
t
(
-
k
3
-
k
a
1
2
).Case 3 (whena
1
=
1, b
0
=
2, and k
=
2).
Consider(37)
u
3
(
x
,
t
)
=
1
+
192
e
-
10
t
+
2
x
48
+
2
(
2
+
2
e
-
10
t
+
2
x
)
2
=
1
+
192
(
cosh
⁡
(
10
t
-
2
x
)
-
sinh
⁡
(
10
t
-
2
x
)
)
48
+
2
(
2
+
2
cosh
⁡
(
10
t
-
2
x
)
-
2
sinh
⁡
(
10
t
-
2
x
)
)
2
.
Description is as shown in Figure 2(a).Evolution of the solutions of (25). (a) Bright soliton: parameters are α
=
1, β
=
1, a
1
=
1, b
0
=
2, and k
=
2. (b) A bright soliton in collision with two dark solitons: parameters are α
=
1, β
=
1, a
1
=
1, b
0
=
-
1, and k
=
-
3.
(a)
(b)Case 4 (whena
1
=
1, b
0
=
-
1, and k
=
-
3).
Consider(38)
u
4
(
x
,
t
)
=
1
-
216
e
30
t
-
3
x
27
+
2
(
-
1
+
2
e
30
t
-
3
x
)
2
=
1
-
216
(
cosh
⁡
(
30
t
-
3
x
)
+
sinh
⁡
(
30
t
-
3
x
)
)
27
+
2
(
-
1
+
2
cosh
⁡
(
30
t
-
3
x
)
+
2
sinh
⁡
(
30
t
-
3
x
)
)
2
.
Description is as shown in Figure 2(b).


### 4.1.3. Whenα
=
1 and β
=
-
1

whenα
=
1 and β
=
-
1, (9) becomes u
t
+
u
2
u
x
-
u
x
x
x
=
0. Then we can obtain bright solitons and dark solitons. The solution formula is shown as follows:
(39)
u
(
x
,
t
)
=
a
1
(
1
-
24
e
y
k
2
b
0
-
3
k
2
b
0
2
+
2
a
1
2
(
2
e
y
+
b
0
)
2
)
=
a
1
(
1
-
24
k
2
(
cosh
⁡
(
y
)
+
sinh
⁡
(
y
)
)
b
0
-
3
k
2
b
0
2
+
2
a
1
2
(
2
cosh
⁡
(
y
)
+
2
sinh
⁡
(
y
)
+
b
0
)
2
)
,
where a
1, b
0, and k are real parameters and y
=
k
x
+
t
(
k
3
-
k
a
1
2
).Case 5 (whena
1
=
-
5, b
0
=
-
2, and k
=
3).
Consider(40)
u
5
(
x
,
t
)
=
-
5
(
1
+
432
e
-
48
t
+
3
x
-
108
+
50
(
-
2
+
2
e
-
48
t
+
3
x
)
2
)
=
-
5
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
1
+
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
432
(
cosh
⁡
(
48
t
-
3
x
)
+
sinh
⁡
(
48
t
-
3
x
)
)
h
h
h
h
h
h
h
h
h
h
×
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
-
108
+
50
(
(
48
t
-
3
x
)
-
2
+
2
cosh
⁡
(
48
t
-
3
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 3(a).Evolution of the solutions of (25). (a) Two-peak bright soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
-
2, and k
=
3. (b) A bright soliton in collision with a dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
3, b
0
=
-
1, and k
=
4.
(a)
(b)Case 6 (whena
1
=
3, b
0
=
-
1, and k
=
4).
Consider(41)
u
7
(
x
,
t
)
=
3
(
1
+
384
e
28
t
+
4
x
-
48
+
18
(
-
1
+
2
e
28
t
+
4
x
)
2
)
=
3
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
)
1
+
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
384
(
cosh
⁡
(
28
t
+
4
x
)
+
sinh
⁡
(
28
t
+
4
x
)
)
h
h
h
h
h
h
h
h
×
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
-
48
+
18
(
(
28
t
+
4
x
)
-
1
+
2
cosh
⁡
(
28
t
+
4
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 3(b).Case 7 (whena
1
=
-
5, b
0
=
2, and k
=
-
3).
Consider(42)
u
6
(
x
,
t
)
=
-
5
(
1
-
432
e
48
t
-
3
x
-
108
+
50
(
2
+
2
e
48
t
-
3
x
)
2
)
=
-
5
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
1
-
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
432
(
cosh
⁡
(
48
t
-
3
x
)
+
sinh
⁡
(
48
t
-
3
x
)
)
h
h
h
h
h
h
h
h
h
×
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
-
108
-
50
(
(
48
t
-
3
x
)
2
+
2
cosh
⁡
(
48
t
-
3
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 4(a).Evolution of the solutions of (25). (a) Dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
2, and k
=
-
3. (b) Two-peak dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
3, and k
=
4.
(a)
(b)Case 8 (whena
1
=
-
5, b
0
=
3, and k
=
4).
Consider(43)
u
8
(
x
,
t
)
=
-
5
(
1
-
1152
e
-
36
t
+
4
x
-
432
+
50
(
3
+
2
e
-
36
t
+
4
x
)
2
)
=
-
5
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
)
1
-
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
1152
(
cosh
⁡
(
36
t
-
4
x
)
-
sinh
⁡
(
36
t
-
4
x
)
)
h
h
h
h
h
h
h
h
h
×
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
-
432
+
50
(
(
36
t
-
4
x
)
3
+
2
cosh
⁡
(
36
t
-
4
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 4(b).


## 4.1.1. Whenα
=
0 and β
=
1

whenα
=
0 and β
=
1, (9) becomes u
t
+
u
x
x
x
=
0. Then we can obtain shock solitons. The solution formula is shown as follows:
(33)
u
(
x
,
t
)
=
a
1
(
1
+
8
e
-
k
3
t
+
k
x
b
0
)
=
a
1
+
8
a
1
cosh
⁡
(
k
3
t
-
k
x
)
b
0
-
8
a
1
sinh
⁡
(
k
3
t
-
k
x
)
b
0
,
where a
1, b
0, and k are real parameters.Case 1 (whena
1
=
1, b
0
=
2, k
=
3).
Consider(34)
u
1
(
x
,
t
)
=
1
+
4
e
-
27
t
+
3
x
=
1
+
4
cosh
⁡
(
27
t
-
3
x
)
-
4
sinh
⁡
(
27
t
-
3
x
)
.
Description is as shown in Figure 1(a).Evolution of the solutions of (25). (a) Shock soliton: parameters are α
=
0, β
=
1, a
1
=
1, b
0
=
2, and k
=
3. (b) A shock soliton in collision with a dark soliton: parameters are α
=
0, β
=
1, a
1
=
3, b
0
=
-
2, and k
=
4.
(a)
(b)Case 2 (whena
1
=
3, b
0
=
-
2, and k
=
4).
Consider(35)
u
2
(
x
,
t
)
=
3
(
1
-
4
e
-
64
t
+
4
x
)
=
3
-
12
cosh
⁡
(
64
t
-
4
x
)
+
12
sinh
⁡
(
64
t
-
4
x
)
.
Description is as shown in Figure 1(b).


## 4.1.2. Whenα
=
1 and β
=
1

whenα
=
1 and β
=
1, (9) becomes u
t
+
u
2
u
x
+
u
x
x
x
=
0. Then we can obtain one-peak bright solitons. The solution formula is shown as follows:
(36)
u
(
x
,
t
)
=
a
1
(
1
+
24
e
y
k
2
b
0
3
k
2
b
0
2
+
2
a
1
2
(
2
e
y
+
b
0
)
2
)
=
a
1
(
1
+
24
k
2
(
cosh
⁡
(
y
)
+
sinh
⁡
(
y
)
)
b
0
3
k
2
b
0
2
+
2
a
1
2
(
2
cosh
⁡
(
y
)
+
2
sinh
⁡
(
y
)
+
b
0
)
2
)
,
where a
1, b
0, and k are real parameters and y
=
k
x
+
t
(
-
k
3
-
k
a
1
2
).Case 3 (whena
1
=
1, b
0
=
2, and k
=
2).
Consider(37)
u
3
(
x
,
t
)
=
1
+
192
e
-
10
t
+
2
x
48
+
2
(
2
+
2
e
-
10
t
+
2
x
)
2
=
1
+
192
(
cosh
⁡
(
10
t
-
2
x
)
-
sinh
⁡
(
10
t
-
2
x
)
)
48
+
2
(
2
+
2
cosh
⁡
(
10
t
-
2
x
)
-
2
sinh
⁡
(
10
t
-
2
x
)
)
2
.
Description is as shown in Figure 2(a).Evolution of the solutions of (25). (a) Bright soliton: parameters are α
=
1, β
=
1, a
1
=
1, b
0
=
2, and k
=
2. (b) A bright soliton in collision with two dark solitons: parameters are α
=
1, β
=
1, a
1
=
1, b
0
=
-
1, and k
=
-
3.
(a)
(b)Case 4 (whena
1
=
1, b
0
=
-
1, and k
=
-
3).
Consider(38)
u
4
(
x
,
t
)
=
1
-
216
e
30
t
-
3
x
27
+
2
(
-
1
+
2
e
30
t
-
3
x
)
2
=
1
-
216
(
cosh
⁡
(
30
t
-
3
x
)
+
sinh
⁡
(
30
t
-
3
x
)
)
27
+
2
(
-
1
+
2
cosh
⁡
(
30
t
-
3
x
)
+
2
sinh
⁡
(
30
t
-
3
x
)
)
2
.
Description is as shown in Figure 2(b).


## 4.1.3. Whenα
=
1 and β
=
-
1

whenα
=
1 and β
=
-
1, (9) becomes u
t
+
u
2
u
x
-
u
x
x
x
=
0. Then we can obtain bright solitons and dark solitons. The solution formula is shown as follows:
(39)
u
(
x
,
t
)
=
a
1
(
1
-
24
e
y
k
2
b
0
-
3
k
2
b
0
2
+
2
a
1
2
(
2
e
y
+
b
0
)
2
)
=
a
1
(
1
-
24
k
2
(
cosh
⁡
(
y
)
+
sinh
⁡
(
y
)
)
b
0
-
3
k
2
b
0
2
+
2
a
1
2
(
2
cosh
⁡
(
y
)
+
2
sinh
⁡
(
y
)
+
b
0
)
2
)
,
where a
1, b
0, and k are real parameters and y
=
k
x
+
t
(
k
3
-
k
a
1
2
).Case 5 (whena
1
=
-
5, b
0
=
-
2, and k
=
3).
Consider(40)
u
5
(
x
,
t
)
=
-
5
(
1
+
432
e
-
48
t
+
3
x
-
108
+
50
(
-
2
+
2
e
-
48
t
+
3
x
)
2
)
=
-
5
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
1
+
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
432
(
cosh
⁡
(
48
t
-
3
x
)
+
sinh
⁡
(
48
t
-
3
x
)
)
h
h
h
h
h
h
h
h
h
h
×
(
+
2
cosh
⁡
(
48
t
-
3
x
)
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
-
108
+
50
(
(
48
t
-
3
x
)
-
2
+
2
cosh
⁡
(
48
t
-
3
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 3(a).Evolution of the solutions of (25). (a) Two-peak bright soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
-
2, and k
=
3. (b) A bright soliton in collision with a dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
3, b
0
=
-
1, and k
=
4.
(a)
(b)Case 6 (whena
1
=
3, b
0
=
-
1, and k
=
4).
Consider(41)
u
7
(
x
,
t
)
=
3
(
1
+
384
e
28
t
+
4
x
-
48
+
18
(
-
1
+
2
e
28
t
+
4
x
)
2
)
=
3
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
)
1
+
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
384
(
cosh
⁡
(
28
t
+
4
x
)
+
sinh
⁡
(
28
t
+
4
x
)
)
h
h
h
h
h
h
h
h
×
(
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
-
48
+
18
(
(
28
t
+
4
x
)
-
1
+
2
cosh
⁡
(
28
t
+
4
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
+
2
sinh
⁡
(
28
t
+
4
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 3(b).Case 7 (whena
1
=
-
5, b
0
=
2, and k
=
-
3).
Consider(42)
u
6
(
x
,
t
)
=
-
5
(
1
-
432
e
48
t
-
3
x
-
108
+
50
(
2
+
2
e
48
t
-
3
x
)
2
)
=
-
5
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
1
-
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
432
(
cosh
⁡
(
48
t
-
3
x
)
+
sinh
⁡
(
48
t
-
3
x
)
)
h
h
h
h
h
h
h
h
h
×
(
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
-
108
-
50
(
(
48
t
-
3
x
)
2
+
2
cosh
⁡
(
48
t
-
3
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
+
2
sinh
⁡
(
48
t
-
3
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 4(a).Evolution of the solutions of (25). (a) Dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
2, and k
=
-
3. (b) Two-peak dark soliton: parameters are α
=
1, β
=
-
1, a
1
=
-
5, b
0
=
3, and k
=
4.
(a)
(b)Case 8 (whena
1
=
-
5, b
0
=
3, and k
=
4).
Consider(43)
u
8
(
x
,
t
)
=
-
5
(
1
-
1152
e
-
36
t
+
4
x
-
432
+
50
(
3
+
2
e
-
36
t
+
4
x
)
2
)
=
-
5
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
)
1
-
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
1152
(
cosh
⁡
(
36
t
-
4
x
)
-
sinh
⁡
(
36
t
-
4
x
)
)
h
h
h
h
h
h
h
h
h
×
(
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
-
432
+
50
(
(
36
t
-
4
x
)
3
+
2
cosh
⁡
(
36
t
-
4
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
2
sinh
⁡
(
36
t
-
4
x
)
)
2
)
-
1
)
)
.
Description is as shown in Figure 4(b).


## 4.2. The Solutions from (32)

Based on the form of (32), if we set α
=
0 and β
=
0 at the same time, the solution should be zero. Therefore, it is not significative to study and analyse. So, in the following, we consider α and β under the circumstance that α
=
0 and β
=
0 are not zero in the meantime.Whenα
=
1 and β
=
1, (9) becomes u
t
+
u
2
u
x
+
u
x
x
x
=
0. Then we can obtain periodic wave solutions. The solution formula is shown as follows:
(44)
u
(
x
,
t
)
=
a
1
-
6
2
K
2
g
2
2
a
1
g
±
cos
⁡
(
y
)
(
-
a
1
+
g
2
(
3
K
2
-
2
a
1
2
)
)
,
where y
=
K
x
-
(
K
a
1
2
-
K
3
)
t and g
=
a
1
2
/
(
-
3
K
2
+
2
a
1
2
).To illustrate the analysis, we obtain the following results.Case 9 (whena
1
=
5, K
=
-
2).
Consider(45)
u
9
(
x
,
t
)
=
5
-
24
19
(
10
/
19
+
6
cos
⁡
(
2
(
21
t
-
x
)
)
)
.
Description is as shown in Figure 5.Figure 5
Evolution of the solutions of (32). Periodic wave solution: parameters are α
=
1, β
=
1, a
1
=
5, and K
=
-
2.


## 5. Conclusions

In this paper, our main attention has been focused on seeking the soliton solutions of (9) via Exp-function method. By applying Exp-function method, we have obtained shock soliton, bright soliton, two-peak bright soliton, dark soliton, two-peak dark soliton, and periodic wave solution. In addition, with figures and symbolic computations, we have described the propagation characteristics of those solitons under different values of those coefficients in the generalized MKdV equation. Furthermore, the problem solving process and the algorithm, by the help of Mathematica, can be easily extended to all kinds of nonlinear equations.



---

*Source: 153974-2014-05-15.xml*