# Application of the Exp-Function and Generalized Kudryashov Methods for Obtaining New Exact Solutions of Certain Nonlinear Conformable Time Partial Integro-Differential Equations

^{1}

^{2}

^{*}

Previous Article in Journal

Previous Article in Special Issue

Previous Article in Special Issue

Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand

Centre of Excellence in Mathematics, CHE, Si Ayutthaya Road, Bangkok 10400, Thailand

Author to whom correspondence should be addressed.

Academic Editors: Yongwimon Lenbury, Ravi P. Agarwal, Philip Broadbridge, Dongwoo Sheen and Shan Zhao

Received: 25 March 2021 / Revised: 16 April 2021 / Accepted: 23 April 2021 / Published: 26 April 2021

(This article belongs to the Special Issue Proceedings of the International Conference in Mathematics and Applications 2020_Mahidol University)

The key objective of this paper is to construct exact traveling wave solutions of the conformable time second integro-differential Kadomtsev–Petviashvili (KP) hierarchy equation using the Exp-function method and the (2 + 1)-dimensional conformable time partial integro-differential Jaulent–Miodek (JM) evolution equation utilizing the generalized Kudryashov method. These two problems involve the conformable partial derivative with respect to time. Initially, the conformable time partial integro-differential equations can be converted into nonlinear ordinary differential equations via a fractional complex transformation. The resulting equations are then analytically solved via the corresponding methods. As a result, the explicit exact solutions for these two equations can be expressed in terms of exponential functions. Setting some specific parameter values and varying values of the fractional order in the equations, their 3D, 2D, and contour solutions are graphically shown and physically characterized as, for instance, a bell-shaped solitary wave solution, a kink-type solution, and a singular multiple-soliton solution. To the best of the authors’ knowledge, the results of the equations obtained using the proposed methods are novel and reported here for the first time. The methods are simple, very powerful, and reliable for solving other nonlinear conformable time partial integro-differential equations arising in many applications.

The study of solutions of nonlinear partial differential equations (NPDEs) attracts the attention of scientists because their solutions can be used to lucidly explain many physical phenomena in various scientific fields, such as fluid mechanics, quantum mechanics, plasma physics, biology, chemistry, fiber optics, and many other branches of engineering. Obtaining solutions for NPDEs is of great significance for analyzing and better understanding the behaviors of the considered problems. There are many robust, stable, and effective methods that have been developed for constructing exact, approximate analytical and numerical solutions for NPDEs. Particularly, the methods have been extensively used to find exact solutions for NPDEs, such as the $(G{}^{\prime}/G,1/G)$-expansion method [1], the modified $\left({G}^{\prime}/{G}^{2}\right)$-expansion method [2], the Jacobi elliptic equation method [3], the sine-Gordon expansion method [4,5], the extended direct algebraic method [6,7], the modified Exp-function method [8,9], and the Kudryashov method [10,11]. However, the focus of this work is to search for exact solutions of certain nonlinear partial integro-differential equations (PIDEs) converted into NPDEs in some ways. Consequently, in the following, we give a brief review of PIDEs and methods through which to find their exact solutions.

A partial integro-differential equation (PIDE) is a mathematical equation involving partial derivatives and integrals of an unknown function of two or more independent variables. In the recent times, partial integro-differential equations (PIDEs) have been of considerable importance because they are widely used to model real-world problems and describe several physical phenomena in engineering, finance, and other fields of science. The applications of PIDEs have been studied in many papers. Sachs and Strauss [12] proposed a new method to efficiently solve Merton’s model written in terms of a PIDE explaining the option price for general Lévy processes. The initial and boundary value problem expressed by the nonlinear weakly singular partial integro-differential equation arising from viscoelasticity was proposed and analyzed using a Legendre wavelet collocation method (LWCM) by Singh et al. [13]. Khaled [14] employed the sinc-Galerkin method to obtain numerical solutions for a parabolic Volterra integro-differential equation presenting the heat transfer of heterogeneous materials. To further grasp other applications of PIDEs, one can refer to [15,16,17,18,19]. As mentioned above, PIDEs are used to represent various real-world problems. Therefore, the investigation of finding solutions to those PIDEs plays an important role in describing mechanisms and physical behaviors of the unknown variables in such problems. Some recent examples of the methods that have been employed to obtain numerical solutions and approximate analytical solutions for PIDEs can be found in [13,14,20,21]. However, the behaviors of the systems formulated using PIDEs can be precisely described via their exact solutions. In consequence, researchers have taken an interest in studying efficient and reliable techniques for constructing exact solutions of PIDEs. Some recent reviews for finding exact solution of the interesting PIDEs are as follows: In 2018, the three shallow water wave models formulated by PIDEs were explored using the improved ansatz method [22]. As a result, their new exact solutions, such as topological soliton solutions, were obtained. In 2019, Liu and Wen [23] constructed N-soliton solutions for the (2 + 1)-dimensional asymmetrical Nizhnik–Novikov–Veselov equation using Hirota’s bilinear method. In 2020, Seadawy et al. [24] used the generalized direct algebraic extended simple equation and modified F-expansion methods to construct some new traveling wave solutions including rational hyperbolic function and trigonometric function solutions for the (1 + 1)-dimensional Ito integro-differential equation. Thus far, there have been large number of studies that focus on computing exact solutions for PIDEs using various powerful techniques [25,26,27,28,29,30,31].

In 2014, Khalil et al. [32] introduced a new interesting derivative definition called the conformable derivative that extends the classical derivative by inserting a fractional order. The advantages of the conformable derivative, beyond classical fractional derivatives such as the Riemann–Liouville fractional order derivative or the Caputo fractional order derivative [33], are as follows: First, the conformable derivative definition satisfies most of the properties that the classical integral derivative has, such as linearity, product rule, quotient rule, power rule, chain rule, vanishing derivatives for constant functions, Rolle’s theorem, the mean value theorem, and Taylor’s theorem. Furthermore, the contribution of this derivative is applied to self-adjointness, the Sturm–Liouville system, and integral transforms such as the Fourier and Laplace transforms. Second, solving differential equations, which are used to model physical problems, with the conformable derivative is easier, as their solutions are simpler than those associated with the classical fractional derivatives. Finally, many researchers have applied the conformable derivative to various real applications. For example, the conformable derivative is employed to develop the Swartzendruber model for description of non-Darcian flow in porous media [34]. The applications of the conformable derivative in quantum mechanics and perturbation theory and their physical interpretation are discussed in [35]. Further details of its mathematical properties and real applications can be found in [35,36,37,38,39,40,41] and the references cited therein.

In this article, we are interested in using the Exp-function method [42,43,44,45] and the generalized Kudryashov method [46,47] to solve the conformable time second integro-differential Kadomtsev–Petviashvili hierarchy equation and the (2 + 1)-dimensional conformable time partial integro-differential Jaulent–Miodek evolution equation, respectively, finding their exact solutions. These equations are the PIDEs in the sense of the conformable partial derivative with respect to time. To the best of the authors’ knowledge, no researchers have found exact solutions for such PIDEs. The mentioned problems are explicitly expressed as follows:

1. The conformable time second integro-differential Kadomtsev–Petviashvili (KP) hierarchy equation of order $0<\alpha \le 1$ is expressed as
where ${\partial}_{t}^{\alpha}\left(\xb7\right)=\frac{{\partial}^{\alpha}}{\partial {t}^{\alpha}}\left(\xb7\right)$ is the conformable partial derivative with respect to t of order $\alpha $ and ${\partial}_{x}^{-n}(\xb7)={\int}_{-\infty}^{x}{\int}_{-\infty}^{{x}_{1}}\cdots {\int}_{-\infty}^{{x}_{n-1}}(\xb7)\phantom{\rule{3.33333pt}{0ex}}d{x}_{n}\cdots d{x}_{2}\phantom{\rule{3.33333pt}{0ex}}d{x}_{1}$ for some positive integer n. When $\alpha =1$, Equation (1) reduces into the second integro-differential Kadomtsev–Petviashvili hierarchy equation [48,49] which is a PIDE in two spatial coordinates and one temporal coordinate and describes the evolution of nonlinear long waves of small amplitude with slow dependence on the transverse coordinate [49].

$$\begin{array}{ccc}\hfill {\partial}_{t}^{\alpha}u& =& \frac{1}{16}{u}_{xxxxx}+\frac{5}{4}{\partial}_{x}^{-1}\left(u{u}_{yy}\right)+\frac{5}{4}{\partial}_{x}^{-1}\left({u}_{y}^{2}\right)+\frac{5}{16}{\partial}_{x}^{-3}\left({u}_{yyyy}\right)+\frac{5}{4}{u}_{x}{\partial}_{x}^{-2}\left({u}_{yy}\right)\hfill \\ & & +\frac{5}{2}u{\partial}_{x}^{-1}\left({u}_{yy}\right)+\frac{5}{2}{u}_{y}{\partial}_{x}^{-1}\left({u}_{y}\right)+\frac{15}{2}{u}^{2}{u}_{x}+\frac{5}{2}{u}_{x}{u}_{xx}+\frac{5}{4}u{u}_{xxx}+\frac{5}{8}{u}_{xyy},\hfill \end{array}$$

2. The (2 + 1)-dimensional conformable time partial integro-differential Jaulent–Miodek (JM) evolution equation of order $0<\alpha \le 1$ is given as
where the symbols ${\partial}_{t}^{\alpha}\left(\xb7\right)$ and ${\partial}_{x}^{-1}\left(\xb7\right)$ can be defined using the notations described above and where ${p}_{1},\phantom{\rule{0.166667em}{0ex}}{p}_{2},\phantom{\rule{0.166667em}{0ex}}{p}_{3},\phantom{\rule{0.166667em}{0ex}}{p}_{4},\phantom{\rule{0.166667em}{0ex}}{p}_{5}$ are arbitrary constants. If $\alpha =1$, then Equation (2) turns out to be the (2 + 1)-dimensional JM equation [50,51], elucidating many branches of physics such as condensed matter physics, fluid dynamics, and optics [16]. In particular, the (2 + 1)-dimensional JM equation associates with energy-dependent Schrödinger potential [18,52,53,54].

$$\begin{array}{c}\hfill {p}_{1}{\partial}_{t}^{\alpha}u+{p}_{2}{u}^{2}{u}_{x}-{u}_{xxx}-{p}_{3}{u}_{x}{\partial}_{x}^{-1}\left({u}_{y}\right)-{p}_{4}u{u}_{y}+{p}_{5}{\partial}_{x}^{-1}\left({u}_{yy}\right)=0,\end{array}$$

The remaining parts of this article are organized as follows: In Section 2, a definition of the conformable derivative and its properties are presented. In Section 3, we provide the main steps of the Exp-function method and the generalized Kudryashov method. The application of each of the methods for the proposed equations is exhibited in Section 4. Graphs and physical explanations of some selected exact solutions to the problems are demonstrated in Section 5. A discussion and conclusions of the results are given in the last section.

In this section, we provide a definition of the conformable derivative and its important properties as established by Khalil et al. [32].

Given a function $f:[0,\infty )\to \mathbb{R}$. Then, the conformable derivative of f of order α is defined as [32,37,39,41,55,56,57,58]
for all $t>0$ and $0<\alpha \le 1.$ If the limit in Equation (3) exists, then we can say that f is α-conformable differentiable at a point $t>0$. Furthermore, if f is α-conformable differentiable in some $(0,a),\phantom{\rule{0.166667em}{0ex}}a>0$ and ${lim}_{t\to {0}^{+}}{D}_{t}^{\alpha}f\left(t\right)$ exists, then one defines ${D}_{t}^{\alpha}f\left(0\right)={lim}_{t\to {0}^{+}}{D}_{t}^{\alpha}f\left(t\right).$

$$\begin{array}{c}\hfill \phantom{\rule{3.33333pt}{0ex}}{D}_{t}^{\alpha}f\left(t\right)=\underset{\epsilon \to 0}{lim}{\displaystyle \frac{f(t+\epsilon {t}^{1-\alpha})-f\left(t\right)}{\epsilon}},\end{array}$$

Let $\alpha \in (0,1]$, and $f\left(t\right),\phantom{\rule{3.33333pt}{0ex}}g\left(t\right)$ be $\alpha $-conformable differentiable functions at a point $t>0$. Then, the properties of the conformable derivative are as follows [32,43,55,56,57,58]:

- (1)
- ${D}_{t}^{\alpha}\left(\lambda \right)=0,\phantom{\rule{3.33333pt}{0ex}}\mathrm{where}\phantom{\rule{3.33333pt}{0ex}}\lambda =\mathrm{constant};$
- (2)
- ${D}_{t}^{\alpha}\left({t}^{\mu}\right)=\mu {t}^{\mu -\alpha},$ for all $\mu \in \mathbb{R};$
- (3)
- ${D}_{t}^{\alpha}(af\left(t\right)+bg\left(t\right))=a{D}_{t}^{\alpha}f\left(t\right)+b{D}_{t}^{\alpha}g\left(t\right),$ for all $a,b\in \mathbb{R};$
- (4)
- ${D}_{t}^{\alpha}\left(f\left(t\right)g\left(t\right)\right)=f\left(t\right){D}_{t}^{\alpha}g\left(t\right)+g\left(t\right){D}_{t}^{\alpha}f\left(t\right);$
- (5)
- ${D}_{t}^{\alpha}\left({\displaystyle \frac{f\left(t\right)}{g\left(t\right)}}\right)={\displaystyle \frac{g\left(t\right){D}_{t}^{\alpha}f\left(t\right)-f\left(t\right){D}_{t}^{\alpha}g\left(t\right)}{g{\left(t\right)}^{2}}};$
- (6)
- If, in addition, f is differentiable, then ${D}_{t}^{\alpha}\left(f\left(t\right)\right)={t}^{1-\alpha}{\displaystyle \frac{df\left(t\right)}{dt}}$.

Using the definition in (3) and the above properties, the conformable derivatives of some interesting functions are defined as follows [32,55,56,57,58]:

- (1)
- ${D}_{t}^{\alpha}\left({e}^{ct}\right)=c{t}^{1-\alpha}{e}^{ct},\phantom{\rule{3.33333pt}{0ex}}c\in \mathbb{R}$;
- (2)
- ${D}_{t}^{\alpha}(sinbt)=b{t}^{1-\alpha}cosbt,b\in \mathbb{R}$;
- (3)
- ${D}_{t}^{\alpha}(cosbt)=-b{t}^{1-\alpha}sinbt,b\in \mathbb{R}$;
- (4)
- ${D}_{t}^{\alpha}\left(\frac{1}{\alpha}{t}^{\alpha}\right)=1$.

([32,43,55,56,57,58]). Let $f:{\mathbb{R}}^{+}\to \mathbb{R}$ be a function such that f is differentiable and α-conformable differentiable. In addition, we assume that g is a differentiable function defined in the range of f. Then, we have
where the prime symbol ${(}^{\prime})$ denotes the classical derivative.

$${D}_{t}^{\alpha}(f\circ g)\left(t\right)={t}^{1-\alpha}{f}^{\prime}\left(g\left(t\right)\right){g}^{\prime}\left(t\right),$$

The Exp-function method and the generalized Kudryashov method described in this section are applied to the proposed PIDEs in the next section. Now, we consider the following general nonlinear conformable partial differential equation:
where ${\partial}_{v}^{\gamma}u=\frac{{\partial}^{\gamma}}{\partial {v}^{\gamma}}u$ is the generic term for the conformable partial derivative of a dependent variable u with respect to an independent variable v of order $\gamma $. Thus, the symbols, for instance, ${\partial}_{t}^{\alpha}u=\frac{{\partial}^{\alpha}}{\partial {t}^{\alpha}}u,\phantom{\rule{0.166667em}{0ex}}{\partial}_{x}^{\beta}u=\frac{{\partial}^{\beta}}{\partial {x}^{\beta}}u$ and ${\partial}_{y}^{\eta}u=\frac{{\partial}^{\eta}}{\partial {y}^{\eta}}u$, that appear in Equation (4) denote the conformable partial derivatives of the dependent variable $u=u(x,y,t)$ with respect to t of order $\alpha $, to x of order $\beta $, and to y of order $\eta $, respectively, with $0<\alpha ,\beta ,\eta \le 1$. The remaining partial derivatives including the mixed partial derivatives in Equation (4) can be defined in a similar manner to the previous ones. The function F in the equation is a polynomial of the unknown function u and its various conformable partial derivatives. Using the definition in (3), we can define the conformable partial derivative ${\partial}_{t}^{\alpha}u(x,y,t)$ as

$$\begin{array}{c}\hfill F(u,{\partial}_{t}^{\alpha}u,{\partial}_{x}^{\beta}u,{\partial}_{y}^{\eta}u,{\partial}_{t}^{2\alpha}u,{\partial}_{t}^{\alpha}{\partial}_{x}^{\beta}u,{\partial}_{t}^{\alpha}{\partial}_{y}^{\eta}u,{\partial}_{x}^{2\beta}u,{\partial}_{x}^{\beta}{\partial}_{y}^{\eta}u,...)=0,\end{array}$$

$$\begin{array}{c}\hfill {\partial}_{t}^{\alpha}u(x,y,t)={\displaystyle \underset{\epsilon \to 0}{lim}}\frac{u(x,y,t+\epsilon {t}^{1-\alpha})-u(x,y,t)}{\epsilon},\phantom{\rule{3.33333pt}{0ex}}t>0.\end{array}$$

The symbols ${\partial}_{x}^{\beta}u,\phantom{\rule{0.166667em}{0ex}}{\partial}_{y}^{\eta}u$ and the other partial derivatives in (4) can be defined in a manner analogous to Equation (5).

The common step of the two methods is to transform the conformable PDE in (4) into an ordinary differential equation (ODE) using the following transformation [59,60,61]:
where $\xi $ is a traveling wave variable with a nonzero arbitrary constant k. Applying the above transformation to Equation (4) and then integrating the resulting equation with respect to $\xi $ (if possible), Equation (4) can be converted into the ODE in the variable $U=U\left(\xi \right)$ as
where P is a polynomial of $U\left(\xi \right)$ and its various integer-order derivatives. The prime notation (${}^{\prime}$) represents an ordinary derivative with respect to $\xi $. Next, the remaining steps of the the Exp-function method [62,63,64] and the generalized Kudryashov method [46,54,57,65] are separately and compactly explained as follows.

$$\begin{array}{c}\hfill U\left(\xi \right)=u(x,y,t),\phantom{\rule{3.33333pt}{0ex}}\phantom{\rule{3.33333pt}{0ex}}\xi =x+y-{\displaystyle \frac{k{t}^{\alpha}}{\alpha}},\end{array}$$

$$\begin{array}{c}\hfill P(U,{U}^{\prime},{U}^{\u2033},{U}^{\u2034},...)=0,\end{array}$$

$$\begin{array}{c}\hfill U\left(\xi \right)=\frac{{\sum}_{n=-c}^{d}{a}_{n}exp\left(n\xi \right)}{{\sum}_{m=-p}^{q}{b}_{m}exp\left(m\xi \right)},\end{array}$$

$$\begin{array}{c}\hfill U\left(\xi \right)=\frac{{a}_{-c}exp(-c\xi )+...+{a}_{d}exp\left(d\xi \right)}{{b}_{-p}exp(-p\xi )+\dots +{b}_{q}exp\left(q\xi \right)}.\end{array}$$

$$\begin{array}{c}\hfill U\left(\xi \right)=\frac{{\sum}_{i=0}^{N}{a}_{i}{Q}^{i}\left(\xi \right)}{{\sum}_{j=0}^{M}{b}_{j}{Q}^{j}\left(\xi \right)},\end{array}$$

$$\begin{array}{c}\hfill {Q}^{\prime}\left(\xi \right)={Q}^{2}\left(\xi \right)-Q\left(\xi \right).\end{array}$$

It is not difficult to ascertain that
is the solution for Equation (11) when $\kappa $ is a constant of integration.

$$\begin{array}{c}\hfill Q\left(\xi \right)=\frac{1}{1+\kappa {e}^{\xi}}\end{array}$$

In this section, we use the Exp-function method and the generalized Kudryashov method to obtain explicit exact solutions for Equations (1) and (2), respectively. Before utilizing the methods as described above, we must convert the conformable PIDEs in (1) and (2) to conformable PDEs and then apply the methods to the resulting equations. These are demonstrated in the following.

Equation (1), which contains at most threefold integral operators, can be converted into a PDE with the conformable time partial derivative of a new variable by differentiating both sides of (1) with respect to x and applying the transformation $u(x,y,t)={v}_{xx}(x,y,t)$ to the resulting equation. We consequently obtain

$$\begin{array}{ccc}\hfill \frac{{\partial}^{\alpha}}{\partial {t}^{\alpha}}{v}_{xxx}& =& \frac{1}{16}{v}_{8x}+\frac{15}{4}{v}_{xx}{v}_{xxyy}+\frac{15}{4}{v}_{xxy}^{2}+\frac{5}{16}{v}_{4y}+\frac{15}{4}{v}_{xxx}{v}_{xyy}+\frac{5}{4}{v}_{xxxx}{v}_{yy}\hfill \\ & & +\frac{5}{2}{v}_{xxxy}{v}_{xy}+\frac{15}{2}{\left({v}_{xx}^{2}{v}_{xxx}\right)}_{x}+\frac{5}{2}{\left({v}_{xxx}{v}_{xxxx}\right)}_{x}+\frac{5}{4}{\left({v}_{xx}{v}_{xxxxx}\right)}_{x}+\frac{5}{8}{v}_{4xyy}.\hfill \end{array}$$

Applying the following transformation:
where $k\ne 0$ is a constant to (13), one can subsequently obtain an ordinary differential equation in the new variable $V=V\left(\xi \right)$. Then, we integrate the resulting equation with respect to $\xi $ by setting a constant of integration to zero. Hence, we have
where the prime notation ${(}^{\prime})$ represents the classical derivative with respect to $\xi $. Letting $W\left(\xi \right)={V}^{\u2033}\left(\xi \right)$ in Equation (15) and then integrating the resulting outcome with respect to $\xi $, we now obtain
where C is a constant of integration.

$$\begin{array}{c}\hfill V\left(\xi \right)=v(x,y,t),\phantom{\rule{3.33333pt}{0ex}}\phantom{\rule{3.33333pt}{0ex}}\xi =x+y-{\displaystyle \frac{k{t}^{\alpha}}{\alpha}},\end{array}$$

$$\begin{array}{ccc}\hfill & & \frac{1}{16}{V}^{\left(7\right)}+\frac{5}{8}{V}^{\left(5\right)}+\left(k+\frac{5}{16}\right){V}^{\u2034}+\frac{15}{2}{V}^{\u2033}{V}^{\u2034}+\frac{15}{2}{\left({V}^{\u2033}\right)}^{2}{V}^{\u2034}+\frac{5}{2}{V}^{\u2034}{V}^{\left(4\right)}\hfill \\ & & +\frac{5}{4}{V}^{\u2033}{V}^{\left(5\right)}=0,\hfill \end{array}$$

$$\begin{array}{c}\hfill \frac{1}{16}{W}^{\left(4\right)}+\left(k+\frac{5}{16}\right)W+\frac{15}{4}{W}^{2}+\frac{5}{2}{W}^{3}+\frac{5}{8}{\left({W}^{\prime}\right)}^{2}+\frac{5}{4}W{W}^{\u2033}+\frac{5}{8}{W}^{\u2033}+C=0,\end{array}$$

Applying the Exp-function method to Equation (16), we assume that the solution for Equation (16) can be expressed in the following form:
where the positive integers $c,\phantom{\rule{0.166667em}{0ex}}p,\phantom{\rule{0.166667em}{0ex}}d,\phantom{\rule{0.166667em}{0ex}}q$ can be determined using the homogeneous balancing principle as follows: In order to determine the values of d and q in (17), we balance the linear term of the highest order derivative in Equation (16) with the highest order nonlinear term. In particular, the terms ${W}^{\left(4\right)}$ and ${W}^{3}$ in (16) are balanced. From (17), we obtain
and
where ${c}_{i},\phantom{\rule{0.166667em}{0ex}}i=1,2,3,4$ are constant coefficients, which are not necessarily computed for their exact values because they are not required in the further steps. On balancing the highest order of the exponential functions in Equations (18) and (19), we get

$$\begin{array}{c}\hfill W\left(\xi \right)=\frac{{a}_{-c}exp(-c\xi )+\dots +{a}_{d}exp\left(d\xi \right)}{{b}_{-p}exp(-p\xi )+\dots +{b}_{q}exp\left(q\xi \right)},\end{array}$$

