Exact Traveling Wave Solutions to Phi-4 Equation and Joseph-Egri (TRLW) Equation and Calogro-Degasperis (CD) Equation by Modified (G'/G²)-Expansion Method

Maha Al-Harbi, Waleed Al-Hamdan, Luwai Wazzan

Department of Mathematics, King Abdulaziz University, Jeddah, Saudi Arabia.

DOI: 10.4236/jamp.2023.117133                

Abstract

In this study, we will introduce the modified (G'/G²)-expansion method to explore some of the exact traveling wave solutions of some nonlinear partial differential equations namely, Phi-4 equation, Joseph-Egri (TRLW) equation, and Calogro-Degasperis (CD) equation. As a result, we have obtained solutions for the equations expressed in terms of trigonometric, hyperbolic and rational functions. Moreover, some selected solutions are plotted using some specific values for the parameters.

Keywords

Exact Solutions, Modified (G'/G²)-Expansion Method, Phi-4 Equation, Joseph-Egri (TRLW) Equation, Calogro-Degasperis (CD) Equation

1. Introduction

The world around us is inherently nonlinear. A large part of equations is nonlinear partial differential equations (NPDEs), which can be utilized to describe many important and sophisticated phenomena in physics, mathematical physics, engineering, and other various scientific fields. It is not easy to solve these equations, but exploration of their exact solutions becomes quite prominent due to the efficiency and reliability of symbolic software packages such as Maple or Mathematica.

In the past decades, authors constructed and developed various powerful and efficient methods for finding exact solutions of (NPDEs) with aid of symbolic software such as, the tanh-function method [1] [2] , sine-cosine function method [3] [4] , Hirota bilinear method [5] [6] , homogeneous balance method [7] [8] , the exp-function method [9] [10] , the sub-equation method, F-expansion method [11] [12] , (1/G')-expansion method [13] , (G'/G, 1/G)-expansion method [14] [15] , (G'/G)-expansion method [16] [17] , and so on.

Recently, the (G'/G²)-expansion method [18] [19] have received a lot of attention for number of authors who work on the method to construct exact solutions of (NPDEs) in applied sciences and engineering. The modified (G'/G²)-expansion method [20] [21] [22] has been used to find solutions for certain mathematical problems.

This work is summarized as follows. In Section 2, the description of the modified (G'/G²)-expansion method is given. In Section 3, we apply the method to Phi-4 equation [23] , Joseph-Egri (TRLW) equation [24] [25] , and Calogro-Degasperis (CD) equation [26] to obtain some of their exact traveling wave solutions. Finally, conclusions and some remarks have been presented.

2. Description of the Method

Consider a general nonlinear partial differential equation (PDE) in the form:

$P\left(u\mathrm{,}{u}_t\mathrm{,}{u}_x\mathrm{,}{u}_y\mathrm{,}{u}_{tt}\mathrm{,}{u}_{xx}\mathrm{,}{u}_{yy}\mathrm{,}{u}_{txy}\mathrm{,}\cdots \right)=\mathrm{0,}$ (1)

where $u=u\left(x,y,t,\cdots \right)$ is an unknown function, P is a polynomial in $u\left(x\mathrm{,}y\mathrm{,}t\mathrm{,}\cdots \right)$ and its partial derivatives.

• First step

Reduction of Equation (1) to a nonlinear ordinary differential equation by introducing the traveling wave variable $\xi =k\left(x+y+\cdots \right)-Vt$ , where t is a real positive number, k is the wave number and V is the velocity of travelling wave, we have:

$Q\left(U\mathrm{,}-V{U}^{\prime }\mathrm{,}k{U}^{\prime }\mathrm{,}k{U}^{\prime }\mathrm{,}{V}^2{U}^{″}\mathrm{,}{k}^2{U}^{″}\mathrm{,}{k}^2{U}^{″}\mathrm{,}-V{k}^2{U}^{‴}\mathrm{,}\cdots \right)=\mathrm{0,}$ (2)

where $U\left(\xi \right)=u\left(x\mathrm{,}y\mathrm{,}t\mathrm{,}\cdots \right)$ and Q is a polynomial of U and it's derivatives with respect to $\xi$ .

• Second step

Assume that the traveling wave solution of Equation (2) can be expressed by a polynomial in $\frac{G^{\prime }}{G^2}$ as follows:

$U\left(\xi \right)={\displaystyle \underset{i=0}{\overset{N}{\sum }}}{a}_i{\left(\frac{G^{\prime }}{G^2}\right)}^i+{\displaystyle \underset{i=1}{\overset{N}{\sum }}}{b}_i{\left(\frac{G^{\prime }}{G^2}\right)}^{-i},$ (3)

where $a_i\left(i=0,1,2,\cdots ,N\right)$ and $b_i\left(i=1,2,\cdots ,N\right)$ are constants which determined later and $G=G\left(\xi \right)$ satisfies the Riccati equation:

${\left(\frac{G^{\prime }}{G^2}\right)}^{\prime }=\rho {\left(\frac{G^{\prime }}{G^2}\right)}^2+\mu \left(\frac{G^{\prime }}{G^2}\right)+\sigma \mathrm{,}$ (4)

$\rho$ , $\mu$ and $\sigma$ are arbitrary constants.

• Third step

Calculate the value of the positive integer N by balancing the highest order nonlinear terms with the highest order derivatives appearing in Equation (2).

• Fourth step

The general solution of Equation (4), has five possible solutions as follows:

$\frac{G^{\prime }}{G^2}=\{\begin{array}{ll}-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }}{2\rho }\left(\frac{A\sinh \left(\frac{1}{2}\sqrt{\Delta }\xi \right)+B\cosh \left(\frac{1}{2}\sqrt{\Delta }\xi \right)}{A\cosh \left(\frac{1}{2}\sqrt{\Delta }\xi \right)+B\sinh \left(\frac{1}{2}\sqrt{\Delta }\xi \right)}\right),\hfill & \text{if}\mu \ne 0,\Delta \ge 0,\hfill \\ -\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }}{2\rho }\left(\frac{-A\sin \left(\frac{1}{2}\sqrt{-\Delta }\xi \right)+B\cos \left(\frac{1}{2}\sqrt{-\Delta }\xi \right)}{A\cos \left(\frac{1}{2}\sqrt{-\Delta }\xi \right)+B\sin \left(\frac{1}{2}\sqrt{-\Delta }\xi \right)}\right),\hfill & \text{if}\mu \ne 0,\Delta <0,\hfill \\ \sqrt{\frac{\sigma }{\rho }}\left(\frac{A\cos \left(\sqrt{\sigma \rho }\xi \right)+B\sin \left(\sqrt{\sigma \rho }\xi \right)}{-A\sin \left(\sqrt{\sigma \rho }\xi \right)+B\cos \left(\sqrt{\sigma \rho }\xi \right)}\right),\hfill & \text{if}\sigma \rho >0,\mu =0,\hfill \\ -\frac{\sqrt{\left|\sigma \rho \right|}}{\rho }\left(\frac{A\sinh \left(2\sqrt{\left|\sigma \rho \right|}\xi \right)+A\cosh \left(2\sqrt{\left|\sigma \rho \right|}\xi \right)+B}{A\sinh \left(2\sqrt{\left|\sigma \rho \right|}\xi \right)+A\cosh \left(2\sqrt{\left|\sigma \rho \right|}\xi \right)-B}\right),\hfill & \text{if}\sigma \rho <0,\mu =0,\hfill \\ -\frac{A}{\rho \left(A\xi +B\right)},\hfill & \text{if}\sigma =0,\mu =0,\rho \ne 0.\hfill \end{array}$ (5)

where A, B are constants and $\Delta ={\mu }^2-4\rho \sigma$ .

• Fifth step

Using Equation (4) to calculate the derivatives then substituting them and

Equation (3) into Equation (2) then collecting all terms with same power of $\frac{G^{\prime }}{G^2}$

and solving the system of algebraic equations by computational software such as Maple. Finally use Equation (5) with the results obtained to derive the exact solutions of Equation (1).

3. Applications

In this section, we will find exact traveling wave solutions to Phi-4 equation, Joseph-Egri (TRLW) equation and Calogro-Degasperis (CD) equation by applying the modified (G'/G²)-expansion method.

3.1. First Application

Consider the Phi-4 equation in the form:

$u_{tt}-u_{xx}+m^{2}u+c{u}^3=0,$ (6)

which has a significant role in the study on the nuclear and particle physics, where m and c are real valued constants. Using the transformation $u\left(x\mathrm{,}t\right)=U\left(\xi \right)$ where $\xi =kx-Vt$ , Equation (6) changed into an ODE equation of the form:

$\left(V^2-k^2\right)U^{″}+m^{2}U+c{U}^3=0.$ (7)

Balancing the highest order derivative term $U^{″}$ and nonlinear term $U^3$ in Equation (7), we obtain $N=1$ . The solution of Equation (7) will be in the following form:

$U\left(\xi \right)=a_0+a_1\left(\frac{G^{\prime }}{G^2}\right)+b_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}$ (8)

$a_0\mathrm{,}{a}_1$ and $b_1$ are constants to be determined later. The derivatives $U^{\prime }$ and $U^{″}$ compute as follows:

$U^{\prime }=a_1\rho {\left(\frac{G^{\prime }}{G^2}\right)}^2+a_1\mu \left(\frac{G^{\prime }}{G^2}\right)-b_1\rho +a_1\sigma -b_1\mu {\left(\frac{G^{\prime }}{G^2}\right)}^{-1}-b_1\sigma {\left(\frac{G^{\prime }}{G^2}\right)}^{-2}\mathrm{,}$ (9)

$\begin{array}{c}{U}^{″}=2a_1{\rho }^2{\left(\frac{G^{\prime }}{G^2}\right)}^3+3a_1\mu \rho {\left(\frac{G^{\prime }}{G^2}\right)}^2+\left(a_1{\mu }^2+2a_1\rho \sigma \right)\left(\frac{G^{\prime }}{G^2}\right)+b_1\mu \rho +a_1\mu \sigma \\ +\left(b_1{\mu }^2+2b_1\sigma \rho \right){\left(\frac{G^{\prime }}{G^2}\right)}^{-1}+3b_1\sigma \mu {\left(\frac{G^{\prime }}{G^2}\right)}^{-2}+2b_1{\sigma }^2{\left(\frac{G^{\prime }}{G^2}\right)}^{-3}\mathrm{.}\end{array}$ (10)

Now, substituting Equation (8) and Equation (10) in Equation (7) we will get a polynomial of degree three in $\left(\frac{G^{\prime }}{G^2}\right)$ , since the polynomial is equal to zero then each of its coefficients are equal to zero, as follows:

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^3:\left(2\left(V^2-k^2\right)a_1{\rho }^2+c{a}_1^3\right)=0,\\ {\left(\frac{G^{\prime }}{G^2}\right)}^2:\left(3\left(V^2-k^2\right)a_1\mu \rho +3c{a}_0{a}_1^2\right)=0,\\ {\left(\frac{G^{\prime }}{G^2}\right)}^1:\left(V^2-k^2\right)\left(a_1{\mu }^2+2a_1\rho \sigma \right)+m^2{a}_1+c\left(b_1{a}_1^2+2a_0^2{a}_1+a_1\left(a_0^2+2b_1{a}_1\right)\right)=0,\\ {\left(\frac{G^{\prime }}{G^2}\right)}^0:\left(V^2-k^2\right)\left(b_1\mu \rho +a_1\mu \sigma \right)+m^2{a}_0+c\left(4b_1{a}_0{a}_1+a_0\left(a_0^2+2b_1{a}_1\right)\right)=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}:\left(V^2-k^2\right)\left(b_1{\mu }^2+2b_1\sigma \rho \right)+m^2{b}_1+c\left(b_1\left(a_0^2+2b_1{a}_1\right)+2a_0^2{b}_1+a_1{b}_1^2\right)=0,\\ {\left(\frac{G^{\prime }}{G^2}\right)}^{-2}:3\left(V^2-k^2\right)b_1\sigma \mu +3c{b}_1^2{a}_0=0,\\ {\left(\frac{G^{\prime }}{G^2}\right)}^{-3}:2\left(V^2-k^2\right)b_1{\sigma }^2+c{b}_1^3=0.\end{array}$ (11)

Solving the above algebraic equations with the assistance of Maple, we obtain the following results:

Result 1.

$V=\pm \sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }},a_0=-\frac{m\mu }{\left(\Delta \right)cs},a_1=0,b_1=2s\sigma m.$ (12)

