Numerical approximation of a variable-order time fractional advection-reaction-diffusion model via shifted Gegenbauer polynomials

Yumei Chen; Jiajie Zhang; Chao Pan; Yumei Chen; Jiajie Zhang; Chao Pan

doi:10.3934/math.2022855

AIMS Mathematics

2022, Volume 7, Issue 8: 15612-15632. doi: 10.3934/math.2022855

Previous Article Next Article

Research article

Numerical approximation of a variable-order time fractional advection-reaction-diffusion model via shifted Gegenbauer polynomials

1.
College of Mathematics Education, China West Normal University, Nanchong 637009, China
2.
School of Mathematics and Information, China West Normal University, Nanchong 637009, China

Received: 11 April 2022 Revised: 11 June 2022 Accepted: 16 June 2022 Published: 23 June 2022
MSC : 65N35, 34A08

The fractional advection-reaction-diffusion equation plays a key role in describing the processes of multiple species transported by a fluid. Different numerical methods have been proposed for the case of fixed-order derivatives, while there are no such methods for the generalization of variable-order cases. In this paper, a numerical treatment is given to solve a variable-order model with time fractional derivative defined in the Atangana-Baleanu-Caputo sense. By using shifted Gegenbauer cardinal function, this approach is based on the application of spectral collocation method and operator matrices. Then the desired problem is transformed into solving a nonlinear system, which can greatly simplifies the solution process. Numerical experiments are presented to illustrate the effectiveness and accuracy of the proposed method.

Keywords:

Citation: Yumei Chen, Jiajie Zhang, Chao Pan. Numerical approximation of a variable-order time fractional advection-reaction-diffusion model via shifted Gegenbauer polynomials[J]. AIMS Mathematics, 2022, 7(8): 15612-15632. doi: 10.3934/math.2022855

Related Papers:

[1]	Ahmed Ghezal, Mohamed Balegh, Imane Zemmouri . Solutions and local stability of the Jacobsthal system of difference equations. AIMS Mathematics, 2024, 9(2): 3576-3591. doi: 10.3934/math.2024175
[2]	Hashem Althagafi, Ahmed Ghezal . Solving a system of nonlinear difference equations with bilinear dynamics. AIMS Mathematics, 2024, 9(12): 34067-34089. doi: 10.3934/math.20241624
[3]	M. T. Alharthi . Correction: On the solutions of some systems of rational difference equations. AIMS Mathematics, 2025, 10(2): 2277-2278. doi: 10.3934/math.2025105
[4]	M. T. Alharthi . On the solutions of some systems of rational difference equations. AIMS Mathematics, 2024, 9(11): 30320-30347. doi: 10.3934/math.20241463
[5]	Eunjung Lee, Dojin Kim . Stability analysis of the implicit finite difference schemes for nonlinear Schrödinger equation. AIMS Mathematics, 2022, 7(9): 16349-16365. doi: 10.3934/math.2022893
[6]	Shulan Kong, Chengbin Wang, Yawen Sun . A recursive filter for a class of two-dimensional nonlinear stochastic systems. AIMS Mathematics, 2025, 10(1): 1741-1756. doi: 10.3934/math.2025079
[7]	Yeyang Jiang, Zhihua Liao, Di Qiu . The existence of entire solutions of some systems of the Fermat type differential-difference equations. AIMS Mathematics, 2022, 7(10): 17685-17698. doi: 10.3934/math.2022974
[8]	Abdulghani R. Alharbi . Traveling-wave and numerical solutions to nonlinear evolution equations via modern computational techniques. AIMS Mathematics, 2024, 9(1): 1323-1345. doi: 10.3934/math.2024065
[9]	Ibraheem M. Alsulami, E. M. Elsayed . On a class of nonlinear rational systems of difference equations. AIMS Mathematics, 2023, 8(7): 15466-15485. doi: 10.3934/math.2023789
[10]	Nan Li, Jiachuan Geng, Lianzhong Yang . Some results on transcendental entire solutions to certain nonlinear differential-difference equations. AIMS Mathematics, 2021, 6(8): 8107-8126. doi: 10.3934/math.2021470

Abstract

1. Introduction

This paper is devoted to study the expressions forms of the solutions and periodic nature of the following third-order rational systems of difference equations

$\begin{equation*} x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}\pm x_{n-2}}, \;y_{n+1} = \dfrac{z_{n-1}x_{n} }{x_{n}\pm y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}\pm z_{n-2}}, \end{equation*}$

with initial conditions are non-zero real numbers.

In the recent years, there has been great concern in studying the systems of difference equations. One of the most important reasons for this is a exigency for some mechanization which can be used in discussing equations emerge in mathematical models characterizing real life situations in economic, genetics, probability theory, psychology, population biology and so on.

Difference equations display naturally as discrete peer and as numerical solutions of differential equations having more applications in ecology, biology, physics, economy, and so forth. For all that the difference equations are quite simple in expressions, it is frequently difficult to realize completely the dynamics of their solutions see ^{[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]} and the related references therein.

There are some papers dealed with the difference equations systems, for example, The periodic nature of the solutions of the nonlinear difference equations system

$\begin{equation*} A_{n+1} = \dfrac{1}{C_{n}}, \;B_{n+1} = \dfrac{B_{n}}{A_{n-1}B_{n-1}}, \;C_{n+1} = \dfrac{1}{A_{n-1}}, \end{equation*}$

has been studied by Cinar in ^[7].

Almatrafi ^[3] determined the analytical solutions of the following systems of rational recursive equations

$\begin{equation*} x_{n+1} = \dfrac{x_{n-1}y_{n-3}}{y_{n-1}(\pm 1\pm x_{n-1}y_{n-3})}, \;y_{n+1} = \dfrac{y_{n-1}x_{n-3}}{x_{n-1}(\pm 1\pm y_{n-1}x_{n-3})}. \end{equation*}$

In ^[20], Khaliq and Shoaib studied the local and global asymptotic behavior of non-negative equilibrium points of a three-dimensional system of two order rational difference equations

$\begin{equation*} x_{n+1} = \dfrac{x_{n-1}}{\varepsilon +x_{n-1}y_{n-1}z_{n-1}}, \;y_{n+1} = \dfrac{ y_{n-1}}{\zeta +x_{n-1}y_{n-1}z_{n-1}}, \ z_{n+1} = \dfrac{z_{n-1}}{\eta +x_{n-1}y_{n-1}z_{n-1}}. \end{equation*}$

In ^[9], Elabbasy et al. obtained the form of the solutions of some cases of the following system of difference equations

$\begin{eqnarray*} x_{n+1} & = &\dfrac{a_{1}+a_{2}y_{n}}{a_{3}z_{n}+a_{4}x_{n-1}z_{n}}, \ y_{n+1} = \dfrac{b_{1}z_{n-1}+b_{2}z_{n}}{b_{3}x_{n}y_{n}+b_{4}x_{n}y_{n-1}}, \\ z_{n+1} & = &\dfrac{c_{1}z_{n-1}+c_{2}z_{n}}{ c_{3}x_{n-1}y_{n-1}+c_{4}x_{n-1}y_{n}+c_{5}x_{n}y_{n}}. \end{eqnarray*}$

In ^[12], Elsayed et al. have got the solutions of the systems of rational higher order difference equations

$\begin{equation*} A_{n+1} = \dfrac{1}{A_{n-p}B_{n-p}}, \;\;B_{n+1} = \dfrac{A_{n-p}B_{n-p}}{ A_{n-q}B_{n-q}}, \end{equation*}$

and

$\begin{equation*} A_{n+1} = \dfrac{1}{A_{n-p}B_{n-p}C_{n-p}}, \;\;B_{n+1} = \dfrac{ A_{n-p}B_{n-p}C_{n-p}}{A_{n-q}B_{n-q}C_{n-q}}, \;C_{n+1} = \dfrac{ A_{n-q}B_{n-q}C_{n-q}}{A_{n-r}B_{n-r}C_{n-r}}. \end{equation*}$

Kurbanli ^[25,26] investigated the behavior of the solutions of the following systems

$\begin{eqnarray*} A_{n+1} & = &\dfrac{A_{n-1}}{A_{n-1}B_{n}-1}, \;B_{n+1} = \dfrac{B_{n-1}}{ B_{n-1}A_{n}-1}, \ \ C_{n+1} = \dfrac{1}{C_{n}B_{n}}, \\ A_{n+1} & = &\dfrac{A_{n-1}}{A_{n-1}B_{n}-1}, \;B_{n+1} = \dfrac{B_{n-1}}{ B_{n-1}A_{n}-1}, \ \ C_{n+1} = \dfrac{C_{n-1}}{C_{n-1}B_{n}-1}. \end{eqnarray*}$

