# On the Solutions of Three-Dimensional Rational Difference Equation Systems

Authors: H. S. Alayachi; A. Q. Khan; M. S. M. Noorani
Journal: Journal of Mathematics
 (2021)

---

## Abstract

In this paper, we are interested in a technique for solving some nonlinear rational systems of difference equations of third order, in three-dimensional case as a special case of the following system:xn+1=ynzn−1/yn±xn−2,yn+1=znxn−1/zn±yn−2,andzn+1=xnyn−1/xn±zn−2 with initial conditions x−2,x−1,x0,y−2, y−1,y0,z−2,z−1,andz0 are nonzero real numbers. Moreover, we study some behavior of the systems such as the boundedness of solutions for such systems. Finally, we present some numerical examples by giving some numerical values for the initial values of each case. Some figures have been given to explain the behavior of the obtained solutions in the case of numerical examples by using the mathematical program MATLAB to confirm the obtained results.


---

## Body

## 1. Introduction

We believe that difference equations, also referred to as recursive sequence, are a hot topic here as there has been increasing interest in the study of qualitative analysis of difference equations and systems of difference equations. Difference equations appear naturally as discrete analogues and as numerical solutions of differential and delay differential equations having applications in biology, ecology, economics, physics, computer sciences, and so on. Although difference equations are very simple in form, it is extremely difficult to understand thoroughly the behaviors of their solution (see [1–9] and the references cited therein). Recently, a great effort has been made in studying the qualitative analysis of rational difference equations and rational difference system (see [10–35]).The study of the nonlinear rational difference equations is quite challenging and rewarding [2, 8]. The results about these equations offer prototypes towards the development of the basic theory of the global behavior of nonlinear difference equations [1]. Recently, many researchers have investigated the behavior of the solution of difference equations. Difference equations arise in the situations in which the discrete values of the independent variable involve. Many practical phenomena are modeled with the help of difference equations [1]. In engineering, difference equations arise in control engineering, digital signal processing, electrical networks, etc. In social sciences, difference equations arise to study the national income of a country and then its variation with time, Cobweb phenomenon in economics, etc.There are many papers related to the difference equation system, for example, the periodicity of the positive solutions of the rational difference equations system:(1)xn+1=1zn,yn+1=ynxn−1yn−1,zn+1=1xn−1,has been obtained by Cinar in [4].Khan et al. [6] studied the equilibrium points, local asymptotic stability of an equilibrium point, instability of equilibrium points, periodicity behavior of positive solutions, and global character of an equilibrium point of a fourth-order system of rational difference equations of the form:(2)xn+1=αxn−3β+γynyn−1yn−2yn−3,yn+1=α1yn−3β1+γ1xnxn−1xn−2xn−3.Elabbasy et al. [7] has obtained the solution of particular cases of the following general system of difference equations:(3)xn+1=a1+a2yna3zn+a4xn−1zn,yn+1=b1zn−1+b2znb3xnyn+b4xnyn−1,zn+1=c1zn−1+c2znc3xn−1yn−1+c4xn−1yn+c5xnyn.In [36], Elsayed et al. dealt with the solutions of the systems of the difference equations:(4)xn+1=1xn−pyn−p,yn+1=xn−pyn−pxn−qyn−q,xn+1=1xn−pyn−pzn−p,yn+1=xn−pyn−pzn−pxn−qyn−qzn−q,zn+1=xn−qyn−qzn−qxn−ryn−rzn−r.Kurbanli [13–15] investigated the behavior of the solutions of the difference equation systems:(5)xn+1=xn−1xn−1yn−1,yn+1=yn−1yn−1xn−1,zn+1=1znyn,xn+1=xn−1xn−1yn−1,yn+1=yn−1yn−1xn−1,zn+1=zn−1zn−1yn−1,xn+1=xn−1xn−1yn−1,yn+1=yn−1yn−1xn−1,zn+1=xnzn−1yn.In [21], Yalcinkaya et al. studied the periodic character of the following two systems of difference equations:(6)xn+11=xn2xn2−1,xn+12=xn3xn3−1,…,xn+1k=xn1xn1−1,xn+11=xnkxnk−1,xn+12=xn1xn1−1,…,xn+1k=xnk−1xnk−1−1,where the initial values are nonzero real numbers for x01,x02,…,x0k≠1.In [37], Zhang et al. studied the boundedness, the persistence, and global asymptotic stability of the positive solutions of the system of difference equations:(7)xn=A+1yn−p,yn=A+yn−1xn−ryn−s.Zkan and Kurbanli [38] have investigated the periodical solutions of the following system of third-order rational difference equations:(8)xn+1=yn−2−1±yn−2xn−1yn,yn+1=xn−2−1±xn−2yn−1xn,zn+1=xn−2+yn−2−1±xn−2yn−1xn.Similar to the references above, this paper is devoted to obtain the form of the solutions of the following third-order systems of rational difference equations:(9)xn+1=ynzn−1yn±xn−2,yn+1=znxn−1zn±yn−2,zn+1=xnyn−1xn±zn−2,where the initial conditions x−2,x−1,x0,y−2, y−1,y0,z−2,z−1,andz0 are nonzero real numbers.


