Investigating One-, Two-, and Triple-Wave Solutions via Multiple Exp-Function Method Arising in Engineering Sciences

Abstract

The multiple Exp-function method is employed for searching the multiple soliton solutions for the new extended ($3+1$)-dimensional Jimbo-Miwa-like (JM) equation, the extended ($2+1$)-dimensional Calogero-Bogoyavlenskii-Schiff (eCBS) equation, the generalization of the ($2+1$)-dimensional Bogoyavlensky-Konopelchenko (BK) equation, and a variable-coefficient extension of the DJKM (vDJKM) equation, which contain one-soliton-, two-soliton-, and triple-soliton-kind solutions. The physical phenomena of these gained multiple soliton solutions are analyzed and indicated in figures by selecting suitable values.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Many nonlinear physical phenomena arise in various fields of engineering and science such as fluid dynamics, nuclear reactor dynamics, hydrodynamics, molecular biology, quantum mechanics, plasma physics, optical fibres, and solid state physics. To describe these complex physical phenomena, nonlinear differential equations play a significant role. Therefore, obtaining the solutions of these nonlinear equations is a topic of great interest in the study of many fields of science. To better understand the working of a physical problem, a mathematical model was brought into the picture in the form of nonlinear PDEs. The solutions of partial differential equations give the detailed summary about the nature of phenomena involved. Many numerical and analytical methods have been derived to deal with these kinds of scientific problems. We need to adopt an effective and powerful method to investigate such type of mathematical model which gives the solutions upholding physical reality. In most of the analytic techniques, linearization of the system is the main topic to focus on, and also, it is assumed that the nonlinearities are relatively insignificant. Sometimes, these assumptions make a strong effect on the solutions in respect to the real physics of the phenomena involved. Thus, finding the solutions of nonlinear ODEs and PDEs is still a significant problem. For this, we need new techniques to develop approximate and exact solutions. Several analytical and numerical techniques have been formulated for tackling these types of nonlinear models, including the Exp-function method [1, 2], the homotopy analysis method [3], the homotopy perturbation method [4], the ($G^{\prime }/G$)-expansion method [5], the improved $\tan \left(\varphi /2\right)$-expansion method ([6–8]), the Hirota bilinear method ([9–13]), the He variational principle [14, 15], the binary Darboux transformation [16], the Lie group analysis [17, 18], and the Bäcklund transformation method [19]. Moreover, many powerful methods have been used to investigate the new properties of mathematical models which are symbolizing serious real world problems ([20–24]). This study is aimed at investigating the following four forms of nonlinear PDEs:

(i)

Extended ($3+1$)-dimensional Jimbo-Miwa-like (JM) equations ([25–34]) is given as $$\begin{array}{}\text{(1)}& {\Psi }_{xxxy}+3{\Psi }_{xx}{\Psi }_y+3{\Psi }_x{\Psi }_{xy}+2{\Psi }_{yt}-3\left({\Psi }_{xz}+{\Psi }_{yz}+{\Psi }_{zz}\right)=0,\end{array}$$and $$\begin{array}{}\text{(2)}& {\Psi }_{xxxy}+3{\Psi }_{xx}{\Psi }_y+3{\Psi }_x{\Psi }_{xy}-3{\Psi }_{xz}+2\left({\Psi }_{xt}+{\Psi }_{yt}+{\Psi }_{zt}\right)=0,\end{array}$$explains some ($3+1$)-dimensional waves in physics [25, 26], in which $\Psi \left(x,y,z,t\right)$ is a function including spaces $x,y,z$ and time $t$. In addition, this equation cannot pass certain conventional integrability tests [27]. Qi et al. [32] studied the extended ($3+1$)-dimensional Jimbo-Miwa-like equation via the Hirota method and obtained the solitary-wave solutions and novel exact soliton solutions.

(ii)

The extended ($2+1$)-dimensional Calogero-Bogoyavlenskii-Schiff (eCBS) equation [35] is given as $$\begin{array}{}\text{(3)}& {\Psi }_{xt}+{\Psi }_{xxxy}+6{\Psi }_{xx}{\Psi }_y+6{\Psi }_x{\Psi }_{xy}+{\delta }_1{\Psi }_{xy}+{\delta }_2{\Psi }_{xx}+{\delta }_3\left({\Psi }_{xxxx}+12{\Psi }_x{\Psi }_{xx}\right)=0,\end{array}$$where ${\delta }_1,{\delta }_2,$ and ${\delta }_3$ are arbitrary values. The eCBS equation will be changed to a ($2+1$)-dimensional CBS equation when ${\delta }_1={\delta }_2={\delta }_3=0$ [36]. Ren et al. [35] investigated the extended ($2+1$)-dimensional Calogero-Bogoyavlenskii-Schiff-like equation by using the generalized bilinear operators based on a prime number $p=3$.

(iii)

The generalization of the ($2+1$)-dimensional Bogoyavlensky-Konopelchenko (BK) equation [37] is given as $$\begin{array}{}\text{(4)}& {\Psi }_{xt}+\alpha \left({\Psi }_{xxxx}+6{\Psi }_x{\Psi }_{xx}\right)+\beta \left({\Psi }_{xxxy}+3{\Psi }_{xx}{\Psi }_y+3{\Psi }_x{\Psi }_{xy}\right)+{\delta }_1{\Psi }_{xx}+{\delta }_2{\Psi }_{xy}+{\delta }_3{\Psi }_{yy}=0.\end{array}$$

 

The ($3+1$)-dimensional VC B-type KP equation is extended from the KP equation and can explain some notable ($3+1$)-dimensional waves in fluid dynamics [37]. Also, Zhang and Pang [38] found the lump and lump-type solutions through the Hirota bilinear form for the variable-coefficient B-type Kadomtsev-Petviashvili equation.

(iv)

The variable-coefficient extension of the DJKM (vDJKM) equation ([39–48]) is given as $$\begin{array}{}\text{(5)}& {\Psi }_{xxxxy}+4{\Psi }_x{\Psi }_{xxy}+2{\Psi }_y{\Psi }_{xxx}+6{\Psi }_{xx}{\Psi }_{xy}-\alpha {\Psi }_{yyy}-2\beta g\left(t\right){\Psi }_{xxt}+h\left(t\right){\Psi }_{xxy}=0,\end{array}$$where the wave amplitude $u\left(x,y,t\right)$ is a function of variables and $x,y,t$, $g\left(t\right)$, and $h\left(t\right)$ are functions of $t$. The vD JKM equation (6) reduces to the DJKM equation when $g\left(t\right)=1$ and $h\left(t\right)=0$ as follows: $$\begin{array}{}\text{(6)}& {\Psi }_{xxxxy}+4{\Psi }_x{\Psi }_{xxy}+2{\Psi }_y{\Psi }_{xxx}+6{\Psi }_{xx}{\Psi }_{xy}-\alpha {\Psi }_{yyy}-2\beta {\Psi }_{xxt}=0.\end{array}$$

The soliton solutions to a few ($3+1$)-dimensional generalized nonlinear integrable equations have been constructed. The adopted approach is the multiple Exp-function method, which was presented earlier in the valuable work by Ma et al. [49]. Recently, special kinds of reductions of soliton solutions to rational functions that are being actively studied are lump solutions to nonlinear partial differential equations by Ma and Zhou [50] and their interactions with solitons to the Hirota-Satsuma-Ito equation in ($2+1$)-dimensions by Ma [51], even for linear PDEs by the same author [52]. Cordero et al. constructed the stability analysis of the fourth-order iterative method for finding multiple roots of nonlinear equations [53]. Gao et al. obtained the optical soliton solutions of the cubic-quartic nonlinear Schrödinger and resonant nonlinear Schrödinger equations with the parabolic law [54]. The truncated Painlevé expansion was employed to derive a Bäcklund transformation of a ($2+1$)-dimensional nonlinear system by Zhao and Han [55]. Also, Gao et al. studied on the conformable ($2+1$)-dimensional Ablowitz-KaupNewell-Segur equation in order to show the existence of complex combined dark-bright soliton solutions [56]. The same authors presented the nonlinear Zoomeron equation by using the newly extended direct algebraic technique [57]. The instability modulation for the ($2+1$)-dimension paraxial wave equation and its new optical soliton solutions in Kerr’s media by utilizing the modified auxiliary expansion method have been studied in [58]. In [59], the authors found new complex solitons to the perturbed nonlinear Schrödinger model with the help of an analytical method. In this paper, we will study the multiple Exp-function method for determining the multiple soliton solutions. The multiple Exp-function method used by some of the powerful authors for various nonlinear equations include the nonlinear evolution equations [60], the ($2+1$)-dimensional Calogero-Bogoyavlenskii-Schiff equation [61], the generalized ($1+1$)-dimensional and ($2+1$)-dimensional Ito equations [62], a new generalization of the associated Camassa-Holm equation [63], the ($3+1$)-dimensional generalized KP and BKP equations [64], and the new ($2+1$)-dimensional Korteweg-de Vries equation [65]. In [65], Liu et al. utilized the multiple Exp-function method for the most well-known equation, namely, the Korteweg-de Vries (KdV) equation, and gained one-soliton-, two-soliton-, and three-soliton-type solutions with interpretations for the obtained soliton solutions. Moreover, they analyzed and illustrated the propagation and interaction of some soliton solutions by selecting appropriate values. Also, the elasticity and being unchanged for some solutions when interacting among three solitons and after the collision were investigated.

The rest of this paper is structured as follows: the multiple Exp-function scheme is summarized in Section 2. In Sections 3–7, the extended JM equations, the extended CBS equation, the generalized BK equation, and the vDJKM equation, respectively, will be investigated to find one-soliton, two-soliton, and triple-soliton solutions. In the last section, the conclusions are given.

2. Multiple Exp-Function Method

This section elucidates a systematic explanation of the multiple Exp-function method [60–64] so that it can be further applied to the nonlinear PDEs in order to furnish its exact solutions:

Step 1.

The following NLPDE are as follows: $$\begin{array}{}\text{(7)}& \mathcal{N}\left(x,y,t,\Psi ,{\Psi }_x,{\Psi }_y,{\Psi }_z,{\Psi }_t,{\Psi }_{xx},{\Psi }_{tt},\cdots \right)=0.\end{array}$$

We commence a sequence of novel variables ${\xi }_i={\xi }_i\left(x,y,z,t\right)$, $1\le i\le n$, by solvable PDEs, for example, the linear ones: $$\begin{array}{c}\text{(8)}& {\xi }_{i,x}={\alpha }_i{\xi }_i,{\xi }_{i,y}={\beta }_i{\xi }_i,\\ {\xi }_{i,z}={\gamma }_i{\xi }_i,\\ {\xi }_{i,t}={\delta }_i{\xi }_i,\\ \hspace{1em}1\le i\le n,\end{array}$$where ${\alpha }_i,{\beta }_i,{\gamma }_i,1\le i\le n$ are the angular wave numbers and ${\delta }_i,1\le i\le n$ are the wave frequencies. It must be pointed out that this is frequently the initiating step for constructing the exact solutions to nonlinear partial differential equations; moreover, solving such linear equations redounds to the exponential function solutions, $$\begin{array}{}\text{(9)}& {\xi }_i={\varpi }_i{e}^{{\theta }_i},\hspace{1em}{\theta }_i={\alpha }_{i}x+{\beta }_{i}y+{\gamma }_{i}z-{\delta }_{i}t,1\le i\le n,\end{array}$$in which ${\varpi }_i,1\le i\le n$, are undetermined values.

Step 2.

Supposing the solution of equation (7) happens to be of the following form in terms of new variables ${\xi }_i,1\le i\le n$: $$\begin{array}{}\text{(10)}& \Psi \left(x,y,z,t\right)=\frac{\Delta \left({\xi }_1,{\xi }_2,\cdots ,{\xi }_n\right)}{\Omega \left({\xi }_1,{\xi }_2,\cdots ,{\xi }_n\right)},\hspace{1em}\Delta =\sum _{r,s=1}^n\sum _{i,j=0}^M{\Delta }_{rs,ij}{\xi }_r^i{\xi }_s^j,\Omega =\sum _{r,s=1}^n\sum _{i,j=0}^N{\Omega }_{rs,ij}{\xi }_r^i{\xi }_s^j,\end{array}$$in which ${\Delta }_{rs,ij}$ and ${\Omega }_{rs,ij}$ are values to be settled. Plugging equation (10) into equation (7) and ordering the numerator of the rational function to zero, we can gain a series of the nonlinear algebraic equations about the variables ${\alpha }_i,{\beta }_i,{\gamma }_i,{\delta }_i,{\Delta }_{rs,ij}$ and ${\Omega }_{rs,ij}$. Solving the solutions for these nonlinear algebraic equations and inserting these solutions into equation (10), the multiple soliton solutions for equation (7) can be achieved as follows: $$\begin{array}{}\text{(11)}& \Psi \left(x,y,z,t\right)=\frac{\Delta \left({\varpi }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}x-{\delta }_{1}t},\cdots ,{\varpi }_n{e}^{{\alpha }_{n}x+{\beta }_{n}y+{\gamma }_{n}x-{\delta }_{n}t}\right)}{\Omega \left({\varpi }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}x-{\delta }_{1}t},\cdots ,{\varpi }_n{e}^{{\alpha }_{n}x+{\beta }_{n}y+{\gamma }_{n}x-{\delta }_{n}t}\right)},\end{array}$$and also we have $$\begin{array}{c}\text{(12)}& {\Delta }_t=\sum _{i=1}^n{\Delta }_{{\xi }_i}{\xi }_{i,t},{\Omega }_t=\sum _{i=1}^n{\Omega }_{{\xi }_i}{\xi }_{i,t},\\ {\Delta }_x=\sum _{i=1}^n{\Delta }_{{\xi }_i}{\xi }_{i,x}\\ {\Omega }_x=\sum _{i=1}^n{\Omega }_{{\xi }_i}{\xi }_{i,x}\\ {\Delta }_y=\sum _{i=1}^n{\Delta }_{{\xi }_i}{\xi }_{i,y},\\ {\Omega }_y=\sum _{i=1}^n{\Omega }_{{\xi }_i}{\xi }_{i,y},\\ \Delta =\Delta \left({\xi }_1,{\xi }_2,\cdots ,{\xi }_n\right),\\ \Omega =\Omega \left({\xi }_1,{\xi }_2,\cdots ,{\xi }_n\right),\\ {\Delta }_{{\xi }_i}=\frac{d\Delta }{d_{{\xi }_i}},\\ {\Omega }_{{\xi }_i}=\frac{d\Omega }{d_{{\xi }_i}},\\ {\xi }_{i,t}=\frac{d{\xi }_i}{dt},\\ {\xi }_{i,x}=\frac{d{\xi }_i}{dx},\\ {\xi }_{i,y}=\frac{d{\xi }_i}{dy},\\ {\Psi }_t=\frac{\Omega {\sum }_{i=1}^n{\Delta }_{{\xi }_i}{\xi }_{i,t}-\Delta {\sum }_{i=1}^n{\Omega }_{{\xi }_i}{\xi }_{i,t}}{{\Omega }^2},\\ {\Psi }_x=\frac{\Omega {\sum }_{i=1}^n{\Delta }_{{\xi }_i}{\xi }_{i,x}-\Delta {\sum }_{i=1}^n{\Omega }_{{\xi }_i}{\xi }_{i,x}}{{\Omega }^2},\\ {\Psi }_y=\frac{\Omega {\sum }_{i=1}^n{\Delta }_{{\xi }_i}{\xi }_{i,y}-\Delta {\sum }_{i=1}^n{\Omega }_{{\xi }_i}{\xi }_{i,y}}{{\Omega }^2},\\ {\Psi }_z=\frac{\Omega {\sum }_{i=1}^n{\Delta }_{{\xi }_i}{\xi }_{i,z}-\Delta {\sum }_{i=1}^n{\Omega }_{{\xi }_i}{\xi }_{i,z}}{{\Omega }^2}.\end{array}$$

3. Multiple Soliton Solutions for the Extended ($3+1$) JM Equation (1)

3.1. Set I: One-Wave Solution

We commence with a one-wave function based on the statement in Step 2 in the previous section; we suppose that equation (1) has the rational function of the one-wave solution as shown in the following form: $$\begin{array}{}\text{(13)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_1}{{\Omega }_1},\hspace{1em}{\Omega }_1=1+{\rho }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t},{\Delta }_1={\sigma }_1+{\rho }_1{\sigma }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t},\end{array}$$in which ${\sigma }_1$ and ${\sigma }_2$ are unfound constants. Plugging (13) into equation (1), we gain to the following case: $$\begin{array}{c}\text{(14)}& {\rho }_1={\rho }_1,{\sigma }_1={\sigma }_1,\\ {\sigma }_2=\frac{{\rho }_2\left(2{\alpha }_1{\rho }_1+{\sigma }_1\right)}{{\rho }_1},\\ {\delta }_1=\frac{{\alpha }_1^3{\beta }_1-3{\alpha }_1{\gamma }_1-3{\beta }_1{\gamma }_1-3{\gamma }_1^2}{2{\beta }_1}.\end{array}$$

Therefore, the resulting one-wave solution reads as $$\begin{array}{}\text{(15)}& \Psi \left(x,y,z,t\right)=\frac{{\sigma }_1+{\rho }_2\left(2{\alpha }_1{\rho }_1+{\sigma }_1\right)e^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-\left(\left({\alpha }_1^3{\beta }_1-3{\alpha }_1{\gamma }_1-3{\beta }_1{\gamma }_1-3{\gamma }_1^2\right)/\left(2{\beta }_1\right)\right)t}}{1+{\rho }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-\left(\left({\alpha }_1^3{\beta }_1-3{\alpha }_1{\gamma }_1-3{\beta }_1{\gamma }_1-3{\gamma }_1^2\right)/\left(2{\beta }_1\right)\right)t}}.\end{array}$$

By selecting the suitable values of parameters, the graphic presentation of the periodic wave solution is presented in Figure 1 including the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-1$, $x=0$, and $x=1$.

[figure omitted; refer to PDF]

3.2. Set II: Two-Wave Solutions

We commence with the two-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (1) has the rational function of two-wave solutions as shown in the following form: $$\begin{array}{c}\text{(16)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_2}{{\Omega }_2},\text{(17)}& {\Omega }_2=1+{\sigma }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t}+{\sigma }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y+{\gamma }_{2}z-{\delta }_{2}t}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y+\left({\gamma }_1+{\gamma }_2\right)z-\left({\delta }_1+{\delta }_2\right)t},{\Delta }_2={\rho }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t}+{\rho }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y+{\gamma }_{2}z-{\delta }_{2}t}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y+\left({\gamma }_1+{\gamma }_2\right)z-\left({\delta }_1+{\delta }_2\right)t}.\end{array}$$

Plugging (16) along with (17) into equation (1), we get to the following case: $$\begin{array}{c}\text{(18)}& {\sigma }_1={\sigma }_1,{\sigma }_2={\sigma }_2,\\ {\delta }_1=\frac{1}{2}{\alpha }_1^3,\\ {\delta }_2=\frac{1}{2}{\alpha }_2^3,\\ {\rho }_1=2{\alpha }_1{\sigma }_1,\\ {\rho }_2=2{\alpha }_2{\sigma }_2,\\ {\rho }_{12}=\frac{\left({\beta }_1-{\beta }_2\right)\left({\alpha }_1-{\alpha }_2\right)}{2{\alpha }_1{\alpha }_2\left({\beta }_1+{\beta }_2\right)},\\ {\sigma }_{12}=\frac{\left({\beta }_1-{\beta }_2\right)\left({\alpha }_1-{\alpha }_2\right)}{\left({\beta }_1+{\beta }_2\right)\left({\alpha }_1+{\alpha }_2\right)}.\end{array}$$