In ^[32], Yalçınkaya has obtained the conditions for the global asymptotically stable of the system

$\begin{equation*} A_{n+1} = \dfrac{B_{n}A_{n-1}+a}{B_{n}+A_{n-1}}, \;B_{n+1} = \dfrac{A_{n}B_{n-1}+a }{A_{n}+B_{n-1}}. \end{equation*}$

Zhang et al. ^[39] investigated the persistence, boundedness and the global asymptotically stable of the solutions of the following system

$\begin{equation*} R_{n} = A+\dfrac{1}{Q_{n-p}}, \ \ \ \;Q_{n} = A+\dfrac{Q_{n-1}}{R_{n-r}Q_{n-s}}. \end{equation*}$

Similar to difference equations and systems were studied see ^{[21,22,23,24,27,28,29,30,31,32,33,34,35,36,37,38]}.

2. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

In this section, we obtain the expressions form of the solutions of the following three dimension system of difference equations

$\begin{equation} x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \;y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}+z_{n-2}}, \end{equation}$

(1)

where $n\in \mathbb{N}_{0}$ and the initial conditions are non-zero real numbers.

Theorem 1. We assume that $\{x_{n}, y_{n}, z_{n}\}$ are solutions of system (1).Then

$\begin{eqnarray*} x_{6n-2} & = &\frac{ak^{3n}}{\prod \limits_{i = 0}^{n-1}(a+(6i)k)(a+(6i+2)k)(a+(6i+4)k)}, \\ x_{6n-1} & = &\frac{bf^{3n}}{\prod \limits_{i = 0}^{n-1}(g+(6i+1)f)(g+(6i+3)f)(g+(6i+5)f)}, \\ x_{6n} & = &\frac{c^{3n+1}}{\prod \limits_{i = 0}^{n-1}(d+(6i+2)c)(d+(6i+4)c)(d+(6i+6)c)}, \\ x_{6n+1} & = &\frac{ek^{3n+1}}{(a+k)\prod \limits_{i = 0}^{n-1}(a+(6i+3)k)(a+(6i+5)k)(a+(6i+7)k)}, \end{eqnarray*}$

$\begin{eqnarray*} x_{6n+2} & = &\frac{f^{3n+2}}{(g+2f)\prod \limits_{i = 0}^{n-1}(g+(6i+4)f)(g+(6i+6)f)(g+(6i+8)f)}, \\ x_{6n+3} & = &\frac{hc^{3n+2}}{(d+c)(d+3c)\prod \limits_{i = 0}^{n-1}(d+(6i+5)c)(d+(6i+7)c)(d+(6i+9)c)}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-2} & = &\frac{dc^{3n}}{\prod \limits_{i = 0}^{n-1}(d+(6i)c)(d+(6i+2)c)(d+(6i+4)c)}, \\ y_{6n-1} & = &\frac{ek^{3n}}{\prod \limits_{i = 0}^{n-1}(a+(6i+1)k)(a+(6i+3)k)(a+(6i+5)k)}, \\ y_{6n} & = &\frac{f^{3n+1}}{\prod \limits_{i = 0}^{n-1}(g+(6i+2)f)(g+(6i+4)f)(g+(6i+6)f)}, \\ y_{6n+1} & = &\frac{hc^{3n+1}}{(d+c)\prod \limits_{i = 0}^{n-1}(d+(6i+3)c)(d+(6i+5)c)(d+(6i+7)c)}, \\ y_{6n+2} & = &\frac{k^{3n+2}}{(a+2k)\prod \limits_{i = 0}^{n-1}(a+(6i+4)k)(a+(6i+6)k)(a+(6i+8)k)}, \\ y_{6n+3} & = &\frac{bf^{3n+2}}{(g+f)(g+3f)\prod \limits_{i = 0}^{n-1}(g+(6i+5)f)(g+(6i+7)f)(g+(6i+9)f)}, \end{eqnarray*}$

and

$\begin{eqnarray*} z_{6n-2} & = &\frac{gf^{3n}}{\prod \limits_{i = 0}^{n-1}(g+(6i)f)(g+(6i+2)f)(g+(6i+4)f)}, \\ z_{6n-1} & = &\frac{hc^{3n}}{\prod \limits_{i = 0}^{n-1}(d+(6i+1)c)(d+(6i+3)c)(d+(6i+5)c)}, \\ z_{6n} & = &\frac{k^{3n+1}}{\prod \limits_{i = 0}^{n-1}(a+(6i+2)k)(a+(6i+4)k)(a+(6i+6)k)}, \\ z_{6n+1} & = &\frac{bf^{3n+1}}{(g+f)\prod \limits_{i = 0}^{n-1}(g+(6i+3)f)(g+(6i+5)f)(g+(6i+7)f)}, \end{eqnarray*}$

$\begin{eqnarray*} z_{6n+2} & = &\frac{c^{3n+2}}{(d+2c)\prod \limits_{i = 0}^{n-1}(d+(6i+4)c)(d+(6i+6)c)(d+(6i+8)c)}, \\ z_{6n+3} & = &\frac{ek^{3n+2}}{(a+k)(a+3k)\prod \limits_{i = 0}^{n-1}(a+(6i+5)k)(a+(6i+7)k)(a+(6i+9)k)}, \end{eqnarray*}$

where $x_{-2} = a, \ x_{-1} = b,$ $x_{0} = c,$ $y_{-2} = d,$ $y_{-1} = e,$ $y_{0} = f,$ $z_{-2} = g,$ $z_{-1} = h$ $and$ $z_{0} = k.$

Proof. For $n = 0$ the result holds. Now assume that $n > 1$ and that our assumption holds for $n-1$ , that is,

$\begin{eqnarray*} x_{6n-8} & = &\frac{ak^{3n-3}}{\prod \limits_{i = 0}^{n-2}(a+(6i)k)(a+(6i+2)k)(a+(6i+4)k)}, \\ x_{6n-7} & = &\frac{bf^{3n-3}}{\prod \limits_{i = 0}^{n-2}(g+(6i+1)f)(g+(6i+3)f)(g+(6i+5)f)}, \\ x_{6n-6} & = &\frac{c^{3n-2}}{\prod \limits_{i = 0}^{n-2}(d+(6i+2)c)(d+(6i+4)c)(d+(6i+6)c)}, \\ x_{6n-5} & = &\frac{ek^{3n-2}}{(a+k)\prod \limits_{i = 0}^{n-2}(a+(6i+3)k)(a+(6i+5)k)(a+(6i+7)k)}, \\ x_{6n-4} & = &\frac{f^{3n-1}}{(g+2f)\prod \limits_{i = 0}^{n-2}(g+(6i+4)f)(g+(6i+6)f)(g+(6i+8)f)}, \\ x_{6n-3} & = &\frac{hc^{3n-1}}{(d+c)(d+3c)\prod \limits_{i = 0}^{n-2}(d+(6i+5)c)(d+(6i+7)c)(d+(6i+9)c)}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-8} & = &\frac{dc^{3n-3}}{\prod \limits_{i = 0}^{n-2}(d+(6i)c)(d+(6i+2)c)(d+(6i+4)c)}, \\ y_{6n-7} & = &\frac{ek^{3n-3}}{\prod \limits_{i = 0}^{n-2}(a+(6i+1)k)(a+(6i+3)k)(a+(6i+5)k)}, \\ y_{6n-6} & = &\frac{f^{3n-2}}{\prod \limits_{i = 0}^{n-2}(g+(6i+2)f)(g+(6i+4)f)(g+(6i+6)f)}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-5} & = &\frac{hc^{3n-2}}{(d+c)\prod \limits_{i = 0}^{n-2}(d+(6i+3)c)(d+(6i+5)c)(d+(6i+7)c)}, \\ y_{6n-4} & = &\frac{k^{3n-1}}{(a+2k)\prod \limits_{i = 0}^{n-2}(a+(6i+4)k)(a+(6i+6)k)(a+(6i+8)k)}, \\ y_{6n-3} & = &\frac{bf^{3n-1}}{(g+f)(g+3f)\prod \limits_{i = 0}^{n-2}(g+(6i+5)f)(g+(6i+7)f)(g+(6i+9)f)}, \end{eqnarray*}$

and

