Degenerate Bogdanov-Takens bifurcation of codimension 4 in Holling-Tanner model with harvesting

In this paper, we revisit the Holling-Tanner model with constant-yield prey harvesting. It is shown that the highest codimension of a nilpotent cusp is 4, and the model can undergo degenerate Bogdanov-Takens bifurcation of codimension 4. Moreover, when the model has a center-type equilibrium, we show that it is a weak focus with order at least 3 and at most 4, and the model can exhibit Hopf bifurcation of codimension 3. Some algebraic methods including resultant elimination and pseudo-division are used to solve the semialgebraic varieties of normal form coefficients or focal values. Our results indicate that constant-yield prey harvesting can cause not only richer dynamics and bifurcations, but also the coextinction of both populations with some positive initial densities. Finally, numerical simulations, including the coexistence of limit cycle and homoclinic cycle, and three limit cycles, are presented to illustrate the theoretical results. (c) 2022 Elsevier Inc. All rights reserved.