On stability analysis of a class of three-dimensional system of exponential difference equations

Abdul Khaliq; Haza Saleh Alayachi; Muhammad Zubair; Muhammad Rohail; Abdul Qadeer Khan; Abdul Khaliq; Haza Saleh Alayachi; Muhammad Zubair; Muhammad Rohail; Abdul Qadeer Khan

doi:10.3934/math.2023251

AIMS Mathematics

2023, Volume 8, Issue 2: 5016-5035. doi: 10.3934/math.2023251

Previous Article Next Article

Research article

On stability analysis of a class of three-dimensional system of exponential difference equations

1.
Department of Mathematics, Riphah International university, Lahore, Pakistan
2.
Department of Mathematics, College of Science, Taibah University, Al-Madinah Al-Munawarah, Saudi Arabia
3.
Government Education Department, Government Elementary School, Walton Road, Lahore, Pakistan
4.
Department of Mathematics, Riphah International university, Lahore, Pakistan
5.
Department of Mathematics, University of Azad Jammu and Kashmir, Muzaffarabad 13100, Pakistan

Received: 08 March 2022 Revised: 22 September 2022 Accepted: 26 September 2022 Published: 12 December 2022
MSC : 39A20, 39A10, 40A05

The boundedness character, persistent nature, and asymptotic conduct of non-negative outcomes of the system of three dimensional exponential form of difference equations were studied in this research:

$\begin{eqnarray*} x_{n+1} & = &ax_{n}+by_{n-1}e^{-x_{n}}, \\ \text{ }y_{n+1} & = &cy_{n}+dz_{n-1}e^{-y_{n}},\ \\ z_{n+1} & = &ez_{n}+fx_{n-1}e^{-z_{n}}, \end{eqnarray*}$

where $a, \ b, \ c$ , $d, \ e$ and $f$ are non-negative real values, and the initial values $x_{-1}, \ x_{0}, \ y_{-1}, \ y_{0}, \ z_{-1}, \ z_{0}$ are non-negative real values.

Keywords:

difference equations in exponential form,
stability character,
periodicity,
rate of convergence

Citation: Abdul Khaliq, Haza Saleh Alayachi, Muhammad Zubair, Muhammad Rohail, Abdul Qadeer Khan. On stability analysis of a class of three-dimensional system of exponential difference equations[J]. AIMS Mathematics, 2023, 8(2): 5016-5035. doi: 10.3934/math.2023251

Related Papers:

[1]	Lv Zhang, Qingbiao Wu . Modified Newton-EHS method for solving nonlinear problems with complex symmetric Jacobian matrices. AIMS Mathematics, 2023, 8(10): 24233-24253. doi: 10.3934/math.20231236
[2]	Lu-Chuan Ceng, Shih-Hsin Chen, Yeong-Cheng Liou, Tzu-Chien Yin . Modified inertial subgradient extragradient algorithms for generalized equilibria systems with constraints of variational inequalities and fixed points. AIMS Mathematics, 2024, 9(6): 13819-13842. doi: 10.3934/math.2024672
[3]	Xiao Guo, Chuanpei Xu, Zhibin Zhu, Benxin Zhang . Nonmonotone variable metric Barzilai-Borwein method for composite minimization problem. AIMS Mathematics, 2024, 9(6): 16335-16353. doi: 10.3934/math.2024791
[4]	Fouad Mohammad Salama, Nur Nadiah Abd Hamid, Umair Ali, Norhashidah Hj. Mohd Ali . Fast hybrid explicit group methods for solving 2D fractional advection-diffusion equation. AIMS Mathematics, 2022, 7(9): 15854-15880. doi: 10.3934/math.2022868
[5]	Mashail M. Al Sobhi . The modified Kies-Fréchet distribution: Properties, inference and application. AIMS Mathematics, 2021, 6(5): 4691-4714. doi: 10.3934/math.2021276
[6]	Zui-Cha Deng, Fan-Li Liu, Liu Yang . Numerical simulations for initial value inversion problem in a two-dimensional degenerate parabolic equation. AIMS Mathematics, 2021, 6(4): 3080-3104. doi: 10.3934/math.2021187
[7]	Xin-Hui Shao, Wan-Chen Zhao . Relaxed modified Newton-based iteration method for generalized absolute value equations. AIMS Mathematics, 2023, 8(2): 4714-4725. doi: 10.3934/math.2023233
[8]	Krunal B. Kachhia, Jyotindra C. Prajapati . Generalized iterative method for the solution of linear and nonlinear fractional differential equations with composite fractional derivative operator. AIMS Mathematics, 2020, 5(4): 2888-2898. doi: 10.3934/math.2020186
[9]	Huiqin Zhang, Yanping Chen . Two-grid finite element method with an H2N2 interpolation for two-dimensional nonlinear fractional multi-term mixed sub-diffusion and diffusion wave equation. AIMS Mathematics, 2024, 9(1): 160-177. doi: 10.3934/math.2024010
[10]	Cui-Xia Li, Long-Quan Yong . Modified BAS iteration method for absolute value equation. AIMS Mathematics, 2022, 7(1): 606-616. doi: 10.3934/math.2022038

Abstract

$\begin{eqnarray*} x_{n+1} & = &ax_{n}+by_{n-1}e^{-x_{n}}, \\ \text{ }y_{n+1} & = &cy_{n}+dz_{n-1}e^{-y_{n}},\ \\ z_{n+1} & = &ez_{n}+fx_{n-1}e^{-z_{n}}, \end{eqnarray*}$

where $a, \ b, \ c$ , $d, \ e$ and $f$ are non-negative real values, and the initial values $x_{-1}, \ x_{0}, \ y_{-1}, \ y_{0}, \ z_{-1}, \ z_{0}$ are non-negative real values.

1. Introduction

Fractional-order calculus has received considerable attention in the engineering and physical sciences over the last few decades to model a number of diverse phenomena in robotic-technology, bio-engineering, control theory, viscoelasticity diffusion model, relaxation processes and signal processing ^[1,2]. The order of derivatives, as well as integrals in the fractional-order calculus, is arbitrary. Therefore, fractional-order NPDEs have developed a fundamental interest in generalising integer-order NPDEs to model complex systems in thermodynamics, engineering, fluid dynamics and optical physics ^[3].