$\begin{eqnarray*} z_{6n-8} & = &\frac{gf^{3n-3}}{\prod \limits_{i = 0}^{n-2}(g+(6i)f)(g+(6i+2)f)(g+(6i+4)f)}, \\ z_{6n-7} & = &\frac{hc^{3n-3}}{\prod \limits_{i = 0}^{n-2}(d+(6i+1)c)(d+(6i+3)c)(d+(6i+5)c)}, \\ z_{6n-6} & = &\frac{k^{3n-2}}{\prod \limits_{i = 0}^{n-2}(a+(6i+2)k)(a+(6i+4)k)(a+(6i+6)k)}, \\ z_{6n-5} & = &\frac{bf^{3n-2}}{(g+f)\prod \limits_{i = 0}^{n-2}(g+(6i+3)f)(g+(6i+5)f)(g+(6i+7)f)}, \\ z_{6n-4} & = &\frac{c^{3n-1}}{(d+2c)\prod \limits_{i = 0}^{n-2}(d+(6i+4)c)(d+(6i+6)c)(d+(6i+8)c)}, \\ z_{6n-3} & = &\frac{ek^{3n-1}}{(a+k)(a+3k)\prod \limits_{i = 0}^{n-2}(a+(6i+5)k)(a+(6i+7)k)(a+(6i+9)k)}. \end{eqnarray*}$

It follows from Eq (1) that

