Global behavior of third order system of rational difference equations

In this paper, our aim is to study the dynamical behavior of third-order system of rational difference equations xn+1 = αxn−2 β + γxnxn−1xn−2 , yn+1 = α1yn−2 β1 + γ1ynyn−1yn−2 , n = 0, 1, · · · . where the parameters α, β, γ, α1, β1, γ1 and initial conditions x0, x−1, x−2, y0, y−1, y−2 are positive real numbers. Some numerical examples are given to verify our theoretical results.