The bifurcation analysis of a predator–prey system of Holling and Leslie type with constant-yield prey harvesting is carried out in this paper. It is shown that the model has a Bogdanov–Takens singularity (cusp case) of codimension at least 4 for some parameter values. Various kinds of bifurcations, such as saddle-node bifurcation, Hopf bifurcation, repelling and attracting Bogdanov–Takens bifurcations of codimensions 2 and 3, are also shown in the model as parameters vary. Hence, there are different parameter values for which the model has ...