$${W}^{\left(4\right)}=\frac{{c}_{1}exp\left[(15q+d)\xi \right]+\dots}{{c}_{2}exp\left(16q\xi \right)+\dots},$$

$${W}^{3}=\frac{{c}_{3}exp\left(3d\xi \right)+\dots}{{c}_{4}exp\left(3q\xi \right)+\dots}=\frac{{c}_{3}exp\left[(3d+13q)\xi \right]+\dots}{{c}_{4}exp\left(16q\xi \right)+\dots},$$

$$15q+d=13q+3d\Rightarrow q=d.$$

Proceeding in a similar manner as illustrated above, we can determine the values of c and p. Using (17) and after manipulating certain calculations, we can write
and
where ${d}_{i},\phantom{\rule{0.166667em}{0ex}}i=1,2,3,4$ are constant coefficients, the exact values of which again are not required. Balancing the lowest order of the exponential functions in Equations (21) and (22), we have

$${W}^{\left(4\right)}=\frac{\dots +{d}_{1}exp[-(15p+c)\xi ]}{\dots +{d}_{2}exp(-16q\xi )},$$

$${W}^{3}=\frac{\dots +{d}_{3}exp(-3c\xi )}{\dots +{d}_{4}exp(-3p\xi )}=\frac{\dots +{d}_{3}exp[-(13p+3c)\xi ]}{\dots +{d}_{4}exp(-16p\xi )},$$

$$-(15p+c)=-(13p+3c)\Rightarrow c=p.$$

For simplicity, we choose $p=c=d=q=1$. Hence, the solution form (17) becomes
where ${a}_{-1},\phantom{\rule{3.33333pt}{0ex}}{a}_{0},\phantom{\rule{3.33333pt}{0ex}}{a}_{1},\phantom{\rule{3.33333pt}{0ex}}{b}_{-1},\phantom{\rule{3.33333pt}{0ex}}{b}_{0},$, and ${b}_{1}$ are unknown constants that are determined in a later step.

$$\begin{array}{c}\hfill W\left(\xi \right)=\frac{{a}_{1}exp\left(\xi \right)+{a}_{0}+{a}_{-1}exp(-\xi )}{{b}_{1}exp\left(\xi \right)+{b}_{0}+{b}_{-1}exp(-\xi )},\end{array}$$

Substituting Equation (24) into Equation (16) and simplifying the resulting equation with the aid of Maple, we have
where $A={\left({b}_{1}exp\left(\xi \right)+{b}_{0}+{b}_{-1}exp(-\xi )\right)}^{5}$. Since each of the coefficients ${C}_{i}$$(i=-5,-4,$ $\dots ,4,5)$ has a very long expression, all of them are transferred to Appendix A. Equating all of the coefficients ${C}_{i}$$(i=-5,-4,\dots ,4,5)$ to be zero, we obtain a system of algebraic equations as follows:

$$\begin{array}{cc}\hfill \phantom{\rule{1.em}{0ex}}& \frac{1}{A}\left[{C}_{5}exp\left(5\xi \right)+{C}_{4}exp\left(4\xi \right)+{C}_{3}exp\left(3\xi \right)+{C}_{2}exp\left(2\xi \right)+{C}_{1}exp\left(\xi \right)+{C}_{0}+{C}_{-1}exp(-\xi )\right.\hfill \\ \hfill \phantom{\rule{1.em}{0ex}}& \left.+{C}_{-2}exp(-2\xi )+{C}_{-3}exp(-3\xi )+{C}_{-4}exp(-4\xi )+{C}_{-5}exp(-5\xi )\right]=0,\hfill \end{array}$$