$V=\pm \sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }},a_0=\frac{m\mu }{\left(\Delta \right)cs},a_1=0,b_1=-2s\sigma m.$ (13)

Result 2.

$V=\pm \sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }},a_0=-\frac{m\mu }{\left(\Delta \right)cs},a_1=2s\rho m,b_1=0.$ (14)

$V=\pm \sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }},a_0=\frac{m\mu }{\left(\Delta \right)cs},a_1=-2s\rho m,b_1=0.$ (15)

where $s=\sqrt{-\frac{1}{c\Delta }}$ . Substituting (12) and (13) into the solution form (8) and using (5), we get the following solutions:

$U_{1,1}\left(\xi \right)=\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}}-\frac{m\mu }{\left(\Delta \right)cs}$

$U_{1,2}\left(\xi \right)=\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}}-\frac{m\mu }{\left(\Delta \right)cs}$

$U_{1,3}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{1,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\rho \left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}{2\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}$

$U_{2,1}\left(\xi \right)=\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}}+\frac{m\mu }{\left(\Delta \right)cs}$

$U_{2,2}\left(\xi \right)=-\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}}+\frac{m\mu }{\left(\Delta \right)cs}$

$U_{2,3}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{2,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\rho \left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}{2\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}$

where $\xi =kx-\sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }}t$ ,