$\begin{eqnarray*} x_{6n-2} & = &\dfrac{y_{6n-4}z_{6n-3}}{z_{6n-3}+x_{6n-5}} \\ & = &\tfrac{\left( \frac{k^{3n-1}}{(a+2k)\prod \limits_{i = 0}^{n-2}(a+(6i+4)k)(a+(6i+6)k)(a+(6i+8)k)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) \left( \frac{ ek^{3n-1}}{(a+k)(a+3k)\prod \limits_{i = 0}^{n-2}(a+(6i+5)k)(a+(6i+7)k)(a+(6i+9)k)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) }{\left( \frac{ ek^{3n-1}}{(a+k)(a+3k)\prod \limits_{i = 0}^{n-2}(a+(6i+5)k)(a+(6i+7)k)(a+(6i+9)k)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) +\left( \frac{ ek^{3n-2}}{(a+k)\prod\limits_{i = 0}^{n-2}(a+(6i+3)k)(a+(6i+5)k)(a+(6i+7)k)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) } \\ & = &\dfrac{\left( \frac{k^{3n}}{(a+2k)\prod \limits_{i = 0}^{n-2}(a+(6i+4)k)(a+(6i+6)k)(a+(6i+8)k)}\right) }{ (a+3k)\prod\limits_{i = 0}^{n-2}(a+(6i+9)k)\left[ \left( \frac{k}{ (a+3k)\prod\limits_{i = 0}^{n-2}(a+(6i+9)k)}\right) +\left( \frac{1}{ \prod\limits_{i = 0}^{n-2}(a+(6i+3)k)}\right) \right] } \\ & = &\dfrac{\left( \frac{k^{3n}}{(a+2k)\prod \limits_{i = 0}^{n-2}(a+(6i+4)k)(a+(6i+6)k)(a+(6i+8)k)}\right) }{\left[ k+\left( \frac{(a+3k)\prod\limits_{i = 0}^{n-2}(a+(6i+9)k)}{ \prod\limits_{i = 0}^{n-2}(a+(6i+3)k)}\right) \right] } \\ & = &\dfrac{\left( \frac{k^{3n}}{(a+2k)\prod \limits_{i = 0}^{n-2}(a+(6i+4)k)(a+(6i+6)k)(a+(6i+8)k)}\right) }{\left[ k+\left( a+(6n-3)k\right) \right] } \\ & = &\frac{ak^{3n}}{a(a+2k)\left( a+(6n-2)k\right) \prod\limits_{i = 0}^{n-2}(a+(6i+4)k)(a+(6i+6)k)(a+(6i+8)k)}. \end{eqnarray*}$

Then we see that

$\begin{equation*} x_{6n-2} = \frac{k^{3n}}{\prod \limits_{i = 0}^{n-1}(a+(6i)k)(a+(6i+2)k)(a+(6i+4)k)}. \end{equation*}$

Also, we see from Eq (1) that

$\begin{eqnarray*} y_{6n-2} & = &\dfrac{z_{6n-4}x_{6n-3}}{x_{6n-3}+y_{6n-5}} \\ & = &\tfrac{\left( \frac{c^{3n-1}}{(d+2c)\prod \limits_{i = 0}^{n-2}(d+(6i+4)c)(d+(6i+6)c)(d+(6i+8)c)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) \left( \frac{ hc^{3n-1}}{(d+c)(d+3c)\prod \limits_{i = 0}^{n-2}(d+(6i+5)c)(d+(6i+7)c)(d+(6i+9)c)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) }{\left( \frac{ hc^{3n-1}}{(d+c)(d+3c)\prod \limits_{i = 0}^{n-2}(d+(6i+5)c)(d+(6i+7)c)(d+(6i+9)c)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) +\left( \frac{ hc^{3n-2}}{(d+c)\prod\limits_{i = 0}^{n-2}(d+(6i+3)c)(d+(6i+5)c)(d+(6i+7)c)} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) } \\ & = &\dfrac{\left( \frac{c^{3n}}{(d+2c)\prod \limits_{i = 0}^{n-2}(d+(6i+4)c)(d+(6i+6)c)(d+(6i+8)c)}\right) }{ (d+3c)\prod\limits_{i = 0}^{n-2}(d+(6i+9)c)\left[ \left( \frac{c}{ (d+3c)\prod\limits_{i = 0}^{n-2}(d+(6i+9)c)}\right) +\left( \frac{1}{ \prod\limits_{i = 0}^{n-2}(d+(6i+3)c)}\right) \right] } \\ & = &\dfrac{\left( \frac{c^{3n}}{(d+2c)\prod \limits_{i = 0}^{n-2}(d+(6i+4)c)(d+(6i+6)c)(d+(6i+8)c)}\right) }{\left[ c+d+(6n-3)c\right] } \\ & = &\dfrac{c^{3n}}{\left[ d+(6n-2)c\right] (d+2c)\prod \limits_{i = 0}^{n-2}(d+(6i+4)c)(d+(6i+6)c)(d+(6i+8)c)}. \end{eqnarray*}$

Then

$\begin{equation*} y_{6n-2} = \dfrac{dc^{3n}}{\prod \limits_{i = 0}^{n-1}(d+(6i)c)(d+(6i+2)c)(d+(6i+4)c)}. \end{equation*}$

Finally from Eq (1), we see that

$\begin{eqnarray*} z_{6n-2} & = &\dfrac{x_{6n-4}y_{6n-3}}{y_{6n-3}+z_{6n-5}} \\ & = &\tfrac{\left( \frac{f^{3n-1}}{(g+2f)\prod \limits_{i = 0}^{n-2}(g+(6i+4)f)(g+(6i+6)f)(g+(6i+8)f)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) \left( \frac{ bf^{3n-1}}{(g+f)(g+3f)\prod \limits_{i = 0}^{n-2}(g+(6i+5)f)(g+(6i+7)f)(g+(6i+9)f)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) }{\left( \frac{ bf^{3n-1}}{(g+f)(g+3f)\prod \limits_{i = 0}^{n-2}(g+(6i+5)f)(g+(6i+7)f)(g+(6i+9)f)}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) +\left( \frac{ bf^{3n-2}}{(g+f)\prod\limits_{i = 0}^{n-2}(g+(6i+3)f)(g+(6i+5)f)(g+(6i+7)f)} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \right) } \\ & = &\dfrac{\left( \frac{f^{3n}}{(g+2f)\prod \limits_{i = 0}^{n-2}(g+(6i+4)f)(g+(6i+6)f)(g+(6i+8)f)}\right) }{ (g+3f)\prod\limits_{i = 0}^{n-2}(g+(6i+9)f)\left[ \left( \frac{f}{ (g+3f)\prod\limits_{i = 0}^{n-2}(g+(6i+9)f)}\right) +\left( \frac{1}{ \prod\limits_{i = 0}^{n-2}(g+(6i+3)f)}\right) \right] } \\ & = &\dfrac{\left( \frac{f^{3n}}{(g+2f)\prod \limits_{i = 0}^{n-2}(g+(6i+4)f)(g+(6i+6)f)(g+(6i+8)f)}\right) }{\left[ f+\left( \frac{(g+3f)\prod\limits_{i = 0}^{n-2}(g+(6i+9)f)}{ \prod\limits_{i = 0}^{n-2}(g+(6i+3)f)}\right) \right] } \\ & = &\dfrac{\left( \frac{f^{3n}}{(g+2f)\prod \limits_{i = 0}^{n-2}(g+(6i+4)f)(g+(6i+6)f)(g+(6i+8)f)}\right) }{\left[ f+(g+(6n-3)f)\right] } \\ & = &\dfrac{f^{3n}}{(g+(6n-2)f)(g+2f)\prod \limits_{i = 0}^{n-2}(g+(6i+4)f)(g+(6i+6)f)(g+(6i+8)f)}. \end{eqnarray*}$

Thus

$\begin{equation*} z_{3n-2} = \dfrac{gf^{3n}}{\prod \limits_{i = 0}^{n-1}(g+(6i)f)(g+(6i+2)f)(g+(6i+4)f)}. \end{equation*}$

By similar way, one can show the other relations. This completes the proof.

Lemma 1. Let $\{x_{n}, y_{n}, z_{n}\}$ be a positive solution of system (1), then all solution of (1) is bounded and approaching to zero.

Proof. It follows from Eq (1) that

$\begin{eqnarray*} x_{n+1} & = &\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}\leq y_{n-1}, \ \ \ \ \ y_{n+1} = \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}}\leq z_{n-1}, \\ z_{n+1} & = &\dfrac{x_{n-1}y_{n}}{y_{n}+z_{n-2}}\leq x_{n-1}, \end{eqnarray*}$

we see that

$\begin{eqnarray*} x_{n+4} &\leq &y_{n+2}, \ \ \ \ \ y_{n+2}\leq z_{n}, \ \ z_{n}\leq x_{n-2}, \ \ \ \Rightarrow \ x_{n+4} < x_{n-2}, \\ y_{n+4} &\leq &z_{n+2}, \ \ \ z_{n+2}\leq x_{n}, \ \ \ x_{n}\leq y_{n-2}, \ \ \ \ \ \Rightarrow \ y_{n+4} < y_{n-2}, \\ z_{n+4} &\leq &x_{n+2}, \ \ \ x_{n+2}\leq y_{n}, \ \ \ y_{n}\leq z_{n-2}, \ \ \ \ \ \Rightarrow \ z_{n+4} < z_{n-2}, \end{eqnarray*}$

Then all subsequences of $\{x_{n}, y_{n}, z_{n}\}$ (i.e., for $\{x_{n}\}$ are $\{x_{6n-2}\}$ , $\{x_{6n-1}\}$ , $\{x_{6n}\},$ $\{x_{6n+1}\}$ , $\{x_{6n+2}\}$ , $\{x_{6n+3}\}\$ are decreasing and at that time are bounded from above by $K, L$ and $M$ since $K = \max \{x_{-2}, x_{-1}, x_{0}, x_{1}, x_{2}, x_{3}\},$ $L = \max \{y_{-2}, y_{-1}, y_{0}, y_{1}, y_{2}, y_{3}\}$ and $M = \max \{z_{-2}, z_{-1}, z_{0}, z_{1}, z_{2}, z_{3}\}$ .

Example 1. We assume an interesting numerical example for the system (1) with $x_{-2} = -.22, \; x_{-1} = -.4, \ x_{0} = .12, \; y_{-2} = -.62,$ $y_{-1} = 4, \ y_{0} = .3, \; z_{-2} = .4, \; z_{-1} = .53\ and\; z_{0} = -2$ . (See Figure 1).

Figure 1. This figure shows the behavior of the solutions of the system (1) with the initial conditions

$x_{-2} = -.22, \; x_{-1} = -.4, \ x_{0} = .12, \; y_{-2} = -.62,$

$y_{-1} = 4, \ y_{0} = .3, \; z_{-2} = .4, \; z_{-1} = .53\ and\; z_{0} = -2.$ (We see from this figure that all solutions converges to zero).

DownLoad: Full-Size Img PowerPoint

3. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}$

In this section, we get the solution's form of the following system of difference equations

$\begin{equation} x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \;y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}-z_{n-2}}, \end{equation}$

(2)

where $n\in \mathbb{N}_{0}$ and the initial values are non-zero real numbers with $x_{-2}\neq \pm z_{0}, \neq -2z_{0},$ $z_{-2}\neq y_{0}, \neq 2y_{0}, \neq 3y_{0}$ and $y_{-2}\neq 2x_{0}, \neq \pm x_{0}.$

Theorem 2. Assume that $\{x_{n}, y_{n}, z_{n}\}$ are solutions of (2). Then for $n = 0, 1, 2, ...,$

$\begin{eqnarray*} x_{6n-2} & = &\frac{(-1)^{n}k^{3n}}{a^{2n-1}(a+2k)^{n}}, \ x_{6n-1} = \frac{ (-1)^{n}bf^{3n}}{(f-g)^{2n}(3f-g)^{n}}, \ x_{6n} = \frac{(-1)^{n}c^{3n+1}}{ d^{2n}(2c-d)^{n}}, \\ x_{6n+1} & = &\frac{ek^{3n+1}}{(a-k)^{n}(a+k)^{2n+1}}, \ x_{6n+2} = \frac{ (-1)^{n}f^{3n+2}}{g^{n}(2f-g)^{2n+1}}, \ x_{6n+3} = \frac{(-1)^{n}hc^{3n+2}}{ (c-d)^{2n+1}(c+d)^{n+1}}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-2} & = &\frac{(-1)^{n}c^{3n}}{d^{2n-1}(2c-d)^{n}}, \ y_{6n-1} = \frac{ ek^{3n}}{(a-k)^{n}(a+k)^{2n}}, \ y_{6n} = \frac{(-1)^{n}f^{3n+1}}{ g^{n}(2f-g)^{2n}}, \\ y_{6n+1} & = &\frac{(-1)^{n}hc^{3n+1}}{(c-d)^{2n}(c+d)^{n+1}}, \ y_{6n+2} = \frac{ (-1)^{n}k^{3n+2}}{a^{2n}(a+2k)^{n+1}}, \ y_{6n+3} = \frac{(-1)^{n}bf^{3n+2}}{ (f-g)^{2n+1}(3f-g)^{n+1}}, \end{eqnarray*}$

and

$\begin{eqnarray*} z_{6n-2} & = &\frac{(-1)^{n}f^{3n}}{g^{n-1}(2f-g)^{2n}}, \ z_{6n-1} = \frac{ (-1)^{n}hc^{3n}}{(c-d)^{2n}(c+d)^{n}}, \ z_{6n} = \frac{(-1)^{n}k^{3n+1}}{ a^{2n}(a+2k)^{n}}, \\ z_{6n+1} & = &\frac{(-1)^{n}bf^{3n+1}}{(f-g)^{2n+1}(3f-g)^{n}}, \ z_{6n+2} = \frac{(-1)^{n+1}c^{3n+2}}{d^{2n+1}(2c-d)^{n}}, \ z_{6n+3} = \frac{-ek^{3n+2}}{ (a-k)^{n}(a+k)^{2n+2}}, \end{eqnarray*}$

where $x_{-2} = a, \ x_{-1} = b,$ $x_{0} = c,$ $y_{-2} = d,$ $y_{-1} = e,$ $y_{0} = f,$ $z_{-2} = g,$ $z_{-1} = h$ $and$ $z_{0} = k.$

Proof. The result is true for $n = 0$ . Now suppose that $n > 0$ and that our claim verified for $n-1$ . That is,

$\begin{eqnarray*} x_{6n-8} & = &\frac{(-1)^{n-1}k^{3n-3}}{a^{2n-3}(a+2k)^{n-1}}, \ x_{6n-7} = \frac{ (-1)^{n-1}bf^{3n-3}}{(f-g)^{2n-2}(3f-g)^{n-1}}, \ x_{6n-6} = \frac{ (-1)^{n-1}c^{3n-2}}{d^{2n-2}(2c-d)^{n-1}}, \\ x_{6n-5} & = &\frac{ek^{3n-2}}{(a-k)^{n-1}(a+k)^{2n-1}}, \ x_{6n-4} = \frac{ (-1)^{n-1}f^{3n-1}}{g^{n-1}(2f-g)^{2n-1}}, \ x_{6n-3} = \frac{ (-1)^{n-1}hc^{3n-1}}{(c-d)^{2n-1}(c+d)^{n}}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-8} & = &\frac{(-1)^{n-1}c^{3n-3}}{d^{2n-3}(2c-d)^{n-1}}, \ y_{6n-7} = \frac{ ek^{3n-3}}{(a-k)^{n-1}(a+k)^{2n-2}}, \ y_{6n-6} = \frac{(-1)^{n-1}f^{3n-2}}{ g^{n-1}(2f-g)^{2n-2}}, \\ y_{6n-5} & = &\frac{(-1)^{n-1}hc^{3n-2}}{(c-d)^{2n-2}(c+d)^{n}}, \ y_{6n-4} = \frac{(-1)^{n-1}k^{3n-1}}{a^{2n-2}(a+2k)^{n}}, \ y_{6n-3} = \frac{ (-1)^{n-1}bf^{3n-1}}{(f-g)^{2n-1}(3f-g)^{n}}, \end{eqnarray*}$

and

$\begin{eqnarray*} z_{6n-8} & = &\frac{(-1)^{n-1}f^{3n-3}}{g^{n-2}(2f-g)^{2n-2}}, \ z_{6n-7} = \frac{ (-1)^{n-1}hc^{3n-3}}{(c-d)^{2n-2}(c+d)^{n-1}}, \ z_{6n-6} = \frac{ (-1)^{n-1}k^{3n-2}}{a^{2n-2}(a+2k)^{n-1}}, \\ z_{6n-5} & = &\frac{(-1)^{n-1}bf^{3n-2}}{(f-g)^{2n-1}(3f-g)^{n-1}}, \ z_{6n-4} = \frac{(-1)^{n}c^{3n-1}}{d^{2n-1}(2c-d)^{n-1}}, \ z_{6n-3} = \frac{-ek^{3n-1}}{ (a-k)^{n-1}(a+k)^{2n}}. \end{eqnarray*}$

Now from Eq (2), it follows that

$\begin{eqnarray*} x_{6n-2} & = &\dfrac{y_{6n-4}z_{6n-3}}{z_{6n-3}+x_{6n-5}} \\ & = &\dfrac{\left( \dfrac{(-1)^{n-1}k^{3n-1}}{a^{2n-2}(a+2k)^{n}}\right) \left( \dfrac{-ek^{3n-1}}{(a-k)^{n-1}(a+k)^{2n}}\right) }{\left( \dfrac{ -ek^{3n-1}}{(a-k)^{n-1}(a+k)^{2n}}\right) +\left( \dfrac{ek^{3n-2}}{ (a-k)^{n-1}(a+k)^{2n-1}}\right) } \\ & = &\dfrac{\left( \dfrac{(-1)^{n}k^{3n}}{a^{2n-2}(a+2k)^{n}}\right) }{\left( -k+a+k\right) } = \dfrac{(-1)^{n}k^{3n}}{a^{2n-1}(a+2k)^{n}}, \\ y_{6n-2} & = &\dfrac{z_{6n-4}x_{6n-3}}{x_{6n-3}+y_{6n-5}} = \dfrac{\left( \frac{ (-1)^{n}c^{3n-1}}{d^{2n-1}(2c-d)^{n-1}}\right) \left( \frac{ (-1)^{n-1}hc^{3n-1}}{(c-d)^{2n-1}(c+d)^{n}}\right) }{\left( \frac{ (-1)^{n-1}hc^{3n-1}}{(c-d)^{2n-1}(c+d)^{n}}\right) +\left( \frac{ (-1)^{n-1}hc^{3n-2}}{(c-d)^{2n-2}(c+d)^{n}}\right) } \\ & = &\dfrac{\left( \frac{(-1)^{n}c^{3n}}{d^{2n-1}(2c-d)^{n-1}}\right) }{c+c-d} = \frac{(-1)^{n}c^{3n}}{d^{2n-1}(2c-d)^{n}}, \\ z_{6n-2} & = &\dfrac{x_{6n-4}y_{6n-3}}{y_{6n-3}-z_{6n-5}} = \dfrac{\left( \frac{ (-1)^{n-1}f^{3n-1}}{g^{n-1}(2f-g)^{2n-1}}\right) \left( \frac{ (-1)^{n-1}bf^{3n-1}}{(f-g)^{2n-1}(3f-g)^{n}}\right) }{\left( \frac{ (-1)^{n-1}bf^{3n-1}}{(f-g)^{2n-1}(3f-g)^{n}}\right) -\left( \frac{ (-1)^{n-1}bf^{3n-2}}{(f-g)^{2n-1}(3f-g)^{n-1}}\right) } \\ & = &\dfrac{\left( \frac{(-1)^{n-1}f^{3n}}{g^{n-1}(2f-g)^{2n-1}}\right) }{ \left( f-3f+g\right) } = \frac{(-1)^{n}f^{3n}}{g^{n-1}(2f-g)^{2n}}. \end{eqnarray*}$

Also, we see from Eq (2) that

$\begin{eqnarray*} x_{6n-1} & = &\dfrac{y_{6n-3}z_{6n-2}}{z_{6n-2}+x_{6n-4}} \\ & = &\dfrac{\left( \frac{(-1)^{n-1}bf^{3n-1}}{(f-g)^{2n-1}(3f-g)^{n}}\right) \left( \frac{(-1)^{n}f^{3n}}{g^{n-1}(2f-g)^{2n}}\right) }{\left( \frac{ (-1)^{n}f^{3n}}{g^{n-1}(2f-g)^{2n}}\right) +\left( \frac{(-1)^{n-1}f^{3n-1}}{ g^{n-1}(2f-g)^{2n-1}}\right) } \\ & = &\dfrac{\left( \frac{(-1)^{n}bf^{3n}}{(f-g)^{2n-1}(3f-g)^{n}}\right) }{ \left( -f+2f-g\right) } = \frac{(-1)^{n}bf^{3n}}{(f-g)^{2n}(3f-g)^{n}}, \\ y_{6n-1} & = &\dfrac{z_{6n-3}x_{6n-2}}{x_{6n-2}+y_{6n-4}} = \dfrac{\left( \frac{ -ek^{3n-1}}{(a-k)^{n-1}(a+k)^{2n}}\right) \left( \frac{(-1)^{n}k^{3n}}{ a^{2n-1}(a+2k)^{n}}\right) }{\left( \frac{(-1)^{n}k^{3n}}{a^{2n-1}(a+2k)^{n}} \right) +\left( \frac{(-1)^{n-1}k^{3n-1}}{a^{2n-2}(a+2k)^{n}}\right) } \\ & = &\dfrac{\left( \frac{ek^{3n}}{(a-k)^{n-1}(a+k)^{2n}}\right) }{-k+a} = \frac{ ek^{3n}}{(a-k)^{n}(a+k)^{2n}}, \\ z_{6n-1} & = &\dfrac{x_{6n-3}y_{6n-2}}{y_{6n-2}-z_{6n-4}} = \dfrac{\left( \frac{ (-1)^{n-1}hc^{3n-1}}{(c-d)^{2n-1}(c+d)^{n}}\right) \left( \frac{ (-1)^{n}c^{3n}}{d^{2n-1}(2c-d)^{n}}\right) }{\left( \frac{(-1)^{n}c^{3n}}{ d^{2n-1}(2c-d)^{n}}\right) -\left( \frac{(-1)^{n}c^{3n-1}}{ d^{2n-1}(2c-d)^{n-1}}\right) } \\ & = &\dfrac{\left( \frac{(-1)^{n-1}hc^{3n}}{(c-d)^{2n-1}(c+d)^{n}}\right) }{ c-\left( 2c-d\right) } = \frac{(-1)^{n}hc^{3n}}{(c-d)^{2n}(c+d)^{n}}. \end{eqnarray*}$

Also, we can prove the other relations.

Example 2. See below for system (2) with the initial conditions $x_{-2} = 11, \; x_{-1} = 5, \ x_{0} = 13, \; y_{-2} = 6,$ $y_{-1} = 7, \ y_{0} = 3, \; z_{-2} = 14,$ $z_{-1} = 9$ $and\; z_{0} = 2.$

Figure 2. This figure shows the behavior of solutions of the systems of rational recursive sequence

$x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1} = \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n} }{y_{n}-z_{n-2}},$ when we take the initial conditions:

