Exact Solutions for a Class of Nonlinear PDE with Variable Coefficients Using ET and ETEM

Abstract

In this article, by using the modified CK direct method, we give a relationship between the generalized fifth-order KDV equations with variable coefficients and the corresponding constant coefficients ones. Then, we construct the abundant travelling solutions by the extended trial equation method (ETEM) in terms of different functions, such as the elliptic functions, rational functions, hyperbolic functions and trigonometric functions. The extended trial equation method is powerful and can be used to other partial differential equations and more research can be done by this method.

Share and Cite:

Sun, S. , Chang, L. and Liu, H. (2019) Exact Solutions for a Class of Nonlinear PDE with Variable Coefficients Using ET and ETEM. Journal of Applied Mathematics and Physics, 7, 2685-2700. doi: 10.4236/jamp.2019.711183.

1. Introduction

It is well known that the exact solutions of partial differential equations are an important problem in nonlinear science; the travelling wave solutions of partial differential equations with variable coefficients are always playing an important role in studying the long-time behavior of solutions and understanding the complex nonlinear fluctuations. Especially in the multifarious real physical background such as the field of nonlinear optical crystal and plasma, nonlinear partial differential equations (PDEs) with variable coefficients can often provide realistic and powerful models than the corresponding constant coefficients ones when the inhomogeneities of media are considered. There are lots of studies have been conducted with different types of the PDEs, such as the modified trigonometric functions series [1] [2], the ${G}^{\prime }/G$ expand method [3] [4], the first integral method [5] [6], the modified CK direct method [7] [8] [9], the $\mathrm{exp}\left(-\varphi \left(\xi \right)\right)$ method [10] [11], modified simple equation method [12] [13], infinite series method [14] [15], the Lie symmetry analysis method [16] - [21] and the extended trial equation method [22] [23]. In this paper, we will use the modified CK direct method and the extended trial equation method (ETEM) to discuss the exact solutions for the following generalized fifth-order nonlinear partial differential KDV equation [24] [25] [26] [27] with time-dependent variable coefficients of the dispersion:

$\frac{\partial }{\partial t}u=-\alpha \left(t\right){u}^{2}\frac{\partial }{\partial x}u-\beta \left(t\right)\frac{{\partial }^{2}u}{\partial {x}^{2}}\frac{\partial u}{\partial x}-F\left(t\right)u\frac{{\partial }^{3}u}{\partial {x}^{3}}-G\left(t\right)\frac{{\partial }^{5}u}{\partial {x}^{5}}.$ (1)

here the derivative ${u}_{t}$ represents the time evolution of a travelling wave as it travels in a certain direction, the function of nonlinear terms ${u}^{2}{u}_{x},{u}_{x}{u}_{xx}$ and $u{u}_{xxx}$ are to collect the waves and the function of the linear term ${u}_{xxxxx}$ is to disperse the waves, respectively, $\alpha \left(t\right),\beta \left(t\right),F\left(t\right),G\left(t\right)$ are arbitrary smooth functions of t. When $\alpha \left(t\right),\beta \left(t\right),F\left(t\right),G\left(t\right)$ be constant, the Equation (1) be the generalized fifth-order KDV equation. What’s more, we can get many other famous fifth-order partial differential equations by taking different values of the parameters:

Sawada-Kotera (SK) equation: when $\beta =F=5,\alpha =\frac{{F}^{2}}{2},G=1$

${u}_{t}+5{u}^{2}{u}_{x}+5{u}_{x}{u}_{xx}+5u{u}_{xxx}+{u}_{xxxxx}=0.$

Caudrey-Dodd-Gibbon (CDG) equation: when $\beta =F=30,\alpha =\frac{{F}^{2}}{5},G=1$

${u}_{t}+180{u}^{2}{u}_{x}+30{u}_{x}{u}_{xx}+30u{u}_{xxx}+{u}_{xxxxx}=0.$

Lax equation: when $\beta =2F,\alpha =\frac{3}{10}{F}^{2},G=1$

${u}_{t}+30{u}^{2}{u}_{x}+20{u}_{x}{u}_{xx}+10u{u}_{xxx}+{u}_{xxxxx}=0.$

Kaup-Kuperschmidt (KK) equation: when $F=10,\beta =\frac{5}{2}F,\alpha =\frac{{F}^{2}}{2},G=1$

${u}_{t}+50{u}^{2}{u}_{x}+25{u}_{x}{u}_{xx}+10u{u}_{xxx}+{u}_{xxxxx}=0.$

Ito equation: when $F=3,\beta =2F,\alpha =\frac{2{F}^{2}}{9},G=1$

${u}_{t}+2{u}^{2}{u}_{x}+6{u}_{x}{u}_{xx}+3u{u}_{xxx}+{u}_{xxxxx}=0.$

These equations can describe motions of waves in shallow water under gravity and in a one-dimensional nonlinear lattice. It is an important mathematical model with wide applications in quantum mechanics and nonlinear optics. Until now, a great number of authors have used lots of methods to solve these PDEs, Alvaro and H. Salas use the exp function method to obtain some exact solutions of Equation (1), Abdul-Majid Wazwaz obtained some new solition solutions of the modified form of fifth-order KDV equation (fKdV). In 2014, F. I. Leite considered the new concepts of self-adjoint equations formulated by Ibragimov and Gandarias are applied to the fifth-order evolution KDV equations. In 1989, Clarkson and Kruskal put forward the CK direct method for the first time, and then Lou modified the CK direct method and proposed a simple method which called the modified CK direct method. Furthermore, many scholars have studied the nonlinear equations with this method. In this paper, a new ETEM method is applied to the variable coefficient KDV equation which is reduced by CK method for the first time, and a series of new exact solutions are obtained.

The aim of this research is to establish an equivalence relation between the variable coefficients PDE and the corresponding constant coefficients one, and use a variety of distinct trail equations to construct some new travelling wave solutions of the constant coefficients equation. In this way, we can obtain a series of new exact solutions of the variable-coefficient equation by the relation between the variable-coefficient equation and the constant-coefficient equation, such as the trigonometric function solution, soliton solution, periodic solution and rational solution. In section 2, the equivalence relation between the fifth-order KDV equation with variable coefficient and the corresponding constant coefficient equation is given by using the modified CK direct method. In section 3, the ETEM is briefly introduced. In section 4, we obtain some new types of travelling wave solutions to Equation (1). In section 5, the more possibility of ETEM method for solving partial differential equations is discussed. In section 6, the conclusion and discussion are given.

2. Equivalence Transformations (ETs) of Generalized Fifth-Order KDV Equation

In this section, we will use the modified CK direct method to look for the equivalence transformations between Equation (1) and the corresponding following equation.

${u}_{t}+a{u}^{2}{u}_{x}+b{u}_{x}{u}_{xx}+cu{u}_{xxx}+d{u}_{xxxxx}=0,$ (2)

where $a,b,c,d$ are arbitrary constants. Suppose that Equation (1) has the following solution:

$u\left(x,t\right)=A+BU\left(X,T\right),$ (3)

