Skip to main content

Theory and Modern Applications

Efficient analytical techniques for solving time-fractional nonlinear coupled Jaulent–Miodek system with energy-dependent Schrödinger potential

Abstract

In this paper, we present analytical-approximate solution to the time-fractional nonlinear coupled Jaulent–Miodek system of equations which comes with an energy-dependent Schrödinger potential by means of a residual power series method (RSPM) and a q-homotopy analysis method (q-HAM). These methods produce convergent series solutions with easily computable components. Using a specific example, a comparison analysis is done between these methods and the exact solution. The numerical results show that present methods are competitive, powerful, reliable, and easy to implement for strongly nonlinear fractional differential equations.

1 Introduction and preliminaries

The term fractional calculus which involves fractional derivatives and fractional integral is nothing new. As stated in the letter L’Hospital wrote to Leibniz, in 1695, he asked him “what is the meaning of the expression \(d^{n} y/dx^{n}\) when \(n=1/2\)?” Leibniz replied to the L’Hospital letter telling him that “\(d^{1/2}x\) will be equal to \(x\sqrt{dx:x}\).” In reality, this is an apparent paradox, from this evident paradox, one day useful consequences will be drawn [1,2,3,4]. Since then, mathematicians have investigated this concept, the like of Riemann–Liouville, Caputo–Hadamard, Erdélyi–Kober, Grünwald–Letnikov, Fourier, Marchaud, Riesz, and Weyl, to mention a few. Most of these derivatives are defined on the basis of the corresponding fractional integral in the Riemann–Liouville sense. Recently, fractional calculus has attracted the attention of researchers in the various field of natural science and engineering due to it wide applications in these various mention fields. The applications can be found in anomalous transport, control theory of dynamical systems, signal and image processing, nanotechnology, financial modeling, viscoelasticity, random walk, nanoprecipitate growth in solid solutions, modeling for shape memory polymers, and anomalous diffusion just to mention a few. We refer the reader to [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20] and the references therein for more details.

Nonlinear partial differential equations (NPDEs) have been attracting attention of engineers, physicists and mathematicians in recent years. Among such NPDEs, we have the Jaulent–Miodek system of equations which comes with energy-dependent Schrödinger potential [21]. These systems of equations are extensively used as a model in solving many real world problems in various fields of engineering and natural sciences; see [22,23,24,25] and the references therein for more details. Due to the momentous and special position these equations have in the above-mentioned fields, it is important to understand the solutions (both the analytic-approximate and the numerical solutions) to these nonlinear partial differential equations (NPDEs). Comprehensive mathematical analysis of the nonlinear fractional-order coupled Jaulent–Miodek equations are still under study and it plays an important role in many parts of science and engineering such as plasma physics [25] and condensed matter physics [26, 27].

There are several methods used in obtaining approximate solutions to linear and nonlinear FPDE such as the Adomian decomposition method (ADM) [28, 29], the homotopy-perturbation method (HPM) [30,31,32,33,34], the variational iteration method (VIM) [35], the q-homotopy analysis transform method (q-HATM) [36, 37], the fractional natural decomposition method (FNDM) [38], the fractional multi-step differential transformed method (FMsDTM) [39], the new iterartive method (NIM) [40,41,42] and the homotopy analysis method (HAM) [43,44,45,46]. In a recent development, [47,48,49,50,51,52,53], a modified homotopy analysis method was established which has potential applications in a wide range of systems of differential equations. This method provides a convenient way to ascertain the convergence of approximation series and even exact solutions. This modification is called a q-homotopy analysis method (q-HAM), we refer this method as one of the most efficient methods of obtaining analytical-approximate and exact solutions for nonlinear partial fractional differential equations as well as the classical type.

Nonlinear time-fractional coupled Jaulent–Miodek sytem of equations where \(0<\alpha\leq1\), is defined as follows:

$$\begin{aligned}& \frac{{{\partial}^{\alpha}}u}{\partial{{t}^{\alpha}}}+\frac {{{\partial} ^{3}}u}{\partial{{x}^{3}}}+\frac{3}{2}v\frac{{{\partial }^{3}}v}{\partial{{x}^{3}}}+ \frac{9}{2}\frac{\partial v}{\partial x}\frac{{{\partial }^{2}}v}{\partial{{x}^{2}}}-6u\frac{\partial u}{\partial x}-6uv \frac{\partial v}{\partial x}-\frac{3}{2}\frac{\partial u}{\partial x}{{v}^{2}} =0, \\& \frac{{{\partial}^{\alpha}}v}{\partial{{t}^{\alpha}}}+\frac {{{\partial}^{3}}v}{\partial{{x}^{3}}}-6\frac{\partial u}{\partial x}v-6u\frac {\partial v}{\partial x}- \frac{15}{2}\frac{\partial v}{\partial x}{{v}^{2}} =0, \end{aligned}$$
(1)

which comes with energy-dependent Schrödinger potential [54,55,56]. Recently, the Sumudu transform homotopy-perturbation method (STHPM) [54], the Hermite wavelets method (HWM) and the optimal homotopy asymptotic method (OHAM) [57], the invariant subspace method [58], the q-homotopy analysis transform method (q-HATM) [59], and others [60,61,62] have been used to obtain approximate solutions of the nonlinear time-fractional Jaulent–Miodek system of equations.

In this paper, we present approximate solutions to time-fractional coupled Jaulent–Miodek equations using the residual power-series method (RSPM) and the q-homotopy analysis method (q-HAM). The paper is organized as follows. In Sect. 2, we explain residual power-series method (RSPM) and the q-homotopy analysis method (q-HAM), describe its convergence analysis, and we present example that shows reliability and efficiency of this method in order to obtain its stable numerical results. In Sect. 3, we obtain approximate solutions of the time-fractional nonlinear coupled Jaulent–Miodek system of equations using both the residual power-series method (RSPM) and the q-homotopy analysis method (q-HAM). In Sect. 4, we discuss the results obtained by RSPM and q-HAM, and Sect. 5 gives our conclusion.

Definition 1.1

The real function \(f(t)\), \(t>0\) is said to be in the space of \(C_{\mu }\) (\(\mu>0\)) when there exists a real number p (>μ) such that \(f(t)=t^{p}f_{1}(t)\) in which \(f_{1}\in C[0,\infty)\) and it is said to be in the space of \(C_{\mu}^{m} \) when \(f^{(m)}\in C_{\mu}\), \(m\in\aleph\) [3, 63].

Definition 1.2

The Riemann–Liouville fractional integral operator (\(J^{\alpha}\)) of order \(\alpha\geq0\) of a function \(f\in C_{\mu}\), \(\mu\geq-1\) is defined as [3, 63]

$$ J^{\alpha}f(t)=\frac{1}{\varGamma(\alpha)} \int_{0}^{t}(t-\tau)^{\alpha-1}f(\tau)\,d\tau, \quad\alpha, t>0, $$
(2)

and \(J^{0}f(t)=f(t)\), where Γ is the well-known Gamma function. Then the following properties hold for the function f:

$$f\in C_{\mu}, \qquad\mu\geq-1,\qquad\alpha, \beta\geq0\quad \text{and}\quad \lambda>-1. $$

Also, these general properties have been itemized as follows:

  1. (a)

    \(J^{\alpha}J^{\beta}f(t)=J^{\alpha+\beta}f(t)\),

  2. (b)

    \(J^{\alpha}J^{\beta}f(t)=J^{\beta}J^{\alpha}f(t)\),

  3. (c)

    \(J^{\alpha}t^{\lambda}=\frac{\varGamma(\lambda +1)}{\varGamma(\lambda +1+\alpha)}t^{\lambda+\alpha}\).

Definition 1.3

The fractional derivative of a function f of order α in the Caputo sense, for \(f\in C_{-1}^{m}\), \(m\in\aleph\cup\{0\}\) is defined as [3]

$$ D^{m\alpha}f(t)=J^{m-\alpha}f^{(m)}(t)=\frac{1}{\varGamma(m-\alpha )}\int_{0}^{t}(t-\tau)^{m-\alpha-1}f^{(m)}( \tau)\,d\tau,\quad \alpha, t>0, $$
(3)

where \(m-1<\alpha<m\) and the function f satisfies some defined properties as follows:

  1. (a)

    \(D^{\alpha} ( af(t)+bg(t) ) =aD^{\alpha }f(t)+bD^{\alpha }g(t)\), \(a,b\in\Re\),

  2. (b)

    \(D^{\alpha}J^{\alpha}f(t)=f(t)\),

  3. (c)

    \(J^{\alpha}D^{\alpha}f(t)=f(t)-\sum_{j=0}^{m-1}f^{(j)}(0) \frac{t^{j}}{j!}\), \(t>0\).

2 Analysis of approximate methods

2.1 Algorithm and convergence of RPSM

Here, using the RPS method, series solutions for the Jaulent–Miodek system of equations are obtained. The RPS method [64,65,66,67,68] consists of expressing the solution of Eq. (1) as a fractional power-series expansion about the initial point \(t=t_{0}\). It is worth mentioning that the proposed method can reduce the computational time and work as compared with other traditional techniques while maintaining the efficiency of the results obtained [69]. We have

$$ u(t,x)=\sum _{m=0}^{\infty}{\frac{D_{t}^{m\alpha }u_{m}(x,{{t}_{0}})}{\varGamma(m\alpha+1)}} {{ ( t-{{t}_{0}} ) }^{m\alpha}},\quad 0< \alpha\leq1, x\in[ a,b ] ,0\leq t< \Re, $$
(4)

and

$$ v(t,x)=\sum _{m=0}^{\infty}{\frac{D_{t}^{m\alpha }v_{m}(x,{{t}_{0}})}{\varGamma(m\alpha+1)}} {{ ( t-{{t}_{0}} ) }^{m\alpha}},\quad 0< \alpha\leq1, x\in[ a,b ] ,0\leq t< \Re. $$
(5)

The zeroth RPS approximate solutions of \(u(t,x)\) and \(v(t,x)\) can be written as follows:

$$\begin{aligned}& {{u}_{0}}(t,x)={{f}_{0}}(x)=f(x), \\& {{v}_{0}}(t,x)={{g}_{0}}(x)=g(x). \end{aligned}$$
(6)

Next, let \(u_{k}(t,x)\) and \(v_{k}(t,x)\) denote, respectively, the kth truncated series of \(u(t,x)\) and \(v(t,x)\) given as

$$\begin{aligned}& {{u}_{k}}(t,x)=\sum _{n=0}^{{k}}{{{f}_{n}}}(x) \frac{{{t}^{n\alpha}}}{\varGamma( 1+n\alpha) } , \end{aligned}$$
(7)
$$\begin{aligned}& {{v}_{k}}(t,x)=\sum _{n=0}^{{k}}{{{g}_{n}}}(x) \frac{{{t}^{n\alpha}}}{\varGamma(1+n\alpha) }. \end{aligned}$$
(8)