## 2. The System:xn+1=ynzn−1/yn+xn−2,yn+1=znxn−1/zn+yn−2,andzn+1=xnyn−1/xn+zn−2

We obtain the form of solutions of the following system of difference equations:(10)xn+1=ynzn−1yn+xn−2,yn+1=znxn−1zn+yn−2,zn+1=xnyn−1xn+zn−2.Theorem 1.
Ifxn,yn,zn be the solution of (10), then(11)x6n−2=x−2y03n∏i=0n−1x−2+6iy0x−2+6i+2y0x−2+6i+4y0,x6n−1=x−1z03n∏i=0n−1y−2+6i+1z0y−2+6i+3z0y−2+6i+5z0,x6n=x03n+1∏i=0n−1z−2+6i+2x0z−2+6i+4x0z−2+6i+6x0,x6n+1=z−1y03n+1x−2+y0∏i=0n−1x−2+6i+3y0x−2+6i+5y0x−2+6i+7y0,x6n+2=z03n+2y−2+2z0∏i=0n−1y−2+6i+4z0y−2+6i+6z0y−2+6i+8z0,x6n+3=y−1x03n+2z−2+x0z−2+3x0∏i=0n−1z−2+6i+5x0z−2+6i+7x0z−2+6i+9x0,y6n−2=y−2z03n∏i=0n−1y−2+6iz0y−2+6i+2z0y−2+6i+4z0,y6n−1=y−1x03n∏i=0n−1z−2+6i+1x0z−2+6i+3x0z−2+6i+5x0,y6n=y03n+1∏i=0n−1x−2+6i+2y0x−2+6i+4y0x−2+6i+6y0,y6n+1=z−1z03n+1y−2+z0∏i=0n−1y−2+6i+3z0y−2+6i+5z0y−2+6i+7z0,y6n+2=x03n+2z−2+2x0∏i=0n−1z−2+6i+4x0z−2+6i+6x0z−2+6i+8x0,y6n+3=z−1y03n+2x−2+y0x−2+3y0∏i=0n−1x−2+6i+5y0x−2+6i+7y0x−2+6i+9y0,z6n−2=z−2x03n∏i=0n−1z−2+6ix0z−2+6i+2x0z−2+6i+4x0,z6n−1=z−1y03n∏i=0n−1x−2+6i+1y0x−2+6i+3y0x−2+6i+5y0,z6n=z03n+1∏i=0n−1y−2+6i+2z0y−2+6i+4z0y−2+6i+6z0,z6n+1=y−1x03n+1z−2+x0∏i=0n−1z−2+6i+3x0z−2+6i+5x0z−2+6i+7x0,z6n+2=y03n+2x−2+2y0∏i=0n−1x−2+6i+4y0x−2+6i+6y0x−2+6i+8y0,z6n+3=x−1z03n+2y−2+z0y−2+3z0∏i=0n−1y−2+6i+5z0y−2+6i+7z0y−2+6i+9z0.Proof.
Obviously, results are true forn=0. Suppose that it is also true for n−1, i.e.,(12)x6n−8=x−2y03n−3∏i=0n−2x−2+6iy0x−2+6i+2y0x−2+6i+4y0,x6n−7=x−1z03n−3∏i=0n−2y−2+6i+1z0y−2+6i+3z0y−2+6i+5z0,x6n−6=x03n−2∏i=0n−2z−2+6i+2x0z−2+6i+4x0z−2+6i+6x0,x6n−5=z−1y03n−2x−2+y0∏i=0n−2x−2+6i+3y0x−2+6i+5y0x−2+6i+7y0,x6n−4=z03n−1y−2+2z0∏i=0n−2y−2+6i+4z0y−2+6i+6z0y−2+6i+8z0,x6n−3=y−1x03n−1z−2+x0z−2+3x0∏i=0n−2z−2+6i+5x0z−2+6i+7x0z−2+6i+9x0,y6n−8=y−2z03n−3∏i=0n−2y−2+6iz0y−2+6i+2z0y−2+6i+4z0,y6n−7=y−1x03n−3∏i=0n−2z−2+6i+1x0z−2+6i+3x0z−2+6i+5x0,y6n−6=y03n−2∏i=0n−2x−2+6i+2y0x−2+6i+4y0x−2+6i+6y0,y6n−5=x−1z03n−2y−2+z0∏i=0n−2y−2+6i+3z0y−2+6i+5z0y−2+6i+7z0,y6n−4=x03n−1z−2+2x0∏i=0n−2z−2+6i+4x0z−2+6i+6x0z−2+6i+8x0,y6n−3=z−1y03n−1x−2+y0x−2+3y0∏i=0n−2x−2+6i+5y0x−2+6i+7y0x−2+6i+9y0,z6n−8=z−2x03n−3∏i=0n−2z−2+6ix0z−2+6i+2x0z−2+6i+4x0,z6n−7=z−1y03n−3∏i=0n−2x−2+6i+1y0x−2+6i+3y0x−2+6i+5y0,z6n−6=z03n−2∏i=0n−2y−2+6i+2z0y−2+6i+4z0y−2+6i+6z0,z6n−5=y−1x03n−2z−2+x0∏i=0n−2z−2+6i+3x0z−2+6i+5x0z−2+6i+7x0,z6n−4=y03n−1x−2+2y0∏i=0n−2x−2+6i+4y0x−2+6i+6y0x−2+6i+8y0,z6n−3=x−1z03n−1y−2+z0y−2+3z0∏i=0n−2y−2+6i+5z0y−2+6i+7z0y−2+6i+9z0.
From equation (10), it follows that(13)x6n−2=z6n−4y6n−3y6n−3+x6n−5=y03n−1/x−2+2y0∏i=0n−2x−2+6i+4y0x−2+6i+6y0x−2+6i+8y0z−1y03n−1/x−2+y0x−2+3y0∏i=0n−2x−2+6i+5y0x−2+6i+7y0x−2+6i+9y0z−1y03n−1/x−2+y0x−2+3y0∏i=0n−2x−2+6i+5y0x−2+6i+7y0x−2+6i+9y0+z−1y03n−2/x−2+y0∏i=0n−2x−2+6i+3y0x−2+6i+5y0x−2+6i+7y0=y03n/x−2+2y0∏i=0n−2x−2+6i+4y0x−2+6i+6y0x−2+6i+8y0y0+x−2+3y0∏i=0n−2x−2+6i+9y0/∏i=0n−2x−2+6i+3y0=y03n/x−2+2y0∏i=0n−2x−2+6i+4y0x−2+6i+6y0x−2+6i+8y0y0+x−2+6n−3y0=y03nx−2+2y0∏i=0n−2x−2+6i+4y0x−2+6i+6y0x−2+6i+8y01x−2+6n−2y0.
Now, one can see that(14)x6n−2=x−2y03n∏i=0n−1x−2+6iy0x−2+6i+2y0x−2+6i+4y0.
Moreover, from (10), one has(15)y6n−2=x6n−4z6n−3z6n−3+y6n−5=z03n−1/y−2+2z0∏i=0n−2y−2+6i+4z0y−2+6i+6z0y−2+6i+8z0x−1z03n−1/y−2+z0y−2+3z0∏i=0n−2y−2+6i+5z0y−2+6i+7z0y−2+6i+9z0x−1z03n−1/y−2+z0y−2+3z0∏i=0n−2y−2+6i+5z0y−2+6i+7z0y−2+6i+9z0+x−1z03n−2/y−2+z0∏i=0n−2y−2+6i+3z0y−2+6i+5z0y−2+6i+7z0=z03n/y−2+2z0∏i=0n−2y−2+6i+4z0y−2+6i+6z0y−2+6i+8z0z0+y−2+3z0∏i=0n−2y−2+6i+9z0/∏i=0n−2y−2+6i+3z0=z03n/y−2+2z0∏i=0n−2y−2+6i+4z0y−2+6i+6z0y−2+6i+8z0z0+y−2+6n−3z0=z03n/y−2+2z0∏i=0n−2y−2+6i+4z0y−2+6i+6z0y−2+6i+8z0y−2+6n−2z0.
Then,(16)y6n−2=y−2z03ny−2+2z0∏i=0n−1y−2+6iz0y−2+6i+2z0y−2+6i+4z0.
Finally, from equation (10),(17)z6n−2=y6n−4x6n−3x6n−3+z6n−5=x03n−1/z−2+2x0∏i=0n−2z−2+6i+4x0z−2+6i+6x0z−2+6i+8x0y−1x03n−1/z−2+x0z−2+3x0∏i=0n−2z−2+6i+5x0z−2+6i+7x0z−2+6i+9x0y−1x03n−1/z−2+x0z−2+3x0∏i=0n−2z−2+6i+5x0z−2+6i+7x0z−2+6i+9x0+y−1x03n−2/z−2+x0∏i=0n−2z−2+6i+3x0z−2+6i+5x0z−2+6i+7x0=x03n/z−2+2x0∏i=0n−2z−2+6i+4x0z−2+6i+6x0z−2+6i+8x0x0+z−2+3x0∏i=0n−2z−2+6i+9x0/∏i=0n−2z−2+6i+3x0=x03n/z−2+2x0∏i=0n−2z−2+6i+4x0z−2+6i+6x0z−2+6i+8x0x0+z−2+6n−3x0.
Thus,(18)z3n−2=z−2x03n∏i=0n−1z−2+6ix0z−2+6i+2x0z−2+6i+4x0.
In a similar way, other relations can also be proved.Lemma 1.
Ifxn,yn,zn is the +ve solution of (10), then it is bounded as well as converge to zero.Proof.
From (10),(19)xn+1=ynzn−1yn+xn−2≤zn−1,yn+1=znxn−1zn+yn−2≤xn−1,zn+1=xnyn−1xn+zn−2≤yn−1,we see that(20)xn+4≤zn+2,zn+2≤yn,yn≤xn−2,⇒xn+4<xn−2,yn+4≤xn+2,xn+2≤zn,zn≤yn−2,⇒yn+4<yn−2,zn+4≤yn+2,yn+2≤xn,xn≤zn−2,⇒zn+4<zn−2.
Therefore, subsequencesx6n−2n=0∞,…,x6n+3n=0∞ are nonincreasing and so are bounded from above by M=maxx−2,x−1,x0,x1,x2,x3. Similarly, the subsequences y6n−2n=0∞,…,y6n+3n=0∞ are nonincreasing and hence bounded above by M=maxy−2,y−1,y0,y1,y2,y3 and z6n−2n=0∞,…,z6n+3n=0∞ are nonincreasing and hence bounded from above by M=maxz−2,z−1,z0,z1,z2,z3.Example 1.
Figure1 represents the dynamics of (10) with x−2=0.18,x−1=0.15,x0=0.13,y−2=0.26, y−1=0.17,y0=−0.13,z−2=0.14,z−1=0.19, and z0=0.15.Figure 1
Difference system (10).