The enormous advantage of using fractional differential equations (FDEs) in modeling real-world problems is their global behavior together with preserving memory ^[4] which is not present in integer-order differential equations. It has also been noted that FDEs fastly converge to ordinary differential equations (ODEs) in a case when fractional-order is equal to one. Moreover, fractional calculus can clarify the basic features of various models and processes them more precisely than integer-order ^[5]. Several techniques have been applied to study analytical as well as numerical solutions of FNPDEs such as, Variational Iteration Method (VIM) ^[6], Laplace transforms Method ^[7], the double Laplace transform ^[8], invariant subspace method ^[9], Integral transform ^[10], Sumudu transform Method (STM) ^[11], natural transform ^[12] and Adomian decomposition method (ADM) ^[13].

The Klein-Gordon equation (KGE) considered herein is a basic non-linear evolution equation that arises in relativistic quantum Mechanics. It was formulated by Erwin Schrödinger for the non-relativistic wave equation in quantum physics, while precisely studied by the famous physicists O. Klein and W. Gordon (as it is named after their work) in 1926 ^[14,15]. The KGE has an extensive variety of applications in classical field theory ^[16] as well as in quantum field theory ^[17]. It has also been extensively used in numerous areas of physical phenomena such as in solid-state physics, dispersive wave-phenomena, nonlinear optics, elementary particles behavior, dislocations propagation in crystals, and different class of soliton solutions ^[18]. Here, we investigate equation of the form ^[19]

$\begin{eqnarray} \frac{\partial^\omega \psi}{\partial t^\omega}-\frac{\partial^{2} \psi}{\partial x^2}+p\psi+qg(\psi) = r(x,t),\ \ \, 1 < \omega\leq 2, \end{eqnarray}$

(1.1)

together with

$\begin{eqnarray*} \psi(x,0) = \mathcal {F}(x),\, \, \psi_t(x,0) = \mathcal {G}(x), \end{eqnarray*}$

where $\psi = \psi(x, t)$ , $g(\psi)$ and $r(x, t)$ represent nonlinear term and external function respectively.

The nonlinear differential equations involve numerous fractional differential operators, such as, Caputo, Hilfer, Riemann-Liouville (R-L), Atangana-Baleanu in Caputo's sense, and Caputo-Fabrizio, ^[20]. The above fractional operators are very useful in FC due to the complexities of fractional-PDEs/ODEs because standard operators cannot handle some equations to obtain explicit solutions. The Caputo fractional derivative is the basic idea of fractional derivatives. All the fractional derivatives will reduce in Caputo or Riemann-Liouville fractional derivatives after some parametric replacement. One can assume that the fractional derivative could provide a power-law of the local behavior of non-differentiable functions. The Caputo fractional derivative was introduced by Michele Caputo in 1967 ^[21] to study initial/boundary value problems in many areas of real-world phenomena. The Caputo's derivative has many advantages as it is the most important tool for dealing with integer order models in a fractional sense with suitable subsidiary conditions ^[23]. Most of the problems have been handled precisely using Caputo operator ^[24].

The integer order KGE has broadly studied by using a variety of methods ^[25]. Time-fractional Klein Gordon equations with Caputo's fractional operator have also been extensively studied using a variety of numerical and analytical techniques ^[26]. Here, we apply double Laplace transform with decomposition method to study the general solutions of the governing model with power law. Some particular examples are also studied numerically with some physical analysis. For preliminaries and some basic definitions of Caputo's derivative, see ^[31] and the reference therein.

Definition 1. Let us suppose $\psi(x, t)$ lies in $x-t\, plane$ , the double Laplace transform (DLT) of $\psi(x, t)$ is defined by ^[32]

$\begin{eqnarray*} \mathcal{L}_x\mathcal{L}_t\left[\psi(x,t)\right] = \int_0^{\infty}e^{-rx}\int_0^{\infty} e^{-s{t}}\,\psi \,dx \,dt, \label{eq2.3} \end{eqnarray*}$

where, $r, s \in \mathcal(C)$ .

Definition 2. Application of DLT on fractional-order operator in Caputo's sense gives

$\begin{equation*} \mathcal{L}_x\mathcal{L}_t\left\{^{C}D_x^\omega\,\psi(x,t)\right\} = r^\omega\,\overline{\psi}(r,s)-\sum\limits_{k = 0}^{n-1}r^{\omega-1-k} \mathcal{L}_t\left\{\frac{\partial^k{\psi(0,t)}}{\partial{x^k}} \right\}, \label{eq2.4} \end{equation*}$

and

$\begin{equation*} \mathcal{L}_x\mathcal{L}_t\left\{^{C}D_t^\beta\,\psi(x,t)\right\} = s^\beta\overline{\psi}(p,s)- \sum\limits_{k = 0}^{m-1}s^{\beta-1-k} \mathcal{L}_{x}\Big\{ \frac{\partial^k{\psi(x,0)}}{\partial{t^{k}}} \Big\}, \label{eq2.5} \end{equation*}$

where, $m = [\beta]+1$ and $n = [\omega]+1$ . Hence, we infer that

$\begin{equation*} \mathcal{L}_x\mathcal{L}_t{\psi(x)\,v(t)} = \overline{\psi}(p)\overline{v}(s) = \mathcal{L}_x{\psi(x)}\mathcal{L}_t{v(t)}. \label{eq2.6} \end{equation*}$

The inverse DLT $\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\{\overline{\psi}\} = \psi$ , is represented by

$\begin{equation*} \mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\Big\{\overline{\psi}(x,t)\Big\} = \frac{1}{2\pi\,i}\int\limits_{c-i\infty}^{c+i\infty}e^{s t}\int\limits_{d-i\infty}^{d+i\infty}e^{px} \overline{\psi }(p,s)\,dp\,ds, \label{eq2.7} \end{equation*}$

where $Re \; (p)\geq c$ and $Re \; (s)\geq d$ , and $c, d \in \mathbb{R}$ to be chosen appropriately.

2. Existence of the solution

It is often more challenging to find the closed form of series solution to a nonlinear FDEs due to their complexity. Therefore question arises about the existence of the solution to such FDEs. For this, we utilize the applications of fixed point theory to study whether the solution of our considered system exists. So far, in the literature there exists no such theory for the existence of our considered system. We use here for the first time the $\beta-l$ -Geraghty type contraction to show that there exists a solution to the considered model. So we progress as follows