$x_{-2} = 11, \; x_{-1} = 5, \ x_{0} = 13, \; y_{-2} = 6,$

$y_{-1} = 7, \ y_{0} = 3, \; z_{-2} = 14,$

$z_{-1} = 9$

$and\; z_{0} = 2.$ (See the figure we can conclude that all the solutions unboundedness solutions).

DownLoad: Full-Size Img PowerPoint

4. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

Here, we obtain the form of solutions of the system

$\begin{equation} x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \;y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}-y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}+z_{n-2}}, \end{equation}$

(3)

where $n\in \mathbb{N}_{0}$ and the initial values are non-zero real numbers with $x_{-2}\neq \pm z_{0}, \neq 2z_{0},$ $z_{-2}\neq \pm y_{0}, \neq -2y_{0}$ and $y_{-2}\neq x_{0}, \neq 2x_{0}, \neq 3x_{0}.$

Theorem 3. If $\ \{x_{n}, y_{n}, z_{n}\}$ are solutions of system (3) where $x_{-2} = a, \ x_{-1} = b,$ $x_{0} = c,$ $y_{-2} = d,$ $y_{-1} = e,$ $y_{0} = f,$ $z_{-2} = g,$ $z_{-1} = h$ $and$ $z_{0} = k.$ Then for $n = 0, 1, 2, ...,$