## 3. The System:xn+1=ynzn−1/yn+xn−2,yn+1=znxn−1/zn+yn−2,andzn+1=xnyn−1/xn−zn−2

Here, we will discuss solutions of the following system:(21)xn+1=ynzn−1yn+xn−2,yn+1=znxn−1zn+yn−2,zn+1=xnyn−1xn−zn−2,where n∈ℕ0 and nonzero initial conditions such that x−2≠±y0,x−2≠2y0, y−2≠±z0,y−2≠−2z0, z−2≠x0,z−2≠2x0, and z−2≠3x0.Theorem 2.
Ifxn,yn,zn are solutions of (21), then solutions of (21) are represented by the following formulas for n=0,1,…:(22)x6n−2=y03nx−22n−1x−2−2y0n,x6n−1=x−1z03ny−2+z02ny−2−z0n,x6n=−1nx03n+1z−2nz−2−2x02n,x6n+1=−1nz−1y03n+1x−2+y0n+1x−2−y02n,x6n+2=−1nz03n+2y−22ny−2+2z0n+1,x6n+3=y−1x03n+2z−2−x02n+1z−2−3x0n+1,y6n−2=−1nz03ny−22n−1y−2+2z0n,y6n−1=y−1x03nx0−z−22nz−2−3x0n,y6n=y03n+1x−22nx−2−2y0n,y6n+1=x−1z03n+1y−2+z02n+1y−2−z0n,y6n+2=−1n+1x03n+2z−2nz−2−2x02n+1,y6n+3=−1n+1z−1y03n+2x−2−y02n+1x−2+y0n+1,z6n−2=−1nx03nz−2n−1z−2−2x02n,z6n−1=−1nz−1y03nx−2−y02nx−2+y0n,z6n=−1nz03n+1y−22ny−2+2z0n,z6n+1=−y−1x03n+1z−2−x02n+1z−2−3x0n,z6n+2=−y03n+2x−22n+1x−2−2y0n,z6n+3=−x−1z03n+2y−2−z0ny−2+z02n+2.Proof.
Obviously, results are true forn=0. Now, one suppose that, for n−1, it is also true, i.e.,(23)x6n−8=y03n−3x−22n−3x−2−2y0n−1,x6n−7=x−1z03n−3y−2+z02n−2y−2−z0n−1,x6n−6=−1n−1x03n−2z−2n−1z−2−2x02n−2,x6n−5=−1n−1z−1y03n−2x−2+y0nx−2−y02n−2,x6n−4=−1n−1z03n−1y−22n−2y−2+2z0n,x6n−3=y−1x03n−1z−2−x02n−1z−2−3x0n,y6n−8=−1n−1z03n−3y−22n−3y−2+2z0n−1,y6n−7=y−1x03n−3x0−z−22n−2z−2−3x0n−1,y6n−6=y03n−2x−22n−2x−2−2y0n−1,y6n−5=x−1z03n−2y−2+z02n−1y−2−z0n−1,y6n−4=−1nx03n−1z−2n−1z−2−2x02n−1,y6n−3=−1nz−1y03n−1x−2−y02n−1x−2+y0n,z6n−8=−1n−1x03n−3z−2n−2z−2−2x02n−2,z6n−7=−1n−1z−1y03n−3x−2−y02n−2x−2+y0n−1,z6n−6=−1n−1z03n−2y−22n−2y−2+2z0n−1,z6n−5=−y−1x03n−2z−2−x02n−1z−2−3x0n−1,z6n−4−y03n−1x−22n−1x−2−2y0n−1,z6n−3=−x−1z03n−1y−2−z0n−1y−2+z02n.
Now, from equation (21), it follows that(24)x6n−2=z6n−4y6n−3y6n−3+x6n−5=−y03n−1/x−22n−1x−2−2y0n−1−1nz−1y03n−1/x−2−y02n−1x−2+y0n−1nz−1y03n−1/x−2−y02n−1x−2+y0n+−1n−1z−1y03n−2/x−2+y0nx−2−y02n−2=−y03n−1/x−22n−1x−2−2y0n−1y0y0−x−2−y0=y03nx−22n−1x−2−2y0n,y6n−1=x6n−3z6n−2z6n−2+y6n−4=y−1x03n−1/z−2−x02n−1z−2−3x0n−1nx03n/z−2n−1z−2−2x02n−1nx03n/z−2n−1z−2−2x02n+−1nx03n−1/z−2n−1z−2−2x02n−1=y−1x03n−1/z−2−x02n−1z−2−3x0nx0x0+z−2−2x0=y−1x03nz−2−x02nz−2−3x0nz6n=y6n−2x6n−1x6n−1−z6n−3=−1nz03n/y−22n−1y−2+2z0nx−1z03n/y−2+z02ny−2−z0nx−1z03n/y−2+z02ny−2−z0n−−x−1z03n−1/y−2−z0n−1y−2+z02n%=−1nz03n/y−22n−1y−2+2z0nz0z0+y−2−z0=−1nz03n+1y−22ny−2+2z0n.
Also, from equation (21), one can see that(25)x6n+1=z6n−1y6ny6n+x6n−2=−1nz−1y03n/x−2−y02nx−2+y0ny03n+1/x−22nx−2−2y0ny03n+1/x−22nx−2−2y0n+y03n/x−22n−1x−2−2y0n=−1nz−1y03n/x−2−y02nx−2+y0ny0y0+x−2=−1nz−1y03n+1x−2−y02nx−2+y0n+1,y6n+2=x6nz6n+1z6n+1+y6n−1=−1nx03n+1/z−2nz−2−2x02n−y−1x03n+1/z−2−x02n+1z−2−3x0n−y−1x03n+1/z−2−x02n+1z−2−3x0n+y−1x03n/x0−z−22nz−2−3x0n=−1nx03n+1/z−2nz−2−2x02nx0x0−z−2−x0=−1n+1x03n+2z−2nz−2−2x02n+1,z6n+3=y6n+1x6n+2x6n+2−z6n=x−1z03n+1/y−2+z02n+1y−2−z0n−1nz03n+2/y−22ny−2+2z0n+1−1nz03n+2/y−22ny−2+2z0n+1−−1nz03n+1/y−22ny−2+2z0n=x−1z03n+1/y−2+z02n+1y−2−z0nz0z0−y−2+2z0=−x−1z03n+2y−2+z02n+2y−2−z0n.
In a similar pattern, one can prove the other relation.Example 2.
Figure2 represents the dynamics of (21) with x−2=0.8,x−1=0.5,x0=0.3, y−2=0.26y−1=0.17,y0=−0.13,z−2=0.14,z−1=0.19, and z0=0.15.
In a similar pattern, one can prove the following cases.Figure 2
Difference system (21).