$\begin{eqnarray} ^{C}_{a}D_{t}^{\omega} \psi-\frac{\partial \psi}{\partial x}+p \psi+ q g(\psi) = r,\; \; \; \; \; \; \; \; 1 < \omega\leq2, \end{eqnarray}$

(2.1)

with

$\begin{eqnarray} \psi(x,0) = \mathcal{F}(x),\; \; \; \; \; \psi_{t}(x,0) = \mathcal{G}(x). \end{eqnarray}$

(2.2)

The above equation can also be expressed in the form

$\begin{eqnarray} ^{C}_{a}D_{t}^{\omega} \psi = \mathcal{H}(x,t,\psi),\; \; \; \; \; \; \; \; 1 < \omega\leq2, \end{eqnarray}$

(2.3)

where

$\begin{eqnarray} \mathcal{H}(x,t,\psi) = \frac{\partial \psi}{\partial x}-p \psi- q g(\psi)+r. \end{eqnarray}$

(2.4)

For the existence of the above model, we use the following notions.

Let $\Omega$ be the family of continuous and increasing functions defined as $\mathtt{l}:[0, \; \infty)\rightarrow [0, \; \infty)$ satisfying

$\begin{eqnarray*} \mathtt{l}(\mathtt{q}x)\leq\mathtt{q}\mathtt{l}(x)\leq\mathtt{q}x,\; \mathtt{q} > 1, \end{eqnarray*}$

and the elements of $\Theta$ are non-decreasing functions, such that

$\begin{eqnarray*} \varepsilon:[0,\; \infty)\rightarrow [0,\; \frac{1}{\rho_{1}^{2}}),\; \; \; \; \; where\; \rho_{1}\geq 1. \end{eqnarray*}$

Definition 3. Suppose that $(\mathcal{M}, \mathtt{d})$ be a complete $b$ -metric space: Let $\mathcal{T}:\mathcal{M}\rightarrow \mathcal{M}$ also consider that $\exists \; \mathbb{F}:\mathcal{M}\times \mathcal{M}\rightarrow [0, \; \infty)$ with $\mathbb{F}(m, n)\mathtt{l}(\rho_{1}^{3}d(\mathcal{T}_{m}, \mathcal{T}_{n}))\leq \varepsilon(\mathtt{l}d(m, n))\mathtt{l}(d(m, n))$ , for $m, n\in\mathcal{M}, \; where\; \rho_{1}\geq 1, \; \varepsilon \in \Theta\; and\; \mathtt{l}\in\Omega$ . Then $\mathcal{T}$ is known a generalized $\mathbb{F}-\mathtt{l}$ -Gergaghty type contraction mapping.

Definition 4. Consider $\mathcal{T}:\mathcal{M}\rightarrow \mathcal{M}$ , where $\mathcal{M}$ is non-empty and $\mathbb{F}:\mathcal{M}\times \mathcal{M}\rightarrow [0, \; \infty)$ , where $\beta(m, n)\geq1\Rightarrow \beta(\mathcal{M}_{m}, \mathcal{M}_{n})\; \forall\; m, n\in\mathcal{M}$ , then $\mathcal{T}$ is called $\beta$ -admissible mapping.

First we show that there exists the fixed point for the considered model Eq (2.3), for this we apply the following theorem.

Theorem 1. ^[33] Let $\mathcal{T}:\mathcal{M}\rightarrow \mathcal{M}$ be a generalized $\mathbb{F}-\mathtt{l}$ -Gergaghty type contraction such that

(1) $\mathcal{T}$ is $\beta$ -admissible.

(2) There exists $\upsilon_{0}\in \mathcal{M}\; with\; \beta(\upsilon_{0}, \mathcal{T}\upsilon_{0})\geq1$ .

(3) $\{\upsilon_{n}\}\subset\mathcal{M}, \lim_{n\rightarrow \infty}\upsilon_{n} = \upsilon$ , where $\upsilon\in \mathcal{M}$ and $\beta(\upsilon_{n}, \upsilon_{n+1})\geq1\Rightarrow \beta(\upsilon_{n}, \upsilon)\geq1$ ,

then $\exists$ a fixed point for $\mathcal{T}$ . Let $\mathcal{M} = \mathcal{C}(\pi, \mathcal{R})$ and $d:\mathcal{M}\times \mathcal{M}\rightarrow [0, \; \infty)$ , where $\pi = [0, \; 1]\times[0, \; 1]$ given by

$\begin{eqnarray*} d(\mathtt{u},\mathtt{v}) = ||(\mathtt{u}-\mathtt{v})^{2}||_{\infty} = \sup\limits_{m\in[0,\; M]\; t\in[0,\mathbb{T}]}(\mathtt{u}-\mathtt{v})^{2}, \end{eqnarray*}$

thus $(\mathcal{M}, d)$ be a complete $b$ -metric space. The following theorem shows the existence of solution of the considered model Eq (2.3).

Theorem 2. Suppose that $\exists$ $\mathcal{J}:\mathbb{R}^{2}\rightarrow \mathbb{R}$ such that

(1) $|\mathcal{H}(x, t, \psi(x, t))-\mathcal{H}(x, t, \phi(x, t))|\leq\frac{\sqrt{\alpha+1}}{3\sqrt{3}}\sqrt{\varepsilon(\mathtt{l}|\mathtt{u}-\mathtt{v}|^{2})\mathtt{l}(|\mathtt{u}-\mathtt{v}|^{2})}$ , for $x \in [0, X], \; t\in [0, T]$ , and $\mathtt{u}, \mathtt{v}\in\mathcal{M}$ with $\mathcal{J}(\mathtt{u}, \mathtt{v})\geq0$ .

(2) there exists $\mathtt{u}_{1}\in \mathcal{M}$ with $\mathcal{J}(\mathtt{u}_{1}, \mathcal{T}\mathtt{u}_{1})\geq0$ , where $\mathcal{T}:\mathcal{M}\rightarrow \mathcal{M}$ is defined by $\mathcal{T}\mathtt{u}_{j} = \mathtt{u}_{0}\mathtt{u}_{0}+\mathtt{u}_{1}\mathtt{u}_{t}+\mathcal{I}_{t}^{\omega}\mathcal{H}(x, t, \psi(x, t))$ .

