Complex dynamics of a prey-predator interaction model with Holling type-II functional response incorporating the effect of fear on prey and non-linear predator harvesting

In the present article, we have investigated the impact of the fear effect and non-linear predator harvesting in a prey-predator interaction model with Holling type-II functional response. Fear for predator population enhances the survival probability of prey population, and it can greatly reduce the reproduction of prey population. We have investigated all the biologically feasible equilibrium points and the positivity and boundedness of the system solutions. We have analyzed the local and global stability of the feasible equilibrium points in terms of the model parameters. Analytically we have established that the intrinsic growth rate of the prey population can stabilize the prey-predator interactions by ignoring the existence of periodic behaviors. The model system undergoes through Transcritical, Saddle-Node, Hopf bifurcation by considering the intrinsic growth rate of the prey population as a bifurcation parameter and Bogdanov–Takens bifurcation with respect to the prey intrinsic growth rate and rate of predator harvesting. Numerically we have identified two parametric planes, which are divided into sub-regions associated with different numbers and nature of the equilibrium points by various bifurcation curves. We also found that the system may exhibit bi-stability behavior by producing stable axial and interior equilibrium points. Lastly, the manuscript is concluded with some recommendations.