$U_{3,1}\left(\xi \right)=\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}}-\frac{m\mu }{\left(\Delta \right)cs}$

$U_{3,2}\left(\xi \right)=\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}}-\frac{m\mu }{\left(\Delta \right)cs}$

$U_{3,3}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{3,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\rho \left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}{2\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}$

$U_{4,1}\left(\xi \right)=\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}}+\frac{m\mu }{\left(\Delta \right)cs}$

$U_{4,2}\left(\xi \right)=\frac{2s\sigma m}{-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}}+\frac{m\mu }{\left(\Delta \right)cs}$

$U_{4,3}\left(\xi \right)=-\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{4,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\sigma m\rho \left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}{2\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}$

where $\xi =kx+\sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }}t$ . Now, substituting (14) and (15) into the solution form (8) and using (5), we get the following solutions:

$U_{5,1}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}+2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)$

$U_{5,2}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}+2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)$

$U_{5,3}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\rho m\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{5,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}m\sqrt{-\rho \sigma }\left(A\sinh \text{}\left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}{2\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}$

$U_{6,1}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}-2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)$

$U_{6,2}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}-2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)$

$U_{6,3}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\rho m\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{6,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}m\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}{2\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}$

where $\xi =kx-\sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }}t$ ,

$U_{7,1}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}+2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)$

$U_{7,2}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}+2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)$

$U_{7,3}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\rho m\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{7,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}m\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}{2\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}$

$U_{8,1}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}-2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)$

$U_{8,2}\left(\xi \right)=\frac{m\mu }{\left(\Delta \right)cs}-2s\rho m\left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)$

$U_{8,3}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}\rho m\sqrt{\frac{\sigma }{\rho }}\left(A\cos \text{}\left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}{2\left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}$

$U_{8,4}\left(\xi \right)=\frac{\sqrt{\frac{4}{c\rho \sigma }}m\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}{2\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}$

where $\xi =kx+\sqrt{-\frac{-k^2{\mu }^2+4k^2\rho \sigma -2m^2}{\Delta }}t$ . Finally, using the fact that

$u_{i\mathrm{,}j}\left(x\mathrm{,}t\right)=U_{i\mathrm{,}j}\left(\xi \right)$ , there are two indices i and j, where i represents the number of solutions of the ODE while j represents the number of solutions of the PDE (As an example of one of the above solutions is shown in Figure 1).

3.2. Second Application

Consider the Joseph-Egri (TRLW) equation in the form:

$u_t+u_x+\alpha u{u}_x+u_{xtt}=0,$ (16)

which an important model to describe the spreading wave, where $\alpha$ is a nonzero constant. Using the transformation $u\left(x,t\right)=U\left(\xi \right)$ where $\xi =kx-Vt$ , Equation (16) changed into an ODE equation of the form:

$\left(k-V\right)U^{\prime }+\alpha kU{U}^{\prime }+k{V}^2{U}^{‴}=0.$ (17)

Figure 1. The solution $u_{\mathrm{1,3}}\left(x\mathrm{,}t\right)$ of Phi-4 equation with $A=1$ , $B=2$ , $k=1$ , $\rho =10$ , $\sigma =0.2$ , $m=1$ and $c=1$ .

Balancing the highest order derivative term $U^{‴}$ and nonlinear term $U{U}^{\prime }$ in Equation (17), we obtain $N=2$ . The solution of Equation (17) will be in the following form:

$U\left(\xi \right)=b_2{\left(\frac{G^{\prime }}{G^2}\right)}^{-2}+b_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}+a_0+a_1\left(\frac{G^{\prime }}{G^2}\right)+a_2{\left(\frac{G^{\prime }}{G^2}\right)}^2$ (18)

$a_0\mathrm{,}{a}_1\mathrm{,}{a}_2\mathrm{,}{b}_1$ and $b_2$ are constants to be determined later. Now, substituting Equation (18) and its derivatives in Equation (17) we will get a polynomial of

degree five in $\left(\frac{G^{\prime }}{G^2}\right)$ , since the polynomial is equal to zero then each of its coefficients are equal to zero, as follows:

