Global dynamics and numerical bifurcation of a bioeconomic system

A predator-prey model was extended to include nonlinear harvesting of the predator guided by its population, such that harvesting is only implemented if the predator population exceeds an economic threshold. Theoretical results showed that the harvesting system undergoes multiple bifurcations, including fold, supercritical Hopf, Bogdanov-Takens and cusp bifurcations. We determine stability and dynamical behaviors of the equilibrium of this system. Numerical simulation results are given to support our theoretical results.