where $A=A\left(x,t\right),B=B\left(x,t\right),X=X\left(x,t\right),T=T\left(x,t\right)$ are functions can be determined by requiring $U=U\left(X,T\right)$ satisfy the transformation $\left\{x,t,u\right\}\to \left\{X,T,U\right\}$. in other words, we let $\left\{X,T,U\right\}$ to satisfied the corresponding constant coefficients equation of Equation (1) as following

${U}_{T}+a{U}^{2}{U}_{X}+b{U}_{X}{U}_{XX}+cU{U}_{XXX}+d{U}_{XXXXX}=0.$ (4)

Then substituting Equation (3) into Equation (1), requiring $U=U\left(X,T\right)$ and ${U}_{T}$ satisfied Equation (4), we can collect the coefficients of U, letting their derivatives to be zero, we have

$GB{X}_{xxxxx}-dB{T}_{t}=0,$

$F{B}^{2}{X}_{x}^{3}-cB{T}_{t}=0,$

$\cdots$

$FB{X}_{x}^{3}A+15GB{X}_{x}{X}_{xx}^{2}+10GB{X}_{x}^{2}{X}_{xxx}=0,$ (5)

$10GB{X}_{xx}{X}_{xxx}+5GB{X}_{x}{X}_{xxxx}+3FB{X}_{x}{X}_{xx}a=0,$

$\beta B{X}_{xx}{A}_{x}+FB{X}_{xxx}A+GB{X}_{xxxxx}+B{X}_{t}=0,$

$G{A}_{xxxxx}+{A}_{x}\alpha {A}^{2}+\beta {A}_{xx}{A}_{x}+F{A}_{xxx}A+{A}_{t}=0.$

Solving above equations, we have

$A=0,\text{\hspace{0.17em}}B={C}_{0},\text{\hspace{0.17em}}F=\frac{c{T}_{t}}{{C}_{2}^{3}{C}_{0}},\text{\hspace{0.17em}}G=d\frac{{T}_{t}}{{C}_{2}^{5}},\text{\hspace{0.17em}}X={C}_{2}x+{C}_{3},$

$\alpha =a\frac{{T}_{t}}{{C}_{2}{C}_{0}^{2}},\text{\hspace{0.17em}}\beta =b\frac{{T}_{t}}{{C}_{2}^{3}{C}_{0}},$ (6)

where ${T}_{t}$ is an arbitrary function of t, ${C}_{1},{C}_{2},{C}_{3}$ are arbitrary constants. Then with Equation (6), we can obtain new travelling wave exact solutions for the Equation (1) as follows:

$u={C}_{0}U\left({C}_{2}x+{C}_{3},T\left(t\right)\right).$ (7)

If $U\left(X,T\right)$ is the solution of fifth-order KDV equation with constant coefficients, then $u\left(x,t\right)$ also is a solution of Equation (1).

3. Description of the Extend Trial Equation Method

In this section, we introduce the extend trial equation method (ETEM) as follows.

Step1. Firstly, we consider the following nonlinear partial differential equation:

${F}_{1}\left(u,{u}_{t},{u}_{x},{u}_{xx},{u}_{xt},{u}_{tt},{u}_{xxx},\cdots \right)=0,$ (8)

In order to look for the solutions of Equation (1), we make the travelling wave transformation

$u\left(x,t\right)=U\left(\xi \right),$

where $\xi =x-ct$, and c is an arbitrary constant.

Then the Equation (8) can be converted to an ordinary differential equation as following

${F}_{2}\left(U,{U}_{\xi },{U}_{\xi \xi },{U}_{\xi \xi \xi },\cdots \right)=0.$ (9)

Step 2. We suppose that the exact solution of Equation (9) is of the form

$u=\underset{i=1}{\overset{\delta }{\sum }}{\tau }_{i}{Y}^{i},$ (10)

where ${\tau }_{i}\left(i=1,2,3,\cdots ,\delta \right)$ are arbitrary constant, and ${Y}^{i}$ satisfied the following condition

${\left({Y}^{\prime }\right)}^{2}=\Lambda \left(Y\right)=\frac{\varphi \left(Y\right)}{\psi \left(Y\right)}=\frac{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+\cdots +{\xi }_{\theta }{Y}^{\theta }}{{\zeta }_{0}+{\zeta }_{1}Y+{\zeta }_{2}{Y}^{2}+\cdots +{\zeta }_{\epsilon }{Y}^{\epsilon }},$ (11)

where ${\xi }_{\theta }$ and ${\zeta }_{\epsilon }$ are arbitrary positive integers to be determined later.

Step 3. According to the balance principle we can determine a relation of $\theta$, $\epsilon$ and $\delta$. We can take some different values of $\theta$, $\epsilon$ and $\delta$ and the solution (10) can be written as a series of infinity solutions.

Step 4. Substituting Equation (10) and Equation (11) into Equation (9), letting the coefficients of Y to be zero, by solving the system of algebraic equations by some software like Maple and Mathematica, we can figure out the values of ${\xi }_{\theta },{\zeta }_{\epsilon }$ and ${\tau }_{\delta }$, then reducing Equation (11) to the elementary form as follows:

$±\left(\xi -{\eta }_{0}\right)=\int \frac{1}{\sqrt{\Lambda \left(Y\right)}}\text{d}Y=\int \sqrt{\frac{\psi \left(Y\right)}{\varphi \left(Y\right)}}\text{d}Y,$ (12)

where ${\eta }_{0}$ is an arbitrary constant.

Using a complete discrimination system for polynomial to classify the different forms of Equation (12), we can write the travelling wave solutions respectively. In this paper, the solutions of PDEs with variable coefficients also can be written respectively by the conclusion we got in the last section. Moreover, this method is appropriate for lots of other PDEs which can be discussed as follows.

4. ETEM to Generalized Fifth-Order KDV Equation

In this section, we discuss the generalized fifth-order KDV equation by using extend trial equation method.

${u}_{t}+a{u}^{2}{u}_{x}+b{u}_{x}{u}_{xx}+cu{u}_{xxx}+d{u}_{xxxxx}=0.$ (13)

Using the travelling wave transformation

$u\left(x,t\right)=U\left(\xi \right)$, where $\xi =x-Ct$,

Equation (13) turns into the following ordinary differential equation, then integrate the equation once:

${C}_{1}-CU+\frac{a}{3}{U}^{3}+\frac{b-c}{2}{\left({U}^{\prime }\right)}^{2}+cU{U}^{″}+d{U}^{\left(4\right)}=0,$ (14)

where ${C}_{1}$ is the integration constant. Substituting Equations (10) and (11) into Equation (14)

$U={\tau }_{\delta }{Y}^{\delta }+{\tau }_{\delta -1}{Y}^{\delta -1}+\cdots +{\tau }_{0},$

${U}^{3}={\tau }_{\delta }^{3}{Y}^{3\delta }+{\tau }_{\delta -1}^{3}{Y}^{3\delta -3}+\cdots +{\tau }_{0}^{3},$ (15)

${U}^{\prime }=\delta {\tau }_{\delta }{Y}^{\delta -1}{Y}^{\prime }+\cdots ={X}_{1}\left(\tau ,\theta ,\delta ,\epsilon \right){Y}^{\frac{2\delta +\theta -\epsilon -2}{2}}+\cdots ,$