(3) for $\mathtt{u}, \mathtt{v}\in \mathcal{M}, \mathcal{J}(\mathtt{u}, \mathtt{v})\geq0 \; \Rightarrow \; \mathcal{J}(\mathcal{T}\mathtt{u}_{n}, \mathcal{T}\mathtt{v})$ .

(4) $\{\mathtt{u}_{n}\}\subset\mathcal{M}$ , $\mathtt{u}_{n}\rightarrow \mathtt{u}$ where $\mathtt{u}\in \mathcal{M}$ and $mathcal{J}(\mathtt{u}_{n}, \mathtt{u}_{n+1}) \; \Rightarrow \; mathcal{J}(\mathtt{u}_{n}, \mathtt{u})\geq0$ , for $n\in N$ .Then there exists a solution of the model Eq (2.3).

Proof. Applying the fractional integral to Eq (2.3), we obtain

$\begin{eqnarray*} \psi(x,t) = C_{0}\psi(x,0)+C_{1}\psi_{t}(x,0)+\mathcal{I}_{t}^{\omega}\mathcal{H}(x,t,\psi(x,t)) = \mathcal{T}\psi(x,t). \label{exi3} \end{eqnarray*}$

Here we prove that $\mathcal{T}$ has a fixed point using the above technique, thus

$\begin{eqnarray*} \left|\mathcal{T}\psi(x,t)-\mathcal{T}\phi(x,t)\right|^{2}& = &\left|\mathcal{I}_{t}^{\omega}\mathcal{H}(x,t,\psi(x,t))-\mathcal{I}_{t}^{\omega}\mathcal{H}(x,t,\phi(x,t))\right|^{2}\\ &\leq&\mathcal{I}_{t}^{\omega}\left\{\left|\mathcal{H}(x,t,\psi(x,t))-\mathcal{H}(x,t,\phi(x,t))\right|\right\}^{2}\\ & = &\left\{\frac{1}{\sqrt{\omega}}\left|\mathcal{H}(x,t,\psi(x,t))-\mathcal{H}(x,t,\phi(x,t))\right|\right\}^{2}\\ &\leq&\left\{\frac{\omega\sqrt{\omega}}{3\sqrt{3}\sqrt{\omega}T^{\omega}}\int_{0}^{t}(t-s)^{\omega-1}\sqrt{\varepsilon(\mathtt{l}(|\mathtt{u}-\mathtt{v}|^{2})\mathtt{l}(|\mathtt{u}-\mathtt{v}|^{2}))}|\right\}^{2}\\ &\leq&\left\{\omega T^{\omega}\int_{0}^{t}(t-s)^{\omega-1}\sqrt{\varepsilon(\mathtt{l}(\sup\limits_{x\in[0,\; X]\; t\in[0,\mathbb{T}]}|\mathtt{u}-\mathtt{v}|^{2})\mathtt{l}(\sup\limits_{x\in[0,\; X]\; t\in[0,\mathbb{T}]}|\mathtt{u}-\mathtt{v}|^{2}))}|\right\}^{2}\\ &\leq&\frac{1}{3\sqrt{3}}\varepsilon (\mathtt{l}(d(\mathtt{u},\mathtt{v}))\mathtt{l}(d(\mathtt{u},\mathtt{v}))). \end{eqnarray*}$