$\begin{eqnarray*} x_{6n-2} & = &\frac{k^{3n}}{a^{2n-1}(a-2k)^{n}}, \ x_{6n-1} = \frac{ (-1)^{n}bf^{3n}}{(f-g)^{n}(f+g)^{2n}}, \ x_{6n} = \frac{(-1)^{n}c^{3n+1}}{ d^{n}(d-2c)^{2n}}, \\ x_{6n+1} & = &\frac{(-1)^{n}ek^{3n+1}}{(a-k)^{2n}(a+k)^{n+1}}, \ x_{6n+2} = \frac{ (-1)^{n}f^{3n+2}}{g^{2n}(2f+g)^{n+1}}, \ x_{6n+3} = \frac{(-1)^{n}hc^{3n+2}}{ (c-d)^{2n+1}(3c-d)^{n+1}}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-2} & = &\frac{(-1)^{n}c^{3n}}{d^{n-1}(d-2c)^{2n}}, \ y_{6n-1} = \frac{ (-1)^{n}ek^{3n}}{(a-k)^{2n}(a+k)^{n}}, \ y_{6n} = \frac{(-1)^{n}f^{3n+1}}{ g^{2n}(2f+g)^{n}}, \\ y_{6n+1} & = &\frac{(-1)^{n}hc^{3n+1}}{(c-d)^{2n+1}(3c-d)^{n}}, \ y_{6n+2} = \frac{-k^{3n+2}}{a^{2n+1}(a-2k)^{n}}, \ y_{6n+3} = \frac{(-1)^{n}bf^{3n+2}}{ (f-g)^{n}(f+g)^{2n+2}}, \end{eqnarray*}$

and

$\begin{eqnarray*} z_{6n-2} & = &\frac{(-1)^{n}f^{3n}}{g^{2n-1}(2f+g)^{n}}, \ z_{6n-1} = \frac{ (-1)^{n}hc^{3n}}{(c-d)^{2n}(3c-d)^{n}}, \ z_{6n} = \frac{k^{3n+1}}{ a^{2n}(a-2k)^{n}}, \\ z_{6n+1} & = &\frac{(-1)^{n}bf^{3n+1}}{(f-g)^{n}(f+g)^{2n+1}}, \ z_{6n+2} = \frac{ (-1)^{n}c^{3n+2}}{d^{n}(2c-d)^{2n+1}}, \ z_{6n+3} = \frac{(-1)^{n+1}ek^{3n+2}}{ (a-k)^{2n+1}(a+k)^{n+1}}. \end{eqnarray*}$

Proof. As the proof of Theorem 2 and so will be left to the reader.

Example 3. We put the initials $x_{-2} = 8, \; x_{-1} = 15, \ x_{0} = 13, \; y_{-2} = 6, \; y_{-1} = 7,$ $y_{0} = 3, \; z_{-2} = 14, \; z_{-1} = 19$ $and\; z_{0} = 2,$ for the system (3), see Figure 3.

Figure 3. This figure shows the unstable of the solutions of the difference equations system (3) with the initial values

$x_{-2} = 8, \; x_{-1} = 15, \ x_{0} = 13, \; y_{-2} = 6, \; y_{-1} = 7,$

$y_{0} = 3, \; z_{-2} = 14, \; z_{-1} = 19$

$and\; z_{0} = 2$ .

DownLoad: Full-Size Img PowerPoint

The following systems can be treated similarly.

5. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

In this section, we deal with the solutions of the following system

$\begin{equation} x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \;y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}+z_{n-2}}, \end{equation}$

(4)

where $n\in \mathbb{N}_{0}$ and the initial values are non-zero real with $x_{-2}\neq z_{0}, \neq 2z_{0}, \neq 3z_{0},$ $z_{-2}\neq \pm y_{0}, \neq 2y_{0}$ and $y_{-2}\neq \pm x_{0}, \neq -2x_{0}.$

Theorem 4. The solutions of system (4) are given by

$\begin{eqnarray*} x_{6n-2} & = &\frac{(-1)^{n}k^{3n}}{a^{n-1}(a-2k)^{2n}}, \ x_{6n-1} = \frac{ (-1)^{n}bf^{3n}}{(f-g)^{2n}(f+g)^{n}}, \ x_{6n} = \frac{(-1)^{n}c^{3n+1}}{ d^{2n}(d+2c)^{n}}, \\ x_{6n+1} & = &\frac{-ek^{3n+1}}{(a-k)^{2n+1}(a-3k)^{n}}, \ x_{6n+2} = \frac{ (-1)^{n+1}f^{3n+2}}{g^{2n+1}(2f-g)^{n}}, \ x_{6n+3} = \frac{(-1)^{n+1}hc^{3n+2} }{(c-d)^{n}(c+d)^{2n+2}}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-2} & = &\frac{(-1)^{n}c^{3n}}{d^{2n-1}(d+2c)^{n}}, \ y_{6n-1} = \frac{ ek^{3n}}{(a-k)^{2n}(a-3k)^{n}}, \ y_{6n} = \frac{(-1)^{n}f^{3n+1}}{ g^{2n}(2f-g)^{n}}, \\ y_{6n+1} & = &\frac{(-1)^{n}hc^{3n+1}}{(c+d)^{2n+1}(c-d)^{n}}, \ y_{6n+2} = \frac{ -k^{3n+2}}{a^{n}(a-2k)^{2n+1}}, \ y_{6n+3} = \frac{(-1)^{n}bf^{3n+2}}{ (f-g)^{2n+1}(f+g)^{n+1}}, \end{eqnarray*}$

and

$\begin{eqnarray*} z_{6n-2} & = &\frac{(-1)^{n}f^{3n}}{g^{2n-1}(2f-g)^{n}}, \ z_{6n-1} = \frac{ (-1)^{n}hc^{3n}}{(c+d)^{2n}(c-d)^{n}}, \ z_{6n} = \frac{(-1)^{n}k^{3n+1}}{ a^{n}(a-2k)^{2n}}, \\ z_{6n+1} & = &\frac{(-1)^{n}bf^{3n+1}}{(f-g)^{2n}(f+g)^{n+1}}, \ z_{6n+2} = \frac{ (-1)^{n}c^{3n+2}}{d^{2n}(2c+d)^{n+1}}, \ z_{6n+3} = \frac{ek^{3n+2}}{ (a-k)^{2n+1}(a-3k)^{n+1}}, \end{eqnarray*}$

where $x_{-2} = a, \ x_{-1} = b,$ $x_{0} = c,$ $y_{-2} = d,$ $y_{-1} = e,$ $y_{0} = f,$ $z_{-2} = g,$ $z_{-1} = h$ $and$ $z_{0} = k.$

Example 4. shows the behavior of the solution of system (4) with $x_{-2} = 18, \; x_{-1} = -15, \ x_{0} = 3, \; y_{-2} = 6,$ $y_{-1} = .7, \ y_{0} = -3,$ $z_{-2} = 4, \; z_{-1} = -9\ and\; z_{0} = 5$ .

Figure 4. This figure shows the behavior of the system

$x_{n+1} = \dfrac{ y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1} = \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}+z_{n-2}}$ with the initial conditions:-

$x_{-2} = 18, \; x_{-1} = -15, \ x_{0} = 3,$

$y_{-2} = 6,$

$y_{-1} = .7,$

$y_{0} = -3, \ z_{-2} = 4, \; z_{-1} = -9\ and\; z_{0} = 5.-0.6, \ x_{-1} = 0.2, \ x_{0} = -5.$ (From the figure, we see that all solutions goes to zero).

DownLoad: Full-Size Img PowerPoint

6. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}$

In this section, we obtain the solutions of the difference system

$\begin{equation} x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \;y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}-y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}-z_{n-2}}, \end{equation}$

(5)

where the initials are arbitrary non-zero real numbers with $x_{-2}\neq z_{0},$ $z_{-2}\neq y_{0}$ and $y_{-2}\neq x_{0}.$

Theorem 5. If $\ \{x_{n}, y_{n}, z_{n}\}$ are solutions of system (5) where $x_{-2} = a, \ x_{-1} = b,$ $x_{0} = c,$ $y_{-2} = d,$ $y_{-1} = e,$ $y_{0} = f,$ $z_{-2} = g,$ $z_{-1} = h$ $and$ $z_{0} = k$ . Then

$\begin{eqnarray*} x_{6n-2} & = &\frac{k^{3n}}{a^{3n-1}}, \ x_{6n-1} = \frac{bf^{3n}}{(f-g)^{3n}}, \ x_{6n} = \frac{c^{3n+1}}{d^{3n}}, \\ x_{6n+1} & = &\frac{ek^{3n+1}}{(k-a)^{3n+1}}, \ x_{6n+2} = \frac{f^{3n+2}}{ g^{3n+1}}, \ x_{6n+3} = \frac{hc^{3n+2}}{(c-d)^{3n+2}}, \end{eqnarray*}$