${U}^{\left(4\right)}={X}_{2}\left(\tau ,\theta ,\delta ,\epsilon \right){Y}^{\delta +2\theta -2\epsilon -4}+\cdots .$ (16)

Using the balance principle, we find the highest degree terms of two Equations (15) and (16) should be the same, so we have

$\delta =\theta -\epsilon -2.$

By assuming and assigning the variables that satisfy the above conditions, we can obtain the traveling wave solution of the Equation (13) in many cases.

If $\epsilon =0,\theta =3$ and $\delta =1$, then

$U={\tau }_{0}+{\tau }_{1}Y,$

${U}^{\prime }={\tau }_{1}\sqrt{\frac{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}}{{\zeta }_{0}}},$

${U}^{″}={\tau }_{1}\frac{\left({\xi }_{1}+2{\xi }_{2}Y+3{\xi }_{3}{Y}^{2}\right)}{2{\zeta }_{0}},$

${U}^{‴}={\tau }_{1}\frac{\sqrt{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}}\left({\xi }_{2}+3{\xi }_{3}Y\right)}{{\zeta }_{0}^{\frac{3}{2}}},$

${U}^{\left(4\right)}={\tau }_{1}\frac{15{\xi }_{3}^{2}{Y}^{3}+15{\xi }_{2}{\xi }_{3}{Y}^{2}+9{\xi }_{1}{\xi }_{3}Y+2{\xi }_{2}^{2}Y+6{\xi }_{0}{\xi }_{3}+{\xi }_{1}{\xi }_{2}}{2{\zeta }_{0}^{2}},$ (17)

where ${\zeta }_{0}{\xi }_{3}\ne 0$.

Substituting Equation (17) into Equation (14), collecting the coefficients of Y and solving these algebraic equation systems respectively, we can obtain:

${\xi }_{3}={L}_{1}{\zeta }_{0}{\tau }_{1},$

${\xi }_{2}=-\frac{{\zeta }_{0}{\tau }_{0}\left(3{L}_{1}c+2a\right)}{{L}_{2}},$

${\xi }_{1}=\frac{{\zeta }_{0}\left(-9{L}_{1}^{2}d{c}^{2}{\tau }_{0}^{2}+3{L}_{1}{L}_{2}{c}^{2}{\tau }_{0}^{2}-12{L}_{1}adc{\tau }_{0}^{2}-{L}_{2}^{2}a{\tau }_{0}^{2}+2{L}_{2}ac{\tau }_{0}^{2}-4d{a}^{2}{\tau }_{0}^{2}+C{L}_{2}^{2}\right)}{{\tau }_{1}{L}_{2}^{2}\left(9d{L}_{1}+b\right)},$

${\xi }_{0}$ (See Appendixes A1), (18)

where ${L}_{1}=±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}$, ${L}_{3}=54{L}_{1}^{2}{d}^{2}+15{L}_{1}db-9{L}_{1}dc+{b}^{2}-bc$, ${L}_{2}=15d{L}_{1}+c+b$. We can use Equation (18) and Equation (12) to get

$±\left(\xi -{\eta }_{0}\right)=\int \frac{1}{\sqrt{\Lambda \left(Y\right)}}\text{d}Y=\int \sqrt{\frac{{\xi }_{3}{Y}^{3}+{\xi }_{2}{Y}^{2}+{\xi }_{1}Y+{\xi }_{0}Y}{{\zeta }_{0}}}\text{d}Y.$ (19)

Then, we will discuss some of the special form of solutions of Equation (14) as follow.

If Equation (19) can be written as following form:

$\begin{array}{l}{\xi }_{3}{Y}^{3}+{\xi }_{2}{Y}^{2}+{\xi }_{1}Y+{\xi }_{0}Y\\ ={L}_{1}{\tau }_{1}{Y}^{3}-\frac{{\tau }_{0}\left(3{L}_{1}c+2a\right)}{{L}_{2}}{Y}^{2}\\ +\frac{\left(-9{L}_{1}^{2}d{c}^{2}{\tau }_{0}^{2}+3{L}_{1}{L}_{2}{c}^{2}{\tau }_{0}^{2}-12{L}_{1}adc{\tau }_{0}^{2}-{L}_{2}^{2}a{\tau }_{0}^{2}+2{L}_{2}ac{\tau }_{0}^{2}-4d{a}^{2}{\tau }_{0}^{2}+C{L}_{2}^{2}\right)}{{\tau }_{1}{L}_{2}^{2}\left(9d{L}_{1}+b\right)}Y\\ \text{ }+2{\zeta }_{0}\left(-6d{\tau }_{0}^{3}{L}_{2}^{2}{a}^{2}-{L}_{2}^{3}ab{\tau }_{0}^{3}+3C{L}_{2}^{3}b{\tau }_{0}-27{C}_{1}{L}_{1}{L}_{2}^{3}d-9{c}^{3}{\tau }_{0}^{3}{L}_{2}^{2}{L}_{1}+3ac{\tau }_{0}^{3}{L}_{2}^{3}\end{array}$

$\begin{array}{l}\text{ }-6a{c}^{2}{\tau }_{0}^{3}{L}_{2}^{2}-3c{\tau }_{0}{L}_{2}^{3}C-81{d}^{2}{\tau }_{0}^{3}{L}_{1}^{3}{c}^{3}-24{d}^{2}{\tau }_{0}^{3}{a}^{3}-3{C}_{1}{L}_{2}^{3}b+6d{\tau }_{0}C{L}_{2}^{2}a\\ \text{ }-9{L}_{1}{L}_{2}^{3}da{\tau }_{0}^{3}+27C{L}_{1}{L}_{2}^{3}d{\tau }_{0}+54{c}^{3}{\tau }_{0}^{3}{L}_{2}{L}_{1}^{2}d+24c{\tau }_{0}^{3}{L}_{2}d{a}^{2}-162{d}^{2}{\tau }_{0}^{3}{L}_{1}^{2}{c}^{2}a\\ \text{ }-108{d}^{2}{\tau }_{0}^{3}{L}_{1}c{a}^{2}+72{c}^{2}{\tau }_{0}^{3}{L}_{1}{L}_{2}ad-9d{\tau }_{0}^{3}{L}_{1}{L}_{2}^{2}ac+9d{\tau }_{0}C{L}_{1}{L}_{2}^{2}c\right)/\left(3{L}_{2}^{3}{\tau }_{1}^{2}{L}_{3}\right)\\ ={\left(Y-{\alpha }_{1}\right)}^{3}.\end{array}$ (20)

where ${\alpha }_{1}$ is an arbitrary constant. Collecting the coefficients of Y of Equation (20), let all coefficients to be zero, we can get the following results

${\tau }_{1}=\frac{1}{{L}_{1}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\tau }_{0}=\frac{3{\alpha }_{1}{L}_{2}}{3c{L}_{1}+3a},$