${\left(\frac{G^{\prime }}{G^2}\right)}^5\mathrm{<mo>\; :\; </mo>}\left(24k{V}^2{a}_2{\rho }^3+2\alpha k{a}_2^2\rho \right)=\mathrm{0,}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^4\mathrm{<mo>\; :\; </mo>2}\alpha k{a}_1{a}_2\rho +\alpha k{a}_2\left(2a_2\mu +a_1\rho \right)\\ +k{V}^2\left(\left(3\left(4a_2\mu +2a_1\rho \right)\rho +18a_2\rho \mu \right)\rho +24a_2{\rho }^2\mu \right)=\mathrm{0,}\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^3\mathrm{<mo>\; :\; </mo>2}\left(k-V\right)a_2\rho +2\alpha k{a}_0{a}_2\rho +\alpha k{a}_1\left(2a_2\mu +a_1\rho \right)+\alpha k{a}_2\left(a_1\mu +2a_2\sigma \right)\\ +k{V}^2\left(\left(2\left(a_1\mu +2a_2\sigma \right)\rho +2\left(4a_2\mu +2a_1\rho \right)\mu +12a_2\rho \sigma \right)\rho \right)\\ +k{V}^2\left(\left(3\left(4a_2\mu +2a_1\rho \right)\rho +18a_2\rho \mu \right)\mu +24a_2{\rho }^2\sigma \right)=\mathrm{0,}\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^2:\left(k-V\right)\left(2a_2\mu +a_1\rho \right)+2\alpha k{b}_1{a}_2\rho +\alpha k{a}_0\left(2a_2\mu +a_1\rho \right)\\ +\alpha k{a}_1\left(a_1\mu +2a_2\sigma \right)+\alpha k{a}_2\left(-b_1\rho +a_1\sigma \right)\\ +k{V}^2\left(\left(\left(a_1\mu +2a_2\sigma \right)\mu +\left(4a_2\mu +2a_1\rho \right)\sigma \right)\rho \right)\\ +k{V}^2\left(\left(2\left(a_1\mu +2a_2\sigma \right)\rho +2\left(4a_2\mu +2a_1\rho \right)\mu +12a_2\rho \sigma \right)\mu \right)\\ +k{V}^2\left(\left(3\left(4a_2\mu +2a_1\rho \right)\rho +18a_2\rho \mu \right)\sigma \right)=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^1\mathrm{<mo>\; :\; </mo>}\left(k-V\right)\left(a_1\mu +2a_2\sigma \right)+2\alpha k{b}_2{a}_2\rho +\alpha k{b}_1\left(2a_2\mu +a_1\rho \right)\\ +\alpha k{a}_0\left(a_1\mu +2a_2\sigma \right)+\alpha k{a}_1\left(-b_1\rho +a_1\sigma \right)+\alpha k{a}_2\left(-b_1\mu -2b_2\rho \right)\\ +k{V}^2\left(\left(\left(a_1\mu +2a_2\sigma \right)\mu +\left(4a_2\mu +2a_1\rho \right)\sigma \right)\mu \right)\\ +k{V}^2\left(\left(2\left(a_1\mu +2a_2\sigma \right)\rho +2\left(4a_2\mu +2a_1\rho \right)\mu +12a_2\rho \sigma \right)\sigma \right)=\mathrm{0,}\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^0:\left(k-V\right)\left(-b_1\rho +a_1\sigma \right)+\alpha k{b}_2\left(2a_2\mu +a_1\rho \right)+\alpha k{b}_1\left(a_1\mu +2a_2\sigma \right)\\ +\alpha k{a}_0\left(-b_1\rho +a_1\sigma \right)+\alpha k{a}_1\left(-b_1\mu -2b_2\rho \right)+\alpha k{a}_2\left(-2b_2\mu -b_1\sigma \right)\\ +k{V}^2\left(\left(-\left(4b_2\mu +2b_1\sigma \right)\rho -\left(b_1\mu +2b_2\rho \right)\mu \right)\rho \right)\\ +k{V}^2\left(\left(\left(a_1\mu +2a_2\sigma \right)\mu +\left(4a_2\mu +2a_1\rho \right)\sigma \right)\sigma \right)=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}\mathrm{<mo>\; :\; </mo>}\left(k-V\right)\left(-b_1\mu -2b_2\rho \right)+\alpha k{b}_2\left(a_1\mu +2a_2\sigma \right)+\alpha k{b}_1\left(-b_1\rho +a_1\sigma \right)\\ +\alpha k{a}_0\left(-b_1\mu -2b_2\rho \right)+\alpha k{a}_1\left(-2b_2\mu -b_1\sigma \right)-2\alpha k{a}_2{b}_2\sigma \\ +k{V}^2\left(\left(-12b_2\sigma \rho -2\left(4b_2\mu +2b_1\sigma \right)\mu -2\left(b_1\mu +2b_2\rho \right)\sigma \right)\rho \right)\\ +k{V}^2\left(\left(-\left(4b_2\mu +2b_1\sigma \right)\rho -\left(b_1\mu +2b_2\rho \right)\mu \right)\mu \right)=\mathrm{0,}\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-2}\mathrm{<mo>\; :\; </mo>}\left(k-V\right)\left(-b_1\mu -2b_2\rho \right)+\alpha k{b}_2\left(a_1\mu +2a_2\sigma \right)+\alpha k{b}_1\left(-b_1\rho +a_1\sigma \right)\\ +\alpha k{a}_0\left(-b_1\mu -2b_2\rho \right)+\alpha k{a}_1\left(-2b_2\mu -b_1\sigma \right)-2\alpha k{a}_2{b}_2\sigma \\ +k{V}^2\left(\left(-12b_2\sigma \rho -2\left(4b_2\mu +2b_1\sigma \right)\mu -2\left(b_1\mu +2b_2\rho \right)\sigma \right)\rho \right)\\ +k{V}^2\left(\left(-\left(4b_2\mu +2b_1\sigma \right)\rho -\left(b_1\mu +2b_2\rho \right)\mu \right)\mu \right)=\mathrm{0,}\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-3}\mathrm{<mo>\; :\; </mo>}-2\left(k-V\right)b_2\sigma +\alpha k{b}_2\left(-b_1\mu -2b_2\rho \right)+\alpha k{b}_1\left(-2b_2\mu -b_1\sigma \right)-2\alpha k{a}_0{b}_2\sigma \\ +k{V}^2\left(-24b_2{\sigma }^2\rho +\left(-18b_2\sigma \mu -3\left(4b_2\mu +2b_1\sigma \right)\sigma \right)\mu \right)\\ +k{V}^2\left(\left(-12b_2\sigma \rho -2\left(4b_2\mu +2b_1\sigma \right)\mu -2\left(b_1\mu +2b_2\rho \right)\sigma \right)\sigma \right)=\mathrm{0,}\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-4}:\alpha k{b}_2\left(-2b_2\mu -b_1\sigma \right)-2\alpha k{b}_1{b}_2\sigma \\ +k{V}^2\left(-24b_2{\sigma }^2\mu +\left(-18b_2\sigma \mu -3\left(4b_2\mu +2b_1\sigma \right)\sigma \right)\sigma \right)=0,\end{array}$

