Bifurcation and stability analysis of a ratio-dependent predator-prey model with predator harvesting rate

In this paper, we study the bifurcation and stability of a ratio-dependent predator-prey model with nonconstant predator harvesting rate. The analysis is carried out both analytically and numerically. We determine stability and dynamical behaviours of the equilibria of this system and characterize codimension 1 and codimension 2 bifurcations of the system analytically. Our bifurcation analysis indicates that the system exhibits numerous types of bifurcation phenomena, including Fold, Hopf, Cusp, and Bogdanov–Takens bifurcations. We use the numerical software MATCONT, to compute curves of equilibria and to compute several bifurcation curves. We especially approximate a family of limit cycles emanating from a Hopf point. Our results generalize and improve some known results and show that the model has more rich dynamics than the ratio-dependent predator-prey model without harvesting rate.