## 4. The System:xn+1=ynzn−1/yn−xn−2, yn+1=znxn−1/zn−yn−2, and zn+1=xnyn−1/xn−zn−2

We discuss the solutions of the following system:(26)xn+1=ynzn−1yn−xn−2,yn+1=znxn−1zn−yn−2,zn+1=xnyn−1xn−zn−2,where n∈ℕ0 and nonzero initial conditions such that x−2≠y0,y−2≠z0, and x0≠z−2.Theorem 3.
Ifxn,yn,zn are solutions of (26), then(27)x6n−2=x−2y0x−23n,x6n−1=x−1z0z0−y−23n,x6n=x0x0z−23n,x6n+1=−1n+1z−1y0x−2−y03n+1,x6n+2=z0z0y−23n+1,x6n+3=y−1x0x0−z−23n+2,y6n−2=y−2z0y−23n,y6n−1=y−1x0x0−z−23n,y6n=y0y0x−23n,y6n+1=−1n+1x−1z0y−2−z03n+1,y6n+2=x0x0z−23n+1,y6n+3=−1nz−1y0x−2−y03n+2,z6n−2=z−2x0z−23n,z6n−1=−1nz−1y0x−2−y03n,z6n=z0z0y−23n,z6n+1=y−1x0x0−z−23n+1,z6n+2=x−2y0x−23n+2,z6n+3=−1nx−1z0y−2−z03n+2.Example 3.
Figure3 represents the behavior of (26) if x−2=2,x−1=0.81,x0=0.49,y−2=0.3, y−1=−0.09,y0=−0.4,z−2=0.32,z−1=0.12,andz0=0.2.Figure 3
Difference system (26).