Therefore, the resulting two-wave solution reads as $$\begin{array}{}\text{(19)}& {\Psi }_2\left(x,y,z,t\right)=\frac{2{\alpha }_1{\sigma }_1{e}^{{\Lambda }_1}+2{\alpha }_2{\sigma }_2{e}^{{\Lambda }_2}+2{\sigma }_1{\sigma }_2\left(\left(\left({\beta }_1-{\beta }_2\right)\left({\alpha }_1-{\alpha }_2\right)\right)/\left({\beta }_1+{\beta }_2\right)\right)e^{{\Lambda }_1+{\Lambda }_2}}{1+{\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_1{\sigma }_2\left(\left(\left({\beta }_1-{\beta }_2\right)\left({\alpha }_1-{\alpha }_2\right)\right)/\left(\left({\beta }_1+{\beta }_2\right)\left({\alpha }_1+{\alpha }_2\right)\right)\right)e^{{\Lambda }_1+{\Lambda }_2}},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y+{\gamma }_{i}z-\left(1/2\right){\alpha }_i^{3}t,i=1,2$. By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 2 including the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

3.3. Set III: Triple-Wave Solutions

We commence with three-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (1) has the rational function of triple-wave solutions as shown in the following form: $$\begin{array}{c}\text{(20)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\text{(21)}& {\Omega }_3=1+{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+{\rho }_3{e}^{{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\rho }_1{\rho }_3{\rho }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\rho }_2{\rho }_3{\rho }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_3{\rho }_{12}{\rho }_{13}{\rho }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},{\Delta }_3={\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_3{e}^{{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\sigma }_1{\sigma }_3{\sigma }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\sigma }_2{\sigma }_3{\sigma }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_3{\sigma }_{12}{\sigma }_{13}{\sigma }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y+{\gamma }_{i}x-{\delta }_{i}t,i=1,2,3$. Plugging (20) along with (21) into equation (1), we obtain the following case: $$\begin{array}{c}\text{(22)}& {\sigma }_i={\sigma }_i,{\delta }_i=\frac{1}{2}{\alpha }_i^3,\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,3,\\ {\rho }_{12}=\frac{\left({\beta }_1-{\beta }_2\right)\left({\alpha }_1-{\alpha }_2\right)}{2{\alpha }_1{\alpha }_2\left({\beta }_1+{\beta }_2\right)},\\ {\sigma }_{12}=\frac{\left({\beta }_1-{\beta }_2\right)\left({\alpha }_1-{\alpha }_2\right)}{\left({\beta }_1+{\beta }_2\right)\left({\alpha }_1+{\alpha }_2\right)},\\ {\rho }_{13}=\frac{\left({\beta }_1-{\beta }_3\right)\left({\alpha }_1-{\alpha }_3\right)}{2{\alpha }_1{\alpha }_3\left({\beta }_1+{\beta }_3\right)},\\ {\sigma }_{13}=\frac{\left({\beta }_1-{\beta }_3\right)\left({\alpha }_1-{\alpha }_3\right)}{\left({\beta }_1+{\beta }_3\right)\left({\alpha }_1+{\alpha }_3\right)},\\ {\rho }_{23}=\frac{\left({\beta }_2-{\beta }_3\right)\left({\alpha }_2-{\alpha }_3\right)}{2{\alpha }_2{\alpha }_3\left({\beta }_2+{\beta }_3\right)},\\ {\sigma }_{23}=\frac{\left({\beta }_2-{\beta }_3\right)\left({\alpha }_2-{\alpha }_3\right)}{\left({\beta }_2+{\beta }_3\right)\left({\alpha }_2+{\alpha }_3\right)}.\end{array}$$

Therefore, the resulting three-wave solution reads as $$\begin{array}{}\text{(23)}& {\Psi }_3\left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\end{array}$$in which ${\Delta }_3$ and ${\Omega }_3$ with their relations are given in (21) and (22). By selecting suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 3 including the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

4. Multiple Soliton Solutions for the Extended ($3+1$) JM Equation (2)

4.1. Set I: One-Wave Solution

We commence with the one-wave function based on the statement in Step 2 in the previous section; we suppose that equation (2) has the rational function of the one-wave solution as shown in the following form: $$\begin{array}{}\text{(24)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_1}{{\Omega }_1},\hspace{1em}{\Omega }_1={\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t},{\Delta }_1={\sigma }_1+{\sigma }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t},\end{array}$$in which ${\rho }_1,{\rho }_2,{\sigma }_1,$ and ${\sigma }_2$ are unfound constants. Plugging (24) into equation (1), we gain the following case: $$\begin{array}{c}\text{(25)}& {\rho }_1={\rho }_1,{\rho }_2={\rho }_2,\\ {\sigma }_1={\sigma }_1,\\ {\sigma }_2=\frac{{\rho }_2\left(2{\alpha }_1{\rho }_1+{\sigma }_1\right)}{{\rho }_1},\\ {\delta }_1=\frac{{\alpha }_1\left({{\alpha }_1}^2{\beta }_1-3{\gamma }_1\right)}{2{\alpha }_1+2{\beta }_1+2{\gamma }_1}.\end{array}$$

Therefore, the resulting one-wave solution reads as $$\begin{array}{}\text{(26)}& \Psi \left(x,y,z,t\right)=\frac{{\sigma }_1+{\rho }_2\left(2{\alpha }_1{\rho }_1+{\sigma }_1\right)e^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-\left(\left({\alpha }_1\left({\alpha }_1^2{\beta }_1-3{\gamma }_1\right)\right)/\left(2{\alpha }_1+2{\beta }_1+2{\gamma }_1\right)\right)t}}{{\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-\left(\left({\alpha }_1\left({\alpha }_1^2{\beta }_1-3{\gamma }_1\right)\right)/\left(2{\alpha }_1+2{\beta }_1+2{\gamma }_1\right)\right)t}}.\end{array}$$

By choosing the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 4 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-1$, $x=0$, and $x=1$.

[figure omitted; refer to PDF]

4.2. Set II: Two-Wave Solutions

We commence with two-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (2) has the rational function of two-wave solutions as shown in the following form: $$\begin{array}{c}\text{(27)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_2}{{\Omega }_2},\text{(28)}& {\Omega }_2=1+{\sigma }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t}+{\sigma }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y+{\gamma }_{2}z-{\delta }_{2}t}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y+\left({\gamma }_1+{\gamma }_2\right)z-\left({\delta }_1+{\delta }_2\right)t},{\Delta }_2={\rho }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y+{\gamma }_{1}z-{\delta }_{1}t}+{\rho }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y+{\gamma }_{2}z-{\delta }_{2}t}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y+\left({\gamma }_1+{\gamma }_2\right)z-\left({\delta }_1+{\delta }_2\right)t}.\end{array}$$