$\begin{array}{l}C=3{\alpha }_{1}^{2}\left(135{L}_{1}^{3}d{c}^{2}-18{L}_{1}^{2}{L}_{2}{c}^{2}+180{L}_{1}^{2}adc+9{L}_{1}^{2}b{c}^{2}+6a{L}_{1}{L}_{2}^{2}-12ac{L}_{1}{L}_{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+60{L}_{1}d{a}^{2}+12{L}_{1}abc+4{a}^{2}b\right)/\left(2{L}_{1}\left(3{L}_{1}c+2a\right)\right),\end{array}$

${C}_{1}$ (See Appendixes A2) (21)

In this family, Equations (18), (19) and (21) lead to get

${\xi }_{3}={\zeta }_{0},\text{\hspace{0.17em}}{\xi }_{2}=-3{\alpha }_{1}{\zeta }_{0},\text{\hspace{0.17em}}{\xi }_{1}=3{\alpha }_{1}^{2}{\zeta }_{0},\text{\hspace{0.17em}}{\xi }_{0}=-{\alpha }_{1}^{3}{\zeta }_{0}.$ (22)

And

$±\left(\xi -{\eta }_{0}\right)=\int \frac{1}{\sqrt{\Lambda \left(Y\right)}}\text{d}Y=\int \frac{\text{d}Y}{{\left(Y-{\alpha }_{1}\right)}^{\frac{3}{2}}}=-\frac{2}{\sqrt{Y-{\alpha }_{1}}}.$ (23)

By Equation (23), (21), we can get the rational exact solution of Equation (14) as following

$U\left(\xi \right)=\frac{3{\alpha }_{1}{L}_{2}}{3c{L}_{1}+3a}+\frac{1}{{L}_{1}}\left[{\alpha }_{1}+\frac{4}{{\left(\xi -{\eta }_{0}\right)}^{2}}\right].$

$u\left(x,t\right)=\frac{3{\alpha }_{1}{L}_{2}}{3c{L}_{1}+3a}+\frac{1}{{L}_{1}}\left[{\alpha }_{1}+\frac{4}{{\left(x-Ct-{\eta }_{0}\right)}^{2}}\right],$ (24)

where

$\begin{array}{l}C=3{\alpha }_{1}^{2}\left(135{L}_{1}^{3}d{c}^{2}-18{L}_{1}^{2}{L}_{2}{c}^{2}+180{L}_{1}^{2}adc+9{L}_{1}^{2}b{c}^{2}+6a{L}_{1}{L}_{2}^{2}-12ac{L}_{1}{L}_{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+60{L}_{1}d{a}^{2}+12{L}_{1}abc+4{a}^{2}b\right)/\left(2{L}_{1}\left(3{L}_{1}c+2a\right)\right),\end{array}$

Hence, with the relation of Equation (7) obtained above, Equation (1) have the following form of rational travelling wave solution

$\begin{array}{l}u\left(x,t\right)={C}_{0}\frac{3{\alpha }_{1}\left(15d\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)+c+b\right)}{3c\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)+3a}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{C}_{0}\frac{1}{\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)}\left({\alpha }_{1}+\frac{4}{{\left(X-CT-{\eta }_{0}\right)}^{2}}\right),\end{array}$ (25)

where $X={C}_{2}x+{C}_{3}$, T is an arbitrary function of t (Figure 1 and Figure 2).

1) In the following, we chose some new coefficients to make it easier to calculate

${\xi }_{3}={L}_{1}{\zeta }_{0}{\tau }_{1},\text{\hspace{0.17em}}{\xi }_{2}={L}_{4}{\zeta }_{0}{\tau }_{0},$

${\xi }_{1}=\left(2\left(-{L}_{4}^{2}d{\tau }_{0}^{2}+{L}_{4}c{\tau }_{0}^{2}-a{\tau }_{0}^{2}+C\right){\zeta }_{0}\right)/\left({\tau }_{1}\left(9d{L}_{4}+b\right)\right),$

$\begin{array}{l}{\xi }_{0}=2{\zeta }_{0}\left(-9ad{L}_{1}{\tau }_{0}^{3}-ab{\tau }_{0}^{3}+27dC{L}_{1}{\tau }_{0}+3bC{\tau }_{0}-27d{L}_{1}{C}_{1}-3b{C}_{1}+6cd{L}_{4}^{2}{\tau }_{0}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-3{c}^{2}{L}_{4}{\tau }_{0}^{3}+3ac{\tau }_{0}^{3}-3cC{\tau }_{0}-3{d}^{2}{L}_{4}^{3}{\tau }_{0}^{3}-3ad{L}_{4}{\tau }_{0}^{3}+3dC{L}_{4}{\tau }_{0}\right)/\left(3{\tau }_{1}^{2}{L}_{3}\right).\end{array}$ (26)

Figure 1. The three-dimensional picture of the exact solution (25) of Equation (1), and its projection at $b=3,c=2,a=d=1,{\eta }_{0}=3$ and $\alpha =11$.

Figure 2. The two-dimensional picture of the equation exact solution (25) of Equation (1), and its projection at $t=0$. When $b=3,c=2,a=d=1,{\eta }_{0}=3$ and $\alpha =11$.

where

${L}_{1}=±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d},$

${L}_{3}=54{L}_{1}^{2}{d}^{2}+15{L}_{1}db-9{L}_{1}dc+{b}^{2}-bc,$

${L}_{4}=\left(3c{L}_{1}+2a\right)/\left(15d{L}_{1}+b+c\right).$

If Equation (19) can be written as following form:

$\begin{array}{l}{L}_{1}{\tau }_{1}{Y}^{3}+{L}_{4}{\tau }_{0}{Y}^{2}+\frac{2\left(-{L}_{4}^{2}d{\tau }_{0}^{2}+{L}_{4}c{\tau }_{0}^{2}-a{\tau }_{0}^{2}+C\right){\zeta }_{0}}{{\tau }_{1}\left(9d{L}_{4}+b\right)}Y+\frac{1}{3{\tau }_{1}^{2}{L}_{5}}\left(-9ad{L}_{1}{\tau }_{0}^{3}\\ \text{ }-ab{\tau }_{0}^{3}+27dC{L}_{1}{\tau }_{0}+3bC{\tau }_{0}-27d{L}_{1}{C}_{1}-3b{C}_{1}+6cd{L}_{4}^{2}{\tau }_{0}^{3}-3{c}^{2}{L}_{4}{\tau }_{0}^{3}+3ac{\tau }_{0}^{3}\\ -3cC{\tau }_{0}-3{d}^{2}{L}_{4}^{3}{\tau }_{0}^{3}-3ad{L}_{4}{\tau }_{0}^{3}+3dC{L}_{4}{\tau }_{0}\right)={\left(Y-{\alpha }_{2}\right)}^{2}\left(Y-{\alpha }_{3}\right),\end{array}$ (27)

where ${\alpha }_{2},{\alpha }_{3}$ are arbitrary constants. Collecting and letting all coefficients of Y of Equation (27) to be zero, we can get the following results

${\tau }_{1}=\frac{1}{{L}_{1}},\text{\hspace{0.17em}}{\tau }_{0}=\frac{2{\alpha }_{2}+{\alpha }_{3}}{{L}_{4}},$

