Abstract
A proposed discretized form of fractional-order prey-predator model is investigated. A sufficient condition for the solution of the discrete system to exist and to be unique is determined. Jury stability test is applied for studying stability of equilibrium points of the discretized system. Then, the effects of varying fractional order and other parameters of the systems on its dynamics are examined. The system undergoes Neimark-Sacker and flip bifurcation under certain conditions. We observe that the model exhibits chaotic dynamics following stable states as the memory parameter alpha decreases and step size h increases. Theoretical results illustrate the rich dynamics and complexity of the model. Numerical simulation validates theoretical results and demonstrates the presence of rich dynamical behaviors include S-asymptotically bounded periodic orbits, quasi-periodicity, and chaos. The system exhibits a wide range of dynamical behaviors for fractional-order alpha key parameter.