Plugging (27) along with (28) into equation (2), we get to the following case: $$\begin{array}{c}\text{(29)}& {\sigma }_i={\sigma }_i,{\delta }_i=\frac{{\alpha }_i^3{\beta }_i}{2{\alpha }_i+2{\beta }_i},\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,\\ {\rho }_{12}=\frac{{\Sigma }_{12}}{2{\alpha }_1{\alpha }_2{\Gamma }_{12}},\\ {\sigma }_{12}=\frac{{\text{Y}}_{12}}{{\Phi }_{12}}.\\ \text{(30)}& {\Sigma }_{12}={\alpha }_1{\alpha }_2\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1+3{\beta }_1^2{\beta }_2-3{\beta }_1{\beta }_2^2\right)\left({\alpha }_1^2-{\alpha }_2^2\right)+{\alpha }_2^2{\beta }_1^2\left({\alpha }_1+{\alpha }_2\right)\left(3{\alpha }_1^2-3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+2{\alpha }_1{\alpha }_2{\beta }_1{\beta }_2\left({\alpha }_1+{\alpha }_2\right)\left({\alpha }_1^2-3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+{\alpha }_1^2{\beta }_2^2\left({\alpha }_1+{\alpha }_2\right)\left({\alpha }_1^2-3{\alpha }_1{\alpha }_2+3{\alpha }_2^2\right),{\Gamma }_{12}={\alpha }_1{\alpha }_2\left({\alpha }_1+{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1+3{\beta }_1^2{\beta }_2+3{\beta }_1{\beta }_2^2\right)+{\alpha }_2^2{\beta }_1^2\left(3{\alpha }_1^2+3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+2{\alpha }_1{\alpha }_2{\beta }_1{\beta }_2\left({\alpha }_1^2+3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+{\alpha }_1^2{\beta }_2^2\left({\alpha }_1^2+3{\alpha }_1{\alpha }_2+3{\alpha }_2^2\right),\\ {\text{Y}}_{12}={\alpha }_1{\alpha }_2\left({\alpha }_1-{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1+3{\beta }_1^2{\beta }_2-3{\beta }_1{\beta }_2^2\right)+{\alpha }_2^2{\beta }_1^2\left(3{\alpha }_1^2-3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+2{\alpha }_1{\alpha }_2{\beta }_1{\beta }_2\left({\alpha }_1^2-3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+{\alpha }_1^2{\beta }_2^2\left({\alpha }_1^2-3{\alpha }_1{\alpha }_2+3{\alpha }_2^2\right),\\ {\Phi }_{12}={\alpha }_1{\alpha }_2\left({\alpha }_1+{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1+3{\beta }_1^2{\beta }_2+3{\beta }_1{\beta }_2^2\right)+{\alpha }_2^2{\beta }_1^2\left(3{\alpha }_1^2+3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+2{\alpha }_1{\alpha }_2{\beta }_1{\beta }_2\left({\alpha }_1^2+3{\alpha }_1{\alpha }_2+{\alpha }_2^2\right)+{\alpha }_1^2{\beta }_2^2\left({\alpha }_1^2+3{\alpha }_1{\alpha }_2+3{\alpha }_2^2\right).\end{array}$$

Therefore, the resulting two-wave solution reads as $$\begin{array}{}\text{(31)}& {\Psi }_2\left(x,y,z,t\right)=\frac{2{\alpha }_1{\sigma }_1{e}^{{\Lambda }_1}+2{\alpha }_2{\sigma }_2{e}^{{\Lambda }_2}+2{\sigma }_1{\sigma }_2\left({\Sigma }_{12}/{\Gamma }_{12}\right)e^{{\Lambda }_1+{\Lambda }_2}}{1+{\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_1{\sigma }_2\left({\text{Y}}_{12}/{\Phi }_{12}\right)e^{{\Lambda }_1+{\Lambda }_2}},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y+{\gamma }_{i}z-\left({\alpha }_i^3{\beta }_i/\left(2{\alpha }_i+2{\beta }_i\right)\right)t,i=1,2$. By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 5 including the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

4.3. Set III: Triple-Wave Solutions

We commence with three-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (2) has the rational function of triple-wave solutions as shown in the following form: $$\begin{array}{c}\text{(32)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\hspace{1em}\text{(33)}& {\Omega }_3=1+{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+{\rho }_3{e}^{{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\rho }_1{\rho }_3{\rho }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\rho }_2{\rho }_3{\rho }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_3{\rho }_{12}{\rho }_{13}{\rho }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},{\Delta }_3={\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_3{e}^{{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\sigma }_1{\sigma }_3{\sigma }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\sigma }_2{\sigma }_3{\sigma }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_3{\sigma }_{12}{\sigma }_{13}{\sigma }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y+{\gamma }_{i}x-{\delta }_{i}t,i=1,2,3$. Plugging (32) along with (33) into equation (2), we obtain the following case: $$\begin{array}{c}\text{(34)}& {\sigma }_i={\sigma }_i,{\delta }_i=\frac{{\alpha }_i^3{\beta }_i}{2{\alpha }_i+2{\beta }_i},\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,3,\\ {\rho }_{12}=\frac{{\Sigma }_{12}}{2{\alpha }_1{\alpha }_2{\Gamma }_{12}},\\ {\sigma }_{12}=\frac{{\text{Y}}_{12}}{{\Phi }_{12}}.\\ \text{(35)}& {\rho }_{13}=\frac{{\Sigma }_{13}}{2{\alpha }_1{\alpha }_3{\Gamma }_{13}},{\sigma }_{13}=\frac{{\text{Y}}_{13}}{{\Phi }_{13}},\\ {\rho }_{23}=\frac{{\Sigma }_{23}}{2{\alpha }_2{\alpha }_3{\Gamma }_{23}},\\ {\sigma }_{23}=\frac{{\text{Y}}_{23}}{{\Phi }_{23}},\\ \hspace{1em}{\Sigma }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\Sigma }_{12},\\ \hspace{1em}{\Sigma }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\Sigma }_{13},\\ \hspace{1em}{\Gamma }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\Gamma }_{12},\\ \hspace{1em}{\Gamma }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\Gamma }_{13},\\ \hspace{1em}{\text{Y}}_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\text{Y}}_{12},\\ \hspace{1em}{\text{Y}}_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\text{Y}}_{13},\\ \hspace{1em}{\Phi }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\Phi }_{12},\\ \hspace{1em}{\Phi }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\Phi }_{13}.\end{array}$$

Therefore, the resulting three-wave solution reads as $$\begin{array}{}\text{(36)}& {\Psi }_3\left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\end{array}$$in which ${\Delta }_3$ and ${\Omega }_3$ with their relations are given in (33) and (34). By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 6 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

5. Multiple Soliton Solutions for the ($2+1$) eCBS Equation

5.1. Set I: One-Wave Solution

We start with one-wave function based on the statement in Step 2 in the previous section; we suppose that equation (3) has the rational function of the one-wave solution as shown in the following form: $$\begin{array}{}\text{(37)}& \Psi \left(x,y,t\right)=\frac{{\Delta }_1}{{\Omega }_1},\hspace{1em}{\Omega }_1={\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t},{\Delta }_1={\sigma }_1+{\sigma }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t},\end{array}$$in which ${\sigma }_1$ and ${\sigma }_2$ are unfound constants. Substituting (37) into equation (3), we get to the following case: $$\begin{array}{c}\text{(38)}& {\rho }_1={\rho }_1,{\rho }_2={\rho }_2,\\ {\sigma }_1={\sigma }_1,\\ {\sigma }_2=\frac{{\rho }_2\left({\alpha }_1{\rho }_1+{\sigma }_1\right)}{{\rho }_1},\\ {\lambda }_1={\alpha }_1^3{\delta }_3+{\alpha }_1^2{\beta }_1+{\alpha }_1{\delta }_2+{\beta }_1{\delta }_1.\end{array}$$

Therefore, the resulting one-wave solution reads as $$\begin{array}{}\text{(39)}& \Psi \left(x,y,t\right)=\frac{{\sigma }_1+\left({\rho }_2\left({\alpha }_1{\rho }_1+{\sigma }_1\right)/{\rho }_1\right)e^{{\alpha }_{1}x+{\beta }_{1}y-\left({\alpha }_1^3{\delta }_3+{\alpha }_1^2{\beta }_1+{\alpha }_1{\delta }_2+{\beta }_1{\delta }_1\right)t}}{{\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-\left({\alpha }_1^3{\delta }_3+{\alpha }_1^2{\beta }_1+{\alpha }_1{\delta }_2+{\beta }_1{\delta }_1\right)t}}.\end{array}$$

By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 7 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-1$, $x=0$, and $x=1$.

[figure omitted; refer to PDF]

5.2. Set II: Two-Wave Solutions

We begin with two-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (3) has the rational function of two-wave solutions as shown in the following form: $$\begin{array}{c}\text{(40)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_2}{{\Omega }_2},\hspace{1em}\text{(41)}& {\Omega }_2=1+{\sigma }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t}+{\sigma }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y-{\lambda }_{2}t}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y+\left({\gamma }_1-\left({\lambda }_1+{\lambda }_2\right)\right.t},{\Delta }_2={\rho }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t}+{\rho }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y-{\lambda }_{2}t}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y-\left({\lambda }_1+{\lambda }_2\right)t}.\end{array}$$

