Global dynamics of a 3 × 6 exponential system of difference equations

We study the global dynamics of following 3×6 exponential system of difference equations: $$x_{n+1}=\frac{{\alpha }_1+{\beta }_1{e}^{?x_n}}{{\gamma }_1+y_{n?1}},y_{n+1}=\frac{{\alpha }_2+{\beta }_2{e}^{?y_n}}{{\gamma }_2+z_{n?1}},z_{n+1}=\frac{{\alpha }_3+{\beta }_3{e}^{?z_n}}{{\gamma }_3+x_{n?1}},n=0,1,\text{?},$$ where parameters α_i,β_i,γ_i (i=1,2,3) and initial conditions x_i,y_i,z_i (i=0,?1) are nonnegative real numbers. The proposed work is considerably extended and improve some existing results in the literature. Finally, theoretical results are verified by numerical simulations.