$\begin{array}{l}C=\frac{1}{2{L}_{4}^{2}{L}_{1}}\left(17{L}_{1}{L}_{4}^{2}d{\alpha }_{2}^{2}+26{L}_{1}{L}_{4}^{2}d{\alpha }_{2}{\alpha }_{3}+2d{L}_{1}{L}_{4}^{2}{\alpha }_{3}^{2}-8{L}_{1}{L}_{4}c{\alpha }_{1}^{2}-8{L}_{1}{L}_{4}c{\alpha }_{2}{\alpha }_{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-2{L}_{1}{L}_{4}c{\alpha }_{4}^{2}+{L}_{4}^{2}b{\alpha }_{1}^{2}+2{L}_{4}^{2}b{\alpha }_{2}{\alpha }_{3}+8{L}_{1}a{\alpha }_{2}^{2}+8{L}_{1}a{\alpha }_{2}{\alpha }_{3}+2{L}_{1}a{\alpha }_{3}^{2}\right),\end{array}$

${C}_{1}$ (See Appendixes A3) (28)

In this family, Equation (18), (19) and (28) lead to get as follows:

${\xi }_{3}={\zeta }_{0},\text{\hspace{0.17em}}{\xi }_{2}=-\left(2{\alpha }_{2}+{\alpha }_{3}\right){\zeta }_{0},\text{\hspace{0.17em}}{\xi }_{1}={\alpha }_{2}\left({\alpha }_{2}+2{\alpha }_{3}\right){\zeta }_{0},\text{\hspace{0.17em}}{\xi }_{0}=-{\alpha }_{3}{\alpha }_{2}^{2}{\zeta }_{0},$ (29)

where ${\zeta }_{0}$ is an arbitrary constant.

When ${\alpha }_{3}-{\alpha }_{2}>0$,

$±\left(\xi -{\eta }_{0}\right)=\int \frac{\text{d}Y}{\left(Y-{\alpha }_{2}\right){\left(Y-{\alpha }_{3}\right)}^{\frac{1}{2}}}=-\frac{2}{\sqrt{{\alpha }_{3}-{\alpha }_{2}}}{\mathrm{tan}}^{-1}\left[\frac{\sqrt{Y-{\alpha }_{3}}}{\sqrt{{\alpha }_{3}-{\alpha }_{2}}}\right],$ (30)

we can solve for Y

$Y={\alpha }_{3}+\left({\alpha }_{3}-{\alpha }_{2}\right){\mathrm{tan}}^{2}\left[\frac{\sqrt{{\alpha }_{3}-{\alpha }_{2}}}{2}\left(\xi -{\eta }_{0}\right)\right].$ (31)

By Equation (31), (29), we can get the rational exact solution of Equation (14) as following

$U\left(\xi \right)=\frac{2{\alpha }_{2}+{\alpha }_{3}}{{L}_{4}}+\frac{1}{{L}_{1}}\left[{\alpha }_{3}+\left({\alpha }_{3}-{\alpha }_{2}\right){\mathrm{tan}}^{2}\left[\frac{\sqrt{{\alpha }_{3}-{\alpha }_{2}}}{2}\left(\xi -{\eta }_{0}\right)\right]\right].$ (32)

With the relation of Equation (7), Equation (1) have the following form of rational travelling wave solution

$\begin{array}{c}u\left(x,t\right)={C}_{0}\left(\frac{\left(2{\alpha }_{2}+{\alpha }_{3}\right)\left(15d\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)+b+c\right)}{3c\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)+2a}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)}\\ \begin{array}{c}\begin{array}{l}\\ \end{array}\\ \\ \\ \end{array}×\left[{\alpha }_{3}+\left({\alpha }_{3}-{\alpha }_{2}\right){\mathrm{tan}}^{2}\left[\frac{\sqrt{{\alpha }_{3}-{\alpha }_{2}}}{2}\left(X-CT-{\eta }_{0}\right)\right]\right]\right),\end{array}$ (33)

where $X={C}_{2}x+{C}_{3}$, T is an arbitrary function of t. (Figure 3 and Figure 4)

When ${\alpha }_{3}-{\alpha }_{2}<0$,

$Y={\alpha }_{2}+\left({\alpha }_{2}-{\alpha }_{3}\right){\text{csch}}^{2}\left[\frac{\sqrt{{\alpha }_{3}-{\alpha }_{2}}}{2}\left(\xi -{\eta }_{0}\right)\right].$ (34)

Similarly, we can get the exact solution of Equation (1) as follow

Figure 3. The three-dimensional picture of the exact solution (33) of Equation (1), When $b=3,c=2,{\eta }_{0}=1,a=d=1$ and ${\alpha }_{2}=1,{\alpha }_{3}=5$.

Figure 4. The two-dimensional picture of the equation exact solution (33) of Equation (1), and its projection at $t=0$. When $b=3,c=2,{\eta }_{0}=1,a=d=1$ and ${\alpha }_{2}=1,{\alpha }_{3}=5$.

$\begin{array}{c}u\left(x,t\right)={C}_{0}\left(\frac{\left(2{\alpha }_{2}+{\alpha }_{3}\right)\left(15d\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)+b+c\right)}{3c\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)+2a}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{\left(±\frac{-2c-b+\sqrt{-40ad+{b}^{2}+4bc+4{c}^{2}}}{30d}\right)}\\ \begin{array}{c}\begin{array}{l}\\ \end{array}\\ \\ \\ \end{array}×\left[{\alpha }_{2}+\left({\alpha }_{2}-{\alpha }_{3}\right){\text{csch}}^{2}\left[\frac{\sqrt{{\alpha }_{2}-{\alpha }_{3}}}{2}\left(X-CT-{\eta }_{0}\right)\right]\right]\right),\end{array}$ (35)

where $X={C}_{2}x+{C}_{3}$, T is an arbitrary function of t (Figure 5 and Figure 6).

Figure 5. The three-dimensional picture of the exact solution (35) of Equation (1), when $b=3,c=2,a=d=1,{\eta }_{0}=1$ and ${\alpha }_{2}=3,{\alpha }_{3}=1$.

Figure 6. The two-dimensional picture of the equation exact solution (35) of Equation (1), and its projection at $t=0$, when $b=3,c=2,a=d=1,{\eta }_{0}=1$ and ${\alpha }_{2}=3,{\alpha }_{3}=1$.

5. More Discussion

In this article, we obtained a series of exact travelling wave solutions of fifth-order KDV equation by the extended trail equation method; according to picking different parameters we can get more exact analytic solutions of nonlinear partial differential equations like fifth-order KDV equation.

The extend trial equation method (ETEM) is proving to play an important role in solving partial differential equations, by using a variety of trail equations, we can construct lots of new types of travelling wave solutions. In this paper, we only considered the following parameters

$\epsilon =0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\theta =3,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\delta =1.$

In a later study we can also study different situations such as $\delta \ne \text{1}$ or $\epsilon \ne 0$. In this paper, the highest degree of Equations (15) and (16) should be the same, according to balance the highest degree terms of these two Equations (15) and (16). So we can take any other parameters that satisfies the following equation

$\delta =\theta -\epsilon -2.$

For example

1) When $\epsilon =0,\theta =4$ and $\delta =2$.