Inserting (40) along with (41) into equation (3), we obtain the following case: $$\begin{array}{c}\text{(42)}& {\sigma }_i={\sigma }_i,{\lambda }_i={\alpha }_i^3{\delta }_3+{\alpha }_i^2{\beta }_i+{\alpha }_i{\delta }_2+{\beta }_i{\delta }_1,\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,\\ {\rho }_{12}=\frac{3{\alpha }_1{\alpha }_2{\delta }_3{\left({\alpha }_1-{\alpha }_2\right)}^2+\left({\alpha }_1-{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)}{{\alpha }_2{\alpha }_1\left(3{\alpha }_1{\alpha }_2{\delta }_3\left({\alpha }_1+{\alpha }_2\right)+2{\alpha }_1{\alpha }_2\left({\beta }_1+{\beta }_2\right)+{\alpha }_1^2{\beta }_2+{\alpha }_2^2{\beta }_1\right)}.\\ \text{(43)}& {\sigma }_{12}=\frac{3{\alpha }_1{\alpha }_2{\delta }_3{\left({\alpha }_1-{\alpha }_2\right)}^2+\left({\alpha }_1-{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)}{3{\alpha }_1{\alpha }_2{\delta }_3{\left({\alpha }_1+{\alpha }_2\right)}^2+\left({\alpha }_1+{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1\right)}.\end{array}$$

Therefore, the resulting two-wave solution reads as $$\begin{array}{}\text{(44)}& {\Psi }_2\left(x,y,t\right)=\frac{{\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}}{{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+4{\alpha }_1{\alpha }_2{\sigma }_1{\sigma }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y-\left({{\alpha }_i}^3{\delta }_3+{{\alpha }_i}^2{\beta }_i+{\alpha }_i{\delta }_2+{\beta }_i{\delta }_1\right)t,i=1,2$. By choosing the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 8 including the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-1$, $x=0$, and $x=1$.

[figure omitted; refer to PDF]

5.3. Set III: Triple-Wave Solutions

We commence with three-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (3) has the rational function of triple-wave solutions as shown in the following form: $$\begin{array}{c}\text{(45)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\text{(46)}& {\Omega }_3=1+{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+{\rho }_3{e}^{{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\rho }_1{\rho }_3{\rho }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\rho }_2{\rho }_3{\rho }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_3{\rho }_{12}{\rho }_{13}{\rho }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},{\Delta }_3={\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_3{e}^{{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\sigma }_1{\sigma }_3{\sigma }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\sigma }_2{\sigma }_3{\sigma }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_3{\sigma }_{12}{\sigma }_{13}{\sigma }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y-{\lambda }_{i}t,i=1,2,3$. Plugging (45) along with (46) into equation (3), we obtain the following case: $$\begin{array}{c}\text{(47)}& {\sigma }_i={\sigma }_i,{\lambda }_i={\alpha }_i^3{\delta }_3+{\alpha }_i^2{\beta }_i+{\alpha }_i{\delta }_2+{\beta }_i{\delta }_1,\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,3,\\ {\rho }_{12}=\frac{3{\alpha }_1{\alpha }_2{\delta }_3{\left({\alpha }_1-{\alpha }_2\right)}^2+\left({\alpha }_1-{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)}{{\alpha }_2{\alpha }_1\left(3{\alpha }_1{\alpha }_2{\delta }_3\left({\alpha }_1+{\alpha }_2\right)+2{\alpha }_1{\alpha }_2\left({\beta }_1+{\beta }_2\right)+{\alpha }_1^2{\beta }_2+{\alpha }_2^2{\beta }_1\right)},\\ {\sigma }_{12}=\frac{3{\alpha }_1{\alpha }_2{\delta }_3{\left({\alpha }_1-{\alpha }_2\right)}^2+\left({\alpha }_1-{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)}{3{\alpha }_1{\alpha }_2{\delta }_3{\left({\alpha }_1+{\alpha }_2\right)}^2+\left({\alpha }_1+{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1\right)},\\ \hspace{1em}{\rho }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\rho }_{12},\\ \hspace{1em}{\sigma }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\sigma }_{12},\\ \hspace{1em}{\rho }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\rho }_{13},\\ \hspace{1em}{\sigma }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\sigma }_{13}.\end{array}$$

Therefore, the resulting three-wave solution reads as $$\begin{array}{}\text{(48)}& {\Psi }_3\left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3}.\end{array}$$

In (48), ${\Delta }_3$ and ${\Omega }_3$ with their relations are given in (46) and (47). By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 9 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-1$, $x=0$, and $x=1$.

[figure omitted; refer to PDF]

6. Multiple Soliton Solutions for the Generalized ($2+1$) BK Equation

6.1. Set I: One-Wave Solution

We start with one-wave function based on the statement in Step 2 in the previous section; we suppose that equation (4) has the rational function of the one-wave solution as shown in the following form: $$\begin{array}{}\text{(49)}& \Psi \left(x,y,t\right)=\frac{{\Delta }_1}{{\Omega }_1},\hspace{1em}{\Omega }_1={\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t},{\Delta }_1={\sigma }_1+{\sigma }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t},\end{array}$$in which ${\sigma }_1$ and ${\sigma }_2$ are unfound constants. Putting (49) into equation (4), we get to the following case: $$\begin{array}{c}\text{(50)}& {\rho }_1={\rho }_1,{\rho }_2={\rho }_2,\\ {\sigma }_1={\sigma }_1,\\ {\sigma }_2=\frac{{\rho }_2\left(2{\alpha }_1{\rho }_1+{\sigma }_1\right)}{{\rho }_1},\\ {\lambda }_1=\frac{{\alpha \alpha }_1^4+\beta {\alpha }_1^3{\beta }_1+{\alpha }_1^2{\delta }_1+{\alpha }_1{\beta }_1{\delta }_2+{\beta }_1^2{\delta }_3}{{\alpha }_1}.\end{array}$$

Therefore, the resulting one-wave solution reads as $$\begin{array}{}\text{(51)}& \Psi \left(x,y,t\right)=\frac{{\sigma }_1+\left({\rho }_2\left(2{\alpha }_1{\rho }_1+{\sigma }_1\right)/{\rho }_1\right)e^{{\alpha }_{1}x+{\beta }_{1}y-\left({\alpha \alpha }_1^4+{\beta \alpha }_1^3{\beta }_1+{\alpha }_1^2{\delta }_1+{\alpha }_1{\beta }_1{\delta }_2+{\beta }_1^2{\delta }_3/{\alpha }_1\right)t}}{{\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-\left({\alpha \alpha }_1^4+{\beta \alpha }_1^3{\beta }_1+{\alpha }_1^2{\delta }_1+{\alpha }_1{\beta }_1{\delta }_2+{\beta }_1^2{\delta }_3/{\alpha }_1\right)t}}.\end{array}$$

By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 10 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

6.2. Set II: Two-Wave Solutions

We begin with two-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (4) has the rational function of two-wave solutions as shown in the following form: $$\begin{array}{c}\text{(52)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_2}{{\Omega }_2},\text{(53)}& {\Omega }_2=1+{\sigma }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t}+{\sigma }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y-{\lambda }_{2}t}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y+\left({\gamma }_1-\left({\lambda }_1+{\lambda }_2\right)\right.t},{\Delta }_2={\rho }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t}+{\rho }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y-{\lambda }_{2}t}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y-\left({\lambda }_1+{\lambda }_2\right)t}.\end{array}$$

Inserting (52) along with (53) into equation (4), we obtain the following case: $$\begin{array}{c}\text{(54)}& {\sigma }_i={\sigma }_i,{\lambda }_i=\frac{{\alpha \alpha }_i^4+\beta {\alpha }_i^3{\beta }_i+{\alpha }_i^2{\delta }_1+{\alpha }_i{\beta }_i{\delta }_2+{\beta }_i^2{\delta }_3}{{\alpha }_i},\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,\\ {\rho }_{12}=\frac{3{\alpha \alpha }_1^2{\alpha }_2^2\left({\alpha }_2+{\alpha }_1\right){\left({\alpha }_1-{\alpha }_2\right)}^2+\left({\alpha }_2+{\alpha }_1\right)\left({\beta \alpha }_1{\alpha }_2\left({\alpha }_1-{\alpha }_2\right)\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2\right)}{2{\alpha }_1{\alpha }_2\left({\alpha }_1{\alpha }_2\right.\left(\left.{\alpha }_2+{\alpha }_1\right)\left(3{\alpha \alpha }_1{\alpha }_2\left({\alpha }_2+{\alpha }_1\right)+\beta \left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1\right)\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2\right)},\\ {\sigma }_{12}=\frac{{\alpha }_1{\alpha }_2\left(-{\alpha }_2+{\alpha }_1\right)\left(3{\alpha \alpha }_1{\alpha }_2\left(-{\alpha }_2+{\alpha }_1\right)+\beta \left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}{{\alpha }_1{\alpha }_2\left({\alpha }_2+{\alpha }_1\right)\left(3{\alpha \alpha }_1{\alpha }_2\left({\alpha }_2+{\alpha }_1\right)+\beta \left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1\right)\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}.\end{array}$$

Therefore, the resulting two-wave solution reads as $$\begin{array}{}\text{(55)}& {\Psi }_2\left(x,y,t\right)=\frac{{\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}}{{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+4{\alpha }_1{\alpha }_2{\sigma }_1{\sigma }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y-\left({\alpha \alpha }_i^4+{\beta \alpha }_i^3{\beta }_i+{\alpha }_i^2{\delta }_1+{\alpha }_i{\beta }_i{\delta }_2+{\beta }_i^2{\delta }_3/{\alpha }_i\right)t,i=1,2$. By choosing suitable values of parameters, the analytical treatment of the periodic wave solution is presented in Figure 11 including the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

6.3. Set III: Triple-Wave Solutions

We commence with three-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (4) has the rational function of triple-wave solutions as shown in the following form: $$\begin{array}{c}\text{(56)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\text{(57)}& {\Omega }_3=1+{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+{\rho }_3{e}^{{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\rho }_1{\rho }_3{\rho }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\rho }_2{\rho }_3{\rho }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_3{\rho }_{12}{\rho }_{13}{\rho }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},{\Delta }_3={\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_3{e}^{{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\sigma }_1{\sigma }_3{\sigma }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\sigma }_2{\sigma }_3{\sigma }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_3{\sigma }_{12}{\sigma }_{13}{\sigma }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y-{\lambda }_{i}t,i=1,2,3$. Plugging (56) along with (57) into equation (4), we obtain the following case: $$\begin{array}{c}\text{(58)}& {\sigma }_i={\sigma }_i,{\lambda }_i=\frac{{\alpha \alpha }_i^4+{\beta \alpha }_i^3{\beta }_i+{\alpha }_i^2{\delta }_1+{\alpha }_i{\beta }_i{\delta }_2+{\beta }_i^2{\delta }_3}{{\alpha }_i},\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,3,\\ {\rho }_{12}=\frac{3{\alpha \alpha }_1^2{\alpha }_2^2\left({\alpha }_2+{\alpha }_1\right){\left({\alpha }_1-{\alpha }_2\right)}^2+\left({\alpha }_2+{\alpha }_1\right)\left({\beta \alpha }_1{\alpha }_2\left({\alpha }_1-{\alpha }_2\right)\right.\left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}{2{\alpha }_1{\alpha }_2\left({\alpha }_1{\alpha }_2\right)\left({\alpha }_2+{\alpha }_1\left(3{\alpha \alpha }_1{\alpha }_2\left({\alpha }_2+{\alpha }_1\right)+\beta \left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1\right)\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2\right)},\\ {\sigma }_{12}=\frac{{\alpha }_1{\alpha }_2\left(-{\alpha }_2+{\alpha }_1\right)\left(3{\alpha \alpha }_1{\alpha }_2\left(-{\alpha }_2+{\alpha }_1\right)+\beta \left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1-2{\alpha }_1{\alpha }_2{\beta }_2-{\alpha }_2^2{\beta }_1\right)\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}{{\alpha }_1{\alpha }_2\left({\alpha }_2+{\alpha }_1\right)\left(3{\alpha \alpha }_1{\alpha }_2\left({\alpha }_2+{\alpha }_1\right)+\beta \left({\alpha }_1^2{\beta }_2+2{\alpha }_1{\alpha }_2{\beta }_1+2{\alpha }_1{\alpha }_2{\beta }_2+{\alpha }_2^2{\beta }_1\right)\right)-{\delta }_3{\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2},\\ \hspace{1em}{\rho }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\rho }_{12},\\ \hspace{1em}{\sigma }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\sigma }_{12},\\ \hspace{1em}{\rho }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\rho }_{13},\\ \hspace{1em}{\sigma }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\sigma }_{13}.\end{array}$$

Therefore, the resulting three-wave solution reads as $$\begin{array}{}\text{(59)}& {\Psi }_3\left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\end{array}$$in which ${\Delta }_3$ and ${\Omega }_3$ with their relations are given in (57) and (58). By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 12 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

7. Multiple Soliton Solutions for the ($2+1$) VC DJKM Equation

7.1. Set I: One-Wave Solution

We start with one-wave function based on the statement in Step 2 in the previous section; we suppose that equation (6) has the rational function of one-wave solutions as shown in the following form: $$\begin{array}{}\text{(60)}& \Psi \left(x,y,t\right)=\frac{{\Delta }_1}{{\Omega }_1},\hspace{1em}{\Omega }_1={\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t},{\Delta }_1={\sigma }_1+{\sigma }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t},\end{array}$$in which ${\rho }_1,{\rho }_2,{\sigma }_1$ and ${\sigma }_2$ are unfound constants. Inserting (60) into equation (6), we gain the following case: $$\begin{array}{c}\text{(61)}& {\rho }_1={\rho }_1,{\rho }_2={\rho }_2,\\ {\sigma }_1=-\frac{{\rho }_1\left(2{\alpha }_1{\rho }_2-{\sigma }_2\right)}{{\rho }_2},\\ {\sigma }_2={\sigma }_2,\\ {\lambda }_1=\frac{{\beta }_1\left(-{\alpha }_1^4+\alpha {\beta }_1^2-h\left(t\right){\alpha }_1^2\right)}{2\beta g\left(t\right){\alpha }_1^2}.\end{array}$$

Therefore, the resulting one-wave solution reads as $$\begin{array}{}\text{(62)}& \Psi \left(x,y,t\right)=\frac{-\left({\rho }_1\left(2{\alpha }_1{\rho }_2-{\sigma }_2\right)/{\rho }_2\right)+{\sigma }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-\left({\beta }_1\left(-{\alpha }_1^4+\alpha {\beta }_1^2-h\left(t\right){\alpha }_1^2\right)/2\beta g\left(t\right){\alpha }_1^2\right)t}}{{\rho }_1+{\rho }_2{e}^{{\alpha }_{1}x+{\beta }_{1}y-\left({\beta }_1\left(-{\alpha }_1^4+\alpha {\beta }_1^2-h\left(t\right){\alpha }_1^2\right)/2\beta g\left(t\right){\alpha }_1^2\right)t}}.\end{array}$$

By choosing the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 13 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-5$, $x=0$, and $x=5$.

[figure omitted; refer to PDF]

7.2. Set II: Two-Wave Solutions

We begin with two-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (6) has the rational function of two-wave solutions as shown in the following form: $$\begin{array}{c}\text{(63)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_2}{{\Omega }_2},\text{(64)}& {\Omega }_2=1+{\sigma }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t}+{\sigma }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y-{\lambda }_{2}t}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y+\left({\gamma }_1-\left({\lambda }_1+{\lambda }_2\right)\right.t},{\Delta }_2={\rho }_1{e}^{{\alpha }_{1}x+{\beta }_{1}y-{\lambda }_{1}t}+{\rho }_2{e}^{{\alpha }_{2}x+{\beta }_{2}y-{\lambda }_{2}t}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{\left({\alpha }_1+{\alpha }_2\right)x+\left({\beta }_1+{\beta }_2\right)y-\left({\lambda }_1+{\lambda }_2\right)t}.\end{array}$$