${\left(\frac{G^{\prime }}{G^2}\right)}^{-5}:-24k{V}^2{b}_2{\sigma }^3-2\alpha k{b}_2^2\sigma =0.$

Solving the above algebraic equations with the assistance of Maple, we obtain the following results:

Result 1.

$a_0=-\frac{V^{2}k{\mu }^2+s}{\alpha k},a_1=-\frac{12V^2\mu \rho }{\alpha },a_2=-\frac{12V^2{\rho }^2}{\alpha },b_1=0,b_2=0.$ (19)

Result 2.

$a_0=-\frac{V^{2}k{\mu }^2+s}{\alpha k},a_1=0,a_2=0,b_1=-\frac{12V^2\mu \sigma }{\alpha },b_2=-\frac{12V^2{\sigma }^2}{\alpha }.$ (20)

where $s=8V^{2}k\rho \sigma -V+k$ . Substituting (19) into the solution form (18) and using (5), we get the following solutions:

$\begin{array}{c}{U}_{1,1}\left(\xi \right)=-\frac{V^{2}k{\mu }^2+s}{\alpha k}-\frac{12V^2\mu \rho \left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)}{\alpha }\\ -\frac{12V^2{\rho }^2{\left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)}^2}{\alpha }\end{array}$

$\begin{array}{c}{U}_{1,2}\left(\xi \right)=-\frac{V^{2}k{\mu }^2+s}{\alpha k}-\frac{12V^2\mu \rho \left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)}{\alpha }\\ -\frac{12V^2{\rho }^2{\left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta \rho }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-{\mu }^2+4\sigma \rho }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-{\mu }^2+4\sigma \rho }\xi }{2}\right)\right)}\right)}^2}{\alpha }\end{array}$

$U_{1,3}\left(\xi \right)=\frac{s}{\alpha k}-\frac{12V^2\rho \sigma {\left(A\cos \left(\sqrt{\sigma \rho }\xi \right)+B\sin \left(\sqrt{\sigma \rho }\xi \right)\right)}^2}{\alpha {\left(-A\sin \left(\sqrt{\sigma \rho }\xi \right)+B\cos \left(\sqrt{\sigma \rho }\xi \right)\right)}^2}$

$U_{\mathrm{1,4}}\left(\xi \right)=\frac{s}{\alpha k}+\frac{12V^2\rho \sigma {\left(A\sinh \left(2\sqrt{-\sigma \rho }\xi \right)+A\cosh \left(2\sqrt{-\sigma \rho }\xi \right)+B\right)}^2}{\alpha {\left(A\sinh \left(2\sqrt{-\sigma \rho }\xi \right)+A\cosh \left(2\sqrt{-\sigma \rho }\xi \right)-B\right)}^2}$

$U_{\mathrm{1,5}}\left(\xi \right)=\frac{k-V}{\alpha k}-\frac{12V^2{A}^2}{\alpha {\left(A\xi +B\right)}^2}$

Now, substituting (20) into the solution form (18) and using (5), we get the following solutions:

$U_{2,1}\left(\xi \right)=\frac{12V^2{\sigma }^2}{\alpha {\left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)}^2}$

$-\frac{12V^2\mu \sigma }{\alpha \left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)}-\frac{V^{2}k{\mu }^2+s}{\alpha k}$

$\begin{array}{c}{U}_{\mathrm{2,2}}\left(\xi \right)=\frac{12V^2{\sigma }^2}{\alpha {\left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)}^2}\\ -\frac{12V^2\mu \sigma }{\alpha \left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)}-\frac{V^{2}k{\mu }^2+s}{\alpha k}\end{array}$

$U_{2,3}\left(\xi \right)=\frac{12V^2\sigma \rho {\left(-A\sin \left(\sqrt{\sigma \rho }\xi \right)+B\cos \left(\sqrt{\sigma \rho }\xi \right)\right)}^2}{\alpha {\left(A\cos \left(\sqrt{\sigma \rho }\xi \right)+B\sin \left(\sqrt{\sigma \rho }\xi \right)\right)}^2}-\frac{s}{\alpha k}$

$U_{2,4}\left(\xi \right)=\frac{12V^2\sigma \rho {\left(A\sinh \left(2\sqrt{-\sigma \rho }\xi \right)+A\cosh \left(2\sqrt{-\sigma \rho }\xi \right)-B\right)}^2}{\alpha {\left(A\sinh \left(2\sqrt{-\sigma \rho }\xi \right)+A\cosh \left(2\sqrt{-\sigma \rho }\xi \right)+B\right)}^2}-\frac{s}{\alpha k}$

Finally, using the fact that $u_{i\mathrm{,}j}\left(x\mathrm{,}t\right)=U_{i\mathrm{,}j}\left(\xi \right)$ (As an example of one of the above solutions is shown in Figure 2).

Figure 2. The solution $u_{\mathrm{1,5}}\left(x\mathrm{,}t\right)$ of TRLW equation with $A=1$ , $B=2$ , $k=1$ , $V=1$ and $\alpha =-1$ .

3.3. Third Application

Consider the Calogro-Degasperis (CD) equation in the form:

$u_{xt}-4u_x{u}_{xx}-2u_y{u}_{xx}+u_{xxxy}=\mathrm{0,}$ (21)