Hence for $\mathtt{u}, \mathtt{v}\in \mathcal{C}(\pi$ with $\mathcal{J}(\mathtt{u}, \mathtt{v})\geq0$ we have $27||(\mathcal{T}\mathtt{u}-\mathcal{T}\mathtt{v})^{2}||\leq\varepsilon (\mathtt{l}(d(\mathtt{u}, \mathtt{v}))\mathtt{l}(d(\mathtt{u}, \mathtt{v})))$ . Now $\mathbb{F}:\mathcal{C}([0, \; X]\times[0, \; T], \mathbb{R})\times\mathcal{C}([0, \; X]\times[0, \; T], \mathbb{R})\rightarrow [0, \; \infty)$ by

$\begin{eqnarray*} \beta(\mathtt{u},\mathtt{v}) = \begin{cases} 1\; \; \; \; \; \; \; \; if\; \mathcal{J}(\mathtt{u},\mathtt{v})\geq0,\\ 0\; \; \; \; \; \; \; \; else, \end{cases} \end{eqnarray*}$

and

$\begin{eqnarray*} \beta(\mathtt{u},\mathtt{v})\mathtt{l}(27d(\mathcal{T}\mathtt{u},\mathcal{T}\mathtt{v}))&\leq&27d(\mathcal{T}\mathtt{u},\mathcal{T}\mathtt{v})\\ &\leq& \varepsilon(\mathtt{l}d(\mathtt{u},\mathtt{v})\mathtt{l}d(\mathtt{u},\mathtt{v})). \end{eqnarray*}$

Thus, $\mathcal{T}$ is an $\beta-\mathtt{l}$ -contraction. Now to show that $\mathcal{T}$ is $\beta$ -admissible, we have from condition $(iii)$

$\begin{eqnarray*} \beta(\mathtt{u},\mathtt{v})\geq1\Rightarrow \mathcal{J}(\mathtt{u},\mathtt{v})\geq0\Rightarrow \mathcal{J}(\mathcal{T}\mathtt{u},\mathcal{T}\mathtt{v})\geq0\Rightarrow \beta(\mathcal{T}\mathtt{u},\mathcal{T}\mathtt{v})\geq1. \end{eqnarray*}$

For $\mathtt{u}, \; \mathtt{v}\in\mathcal{C}(\pi, \mathbb{R})$ , from condition $(ii)$ we have $\mathtt{u}\in\mathcal{C}(\pi, \mathbb{R})$ . Such that $\beta(\mathtt{u}_{0}, \mathcal{T}\mathtt{u}_{0})\geq1$ . Similarly from $(iv)$ and Theorem 1, there exists $\mathtt{u}^{*}\in\mathcal{C}(\pi, \mathbb{R})$ , such that $\mathtt{u}^{*} = \mathcal{T}\mathtt{u}^{*}$ . Therefore we proved that the model Eq (2.3) has a solution.

3. Modified double Laplace decomposition method (MDLDM)

Here, we study the above technique, which is a composition of DLT with the decomposition method. This method can be applied to find the general series solutions for various PDEs/ODEs. This is an efficient technique to study the analytical solutions of several nonlinear systems ^[34]. Let us consider the general non-linear system

$\begin{eqnarray} L\,\psi+R\,\psi+N\,\psi = r(x,t). \end{eqnarray}$

(3.1)

Here, $L$ and $R$ is linear and nonlinear operators, $r(x, t)$ is some particular external function and $N$ is nonlinearity in the system. The convergence analysis of the considered technique can be seen in ^[35].

General solution of proposed model in Caputo's sense

Using the technique defined above and expressing Eq (1.1) in the form

$\begin{eqnarray} ^{C}D_t^\omega\,\psi(x,t)-\frac{\partial^{2} \psi}{\partial x^2}+p\psi+qg(\psi) = r(x,t),\, \,1 < \omega\leq 2, \end{eqnarray}$

(3.2)

with

$\begin{eqnarray} \psi(x,0) = \mathcal {F}(x),\, \, \psi_t(x,t) = \mathcal {G}(x). \end{eqnarray}$

(3.3)

Applying DLT to above equation, we obtain

$\begin{eqnarray} \mathcal{L}_x\mathcal{L}_t\{^CD_t^\omega\,\psi\}-\mathcal{L}_x\mathcal{L}_t\{\frac{\partial^2}{\partial x^2}\psi\}+p\mathcal{L}_x\mathcal{L}_t\{ \psi\} +q\mathcal{L}_x\mathcal{L}_t\{g\psi\} = \mathcal{L}_x\mathcal{L}_t\{r(x,t)\}. \end{eqnarray}$

(3.4)

Applying DLT on fractional order, gives

$\begin{eqnarray} \mathcal{L}_x\mathcal{L}_t\{\psi\}& = &\frac{1}{s}\mathcal{L}_x\{\psi(x,0)\}+\frac{1}{s^{2}}\mathcal{L}_x\{\psi_t(x,0)\}+\frac{1}{s^{\omega}}\mathcal{L}_x\mathcal{L}_t\{\frac{\partial^2}{\partial x^2}\psi\}+p\frac{1}{s^{\omega}}\mathcal{L}_x\mathcal{L}_t\{ \psi \} \\&&+q\frac{1}{s^{\omega}}\mathcal{L}_x\mathcal{L}_t\{g(\psi)\}+\mathcal{L}_x\mathcal{L}_t\{r(x,t)\}. \end{eqnarray}$

(3.5)

Similarly, applying Laplace transform on Eq (3.3), gives

$\begin{eqnarray} \mathcal{L}_x\{\psi(x,0)\} = \overline{\mathcal{F}}(p),\,\,\,\,\mathcal{L}_x\{\psi_t(x,0)\} = t\overline{\mathcal{G}}(p). \end{eqnarray}$

(3.6)

Now consider

$\begin{eqnarray} \psi = \sum\limits_{n = 0}^\infty \psi_n, \end{eqnarray}$

(3.7)

where the non-linear term can be degraded as

$\begin{eqnarray} g(\psi) = \sum\limits_{i = 0}^\infty A_n, \end{eqnarray}$

(3.8)

where $A_n$ , is given by ^[36]

$\begin{eqnarray} A_n = \frac{1}{n!}\frac{d^n}{d\lambda^n}\left[\sum\limits_{k = 0}^n \lambda^{k}g(\psi_k)\right]_{\lambda = 0}. \end{eqnarray}$

(3.9)

Finally, applying inverse DLT to Eq (3.2), using Eq (3.6) and Eq (3.9), gives

$\begin{eqnarray*} \psi_0& = &\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s}\overline{\mathscr{F}}(p,0)\right]+t{\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s^{2}}\overline{\mathscr{G}}(p,0)\right] = \psi(x,0)},\\ \psi_1& = &\mathcal{L}_x^{-1}\mathcal{L}_t^{-1} \left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{\psi_{0xx}\}\right] -p\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{\psi_{0}\}\right]-q\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{A_0\}\right]\\ &+&\left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{r(x,t)\}\right],\\ \psi_2& = &\mathcal{L}_x^{-1}\mathcal{L}_t^{-1} \left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{\psi_{1xx}\}\right] -p\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{\psi_{1}\}\right]-q\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{A_1\}\right],\\ \psi_3& = &\mathcal{L}_x^{-1}\mathcal{L}_t^{-1} \left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{\psi_{2xx}\}\right] -p\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{\psi_{2}\}\right]-q\mathcal{L}_x^{-1}\mathcal{L}_t^{-1}\left[\frac{1}{s^\omega}\mathcal{L}_x\mathcal{L}_t\{A_2\}\right].\\ \end{eqnarray*}$

In a similar manner, other terms can be computed. Final result can be obtained as

$\begin{eqnarray} \psi(x,t) = \sum\limits_{n = 0}^\infty \psi_n(x,t). \end{eqnarray}$

(3.10)

which is the general solution of Eq (3.2) in series form by using the proposed method as discussed above.

4. Applications of MDLDM

Here, we present numerical examples on the TFKG equation in Caputo's sense given as Eq (3.2) and discuss the behaviour of each example. We apply the aforesaid technique discussed in Section 3, to obtain the approximate solution of the problems.

Example 1. Consider the nonlinear TFKG equation

$\begin{eqnarray} D_{t}^{\omega}\psi-\frac{\partial^{2}\psi}{\partial x^{2}}+\frac{3}{4}\psi-\frac{3}{2}\psi^{3} = 0,\,\,\,1 < \omega\leq 2, \,\,\,p = \frac{3}{4},\,\,\, q = \frac{3}{2}, \end{eqnarray}$

(4.1)

$\begin{eqnarray} g(\psi) = \psi^{3}, \,r(x,t) = 0, \end{eqnarray}$

(4.2)

with

$\begin{eqnarray} \psi(x,0) = - \operatorname{sech}(x),\,\,\,\,\,\,\,\,\psi_{t}(x,0) = \frac{1}{2} \operatorname{sech}(x)\tanh(x). \end{eqnarray}$

