On Some Solvable Systems of Some Rational Difference Equations of Third Order

: Our aim in this paper is to obtain formulas for solutions of rational difference equations such as x n + 1 = 1 ± ( x n − 1 y n ) / ( 1 − y n ) , y n + 1 = 1 ± ( y n − 1 x n ) / ( 1 − x n ) , and x n + 1 = 1 ± ( x n − 1 y n − 2 ) / ( 1 − y n ) , y n + 1 = 1 ± ( y n − 1 x n − 2 ) / ( 1 − x n ) , where the initial conditions x − 2 , x − 1 , x 0 , y − 2 , y − 1 , y 0 are non-zero real numbers. In addition, we show that the some of these systems are periodic with different periods. We also verify our theoretical outcomes at the end with some numerical applications and draw it by using some mathematical programs to illustrate the results.


Introduction
Nonlinear difference equations have recently captured the attention of numerous scholars.In fact, during the past ten years, we have encouraged a rapid rise in interest in these kinds of equations.The fact that these types of equations have several applications not just in mathematics but also in related areas, particularly in biological sciences, engineering, ecology, discrete temporal systems, economics, physics, and so on, may have contributed to the desire.We believe that, as more appealing and engaging results are achieved and communicated in studies, this area of research will continue to captivate the minds of more scholars in the years to come.The challenge of solving nonlinear difference equations in closed form has emerged as a common theme in this research area.In reality, a large number of articles attempt to solve nonlinear difference equations in any way they can; for an example, see [1][2][3][4][5][6].Evidently, it can be very difficult to obtain the solution form for these kinds of equations in general.However, a number of approaches have lately been proposed to simplify challenging nonlinear difference equations into linear forms with established solution forms.For instance, a sizable class of nonlinear difference equations were solved in closed-form by converting into linear types (see, e.g., [7][8][9][10][11][12]).
Numerous academics have examined how systems of solved difference equations behave, for instance: Cinar examined the answers to the following system of difference equations in [13].
El-Metwally [14] found the solutions form for the following systems of rational difference equations: Kara and Yazlik in [15] showed that the following three-dimensional system of difference equations , can be solved.Furthermore, they determined the forbidden set of the initial conditions by employing acquired formulas.Finally, they provided various applications involving the difference equation system discussed before.Mansour et al. [16] examined the behavior of the difference equations systems' solutions In [17], Ozban studied the positive solutions of the following system of rational difference equations Sroysang [18] focused on a system of a rational higher-order difference equation Touafek et al. [19] investigated the periodic nature and provided the form of the solutions of the following systems of rational difference equations Furthermore, Yalçınkaya [20] has obtained the sufficient conditions for the global asymptotic stability of the following system of two nonlinear difference equations In [21], Zhang et al. studied the boundedness, the persistence, and global asymptotic stability of the positive solutions of the following system Zhang et al. [22] studied the dynamics of a system of the rational third-order difference equation For more studies for nonlinear difference equations and systems of rational difference equations see [2,[23][24][25][26][27][28][29][30][31][32].
Furthermore, difference equations are appropriate models for describing situations where population growth is not continuous but seasonal with overlapping generations.
Researchers have looked at the generalized Beverton-Holt stock recruitment model in [33] Khaliq et al. [34] studied the dynamical analysis of the following system of discretetime two-predators and the one-prey Lotka-Volterra model The boundedness character, persistence, local, and global behavior of the following two-directional interacting and invasive species model were examined by Din and Elsayed [35] The authors in [36] explored local dynamics with topological classifications, bifurcation analysis, and chaos control in a discrete-time COVID-19 epidemic model.See also [37][38][39].
In this paper, we deal with the existence of the form for the solutions of the following systems of difference equations and with the initial conditions x −2 , x −1 , x 0 , y −2 , y −1 , and y 0 are arbitrary non zero real numbers.

Main Results
2.1.System (1) When δ = +1 and γ = +1 In this section, we investigate the solutions of the following system of two difference equations where n ∈ N 0 and the initial conditions are arbitrary non zero real numbers with x 0 = 1 and y 0 = 1.

2.
{x n , y n } ∞ n=−2 has the following form Proof.From Equation (1), we have Furthermore, we see from Equation (1) that Finally we obtain This completes the proof.

System (1) When δ = −1 and γ = +1
The system of difference equations' solutions are provided in this section since n ∈ N 0 and the initial values are arbitrary nonzero real numbers such that x 0 = 1 and y 0 = 1.
Theorem 2. Assume {x n , y n } is a solution of System (2).Then, for n = 0, 1, 2, . . ., Proof.The result is true for n = 0. Assume n exceeds 1 and that n − 1 is consistent with our premise, which is , x −1 , It follows from System (2) that . Furthermore, we obtain from System (2) that x −1 .
Similarly, we can prove the other relations.This completes the proof.
Lemma 1.Let {x n , y n } be a solution of System (2), then {x n }, {y n } are unbounded solutions.
Proof.The proof follows from the expressions of solutions of System (2).
The following theorems can be proved similar to the previous theorem so it will be omitted.Theorem 3. Assume that the sequences {x n , y n } are a solution of the system Then, all solutions are unbounded and provided by the following formulas where n = 0, 1, 2, . . .and x 0 , y 0 = 1.

System (2) When δ = +1 and γ = +1
In this section, we obtain the form of the solutions of the following system of difference equations where n ∈ N 0 and the initial conditions are arbitrary non-zero real numbers such that x 0 , y 0 = 1.
Proof.The conclusion is true for n = 0. Assume n exceeds 1 and that n − 1 is covered by our supposition, so that Now, it follows from System (3) that We can prove the other relations similarly.The proof is complete.
Lemma 2. Let {x n , y n } be a solution of System (3), then {x n }, {y n } are unbounded solutions.
The following theorems can be treated similarly to the previous results.Theorem 6.The solutions of the system where x 0 , y 0 = 1 are provided as follows Lemma 3. If x 0 + x −2 = 1, and y 0 + y −2 = 1, then the solutions are bounded and periodic with period four and adopt the form Otherwise, every solution of System (4) is unbounded.
Theorem 7. Every solution {x n , y n } of the following system with non-zero real numbers, the initial conditions satisfies x 0 , y 0 = 1 adopts the form Lemma 4. If x 0 + x −2 = 1 and y 0 + y −2 = 1, then the solutions of System ( 5) are bounded and periodic with period eight as follows Otherwise, every solution of System ( 5) is unbounded.
Remark 1.As in the previous theorem, we can obtain the solutions of the following system

Numerical Examples
In this section, we present some numerical examples that support the above theoretical results.Example 4. Suppose the initial conditions for the system x n+1 = 1 − x , and y 0 = 0.7.See

Conclusions
In this paper, we obtained the expressions of the solutions of different classes of third-order rational systems of difference equations.In Section 1, the work of the authors on the same side of the difference equations, whether they are equations, systems of equations, or some applications of difference equations, is presented.After the introduction in Section 2, we have solved the first system of second order rational difference equations x n+1 = 1 + x n−1 y n 1 − y n , y n+1 = 1 + y n−1 x n 1 − x n .After finding the solutions, we provided numerical examples to illustrate the results.In Section 3, we obtained the form of the solution of the second system x n+1 = 1 − x n−1 y n 1 − y n , y n+1 = 1 + y n−1 x n 1 − x n ; we also mentioned the solutions of the other systems x n+1 = 1 + x n−1 y n 1 − y n , y n+1 = 1 − y n−1 x n 1 − x n , and x n+1 = 1 − x n−1 y n 1 − y n , y n+1 = 1 − y n−1 x n 1 − x n .Finally, Section 4 was devoted to the form of the solutions of the main system of third order fractional difference equations x n+1 = 1 + x n−1 y n−2 1 − y n , y n+1 = 1 + y n−1 x n−2 1 − x n and some other systems that we obtained the expressions of the solutions for and studied the periodicity nature of the solutions.Moreover, we confirmed our results using numerical simulations and drew them using Matlab program.