Substitution of the kth truncated series \(u_{k}(t,x)\) and \(v_{k}(t,x)\) of Eqs. (7) and (8) into the main equation Eq. (1) leads to the kth residual function denoted by \(\operatorname{Res} u_{k}(t,x)\) and \(\operatorname{Res} v_{k}(t,x)\) given by

$$\begin{aligned}& \begin{aligned}[b] \operatorname{Res}u_{k}(t,x) ={}&D_{t}^{\alpha }u(t,x)+(u_{k})_{xxx}(t,x)+ \frac{3 }{2}(v_{k})_{xxx}(t,x)v_{k}(t,x)+ \frac {9}{2}(v_{k})_{x}(t,x) (v_{k})_{xx}(t,x) \\ &- 6u_{k}(t,x) (u_{k})_{x}(t,x)-6{{u}_{k}}(t,x) (u_{k})_{x}(t,x)-6{{u}_{k}}(t,x){{v}_{k}}(t,x) (v_{k})_{x}(t,x) \\ &- \frac{3}{2}(u_{k})_{x}(t,x)v_{k}^{2}(t,x), \end{aligned} \hspace{-6pt} \\& \begin{aligned} \operatorname{Res}v_{k}(t,x) ={}&D_{t}^{\alpha }v(t,x)+(v_{k})_{xxx}(t,x)-6(u_{k})_{x}(t,x)v_{k}(t,x)-6u_{k}(t,x) (v_{k})_{x}(t,x) \\ &- \frac{15}{2}(v_{k})_{x}(t,x)v_{k}^{2}(t,x). \end{aligned} \end{aligned}$$
(9)

Also, we obtain

$$ \operatorname{Res}u(t,x)=\lim_{k\rightarrow\infty} \operatorname{Res}_{k}(t,x). $$
(10)

It is noticed that \(\operatorname{Res}u(t,x)=0\) for all values of \(x\in[ a,b]\). This means that \(\operatorname{Res}u(t,x)\) is infinitely many times differentiable at \(x=a\). Besides, \(\frac {{{d}^{k}}}{d{{x}^{k-1}}}\operatorname{Res}u(0,x)=\frac {{{d}^{k}}}{d{{x}^{k-1}}}\operatorname{Res}u_{k}(0,x)=0\). In fact, this equation is a fundamental rule in RPS method in order to apply it on many linear and nonlinear problems. To obtain the kth approximate solutions, we consider Eqs. (7) and (8) and then we differentiate both sides of these equations with respect to independent variables x and t and then we substitute \(t=0\) in order to find f and g constant parameters. After substituting these constant parameters in \(u_{k}(t,x)\), we can obtain the kth truncated series and by putting it in Eq. (1), we get our favorite approximate solution. This procedure can be iterated for other arbitrary order of coefficients with RPS solutions of Eq. (1). Now, regarding the convergence of the above iteration scheme, we present the following theorem.

Theorem 2.1

If there exists a fixed constant \(0< K<1\)such that

$$\bigl\Vert {{u}_{n+1}} ( t,x ) \bigr\Vert \leq K \bigl\Vert {{u}_{n}} ( t,x ) \bigr\Vert $$

for all \(n\in\mathbb{N}\)and \(0< t< R<1\), then the sequence of approximate solution converges to an exact solution.

Proof

For all \(0< t< R<1\), we have

$$\begin{aligned} \bigl\Vert u(t,x)-{{u}_{n}}(t,x) \bigr\Vert =& \Biggl\Vert \sum _{i=n+1}^{\infty}{{{u}_{i}}(t,x)} \Biggr\Vert \\ \leq&\sum _{i=n+1}^{\infty}{ \bigl\Vert {{u}_{i}}(t,x) \bigr\Vert } \\ \leq& \bigl\Vert {f(x)} \bigr\Vert \sum _{i=n+1}^{\infty }{{{K}^{i}}} \\ =&\frac{{{K}^{n+1}}}{1-K} \bigl\Vert {f(x)} \bigr\Vert \overset{n\rightarrow \infty}{\longrightarrow}0. \end{aligned}$$
(11)

 □

2.2 Fundamentals of the q-HAM

Here, we give a brief analysis of the q-homotopy analysis method as applied to differential equations. Generally, we consider the following differential equation:

$$ \mathcal{N} \bigl[\mathcal {D}_{t}^{\alpha}w(t,x) \bigr]-f(t,x)=0, $$
(12)

where \(\mathcal {D}_{t}^{\alpha}\) is the fractional derivative in time, \(\mathcal{N}\) denotes the nonlinear operator, f is a given function and \(w(t,x)\) is an unknown function. This is the the zeroth-order deformation equation

$$\begin{aligned} (1-nq)\mathcal{L} \bigl(\varphi(t,x;q)-w_{0}(t,x) \bigr)=qhH(t,x) \bigl(\mathcal{N}\bigl[\mathcal {D}_{t}^{\alpha}\varphi(t,x;q) \bigr]-f(t,x) \bigr), \end{aligned}$$
(13)

where \(n\geqslant1\), \(q\in[0,\frac{1}{n} ]\) denotes the so-called embedded parameter, \(h\neq0\) is an auxiliary parameter, \(H(t,x)\) is a non-zero auxiliary function, L is an auxiliary linear operator.

The following equations are obtained for \(q=0\) and \(q=\frac{1}{n}\), respectively:

$$ \varphi(t,x;0)=w_{0}(t,x) \quad\text{and}\quad\varphi\biggl(t,x; \frac {1}{n} \biggr)=w(t,x) . $$
(14)

So, starting from 0, as q approaches \(\frac{1}{n}\), the solution \(\varphi(t,x;q)\) varies from the initial guess \(w_{0}(t,x)\) to the solution \(w(t,x)\). We need to carefully choose \(w_{0}(t,x)\), \(\mathcal {L}\), h, \(H(t,x)\) to ascertain the existence of the solution \(\varphi (t,x;q)\) of Eq. (12) for \(q\in[0,\frac{1}{n} ]\). The Taylor series expansion of \(\varphi(t,x;q)\) gives

$$ \varphi(x,t;q)=w_{0}(t,x)+\sum_{m=1}^{\infty}w_{m}(t,x)q^{m}, $$
(15)

where

$$ w_{m}(t,x)=\frac{1}{m!}\frac{\partial^{m}\varphi (t,x;q)}{\partial q^{m}} \bigg|_{q=0}. $$
(16)

Assume that the auxiliary linear operator L, the initial guess \(c_{0}\), the auxiliary parameter h and \(H(t,x)\) are properly chosen such that the series Eq. (15) converges at \(q=\frac{1}{n}\), then we have

$$ w(t,x)=w_{0}(t,x)+\sum_{m=1}^{\infty}w_{m}(t,x) \biggl(\frac{1}{n} \biggr)^{m}. $$
(17)

Let the vector \(\vec{w}_{n}\) be defined as follows:

$$ \vec{w}_{n}= \bigl\{ w_{0}(t,x),w_{1}(t,x), \ldots,w_{n}(t,x) \bigr\} . $$
(18)

First, differentiate Eq. (13) m times with respect to the parameter q, then evaluate at \(q=0\) and finally divide them by m!. We have what is known as the mth-order deformation equation [52, 70]

$$ \mathcal{L} \bigl[w_{m}(t,x)-\chi^{*}_{m}w_{m-1}(t,x) \bigr]=hH(t,x)\mathcal{R}_{m} (\vec{w}_{m-1} ), $$
(19)

with initial conditions

$$ w_{m}^{(k)}(0,x)=0,\quad k=0,1,2,\ldots,m-1, $$
(20)

where

$$ \mathcal{R}_{m} (\vec{w}_{m-1} )=\frac {1}{(m-1)!} \frac{\partial^{m-1} (\mathcal{N}[\mathcal {D}_{t}^{\alpha}\varphi(t,x;q)]-g(t,x) )}{\partial q^{m-1}} \bigg|_{q=0} $$
(21)

and