(4.3)

For $\alpha = 2$ , the exact solution of Eq (4.1) can be obtained in the form ^[19]

$\begin{eqnarray} \psi(x,t) = -sech\left( x+\frac{t}{2}\right). \end{eqnarray}$

(4.4)

Consider TFKG Eq (4.1) in Caputo's sense

$\begin{eqnarray} ^{C}D_{t}^{\omega}\psi-\frac{\partial^{2}\psi}{\partial x^{2}}+\frac{3}{4}\psi-\frac{3}{2}\psi^{3} = 0,\,\,\,\,\,1 < \omega\leq 2. \end{eqnarray}$

(4.5)

Applying MDLDM scheme discussed in Section 3, we obtain

$\begin{eqnarray*} \psi_0& = &- \operatorname{sech}(x)+\frac{t}{2} \operatorname{sech}(x)\tanh(x),\\ \psi_1& = &\frac{t^{\omega}}{\Gamma(\omega+1)}\left[-\frac{1}{4}-\frac{3}{2} \operatorname{sech}^{2}(x)\right] \operatorname{sech}(x) +\frac{t^{\omega+1}}{\Gamma(\omega+2)}\left[\frac{11}{8}-\frac{7}{4} \operatorname{sech}^{2}(x)\right] \operatorname{sech}(x)\tanh(x)\\&& -\frac{2!t^{\omega+2}}{\Gamma(\omega+3)}\left[\frac{9}{8}sech^{3}(x)tanh^{2}(x)\right] +\frac{3!t^{\omega+3}}{\Gamma(\omega+4)}\left[\frac{3}{16}sech^{3}(x)tanh^{3}(x)\right]. \end{eqnarray*}$

In a similar manner, other terms can be computed. Final result can be obtained as

$\begin{eqnarray} \psi(x,t) = \sum\limits_{n = 0}^\infty \psi_n(x,t). \end{eqnarray}$

(4.6)

Discussion

The error analysis between series solution Eq (4.6) and the exact solution Eq (4.4) is shown in Table 1, while the surface behaviour is shown in reveals that Eq (4.1) is depends on time (t). It should be noted that when the time (t) is small enough, there is less extent of error exist amongst the approximate and exact solutions obtained by the MDLDM method. [left panel] shows the absolute of wave solution Eq (4.6) with deviations in $(\alpha)$ with $t = 0.65$ in comparison with exact solution Eq (4.4). Notice that numerical result, Eq (4.6) exactly matches to the exact solution Eq (4.4). This shows that the governing equation admits a soliton solution. [right panel] represents Eq (4.6) reveals that the amplitude of the solitary potentials blow-up as t increases. The $3D$ profiles for Eq (4.6) is shown in versus x for $t = 0.65$ . It reveals the progression of localized mode in the governing system. The solution obtained as Eq (4.6) versus x for $t = 1\, (dashed\, line), \, 0.8\, (solid\, curve), \, 0.6\, (dotted\, curve)$ , when $\omega = 2\, and \, 1.7$ are also depicted. Clearly, one can see the wave amplitude enhancement with variations in t that concludes that coefficient $(\omega)$ considerably increases the wave amplitudes.

Table 1. Comparison between the exact solution (4.4) with approximate solution obtained in the form (4.6).

(x, t)	$\mid$ Exact $\mid$	$\mid \; \psi \; \mid$	$\mid$ Exact $-$ $\psi$ $\mid$	(x, t)	$\mid$ Exact $\mid$	$\mid \; \psi \; \mid$	$\mid$ Exact $-$ $\psi$ $\mid$
(-6, 0.6)	0.0067	0.0067	2.4083 $\times 10^{-5}$	(-4, 0.6)	0.0494	0.0492	1.7609 $\times 10^{-4}$
(-2, 0.6)	0.3536	0.3529	6.2194 $\times 10^{-4}$	(0, 0.6)	0.9566	0.9550	1.6000 $\times 10^{-3}$
(2, 0.6)	0.1985	0.1992	6.7883 $\times 10^{-4}$	(4, 0.6)	0.0271	0.0273	1.5198 $\times 10^{-4}$
(6, 0.6)	0.0037	0.0037	2.0728 $\times 10^{-5}$	(-6, 0.2)	0.0055	0.0055	8.4718 $\times 10^{-7}$
(-4, 0.2)	0.0405	0.0405	6.2013 $\times 10^{-6}$	(-2, 0.2)	0.2926	0.2926	2.4219 $\times 10^{-5}$
(0, 0.2)	0.9950	0.9950	2.0749 $\times 10^{-5}$	(2, 0.2)	0.2413	0.2413	2.4874 $\times 10^{-5}$
(4, 0.2)	0.0331	0.0331	5.9043 $\times 10^{-6}$	(6, 0.2)	0.0045	0.0045	8.0587 $\times 10^{-7}$
(-6, 0.05)	0.0051	0.0051	1.2989 $\times 10^{-8}$	(-4, 0.05)	0.0375	0.0375	9.5114 $\times 10^{-8}$
(-2, 0.05)	0.2723	0.2723	3.8310 $\times 10^{-7}$	(0, 0.05)	0.9997	0.9997	8.1360 $\times 10^{-8}$
(2, 0.05)	0.2595	0.2595	3.8564 $\times 10^{-7}$	(4, 0.05)	0.0357	0.0357	9.3954 $\times 10^{-8}$
(6, 0.05)	0.0048	0.0048	1.2828 $\times 10^{-8}$

| Show Table

DownLoad: CSV

Figure 1. The surface plot of the error analysis given in Table 1.

DownLoad: Full-Size Img PowerPoint

Figure 2. The left plot portrays comparison between Eq (4.4) and Eq (4.6) for various values of

$\omega$ , while the right panel portrays solution profiles of

$\psi(x, t)$ vs t for various values of

$\omega$ .

DownLoad: Full-Size Img PowerPoint

Figure 3. The surface plot of Eq (4.6) for the parameters used in the Figure 2's left panel.

DownLoad: Full-Size Img PowerPoint

Figure 4. The solution profiles of Eq (4.6) for desperate values of

$\omega$ with desperate values of t.

DownLoad: Full-Size Img PowerPoint