which arises in fluid mechanics, solid state physics and plasma physics. Using the transformation $u\left(x\mathrm{,}y\mathrm{,}t\right)=U\left(\xi \right)$ where $\xi =k\left(x+y-Vt\right)$ , Equation (21) changed into an ODE equation of the form:

$-V{U^{\prime }}^{\prime }-4k{U}^{\prime }{U^{\prime }}^{\prime }-2k{U}^{\prime }{U^{\prime }}^{\prime }+k^2{{{U^{\prime }}^{\prime }}^{\prime }}^{\prime }=0.$ (22)

Balancing the highest order derivative term ${{{U^{\prime }}^{\prime }}^{\prime }}^{\prime }$ and nonlinear term $U^{\prime }{U^{\prime }}^{\prime }$ in Equation (22), we obtain $N=1$ . The solution of Equation (22) will be in the following form:

$U\left(\xi \right)=a_0+a_1\left(\frac{G^{\prime }}{G^2}\right)+b_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}$ (23)

$a_0\mathrm{,}{a}_1$ and $b_1$ are constants to be determined later. Now, substituting Equation (23) and its derivatives in Equation (22) we will get a polynomial of degree five

in $\left(\frac{G^{\prime }}{G^2}\right)$ , since the polynomial is equal to zero then each of its coefficients are equal to zero, as follows:

${\left(\frac{G^{\prime }}{G^2}\right)}^5:24a_1{\rho }^4{k}^2-12a_1^2{\rho }^{3}k=0,$

${\left(\frac{G^{\prime }}{G^2}\right)}^4:60a_1\mu {\rho }^3{k}^2-30a_1^2\mu {\rho }^{2}k=0,$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^3:-6\left(2\left(-b_1\rho +a_1\sigma \right)a_1{\rho }^2+3a_1^2{\mu }^2\rho +a_1\rho \left(a_1{\mu }^2+2a_1\rho \sigma \right)\right)k\\ -2a_1{\rho }^{2}V+(\left(2\left(a_1{\mu }^2+2a_1\rho \sigma \right)\rho +12a_1{\mu }^2\rho +12a_1{\rho }^2\sigma \right)\rho \\ {}_{{}^{}}{}^{}+36a_1{\mu }^2{\rho }^2+24a_1{\rho }^3\sigma )k^2=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^2:-6(-2b_1\mu a_1{\rho }^2+3\left(-b_1\rho +a_1\sigma \right)a_1\mu \rho +a_1\mu \left(a_1{\mu }^2+2a_1\rho \sigma \right)\\ +a_1\rho \left(b_1\mu \rho +a_1\mu \sigma \right))k+(\left(\left(a_1{\mu }^2+2a_1\rho \sigma \right)\mu +6a_1\mu \rho \sigma \right)\rho \\ +\left(2\left(a_1{\mu }^2+2a_1\rho \sigma \right)\rho +12a_1{\mu }^2\rho +12a_1{\rho }^2\sigma \right)\mu +36a_1\mu {\rho }^2\sigma )k^2\\ -3a_1\mu \rho V=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^1:-6(-2b_1\sigma a_1{\rho }^2-3b_1{\mu }^2{a}_1\rho +\left(-b_1\rho +a_1\sigma \right)\left(a_1{\mu }^2+2a_1\rho \sigma \right)\\ +a_1\mu \left(b_1\mu \rho +a_1\mu \sigma \right)+a_1\rho \left(b_1{\mu }^2+2b_1\sigma \rho \right))k+((\left(a_1{\mu }^2+2a_1\rho \sigma \right)\mu \\ +6a_1\mu \rho \sigma )\mu +\left(2\left(a_1{\mu }^2+2a_1\rho \sigma \right)\rho +12a_1{\mu }^2\rho +12a_1{\rho }^2\sigma \right)\sigma )k^2\\ -\left(a_1{\mu }^2+2a_1\rho \sigma \right)V=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^0:6(-b_1\mu \left(a_1{\mu }^2+2a_1\rho \sigma \right)+\left(-b_1\rho +a_1\sigma \right)\left(b_1\mu \rho +a_1\mu \sigma \right)\\ +a_1\mu \left(b_1{\mu }^2+2b_1\sigma \rho \right))k+(\left(6b_1\sigma \mu \rho -\left(-b_1{\mu }^2-2b_1\sigma \rho \right)\mu \right)\rho \\ +\left(\left(a_1{\mu }^2+2a_1\rho \sigma \right)\mu +6a_1\mu \rho \sigma \right)\sigma )k^2-\left(b_1\mu \rho +a_1\mu \sigma \right)V=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}:-6(-b_1\sigma \left(a_1{\mu }^2+2a_1\rho \sigma \right)-b_1\mu \left(b_1\mu \rho +a_1\mu \sigma \right)+\left(-b_1\rho +a_1\sigma \right)(b_1{\mu }^2\\ +2b_1\sigma \rho )+3a_1{\mu }^2{b}_1\sigma +2a_1\rho b_1{\sigma }^2)k+((12b_1{\sigma }^2\rho +12b_1\sigma {\mu }^2\\ -2\left(-b_1{\mu }^2-2b_1\sigma \rho \right)\sigma )\rho +\left(6b_1\sigma \mu \rho -\left(-b_1{\mu }^2-2b_1\sigma \rho \right)\mu \right)\mu )k^2\\ -\left(b_1{\mu }^2+2b_1\sigma \rho \right)V=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-2}:-6(-b_1\sigma \left(b_1\mu \rho +a_1\mu \sigma \right)-b_1\mu \left(b_1{\mu }^2+2b_1\sigma \rho \right)+3\left(-b_1\rho +a_1\sigma \right)b_1\sigma \mu \\ +2a_1\mu b_1{\sigma }^2)k+(36b_1{\sigma }^2\mu \rho +(12b_1{\sigma }^2\rho +12b_1\sigma {\mu }^2\\ -2\left(-b_1{\mu }^2-2b_1\sigma \rho \right)\sigma )\mu +\left(6b_1\sigma \mu \rho -\left(-b_1{\mu }^2-2b_1\sigma \rho \right)\mu \right)\sigma )k^2\\ -3b_1\sigma \mu V=0,\end{array}$