Inserting (63) along with (64) into equation (6), we obtain the following case: $$\begin{array}{c}\text{(65)}& {\sigma }_i={\sigma }_i,{\lambda }_i=\frac{{\beta }_i\left(-{\alpha }_i^4+\alpha {\beta }_i^2-h\left(t\right){\alpha }_i^2\right)}{2\beta g\left(t\right){\alpha }_i^2},\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2.\\ \text{(66)}& {\rho }_{12}=\frac{1}{2}\frac{{\alpha }_1^2{\alpha }_2^2\left({\alpha }_1+{\alpha }_2\right){\left(-{\alpha }_2+{\alpha }_1\right)}^2+\alpha \left({\alpha }_1+{\alpha }_2\right){\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}{{\alpha }_2{\alpha }_1\left({\alpha }_1^2{\alpha }_2^2{\left({\alpha }_1+{\alpha }_2\right)}^2+\alpha {\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2\right)},{\sigma }_{12}=\frac{{\alpha }_1^2{\alpha }_2^2{\left(-{\alpha }_2+{\alpha }_1\right)}^2+\alpha {\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}{{\alpha }_1^2{\alpha }_2^2{\left({\alpha }_1+{\alpha }_2\right)}^2+\alpha {\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}.\end{array}$$

Therefore, the resulting two-wave solution reads as $$\begin{array}{}\text{(67)}& {\Psi }_2\left(x,y,t\right)=\frac{{\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}}{{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+4{\alpha }_1{\alpha }_2{\sigma }_1{\sigma }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y-\left({\beta }_i\left(-{\alpha }_i^4+{\alpha \beta }_i^2-h\left(t\right){\alpha }_i^2\right)/2\beta g\left(t\right){\alpha }_i^2\right)t,i=1,2$. By choosing the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 14 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

7.3. Set III: Triple-Wave Solutions

We commence with three-wave functions based on the statement in Step 2 in the previous section; we suppose that equation (6) has the rational function of triple-wave solutions as shown in the following form: $$\begin{array}{c}\text{(68)}& \Psi \left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3},\text{(69)}& {\Omega }_3=1+{\rho }_1{e}^{{\Lambda }_1}+{\rho }_2{e}^{{\Lambda }_2}+{\rho }_3{e}^{{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\rho }_1{\rho }_3{\rho }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\rho }_2{\rho }_3{\rho }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\rho }_1{\rho }_2{\rho }_3{\rho }_{12}{\rho }_{13}{\rho }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},{\Delta }_3={\sigma }_1{e}^{{\Lambda }_1}+{\sigma }_2{e}^{{\Lambda }_2}+{\sigma }_3{e}^{{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_{12}{e}^{{\Lambda }_1+{\Lambda }_2}+{\sigma }_1{\sigma }_3{\sigma }_{13}{e}^{{\Lambda }_1+{\Lambda }_3}+{\sigma }_2{\sigma }_3{\sigma }_{23}{e}^{{\Lambda }_2+{\Lambda }_3}+{\sigma }_1{\sigma }_2{\sigma }_3{\sigma }_{12}{\sigma }_{13}{\sigma }_{23}{e}^{{\Lambda }_1+{\Lambda }_2+{\Lambda }_3},\end{array}$$in which ${\Lambda }_i={\alpha }_{i}x+{\beta }_{i}y-{\lambda }_{i}t,i=1,2,3$. Plugging (68) along with (69) into equation (6), we obtain the following case: $$\begin{array}{c}\text{(70)}& {\sigma }_i={\sigma }_i,{\lambda }_i=\frac{{\beta }_i\left(-{\alpha }_i^4+{\alpha \beta }_i^2-h\left(t\right){\alpha }_i^2\right)}{2\beta g\left(t\right){\alpha }_i^2},\\ {\rho }_i=2{\alpha }_i{\sigma }_i,\\ \hspace{1em}i=1,2,3,\\ {\rho }_{12}=\frac{1}{2}\frac{{\alpha }_1^2{\alpha }_2^2\left({\alpha }_1+{\alpha }_2\right){\left(-{\alpha }_2+{\alpha }_1\right)}^2+\alpha \left({\alpha }_1+{\alpha }_2\right){\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}{{\alpha }_2{\alpha }_1\left({\alpha }_1^2{\alpha }_2^2{\left({\alpha }_1+{\alpha }_2\right)}^2+\alpha {\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2\right)},\\ {\sigma }_{12}=\frac{{\alpha }_1^2{\alpha }_2^2{\left(-{\alpha }_2+{\alpha }_1\right)}^2+\alpha {\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2}{{\alpha }_1^2{\alpha }_2^2{\left({\alpha }_1+{\alpha }_2\right)}^2+\alpha {\left({\alpha }_1{\beta }_2-{\alpha }_2{\beta }_1\right)}^2},\\ \hspace{1em}{\rho }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\rho }_{12},\\ \hspace{1em}{\sigma }_{13}=\text{substitute}\left\{{\alpha }_2={\alpha }_3,{\beta }_2={\beta }_3\right\},\text{at}{\sigma }_{12},\\ \hspace{1em}{\rho }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\rho }_{13},\\ \hspace{1em}{\sigma }_{23}=\text{substitute}\left\{{\alpha }_1={\alpha }_2,{\beta }_1={\beta }_2\right\},\text{at}{\sigma }_{13}.\end{array}$$

Therefore, the resulting three-wave solution reads as $$\begin{array}{}\text{(71)}& {\Psi }_3\left(x,y,z,t\right)=\frac{{\Delta }_3}{{\Omega }_3}.\end{array}$$

In (71), ${\Delta }_3$ and ${\Omega }_3$ with their relations are given in (69) and (70). By selecting the suitable values of parameters, the graphic representation of the periodic wave solution is presented in Figure 15 containing the 3D plot, the contour plot, the density plot, and the 2D plot when three spaces arise at spaces $x=-10$, $x=0$, and $x=10$.

[figure omitted; refer to PDF]

This paper finds many novel one-soliton-, two-soliton-, and triple-soliton-type solutions to governing models. With the help of some computations, surfaces of results reported have been observed in Figures 1–15. These figures are dependent on the family conditions which are of importance physically. It has been investigated that all figures plotted have been symbolized for the four types of the nonlinear PDEs containing the extended ($3+1$)-dimensional Jimbo-Miwa-like equation, the extended ($2+1$)-dimensional Calogero-Bogoyavlenskii-Schiff equation, the generalization of the ($2+1$)-dimensional Bogoyavlensky-Konopelchenko equation, and a variable-coefficient extension of the DJKM equation. These mathematical properties come from exponential function properties. In this sense, from the mathematical and physical points of view, these results play an important role in explaining wave propagation of nonlinear phenomena. Hence, we consider that surfaces plotted in this paper have proven such physical meaning of the obtained solutions.

8. Conclusion

In this article, we obtained the multiple soliton solutions of the novel extended ($3+1$)-dimensional Jimbo-Miwa-like equation, the extended ($2+1$)-dimensional Calogero-Bogoyavlenskii-Schiff equation, the generalization of the ($2+1$)-dimensional Bogoyavlensky-Konopelchenko equation, and a variable-coefficient extension of the DJKM equation via operating the multiple Exp-function method, including one-soliton-, two-soliton-, and triple-soliton-type solutions. It is quite visible that this novel scheme has plenty of family of solutions containing rational exponential functions by selecting particular parameters. Thus, this paper provides a lot of encouragement for future research in soliton topics. The behaviours of the solutions for the known nonlinear equations obtained by the multiple Exp-function method by choosing the suitable values are cited in Figures 1–15. Moreover, the numerical simulations include gaining the coefficients${\rho }_i,{\sigma }_i,\text{and}{\lambda }_i$ and have been carried out to show that the projected algorithm is applicable and efficient. The analytical study has been conducted for the solutions $\Psi \left(x,y,t\right)$ obtained by employing the aforementioned method, in which obtained solutions were as one-soliton, two-soliton, and triple-soliton. The results are beneficial to the study of wave propagation. All computations in this paper have been made quickly with the aid of Maple.

References

[1] S. H. Seyedi, B. N. Saray, A. J. Chamkha, "Heat and mass transfer investigation of MHD Eyring-Powell flow in a stretching channel with chemical reactions," Physica A: Statistical Mechanics and its Applications, vol. 544,DOI: 10.1016/j.physa.2019.124109, 2020.

[2] M. Dehghan, J. Manafian, A. Saadatmandi, "Analytical treatment of some partial differential equations arising in mathematical physics by using the Exp-function method," International Journal of Modern Physics B, vol. 25, pp. 2965-2981, 2012.

[3] M. Dehghan, J. Manafian, A. Saadatmandi, "Solving nonlinear fractional partial differential equations using the homotopy analysis method," Numerical Methods for Partial Differential Equations Journal, vol. 26, pp. 448-479, 2010.

[4] M. Dehghan, J. Manafian, "The solution of the variable coefficients fourth-order parabolic partial differential equations by the homotopy perturbation method," Zeitschrift für Naturforschung A, vol. 64 no. 7-8, pp. 420-430, DOI: 10.1515/zna-2009-7-803, 2009.

[5] J. Manafian, M. Lakestani, "Solitary wave and periodic wave solutions for Burgers, Fisher, Huxley and combined forms of these equations by the ( G ′ / G )-expansion method," Pramana, vol. 130, pp. 31-52, 2015.

[6] J. Manafian, M. Lakestani, "Application of tan( ϕ /2)-expansion method for solving the Biswas–Milovic equation for Kerr law nonlinearity," Optik, vol. 127 no. 4, pp. 2040-2054, DOI: 10.1016/j.ijleo.2015.11.078, 2016.

[7] J. Manafian, M. Lakestani, "Dispersive dark optical soliton with Tzitzéica type nonlinear evolution equations arising in nonlinear optics," Optical and Quantum Electronics, vol. 48 no. 2,DOI: 10.1007/s11082-016-0371-y, 2016.

[8] A. R. Seadawy, J. Manafian, "New soliton solution to the longitudinal wave equation in a magneto-electro-elastic circular rod," Results in Physics, vol. 8, pp. 1158-1167, DOI: 10.1016/j.rinp.2018.01.062, 2018.

[9] J. Manafian, "Novel solitary wave solutions for the ( 3 + 1 )-dimensional extended Jimbo-Miwa equations," Computers & Mathematcs with Applications, vol. 76 no. 5, pp. 1246-1260, DOI: 10.1016/j.camwa.2018.06.018, 2018.

[10] J. Manafian, B. Mohammadi-Ivatlo, M. Abapour, "Lump-type solutions and interaction phenomenon to the ( 2 + 1 )-dimensional breaking soliton equation," Applied Mathematics and Computation, vol. 13, pp. 13-41, 2019.

[11] W. X. Ma, Y. Zhou, R. Dougherty, "Lump-type solutions to nonlinear differential equations derived from generalized bilinear equations," International Journal of Modern Physics B, vol. 30 no. 28n29, article 1640018,DOI: 10.1142/S021797921640018X, 2016.

[12] J. Lü, S. Bilige, X. Gao, Y. Bai, R. Zhang, "Abundant lump solutions and interaction phenomena to the Kadomtsev-Petviashvili-Benjamin-Bona-Mahony equation," Journal of Applied Mathematics and Physics, vol. 6 no. 8, pp. 1733-1747, DOI: 10.4236/jamp.2018.68148, 2018.

[13] C. Wang, "Spatiotemporal deformation of lump solution to ( 2 + 1 )-dimensional KdV equation," Nonlinear Dynamics, vol. 84 no. 2, pp. 697-702, DOI: 10.1007/s11071-015-2519-x, 2016.

[14] J. H. He, "A modified Li-He’s variational principle for plasma," International Journal of Numerical Methods for Heat and Fluid Flow,DOI: 10.1108/HFF-06-2019-0523, 2019.

[15] J. H. He, "Lagrange crisis and generalized variational principle for 3D unsteady flow," International Journal of Numerical Methods for Heat and Fluid Flow, vol. 30 no. 3, pp. 1189-1196, DOI: 10.1108/HFF-07-2019-0577, 2019.

[16] S. S. Chen, B. Tian, L. Liu, Y. Q. Yuan, C. R. Zhang, "Conservation laws, binary Darboux transformations and solitons for a higher-order nonlinear Schrödinger system," Chaos, Solitons and Fractals, vol. 118, pp. 337-346, DOI: 10.1016/j.chaos.2018.11.010, 2019.

[17] X. X. Du, B. Tian, X. Y. Wu, H. M. Yin, C. R. Zhang, "Lie group analysis, analytic solutions and conservation laws of the ( 3 + 1 )-dimensional Zakharov-Kuznetsov-Burgers equation in a collisionless magnetized electron-positron-ion plasma," The European Physical Journal Plus, vol. 133 no. 9,DOI: 10.1140/epjp/i2018-12239-y, 2018.

[18] S. Saha Ray, "On conservation laws by Lie symmetry analysis for ( 2 + 1 )-dimensional BogoyavlenskyKonopelchenko equation in wave propagation," Computers & Mathematics with Applications, vol. 74, pp. 1158-1165, 2017.

[19] X. H. Zhao, B. Tian, X. Y. Xie, X. Y. Wu, Y. Sun, Y. J. Guo, "Solitons, Bäcklund transformation and Lax pair for a ( 2 + 1 )-dimensional Davey-Stewartson system on surface waves of finite depth," Wave Random Complex, vol. 28, pp. 356-366, 2017.

[20] H. M. Baskonus, H. Bulut, "Exponential prototype structures for ( 2 + 1 )-dimensional Boiti-Leon-Pempinelli systems in mathematical physics," Waves in Random and Complex Media, vol. 26, pp. 201-208, 2016.

[21] H. M. Baskonus, D. A. Koç, H. Bulut, "New travelling wave prototypes to the nonlinear Zakharov-Kuznetsov equation with power law nonlinearity," Nonlinear Analysis: Real World Applications, vol. 7, pp. 67-76, DOI: 10.1016/j.nonrwa.2009.08.007, 2010.

[22] Q. Zhou, "Optical solitons in medium with parabolic law nonlinearity and higher order dispersion," Waves in Random and Complex Media, vol. 25, pp. 52-59, 2016.

[23] T. A. Sulaiman, R. I. Nuruddeen, E. Zerrad, B. B. Mikail, "Dark and singular solitons to the two nonlinear Schrödinger equations," Optik, vol. 186, pp. 423-430, DOI: 10.1016/j.ijleo.2019.04.023, 2019.

[24] M. Inc, A. I. Aliyu, A. Yusuf, D. Baleanu, "Optical solitary waves, conservation laws and modulation instability analysis to the nonlinear Schrödinger’s equation in compressional dispersive Alvèn waves," Optik, vol. 155, pp. 257-266, DOI: 10.1016/j.ijleo.2017.10.109, 2018.

[25] M. Jimbo, "Solitons and infinite-dimensional lie algebras," Kyoto University, vol. 19 no. 3, pp. 943-1001, DOI: 10.2977/prims/1195182017, 1983.

[26] W. Tan, Z. D. Dai, J. L. Xie, L. L. Hu, "Emergence and interaction of the lump-type solution with the ( 3 + 1 )-D Jimbo-Miwa equation," Zeitschrift für Naturforschung, vol. 73 no. 1, pp. 43-49, DOI: 10.1515/zna-2017-0293, 2017.

[27] J. Y. Yang, W. X. Ma, "Abundant lump-type solutions of the Jimbo-Miwa equation in ( 3 + 1 )-dimensions," Computers & Mathematcs with Applications, vol. 73 no. 2, pp. 220-225, DOI: 10.1016/j.camwa.2016.11.007, 2017.

[28] W. X. Ma, "Lump-type solutions to the ( 3 + 1 )-dimensional Jimbo-Miwa equation," International Journal of Nonlinear Sciences and Numerical Simulation, vol. 17 no. 7-8,DOI: 10.1515/ijnsns-2015-0050, 2016.

[29] W. X. Ma, J. H. Lee, "A transformed rational function method and exact solutions to the 3 + 1 dimensional Jimbo-Miwa equation," Chaos, Solitons & Fractals, vol. 42 no. 3, pp. 1356-1363, DOI: 10.1016/j.chaos.2009.03.043, 2009.

[30] Y. N. Tang, W. X. Ma, W. Xu, L. Gao, "Wronskian determinant solutions of the ( 3 + 1 )-dimensional Jimbo-Miwa equation," Applied Mathematics and Computation, vol. 217 no. 21, pp. 8722-8730, DOI: 10.1016/j.amc.2011.03.120, 2011.

[31] S. Batwa, W. X. Ma, "A study of lump-type and interaction solutions to a (3+1)-dimensional Jimbo-Miwa-like equation," Computers & Mathematcs with Applications, vol. 76 no. 7, pp. 1576-1582, DOI: 10.1016/j.camwa.2018.07.008, 2018.

[32] F. H. Qi, Y. H. Huang, P. Wang, "Solitary-wave and new exact solutions for an extended ( 3 + 1 )-dimensional Jimbo-Miwa-like equation," Applied Mathematics Letters, vol. 100,DOI: 10.1016/j.aml.2019.106004, 2020.

[33] H. Q. Sun, A. H. Chen, "Lump and lump-kink solutions of the ( 3 + 1 )-dimensional Jimbo-Miwa and two extended Jimbo-Miwa equations," Applied Mathematics Letters, vol. 68, pp. 55-61, DOI: 10.1016/j.aml.2016.12.008, 2017.

[34] Y. F. Yue, L. L. Huang, Y. Chen, "Localized waves and interaction solutions to an extended ( 3 + 1 )-dimensional Jimbo-Miwa equation," Applied Mathematics Letters, vol. 89, pp. 70-77, DOI: 10.1016/j.aml.2018.09.020, 2019.

[35] B. Ren, J. Lin, Z. M. Lou, "A new nonlinear equation with lump-soliton, lump-periodic, and lump-periodic-soliton solutions," Complexity, vol. 2019,DOI: 10.1155/2019/4072754, 2019.

[36] A. M. Wazwaz, "The ( 2 + 1 ) and ( 3 + 1 )-dimensional CBS equations: multiple soliton solutions and multiple singular soliton solutions," Zeitschrift für Naturforschung A, vol. 65 no. 3, pp. 173-181, DOI: 10.1515/zna-2010-0304, 2010.

[37] W. R. Sun, B. Tian, H. L. Zhen, Y. Sun, "Breathers and rogue waves of the fifth-order nonlinear Schrödinger equation in the Heisenberg ferromagnetic spin chain," Nonlinear Dynamics, vol. 81 no. 1-2, pp. 725-732, DOI: 10.1007/s11071-015-2022-4, 2015.

[38] Y. Zhang, J. Pang, "Lump and lump-type solutions of the generalized ( 3 + 1 )-dimensional variable-coefficient B-type Kadomtsev-Petviashvili equation," Journal of Applied Mathematics, vol. 2019,DOI: 10.1155/2019/7172860, 2019.

[39] H. Leblond, D. Mihalache, "Models of few optical cycle solitons beyond the slowly varying envelope approximation," Physics Reports, vol. 523 no. 2, pp. 61-126, DOI: 10.1016/j.physrep.2012.10.006, 2013.

[40] D. Mihalache, "Multidimensional localized structures in optical and matter-wave media: a topical survey of recent literature," Romanian Reports in Physics, vol. 69, 2017.

[41] P. A. Clarkson, M. D. Kruskal, "New similarity reductions of the Boussinesq equation," Journal of Mathematical Physics, vol. 30 no. 10, pp. 2201-2213, DOI: 10.1063/1.528613, 1989.

[42] J. Weiss, M. Tabor, G. Carnevale, "The Painlevé property for partial differential equations," Journal of Mathematical Physics, vol. 24 no. 3, pp. 522-526, DOI: 10.1063/1.525721, 1983.

[43] K. R. Adem, C. M. Khalique, "Exact solutions and conservation laws of Zakharov-Kuznetsov modified equal width equation with power law nonlinearity," Nonlinear Analysis: Real World Applications, vol. 13 no. 4, pp. 1692-1702, DOI: 10.1016/j.nonrwa.2011.12.001, 2012.

[44] A. M. Wazwaz, Partial Differential Equations and Solitary Waves Theorem, 2009.

[45] A. M. Wazwaz, "Negative-order forms for the Calogero-Bogoyavlenskii-Schiff equation and the modified Calogero-Bogoyavlenskii-Schiff equation," Proceedings of the Romanian Academy. Series A, vol. 18 no. 4, pp. 337-344, 2017.

[46] S. Chen, F. Baronio, J. M. Soto-Crespo, P. Grelu, D. Mihalache, "Versatile rogue waves in scalar, vector, and multidimensional nonlinear systems," Journal of Physics A, vol. 50 no. 46,DOI: 10.1088/1751-8121/aa8f00, 2017.

[47] Y. V. Kartashov, G. E. Astrakharchik, B. A. Malomed, L. Torner, "Frontiers in multidimensional self-trapping of nonlinear fields and matter," Nature Reviews Physics, vol. 1 no. 3, pp. 185-197, DOI: 10.1038/s42254-019-0025-7, 2019.

[48] A. M. Wazwaz, "A ( 2 + 1 )-dimensional time-dependent Date-Jimbo-Kashiwara-Miwa equation: Painlevé integrability and multiple soliton solutions," Computers & Mathematcs with Applications, vol. 79 no. 4, pp. 1145-1149, DOI: 10.1016/j.camwa.2019.08.025, 2020.

[49] W. X. Ma, T. Huang, Y. Zhang, "A multiple exp-function method for nonlinear differential equations and its application," Physica Scripta, vol. 82 no. 6, article 065003,DOI: 10.1088/0031-8949/82/06/065003, 2010.

[50] W. X. Ma, Y. Zhou, "Lump solutions to nonlinear partial differential equations via Hirota bilinear forms," Journal of Differential Equations, vol. 264 no. 4, pp. 2633-2659, DOI: 10.1016/j.jde.2017.10.033, 2018.

[51] W. X. Ma, "Interaction solutions to Hirota-Satsuma-Ito equation in ( 2 + 1 )-dimensions," Frontiers of Mathematics in China, vol. 14 no. 3, pp. 619-629, DOI: 10.1007/s11464-019-0771-y, 2019.

[52] W. X. Ma, "Lump and interaction solutions to linear PDEs in 2 + 1 dimensions via symbolic computation," Modern Physics Letters B, vol. 33 no. 36,DOI: 10.1142/S0217984919504578, 2019.

[53] A. Cordero, J. P. Jaiswal, J. R. Torregrosa, "Stability analysis of fourth-order iterative method for finding multiple roots of non-linear equations," Applied Mathematics and Nonlinear Sciences, vol. 4 no. 1, pp. 43-56, DOI: 10.2478/AMNS.2019.1.00005, 2019.

[54] W. Gao, H. F. Ismael, A. M. Husien, H. Bulut, H. M. Baskonus, "Optical soliton solutions of the cubic-quartic nonlinear Schrödinger and resonant nonlinear Schrödinger equation with the parabolic law," Applied Sciences, vol. 10 no. 1,DOI: 10.3390/app10010219, 2020.

[55] Z. Zhao, B. Han, "Residual symmetry, Bäcklund transformation and CRE solvability of a ($$\mathbf{2}{\varvec{+}}{} \mathbf{1}$$ 2 + 1)-dimensional nonlinear system," Nonlinear Dynamics, vol. 94 no. 1, pp. 461-474, DOI: 10.1007/s11071-018-4371-2, 2018.

[56] W. Gao, G. Yel, H. M. Baskonus, C. Cattani, "Complex solitons in the conformable ( 2 + 1 )-dimensional Ablowitz-Kaup-Newell-Segur equation," Aims Mathematics, vol. 5 no. 1, pp. 507-521, DOI: 10.3934/math.2020034, 2020.

[57] W. Gao, H. Rezazadeh, Z. Pinar, H. M. Baskonus, S. Sarwar, G. Yel, "Novel explicit solutions for the nonlinear Zoomeron equation by using newly extended direct algebraic technique," Optical and Quantum Electronics, vol. 52 no. 52, 2020.

[58] W. Gao, H. F. Ismael, H. Bulut, H. M. Baskonus, "Instability modulation for the ( 2 + 1 )-dimension paraxial wave equation and its new optical soliton solutions in Kerr media," Physica Scripta, vol. 95 no. 35207, 2020.

[59] W. Gao, B. Ghanbari, H. Günerhan, H. M. Baskonus, "Some mixed trigonometric complex soliton solutions to the perturbed nonlinear Schrödinger equation," Modern Physics Letters B, vol. 34 no. 3, article 2050034,DOI: 10.1142/s0217984920500347, 2020.

[60] Y. Yakup, Y. Emrullah, "Multiple exp-function method for soliton solutions of nonlinear evolution equations," Chinese Physics B, vol. 26 no. 7, pp. 20-26, 2017.

[61] A. G. Alnowehy, "The multiple exp-function method and the linear superposition principle for solving the ( 2 + 1 )-dimensional Calogero-Bogoyavlenskii-Schiff equation," Zeitschrift für Naturforschung A, vol. 70 no. 9, pp. 775-779, DOI: 10.1515/zna-2015-0151, 2015.

[62] R. A. Abdullahi, "The generalized ( 1 + 1 )-dimensional and ( 2 + 1 )-dimensional Ito equations: multiple exp-function algorithm and multiple wave solutions," Computers & Mathematcs with Applications, vol. 71, pp. 1248-1258, 2016.

[63] Y. Long, Y. He, S. Li, "Multiple soliton solutions for a new generalization of the associated Camassa-Holm equation by exp-function method," Mathematical Problems in Engineering, vol. 2014,DOI: 10.1155/2014/418793, 2014.

[64] W. X. Ma, Z. Zhu, "Solving the ( 3 + 1 )-dimensional generalized KP and BKP equations by the multiple exp-function algorithm," Applied Mathematics and Computation, vol. 218 no. 24, pp. 11871-11879, DOI: 10.1016/j.amc.2012.05.049, 2012.

[65] J. G. Liu, L. Zhou, Y. He, "Multiple soliton solutions for the new ( 2 + 1 )-dimensional Korteweg-de Vries equation by multiple exp-function method," Applied Mathematics Letters, vol. 80, pp. 71-78, DOI: 10.1016/j.aml.2018.01.010, 2018.