Example 2. Consider the time-fractional nonlinear KGE in the form

$\begin{eqnarray} D_t^\omega\,\phi-\frac{\partial^{2} \phi}{\partial x^2}+q\phi^{3} = 0,\,\,\,1 < \omega\leq 2,\,\,\, p = 0\,\,\, q = 1, \end{eqnarray}$

(4.7)

$\begin{eqnarray} y(\phi) = \phi^{3},\,\,\,r(x,t) = 0, \end{eqnarray}$

(4.8)

with

$\begin{eqnarray} \phi(x,0) = Rtan(\lambda\,x),\,\,\,\,\,\phi_{t}(x,0) = R\eta\lambda sec^{2}(\lambda\,x), \end{eqnarray}$

(4.9)

where

$\begin{eqnarray*} R = \sqrt{\frac{\rho}{\kappa}}\,\,\,and\,\,\,\,\,\lambda = \sqrt{\frac{-\rho}{2(\sigma+\eta^{2})}}. \end{eqnarray*}$

The parameters $\rho$ , $\kappa$ , $\sigma$ , and $\eta \; \in \mathbb{R}$ . It should be noted that, for $\alpha = 2$ , an exact solution of Eq (4.7) can be obtained in the form ^[37]

$\begin{eqnarray} \phi(x,t) = R\,\tan[\lambda(x+\eta\,t)]. \end{eqnarray}$

(4.10)

Writing Eq (4.7) in Caputo's sense gives

$\begin{eqnarray} ^{C}D_t^{\omega}\,\phi-\frac{\partial^{2} \phi}{\partial x^2}+\phi^{3} = 0,\,\,\,1 < \omega\leq2. \end{eqnarray}$

(4.11)

The series solution of Eq (4.11) with conditions (4.9) gives

$\begin{eqnarray*} \phi_0& = &Rtan(\lambda\,x)+tR\eta\lambda sec^{2}(\lambda\,x),\\ \phi_1& = &\left(\frac{t^{\omega}}{\Gamma(\omega+1)}\right) \left[2R\lambda^{2}sec^{2}(\lambda\,x)tan(\lambda\,x)-R^{2}tan^{3}(\lambda\,x)\right]\\&& +R\eta\lambda\left(\frac{t^{\omega+1}}{\Gamma(\omega+2)}\right) \left[4\lambda^{2}sec^{2}(\lambda\,x)tan^{2}(\lambda\,x)+2\lambda^{2}sec^{4}(\lambda\,x)-3R^{2}tan^{2}(\lambda\,x)sec^{2}(\lambda\,x)\right]\\&& -3R^{3}\lambda^{2}\eta^{2}\left(\frac{2!t^{\omega+2}}{\Gamma(\omega+3)}\right) \left[sec^{4}(\lambda\,x)tan(\lambda\,x)\right]-R^{3}\lambda^{3}\eta^{3}\left(\frac{3!t^{\omega+3}}{\Gamma(\omega+4)}\right) \left[sech^{6}(\lambda\,x)\right]. \end{eqnarray*}$

In a similar manner, other terms can be computed. The final result can be written in the form

$\begin{eqnarray} \phi & = & \sum\limits_{n = 0}^\infty \phi_n. \end{eqnarray}$

(4.12)

Discussion

The parameters as $\kappa = 1, \sigma = -8.5, \, \eta = 0.05$ , and $\rho = -1$ are considered for numerical illustration. The error amongst the approximate and exact solutions of Eq (4.7) is shown in Table 2 and its corresponding surface plot is presented in . The numerical solution, Eq (4.12) and exact solution Eq (4.10) is depicted in [left panel], for $t = 7$ and for different values of time-fractional coefficient $(\omega)$ . It is noted that TFKG Eq (4.12) may admits the excitation in the system. This amount enrichment in $\omega$ overturned the wave amplitude as it interrupt the dispersion/nonlinearity effects. To see the effect of a temporal variable (t) on the wave solution, Eq (4.12) is displayed in [right panel] which shows that $\phi(x, t)$ increases with time. Further, the $3D$ profiles for Eq (4.12) is shown versus x with $t = 7$ in for $(\omega = 2)$ , which represents the physical behaviour of Eq (4.12). It shows the advancement of localized mode excitations in the governing equation. The solution of Eq (4.12) versus x with $t = 7$ (dashed curve), 6 (solid green curve), 5 (dotted curve), in , with $\omega = 2\, and\, 1.5$ respectively is depicted. Clearly, the wave amplitude increases with deviations in t. It infers that the fractional order $(\omega)$ significantly increases the wave amplitudes.

Table 2. Comparison between approximate solution obtained in the form (4.12) with exact solution (4.7).

(x, t)	$\mid$ Exact $\mid$	$\mid \; \phi \; \mid$	$\mid$ Exact $-$ $\phi$ $\mid$	(x, t)	$\mid$ Exact $\mid$	$\mid \; \phi \; \mid$	$\mid$ Exact $-$ $\phi$ $\mid$
(-10, 0.02)	0.9845	0.9845	4 $\times 10^{-7}$	(-8, 0.02)	0.9596	0.9596	4 $\times 10^{-6}$
(-6, 0.02)	0.8968	0.8965	3 $\times 10^{-5}$	(-4, 0.02)	0.7488	0.7486	2 $\times 10^{-5}$
(-2, 0.02)	0.4504	0.4503	1 $\times 10^{-5}$	(0, 0.02)	0	0	0
(2, 0.02)	0.4504	0.4503	1 $\times 10^{-5}$	(4, 0.02)	0.7488	0.7486	2 $\times 10^{-5}$
(6, 0.02)	0.8968	0.8965	3 $\times 10^{-5}$	(8, 0.02)	0.9596	0.9596	4 $\times 10^{-6}$
(10, 0.02)	0.9845	0.9845	4 $\times 10^{-7}$	(-10, 0.05)	0.9845	0.9845	1.8751 $\times 10^{-5}$
(-8, 0.05)	0.9595	0.9596	4.8543 $\times 10^{-5}$	(-6, 0.05)	0.8966	0.8968	1.2196 $\times 10^{-4}$
(-4, 0.05)	0.7486	0.7488	2.8545 $\times 10^{-4}$	(-2, 0.05)	0.4499	0.4499	5.8734 $\times 10^{-4}$
(0, 0.05)	0.0012	0.0012	1.1300 $\times 10^{-6}$	(2, 0.05)	0.4508	0.4505	3.7952 $\times 10^{-4}$
(4, 0.05)	0.7491	0.7488	2.4730 $\times 10^{-4}$	(6, 0.05)	0.8969	0.8968	1.1554 $\times 10^{-4}$
(8, 0.05)	0.9596	0.9596	4.7525 $\times 10^{-5}$	(10, 0.05)	0.9845	0.9845	1.8587 $\times 10^{-5}$

