Codimension-2 bifurcation in a discrete predator–prey system with constant yield predator harvesting

This work investigates the bifurcation analysis in a discrete-time Leslie–Gower predator–prey model with constant yield predator harvesting. The stability analysis for the fixed points of the discretized model is shown briefly. In this study, the model undergoes codimension-1 bifurcation such as fold bifurcation (limit point), flip bifurcation (period-doubling) and Neimark–Sacker bifurcation at a positive fixed point. Further, the model exhibits codimension-2 bifurcations, including Bogdanov–Takens bifurcation and generalized flip bifurcation at the fixed point. For each bifurcation, by using the critical normal form coefficient method, various critical states are calculated. To validate our analytical findings, the bifurcation curves of fixed points are drawn by using MATCONTM. The system exhibits interesting rich dynamics including limit cycles and chaos. Moreover, it has been shown that the predator harvesting may control the chaos in the system.