$$\left\{\begin{array}{ccccc}{C}_{-5}=0,& {C}_{-4}=0,& {C}_{-3}=0,& {C}_{-2}=0,& {C}_{-1}=0,\\ {C}_{0}=0,\\ {C}_{1}=0,& {C}_{2}=0,& {C}_{3}=0,& {C}_{4}=0,& {C}_{5}=0.\end{array}\right.$$

With the help of Maple, we can simultaneously solve system (26) for ${a}_{-1},\phantom{\rule{3.33333pt}{0ex}}{a}_{0},\phantom{\rule{3.33333pt}{0ex}}{a}_{1},\phantom{\rule{3.33333pt}{0ex}}{b}_{-1},\phantom{\rule{3.33333pt}{0ex}}{b}_{0},{b}_{1},\phantom{\rule{3.33333pt}{0ex}}k$, and C. The only result obtained from the above process is as follows:
where ${a}_{-1},\phantom{\rule{3.33333pt}{0ex}}{b}_{0}$ and ${b}_{1}$ are free parameters. Substituting the above result into (24), we get the solution for Equation (16) as

$$\begin{array}{ccc}\hfill & & {a}_{0}=\frac{8{a}_{-1}{b}_{1}+{b}_{0}^{2}}{2{b}_{0}},\phantom{\rule{3.33333pt}{0ex}}{a}_{1}=\frac{8{a}_{-1}{b}_{1}+{b}_{0}^{2}}{2{b}_{0}},\phantom{\rule{3.33333pt}{0ex}}{b}_{-1}=\frac{{b}_{0}^{2}}{4{b}_{1}},\hfill \\ & & C=\frac{{a}_{-1}{b}_{1}\left(1280{a}_{-1}^{2}{b}_{1}^{2}+320{a}_{-1}{b}_{0}^{2}{b}_{1}+11{b}_{0}^{4}\right)}{4{b}_{0}^{6}},\phantom{\rule{3.33333pt}{0ex}}k=-\frac{120{a}_{-1}^{2}{b}_{1}^{2}+35{a}_{-1}{b}_{0}^{2}{b}_{1}+{b}_{0}^{4}}{{b}_{0}^{4}},\hfill \end{array}$$

$$\begin{array}{ccc}\hfill W\left(\xi \right)& =& \frac{{b}_{1}\left(16{a}_{-1}{b}_{1}^{2}{\mathrm{e}}^{2\xi}+\left(16{a}_{-1}{b}_{0}{b}_{1}+2{b}_{0}^{3}\right){\mathrm{e}}^{\xi}+4{a}_{-1}{b}_{0}^{2}\right)}{{b}_{0}^{2}\left(4{b}_{1}^{2}{\mathrm{e}}^{2\xi}+4{b}_{0}{b}_{1}{\mathrm{e}}^{\xi}+{b}_{0}^{2}\right)}.\hfill \end{array}$$

Since $W={V}^{\u2033}$, we then have

$$\begin{array}{c}\hfill V\left(\xi \right)=\int \int W\left(\xi \right)\phantom{\rule{3.33333pt}{0ex}}d\xi \phantom{\rule{3.33333pt}{0ex}}d\xi =\frac{2{b}_{1}{a}_{-1}{\xi}^{2}+{b}_{0}^{2}ln\left({b}_{0}+2{b}_{1}{\mathrm{e}}^{\xi}\right)-{b}_{0}^{2}\xi}{{b}_{0}^{2}}.\end{array}$$

Inserting $\xi =x+y-\frac{k{t}^{\alpha}}{\alpha}$ into the above result to obtain $v(x,y,t)$ and then taking the partial derivative of v with respect to x twice, we obtain the exact solution $u(x,y,t)$ of (1) as
where $\theta (x,y,t)=\frac{\left(120{a}_{-1}^{2}{b}_{1}^{2}+35{a}_{-1}{b}_{0}^{2}{b}_{1}+{b}_{0}^{4}\right){t}^{\alpha}+\alpha {b}_{0}^{4}\left(x+y\right)}{\alpha {b}_{0}^{4}}$.

$$\begin{array}{ccc}\hfill u(x,y,t)& =& \frac{2{b}_{1}\left(8{a}_{-1}{b}_{1}^{2}{\mathrm{e}}^{2\theta (x,y,t)}+\left(8{a}_{-1}{b}_{0}{b}_{1}+{b}_{0}^{3}\right){\mathrm{e}}^{\theta (x,y,t)}+2{a}_{-1}{b}_{0}^{2}\right)}{{b}_{0}^{2}{\left(2{b}_{1}{\mathrm{e}}^{\theta (x,y,t)}+{b}_{0}\right)}^{2}},\hfill \end{array}$$

Using the transformation $u(x,y,t)={v}_{x}(x,y,t)$, Equation (2) can be converted into the following conformable PDE:

$$\begin{array}{c}\hfill {p}_{1}\left(\frac{{\partial}^{\alpha}}{\partial {t}^{\alpha}}{v}_{x}\right)+{p}_{2}{v}_{x}^{2}{v}_{xx}-{v}_{xxxx}-{p}_{3}{v}_{xx}{v}_{y}-{p}_{4}{v}_{x}{v}_{xy}+{p}_{5}{v}_{yy}=0.\end{array}$$

Applying the traveling wave transform (14) to the above equation and integrating the resulting equation once, we then obtain the following ODE in $V\left(\xi \right)$ as
where k is a constant of the traveling wave transform and C is a constant of integration. These constants are determined in a further step. Letting $W={V}^{\prime}$, Equation (32) becomes

$$\begin{array}{c}\hfill {V}^{\u2034}-({p}_{5}-k{p}_{1}){V}^{\prime}+\frac{({p}_{3}+{p}_{4})}{2}{\left({V}^{\prime}\right)}^{2}-\frac{{p}_{2}}{3}{\left({V}^{\prime}\right)}^{3}+C=0,\end{array}$$

$$\begin{array}{c}\hfill {W}^{\u2033}-({p}_{5}+k{p}_{1})W+\frac{({p}_{3}+{p}_{4})}{2}{W}^{2}-\frac{{p}_{2}}{3}{W}^{3}+C=0.\end{array}$$

On the basis of Equation (10) of the generalized Kudryashov method, we assume that the solution form for Equation (33) is $W\left(\xi \right)=\frac{{\sum}_{i=0}^{N}{a}_{i}{Q}^{i}\left(\xi \right)}{{\sum}_{j=0}^{M}{b}_{j}{Q}^{j}\left(\xi \right)}$. Then, $\mathrm{Deg}\left[W\right]=N-M$. Balancing the highest order derivative ${W}^{\u2033}$ in Equation (33) with the nonlinear term ${W}^{3}$ via formulas (10) of [57], we get

$$\begin{array}{c}\hfill N-M+2=3(N-M)\Rightarrow N=M+1.\end{array}$$

Choosing $M=1$, we obtain $N=2$. Hence, the solution for (33) takes the following particular form:
where ${a}_{0},\phantom{\rule{3.33333pt}{0ex}}{a}_{1},\phantom{\rule{3.33333pt}{0ex}}{a}_{2},\phantom{\rule{3.33333pt}{0ex}}{b}_{0}$, and ${b}_{1}$ are constants to be determined at a later step and where $Q=Q\left(\xi \right)$ expressed in (12) satisfies Equation (11). Inserting solution (35) along with Equation (11) into Equation (33), we obtain a polynomial in $Q\left(\xi \right)$. Equating all the coefficients of like power of Q in the polynomial to be zero, the system of the algebraic equations expressed in terms of ${a}_{0},\phantom{\rule{3.33333pt}{0ex}}{a}_{1},\phantom{\rule{3.33333pt}{0ex}}{a}_{2},\phantom{\rule{3.33333pt}{0ex}}{b}_{0},\phantom{\rule{3.33333pt}{0ex}}{b}_{1},\phantom{\rule{3.33333pt}{0ex}}C$ and k is obtained, as shown below:

$$\begin{array}{c}\hfill W\left(\xi \right)=\frac{{a}_{0}+{a}_{1}Q+{a}_{2}{Q}^{2}}{{b}_{0}+{b}_{1}Q},\end{array}$$

$$\begin{array}{cc}\hfill {Q}^{0}:& {a}_{0}{b}_{0}^{2}{p}_{1}-C{{b}_{0}}^{3}-{a}_{0}{b}_{0}^{2}{p}_{5}+\frac{{a}_{0}^{2}{b}_{0}{p}_{3}}{2}+\frac{{a}_{0}^{2}{b}_{0}{p}_{4}}{2}-\frac{{p}_{2}{a}_{0}^{3}}{3}=0,\hfill \\ \hfill {Q}^{1}:& -{a}_{0}^{2}{p}_{2}{a}_{1}+\frac{{a}_{0}^{2}{b}_{1}{p}_{3}}{2}+\frac{{a}_{0}^{2}{b}_{1}{p}_{4}}{2}-3C{{b}_{0}}^{2}{b}_{1}-{a}_{0}{b}_{0}{b}_{1}-{a}_{1}{b}_{0}^{2}{p}_{5}+2k{a}_{0}{b}_{0}{b}_{1}{p}_{1}+{b}_{0}^{2}{a}_{1}\hfill \\ \hfill \phantom{\rule{1.em}{0ex}}& +{a}_{0}{a}_{1}{b}_{0}{p}_{3}+{a}_{0}{a}_{1}{b}_{0}{p}_{4}-2{a}_{0}{b}_{0}{b}_{1}{p}_{5}+k{a}_{1}{b}_{0}^{2}{p}_{1}=0,\hfill \\ \hfill {Q}^{2}:& 2{b}_{1}^{2}{a}_{2}-\frac{{p}_{2}{a}_{2}^{3}}{3}=0,\hfill \\ \hfill {Q}^{3}:& {\frac{{a}_{2}}{2}}^{2}{b}_{1}{p}_{4}-{p}_{2}{a}_{1}{a}_{2}^{2}+6{a}_{2}{b}_{0}{b}_{1}+\frac{{a}_{2}^{2}{b}_{1}{p}_{3}}{2}-3{b}_{1}^{2}{a}_{2}=0,\hfill \\ \hfill \phantom{\rule{4pt}{0ex}}{Q}^{4}:& \frac{{a}_{2}^{2}{b}_{0}{p}_{4}}{2}-{p}_{2}{a}_{0}{{a}_{2}}^{2}-{p}_{2}{{a}_{1}}^{2}{a}_{2}-9\phantom{\rule{0.166667em}{0ex}}{a}_{2}{b}_{0}{b}_{1}-{a}_{2}{{b}_{1}}^{2}{p}_{5}+\frac{{{a}_{2}}^{2}{b}_{0}{p}_{3}}{2}+6{{b}_{0}}^{2}{a}_{2}+{a}_{1}{a}_{2}{b}_{1}{p}_{3}\hfill \\ \hfill \phantom{\rule{1.em}{0ex}}& +{a}_{1}{a}_{2}{b}_{1}{p}_{4}+k{a}_{2}{{b}_{1}}^{2}{p}_{1}+{{b}_{1}}^{2}{a}_{2}=0,\hfill \\ \hfill {Q}^{5}:& -2{a}_{0}{b}_{0}{b}_{1}+{a}_{1}{b}_{0}{b}_{1}+3{a}_{2}{b}_{0}{b}_{1}-{a}_{1}{b}_{1}^{2}{p}_{5}+\frac{{b}_{1}{a}_{1}^{2}{p}_{3}}{2}+\frac{{b}_{1}{a}_{1}^{2}{p}_{4}}{2}+2k{a}_{2}{b}_{0}{b}_{1}{p}_{1}+2{b}_{0}^{2}{a}_{1}\hfill \\ \hfill \phantom{\rule{1.em}{0ex}}& -10{b}_{0}^{2}{a}_{2}-C{b}_{1}^{3}+{b}_{1}{a}_{0}{a}_{2}{p}_{3}+{a}_{1}{a}_{2}{b}_{0}{p}_{3}+{b}_{1}{a}_{0}{a}_{2}{p}_{4}+{a}_{1}{a}_{2}{b}_{0}{p}_{4}-2{a}_{2}{b}_{0}{b}_{1}{p}_{5}+k{a}_{1}{b}_{1}^{2}{p}_{1}\hfill \\ \hfill \phantom{\rule{1.em}{0ex}}& -\frac{{p}_{2}{a}_{1}^{3}}{3}-{a}_{0}{b}_{1}^{2}-2{p}_{2}{a}_{0}{a}_{1}{a}_{2}=0,\hfill \\ \hfill {Q}^{6}:& \frac{{a}_{1}^{2}{b}_{0}{p}_{4}}{2}-{a}_{0}^{2}{a}_{2}{p}_{2}-{a}_{0}{a}_{1}^{2}{p}_{2}-{a}_{0}{b}_{1}^{2}{p}_{5}-{a}_{1}{b}_{0}{b}_{1}-3C{b}_{0}{b}_{1}^{2}+3{a}_{0}{b}_{0}{b}_{1}-{a}_{2}{b}_{0}^{2}{p}_{5}+\frac{{a}_{1}^{2}{b}_{0}{p}_{3}}{2}\hfill \\ \hfill \phantom{\rule{1.em}{0ex}}& +2k{a}_{1}{b}_{0}{b}_{1}{p}_{1}-3{b}_{0}^{2}{a}_{1}+4{b}_{0}^{2}{a}_{2}+{a}_{0}{a}_{1}{b}_{1}{p}_{3}+{a}_{0}{a}_{2}{b}_{0}{p}_{3}+{a}_{0}{a}_{1}{b}_{1}{p}_{4}+{a}_{0}{a}_{2}{b}_{0}{p}_{4}+k{a}_{0}{b}_{1}^{2}{p}_{1}\hfill \\ \hfill \phantom{\rule{1.em}{0ex}}& -2{a}_{1}{b}_{0}{b}_{1}{p}_{5}+k{a}_{2}{b}_{0}^{2}{p}_{1}+{a}_{0}{b}_{1}^{2}=0.\hfill \end{array}$$

Solving the above system with the help of Maple, we obtain two different results for the exact solutions for Equation (2) as follows:

Result 1:

$$\begin{array}{ccc}\hfill & & {a}_{0}=\frac{{b}_{0}\left(\mp \sqrt{6{p}_{2}}+{p}_{3}+{p}_{4}\right)}{2{p}_{2}},\phantom{\rule{3.33333pt}{0ex}}{a}_{1}=\frac{\pm \sqrt{6{p}_{2}}\left(2{b}_{0}-{b}_{1}\right)+{b}_{1}\left({p}_{3}+{p}_{4}\right)}{2{p}_{2}},\hfill \\ & & {a}_{2}=\pm {b}_{1}\sqrt{\frac{6}{{p}_{2}}},\phantom{\rule{3.33333pt}{0ex}}C=\frac{\left(-{p}_{3}^{2}-2{p}_{3}{p}_{4}-{p}_{4}^{2}+6{p}_{2}\right)\left({p}_{3}+{p}_{4}\right)}{24{p}_{2}^{2}},\hfill \\ & & k=\frac{4{p}_{2}{p}_{5}-{p}_{3}^{2}-2{p}_{3}{p}_{4}-{p}_{4}^{2}+2{p}_{2}}{4{p}_{1}{p}_{2}},\hfill \end{array}$$

where ${b}_{0}$ and ${b}_{1}$ are arbitrary constants. As a result of substituting the above results into solution (35) with the replacement of $Q\left(\xi \right)$ in (12) and then simplifying the resulting outcome, the solution for Equation (33) is obtained:

$$\begin{array}{ccc}\hfill {W}_{1}\left(\xi \right)& =& \frac{\left(\sqrt{6}{p}_{2}\left(\kappa {\mathrm{e}}^{\xi}-1\right)\pm \sqrt{{p}_{2}}\left({p}_{3}+{p}_{4}\right)\left(1+\kappa {\mathrm{e}}^{\xi}\right)\right)}{2{p}_{2}^{3/2}\left(\kappa {\mathrm{e}}^{\xi}+1\right)}.\hfill \end{array}$$

As a result of ${W}_{1}={V}_{1}^{\prime}$, we then obtain

$$\begin{array}{ccc}\hfill {V}_{1}\left(\xi \right)& =& \int {W}_{1}\left(\xi \right)\phantom{\rule{3.33333pt}{0ex}}d\xi ,\hfill \\ & =& \frac{\left(\sqrt{{p}_{2}}\left({p}_{3}+{p}_{4}\right)\pm \sqrt{6}{p}_{2}\right)\xi \mp 2\sqrt{6}{p}_{2}ln\left(1+\kappa {\mathrm{e}}^{\xi}\right)}{2{p}_{2}^{3/2}}.\hfill \end{array}$$

By replacing $\xi =x+y-{\displaystyle \frac{k{t}^{\alpha}}{\alpha}}$ in the above result to get ${v}_{1}(x,y,t)$ and then taking the partial derivative of ${v}_{1}$ with respect to x, the exact solution of (2) is obtained:
where $\omega (x,y,t)=\frac{\left({\left({p}_{3}+{p}_{4}\right)}^{2}-\left(4{p}_{5}+2\right){p}_{2}\right){t}^{\alpha}+4\alpha {p}_{1}{p}_{2}\left(x+y\right)}{4\alpha {p}_{1}{p}_{2}}$.

$$\begin{array}{ccc}\hfill {u}_{1}(x,y,t)& =& \frac{\sqrt{{p}_{2}}\left({p}_{3}+{p}_{4}\right)\left(\kappa {\mathrm{e}}^{\omega (x,y,t)}+1\right)\mp \sqrt{6}{p}_{2}\left(\kappa {\mathrm{e}}^{\omega (x,y,t)}-1\right)}{{p}_{2}^{3/2}\left(2\kappa {\mathrm{e}}^{\omega (x,y,t)}+2\right)},\hfill \end{array}$$

Result 2:

$$\begin{array}{ccc}\hfill & & {a}_{0}=\frac{\pm {b}_{0}\left(2{b}_{0}\sqrt{6{p}_{2}}+{b}_{1}\left(\sqrt{6{p}_{2}}+{p}_{3}+{p}_{4}\right)\right)}{2{p}_{2}{b}_{1}},\hfill \\ & & {a}_{1}=\pm \frac{2{b}_{0}\sqrt{6{p}_{2}}+{b}_{1}\left({p}_{3}+{p}_{4}-\sqrt{6{p}_{2}}\right)}{2{p}_{2}},\hfill \\ & & {a}_{2}=\pm {b}_{1}\sqrt{\frac{6}{{p}_{2}}},\phantom{\rule{3.33333pt}{0ex}}C=\frac{1}{24{p}_{2}^{2}{b}_{1}^{3}}(\pm 96\sqrt{6}{b}_{0}{p}_{2}^{3/2}\left({b}_{0}+\frac{{b}_{1}}{2}\right)\left({b}_{0}+{b}_{1}\right)\hfill \\ & & +72{b}_{1}\left({p}_{3}+{p}_{4}\right)\left({b}_{1}^{2}\left(-\frac{{p}_{3}^{2}}{72}-\frac{{p}_{4}{p}_{3}}{36}-\frac{{p}_{4}^{2}}{72}+\frac{{p}_{2}}{12}\right)+{b}_{0}{b}_{1}{p}_{2}+{b}_{0}^{2}{p}_{2}\right)),\hfill \\ & & k=\frac{{b}_{1}^{2}\left({p}_{2}\left(4{p}_{5}+2\right)-{\left({p}_{3}+{p}_{4}\right)}^{2}\right)+24{p}_{2}{b}_{0}\left({b}_{1}+{b}_{0}\right)}{4{b}_{1}^{2}{p}_{1}{p}_{2}},\hfill \end{array}$$

where ${b}_{0}$ and ${b}_{1}$ are arbitrary constants. Substituting these values into solution (35), using the function $Q\left(\xi \right)$ in (12), and then simplifying the resulting equation, the solution of Equation (33) is expressed as

$$\begin{array}{ccc}\hfill {W}_{2}\left(\xi \right)& =& \frac{\pm 1}{12{p}_{2}^{3/2}{b}_{1}\left({b}_{0}\kappa {\mathrm{e}}^{\xi}+{b}_{0}+{b}_{1}\right)\left(1+\kappa {\mathrm{e}}^{\xi}\right)}[6{b}_{0}{\kappa}^{2}{\mathrm{e}}^{2\xi}\left({b}_{1}\sqrt{{p}_{2}}\left({p}_{3}+{p}_{4}\right)\right.\hfill \\ & & \left.+\sqrt{6}{p}_{2}\left(2{b}_{0}+{b}_{1}\right)\right)+{b}_{1}\sqrt{{p}_{2}}\left(\kappa {\mathrm{e}}^{\xi}\left(12{b}_{0}+6{b}_{1}\right)+6{b}_{0}+6{b}_{1}\right)\left({p}_{3}+{p}_{4}\right)\hfill \\ & & +24\sqrt{6}{p}_{2}\left(\kappa {\mathrm{e}}^{\xi}\left({b}_{0}^{2}+{b}_{0}{b}_{1}-\frac{{b}_{1}^{2}}{4}\right)+\frac{1}{2}\left({b}_{0}+\frac{{b}_{1}}{2}\right)\left({b}_{0}+{b}_{1}\right)\right)].\hfill \end{array}$$

Since ${W}_{2}={V}_{2}^{\prime}$, we have

$$\begin{array}{ccc}\hfill {V}_{2}\left(\xi \right)& =& \int {W}_{2}\left(\xi \right)\phantom{\rule{3.33333pt}{0ex}}d\xi ,\hfill \\ & =& \frac{\pm 1}{2{p}_{2}^{3/2}{b}_{1}}[2\sqrt{6}{b}_{1}{p}_{2}ln\left({b}_{0}\kappa {\mathrm{e}}^{\xi}+{b}_{0}+{b}_{1}\right)-2\sqrt{6}{b}_{1}{p}_{2}ln\left(1+\kappa {\mathrm{e}}^{\xi}\right)\hfill \\ & & +\xi \left({b}_{1}\sqrt{{p}_{2}}\left({p}_{3}+{p}_{4}\right)+2\sqrt{6}{p}_{2}\left({b}_{0}+\frac{{b}_{1}}{2}\right)\right)].\hfill \end{array}$$

Inserting $\xi =x+y-{\displaystyle \frac{k{t}^{\alpha}}{\alpha}}$ into Equation (42) to get ${v}_{2}(x,y,t)$ and then the taking partial derivative of ${v}_{2}$ with respect to x, the other exact solution for (2) can be written as
where $\delta (x,y,t)=\frac{\left(\left({\left({p}_{3}+{p}_{4}\right)}^{2}-\left(4{p}_{5}+2\right){p}_{2}\right){b}_{1}^{2}-24{b}_{0}{b}_{1}{p}_{2}-24{b}_{0}^{2}{p}_{2}\right){t}^{\alpha}+4\alpha {b}_{1}^{2}{p}_{1}{p}_{2}\left(x+y\right)}{4\alpha {b}_{1}^{2}{p}_{1}{p}_{2}}.$

$$\begin{array}{ccc}\hfill {u}_{2}(x,y,t)& =& \frac{\pm 1}{8\alpha {p}_{2}^{5/2}{p}_{1}{b}_{1}^{3}}[\frac{8\sqrt{6}\alpha \kappa {b}_{0}{b}_{1}^{3}{p}_{1}{p}_{2}^{2}{\mathrm{e}}^{\delta (x,y,t)}}{\kappa {b}_{0}{\mathrm{e}}^{\delta (x,y,t)}+{b}_{0}+{b}_{1}}-\frac{8\sqrt{6}\alpha \kappa {b}_{1}^{3}{p}_{1}{p}_{2}^{2}{\mathrm{e}}^{\delta (x,y,t)}}{\kappa {\mathrm{e}}^{\delta (x,y,t)}+1}\hfill \\ & & +4\alpha {{b}_{1}}^{3}{p}_{1}{p}_{2}^{3/2}\left({p}_{3}+{p}_{4}\right)+8\sqrt{6}\alpha {b}_{1}^{2}{p}_{1}{p}_{2}^{2}\left({b}_{0}+\frac{{b}_{1}}{2}\right)],\hfill \end{array}$$

In this section, we provide interesting graphical representations of the exact solutions of the conformable time second integro-differential KP hierarchy Equation (1) obtained using the Exp-function method and of the (2 + 1)-dimensional conformable time partial integro-differential Jaulent–Miodek evolution Equation (2) obtained utilizing the generalized Kudryashov method. The time-fractional order $\alpha $ for the equations is varied in order to explore graphical behaviors of the exact solutions reported in the previous section. In particular, the values of the time-fractional order used for all simulations are $\alpha =1,\phantom{\rule{3.33333pt}{0ex}}0.8$ and $0.2$. The exact solution (30) for Equation (1) and the exact solutions (39) and (43) for Equation (2) are graphically portrayed as 3D, 2D, and contour plots according to the used values of $\alpha $. All of the 3D plots are drawn on the domain $0\le x,t\le 10$ with the fixed value $y=1$. The 2D graphs, showing a relation of $u\left(x\right)$ and x, are depicted on $0\le x\le 10$ with $y=t=1.$ Here, the contour plots, presenting a 3D surface by plotting constant u slices on a 2D plane, are plotted to connect the $(x,t)$ coordinates, where the given values of u occur and y is fixed at $y=1$. Furthermore, the physical explanations of the displayed graphs are discussed in this section.

In Figure 1, the graphs of the exact solution $u(x,y,t)$ in (30), constructed using the Exp-function method, for problem (1) are unfolded in different aspects. Using the values of $\alpha $ as shown above and setting the free parameters ${a}_{-1}=-1,\phantom{\rule{3.33333pt}{0ex}}{b}_{0}=5$ and ${b}_{1}=0.001$, the solution (30) is evaluated and plotted on the domains for the 3D, 2D, and contour plots in Figure 1. In particular, Figure 1a–c shows the 3D, 2D, and contour plots for solution (30), respectively, when $\alpha =1$. Figure 1d–i are portrayed in the same manner as before except using $\alpha =0.8$ and $\alpha =0.2,$ respectively. By classifying the shapes of the 3D and 2D graphs in Figure 1, one can identify that solution (30) is a bell-shaped solitary wave solution.

Figure 2 and Figure 3 show the solution graphs for problem (2), which are generated using the generalized Kudryashov method. Specifically, the solution ${u}_{1}(x,y,t)$ in (39), in which the top sign of ∓ is chosen, is graphically plotted in Figure 2, and the solution ${u}_{2}(x,y,t)$ in (43), in which the top sign of ± is taken, is drawn in Figure 3. Apart from using $\alpha =1,\phantom{\rule{3.33333pt}{0ex}}0.8,\phantom{\rule{3.33333pt}{0ex}}0.2$ and the domains mentioned above, the following two sets of the parameter values $\{\kappa =1,\phantom{\rule{0.166667em}{0ex}}{p}_{1}=-3,\phantom{\rule{0.166667em}{0ex}}{p}_{2}=3,\phantom{\rule{0.166667em}{0ex}}{p}_{3}=3,\phantom{\rule{0.166667em}{0ex}}{p}_{4}=5,\phantom{\rule{0.166667em}{0ex}}{p}_{5}=1\}$ and $\{\kappa =1,\phantom{\rule{0.166667em}{0ex}}{p}_{1}=3,\phantom{\rule{0.166667em}{0ex}}{p}_{2}=1,\phantom{\rule{0.166667em}{0ex}}{p}_{3}=-1,\phantom{\rule{0.166667em}{0ex}}{p}_{4}=-0.1,\phantom{\rule{0.166667em}{0ex}}{p}_{5}=3,\phantom{\rule{0.166667em}{0ex}}{b}_{0}=-0.05,{b}_{1}=0.6\}$ are inserted into solution (39) and solution (43), respectively, in order to depict their 3D, 2D, and contour plots. Figure 2a–c demonstrates the 3D, 2D, and contour graphs of the solution ${u}_{1}(x,y,t)$ when $\alpha =1$ is used. However, the 3D, 2D, and contour plots of the solution ${u}_{1}(x,y,t)$ when $\alpha =0.8$ and $\alpha =0.2$ are shown in Figure 2d–i. The physical behavior of these graphs is characterized as a kink-type solution.

Furthermore, Figure 3a–c displays the 3D, 2D, and contour plots of the solution ${u}_{2}(x,y,t)$ in (43) when $\alpha =1$. Proceeding in a similar manner to the above plots, except utilizing $\alpha =0.8$ and $\alpha =0.2$, the corresponding graphs are portrayed in Figure 3d–i. As can be seen in the graph structures shown in Figure 3, their physical behavior is considered as a singular multiple-soliton solution.

We can observe some effects of the variation of the time-fractional order $\alpha $ on the 3D and 2D graphs in Figure 1, Figure 2 and Figure 3. It can be possibly concluded that when the value of $\alpha $ is varied to depict the 3D and 2D solutions in each figure, their graphical structures remain the same, but their magnitudes change and translations are detected. In addition, how the magnitude and translation of the 3D and 2D graphs alter depends on the change in $\alpha $.

In summary, the Exp-function method, the generalized Kudryashov method, the use of transformation (6), and symbolic software packages such as Maple were successfully employed to obtain exact traveling wave solutions for the conformable time second integro-differential KP hierarchy Equation (1) and the (2 + 1)-dimensional conformable time partial integro-differential JM evolution Equation (2), respectively. Utilizing the Exp-function method to analytically solve Equation (1) and removing the trivial and disqualified solutions, only one exact solution (30) was obtained. As a result, solution (30) is expressed in terms of the exponential functions, or equivalently, the hyperbolic functions, and it behaves as a bell-shaped solitary wave solution under a certain set of parameter values. In order to compare our solution (30) with certain solutions from the relevant literature, one can consider the mathematical expressions between our solution (when $\alpha =1$) and the solutions in [28,30]. In [28], the generalized Kudryashov method was utilized to solve the second integro-differential KP hierarchy equation for which the exponential function solutions were generated. Their solutions displayed the behavior of a bell-shaped solution. In [30], one of the solutions of the second integro-differential KP hierarchy equation, established by the $({G}^{\prime}/G,1/G)$-expansion method, was expressed in terms of the hyperbolic functions and characterized as an anti-bell soliton solution. Generally speaking, solution (30) when $\alpha =1$ has a mathematical structure similar to that of the solutions in the aforementioned literature.

Applying the generalized Kudryashov method to obtain exact solutions for Equation (2) and after eliminating the trivial solutions, we obtained two different explicit exact solutions, as shown in Equations (39) and (43). The solutions are written in terms of the exponential functions. By plotting these solutions computed using the specific sets of parameter values, a kink-type solution and a singular multiple-soliton solution were obtained. A comparison between our solutions for (2) when $\alpha =1$ and the exact solutions of the (2 + 1)-dimensional Jaulent–Miodek equation obtained in [51] is as follows: In [51], the complex method was used to solve the (2 + 1)-dimensional JM equation for which the hyperbolic cotangent function solutions were established. The solutions were characterized as a solitary wave solution of singular kink-type. It can be seen that the structures of their solution and our solution are found in many physical phenomena, such optics and nonlinear waves.

All of the obtained exact solutions discussed in this paper were verified by substituting them back into their corresponding equations with the help of the Maple package program. The Exp-function method and the generalized Kudryashov method are straightforward, reliable, efficient, and pragmatic mathematical tools for solving the proposed equations because they produce uncomplicated exact solution forms. These two methods could be effectively applied to solve a wide range of nonlinear partial integro-differential equations arising in natural phenomena, giving their analytically extracted exact solutions.

Conceptualization, S.S. (Sekson Sirisubtawee) and S.S. (Surattana Sungnul); methodology, S.K. and S.S. (Sekson Sirisubtawee); software, S.K. and S.S. (Sekson Sirisubtawee); validation, S.K., S.S. (Sekson Sirisubtawee) and S.S. (Surattana Sungnul); formal analysis, S.K. and S.S. (Sekson Sirisubtawee); investigation, S.K. and S.S. (Sekson Sirisubtawee); resources, S.S. (Sekson Sirisubtawee) and S.S. (Surattana Sungnul); data curation, S.K.; writing—original draft preparation, S.K. and S.S. (Sekson Sirisubtawee); writing—review and editing, S.S. (Sekson Sirisubtawee) and S.S. (Surattana Sungnul); visualization, S.K., S.S. (Sekson Sirisubtawee) and S.S. (Surattana Sungnul); supervision, S.S. (Sekson Sirisubtawee) and S.S. (Surattana Sungnul); project administration, S.S. (Sekson Sirisubtawee) and S.S. (Surattana Sungnul); funding acquisition, S.K. and S.S. (Sekson Sirisubtawee). All authors have read and agreed to the published version of the manuscript.

The first author was funded by King Mongkut’s University of Technology North Bangkok under contract no. KMUTNB-61-PHD-014. The second and third authors were financially supported by the Centre of Excellence in Mathematics, the Commission on Higher Education, Thailand.

Not applicable.

Not applicable.

Not applicable.

The authors are grateful to anonymous referees for the valuable comments which have significantly improved this article. In addition, the first author would like to acknowledge the partial support from the Graduate College, King Mongkut’s University of Technology North Bangkok.

The authors declare no conflict of interest.

This part demonstrates the expression of each coefficient ${C}_{i}$ of $exp\left(i\xi \right),$$i=-5,-4,$$\dots ,4,5$ in Equation (25) as follows:

$$\begin{array}{ccc}\hfill {C}_{-5}& =& k{a}_{-1}{b}_{-1}^{4}+\frac{5{a}_{-1}^{3}{b}_{-1}^{2}}{2}+\frac{5{a}_{-1}{b}_{-1}^{4}}{16}+C{b}_{-1}^{5}+\frac{15{a}_{-1}^{2}{b}_{-1}^{3}}{4},\hfill \\ \hfill {C}_{-4}& =& k{a}_{0}{b}_{-1}^{4}+\frac{35{a}_{0}{a}_{-1}{b}_{-1}^{3}}{4}+\frac{15{a}_{0}{a}_{-1}^{2}{b}_{-1}^{2}}{2}+10{a}_{-1}^{2}{b}_{-1}^{2}{b}_{0}+5C{b}_{0}{b}_{-1}^{4}+5{a}_{-1}^{3}{b}_{0}{b}_{-1}\hfill \\ & & +\frac{9{a}_{-1}{b}_{-1}^{3}{b}_{0}}{16}+{a}_{0}{b}_{-1}^{4}+4k{a}_{-1}{b}_{0}{b}_{-1}^{3},\hfill \\ \hfill {C}_{-3}& =& 10C{b}_{0}^{2}{b}_{-1}^{3}+5{a}_{-1}^{3}{b}_{1}{b}_{-1}+\frac{25{a}_{1}{a}_{-1}{b}_{-1}^{3}}{2}+k{a}_{1}{b}_{-1}^{4}-\frac{94{a}_{-1}{b}_{1}{b}_{-1}^{3}}{4}+\frac{31{a}_{-1}{b}_{-1}^{2}{b}_{0}^{2}}{16}\hfill \\ & & +\frac{19{a}_{0}{b}_{-1}^{3}{b}_{0}}{16}+\frac{15{a}_{1}{a}_{-1}^{2}{b}_{-1}^{2}}{2}+\frac{15{a}_{0}^{2}{a}_{-1}{b}_{-1}^{2}}{2}+\frac{25{a}_{-1}^{2}{b}_{1}{b}_{-1}^{2}}{4}+\frac{95{a}_{-1}^{2}{b}_{-1}{b}_{0}^{2}}{8}+5C{b}_{1}{b}_{-1}^{4}\hfill \\ & & +\frac{5{a}_{-1}^{3}{b}_{0}^{2}}{2}+\frac{45{a}_{0}^{2}{b}_{-1}^{3}}{8}+\frac{61{a}_{1}{b}_{-1}^{4}}{16}+4k{a}_{0}{b}_{0}{b}_{-1}^{3}+15{a}_{0}{a}_{-1}^{2}{b}_{0}{b}_{-1}+4k{a}_{-1}{b}_{1}{b}_{-1}^{3}\hfill \\ & & +6k{a}_{-1}{b}_{0}^{2}{b}_{-1}^{2}+20{a}_{0}{a}_{-1}{b}_{0}{b}_{-1}^{2},\hfill \\ \hfill {C}_{-2}& =& \frac{19{a}_{-1}{b}_{-1}{b}_{0}^{3}}{16}+\frac{129{a}_{1}{b}_{-1}^{3}{b}_{0}}{16}+\frac{31{a}_{0}{b}_{-1}^{2}{b}_{0}^{2}}{16}+\frac{65{a}_{1}{a}_{0}{b}_{-1}^{3}}{4}+\frac{15{a}_{0}{{a}_{-1}}^{2}{b}_{0}^{2}}{2}+\frac{95{a}_{0}^{2}{b}_{-1}^{2}{b}_{0}}{8}\hfill \\ & & +10C{b}_{0}^{3}{b}_{-1}^{2}+5{a}_{-1}^{3}{b}_{1}{b}_{0}-6{a}_{0}{b}_{1}{b}_{-1}^{3}+\frac{55{a}_{1}{a}_{-1}{b}_{0}{b}_{-1}^{2}}{2}+4k{a}_{0}{b}_{1}{b}_{-1}^{3}+6k{a}_{0}{b}_{0}^{2}{b}_{-1}^{2}\hfill \\ & & +15{a}_{1}{a}_{0}{a}_{-1}{b}_{-1}^{2}+15{a}_{1}{a}_{-1}^{2}{b}_{0}{b}_{-1}+15{a}_{0}^{2}{a}_{-1}{b}_{0}{b}_{-1}+15{a}_{0}{a}_{-1}^{2}{b}_{1}{b}_{-1}+\frac{45{a}_{-1}^{2}{b}_{1}{b}_{0}{b}_{-1}}{2}\hfill \\ & & +20C{b}_{1}{b}_{0}{b}_{-1}^{3}+4k{a}_{1}{b}_{0}{b}_{-1}^{3}+4k{a}_{-1}{b}_{0}^{3}{b}_{-1}+\frac{67{a}_{-1}{b}_{-1}^{2}{b}_{1}{b}_{0}}{16}+\frac{35{a}_{0}{a}_{-1}{b}_{1}{b}_{-1}^{2}}{4}\hfill \\ & & +20{a}_{0}{a}_{-1}{b}_{0}^{2}{b}_{-1}+12k{a}_{-1}{b}_{1}{b}_{0}{b}_{-1}^{2}+\frac{5{a}_{0}^{3}{b}_{-1}^{2}}{2}+\frac{45{{a}_{-1}}^{2}{b}_{0}^{3}}{8},\hfill \\ \hfill {C}_{-1}& =& \frac{25}{2}{a}_{1}{a}_{-1}{b}_{1}{b}_{-1}^{2}+\frac{95{a}_{1}{a}_{-1}{b}_{0}^{2}{b}_{-1}}{4}+12k{a}_{0}{b}_{1}{b}_{0}{b}_{-1}^{2}+12k{a}_{-1}{b}_{1}{b}_{0}^{2}{b}_{-1}+30{a}_{0}{a}_{-1}{b}_{1}{b}_{0}{b}_{-1}\hfill \\ & & +30{a}_{1}{a}_{0}{a}_{-1}{b}_{0}{b}_{-1}+4k{a}_{0}{b}_{0}^{3}{b}_{-1}+15{a}_{1}{a}_{-1}^{2}{b}_{1}{b}_{-1}+15{a}_{0}^{2}{a}_{-1}{b}_{1}{b}_{-1}+15{a}_{0}{a}_{-1}^{2}{b}_{1}{b}_{0}\hfill \\ & & +30C{b}_{1}{b}_{0}^{2}{b}_{-1}^{2}+4k{a}_{1}{b}_{1}{b}_{-1}^{3}+6k{a}_{1}{b}_{0}^{2}{b}_{-1}^{2}+6k{a}_{-1}{b}_{1}^{2}{b}_{-1}^{2}+\frac{11{a}_{-1}{b}_{-1}{b}_{1}{b}_{0}^{2}}{8}\hfill \\ & & -\frac{23{a}_{0}{b}_{1}{b}_{-1}^{2}{b}_{0}}{16}+35{a}_{1}{a}_{0}{b}_{0}{b}_{-1}^{2}+\frac{15{a}_{1}^{2}{a}_{-1}{b}_{-1}^{2}}{2}+\frac{35{a}_{0}{a}_{-1}{b}_{0}^{3}}{4}+\frac{151{a}_{1}{b}_{-1}^{2}{b}_{0}^{2}}{16}\hfill \\ & & +\frac{15{a}_{1}{a}_{0}^{2}{b}_{-1}^{2}}{2}+\frac{15{a}_{0}^{2}{a}_{-1}{b}_{0}^{2}}{2}+\frac{15{a}_{1}{a}_{-1}^{2}{b}_{0}^{2}}{2}+\frac{9{a}_{0}{b}_{-1}{b}_{0}^{3}}{16}+\frac{55{a}_{-1}^{2}{b}_{1}^{2}{b}_{-1}}{4}+\frac{155{a}_{-1}^{2}{b}_{1}{b}_{0}^{2}}{8}\hfill \\ & & +\frac{35{a}_{0}^{2}{b}_{-1}^{2}{b}_{1}}{8}+10{a}_{0}^{2}{b}_{-1}{b}_{0}^{2}+10C{b}_{1}^{2}{b}_{-1}^{3}+5C{{b}_{0}}^{4}{b}_{-1}+5{a}_{0}^{3}{b}_{0}{b}_{-1}+k{a}_{-1}{b}_{0}^{4}\hfill \\ & & +\frac{83{a}_{-1}{b}_{1}^{2}{b}_{-1}^{2}}{8}-\frac{29{a}_{1}{b}_{1}{b}_{-1}^{3}}{4}+\frac{45{a}_{1}^{2}{b}_{-1}^{3}}{4}+\frac{5{a}_{-1}^{3}{b}_{1}^{2}}{2}+{a}_{-1}{b}_{0}^{4},\hfill \end{array}$$

$$\begin{array}{ccc}\hfill {C}_{0}& =& -\frac{25{a}_{0}{b}_{1}{b}_{-1}{b}_{0}^{2}}{8}+\frac{15{a}_{1}^{2}{a}_{0}{b}_{-1}^{2}}{2}+15{a}_{1}{a}_{0}{a}_{-1}{b}_{0}^{2}+\frac{15{a}_{0}{a}_{-1}^{2}{b}_{1}^{2}}{2}+25{a}_{1}^{2}{b}_{-1}^{2}{b}_{0}+25{a}_{-1}^{2}{b}_{1}^{2}{b}_{0}\hfill \\ & & +30C{b}_{1}^{2}{b}_{0}{b}_{-1}^{2}+20C{b}_{1}{b}_{0}^{3}{b}_{-1}+\frac{75{a}_{1}{b}_{0}^{3}{b}_{-1}}{16}+10{a}_{0}{b}_{1}^{2}{b}_{-1}^{2}+12k{a}_{1}{b}_{1}{b}_{0}{b}_{-1}^{2}\hfill \\ & & +12k{a}_{-1}{b}_{1}^{2}{b}_{0}{b}_{-1}+6k{a}_{0}{b}_{1}^{2}{b}_{-1}^{2}+k{a}_{0}{b}_{0}^{4}+\frac{35{a}_{1}{a}_{-1}{b}_{0}^{3}}{4}+\frac{75{a}_{-1}{b}_{1}{b}_{0}^{3}}{16}+5{a}_{0}^{3}{b}_{1}{b}_{-1}\hfill \\ & & +C{b}_{0}^{5}+\frac{15{a}_{0}^{2}{b}_{0}^{3}}{4}+\frac{5{a}_{0}^{3}{b}_{0}^{2}}{2}+\frac{5{a}_{0}{b}_{0}^{4}}{16}+30{a}_{1}{a}_{-1}{b}_{1}{b}_{0}{b}_{-1}+12k{a}_{0}{b}_{1}{b}_{0}^{2}{b}_{-1}+15{a}_{1}^{2}{a}_{-1}{b}_{0}{b}_{-1}\hfill \\ & & +15{a}_{1}{a}_{0}^{2}{b}_{0}{b}_{-1}+30{a}_{1}{a}_{0}{a}_{-1}{b}_{1}{b}_{-1}+15{a}_{1}{a}_{-1}^{2}{b}_{1}{b}_{0}+15{a}_{0}^{2}{a}_{-1}{b}_{1}{b}_{0}+\frac{45{a}_{0}^{2}{b}_{1}{b}_{0}{b}_{-1}}{4}\hfill \\ & & +4k{a}_{1}{b}_{0}^{3}{b}_{-1}+4k{a}_{-1}{b}_{1}{b}_{0}^{3}-\frac{5{a}_{1}{b}_{1}{b}_{0}{b}_{-1}^{2}}{16}-\frac{5{a}_{-1}{b}_{-1}{b}_{1}^{2}{b}_{0}}{16}+\frac{65{a}_{1}{a}_{0}{b}_{1}{b}_{-1}^{2}}{4}\hfill \\ & & +\frac{55{a}_{1}{a}_{0}{b}_{0}^{2}{b}_{-1}}{2}+\frac{65{a}_{0}{a}_{-1}{b}_{1}^{2}{b}_{-1}}{4}+\frac{55{a}_{0}{a}_{-1}{b}_{1}{b}_{0}^{2}}{2},\hfill \\ \hfill {C}_{1}& =& 30C{b}_{1}^{2}{b}_{0}^{2}{b}_{-1}+4k{a}_{-1}{b}_{1}^{3}{b}_{-1}+6k{a}_{-1}{b}_{1}^{2}{b}_{0}^{2}+\frac{11{a}_{1}{b}_{1}{b}_{0}^{2}{b}_{-1}}{8}-\frac{23{a}_{0}{b}_{1}^{2}{b}_{0}{b}_{-1}}{16}+35{a}_{0}{a}_{-1}{b}_{1}^{2}{b}_{0}\hfill \\ & & +6k{a}_{1}{b}_{1}^{2}{b}_{-1}^{2}+\frac{25{a}_{1}{a}_{-1}{b}_{1}^{2}{b}_{-1}}{2}+\frac{95{a}_{1}{a}_{-1}{b}_{1}{b}_{0}^{2}}{4}+4k{a}_{0}{b}_{1}{b}_{0}^{3}+15{a}_{1}^{2}{a}_{0}{b}_{0}{b}_{-1}+15{a}_{1}^{2}{a}_{-1}{b}_{1}{b}_{-1}\hfill \\ & & +15{a}_{1}{a}_{0}^{2}{b}_{1}{b}_{-1}+30{a}_{1}{a}_{0}{b}_{1}{b}_{0}{b}_{-1}+30{a}_{1}{a}_{0}{a}_{-1}{b}_{1}{b}_{0}+12k{a}_{0}{b}_{1}^{2}{b}_{0}{b}_{-1}+12k{a}_{1}{b}_{1}{b}_{0}^{2}{b}_{-1}\hfill \\ & & +k{a}_{1}{b}_{0}^{4}+\frac{83{a}_{1}{b}_{1}^{2}{b}_{-1}^{2}}{8}-\frac{29{a}_{-1}{b}_{1}^{3}{b}_{-1}}{4}+\frac{151{a}_{-1}{b}_{1}^{2}{b}_{0}^{2}}{16}+\frac{9{a}_{0}{b}_{1}{b}_{0}^{3}}{16}+\frac{35{a}_{1}{a}_{0}{b}_{0}^{3}}{4}+\frac{15{a}_{1}^{2}{a}_{-1}{{b}_{0}}^{2}}{2}\hfill \\ & & +\frac{15{a}_{1}{{a}_{0}}^{2}{{b}_{0}}^{2}}{2}+\frac{15{a}_{1}{{a}_{-1}}^{2}{{b}_{1}}^{2}}{2}+\frac{15{a}_{0}^{2}{a}_{-1}{{b}_{1}}^{2}}{2}+\frac{55{a}_{1}^{2}{b}_{1}{b}_{-1}^{2}}{4}+\frac{155{a}_{1}^{2}{b}_{-1}{b}_{0}^{2}}{8}+\frac{35{a}_{0}^{2}{b}_{1}^{2}{b}_{-1}}{8}\hfill \\ & & +10{a}_{0}^{2}{b}_{1}{b}_{0}^{2}+5C{b}_{1}{b}_{0}^{4}+10C{b}_{1}^{3}{b}_{-1}^{2}+5{a}_{0}^{3}{b}_{1}{b}_{0}+{a}_{1}{b}_{0}^{4}+\frac{45{a}_{-1}^{2}{b}_{1}^{3}}{4}+\frac{5{a}_{1}^{3}{{b}_{-1}}^{2}}{2},\hfill \\ \hfill {C}_{2}& =& \frac{55{a}_{1}{a}_{-1}{b}_{1}^{2}{b}_{0}}{2}+12k{a}_{1}{b}_{1}^{2}{b}_{0}{b}_{-1}+4k{a}_{0}{b}_{1}^{3}{b}_{-1}+6k{a}_{0}{b}_{1}^{2}{b}_{0}^{2}+15{a}_{1}^{2}{a}_{0}{b}_{1}{b}_{-1}+15{a}_{1}^{2}{a}_{-1}{b}_{1}{b}_{0}\hfill \\ & & +15{a}_{1}{a}_{0}^{2}{b}_{1}{b}_{0}+15{a}_{1}{a}_{0}{a}_{-1}{b}_{1}^{2}+\frac{45{a}_{1}^{2}{b}_{1}{b}_{0}{b}_{-1}}{2}+20C{b}_{1}^{3}{b}_{0}{b}_{-1}+4k{a}_{1}{b}_{1}{b}_{0}^{3}+4k{a}_{-1}{b}_{1}^{3}{b}_{0}\hfill \\ & & +\frac{67{a}_{1}{b}_{1}^{2}{b}_{0}{b}_{-1}}{16}+\frac{35{a}_{1}{a}_{0}{b}_{1}^{2}{b}_{-1}}{4}+20{a}_{1}{a}_{0}{b}_{1}{b}_{0}^{2}+\frac{31{a}_{0}{b}_{1}^{2}{b}_{0}^{2}}{16}+\frac{15{a}_{1}^{2}{a}_{0}{{b}_{0}}^{2}}{2}+\frac{129{a}_{-1}{b}_{1}^{3}{b}_{0}}{16}\hfill \\ & & +\frac{65{a}_{0}{a}_{-1}{b}_{1}^{3}}{4}+\frac{95{a}_{0}^{2}{b}_{1}^{2}{b}_{0}}{8}+10C{b}_{1}^{2}{b}_{0}^{3}+5{a}_{1}^{3}{b}_{0}{b}_{-1}+\frac{19{a}_{1}{b}_{1}{b}_{0}^{3}}{16}-6{a}_{0}{b}_{1}^{3}{b}_{-1}+\frac{5{a}_{0}^{3}{{b}_{1}}^{2}}{2}\hfill \\ & & +\frac{45{a}_{1}^{2}{b}_{0}^{3}}{8},\hfill \\ \hfill {C}_{3}& =& k{a}_{-1}{b}_{1}^{4}+\frac{31{a}_{1}{b}_{1}^{2}{b}_{0}^{2}}{16}-\frac{9{a}_{1}{{b}_{1}}^{3}{b}_{-1}}{4}+{\frac{15{a}_{1}}{2}}^{2}{a}_{-1}{b}_{1}^{2}+\frac{15{a}_{1}{{a}_{0}}^{2}{{b}_{1}}^{2}}{2}+\frac{25{a}_{1}^{2}{b}_{1}^{2}{b}_{-1}}{4}+\frac{95{a}_{1}^{2}{b}_{1}{b}_{0}^{2}}{8}\hfill \\ & & +5C{b}_{1}^{4}{b}_{-1}+10C{b}_{1}^{3}{b}_{0}^{2}+5{a}_{1}^{3}{b}_{1}{b}_{-1}+\frac{25}{2}{a}_{1}{a}_{-1}{b}_{1}^{3}+20{a}_{1}{a}_{0}{b}_{1}^{2}{b}_{0}+4k{a}_{1}{b}_{1}^{3}{b}_{-1}+6k{a}_{1}{b}_{1}^{2}{b}_{0}^{2}\hfill \\ & & +4k{a}_{0}{b}_{1}^{3}{b}_{0}+15{a}_{1}^{2}{a}_{0}{b}_{1}{b}_{0}+\frac{19{a}_{0}{b}_{1}^{3}{b}_{0}}{16}+{\frac{5{a}_{1}}{2}}^{3}{b}_{0}^{2}+\frac{61{a}_{-1}{b}_{1}^{4}}{16}+\frac{45{a}_{0}^{2}{b}_{1}^{3}}{8},\hfill \\ \hfill {C}_{4}& =& 5C{b}_{1}^{4}{b}_{0}+5{a}_{1}^{3}{b}_{1}{b}_{0}+\frac{9{a}_{1}{b}_{1}^{3}{b}_{0}}{16}+k{a}_{0}{b}_{1}^{4}+\frac{35{a}_{1}{a}_{0}{b}_{1}^{3}}{4}+{\frac{15{a}_{1}}{2}}^{2}{a}_{0}{b}_{1}^{2}+10{a}_{1}^{2}{b}_{1}^{2}{b}_{0}+4k{a}_{1}{b}_{1}^{3}{b}_{0}\hfill \\ & & +{a}_{0}{b}_{1}^{4},\hfill \\ \hfill {C}_{5}& =& \frac{15{a}_{1}^{2}{b}_{1}^{3}}{4}+\frac{5{a}_{1}{b}_{1}^{4}}{16}+{\frac{5{a}_{1}}{2}}^{3}{b}_{1}^{2}+C{b}_{1}^{5}+k{a}_{1}{b}_{1}^{4}.\hfill \end{array}$$

- Ghafoor, A.; Firdous, S.; Zubair, T.; Iftikhar, A.; Zainab, S.; Mohyud-Din, S.T. (G
^{′}/G,1/G)-Expansion method for generalized ZK, Sharma-Tasso-Olver (STO) and modified ZK equations. Qscience Connect**2013**, 2013, 24. [Google Scholar] [CrossRef] - Aljahdaly, N.H. Some applications of the modified (G
^{′}/G^{2})-expansion method in mathematical physics. Results Phys.**2019**, 13, 102272. [Google Scholar] [CrossRef] - Zheng, B.; Feng, Q. The Jacobi elliptic equation method for solving fractional partial differential equations. Abstr. Appl. Anal.
**2014**, 2014, 249071. [Google Scholar] [CrossRef] - Baskonus, H.M.; Sulaiman, T.A.; Bulut, H. Novel complex and hyperbolic forms to the strain wave equation in microstructured solids. Opt. Quantum Electron.
**2018**, 50, 14. [Google Scholar] [CrossRef] - Eskitasçıoglu, E.; Aktas, M.B.; Baskonus, H.M. New complex and hyperbolic forms for Ablowitz-Kaup-Newell-Segur wave equation with fourth order. Appl. Math. Nonlinear Sci.
**2019**, 4, 105–112. [Google Scholar] [CrossRef] - Seadawy, A.R. Three-dimensional nonlinear modified Zakharov-Kuznetsov equation of ion-acoustic waves in a magnetized plasma. Comput. Math. Appl.
**2016**, 71, 201–212. [Google Scholar] [CrossRef] - Seadawy, A.R. Solitary wave solutions of two-dimensional nonlinear Kadomtsev-Petviashvili dynamic equation in dust-acoustic plasmas. Pramana
**2017**, 89, 49. [Google Scholar] [CrossRef] - Shakeel, M.; Mohyud-Din, S.T.; Iqbal, M.A. Modified extended exp-function method for a system of nonlinear partial differential equations defined by seismic sea waves. Pramana
**2018**, 91, 1–8. [Google Scholar] [CrossRef] - Yel, G.; Akturk, T. A new approach to (3 + 1)-dimensional Boiti-Leon-Manna-Pempinelli equation. Appl. Math. Nonlinear Sci.
**2020**, 5, 309–316. [Google Scholar] [CrossRef] - Ryabov, P.N.; Sinelshchikov, D.I.; Kochanov, M.B. Application of the Kudryashov method for finding exact solutions of the high order nonlinear evolution equations. Appl. Math. Comput.
**2011**, 218, 3965–3972. [Google Scholar] [CrossRef] - Kudryashov, N.A.; Ryabov, P.N.; Fedyanin, T.E.; Kutukov, A.A. Evolution of pattern formation under ion bombardment of substrate. Phys. Lett. A
**2013**, 377, 753–759. [Google Scholar] [CrossRef] - Sachs, E.; Strauss, A. Efficient solution of a partial integro-differential equation in finance. Appl. Numer. Math.
**2008**, 58, 1687–1703. [Google Scholar] [CrossRef] - Singh, S.; Patel, V.K.; Singh, V.K. Convergence rate of collocation method based on wavelet for nonlinear weakly singular partial integro-differential equation arising from viscoelasticity. Numer. Methods Partial. Differ. Eqs.
**2018**, 34, 1781–1798. [Google Scholar] [CrossRef] - Al-Khaled, K.; Darweesh, A.; Yousef, M.H. Convergence of numerical schemes for the solution of partial integro-differential equations used in heat transfer. J. Appl. Math. Comput.
**2019**, 61, 657–675. [Google Scholar] [CrossRef] - Zadeh, K.S. An integro-partial differential equation for modeling biofluids flow in fractured biomaterials. J. Theor. Biol.
**2011**, 273, 72–79. [Google Scholar] [CrossRef] - Li, Y.; Liu, X.; Xin, X. Explicit solutions, conservation laws of the extended (2 + 1)-dimensional Jaulent-Miodek equation. arXiv
**2015**, arXiv:1512.09196. [Google Scholar] - Zheng, X.; Qiu, W.; Chen, H. Three semi-implicit compact finite difference schemes for the nonlinear partial integro-differential equation arising from viscoelasticity. Int. J. Model. Simul.
**2020**. [Google Scholar] [CrossRef] - Kaewta, S.; Sirisubtawee, S.; Khansai, N. Explicit exact solutions of the (2 + 1)-dimensional integro-differential Jaulent-Miodek evolution equation using the reliable methods. Int. J. Math. Math. Sci.
**2020**, 2020, 2916395. [Google Scholar] [CrossRef] - Yapman, Ö.; Amiraliyev, G.M.; Amirali, I. Convergence analysis of fitted numerical method for a singularly perturbed nonlinear Volterra integro-differential equation with delay. J. Comput. Appl. Math.
**2019**, 355, 301–309. [Google Scholar] [CrossRef] - Hendi, F.; Al-Qarni, M. The variational Adomian decomposition method for solving nonlinear two-dimensional Volterra-Fredholm integro-differential equation. J. King Saud Univ. Sci.
**2019**, 31, 110–113. [Google Scholar] [CrossRef] - Hasseine, A.; Attarakih, M.; Belarbi, R.; Bart, H.J. On the semi-analytical solution of integro-partial differential equations. Energy Procedia
**2017**, 139, 358–366. [Google Scholar] [CrossRef] - Liu, J.; Zhang, Y.; Wang, Y. Topological soliton solutions for three shallow water waves models. Waves Random Complex Media
**2018**, 28, 508–515. [Google Scholar] [CrossRef] - Liu, Y.; Wen, X.Y. Soliton, breather, lump and their interaction solutions of the (2 + 1)-dimensional asymmetrical Nizhnik-Novikov-Veselov equation. Adv. Differ. Eqs.
**2019**, 2019, 332. [Google Scholar] [CrossRef] - Seadawy, A.; Ali, A.; Aljahdaly, N. The nonlinear integro-differential Ito dynamical equation via three modified mathematical methods and its analytical solutions. Open Phys.
**2020**, 18, 24–32. [Google Scholar] [CrossRef] - Bhrawy, A.; Alhuthali, M.S.; Abdelkawy, M. New solutions for (1 + 1)-dimensional and (2 + 1)-dimensional Ito equations. Math. Probl. Eng.
**2012**, 1–24. [Google Scholar] [CrossRef] - Lin, L.; Zhu, S.; Xu, Y.; Shi, Y. Exact solutions of Gardner equations through tanh-coth method. Appl. Math.
**2016**, 7, 2374–2381. [Google Scholar] [CrossRef] - Akbari, M.; Rasht, I. Application of Kudryashov method for the Ito equations. Appl. Math.
**2017**, 12, 136–142. [Google Scholar] - Gepreel, K.A.; Nofal, T.A.; Alasmari, A.A. Exact solutions for nonlinear integro–partial differential equations using the generalized Kudryashov method. J. Egypt. Math. Soc.
**2017**, 25, 438–444. [Google Scholar] [CrossRef] - Yin, Y.H.; Ma, W.X.; Liu, J.G.; Lü, X. Diversity of exact solutions to a (3 + 1)-dimensional nonlinear evolution equation and its reduction. Comput. Math. Appl.
**2018**, 76, 1275–1283. [Google Scholar] [CrossRef] - Miah, M.M.; Ali, H.S.; Akbar, M.A.; Seadawy, A.R. New applications of the two variable (G
^{′}/G,1/G)-expansion method for closed form traveling wave solutions of integro-differential equations. J. Ocean Eng. Sci.**2019**, 4, 132–143. [Google Scholar] [CrossRef] - Alharbi, A.R.; Almatrafi, M.B. New exact and numerical solutions with their stability for Ito integro-differential equation via Riccati–Bernoulli sub-ODE method. J. Taibah Univ. Sci.
**2020**, 14, 1447–1456. [Google Scholar] [CrossRef] - Khalil, R.; Al Horani, M.; Yousef, A.; Sababheh, M. A new definition of fractional derivative. J. Comput. Appl. Math.
**2014**, 264, 65–70. [Google Scholar] [CrossRef] - Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
- Yang, S.; Wang, L.; Zhang, S. Conformable derivative: Application to non-Darcian flow in low-permeability porous media. Appl. Math. Lett.
**2018**, 79, 105–110. [Google Scholar] [CrossRef] - Anderson, D.R.; Camrud, E.; Ulness, D.J. On the nature of the conformable derivative and its applications to physics. arXiv
**2018**, arXiv:1810.02005. [Google Scholar] - Anderson, D.R.; Ulness, D.J. Newly defined conformable derivatives. Adv. Dyn. Syst. Appl.
**2015**, 10, 109–137. [Google Scholar] - Chung, W.S. Fractional Newton mechanics with conformable fractional derivative. J. Comput. Appl. Math.
**2015**, 290, 150–158. [Google Scholar] [CrossRef] - Lazo, M.J.; Torres, D.F. Variational calculus with conformable fractional derivatives. IEEE/CAA J. Autom. Sin.
**2016**, 4, 340–352. [Google Scholar] [CrossRef] - Çenesiz, Y.; Baleanu, D.; Kurt, A.; Tasbozan, O. New exact solutions of Burgers’ type equations with conformable derivative. Waves Random Complex Media
**2017**, 27, 103–116. [Google Scholar] [CrossRef] - Birgani, O.T.; Chandok, S.; Dedovic, N.; Radenovic, S. A note on some recent results of the conformable fractional derivative. Adv. Theory Nonlinear Anal. Appl.
**2019**, 3, 11–17. [Google Scholar] - Gao, F.; Chi, C. Improvement on conformable fractional derivative and its applications in fractional differential equations. J. Funct. Spaces
**2020**, 2020. [Google Scholar] [CrossRef] - Tasbozan, O.; Çenesiz, Y.; Kurt, A.; Baleanu, D. New analytical solutions for conformable fractional PDEs arising in mathematical physics by exp-function method. Open Phys.
**2017**, 15, 647–651. [Google Scholar] [CrossRef] - Hosseini, K.; Manafian, J.; Samadani, F.; Foroutan, M.; Mirzazadeh, M.; Zhou, Q. Resonant optical solitons with perturbation terms and fractional temporal evolution using improved tan (ϕ (η)/2)-expansion method and exp function approach. Optik
**2018**, 158, 933–939. [Google Scholar] [CrossRef] - Yaslan, H.Ç.; Girgin, A. Exp-function method for the conformable space-time fractional STO, ZKBBM and coupled Boussinesq equations. Arab. J. Basic Appl. Sci.
**2019**, 26, 163–170. [Google Scholar] [CrossRef] - Ghany, H.A.; Hyder, A.A.; Zakarya, M. Exact solutions of stochastic fractional Korteweg de-Vries equation with conformable derivatives. Chin. Phys. B
**2020**, 29, 030203. [Google Scholar] [CrossRef] - Biswas, A.; Sonmezoglu, A.; Ekici, M.; Mirzazadeh, M.; Zhou, Q.; Moshokoa, S.P.; Belic, M. Optical soliton perturbation with fractional temporal evolution by generalized Kudryashov’s method. Optik
**2018**, 164, 303–310. [Google Scholar] [CrossRef] - El-Ganaini, S.; Al-Amr, M.O. New abundant wave solutions of the conformable space-time fractional (4 + 1)-dimensional Fokas equation in water waves. Comput. Math. Appl.
**2019**, 78, 2094–2106. [Google Scholar] [CrossRef] - Wazwaz, A.M. Kadomtsev-Petviashvili hierarchy: N-soliton solutions and distinct dispersion relations. Appl. Math. Lett.
**2016**, 52, 74–79. [Google Scholar] [CrossRef] - Gepreel, K.A. Exact solutions for nonlinear integral member of Kadomtsev-Petviashvili hierarchy differential equations using the modified (w/g)-expansion method. Comput. Math. Appl.
**2016**, 72, 2072–2083. [Google Scholar] [CrossRef] - Zhang, Y.Y.; Liu, X.Q.; Wang, G.W. Symmetry reductions and exact solutions of the (2 + 1)-dimensional Jaulent-Miodek equation. Appl. Math. Comput.
**2012**, 219, 911–916. [Google Scholar] [CrossRef] - Gu, Y.; Deng, B.; Lin, J. Exact traveling wave solutions to the (2 + 1)-dimensional Jaulent-Miodek equation. Adv. Math. Phys.
**2018**, 2018. [Google Scholar] [CrossRef] - Ma, H.C.; Deng, A.P.; Yu, Y.D. Lie symmetry group of (2 + 1)-dimensional Jaulent-Miodek equation. Therm. Sci.
**2014**, 18, 1547–1552. [Google Scholar] [CrossRef] - Matinfar, M.; Eslami, M.; Roshandel, S. The first integral method to study the (2 + 1)-dimensional Jaulent-Miodek equations. Pramana
**2015**, 85, 593–603. [Google Scholar] [CrossRef] - Kaplan, M.; Bekir, A.; Akbulut, A. A generalized Kudryashov method to some nonlinear evolution equations in mathematical physics. Nonlinear Dyn.
**2016**, 85, 2843–2850. [Google Scholar] [CrossRef] - Sirisubtawee, S.; Koonprasert, S.; Sungnul, S. Some Applications of the (G
^{′}/G,1/G)-expansion method for finding exact traveling wave solutions of nonlinear fractional evolution equations. Symmetry**2019**, 11, 952. [Google Scholar] [CrossRef] - Sirisubtawee, S.; Koonprasert, S.; Sungnul, S.; Leekparn, T. Exact traveling wave solutions of the space-time fractional complex Ginzburg-Landau equation and the space-time fractional Phi-4 equation using reliable methods. Adv. Differ. Eqs.
**2019**, 2019, 219. [Google Scholar] [CrossRef] - Sirisubtawee, S.; Koonprasert, S.; Sungnul, S. New exact solutions of the conformable space-time Sharma-Tasso-Olver equation using two reliable methods. Symmetry
**2020**, 12, 644. [Google Scholar] [CrossRef] - Al-Ghafri, K.S. Soliton behaviours for the conformable space-time fractional complex Ginzburg-Landau equation in optical fibers. Symmetry
**2020**, 12, 219. [Google Scholar] [CrossRef] - He, J.H.; Elagan, S.; Li, Z. Geometrical explanation of the fractional complex transform and derivative chain rule for fractional calculus. Phys. Lett. A
**2012**, 376, 257–259. [Google Scholar] [CrossRef] - Su, W.H.; Yang, X.J.; Jafari, H.; Baleanu, D. Fractional complex transform method for wave equations on Cantor sets within local fractional differential operator. Adv. Differ. Eqs.
**2013**, 2013, 97. [Google Scholar] [CrossRef] - Güner, O.; Bekir, A.; Cevikel, A.C. A variety of exact solutions for the time fractional Cahn-Allen equation. Eur. Phys. J. Plus
**2015**, 130, 146. [Google Scholar] [CrossRef] - Rahmatullah; Mohyud-Din, S.T.; Khan, U. Exact traveling wave solutions of fractional order Boussinesq-like equations by applying Exp-function method. Results Phys.
**2018**, 8, 114–120. [Google Scholar] [CrossRef] - Guner, O.; Aksoy, E.; Bekir, A.; Cevikel, A.C. Different methods for (3 + 1)-dimensional space-time fractional modified KdV-Zakharov-Kuznetsov equation. Comput. Math. Appl.
**2016**, 71, 1259–1269. [Google Scholar] [CrossRef] - Guner, O.; Bekir, A. The Exp-function method for solving nonlinear space-time fractional differential equations in mathematical physics. J. Assoc. Arab. Univ. Basic Appl. Sci.
**2017**, 24, 277–282. [Google Scholar] [CrossRef] - Demiray, S.T.; Bulut, H. Soliton solutions of some nonlinear evolution problems by GKM. Neural Comput. Appl.
**2019**, 31, 287–294. [Google Scholar] [CrossRef] - He, J.H.; Wu, X.H. Exp-function method for nonlinear wave equations. Chaos Solitons Fractals
**2006**, 30, 700–708. [Google Scholar] [CrossRef]

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |

© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).