Solution Form of a Higher Order System of Difference Equation and Dynamical Behavior of Its Special Case

Abstract

The solution form of the system of nonlinear difference equations \begin{equation*} x_{n+1} = \frac{x_{n-k+1}^{p}y_{n}}{a y_{n-k}^{p}+b y_{n}},\ y_{n+1} = \frac{y_{n-k+1}^{p}x_{n}}{αx_{n-k}^{p}+βx_{n}}, \quad n, p \in \mathbb{N}_{0},\ k\in \mathbb{N}, \end{equation*} where the coefficients $a,b,\alpha ,\beta$ and the initial values $x_{-i},y_{-i},i\in \{0,1,\dots ,k\}$ are real numbers, is obtained. Furthermore, the behavior of solutions of the above system when $p=1$ is examined. Numerical examples are presented to illustrate the results exhibited in the paper.

A preprint of this manuscript, which is of 13 pages with 8 figures, is submitted for publication