$\begin{array}{l}{\left(\frac{G^{\prime }}{G^2}\right)}^{-3}:-6\left(-b_1\sigma \left(b_1{\mu }^2+2b_1\sigma \rho \right)-3b_1^2{\mu }^2\sigma +2\left(-b_1\rho +a_1\sigma \right)b_1{\sigma }^2\right)k\\ +(24b_1{\sigma }^3\rho +36b_1{\sigma }^2{\mu }^2+(12b_1{\sigma }^2\rho +12b_1\sigma {\mu }^2\\ -2\left(-b_1{\mu }^2-2b_1\sigma \rho \right)\sigma )\sigma )k^2-2b_1{\sigma }^{2}V=0,\end{array}$

${\left(\frac{G^{\prime }}{G^2}\right)}^{-4}:60b_1{\sigma }^3\mu k^2+30b_1^2{\sigma }^2\mu k=0,$

${\left(\frac{G^{\prime }}{G^2}\right)}^{-5}:24b_1{\sigma }^4{k}^2+12b_1^2{\sigma }^{3}k=0.$ (24)

Solving the above algebraic equations with the assistance of Maple, we obtain the following results:

Result 1.

$V=k^2{\mu }^2-4k^2\rho \sigma ,a_0=a_0,a_1=2k\rho ,b_1=0,$ (25)

Result 2.

$V=k^2{\mu }^2-4k^2\rho \sigma ,a_0=a_0,a_1=0,b_1=-2k\sigma .$ (26)

Substituting (25) into the solution form (23) and using (5), we get the following solutions:

$U_{\mathrm{1,1}}\left(\xi \right)=2k\rho \left(-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}\right)+a_0$

$U_{\mathrm{1,2}}\left(\xi \right)=2k\rho \left(-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}\right)+a_0$

$U_{1,3}\left(\xi \right)=\frac{2k\rho \sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}{-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)}+a_0,\sigma ,\rho >0$

$U_{\mathrm{1,4}}\left(\xi \right)=\frac{2k\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}{A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B}+a_0$

$U_{\mathrm{1,5}}\left(\xi \right)=\frac{2kA}{A\xi +B}+a_0$

Now, substituting (26) into the solution form (23) and using (5), we get the following solutions:

$U_{\mathrm{2,1}}\left(\xi \right)=\frac{2k\sigma }{-\frac{\mu }{2\rho }-\frac{\sqrt{\Delta }\left(A\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+B\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}{2\rho \left(B\cosh \left(\frac{\sqrt{\Delta }\xi }{2}\right)+A\sinh \left(\frac{\sqrt{\Delta }\xi }{2}\right)\right)}}+a_0$

$U_{\mathrm{2,2}}\left(\xi \right)=-\frac{2k\sigma }{-\frac{\mu }{2\rho }-\frac{\sqrt{-\Delta }\left(A\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)-B\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}{2\rho \left(A\sin \left(\frac{\sqrt{-\Delta }\xi }{2}\right)+B\cos \left(\frac{\sqrt{-\Delta }\xi }{2}\right)\right)}}+a_0$

$U_{2,3}\left(\xi \right)=-\frac{2k\sigma \left(-A\sin \left(\sqrt{\rho \sigma }\xi \right)+B\cos \left(\sqrt{\rho \sigma }\xi \right)\right)}{\sqrt{\frac{\sigma }{\rho }}\left(A\cos \left(\sqrt{\rho \sigma }\xi \right)+B\sin \left(\sqrt{\rho \sigma }\xi \right)\right)}+a_0,\sigma ,\rho >0$

$U_{\mathrm{2,4}}\left(\xi \right)=\frac{2k\sigma \rho \left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)-B\right)}{\sqrt{-\rho \sigma }\left(A\sinh \left(2\sqrt{-\rho \sigma }\xi \right)+A\cosh \left(2\sqrt{-\rho \sigma }\xi \right)+B\right)}+a_0$

Finally, using the fact that $u_{i\mathrm{,}j}\left(x\mathrm{,}y\mathrm{,}t\right)=U_{i\mathrm{,}j}\left(\xi \right)$ , where $\xi =k\left(x+y-\left(k^2{\mu }^2-4k^2\right)t\right)$ (As an example of one of the above solutions is shown in Figure 3).

Figure 3. The solution $u_{\mathrm{1,5}}\left(x\mathrm{,}y\right)$ of CD equation with $A=1$ , $B=1$ , $k=1$ , and $a_0=1$ .

4. Conclusion and Future Work