$\begin{eqnarray*} y_{6n-2} & = &\frac{c^{3n}}{d^{3n-1}}, \ y_{6n-1} = \frac{ek^{3n}}{(k-a)^{3n}}, \ y_{6n} = \frac{f^{3n+1}}{g^{3n}}, \\ y_{6n+1} & = &\frac{hc^{3n+1}}{(c-d)^{3n+1}}, \ y_{6n+2} = \frac{k^{3n+2}}{ a^{3n+1}}, \ y_{6n+3} = \frac{bf^{3n+2}}{(f-g)^{3n+2}}, \end{eqnarray*}$

and

$\begin{eqnarray*} z_{6n-2} & = &\frac{f^{3n}}{g^{3n-1}}, \ z_{6n-1} = \frac{hc^{3n}}{(c-d)^{3n}}, \ z_{6n} = \frac{k^{3n+1}}{a^{3n}}, \\ z_{6n+1} & = &\frac{bf^{3n+1}}{(f-g)^{3n+1}}, \ z_{6n+2} = \frac{c^{3n+2}}{ d^{3n+1}}, \ z_{6n+3} = \frac{ek^{3n+2}}{(k-a)^{3n+2}}. \end{eqnarray*}$

Example 5. shows the dynamics of the solution of system (5) with $x_{-2} = 18, \; x_{-1} = -15, \; x_{0} = 3, y_{-2} = 6, y_{-1} = .7,$ $y_{0} = -3, z_{-2} = 4, z_{-1} = -9\ and\; z_{0} = 5.$

Figure 5. This figure shows the behavior of the system of nonlinear difference equations (5) with the initial conditions considered as follows:-

$x_{-2} = 18, \; x_{-1} = -15, \ x_{0} = 3, \; y_{-2} = 6,$

$y_{-1} = .7, \ y_{0} = -3, \; z_{-2} = 4, \; z_{-1} = -9\ and\; z_{0} = 5$ .

DownLoad: Full-Size Img PowerPoint

7. Conclusions

This paper discussed the expression's form and boundedness of some systems of rational third order difference equations. In Section 2, we studied the qualitative behavior of system $x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1} = \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n} }{y_{n}+z_{n-2}},$ first we have got the form of the solutions of this system, studied the boundedness and gave numerical example and drew it by using Matlab. In Section 3, we have got the solution's of the system $x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}-z_{n-2}},$ and take a numerical example. In Sections 4–6, we obtained the solution of the following systems respectively, $x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}-y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}+z_{n-2}},$ $x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}+y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}+z_{n-2}},$ and $x_{n+1} = \dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1} = \dfrac{z_{n-1}x_{n}}{ x_{n}-y_{n-2}}, \ z_{n+1} = \dfrac{x_{n-1}y_{n}}{y_{n}-z_{n-2}}.$ Also, in each case we take a numerical example to illustrates the results.

Acknowledgements

This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under grant no. (G: 233–130–1441). The authors, therefore, acknowledge with thanks DSR for technical and financial support.

Conflict of interest

All authors declare no conflicts of interest in this paper.

References