${\left({Y}^{\prime }\right)}^{2}=\frac{\Phi \left(Y\right)}{\Psi \left(Y\right)}=\frac{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}}{{\zeta }_{0}},$ (36)

and

$U={\tau }_{0}+{\tau }_{1}Y+{\tau }_{2}{Y}^{2},$

${U}^{\prime }=\frac{\left({\tau }_{1}+2{\tau }_{2}Y\right)\sqrt{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}}}{\sqrt{{\zeta }_{0}}},$

$\begin{array}{l}{U}^{″}=\\ \frac{12{Y}^{4}{\tau }_{2}{\xi }_{4}+4{Y}^{3}{\tau }_{1}{\xi }_{4}+10{Y}^{3}{\tau }_{2}{\xi }_{3}+3{Y}^{2}{\tau }_{1}{\xi }_{3}+8{Y}^{2}{\tau }_{2}{\xi }_{3}+2Y{\tau }_{1}{\xi }_{2}+6Y{\tau }_{2}{\xi }_{1}+{\tau }_{1}{\xi }_{1}+4{\tau }_{2}{\xi }_{0}}{2{\zeta }_{0}}.\end{array}$

$\begin{array}{l}{U}^{‴}=\frac{1}{\sqrt{{\zeta }_{0}^{3}}}\sqrt{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}}\left(24{Y}^{3}{\tau }_{2}{\xi }_{4}+6{Y}^{2}{\tau }_{1}{\xi }_{4}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+15{Y}^{2}{\tau }_{2}{\xi }_{3}+3Y{\tau }_{1}{\xi }_{3}+8Y{\tau }_{2}{\xi }_{3}+{\tau }_{1}{\xi }_{2}+3{\tau }_{2}{\xi }_{1}\right),\end{array}$ (37)

where ${\zeta }_{0}{\xi }_{4}\ne 0$.

2) When $\epsilon =1,\theta =4$ and $\delta =1$.

${\left({Y}^{\prime }\right)}^{2}=\frac{\Phi \left(Y\right)}{\Psi \left(Y\right)}=\frac{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}}{{\zeta }_{0}+{\zeta }_{1}Y},$ (38)

then

$U={\tau }_{0}+{\tau }_{1}Y,$

${U}^{\prime }={\tau }_{1}\frac{\sqrt{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}}}{\sqrt{{\zeta }_{0}+{\zeta }_{1}Y}},$

${U}^{″}={\tau }_{1}\left(\frac{{\xi }_{1}+2{\xi }_{2}Y+3{\xi }_{3}{Y}^{2}+4{\xi }_{4}{Y}^{3}}{2{\left({\zeta }_{0}+{\zeta }_{1}Y\right)}^{2}}-\frac{{\zeta }_{1}\left({\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}\right)}{2{\left({\zeta }_{0}+{\zeta }_{1}Y\right)}^{3}}\right),$

$\begin{array}{c}{U}^{‴}={\tau }_{1}\left(\frac{1}{2}{\left({\zeta }_{0}+{\zeta }_{1}Y\right)}^{-\frac{5}{2}}\left(12{\xi }_{4}{Y}^{2}+6{\xi }_{3}Y+2{\xi }_{2}\right)\sqrt{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{3}{2}{\zeta }_{1}{\left({\zeta }_{0}+{\zeta }_{1}Y\right)}^{-\frac{7}{2}}\left({\xi }_{1}+2{\xi }_{2}Y+3{\xi }_{3}{Y}^{2}+4{\xi }_{4}{Y}^{3}\right)\sqrt{{\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{3}{2}{\zeta }_{1}^{2}{\left({\zeta }_{0}+{\zeta }_{1}Y\right)}^{-\frac{9}{2}}{\left({\xi }_{0}+{\xi }_{1}Y+{\xi }_{2}{Y}^{2}+{\xi }_{3}{Y}^{3}+{\xi }_{4}{Y}^{4}\right)}^{\frac{3}{2}}\right),\end{array}$ (39)

where ${\xi }_{\theta }$ and ${\zeta }_{\epsilon }$ are arbitrary positive integers to be determined in later calculations.

In later studies, when studying other nonlinear partial differential equations, we can also obtain the relationship of parameters through balancing the highest degree terms of these equations, select parameters and apply the Extend trial equation method (ETEM) to solve the exact solution of the equation.

6. Conclusion

In this letter, the ETEM has been successfully applied to construct exact traveling wave solutions for fifth-order KDV equation. Then, the solutions of corresponding nonlinear partial differential equations with variable coefficients are obtained by the equivalence transformation given in Section 2. In later studies, many solutions of variable coefficient PDEs can be considered in the same procedure. Generally, for tackling exact solutions to vc-PDEs are difficult, the results in this paper provide a useful supplement to the existing literature. Moreover, the equivalence transformation and improved ETEM can be used to other types of vc-PDEs in mathematical physics.

Appendixes

$\begin{array}{c}{\xi }_{0}=2{\zeta }_{0}\left(-6d{\tau }_{0}^{3}{L}_{2}^{2}{a}^{2}-{L}_{2}^{3}ab{\tau }_{0}^{3}+3C{L}_{2}^{3}b{\tau }_{0}-27{C}_{1}{L}_{1}{L}_{2}^{3}d-9{c}^{3}{\tau }_{0}^{3}{L}_{2}^{2}{L}_{1}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+3ac{\tau }_{0}^{3}{L}_{2}^{3}-6a{c}^{2}{\tau }_{0}^{3}{L}_{2}^{2}-3c{\tau }_{0}{L}_{2}^{3}C-81{d}^{2}{\tau }_{0}^{3}{L}_{1}^{3}{c}^{3}-24{d}^{2}{\tau }_{0}^{3}{a}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-3{C}_{1}{L}_{2}^{3}b+6d{\tau }_{0}C{L}_{2}^{2}a-9{L}_{1}{L}_{2}^{3}da{\tau }_{0}^{3}+27C{L}_{1}{L}_{2}^{3}d{\tau }_{0}+54{c}^{3}{\tau }_{0}^{3}{L}_{2}{L}_{1}^{2}d\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+24c{\tau }_{0}^{3}{L}_{2}d{a}^{2}-162{d}^{2}{\tau }_{0}^{3}{L}_{1}^{2}{c}^{2}a-108{d}^{2}{\tau }_{0}^{3}{L}_{1}c{a}^{2}+72{c}^{2}{\tau }_{0}^{3}{L}_{1}{L}_{2}ad\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-9d{\tau }_{0}^{3}{L}_{1}{L}_{2}^{2}ac+9d{\tau }_{0}C{L}_{1}{L}_{2}^{2}c\right)/\left(3{L}_{2}^{3}{\tau }_{1}^{2}{L}_{3}\right).\end{array}$ (A1)

$\begin{array}{c}{C}_{1}={\alpha }_{1}^{3}\left(10935{L}_{1}^{5}{L}_{2}{d}^{2}{c}^{2}+2187{L}_{1}^{5}{d}^{2}{c}^{3}-1458{L}_{2}^{4}{L}_{2}^{2}d{c}^{2}+14580{L}_{1}^{4}{L}_{2}{d}^{2}ac\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+1944{L}_{1}^{4}{L}_{2}db{c}^{2}-729{L}_{1}^{4}{L}_{2}d{c}^{3}+4374{L}_{1}^{4}{d}^{2}a{c}^{2}+243{L}_{1}^{4}db{c}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+324{L}_{1}^{3}{L}_{2}^{3}da-972{L}_{1}^{3}{L}_{2}^{2}dac-162{L}_{1}^{3}{L}_{2}^{2}b{c}^{2}+4860{L}_{1}^{3}{L}_{2}{d}^{2}{a}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+2592{L}_{1}^{3}{L}_{2}abcd-972{L}_{1}^{3}{L}_{2}ad{c}^{2}+81{L}_{1}^{3}{L}_{2}{b}^{2}{c}^{2}-81{L}_{1}^{3}{L}_{2}b{c}^{3}\end{array}$

$\begin{array}{c}\text{ }+2916{L}_{1}^{3}{a}^{2}{d}^{2}c+486{L}_{1}^{3}adb{c}^{2}+27{L}_{1}^{3}{L}_{3}{c}^{3}+36{L}_{1}^{2}{L}_{2}^{3}ab-108{L}_{1}^{2}{L}_{2}^{2}abc\\ \text{ }+864{L}_{1}^{2}{L}_{2}{a}^{2}bd-324{L}_{1}^{2}{L}_{2}{a}^{2}dc+108{L}_{1}^{2}{L}_{2}a{b}^{2}c-108{L}_{1}^{2}{L}_{2}ab{c}^{2}\\ \text{ }+648{L}_{1}^{2}{a}^{3}{d}^{2}+324{L}_{1}^{2}{a}^{2}bcd+54{L}_{1}^{2}{L}_{3}a{c}^{2}+36{L}_{1}{L}_{2}{a}^{2}{b}^{2}-36{L}_{1}{L}_{2}{a}^{2}bc\\ \text{ }+72{L}_{1}d{a}^{3}b+36{L}_{1}{L}_{3}{a}^{2}c+8{L}_{3}{a}^{3}\right)/\left(2{L}_{1}^{2}{\left(3{L}_{1}c+2a\right)}^{3}\left(9d{L}_{1}+b\right)\right).\end{array}$ (A2)

$\begin{array}{c}{C}_{1}=\frac{1}{6{L}_{1}^{2}{L}_{4}^{3}\left(9d{L}_{1}+b\right)}\left(156{L}_{1}^{2}d{L}_{4}^{2}b{\alpha }_{2}^{3}+432{L}_{1}^{3}ad{\alpha }_{2}^{2}{\alpha }_{3}+135{L}_{4}^{3}{d}^{2}{L}_{1}^{2}{\alpha }_{2}^{2}{\alpha }_{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+48{L}_{1}^{2}ab{\alpha }_{2}^{2}{\alpha }_{3}+1863{L}_{1}^{3}{d}^{2}{L}_{4}^{2}{\alpha }_{2}^{2}{\alpha }_{3}+810{\alpha }_{2}{L}_{1}^{3}{d}^{2}{L}_{4}^{2}{\alpha }_{3}^{2}+216{\alpha }_{2}{L}_{1}^{3}ad{\alpha }_{3}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+24{\alpha }_{2}{L}_{1}^{2}ab{\alpha }_{3}^{2}+54{L}_{4}^{3}{d}^{2}{L}_{1}^{2}{\alpha }_{3}^{2}{\alpha }_{2}+6{L}_{1}^{2}{L}_{4}^{2}bd{\alpha }_{3}^{3}+6{L}_{1}{L}_{4}^{3}db{\alpha }_{2}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+15{L}_{1}{L}_{4}^{2}{b}^{2}{\alpha }_{2}^{2}{\alpha }_{3}+6{L}_{1}{L}_{4}^{2}{b}^{2}{\alpha }_{3}^{2}{\alpha }_{2}-432{L}_{1}^{3}dc{L}_{4}{\alpha }_{2}^{3}-54{L}_{4}^{2}dc{L}_{1}^{2}{\alpha }_{2}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-48{L}_{1}^{2}{L}_{4}bc{\alpha }_{2}^{3}-54{L}_{1}^{3}dc{L}_{4}{\alpha }_{3}^{3}-6{L}_{1}^{2}{L}_{4}bc{\alpha }_{3}^{3}-6{L}_{1}{L}_{4}^{2}bc{\alpha }_{2}^{3}+3{\alpha }_{2}^{2}{\alpha }_{3}{L}_{4}^{3}{L}_{3}\end{array}$

$\begin{array}{c}\text{ }+6{L}_{1}{L}_{4}^{2}{b}^{2}{\alpha }_{2}^{3}+918{L}_{1}^{3}{d}^{2}{L}_{4}^{2}{\alpha }_{2}^{3}+288{L}_{1}^{3}ad{\alpha }_{2}^{3}+54{L}_{4}^{3}{d}^{2}{L}_{1}^{2}{\alpha }_{2}^{3}+32{L}_{1}^{2}ab{\alpha }_{2}^{3}\\ \text{ }+54{L}_{1}^{3}{d}^{2}{L}_{4}^{2}{\alpha }_{3}^{3}+36{L}_{1}^{3}ad{\alpha }_{3}^{3}+4{L}_{1}^{2}ab{\alpha }_{3}^{3}+144{\alpha }_{2}{L}_{1}^{2}{L}_{4}^{2}bd{\alpha }_{3}^{2}\\ \text{ }+342{L}_{1}^{2}bd{L}_{1}^{2}{L}_{4}^{2}{\alpha }_{2}^{2}{\alpha }_{3}-36{\alpha }_{2}{L}_{1}^{2}{L}_{4}bc{\alpha }_{3}^{2}-648{L}_{1}^{3}dc{L}_{4}{\alpha }_{2}^{2}{\alpha }_{3}\\ \text{ }-135{L}_{4}^{2}dc{L}_{1}^{2}{\alpha }_{2}^{2}{\alpha }_{3}-72{L}_{1}^{2}{L}_{4}bc{\alpha }_{2}^{2}{\alpha }_{3}-54{L}_{4}^{2}dc{L}_{1}^{2}{\alpha }_{3}^{2}{\alpha }_{2}-15{L}_{1}{L}_{4}^{2}bc{\alpha }_{2}^{2}{\alpha }_{3}\\ \text{ }-6{L}_{1}{L}_{4}^{2}bc{\alpha }_{3}^{2}{\alpha }_{2}-324{\alpha }_{3}^{2}{\alpha }_{2}{L}_{1}^{3}{L}_{4}d+15{L}_{1}{L}_{4}^{3}bd{\alpha }_{2}^{2}{\alpha }_{3}+6{L}_{1}{L}_{4}^{3}bd{\alpha }_{3}^{2}{\alpha }_{2}\right).\end{array}$ (A3)

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

 [1] Qureshi, M.I., Quraishi, K.A. and Srivastava, H.M. (2008) Some Hypergeometric Summation Formulas and Series Identities Associated with Exponential and Trigonometric Functions. Integral Transforms & Special Functions, 19, 267-276.https://doi.org/10.1080/10652460801896024 [2] Yu, H. (2016) On Generalized Trigonometric Functions and Series of Rational Functions. Journal of Number Theory, 180, 512-532. https://doi.org/10.1016/j.jnt.2017.05.015 [3] Zhang, Z., Zhong, J., Dou, S., Liu, J., Peng, D. and Gao, T. (2013) A New Method to Construct Travelling Wave Solutions for the Klein-Gord-Zakharov Equations. Romanian Journal of Physics, 58, 766-777. [4] Miao, X. and Zhang, X. (2011) The Modified (G’/G)-Expansion Method and Traveling Wave Solutions of Nonlinear the Perturbed Nonlinear Schrödinger’s Equation with Kerr Law Nonlinearity. Communications in Nonlinear Science and Numerical Simulation, 16, 4259-4267. https://doi.org/10.1016/j.cnsns.2011.03.032 [5] Feng, Z. and Wang, X. (2003) The First Integral Method to the Two-Dimensional Burgers-Korteweg-de Vries Equation. Physics Letters A, 308, 173-178.https://doi.org/10.1016/S0375-9601(03)00016-1 [6] Eslami, M. and Rezazadeh, H. (2015) The First Integral Method for Wu-Zhang System with Conformable Time-Fractional Derivative. Calcolo, 53, 475-485.https://doi.org/10.1007/s10092-015-0158-8 [7] Qian, S. and Tian, L. (2007) Modification of the Clarkson-Kruskal Direct Method for a Coupled System. Chinese Physics Letters, 24, 2720-2723.https://doi.org/10.1088/0256-307X/24/10/002 [8] Chen, M. and Liu, X.Q. (2011) Exact Solutions and Conservation Laws of the Konopelchenko-Dubrovsky Equations. Pure & Applied Mathematics, 27, 533-532. [9] Yuan, Q. and Lei, Y. (2011) Symmetry Groups and New Exact Solutions of (2+1)-Dimensional Broer-Kaup-Kupershmidt (BKK) System. International Conference on Multimedia Technology. [10] He, J. and Abdou, M. (2007) New Periodic Solutions for Nonlinear Evolution Equations Using Exp-Function Method. Chaos, Solitons & Fractals, 34, 1421-1429.https://doi.org/10.1016/j.chaos.2006.05.072 [11] Kudryashov, N.A. and Loguinova, N.B. (2009) Be Careful with the Exp-Function Method. Communications in Nonlinear Science & Numerical Simulation, 14, 1881-1890. https://doi.org/10.1016/j.cnsns.2008.07.021 [12] Hossain, A.K.M.K.S., Akbar, M.A. and Wazwaz, A.M. (2017) Closed Form Solutions of Complex Wave Equations via the Modified Simple Equation Method. Cogent Physics, 4, Article ID: 1312751. https://doi.org/10.1080/23311940.2017.1312751 [13] Akter, J. and Ali Akbar, M. (2015) Exact Solutions to the Benney-Luke Equation and the Phi-4 Equations by Using Modified Simple Equation Method. Results in Physics, 5, 125-130. https://doi.org/10.1016/j.rinp.2015.01.008 [14] Jiao, X.Y. and Lou, S.Y. (2009) Approximate Direct Reduction Method: Infinite Series Reductions to the Perturbed mKdV Equation. Chinese Physics B, 18, 3611-3615. https://doi.org/10.1088/1674-1056/18/9/001 [15] Taghizadeh, N., Mirzazadeh, M. and Farahrooz, F. (2011) Exact Travelling Wave Solutions of the Coupled Klein-Gordon Equation by the Infinite Series Method. Applications and Applied Mathematics, 6, 1964-1972. [16] Liu, H.Z., Li, J.B. and Liu, L. (2010) Lie Symmetry Analysis, Optimal Systems and Exact Solutions to the Fifth-Order KdV Types of Equations. Journal of Mathematical Analysis and Applications, 368, 551-558.https://doi.org/10.1016/j.jmaa.2010.03.026 [17] Wang, G.W., Liu, X.Q. and Zhang, Y.Y. (2013) Symmetry Reduction, Exact Solutions and Conservation Laws of a New Fifth-Order Nonlinear Integrable Equation. Communications in Nonlinear Science and Numerical Simulation, 18, 2313-2320.https://doi.org/10.1016/j.cnsns.2012.12.003 [18] Liu, H.Z. (2015) Painlevé Test, Generalized Symmetries, Bäcklund Transformations and Exact Solutions to the Third-Order Burgers’ Equations. Journal of Statistical Physics, 158, 433-446. https://doi.org/10.1007/s10955-014-1130-8 [19] Liu, H.Z. and Li, J.B. (2012) Painlevé Analysis, Complete Lie Group Classifications and Exact Solutions to the Time-Dependent Coefficients Gardner Types of Equations. Nonlinear Dynamics, 80, 515-527. https://doi.org/10.1007/s11071-014-1885-0 [20] Sahoo, S. and Ray, S.S. (2017) Analysis of Lie Symmetries with Conservation Laws for the (3+1) Dimensional Time-Fractional mKdV–ZK Equation in Ion-Acoustic Waves. Nonlinear Dynamics, 90, 1105-1113. https://doi.org/10.1007/s11071-017-3712-x [21] Xin, X.P. (2018) Non-Local Symmetries and Exact Solutions of Nonlinear Development Equations. Journal of Liaocheng University (Natural Science), 32, 15-20. (In Chinese) [22] Ekici, M., Mirzazadeh, M., Sonmezoglu, A., et al. (2017) Nematicons in Liquid Crystals by Extended Trial Equation Method. Journal of Nonlinear Optical Physics & Materials, 26, Article ID: 1750005. https://doi.org/10.1142/S0218863517500059 [23] Ekici, M., Mirzazadeh, M., Sonmezoglu, A., et al. (2017) Optical Solitons with Anti-Cubic Nonlinearity by Extended Trial Equation Method. Optik—International Journal for Light and Electron Optics, 136, 368-373.https://doi.org/10.1016/j.ijleo.2017.02.004 [24] Abbasbandy, S. and Zakaria, F.S. (2008) Soliton Solutions for the Fifth-Order KdV Equation with the Homotopy Analysis Method. Nonlinear Dynamics, 51, 83-87.https://doi.org/10.1007/s11071-006-9193-y [25] Bridges, T.J., Derks, G. and Gottwald, G. (2002) Stability and Instability of Solitary Waves of the Fifth-Order KdV Equation: A Numerical Framework. Physica D—Nonlinear Phenomena, 172, 190-216. https://doi.org/10.1016/S0167-2789(02)00655-3 [26] Leite Freire, I. and Santos Sampaio, J.C. (2012) Corrigendum: Nonlinear Self-Adjointness of a Generalized Fifth-Order KdV Equation. Journal of Physics A Mathematical & Theoretical, 45, Article ID: 119502. https://doi.org/10.1088/1751-8113/45/11/119502 [27] Triki, H. and Biswas, A. (2011) Soliton Solutions for a Generalized Fifth-Order KdV Equation with T-Dependent Coefficients. Waves in Random & Complex Media, 21, 151-160. https://doi.org/10.1080/17455030.2010.539632