A modified (G'/G²)-expansion method has been proposed in this paper to obtain exact traveling wave solutions to the Phi-4 equation, the Joseph-Egri (TRLW) equation, and the Calogro-Degasperis (CD) equation. All solutions obtained to these equations have been verified by substitution into the original equation. With the help of the symbolic computation software Maple, this method is easy to implement and powerful. For future work, we will try to modify the auxiliary equation in order to get different new solutions.

Conflicts of Interest

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

References

[1]	Hereman, W. and Malfliet, W. (2005) The Tanh Method: A Tool to Solve Nonlinear Partial Differential Equations with Symbolic Software. Proceedings 9th World Multi-Conference on Systemics, Cybernetics and Informatics, Orlando, 10-13 July 2005, 165-168.
[2]	Zahran, E.H.M. and Khater, M.M.A. (2016) Modified Extended Tanh-Function Method and Its Applications to the Bogoyavlenskii Equation. Applied Mathematical Modelling, 40, 1769-1775. https://doi.org/10.1016/j.apm.2015.08.018
[3]	Jawad, A.J.M. (2012) The Sine-Cosine Function Method for the Exact Solutions of Nonlinear Partial Differential Equations. System, 2, 3.
[4]	Wazwaz, A.-M. (2004) A Sine-Cosine Method for Handling Nonlinear Wave Equations. Mathematical and Computer Modelling, 40, 499-508. https://doi.org/10.1016/j.mcm.2003.12.010
[5]	Hietarinta, J. (2005) Hirota’s Bilinear Method and Soliton Solutions. Physics AUC, 15, 31-37.
[6]	Wazwaz, A.-M. (2007) Multiple-Soliton Solutions for the KP Equation by Hirota’s Bilinear Method and by the Tanh-Coth Method. Applied Mathematics and Computation, 190, 633-640. https://doi.org/10.1016/j.amc.2007.01.056
[7]	Wang, M.L., Zhou, Y.B. and Li, Z.B. (1996) Application of a Homogeneous Balance Method to Exact Solutions of Nonlinear Equations in Mathematical Physics. Physics Letters A, 216, 67-75. https://doi.org/10.1016/0375-9601(96)00283-6
[8]	Abdel Rady, A.S., Osman, E.S. and Khalfallah, M. (2010) The Homogeneous Balance Method and Its Application to the Benjamin-Bona-Mahoney (BBM) Equation. Applied Mathematics and Computation, 217, 1385-1390. https://doi.org/10.1016/j.amc.2009.05.027
[9]	Yusufoglu, E. (2008) New Solitonary Solutions for the MBBM Equations Using Expfunction Method. Physics Letters A, 372, 442-446. https://doi.org/10.1016/j.physleta.2007.07.062
[10]	He, J.-H. and Wu, X.-H. (2006) Exp-Function Method for Nonlinear Wave Equations. Chaos, Solitons & Fractals, 30, 700-708. https://doi.org/10.1016/j.chaos.2006.03.020
[11]	Akinyemi, L., Senol, M. and Iyiola, O.S. (2021) Exact Solutions of the Generalized Multidimensional Mathematical Physics Models via Sub-Equation Method. Mathematics and Computers in Simulation, 182, 211-233. https://doi.org/10.1016/j.matcom.2020.10.017
[12]	Durur, H., Kurt, A. and Tasbozan, O. (2020) New Travelling Wave Solutions for KdV6 Equation Using Sub Equation Method. Applied Mathematics and Nonlinear Sciences, 5, 455-460. https://doi.org/10.2478/amns.2020.1.00043
[13]	Demiray, S., Unsal, O. and Bekir, A. (2014) New Exact Solutions for Boussinesq Type Equations by Using (G’/G, 1/G) and (1/G’)-Expansion Methods. Acta Physica Polonica A, 125, 1093-1098. https://doi.org/10.12693/APhysPolA.125.1093
[14]	Mamun Miah, M., Shahadat Ali, H.M., Ali Akbar, M. and Wazwaz, A.M. (2017) Some Applications of the (G’/G, 1/G)-Expansion Method to Find New Exact Solutions of NLEEs. The European Physical Journal Plus, 132, 1-15. https://doi.org/10.1140/epjp/i2017-11571-0
[15]	Li, L.-X., Li, E.-Q. and Wang, M.-L. (2010) The (G’/G, 1/G)-Expansion Method and Its Application to Travelling Wave Solutions of the Zakharov Equations. Applied Mathematics—A Journal of Chinese Universities, 25, 454-462. https://doi.org/10.1007/s11766-010-2128-x
[16]	Bekir, A. (2008) Application of the (G’/G)-Expansion Method for Nonlinear Evolution Equations. Physics Letters A, 372, 3400-3406. https://doi.org/10.1016/j.physleta.2008.01.057
[17]	Wang, M.L., Li, X.Z. and Zhang, J.L. (2008) The (G’/G)-Expansion Method and Travelling Wave Solutions of Nonlinear Evolution Equations in Mathematical Physics. Physics Letters A, 372, 417-423. https://doi.org/10.1016/j.physleta.2007.07.051
[18]	Bibi, S., Mohyud-Din, S.T., Ullah, R., Ahmed, N. and Khan, U. (2017) Exact Solutions for STO and (3+1)-Dimensional KdV-ZK Equations Using (G’/G2)-Expansion Method. Results in Physics, 7, 4434-4439. https://doi.org/10.1016/j.rinp.2017.11.009
[19]	Mohyud-Din, S.T. and Bibi, S. (2018) Exact Solutions for Nonlinear Fractional Differential Equations Using (G’/G2)-Expansion Method. Alexandria Engineering Journal, 57, 1003-1008. https://doi.org/10.1016/j.aej.2017.01.035
[20]	Aljahdaly, N.H. (2019) Some Applications of the Modified (G’/G2)-Expansion Method in Mathematical Physics. Results in Physics, 13, Article ID: 102272. https://doi.org/10.1016/j.rinp.2019.102272
[21]	Behera, S., Aljahdaly, N.H. and Virdi, J.P.S. (2022) On the Modified (G’/G2)-Expansion Method for Finding Some Analytical Solutions of the Traveling Waves. Journal of Ocean Engineering and Science, 7, 313-320. https://doi.org/10.1016/j.joes.2021.08.013
[22]	Zayed, E.M. and Arnous, A.H. (2014) The Modified (w/g)-Expansion Method and Its Applications for Solving the Modified Generalized Vakhnenko Equation. Italian Journal of Pure and Applied Mathematics, 32, 477-492.
[23]	Akter, J. and Akbar, M.A. (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
[24]	Taghizadeha, N., Mirzazadeha, M. and Paghaleh, A.S. (2012) Exact Travelling Wave Solutions of Joseph-Egri(TRLW) Equation by the Extended Homogeneous Balance Method. International Journal of Applied Mathematics and Computation, 4, 96-104.
[25]	Fan, F.T., Zhou, Y.Q. and Liu, Q. (2019) Bifurcation of Traveling Wave Solutions for the Joseph-Egri Equation. Reports on Mathematical Physics, 83, 175-190. https://doi.org/10.1016/S0034-4877(19)30038-2
[26]	Mohyud-Din, S.T., Nawaz, T., Azhar, E. and Akbar, M.A. (2017) Fractional Sub-Equation Method to Space-Time Fractional Calogero-Degasperis and Potential Kadomtsev-Petviashvili Equations. Journal of Taibah University for Science, 11, 258-263. https://doi.org/10.1016/j.jtusci.2014.11.010