## 5. The System:xn+1=ynzn−1/yn+xn−2, yn+1=znxn−1/zn−yn−2, and zn+1=xnyn−1/xn+zn−2

We discuss the solutions of the following system:(28)xn+1=ynzn−1yn+xn−2,yn+1=znxn−1zn−yn−2,zn+1=xnyn−1xn+zn−2,where n∈ℕ0 and nonnegative initial conditions such that x−2≠±y0,x−2≠−2y0, y−2≠z0,y−2≠2z0, y−2≠3z0,z−2≠±x0, and z−2≠2x0.Theorem 4.
Ifxn,yn,zn are solutions of (28), then(29)x6n−2=−1ny03nx−22n−1x−2+2y0n,x6n−1=x−1z03ny−2−z02ny−2−3z0n,x6n=−1nx03n+1z−22n2x0−z−2n,x6n+1=z−1y03n+1x−2−y0nx−2+y02n+1,x6n+2=−1n+1z03n+2y−2ny−2−2z02n+1,x6n+3=−1ny−1x03n+2x0−z−22n+1x0+z−2n+1,y6n−2=−1nz03ny−2n−1y−2−2z02n,y6n−1=−1ny−1x03nx0−z−22nx0+z−2n,y6n=−1ny03n+1x−22nx−2+2y0n,y6n+1=−x−1z03n+1y−2−z02n+1y−2−3z0n,y6n+2=−1n+1x03n+2z−22n+12x0−z−2n,y6n+3=−z−1y03n+2x−2+y02n+2x−2−y0n,z6n−2=−1nx03nz−22n−12x0−z−2n,z6n−1=z−1y03nx−2−y0nx−2+y02n,z6n=−1nz03n+1y−2ny−2−2z02n,z6n+1=−1ny−1x03n+1x0−z−22nx0+z−2n+1,z6n+2=−1ny03n+2x−22nx−2+2y0n+1,z6n+3=x−1z03n+2y−2−3z0n+1y−2−z02n+1.Example 4.
Figure4 represents the dynamics of (28) with x−2=0.8,x−1=0.25,x0=−0.31,y−2=0.2, y−1=0.7,y0=−0.13,z−2=0.4,z−1=0.9, and z0=0.71.Figure 4
Difference system (28).


