In this paper we study the complex dynamics of predator-prey systems with nonmonotonicfunctional response and harvesting. When the harvesting is constant-yield for prey, it isshown that various kinds of bifurcations, such as saddle-node bifurcation, degenerate Hopfbifurcation, and Bogdanov-Takens bifurcation, occur in the model as parameters vary. Theexistence of two limit cycles and a homoclinic loop is established by numericalsimulations. When the harvesting is seasonal for both species, sufficient conditions forthe existence of an asymptotically stable periodic solution and bifurcation of a stableperiodic orbit into a stable invariant torus of the model are given. Numerical simulationsare carried out to demonstrate the existence of bifurcation of a stable periodic orbitinto an invariant torus and transition from invariant tori to periodic solutions,respectively, as the amplitude of seasonal harvesting increases.