$$\begin{aligned} \chi^{*}_{m}=\left \{ \textstyle\begin{array}{l@{\quad}l} 0,& m\leqslant1, \\ n,&\text{otherwise}. \end{array}\displaystyle \right . \end{aligned}$$
(22)

3 Solutions of time-fractional coupled Jaulent–Miodek system of equations

This section presents application of the above approximate methods for obtaining solutions of the time-fractional coupled Jaulent–Miodek system of equations.

3.1 RPSM solution

Consider the following time-fractional coupled Jaulent–Miodek (JM) system of equations:

$$\begin{aligned}& \frac{\partial^{\alpha}u}{\partial t^{\alpha}}+\frac {\partial ^{3}u}{\partial x^{3}}+\frac{3}{2}v \frac{\partial^{3}v}{\partial x^{3}}+\frac{9}{2}\frac{\partial v}{\partial x} \frac{\partial^{2}v}{\partial x^{2}}-6u\frac{\partial u}{\partial x}-6uv\frac{\partial v}{\partial x}-\frac {3}{2} \frac{\partial u}{\partial x}v^{2}=0, \\& \frac{\partial^{\alpha}v}{\partial t^{\alpha}}+\frac{\partial ^{3}v}{\partial x^{3}}-6\frac{\partial u}{\partial x}v-6u \frac{\partial v}{\partial x}-\frac{15}{2}\frac{\partial v}{\partial x}v^{2}=0, \end{aligned}$$
(23)

subject to the initial conditions

$$ u(0,x)=f(x),\qquad v(0,x)=g(x). $$
(24)

The method reveals the solution of the problem as a fractional power series around the initial value \(t=0\) as

$$\begin{aligned}& u(t,x)=\sum_{n=o}^{\infty}f_{n}(x) \frac{t^{n\alpha}}{\varGamma (1+n\alpha)}, \end{aligned}$$
(25)
$$\begin{aligned}& v(t,x)=\sum_{n=o}^{\infty}g_{n}(x) \frac{t^{n\alpha}}{\varGamma (1+n\alpha)}. \end{aligned}$$
(26)

Next, let \(u_{k}(t,x)\) and \(v_{k}(t,x)\) denote the kth truncated series of \(u(t,x)\) and \(v(t,x)\), respectively, we have

$$\begin{aligned}& u_{k}(t,x)=\sum_{n=o}^{k}f_{n}(x) \frac{t^{n\alpha}}{\varGamma (1+n\alpha)}, \end{aligned}$$
(27)
$$\begin{aligned}& v_{k}(t,x)=\sum_{n=o}^{k}g_{n}(x) \frac{t^{n\alpha}}{\varGamma (1+n\alpha)} . \end{aligned}$$
(28)

We define the kth residual functions as

$$\begin{aligned}& \begin{aligned}[b] \operatorname{Res} u_{k}(t,x) ={}& D_{t}^{\alpha }u(t,x)+(u_{k})_{xxx}(t,x)+ \frac{3}{2}(v_{k})_{xxx}(t,x)v_{k}(t,x)+ \frac {9}{2}(v_{k})_{x}(t,x) (v_{k})_{xx}(t,x)\hspace{-12pt} \\ &-6u_{k}(t,x) (u_{k})_{x}(t,x)-6u_{k}(t,x) (v_{k})_{x}(t,x)v_{k}(t,x)\\&- \frac {3}{2}(u_{k})_{x}(t,x)v_{k}^{2}(t,x), \end{aligned} \end{aligned}$$
(29)
$$\begin{aligned}& \begin{aligned}[b] \operatorname{Res} v_{k}(t,x) ={}& D_{t}^{\alpha }v(t,x)+(v_{k})_{xxx}(t,x)-6(u_{k})_{x}(t,x)v_{k}(t,x)-6u_{k}(t,x) (v_{k})_{x}(t,x) \\ &-\frac{15}{2}(v_{k})_{x}(t,x)v_{k}^{2}(t,x). \end{aligned} \end{aligned}$$
(30)

Since \(\operatorname{Res}u(t,x)=0\) and \(\lim_{k\rightarrow\infty }\operatorname{Res} u_{k}(t,x)=\operatorname{Res}u(t,x)\) for all \(x\in I\) and \(t\geq0\) [71, 72], \(D_{t}^{m\alpha}\operatorname {Res}u(t,x)=0\), because the fractional derivative of a constant is 0 in the Caputo sense. Also, the fractional derivative \(D_{t}^{m\alpha}\operatorname{Res}u(t,x)=0\) and \(\operatorname {Res}u_{k}(t,x)\) is matching at \(t=0\) for each \(m=0,1,2,\ldots,k\). We first substitute \(u_{k}(t,x)\) and \(v_{k}(t,x)\) into Eq. (1) and find the fractional derivative formula \(D_{t}^{(k-1)\alpha}\operatorname{Res}u(0,x)=0\) for \(k=1,2,3\). Solving these algebraic equations gives the \(f_{n}(x)\) and \(g_{n}(x)\) coefficients.

For the first step, we consider \(k=1\),

$$\begin{aligned}& u_{1}(t,x)=f(x)+f_{1}(x)\frac{t^{\alpha}}{\varGamma(\alpha +1)}, \end{aligned}$$
(31)
$$\begin{aligned}& v_{1}(t,x)=g(x)+g_{1}(x)\frac{t^{\alpha}}{\varGamma(\alpha+1)} , \end{aligned}$$
(32)

and

$$\begin{aligned}& \begin{aligned}[b]\operatorname{Res} u_{1}(0,x) ={}&f_{1}(x)-6f^{\prime}(x)f(x)- \frac {3}{2}f^{\prime }(x)g^{2}(x)-6f(x)g(x)g^{\prime}(x)+ \frac{9}{2}g^{\prime }(x)g^{\prime \prime}(x) \\ &+f^{(3)}(x)+\frac{3}{2}g^{(3)}(x)g(x), \end{aligned} \end{aligned}$$
(33)
$$\begin{aligned}& \operatorname{Res} v_{1}(0,x) =g_{1}(x)-6f^{\prime }(x)g(x)-6f(x)g^{\prime }(x)-\frac{15}{2}g^{2}(x)g^{\prime}(x)+g^{(3)}(x). \end{aligned}$$
(34)

Therefore

$$\begin{aligned}& f_{1}(x) =\frac{1}{2} \bigl( -2f^{(3)}+3g^{2}f^{\prime }+12ff^{\prime }+12fgg^{\prime(3)}g-9g^{\prime}g^{\prime\prime} \bigr) , \end{aligned}$$
(35)
$$\begin{aligned}& g_{1}(x) =\frac{1}{2} \bigl( 12gf^{\prime}+12fg^{\prime (3)}+15g^{2}g^{\prime} \bigr) , \end{aligned}$$
(36)

and

$$\begin{aligned}& \begin{aligned}[b] u_{1}(t,x) ={}&f(x)+\frac{t^{\alpha} ( 12f(x)f^{\prime }(x)+3f^{\prime }(x)g^{2}(x)+12f(x)g(x)g^{\prime}(x) ) }{2\varGamma(\alpha+1)} \\ &-\frac{t^{\alpha} (9g^{\prime}(x)g^{\prime\prime }(x)+2f^{(3)}(x)+3g(x)g^{(3)}(x) ) }{2\varGamma(\alpha+1)}, \end{aligned} \\& v_{1}(t,x) =g(x)+\frac{t^{\alpha} ( 12f^{\prime }(x)g(x)+12f(x)g^{\prime}(x)+15g^{2}(x)g^{\prime }(x)-2g^{(3)}(x) ) }{2\varGamma(\alpha+1)} . \end{aligned}$$
(37)

To obtain \(f_{2}(x)\) and \(g_{2}(x)\), we substitute the second truncated series,

$$\begin{aligned}& u_{2}(t,x) =f(x)+f_{2}(x)\frac{t^{2\alpha}}{\varGamma(2\alpha +1)}+f_{1}(x)\frac{t^{\alpha}}{\varGamma(\alpha+1)} , \end{aligned}$$
(38)
$$\begin{aligned}& v_{2}(t,x) =g(x)+g_{2}(x)\frac{t^{2\alpha}}{\varGamma(2\alpha +1)}+g_{1}(x)\frac{t^{\alpha}}{\varGamma(\alpha+1)} , \end{aligned}$$
(39)

into the second residual functions \(\operatorname{Res}u_{2}(t,x)\) and \(\operatorname{Res}v_{2}(t,x) \) and applying \(D_{t}^{\alpha}\) on both sides for \(t=0\) yields

$$\begin{aligned}& \begin{aligned}[b] \operatorname{Res} u_{2}(0,x) ={}&f_{2}(x)-f^{(6)}(x)+12f^{(4)}(x)g(x)^{2}+12f(x)f^{(4)}(x)+84f^{(3)}(x)g(x)g^{\prime }(x) \\ &+108f^{\prime\prime}(x)g^{\prime\prime}(x)g(x)+108f^{\prime \prime }(x)g^{\prime2}(x)- \frac{9}{4}^{\prime\prime2}f^{\prime\prime }(x)g^{4}(x) \\ &+18f^{\prime\prime2}-36f(x)^{2}f^{\prime\prime(3)}(x)f^{\prime 3}f^{\prime}(x)g^{\prime}(x)g(x)-252f(x)f^{\prime}(x)g^{\prime}(x)g(x) \\ &+225f^{\prime}(x)f^{\prime2}(x)g^{\prime}(x)g^{\prime\prime 2}-72f(x)f^{\prime2}+30f^{(3)}(x)f^{\prime(4)}(x)g(x) \\ &-54f(x)g(x)^{3}g^{\prime\prime2}g(x)g^{\prime\prime }(x)+72f(x)g^{\prime\prime2}-180f(x)g(x)^{2}g^{\prime 2}(x) \\ &-72f(x)^{2}g^{\prime2}(x)+102f(x)g^{(3)}(x)g^{\prime }(x)-3g^{(6)}(x)g(x)+ \frac{27}{2}g^{(4)}(x)g^{3}(x)\hspace{-12pt} \\ &-\frac{33}{2}g^{(3)}(x)^{2}+108g(x)^{2}g^{\prime\prime2}(x)+ \frac {135}{2}g^{\prime4}(x)+432g(x)g^{\prime2}(x)g^{\prime\prime}(x) \\ &-\frac{27}{2}g^{(5)}(x)g^{\prime(4)}(x)g^{\prime\prime (3)}(x)g(x)^{2}g^{\prime}(x), \end{aligned} \end{aligned}$$
(40)
$$\begin{aligned}& \begin{aligned}[b] \operatorname{Res} v_{2}(0,x) ={}&g_{2}(x)+12f^{(4)}(x)g(x)+30f^{(3)}(x)g^{\prime }(x)+36f^{\prime\prime}(x)f^{\prime\prime}(x)g^{\prime\prime 3}(x) \\ &-72f(x)f^{\prime\prime(3)}(x)f^{\prime2}g^{\prime\prime 3}(x)f^{\prime}(x)g^{\prime}(x)g(x)-144f(x)f^{\prime}(x)g^{\prime }(x) \\ &-72f^{\prime 2}(x)g(x)+12f(x)g^{(4)}(x)-126f(x)g(x)^{2}g^{\prime\prime 2}(x)g^{\prime \prime}(x) \\ &-252f(x)g(x)g^{\prime2}(x)-g^{(6)}(x)+24g^{(4)}(x)g(x)^{2}- \frac {225}{4}g(x)^{4}g^{\prime\prime }(x) \\ &+72g(x)g^{\prime\prime2}(x)-225g(x)^{3}g^{\prime 2}(x)+120g^{(3)}(x)g(x)g^{\prime}(x)+117g^{\prime2}g^{\prime\prime }(x), \end{aligned} \end{aligned}$$
(41)

and

$$\begin{aligned}& \begin{aligned}[b] f_{2}(x) ={}&f^{(6)}-12f^{(4)}g^{2}-12ff^{(4)}-84f^{(3)}gg^{\prime }+ \frac{9g^{4}f^{\prime\prime}}{4}+54fg^{2}f^{\prime\prime2}gg^{\prime \prime } \\ &-108gf^{\prime\prime}g^{\prime\prime}-108f^{\prime\prime} \bigl( g^{\prime} \bigr) ^{2} +72f^{2} \bigl( g^{\prime} \bigr) ^{2}-18 \bigl( f^{\prime\prime} \bigr) ^{2}+36f^{2}f^{\prime\prime(3)}gf^{\prime2} \bigl( f^{\prime} \bigr) ^{2} \\ &+252fgf^{\prime}g^{\prime3}f^{\prime}g^{\prime}-225f^{\prime }g^{\prime}g^{\prime\prime}+72f \bigl( f^{\prime} \bigr) ^{2}-30f^{(3)}f^{\prime(4)}g-72f \bigl( g^{\prime\prime} \bigr) ^{2} \\ &+54fg^{3}g^{\prime\prime(3)}g^{\prime2} \bigl( g^{\prime} \bigr) ^{2}+3g^{(6)}g+\frac{33}{2} \bigl( g^{(3)} \bigr) ^{2}-108g^{2} \bigl( g^{\prime\prime} \bigr) ^{2}, \end{aligned} \end{aligned}$$
(42)
$$\begin{aligned}& \begin{aligned}[b] g_{{2}}(x) ={}&{-}30f^{(3)}g^{\prime2}g^{\prime\prime3} \bigl( 6f^{\prime \prime}+25 \bigl( g^{\prime} \bigr) ^{2} \bigr) -9g^{\prime\prime } \bigl( 4f^{\prime\prime}+13 \bigl( g^{\prime} \bigr) ^{2} \bigr) -30g^{(3)}f^{\prime} \\ &+12f \bigl( 12f^{\prime}g^{\prime(4)} \bigr)+6g^{2} \bigl( 48f^{\prime}g^{\prime}+21fg^{\prime\prime(4)} \bigr) \\ &+12g \bigl( -f^{(4)}+6ff^{\prime\prime}+6 \bigl( f^{\prime} \bigr) ^{2}+21f \bigl( g^{\prime} \bigr) ^{2}-6 \bigl( g^{\prime\prime } \bigr) ^{2}-10g^{(3)}g^{\prime} \bigr) \\ &+g^{(6)}+\frac{225g^{4}g^{\prime\prime}}{4} . \end{aligned} \end{aligned}$$
(43)

Therefore

$$\begin{aligned}& \begin{aligned}[b] u_{{2}}(t,x) ={}&f(x)\\&+\frac{t^{\alpha} ( -2f^{(3)}(x)+12f(x) ( f^{\prime}(x)+g(x)g^{\prime}(x) ) +3g(x)^{2}f^{\prime (3)}(x)g(x)-9g^{\prime}(x)g^{\prime\prime}(x) ) }{2\varGamma (\alpha+1)}\hspace{-24pt} \\ &+\frac{1}{4\varGamma(2\alpha+1)}(4f^{(6)}(x)t^{2\alpha }-48f^{(4)}(x)g(x)^{2}t^{2\alpha}-48f(x)f^{(4)}(x)t^{2\alpha} \\ &-336f^{(3)}(x)g(x)t^{2\alpha}g^{\prime2\alpha}f^{\prime\prime }(x)g^{\prime\prime2\alpha}f^{\prime\prime}(x)g^{\prime 2}+9g(x)^{4}t^{2\alpha}f^{\prime\prime2}t^{2\alpha}f^{\prime\prime 2\alpha}f^{\prime\prime2} \\ &+144f(x)^{2}t^{2\alpha}f^{\prime\prime(3)}(x)t^{2\alpha}f^{\prime 3}t^{2\alpha}f^{\prime}(x)g^{\prime2\alpha}f^{\prime}(x)g^{\prime 2\alpha}f^{\prime}(x)g^{\prime}(x)g^{\prime\prime2}t^{2\alpha }f^{\prime2} \\ &+288f(x)t^{2\alpha}f^{\prime2}-120f^{(3)}(x)t^{2\alpha}f^{\prime (4)}(x)t^{2\alpha}-288f(x)t^{2\alpha}g^{\prime\prime2} \\ &+216f(x)g(x)^{3}t^{2\alpha}g^{\prime\prime2}g(x)t^{2\alpha }g^{\prime \prime2}t^{2\alpha}g^{\prime2}+720f(x)g(x)^{2}t^{2\alpha}g^{\prime 2} \\ &-408f(x)g^{(3)}(x)t^{2\alpha}g^{\prime(6)}(x)t^{2\alpha }-54g(x)^{3}g^{(4)}(x)t^{2\alpha}+66g^{(3)}(x)^{2}t^{2\alpha } \\ &-432g(x)^{2}t^{2\alpha}g^{\prime\prime2}-270t^{2\alpha}g^{\prime 4}\\&+54g^{(5)}(x)t^{2\alpha}g^{\prime(4)}(x)t^{2\alpha}g^{\prime \prime 2}g^{(3)}(x)t^{2\alpha}g^{\prime2\alpha}g^{\prime2}g^{\prime\prime}(x), \end{aligned} \end{aligned}$$
(44)
$$\begin{aligned}& v_{2}(t,x) =g(x)+\frac{t^{\alpha} ( 12g(x)f^{\prime }(x)+12f(x)g^{\prime(3)}(x)+15g(x)^{2}g^{\prime}(x) ) }{2\varGamma (\alpha+1)} \\& \phantom{v_{2}(t,x) =}+\frac{1}{4\varGamma(2\alpha +1)}(-48f^{(4)}(x)g(x)t^{2\alpha }-120f^{(3)}(x)g^{\prime}(x)t^{2\alpha}-144f^{\prime\prime }(x)g^{\prime \prime3}f^{\prime\prime}(x)t^{2\alpha} \\& \phantom{v_{2}(t,x) =}+288f(x)g(x)f^{\prime\prime(3)}(x)f^{\prime 2}f^{\prime }(x)g^{\prime }(x)t^{2\alpha}+576f(x)f^{\prime}(x)g^{\prime}(x)t^{2\alpha }+288g(x)f^{\prime2}t^{2\alpha} \\& \phantom{v_{2}(t,x) =}-48f(x)g^{(4)}(x)t^{2\alpha }+504f(x)g(x)^{2}g^{\prime\prime 2}g^{\prime \prime}(x)t^{2\alpha}+1008f(x)g(x)g^{\prime2}t^{2\alpha} \\& \phantom{v_{2}(t,x) =}+4g^{(6)}(x)t^{2\alpha }-96g^{(4)}(x)g(x)^{2}t^{2\alpha }+225g(x)^{4}g^{\prime\prime}(x)t^{2\alpha}-288g(x)g^{\prime\prime 2}t^{2\alpha} \\& \phantom{v_{2}(t,x) =}+900g(x)^{3}g^{\prime2}t^{2\alpha }-480g^{(3)}(x)g(x)g^{\prime }(x)t^{2\alpha}-468g^{\prime2}g^{\prime\prime}(x)t^{2\alpha}) . \end{aligned}$$
(45)

Applying the same procedure we finally calculated the following values:

$$\begin{aligned}& \begin{aligned}[b] u_{3}(t,x) ={}&f(x)+\frac{f_{1}(x)t^{\alpha}}{\varGamma(\alpha +1)}+\frac{f_{2}(x)t^{2\alpha}}{\varGamma(2\alpha+1)}+ \frac{t^{3\alpha}}{2\varGamma (\alpha+1)^{2}\varGamma(3\alpha+1)} \\ &\times\bigl(6\varGamma(\alpha+1)^{2}g(x)g_{2}(x)f^{\prime}(x)\\&+6 \alpha\varGamma(2\alpha)g_{1}(x)^{2}f^{\prime2}f_{2}(x)f^{\prime 2}f_{2}(x)g(x)g^{\prime 2}f(x)g_{2}(x)g^{\prime}(x) \\ &+24\alpha\varGamma(2\alpha)f_{1}(x)g_{1}(x)g^{\prime 2}g(x)^{2}f_{2}^{\prime2}f(x)g(x)g_{2}^{\prime}(x)+12 \alpha\varGamma(2\alpha)g(x)g_{1}(x)f_{1}^{\prime}(x)\hspace{-12pt} \\ &+24\alpha\varGamma(2\alpha)f_{1}(x)g(x)g_{1}^{\prime}(x)+24 \alpha\varGamma(2\alpha)f(x)g_{1}(x)g_{1}^{\prime 2}f(x)f_{2}^{\prime}(x) \\ &+24\alpha\varGamma(2\alpha)f_{1}(x)f_{1}^{\prime2}f_{2}^{(3)}(x)-3 \varGamma(\alpha+1)^{2}g_{2}(x)g^{(3)}(x) \\ &-9\varGamma(\alpha+1)^{2}g_{2}^{\prime}(x)g^{\prime\prime 2}g^{\prime}(x)g_{2}^{\prime \prime}(x)-18 \alpha\varGamma(2\alpha)g_{1}^{\prime}(x)g_{1}^{\prime \prime 2}g(x)g_{2}^{(3)}(x) \\ &-6\alpha\varGamma(2\alpha)g_{1}(x)g_{1}^{(3)}(x)\bigr), \end{aligned} \end{aligned}$$
(46)
$$\begin{aligned}& \begin{aligned}[b] v_{{3}}(t,x) ={}&g(x)+\frac{g_{1}(x)t^{\alpha}}{\varGamma(\alpha +1)}+\frac{g_{2}(x)t^{2\alpha}}{\varGamma(2\alpha+1)}+ \frac{t^{3\alpha}}{2\varGamma (\alpha+1)^{2}\varGamma(3\alpha+1)} \\ &\times\bigl(12\varGamma(\alpha+1)^{2}g_{2}(x)f^{\prime 2}f_{2}(x)g^{\prime 2}g(x)f_{2}^{\prime2}f(x)g_{2}^{\prime}(x)+24 \alpha\varGamma(2\alpha)g_{1}(x)f_{1}^{\prime}(x) \\ &+24\alpha\varGamma(2\alpha)f_{1}(x)g_{1}^{\prime 2}g(x)g_{2}(x)g^{\prime }(x)+30 \alpha\varGamma(2\alpha)g_{1}(x)^{2}g^{\prime 2}g(x)^{2}g_{2}^{\prime}(x) \\ &+60\alpha\varGamma(2\alpha)g(x)g_{1}(x)g_{1}^{\prime 2}g_{2}^{(3)}(x)\bigr). \end{aligned} \end{aligned}$$
(47)

3.2 q-HAM solution

Consider the same Eq. (23), we use initial approximations

$$\begin{aligned} u_{0}(t,x):=u(0,x)=f(x) \end{aligned}$$

and

$$\begin{aligned} v_{0}(t,x):=v(0,x)=g(x). \end{aligned}$$

We apply q-HAM and obtain the following:

$$\begin{aligned}& \begin{aligned}[b] \mathcal{R}_{m} (\vec{u}_{m-1} )={}& \mathcal {D}_{t}^{\alpha}u_{m-1}+u_{(m-1)xxx}+1.5\sum _{k=0}^{m-1}v_{k}v_{(m-1-k)xxx}+4.5 \sum_{k=0}^{m-1}v_{kx}v_{(m-1-k)xx} \\ &-6\sum_{k=0}^{m-1}u_{k}u_{(m-1-k)x}-6 \sum_{k=0}^{m-1}\sum _{j=0}^{k}u_{j}v_{k-j}v_{(m-1-k)x} \\ &-1.5\sum_{k=0}^{m-1}\sum _{j=0}^{k}v_{j}v_{k-j}u_{(m-1-k)x}, \end{aligned} \end{aligned}$$
(48)
$$\begin{aligned}& \begin{aligned}[b] \hat{\mathcal{R}}_{m} (\vec{v}_{m-1} )={}& \mathcal {D}_{t}^{\alpha}v_{m-1}+v_{(m-1)xxx}-6\sum _{k=0}^{m-1}v_{k}u_{(m-1-k)x}-6 \sum_{k=0}^{m-1}u_{kx}v_{(m-1-k)x} \\ &-7.5\sum_{k=0}^{m-1}\sum _{j=0}^{k}v_{j}v_{k-j}v_{(m-1-k)x}. \end{aligned} \end{aligned}$$
(49)

Therefore, we have the solution to the system for \(m\geq1\)

$$\begin{aligned}& u_{m}(t,x)=\chi^{*}_{m} u_{m-1}(t,x)+h J^{\alpha}\bigl[\mathcal{R}_{m} \bigl(\vec{u}_{m-1}(t,x) \bigr) \bigr], \end{aligned}$$
(50)
$$\begin{aligned}& v_{m}(t,x)=\chi^{*}_{m} v_{m-1}(t,x)+h J^{\alpha}\bigl[\hat{\mathcal{R}}_{m} \bigl( \vec{v}_{m-1}(t,x) \bigr) \bigr]. \end{aligned}$$
(51)

Hence, the expression of the series solutions by q-HAM are

$$\begin{aligned}& u(t,x;n;h)=U^{(M)}(t,x,n,h)=u_{0}(t,x)+\sum _{i=1}^{M}u_{i}(t,x;n;h) \biggl( \frac{1}{n} \biggr)^{i}, \end{aligned}$$
(52)
$$\begin{aligned}& v(t,x;n;h)=V^{(M)}(t,x,n,h)=v_{0}(t,x)+\sum _{i=1}^{M}v_{i}(t,x;n;h) \biggl( \frac{1}{n} \biggr)^{i}. \end{aligned}$$
(53)

The series solutions obtained in Eqs. (52) and (53) are appropriate solutions to the system Eq. (23) in terms of the convergence parameter h and n.

4 Numerical comparison

For numerical comparison purposes, we consider the system of equations 23 with different initial data.

4.1 Case I

We begin with a case where exact solutions are known when \(\alpha=1\). Consider the following initial conditions:

$$\begin{aligned}& u(0,x) =f(x)=\frac{1}{8}\lambda^{2} \biggl( 1-4\sec h^{2} \biggl( \frac{ \lambda x}{2} \biggr) \biggr), \\& v(0,x) =g(x)=\lambda\sec h \biggl( \frac{\lambda x}{2} \biggr). \end{aligned}$$
(54)

The exact solutions of the problem (\(\alpha=1\)) are given by

$$\begin{aligned}& u(t,x) =\frac{1}{8}\lambda^{2} \biggl( 1-4\sec h^{2} \biggl( \frac{1}{2} \lambda\biggl( x+\frac{1}{2} \lambda^{2}t \biggr) \biggr) \biggr), \\& v(t,x) =\lambda\sec h \biggl( \frac{1}{2}\lambda\biggl( x+ \frac {1}{2}\lambda^{2}t \biggr) \biggr), \end{aligned}$$
(55)

where λ is an arbitrary constant. The absolute errors for this Case I are reported in Tables 1, 2. In Figs. 1, 2, 3, 4, we present the graphical representation of the obtained results by RPSM, q-HAM and the exact solutions. We compare solutions, \((U^{(3)},V^{(3)})\), obtained by RPSM and q-HAM with the exact solutions by using different values of the parameters x and t.

Figure 1
figure 1

Comparison for Case I: \(u(t,x)\), \(\lambda=0.5\), \(\alpha =1\)

Figure 2
figure 2

Comparison for Case I: \(v(t,x)\), \(\lambda=0.5\), \(\alpha =1\)

Figure 3
figure 3

RPSM vs. exact solution for Case I in 2D with \(\lambda=0.5\), \(\alpha=1\)

Figure 4
figure 4

q-HAM vs. exact solution for Case I in 2D with \(\lambda =0.5\), \(\alpha=1\)

Table 1 Case I Comparison (\(U^{(3)}\)): RPSM, q-HAM with the exact solution
Table 2 Case I Comparison (\(V^{(3)}\)): RPSM, q-HAM with the exact solution

4.2 Case II

$$\begin{aligned}& u(0,x) =f(x)=e^{\gamma x}, \\& v(0,x) =g(x)=e^{\beta x}, \end{aligned}$$
(56)

where γ and β are arbitrary constant. Applying a similar procedure, the series solutions of RPSM are given by

$$\begin{aligned}& u_{1}(t,x)=\frac{t^{\alpha} (-12 \beta^{3} e^{2 \beta x}+3 (4 \beta +\gamma) e^{x (2 \beta+\gamma)}-2 \gamma^{3} e^{\gamma x}+12 \gamma e^{2 \gamma x} )}{2 \varGamma(\alpha+1)}+e^{\gamma x}, \\& v_{1}(t,x)=\frac{1}{2} e^{\beta x} \biggl( \frac{t^{\alpha} (-2 \beta ^{3}+12 (\beta+\gamma) e^{\gamma x}+15 \beta e^{2 \beta x} )}{\varGamma (\alpha+1)}+2 \biggr), \\& u_{2}(t,x) =t^{2\alpha} \bigl( 9 \bigl( 104\beta^{2}+20\beta \gamma +\gamma^{2} \bigr) e^{x(4\beta+\gamma)}\\& \phantom{u_{2}(t,x) =}+4 \bigl( 60\beta ^{6}e^{2\beta x}-774\beta^{4}e^{4\beta x}+\gamma^{6}e^{\gamma x}-60\gamma ^{4}e^{2\gamma x} \bigr) +108\gamma^{2}e^{3\gamma x} \bigr) \\& \phantom{u_{2}(t,x) =}/{\bigl(4\varGamma(2\alpha +1)\bigr)} \\& \phantom{u_{2}(t,x) =}+t^{2\alpha} \bigl( 18(2\beta+\gamma)(4\beta +5\gamma )e^{2x(\beta +\gamma)}\\& \phantom{u_{2}(t,x) =}-3 \bigl( 68\beta^{4}+101\beta^{3}\gamma +72\beta ^{2}\gamma ^{2}+28\beta\gamma^{3}+4\gamma^{4} \bigr) e^{x(2\beta+\gamma )} \bigr) \\& \phantom{u_{2}(t,x) =}/\bigl(4\varGamma(2\alpha+1)\bigr) \\& \phantom{u_{2}(t,x) =}+\frac{t^{\alpha}}{4\varGamma(\alpha +1)}-24\beta^{3}e^{2\beta x}+6(4\beta+ \gamma)e^{x(2\beta+\gamma)}-4\gamma^{3}e^{\gamma x}+e^{\gamma x}+24 \gamma e^{2\gamma x}, \\& v_{2}(t,x) =\frac{t^{2\alpha}e^{\beta x} ( 72(3\beta+\gamma )(7\beta +3\gamma)e^{x(2\beta+\gamma)}-24(\beta+\gamma)^{2} ( 2\beta ^{2}+\beta\gamma+2\gamma^{2} ) e^{\gamma x} ) }{4\varGamma (2\alpha+1)} \\& \phantom{v_{2}(t,x) =}+\frac{t^{2\alpha}e^{\beta x} ( 4\beta ^{6}-1332\beta ^{4}e^{2\beta x}+1125\beta^{2}e^{4\beta x}+144(\beta+2\gamma)^{2}e^{2\gamma x} ) }{4\varGamma(2\alpha+1)} \\& \phantom{v_{2}(t,x) =}+\frac{t^{\alpha}e^{\beta x} ( -4\beta ^{3}+24(\beta+\gamma )e^{\gamma x}+30\beta e^{2\beta x} ) }{4\varGamma(\alpha +1)}+e^{\beta x}. \end{aligned}$$

5 Conclusion

In this paper, we have employed efficient analytical techniques, called the residual power series method and the q-homotopy analysis method, to obtain approximate series solutions to the time-fractional coupled Jaulent–Miodek system of equations which comes with energy-dependent Schrödinger potential with different initial conditions. The numerical results are compared with the known exact solutions when \(\alpha=1\). From our numerical results, see Tables 12 and Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1112, we demonstrate the fast convergence rate of the present methods even after computing a few iterations in solving a system of strongly nonlinear fractional differential equations. Although these methods are different, the results are similar and the resulting errors are comparable. Both methods are elegant and do not require any transformations, perturbations or discretization. Our numerical results further show that the methods are reliable, powerful and easy to implement when compared to other numerical and approximate methods. Since the time-fractional coupled Jaulent–Miodek system of equations is a complex system, the results prove that the present methods could be applied to various complex fractional linear and nonlinear models occurring in various fields of science and engineering, respectively.

Figure 5
figure 5

RPSM vs. q-HAM for Case II for \(U^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=0.35\)

Figure 6
figure 6

RPSM vs. q-HAM for Case II for \(U^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=0.45\)

Figure 7
figure 7

RPSM vs. q-HAM for Case II for \(U^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=0.55\)

Figure 8
figure 8

RPSM vs. q-HAM for Case II for \(U^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=1\)

Figure 9
figure 9

RPSM vs. q-HAM for Case II for \(V^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=0.35\)

Figure 10
figure 10

RPSM vs. q-HAM for Case II for \(V^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=0.45\)

Figure 11
figure 11

RPSM vs. q-HAM for Case II for \(V^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=0.55\)

Figure 12
figure 12

RPSM vs. q-HAM for Case II for \(V^{(2)}(t,x)\) with \(\beta =0.1\), \(\gamma=0.1\) and \(\alpha=1\)

References

  1. Leibniz, G.W.: Letter from Hanover, Germany to G.F.A. L’Hospital, September 30, 1695. In: Leibniz Mathematische Schriften, pp. 301–302. Olms-Verlag, Hildesheim (1962) (first published in 1849)

    Google Scholar 

  2. Leibniz, G.W.: Letter from Hanover, Germany to Johann Bernoulli, December 28, 1695. In: Leibniz Mathematische Schriften, p. 226. Olms-Verlag, Hildesheim (1962) (first published in 1849)

    Google Scholar 

  3. Podlubny, I.: Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications. Mathematics in Science and Engineering, vol. 198. Academic Press, San Diego (1998)

    MATH  Google Scholar 

  4. Ross, B.: The development of fractional calculus 1695–1900. Hist. Math. 4(1), 75–89 (1977)

    Article  MathSciNet  MATH  Google Scholar 

  5. Baleanu, D., Guvenc, Z.B., Machado, J.T.: New Trends in Nanotechnology and Fractional Calculus Applications. Springer, Berlin (2010)

    Book  MATH  Google Scholar 

  6. Chen, D., Chen, Y., Xue, D.: Three fractional-order TV-models for image de-noising. J. Comput. Inf. Syst. 9(12), 4773–4780 (2013)

    Google Scholar 

  7. Ullah, A., Chen, W., Sun, H.G., Khan, M.A.: An efficient variational method for restoring images with combined additive and multiplicative noise. Int. J. Appl. Comput. Math. 3(3), 1999–2019 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  8. Hilfer, R., Anton, L.: Fractional master equations and fractal time random walks. Phys. Rev. E 51, R848–R851 (1995)

    Article  Google Scholar 

  9. Klages, R., Radons, G., Sokolov, I.: Anomalous Transport: Foundations and Applications. Wiley, New York (2008)

    Book  Google Scholar 

  10. Li, Z., Wang, H., Xiao, R., Yang, S.: A variable-order fractional differential equation model of shape memory polymers. Chaos Solitons Fractals 102, 473–485 (2017)

    Article  MathSciNet  Google Scholar 

  11. Mainardi, F.: Fractional Calculus and Waves in Linear Viscoelasticity. Imperial College Press, London (2010)

    Book  MATH  Google Scholar 

  12. Monje, C.A., Chen, Y., Vinagre, B.M., Xue, D., Feliu, V.: Fractional-Order Systems and Controls, Advances in Industrial Control. Springer, Berlin (2010)

    MATH  Google Scholar 

  13. Sibatov, R.T., Svetukhin, V.V.: Subdiffusion kinetics of nanoprecipitate growth and destruction in solid solutions. Theor. Math. Phys. 183(3), 846–859 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  14. Uchaikin, V.V.: Fractional Derivatives for Physicists and Engineers, vol. I. Nonlinear Physical Science. Background and Theory. Higher Education Press, Beijing; Springer, Heidelberg (2013)

    Book  MATH  Google Scholar 

  15. Uchaikin, V.V.: Fractional Derivatives for Physicists and Engineers, vol. II. Nonlinear Physical Science. Applications. Higher Education Press, Beijing; Springer, Heidelberg (2013)

    Book  MATH  Google Scholar 

  16. Tarasov, V.E., Tarasova, V.V.: Time-dependent fractional dynamics with memory in quantum and economic physics. Ann. Phys. 383, 579–599 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  17. Zhang, Y., Benson, D.A., Meerschaert, M.M., LaBolle, E.M., Scheffler, H.P.: Random walk approximation of fractional order multiscaling anomalous diffusion. Phys. Rev. E 74, 026706 (2006)

    Article  Google Scholar 

  18. Zhang, J., Wei, Z., Xiao, L.: Adaptive fractional multiscale method for image de-noising. J. Math. Imaging Vis. 43, 39–49 (2012)

    Article  Google Scholar 

  19. Zhang, Y., Pu, Y.F., Hu, J.R., Zhou, J.L.: A class of fractional-order variational image in-painting models. Appl. Math. Inf. Sci. 6(2), 299–306 (2012)

    MathSciNet  Google Scholar 

  20. Pu, Y.F.: Fractional differential analysis for texture of digital image. J. Algorithms Comput. Technol. 1(3), 357–380 (2007)

    Article  MathSciNet  Google Scholar 

  21. Jaulent, M., Miodek, I.: Nonlinear evolution equations associated with energy-dependent Schrödinger potential. Lett. Math. Phys. 1, 243–250 (1976)

    Article  MATH  Google Scholar 

  22. Atangana, A., Alabaraoye, E.: Solving system of fractional partial differential equations arisen in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller–Segel equation. Adv. Differ. Equ. 2013, Article ID 94 (2013)

    Article  MATH  Google Scholar 

  23. Atangana, A., Kilicman, A.: Analytical solutions of the spacetime fractional derivative of advection dispersion equation. Math. Probl. Eng. 2013, Article ID 853127 (2013)

    MATH  Google Scholar 

  24. Atangana, A., Botha, J.F.: Analytical solution of the ground water flow equation obtained via homotopy decomposition method. J. Earth Sci. Clim. Change 3, Article ID 115 (2012)

    Article  Google Scholar 

  25. Das, G.C., Sarma, J., Uberoi, C.: Explosion of a soliton in a multicomponent plasma. Phys. Plasmas 4(6), 2095–2100 (1997)

    Article  MathSciNet  Google Scholar 

  26. Hong, T., Wang, Y.Z., Huo, Y.S.: Bogoliubov quasiparticles carried by dark solitonic excitations in non-uniform Bose–Einstein condensates. Chin. Phys. Lett. 15, 550–552 (1998)

    Article  Google Scholar 

  27. Ma, W.X., Li, C.X., He, J.: A second Wronskian formulation of the Boussinesq equation. Nonlinear Anal., Theory Methods Appl. 70(12), 4245–4258 (2009)

    Article  MathSciNet  MATH  Google Scholar 

  28. Momani, S., Odibat, Z.: Numerical approach to differential equations of fractional order. J. Comput. Appl. Math. 207(1), 96–110 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  29. Ray, S.S., Bera, R.: Analytical solution of a fractional diffusion equation by Adomian decomposition method. Appl. Math. Comput. 174(1), 329–336 (2006)

    MathSciNet  MATH  Google Scholar 

  30. Abdulaziz, O., Hashim, I., Momani, S.: Application of homotopy-perturbation method to fractional IVPs. J. Comput. Appl. Math. 216(2), 574–584 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  31. Abdulaziz, O., Hashim, I., Momani, S.: Solving systems of fractional differential equations by homotopy-perturbation method. Phys. Lett. A 372, 451–459 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  32. Ganji, Z., Ganji, D., Jafari, H., Rostamian, M.: Application of the homotopy perturbation method to coupled system of partial differential equations with time fractional derivatives. Topol. Methods Nonlinear Anal. 31(2), 341–348 (2008)

    MathSciNet  MATH  Google Scholar 

  33. Hosseinnia, S.H., Ranjbar, A., Momani, S.: Using an enhanced homotopy perturbation method in fractional equations via deforming the linear part. Comput. Math. Appl. 56, 3138–3149 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  34. Yildirim, A., Gulkanat, Y.: Analytical approach to fractional Zakharov–Kuznetsov equations by He’s homotopy perturbation method. Commun. Theor. Phys. 53(6), 1005–1010 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  35. Neamaty, A., Agheli, B., Darzi, R.: Variational iteration method and He’s polynomials for time-fractional partial differential equations. Prog. Fract. Differ. Appl. 1(1), 47–55 (2015)

    MATH  Google Scholar 

  36. Veeresha, P., Prakasha, D.G., Baskonus, H.M.: Solving smoking epidemic model of fractional order using a modified homotopy analysis transform method. Math. Sci. 13(2), 115–128 (2019)

    Article  MathSciNet  Google Scholar 

  37. Veeresha, P., Prakasha, D.G.: Solution for fractional Zakharov–Kuznetsov equations by using two reliable techniques. Chin. J. Phys. 60, 313–330 (2019)

    Article  MathSciNet  Google Scholar 

  38. Prakasha, D.G., Veeresha, P., Rawashdeh, M.S.: Numerical solution for \((2 + 1)\)-dimensional time-fractional coupled Burger equations using fractional natural decomposition method. Math. Methods Appl. Sci. 42(10), 3409–3427 (2019)

    Article  MathSciNet  MATH  Google Scholar 

  39. Momani, S., Abu Arqub, O., Freihat, A., Al-Smadi, M.: Analytical approximations for Fokker–Planck equations of fractional order in multistep schemes. Appl. Comput. Math. 15(3), 319–330 (2016)

    MathSciNet  MATH  Google Scholar 

  40. Bhalekar, S., Daftardar-Gejji, V.: New iterative method: application to partial differential equations. Appl. Math. Comput. 203(2), 778–783 (2008)

    MathSciNet  MATH  Google Scholar 

  41. Daftardar-Gejji, V., Bhalekar, S.: Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method. Comput. Math. Appl. 59(5), 1801–1809 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  42. Akinyemi, L., Iyiola, O.S., Akpan, U.: Iterative methods for solving fourth and sixth order time-fractional Cahn–Hillard equation (2019). arXiv:1903.10337

  43. Abu Arqub, O., Abo-Hammour, Z., Al-Badarneh, R., Momani, S.: A reliable analytical method for solving higher-order initial value problems. Discrete Dyn. Nat. Soc. 2013, Article ID 673829 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  44. Abdulaziz, O., Hashim, I., Saif, A.: Series solutions of time-fractional PDEs by homotopy analysis method. Int. J. Differ. Equ. 2008, Article ID 686512 (2008)

    MathSciNet  MATH  Google Scholar 

  45. Rashidi, M.M., Domairry, G., Dinarvand, S.: The homotopy analysis method for explicit analytical solutions of Jaulent–Miodek equations. Numer. Methods Partial Differ. Equ. 25(2), 430–439 (2009)

    Article  MathSciNet  MATH  Google Scholar 

  46. Abbasbandy, S., Shirzadi, A.: Homotopy analysis method for multiple solutions of the fractional Sturm–Liouville problems. Numer. Algorithms 54(4), 521–532 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  47. El-Tawil, M.A., Huseen, S.N.: The Q-homotopy analysis method (qHAM). Int. J. Appl. Math. Mech. 8(15), 51–75 (2012)

    Google Scholar 

  48. Iyiola, O.S.: Exact and approximate solutions of fractional diffusion equations with fractional reaction terms. Prog. Fract. Differ. Appl. 2(1), 21–30 (2016)

    Article  Google Scholar 

  49. Iyiola, O.S.: On the solutions of nonlinear time-fractional gas dynamic equations: an analytical approach. Int. J. Pure Appl. Math. 98(4), 491–502 (2015)

    Article  Google Scholar 

  50. Iyiola, O.S., Zaman, F.D.: A fractional diffusion equation model for cancer tumor. AIP Adv. 4, 107121 (2014)

    Article  Google Scholar 

  51. Iyiola, O.S., Ojo, G.O.: Analytical solutions of time-fractional models for homogeneous Gardner equation and nonhomogeneous differential equations. Ain Shams Eng. J. 5, 999–1004 (2014)

    Article  Google Scholar 

  52. Iyiola, O.S., Soh, M.E., Enyi, C.D.: Generalised homotopy analysis method (qHAM) for solving foam drainage equation of time fractional type. Math. Eng. Sci. Aerosp. 4(4), 429–440 (2013)

    MATH  Google Scholar 

  53. Akinyemi, L.: Q-homotopy analysis method for solving seventh-order time-fractional Lax’s Korteweg–de Vries and Sawada–Kotera equations. Comput. Appl. Math. 38, Article ID 191 (2019). https://doi.org/10.1007/s40314-019-0977-3

    Article  MathSciNet  MATH  Google Scholar 

  54. Atangana, A., Baleanu, D.: Nonlinear fractional Jaulent–Miodek and Whitham–Broer–Kaup equations with Sumudu transform. Abstr. Appl. Anal. 2013, Article ID 160681 (2013)

    MathSciNet  MATH  Google Scholar 

  55. Lou, S.Y.: A direct perturbation method: nonlinear Schrödinger equation with loss. Chin. Phys. Lett. 16, 659–661 (1999)

    Article  Google Scholar 

  56. Ozer, H.T., Salihoglu, S.: Nonlinear Schrödinger equations and \(N=1\) super-conformal algebra. Chaos Solitons Fractals 33, 1417–1423 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  57. Gupta, A.K., Ray, S.S.: An investigation with Hermite wavelets for accurate solution of fractional Jaulent–Miodek equation associated with energy-dependent Schrödinger potential. Appl. Math. Comput. 270, 458–471 (2015)

    MathSciNet  MATH  Google Scholar 

  58. Majlesi, A., Ghehsareha, H.R., Zaghian, A.: On the fractional Jaulent–Miodek equation associated with energy-dependent Schrodinger potential: Lie symmetry reductions, explicit exact solutions and conservation laws. Eur. Phys. J. Plus 132, Article ID 516 (2017)

    Article  Google Scholar 

  59. Veeresha, P., Prakasha, D.G., Magesh, N., Nandeppanavar, M.M., Christopher, A.J.: Numerical simulation for fractional Jaulent–Miodek equation associated with energy-dependent Schrödinger potential using two novel techniques (2019). arXiv:1810.06311

  60. He, J.H., Zhang, L.N.: Generalized solitary solution and compacton-like solution of the Jaulent–Miodek equations using the Exp-function method. Phys. Lett. A 372(7), 1044–1047 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  61. Yildirim, A., Kelleci, A.: Numerical simulation of the Jaulent–Miodek equation by He’s homotopy perturbation method. World Appl. Sci. J. 7, 84–89 (2009)

    Google Scholar 

  62. Rashidi, M.M., Domairry, G., Dinarvand, S.: The homotopy analysis method for explicit analytical solutions of Jaulent–Miodek equations. Numer. Methods Partial Differ. Equ. 25(2), 430–439 (2009)

    Article  MathSciNet  MATH  Google Scholar 

  63. Luchko, Y.F., Srivastava, H.: The exact solution of certain differential equations of fractional order by using operational calculus. Comput. Math. Appl. 29(8), 73–85 (1995)

    Article  MathSciNet  MATH  Google Scholar 

  64. El-Ajou, A., Abu Arqub, O., Momani, S.: Approximate analytical solution of the nonlinear fractional KdV–Burgers equation: a new iterative algorithm. J. Comput. Phys. 293, 81–95 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  65. Alquran, M.: Analytical solution of time-fractional two-component evolutionary system of order 2 by residual power series method. J. Appl. Anal. Comput. 5(4) 589–599 (2015)

    MathSciNet  Google Scholar 

  66. Alquran, M., Jaradat, H.M., Syam, M.I.: Analytical solution of the time-fractional Phi-4 equation by using modified residual power series method. Nonlinear Dyn. 90(4), 2525–2529 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  67. Freihet, A., Hasan, S., Al-Smadi, M., Gaith, M., Momani, S.: Construction of fractional power series solutions to fractional stiff system using residual functions algorithm. Adv. Differ. Equ. 2019, Article ID 95 (2019)

    Article  MathSciNet  MATH  Google Scholar 

  68. Hasan, S., Al-Smadi, M., Freihet, A., Momani, S.: Two computational approaches for solving a fractional obstacle system in Hilbert space. Adv. Differ. Equ. 2019(1), 55 (2019)

    Article  MathSciNet  MATH  Google Scholar 

  69. Prakasha, D.G., Veeresha, P., Baskonus, H.M.: Residual power series method for fractional Swift–Hohenberg equation. Fractal Fract. 3(1), Article ID 9 (2019)

    Article  Google Scholar 

  70. Liao, S.J.: An approximate solution technique not depending on small parameters: a special example. Int. J. Non-Linear Mech. 30(3), 371–380 (1995)

    Article  MathSciNet  MATH  Google Scholar 

  71. Abu Arqub, O.: Series solution of fuzzy differential equations under strongly generalized differentiability. J. Adv. Res. Appl. Math. 5(1), 31–52 (2013)

    Article  MathSciNet  Google Scholar 

  72. Abu Arqub, O., El-Ajou, A., Bataineh, A.S., Hashim, I.: A representation of the exact solution of generalized Lane–Emden equations using a new analytical method. Abstr. Appl. Anal. 2013, Article ID 378593 (2013)

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

The authors thank the referees for a number of suggestions which have improved many aspects of this article.

Availability of data and materials

Not applicable.

Funding

No funding available for this project.

Author information

Authors and Affiliations

Authors

Contributions

The authors read and approved the final manuscript.

Corresponding author

Correspondence to Olaniyi S. Iyiola.

Ethics declarations

Competing interests

The authors declare to have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix

Appendix

The components of the series solutions obtained by qHAM given in Eqs. (52) and (53) are obtained for Case I and Case II, successively below.

1.1 Case I

u 1 ( t , x ) = χ 1 u 0 ( t , x ) + h J α [ R 1 ( u 0 ( t , x ) ) ] = h J α [ D t α u 0 + u 0 x x x + 1.5 v 0 v 0 x x x + 4.5 v 0 x v 0 x x ] h J α [ 6 u 0 u 0 x + 6 u 0 v 0 v 0 x + 1.5 v 0 v 0 u 0 x ] = h λ 5 tanh ( λ x 2 ) sech 2 ( λ x 2 ) 4 Γ ( α + 1 ) t α , v 1 ( t , x ) = χ 1 v 0 ( t , x ) + h J α [ R ˆ 1 ( v 0 ( t , x ) ) ] = h J α [ D t α v 0 + v 0 x x x 6 v 0 u 0 x 6 u 0 v 0 x 7.5 v 0 v 0 v 0 x ] = h λ 4 tanh ( λ x 2 ) sech ( λ x 2 ) 4 Γ ( α + 1 ) t α , u 2 ( t , x ) = χ 2 u 1 ( t , x ) + h J α [ R 2 ( u 1 ( t , x ) ) ] u 2 ( t , x ) = n u 1 + h J α [ D t α u 1 + u 1 x x x + 1.5 v 0 v 1 x x x + 1.5 v 1 v 0 x x x + 4.5 v 0 x v 1 x x + 4.5 v 1 x v 0 x x ] u 2 ( t , x ) = h J α [ 6 u 0 u 1 x + 6 u 1 u 0 x + 6 u 0 v 0 v 1 x + 6 u 0 v 1 v 0 x + 6 u 1 v 0 v 0 x u 2 ( t , x ) = + 1.5 v 0 v 0 u 1 x + 1.5 v 0 v 1 u 0 x + 1.5 v 1 v 0 u 0 x ] u 2 ( t , x ) = ( n + h ) u 1 + h J α [ u 1 x x x + 1.5 v 0 v 1 x x x + 1.5 v 1 v 0 x x x + 4.5 v 0 x v 1 x x + 4.5 v 1 x v 0 x x ] u 2 ( t , x ) = h J α [ 6 u 0 u 1 x + 6 u 1 u 0 x + 6 u 0 v 0 v 1 x + 6 u 0 v 1 v 0 x + 6 u 1 v 0 v 0 x u 2 ( t , x ) = + 1.5 v 0 v 0 u 1 x + 1.5 v 0 v 1 u 0 x + 1.5 v 1 v 0 u 0 x ] u 2 ( t , x ) = h λ 5 ( n + h ) tanh ( λ x 2 ) sech 2 ( λ x 2 ) 4 Γ ( α + 1 ) t α h 2 λ 8 ( cosh ( λ x ) 2 ) sech 4 ( λ x 2 ) 16 Γ ( 2 α + 1 ) t 2 α , v 2 ( t , x ) = χ 2 v 1 ( t , x ) + h J α [ R ˆ 2 ( v 1 ( t , x ) ) ] = n v 1 + h J α [ D t α v 1 + v 1 x x x 6 v 0 u 1 x 6 v 1 u 0 x 6 u 0 v 1 x 6 u 1 v 0 x 7.5 v 0 v 0 v 1 x 7.5 v 0 v 1 v 0 x 7.5 v 1 v 0 v 0 x ] = ( n + h ) v 1 + h J α [ v 1 x x x 6 v 0 u 1 x 6 v 1 u 0 x 6 u 0 v 1 x 6 u 1 v 0 x 7.5 v 0 v 0 v 1 x 7.5 v 0 v 1 v 0 x 7.5 v 1 v 0 v 0 x ] = h λ 4 ( h + n ) tanh ( λ x 2 ) sech ( λ x 2 ) 4 Γ ( α + 1 ) t α + λ 7 h 2 ( cosh ( λ x ) 3 ) sech 3 ( λ x 2 ) 32 Γ ( 2 α + 1 ) t 2 α , u 3 ( t , x ) = χ 3 u 2 ( t , x ) + h J α [ R 3 ( u 2 ( t , x ) ) ] = n u 2 + h J α [ D t α u 2 + u 2 x x x + 1.5 v 0 v 2 x x x + 1.5 v 1 v 1 x x x + 1.5 v 2 v 0 x x x + 4.5 v 0 x v 2 x x + 4.5 v 1 x v 1 x x + 4.5 v 2 x v 0 x x ] h J α [ 6 u 0 u 2 x + 6 u 1 u 1 x + 6 u 2 u 0 x + 6 u 0 v 0 v 2 x + 6 u 0 v 1 v 1 x + 6 u 1 v 0 v 1 x + 6 u 0 v 2 v 0 x + 6 u 1 v 1 v 0 x + 6 u 2 v 0 v 0 x ] h J α [ 1.5 v 0 v 0 u 2 x + 1.5 v 0 v 1 u 1 x + 1.5 v 1 v 0 u 1 x + 1.5 v 0 v 2 u 0 x + 1.5 v 1 v 1 u 0 x + 1.5 v 2 v 0 u 0 x ] = ( n + h ) u 2 + h J α [ u 2 x x x + 1.5 v 0 v 2 x x x + 1.5 v 1 v 1 x x x + 1.5 v 2 v 0 x x x + 4.5 v 0 x v 2 x x + 4.5 v 1 x v 1 x x + 4.5 v 2 x v 0 x x ] h J α [ 1.5 v 0 v 0 u 2 x + 1.5 v 1 v 1 u 1 x + 1.5 v 1 v 0 u 1 x + 1.5 v 0 v 2 u 0 x + 1.5 v 1 v 1 u 0 x + 1.5 v 2 v 0 u 0 x ] = h λ 5 ( h + n ) 2 sinh ( λ x ) sech 4 ( λ x 2 ) 8 Γ ( α + 1 ) t α h 2 λ 8 ( h + n ) ( cosh ( λ x ) 2 ) sech 4 ( λ x 2 ) 8 Γ ( 2 α + 1 ) t 2 α + [ ( h 3 λ 11 [ 3 Γ ( 2 α + 1 ) ( 20 cosh ( λ x ) cosh ( 2 λ x ) 43 ) 2 Γ ( α + 1 ) 2 ( cosh ( 2 λ x ) + 28 cosh ( λ x ) 165 ) ] tanh ( λ x 2 ) sech 6 ( λ x 2 ) ) / ( 1024 Γ ( α + 1 ) 2 Γ ( 3 α + 1 ) ) ] t 3 α , v 3 ( t , x ) = χ 3 v 2 ( t , x ) + h J α [ R ˆ 3 ( v 2 ( t , x ) ) ] v 3 ( t , x ) = n v 2 + h J α [ D t α v 2 + v 2 x x x 6 v 0 u 2 x 6 v 1 u 1 x 6 v 2 u 0 x v 3 ( t , x ) = 6 u 0 x v 2 x 6 u 1 x v 1 x 6 u 2 x v 0 x ] v 3 ( t , x ) = h J α [ 7.5 v 0 v 0 v 2 x + 7.5 v 0 v 1 v 1 x + 7.5 v 1 v 0 v 1 x v 3 ( t , x ) = + 7.5 v 0 v 2 v 0 x + 7.5 v 1 v 1 v 0 x + 7.5 v 2 v 0 v 0 x ] v 3 ( t , x ) = ( n + h ) v 2 + h J α [ v 2 x x x 6 v 0 u 2 x 6 v 1 u 1 x 6 v 2 u 0 x v 3 ( t , x ) = 6 u 0 x v 2 x 6 u 1 x v 1 x 6 u 2 x v 0 x ] v 3 ( t , x ) = h J α [ 7.5 v 0 v 0 v 2 x + 7.5 v 0 v 1 v 1 x + 7.5 v 1 v 0 v 1 x v 3 ( t , x ) = + 7.5 v 0 v 2 v 0 x + 7.5 v 1 v 1 v 0 x + 7.5 v 2 v 0 v 0 x ] v 3 ( t , x ) = h λ 4 ( h + n ) 2 tanh ( λ x 2 ) sech ( λ x 2 ) 4 Γ ( α + 1 ) t α + λ 7 h 2 ( h + n ) ( cosh ( λ x ) 3 ) sech 3 ( λ x 2 ) 16 Γ ( 2 α + 1 ) t 2 α v 3 ( t , x ) = + [ ( h 3 λ 10 [ 12 Γ ( 2 α + 1 ) ( 3 cosh ( λ x ) 7 ) v 3 ( t , x ) = + Γ ( α + 1 ) 2 ( cosh ( 2 λ x ) 92 cosh ( λ x ) + 147 ) ] tanh ( λ x 2 ) sech 5 ( λ x 2 ) ) v 3 ( t , x ) = / ( 512 Γ ( α + 1 ) 2 Γ ( 3 α + 1 ) ) ] t 3 α .

1.2 Case II

$$\begin{aligned}& \begin{aligned} u_{1}(t,x)={}&\chi^{*}_{1} u_{0}(t,x)+h J^{\alpha}\bigl[\mathcal{R}_{1} \bigl(\vec{u}_{0}(t,x) \bigr) \bigr] \\ ={}&h J^{\alpha}\bigl[\mathcal{D}_{t}^{\alpha}u_{0} + u_{0xxx}+1.5v_{0}v_{0xxx}+4.5v_{0x}v_{0xx} \bigr] \\ &-h J^{\alpha}[6u_{0}u_{0x}+6u_{0}v_{0}v_{0x}+1.5v_{0}v_{0}u_{0x} ] \\ ={}&{-}\frac{h[3(\gamma+4\beta)e^{(\gamma+2\beta)x}-2(\gamma ^{3}e^{\gamma x}-6\gamma e^{2\gamma x}+6\beta^{3}e^{2\beta x})]}{2\varGamma (\alpha+ 1)} t^{\alpha}, \end{aligned} \\& \begin{aligned}v_{1}(t,x)&=\chi^{*}_{1} v_{0}(t,x)+h J^{\alpha}\bigl[\hat{\mathcal{R}}_{1} \bigl( \vec{v}_{0}(t,x) \bigr) \bigr] \\ &=h J^{\alpha}\bigl[\mathcal{D}_{t}^{\alpha}v_{0}+ v_{0xxx}-6v_{0}u_{0x}-6u_{0}v_{0x}-7.5v_{0}v_{0}v_{0x} \bigr] \\ &=-\frac{he^{\beta x}[12(\gamma+\beta)e^{\gamma x}-2\beta ^{3}+15\beta e^{2\beta x}]}{ 2\varGamma(\alpha+1)} t^{\alpha}, \end{aligned} \\& u_{2}(t,x)=\chi^{*}_{2} u_{1}(t,x)+h J^{\alpha}\bigl[\mathcal{R}_{2} \bigl(\vec{u}_{1}(t,x) \bigr) \bigr] \\& \phantom{u_{2}(t,x)}=n u_{1}+h J^{\alpha}\bigl[\mathcal{D}_{t}^{\alpha}u_{1} + u_{1xxx}+1.5v_{0}v_{1xxx}+1.5v_{1}v_{0xxx}+4.5v_{0x}v_{1xx}+4.5v_{1x}v_{0xx} \bigr] \\& \phantom{u_{2}(t,x)=}-h J^{\alpha}[6u_{0}u_{1x}+6u_{1}u_{0x}+6u_{0}v_{0}v_{1x}+6u_{1}v_{1}v_{0x}+6u_{1}v_{0}v_{0x}\\& \phantom{u_{2}(t,x)=}+1.5v_{0}v_{0}u_{1x}+1.5v_{1}v_{1}u_{0x}+1.5v_{1}v_{0}u_{0x} ] \\& \phantom{u_{2}(t,x)}=(n+h) u_{1}+h J^{\alpha}[u_{1xxx}+1.5v_{0}v_{1xxx}+1.5v_{1}v_{0xxx}+4.5v_{0x}v_{1xx}+4.5v_{1x}v_{0xx} ] \\& \phantom{u_{2}(t,x)=}-h J^{\alpha}[6u_{0}u_{1x}+6u_{1}u_{0x}+6u_{0}v_{0}v_{1x}+6u_{1}v_{1}v_{0x}+6u_{1}v_{0}v_{0x}\\& \phantom{u_{2}(t,x)=}+1.5v_{0}v_{0}u_{1x}+1.5v_{1}v_{1}u_{0x}+1.5v_{1}v_{0}u_{0x} ] \\& \phantom{u_{2}(t,x)=}=-\frac{h(n+h) [3(\gamma+4\beta)e^{(\gamma +2\beta )x}-2(\gamma^{3}e^{\gamma x}-6\gamma e^{2\gamma x}+6\beta^{3}e^{2\beta x}) ]}{2\varGamma(\alpha+ 1)} t^{\alpha}\\& \phantom{u_{2}(t,x)=}-\Bigl(3h^{2}e^{(\gamma+2\beta)x} \bigl[4\gamma ^{4}+28\gamma^{3}\beta +72\gamma^{2}\beta^{2}+101\gamma\beta^{3}-84\gamma\beta e^{(\gamma x)}\\& \phantom{u_{2}(t,x)=}-15\gamma\beta e^{(2\beta x)}+68\beta^{4}-78\beta ^{2}e^{(2\beta x)} \bigr]\\& \phantom{u_{2}(t,x)=}/\bigl(\varGamma(2\alpha+ 1)\bigr) \Bigr)t^{2\alpha} \\& \phantom{u_{2}(t,x)=}+\Bigl(h^{2} \bigl[4\gamma^{6}e^{(\gamma x)}-240\gamma^{4}e^{(2\gamma x)}+9\gamma^{2} \bigl(40e^{2(\gamma+\beta)x}+e^{(\gamma+4\beta )x}+48e^{3\gamma x} \bigr)\\& \phantom{u_{2}(t,x)=}+24\beta^{2} e^{2\beta x } \bigl(24e^{2\gamma x}+10\beta^{4}-129\beta^{2}e^{2\beta x } \bigr) \bigr]\\& \phantom{u_{2}(t,x)=}/\bigl(4\varGamma(2\alpha+ 1)\bigr)\Bigr) t^{2\alpha}, \\& \begin{aligned} v_{2}(t,x)={}&\chi^{*}_{2} v_{1}(t,x)+h J^{\alpha}\bigl[\hat{\mathcal{R}}_{2} \bigl( \vec{v}_{1}(t,x) \bigr) \bigr] \\ ={}&n v_{1}+h J^{\alpha}\bigl[\mathcal{D}_{t}^{\alpha}v_{1}+ v_{1xxx}-6v_{0}u_{1x}-6v_{1}u_{0x}-6u_{0}v_{1x}-6u_{1}v_{0x}\\ &-7.5v_{0}v_{0}v_{1x}-7.5v_{1}v_{1}v_{0x}-7.5v_{1}v_{0}v_{0x} \bigr] \\ ={}&(n+h) v_{1}+hJ^{\alpha}[v_{1xxx}-6v_{0}u_{1x}-6v_{1}u_{0x}-6u_{0}v_{1x}-6u_{1}v_{0x}\\ &-7.5v_{0}v_{0}v_{1x}-7.5v_{1}v_{1}v_{0x}-7.5v_{1}v_{0}v_{0x} ] \\ ={}&{-}\frac{h(n+h)e^{\beta x}[12(\gamma+\beta)e^{\gamma x}-2\beta ^{3}+15\beta e^{2\beta x}]}{ 2\varGamma(\alpha+1)} t^{\alpha}\\ &+\frac{h^{2}e^{\beta x}[-48e^{\gamma x}(\gamma^{4}+\beta^{4})+576\gamma ^{2}e^{2\gamma x}+4\beta^{6}-1332\beta^{4}e^{2\beta x}+1125\beta ^{2}e^{4\beta x}] }{4\varGamma(2\alpha+1)} t^{2\alpha} \\ &+\Bigl(6h^{2}e^{(\gamma+\beta)x}\bigl[6\beta e^{\gamma x}(4\gamma +\beta )-\gamma{\beta}\bigl(5\gamma^{2}+6\gamma{\beta}+5\beta^{2}\bigr)\\ &+3e^{2\beta x}(\gamma+3\beta)(3\gamma+7\beta)\bigr]\\ &/\varGamma(2\alpha+1) \Bigr)t^{2\alpha}. \end{aligned} \end{aligned}$$

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Şenol, M., Iyiola, O.S., Daei Kasmaei, H. et al. Efficient analytical techniques for solving time-fractional nonlinear coupled Jaulent–Miodek system with energy-dependent Schrödinger potential. Adv Differ Equ 2019, 462 (2019). https://doi.org/10.1186/s13662-019-2397-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13662-019-2397-5

Keywords