## 6. The System:xn+1=ynzn−1/yn−xn−2, yn+1=znxn−1/zn+yn−2, and zn+1=xnyn−1/xn+zn−2

We will obtain the solutions of the following system:(30)xn+1=ynzn−1yn−xn−2,yn+1=znxn−1zn+yn−2,zn+1=xnyn−1xn+zn−2,where n∈ℕ0 and nonzero initial conditions such that x−2≠y0,x−2≠2y0, x−2≠3y0, y−2≠±z0,y−2≠2z0, z−2≠±x0,and z−2≠−2x0.Theorem 5.
Ifxn,yn,zn are solutions of (30), then(31)x6n−2=−1ny03nx−2n−1x−2−2y02n,x6n−1=−1nx−1z03ny−2+z0ny−2−z02n,x6n=−1nx03n+1z−22nz−2+2x0n,x6n+1=−z−1y03n+1x−2−3y0nx−2−y02n+1,x6n+2=−z03n+2y−22n+1y−2−2z0n,x6n+3=−1n+1y−1x03n+2x0−z−2nx0+z−22n+2,y6n−2=z03ny−22n−1y−2−2z0n,y6n−1=−1ny−1x03nx0−z−2nx0+z−22n,y6n=−1ny03n+1x−2nx−2−2y02n,y6n+1=−1nx−1z03n+1y−2−z02ny−2+z0n+1,y6n+2=−1nx03n+2z−22nz−2+2x0n+1,y6n+3=z−1y03n+2x−2−3y0n+1x−2−y02n+1,z6n−2=−1nx03nz−22n−1z−2+2x0n,z6n−1=z−1y03nx−2−y02nx−2−3y0n,z6n=z03n+1y−22ny−2−2z0n,z6n+1=−1ny−1x03n+1x0−z−2nx0+z−22n+1,z6n+2=−1n+1y03n+2x−2nx−2−2y02n+1,z6n+3=−1n+1x−1z03n+2y−2−z02n+1y−2+z0n+1.Example 5.
Figure5 represents the dynamics of (30) with x−2=−0.2,x−1=0.05,x0=−0.39,y−2=2, y−1=1.7,y0=−0.3,z−2=1.1,z−1=8,andz0=−0.16.Figure 5
Difference system (30).


## 7. Conclusion

This paper discussed boundedness and the form expressions of the solutions of some special cases of difference equation systemsxn+1=ynzn−1/yn±xn−2,yn+1=znxn−1/zn±yn−2,andzn+1=xnyn−1/xn±zn−2, after the introduction in Section 1, we investigated in Section 2, the boundedness behavior of equation (10), and we have obtained the solutions of this system. In Section 3, we obtained the form of the solutions of the system xn+1=ynzn−1/yn+xn−2,yn+1=znxn−1/zn+yn−2,andzn+1=xnyn−1/xn−zn−2. Also, in Section 4, we gave the solution form of the system xn+1=ynzn−1/yn−xn−2,yn+1=znxn−1/zn−yn−2,andzn+1=xnyn−1/xn−zn−2. In Section 5, we dealt with the solution of system (28). Finally, in Section 6, we have the expressions of the form of the solutions of the system xn+1=ynzn−1/yn−xn−2,yn+1=znxn−1/zn+yn−2,andzn+1=xnyn−1/xn+zn−2, and in each case, we presented some numerical examples of each case and drew it by using MATLAB.



---

*Source: 2480294-2021-04-16.xml*