[1]	M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Progr. Fract. Differ. Appl., 1 (2015), 73–85. http://dx.doi.org/10.12785/pfda/010201 doi: 10.12785/pfda/010201
[2]	A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model, Therm. Sci., 20 (2016), 763–769. https://doi.org/10.2298/TSCI160111018A doi: 10.2298/TSCI160111018A
[3]	X. J. Yang, H. M. Srivastava, J. A. Tenreiro Machado, A new fractional derivative without singular kernel: Application to the modelling of the steady heat flow, Therm. Sci., 20 (2016), 753–756. https://doi.org/10.2298/TSCI151224222Y doi: 10.2298/TSCI151224222Y
[4]	S. T. Sutar, K. D. Kucche, On nonlinear hybrid fractional diferential equations with Atangana-Baleanu-Caputo derivative, Chaos Soliton. Fract., 143 (2021), 110557. https://doi.org/10.1016/j.chaos.2020.110557 doi: 10.1016/j.chaos.2020.110557
[5]	X. M. Gu, H. W. Sun, Y. L. Zhao, X. Zheng, An implicit difference scheme for time-fractional diffusion equations with a time-invariant type variable order, Appl. Math. Lett., 120 (2021), 107270. https://doi.org/10.1016/j.aml.2021.107270 doi: 10.1016/j.aml.2021.107270
[6]	M. Hassouna, E. H. El Kinani, A. Ouhadan, Global existence and uniqueness of solution of Atangana-Baleanu-Caputo fractional differential equation with nonlinear term and approximate solutions, Int. J. Differ. Equations, 2021 (2021), 5675789. https://doi.org/10.1155/2021/5675789 doi: 10.1155/2021/5675789
[7]	J. Gómez-Aguilar, R. Escobar-Jiménez, M. López-López, V. Alvarado-Martínez, Atangana-Baleanu fractional derivative applied to electromagnetic waves in dielectric media, J. Electromagn. Waves Appl., 30 (2016), 1937–1952. https://doi.org/10.1080/09205071.2016.1225521 doi: 10.1080/09205071.2016.1225521
[8]	S. Ullah, M. A. Khan, M. Farooq, Modeling and analysis of the fractional HBV model with Atangana-Baleanu derivative, Eur. Phys. J. Plus, 133 (2018), 313. https://doi.org/10.1140/epjp/i2018-12120-1 doi: 10.1140/epjp/i2018-12120-1
[9]	O. J. Peter, A. S. Shaikh, M. O. Ibrahim, K. S. Nisar, D. Baleanu, I. Khan, et al., Analysis and dynamics of fractional order mathematical model of covid-19 in Nigeria using Atangana-Baleanu operator, Comput. Mater. Con., 66 (2020), 1823–1848. http://dx.doi.org/10.32604/cmc.2020.012314 doi: 10.32604/cmc.2020.012314
[10]	C. F. Lorenzo, T. T. Hartley, Variable order and distributed order fractional operators, Nonlinear Dyn., 29 (2002), 57–98. https://doi.org/10.1023/A:1016586905654 doi: 10.1023/A:1016586905654
[11]	H. Sun, W. Chen, H. Wei, A comparative study of constant-order and variable-order fractional models in characterizing memory property of systems, Eur. Phys. J. Spec. Top., 193 (2011), 185. https://doi.org/10.1140/epjst/e2011-01390-6 doi: 10.1140/epjst/e2011-01390-6
[12]	M. A. Abdelkawy, M. A. Zaky, A. H. Bhrawy, D. Baleanu, Numerical simulation of time variable fractional order mobile-immobile advection-dispersion model, Rom. Rep. Phys., 67 (2015), 773–791.
[13]	C. Canuto, M. Y. Hussaini, A. Quarteroni, T. A. Zang, Spectral methods in fluid dynamics, Springer, Berlin, 1987.
[14]	J. Solís-Pérez, J. Gómez-Aguilar, A. Atangana, Novel numerical method for solving variable-order fractional differential equations with power, exponential and Mittag-Leffler laws, Chaos Soliton. Fract., 114 (2018), 175–185. https://doi.org/10.1016/j.chaos.2018.06.032 doi: 10.1016/j.chaos.2018.06.032
[15]	X. Li, Y. Gao, B. Wu, Approximate solutions of Atangana-Baleanu variable order fractional problems, AIMS Math., 5 (2020), 2285–2294. https://doi.org/10.3934/math.2020151 doi: 10.3934/math.2020151
[16]	M. H. Heydari, Z. Avazzadeh, A. Atangana, Shifted Jacobi polynomials for nonlinear singular variable-order time fractional Emden-Fowler equation generated by derivative with non-singular kernel, Adv. Differ. Equations, 2021 (2021), 188. https://doi.org/10.1186/s13662-021-03349-1 doi: 10.1186/s13662-021-03349-1
[17]	A. H. Bhrawy, M. A. Zaky, Numerical simulation for two-dimensional variable-order fractional nonlinear cable equation, Nonlinear Dyn., 80 (2015), 101–116. https://doi.org/10.1007/s11071-014-1854-7 doi: 10.1007/s11071-014-1854-7
[18]	T. El-Gindy, H. Ahmed, M. Melad, Shifted Gegenbauer operational matrix and its applications for solving fractional differential equations, J. Egypt. Math. Soc., 26 (2018), 72–90. https://doi.org/10.21608/JOMES.2018.9463 doi: 10.21608/JOMES.2018.9463
[19]	M. Usman, M. Hamid, T. Zubair, R. U. Haq, W. Wang, M. Liu, Novel operational matrices-based method for solving fractional-order delay differential equations via shifted Gegenbauer polynomials, Appl. Math. Comput., 372 (2020), 124985. https://doi.org/10.1016/j.amc.2019.124985 doi: 10.1016/j.amc.2019.124985
[20]	F. Soufivand, F. Soltanian, K. Mamehrashi, An operational matrix method based on the Gegenbauer polynomials for solving a class of fractional optimal control problems, Int. J. Industrial Electron. Control Optim., 4 (2021), 475–484. https://doi.org/10.22111/IECO.2021.39546.1371 doi: 10.22111/IECO.2021.39546.1371
[21]	M. Heydari, A. Atangana, A cardinal approach for nonlinear variable-order time fractional schr $\ddot{o}$ dinger equation defined by Atangana-Baleanu-Caputo derivative, Chaos Soliton. Fract., 128 (2019), 339–348. https://doi.org/10.1016/j.chaos.2019.08.009 doi: 10.1016/j.chaos.2019.08.009
[22]	Y. V. Mukhartova, M. A. Davydova, N. F. Elansky, O. V. Postylyakov, S. A. Zakharova, A. N. Borovski, On application of nonlinear reaction-diffusion-advection models to simulation of transport of chemically-active impurities, Remote Sensing Technologies and Applications in Urban Environments IV, 11157 (2019), 180–187. https://doi.org/10.1117/12.2535489 doi: 10.1117/12.2535489
[23]	F. Heße, F. A. Radu, M. Thullner, S. Attinger, Upscaling of the advection-diffusion-reaction equation with Monod reaction, Adv. Water Resour., 32 (2009), 1336–1351. https://doi.org/10.1016/j.advwatres.2009.05.009 doi: 10.1016/j.advwatres.2009.05.009
[24]	A. Hamdi, Identification of point sources in two-dimensional advection-diffusion-reaction equation: Application to pollution sources in a river. Stationary case, Inverse Probl. Sci. Eng., 15 (2007), 855–870. https://doi.org/10.1080/17415970601162198 doi: 10.1080/17415970601162198
[25]	A. Rubio, A. Zalts, C. El Hasi, Numerical solution of the advection-reaction-diffusion equation at different scales, Environ. Modell. Softw., 23 (2008), 90–95. https://doi.org/10.1016/j.envsoft.2007.05.009 doi: 10.1016/j.envsoft.2007.05.009
[26]	K. Issa, B. M. Yisa, J. Biazar, Numerical solution of space fractional diffusion equation using shifted Gegenbauer polynomials, Comput. Methods Differ. Equations, 10 (2022), 431–444. https://dx.doi.org/10.22034/cmde.2020.42106.1818 doi: 10.22034/cmde.2020.42106.1818
[27]	U. Ali, A. Iqbal, M. Sohail, F. A. Abdullah, Z. Khan, Compact implicit difference approximation for time-fractional diffusion-wave equation, Alex. Eng. J., 61 (2022), 4119–4126. https://doi.org/10.1016/j.aej.2021.09.005 doi: 10.1016/j.aej.2021.09.005
[28]	M. A. Zaky, S. S. Ezz-Eldien, E. H. Doha, J. A. Tenreiro Machado, A. H. Bhrawy, An efficient operational matrix technique for multidimensional variable-Order time fractional diffusion equations, ASME J. Comput. Nonlinear Dyn., 11 (2016), 061002. https://doi.org/10.1115/1.4033723 doi: 10.1115/1.4033723
[29]	M. M. Izadkhah, J. Saberi-Nadjafi, Gegenbauer spectral method for time-fractional convection-difffusion equations with variable coefficients, Math. Methods Appl. Sci., 38 (2015), 3183–3194. https://doi.org/10.1002/mma.3289 doi: 10.1002/mma.3289
[30]	M. H. Heydari, A. Atangana, Z. Avazzadeh, M. R. Mahmoudi, An operational matrix method for nonlinear variable-order time fractional reaction-diffusion equation involving Mittag-Leffler kernel, Eur. Phys. J. Plus, 135 (2020), 237. https://doi.org/10.1140/epjp/s13360-020-00158-5 doi: 10.1140/epjp/s13360-020-00158-5
[31]	P. Pandey, S. Kumar, J. G $\acute{o}$ mez-Aguilar, Numerical solution of the time fractional reaction-advection-diffusion equation in porous media, J. Appl. Comput. Mech., 8 (2022), 84–96. https://doi.org/10.22055/JACM.2019.30946.1796 doi: 10.22055/JACM.2019.30946.1796
[32]	S. Kumar, D. Zeidan, An efficient Mittag-Leffler kernel approach for time-fractional advection-reaction-diffusion equation, Appl. Numer. Math., 170 (2021), 190–207. https://doi.org/10.1016/j.apnum.2021.07.025 doi: 10.1016/j.apnum.2021.07.025
[33]	K. S. Miller, B. Ross, An introduction to the fractional calculus and fractional differential equations, New York: Wiley, 1993.
[34]	M. Hosseininia, M. H. Heydari, Legendre wavelets for the numerical solution of nonlinear variable-order time fractional 2d reaction-diffusion equation involving Mittag-Leffler non-singular kernel, Chaos Soliton. Fract., 127 (2019), 400–407. https://doi.org/10.1016/j.chaos.2019.07.017 doi: 10.1016/j.chaos.2019.07.017
[35]	F. R. Lin, H. Qu, A Runge-Kutta Gegenbauer spectral method for nonlinear fractional differential equations with Riesz fractional derivatives, Int. J. Comput. Math., 96 (2018), 417–435. https://doi.org/10.1080/00207160.2018.1487059 doi: 10.1080/00207160.2018.1487059
[36]	H. Tajadodi, A numerical approach of fractional advection-diffusion equation with Atangana-Baleanu derivative, Chaos Soliton. Fract., 130 (2020), 109527. https://doi.org/10.1016/j.chaos.2019.109527 doi: 10.1016/j.chaos.2019.109527
[37]	S. Yadav, R. K. Pandey, Numerical approximation of fractional Burgers equation with Atangana-Baleanu derivative in Caputo sense, Chaos Soliton. Fract., 133 (2020), 109630. https://doi.org/10.1016/j.chaos.2020.109630 doi: 10.1016/j.chaos.2020.109630

This article has been cited by:

1.	Khalil S. Al-Basyouni, Elsayed M. Elsayed, On Some Solvable Systems of Some Rational Difference Equations of Third Order, 2023, 11, 2227-7390, 1047, 10.3390/math11041047
2.	Ibraheem M. Alsulami, E. M. Elsayed, On a class of nonlinear rational systems of difference equations, 2023, 8, 2473-6988, 15466, 10.3934/math.2023789
3.	E.M. Elsayed, B.S. Alofi, The periodic nature and expression on solutions of some rational systems of difference equations, 2023, 74, 11100168, 269, 10.1016/j.aej.2023.05.026
4.	Hashem Althagafi, Dynamics of difference systems: a mathematical study with applications to neural systems, 2025, 10, 2473-6988, 2869, 10.3934/math.2025134

Reader Comments

Your name:*

Email:*
© 2022 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(1683) PDF downloads(70) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(13) / Tables(5)

AIMS Mathematics

Numerical approximation of a variable-order time fractional advection-reaction-diffusion model via shifted Gegenbauer polynomials

Related Papers:

Abstract

1. Introduction

2. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

3. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}$

4. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

5. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

6. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}$

7. Conclusions

Acknowledgements

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

Numerical approximation of a variable-order time fractional advection-reaction-diffusion model via shifted Gegenbauer polynomials

Related Papers:

Abstract

1. Introduction

2. The system: xn+1=yn−1znzn+xn−2,yn+1=zn−1xnxn+yn−2, x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}}, zn+1=xn−1ynyn+zn−2z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}

3. The system: xn+1=yn−1znzn+xn−2,yn+1=zn−1xnxn+yn−2, x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}}, zn+1=xn−1ynyn−zn−2z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}

4. The system: xn+1=yn−1znzn+xn−2,yn+1=zn−1xnxn−yn−2, x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}}, zn+1=xn−1ynyn+zn−2z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}

5. The system: xn+1=yn−1znzn−xn−2,yn+1=zn−1xnxn+yn−2, x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}}, zn+1=xn−1ynyn+zn−2z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}

6. The system: xn+1=yn−1znzn−xn−2,yn+1=zn−1xnxn−yn−2, x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}}, zn+1=xn−1ynyn−zn−2z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}

7. Conclusions

Acknowledgements

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

2. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

3. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}$

4. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}+x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

5. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}+y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}+z_{n-2}}$

6. The system: $x_{n+1}=\dfrac{y_{n-1}z_{n}}{z_{n}-x_{n-2}}, \; y_{n+1}= \dfrac{z_{n-1}x_{n}}{x_{n}-y_{n-2}},$ $z_{n+1}=\dfrac{x_{n-1}y_{n}}{ y_{n}-z_{n-2}}$