| Show Table

DownLoad: CSV

Figure 5. The surface plot for Table 2.

DownLoad: Full-Size Img PowerPoint

Figure 6. Comparison between Eq (4.10) and Eq (4.12) for different values of

$\omega$ [left panel]. The solution profiles of

$\psi(x, t)$ against time (t) interval with various values of

$\omega$ [right panel].

DownLoad: Full-Size Img PowerPoint

Figure 7. The surface plot for the parameters used for the left panel of Figure 6.

DownLoad: Full-Size Img PowerPoint

Figure 8. The solution profiles of Eq (4.12) for different

$\omega$ with different values of time(t).

DownLoad: Full-Size Img PowerPoint

5. Conclusions

We have studied the TFKG equation using double Laplace transforms with the decomposition method. The general solution of proposed system is obtained as a class of general series solution. It is relevant to note that following only two iterations, fairly precise results are obtained that converges to the exact solution of the governing equation. The proposed method offers perfect numerical results without any alteration and complicated numerical methods for the governing equation in the fractional case. The numerical results obtained for particular examples are compared with the exact solutions at the classical order. The result profiles with physical interpretations for different fraction orders were revealed explicitly.

Acknowledgments

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R8). Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflict of interest

It is declared that all the authors have no conflict of interest regarding this manuscript.

References

[1]	H. El-Metwally, E. A. Grove, G. Ladas, R. Levins, M. Radin, On the difference equation $x_{n+1} = \alpha +\beta x_{n-1}e^{-x_{n}}$ , Nonlinear Anal., 47 (2001), 4623–4634. https://doi.org/10.1016/S0362-546X(01)00575-2 doi: 10.1016/S0362-546X(01)00575-2
[2]	E. A. Grove, G. Ladas, N. R. Prokup, R. Levins, On the global behavior of solutions of a biological model, Commun. Appl. Nonlinear Anal, 7(2000), 33–46.
[3]	G. Papaschinopoulos, G. Ellina, K. B. Papadopoulos, Asymptotic behavior of the positive solutions of an exponential type system of difference equations, Appl. Math. Comput., 245 (2014), 181–190. https://doi.org/10.1016/j.amc.2014.07.074 doi: 10.1016/j.amc.2014.07.074
[4]	W. J. Wang, H. Feng, On the dynamics of positive solutions for the difference equation in a new population model, J. Nonlinear Sci. Appl., 9 (2016), 1748–1754. http://dx.doi.org/10.22436/jnsa.009.04.30 doi: 10.22436/jnsa.009.04.30
[5]	G. Stefanidou, G. Papaschinopoulos, C. J. Schinas, On a system of two exponential type difference equations, Commun. Appl. Nonlinear Anal., 17 (2010), 1–13.
[6]	G. Papaschinopoulos, M. A. Radin, C. J. Schinas, On the system of two difference equations of exponential form: $x_{n+1} = a+bx_{n-1}e^{-y_{n}}$ , $y_{n+1} = c+dy_{n-1}e^{-x_{n}}$ , Math. Comput. Model., 54 (2011), 2969–2977. https://doi.org/10.1016/j.mcm.2011.07.019 doi: 10.1016/j.mcm.2011.07.019
[7]	G. Papaschinopoulos, M. A. Radin, C. J. Schinas, Study of the asymptotic behavior of the solutions of three systems of difference equations of exponential form, Appl. Math. Comput., 218 (2012), 5310–5318. https://doi.org/10.1016/j.amc.2011.11.014 doi: 10.1016/j.amc.2011.11.014
[8]	I. Ozturk, F. Bozkurt, S. Ozen, On the difference equation $y_{n+1} = (\alpha +\beta e^{-y_{n}})/(\gamma +y_{n-1}$ ), Appl. Math. Comput., 181 (2006), 1387–1393. https://doi.org/10.1016/j.amc.2006.03.007 doi: 10.1016/j.amc.2006.03.007
[9]	G. Papaschinopoulos, C. J. Schinas, On the dynamics of two exponential type systems of difference equations, Comput. Math. Appl., 64 (2012), 2326–2334. https://doi.org/10.1016/j.camwa.2012.04.002 doi: 10.1016/j.camwa.2012.04.002
[10]	G. Papaschinopoulos, N. Fotiades, C. J. Schinas, On a system of difference equations including negative exponential terms, J. Differ. Equ. Appl., 20 (2014), 717–732. https://doi.org/10.1080/10236198.2013.814647 doi: 10.1080/10236198.2013.814647
[11]	M. Pituk, More on Poincaré's and Peron's theorems for difference equations, J. Differ. Equ. Appl., 8 (2002), 201–216. https://doi.org/10.1080/10236190211954 doi: 10.1080/10236190211954

This article has been cited by:

Yuan Xu, Siu-Long Lei, Hai-Wei Sun, An efficient red–black skewed extrapolation cascadic multigrid method for two-dimensional Poisson equation, 2024, 43, 2238-3603, 10.1007/s40314-023-02514-4

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Mathematics

1.8 3.4

Metrics

Article views(1322) PDF downloads(68) Cited by(4)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(2)

AIMS Mathematics

On stability analysis of a class of three-dimensional system of exponential difference equations

Related Papers:

Abstract

1. Introduction

2. Existence of the solution

3. Modified double Laplace decomposition method (MDLDM)

4. Applications of MDLDM

5. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

AIMS Mathematics

On stability analysis of a class of three-dimensional system of exponential difference equations

Related Papers:

Abstract

1. Introduction

2. Existence of the solution

3. Modified double Laplace decomposition method (MDLDM)

4. Applications of MDLDM

5. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog