Identifying numerical bifurcation structures of codimensions 1 and 2 in interacting species system

Abstract

The present study has focused on the examination of a Holling type III nonlinear mathematical model that incorporates the influence of fear and Michaelis–Menten‐type predator harvesting. The incorporation of fear with respect to the prey species has been observed to result in a reduction in the survival probability of the prey population and a concurrent reduction in the reproduction rate of the prey species. The existence and stability of ecologically significant equilibria have been ascertained through mathematical analysis. Emphasis within the proposed model primarily centers on numerical bifurcations of codimensions 1 and 2. Numerical validation has been performed on all simulated outcomes within the feasible range of parametric values. Dynamical characteristics of the model have subsequently undergone investigation through a series of numerical simulations, successfully revealing various forms of local and global bifurcations. In addition to the identification of saddle‐node, Hopf, Bogdanov–Takens, transcritical, cusp, homoclinic, and limit point cycle (LPC) bifurcations, the model has also demonstrated bistability and global asymptotic stability. These bifurcation phenomena serve as illustrative examples of the intricate dynamical behavior inherent to the model. Numerical validation through graphical representations has been utilized to elucidate the effects of factors such as fear, nonlinear predator harvesting, and predation rate on the dynamics of the interacting species under different parametric conditions.