Dynamics of ratio-dependent Predator-Prey models with nonconstant harvesting

The dynamics of constant harvesting of a single species has been studied extensively within the framework of ratio-dependent predator-prey models. In this work, we investigate the properties of a Michaelis-Menten ratio-dependent predator-prey model with two nonconstant harvesting functions depending on the prey population. Equilibria and periodic orbits are computed and their stability properties are analyzed. Several bifurcations are detected as well as connecting orbits, with an emphasis on analyzing the equilibrium points at which the species coexist. Smooth numerical continuation is performed that allows computation of branches of solutions.