A Simple Predator-Prey Population Model with Rich Dynamics

A non-smooth switched harvest on predators is introduced into a simple predator-prey model with logistical growth of the prey and a bilinear functional response. If the density of the predator is below a switched value, the harvesting rate is linear; otherwise, it is constant. The model links the well studied predator-prey model with constant harvesting to that with a proportional harvesting rate. It is shown that when the net reproductive number for the predator is greater than unity, the system is permanent and there may exist multiple positive equilibria due to the effects of the switched harvest, a saddle-node bifurcation, a limit cycle, and the coexistence of a stable equilibrium and a unstable circled inside limit cycle and a stable circled outside limit cycle. When the net reproductive number is less than unity, a backward bifurcation from a positive equilibrium occurs, which implies that the stable predator-extinct equilibrium may coexist with two coexistence equilibria. In this situation, reducing the net reproductive number to less than unity is not enough to enable the predator to go extinct. Numerical simulations are provided to illustrate the theoretical results. It seems that the model possesses new complex dynamics compared to the existing harvesting models.


Introduction
Mathematical modeling of predator-prey interactions have attracted wide attention since the original work by Lotka and Volterra in 1920s, and there have been extensively studied for their rich dynamics [1][2][3].Since the rich and complex dynamics for interactive species are common in the real world, many researchers have investigated the processes that affect the dynamics of prey-predator models and wanted to know what models can best represent species interactions.
As a simplest form, the interaction between a predator and prey may be modeled by a pair of differential equations [1,[3][4][5], where N and P represent the prey and predator species, respectively; r, K, a, c, and d are positive constants.In the absence of the predation, the prey grows logistically with intrinsic growth rate r and carrying capacity K.In the presence of the predator, the prey species decreases at a rate proportional to the functional response aN, where a presents the rate of predation.The factor c denotes the efficiency of predation which divides a maximum per capita birth rate of the predators into a maximum per capita consumption rate.Without the prey, no predation occurs and the predator species decreases exponentially with mortality rate d.
To consider the dynamics of model Equation (1), it is shown that the origin (0, 0) is a saddle point.Define the net reproductive number of the predator population n 0 , i.e., the expected number of a predator individual producing as the predator population is introduced into a stable prey population [6,7], as Then if n 0 < 1, the boundary equilibrium (K, 0) with the predator going extinct is globally asymptotically stable and there exists no positive equilibrium for the prey-predator interaction.If n 0 > 1, the boundary equilibrium (K, 0) is a saddle point and there exists a positive coexistence equilibrium (N * , P * ) = K n 0 , r a 1 − 1 n 0 , which is globally asymptotically stable.The model dynamics are relatively simple.Using n 0 as bifurcation parameter, we have a transcritical bifurcation at n 0 = 1 as (N * , P * ) is bifurcated.Notice that model Equation ( 1) is a special case of the prey-predator model in [8].
In the real world, from the point of view of predators' needs, the exploitation of biological resources and harvest are commonly practiced in fishery, forestry, and wildlife management.There is an interest in the use of bioeconomic models to gain insight into the scientific management of renewable resources [16].Moreover, harvesting is an important and effective method to prevent and control the explosive growth of predators or prey when they are enough.So, generally speaking, it is reasonable and necessary for one to introduce the harvest of populations into models.Taking the above reasons into a consideration, we focus on the predator-prey model with harvest [12,13,16,17,[20][21][22][23][24][25][26][27][28].
Normally, harvesting has several forms in predator-prey models.The most common one of these harvesting forms is a nonzero constant [16,17,20,[22][23][24]26] or a linear harvesting rate [12,13,22,[25][26][27][28].In Ref. [22], a two-prey-one-predator model with predator harvested was studied.The authors are particularly interested in the stability properties of different harvest strategies.Two types of harvest strategies are: with a nonzero constant and a linear harvesting rate.The choice of idealized harvest strategies will contribute to a qualitative understanding of the properties of different harvesting strategies.Xiao and Jennings [16] considered the dynamical properties of the ratio-dependent predator-prey model with constant prey harvesting.There existed numerous kinds of bifurcations, such as the saddle-node bifurcation, the subcritical and supercritical Hopf bifurcations.There also existed a limit cycle, a homoclinic or heteroclinic orbit satisfying different parameter values.In Ref. [17], the ratio-dependent predator-prey model with constant predator harvesting was focused on.Philip et al. [26] also discussed two predator-prey models with linear or nonzero constant predator harvesting.
The above two types of harvesting rates seemingly have their own advantages as well as disadvantages in fitting the harvest in the real world.When the density of the predator or prey is rather low, the nonzero constant harvesting rate is not as reasonable as of the proportional type [13,22,26]; while if the predator or prey is abundant, linear harvesting rate is less possible than the constant harvesting rate [16].In Ref. [22], to compare the stability properties of the system with two harvest strategies, they applied linear or constant harvesting rate, respectively.In that comparative study, the authors demonstrated that switching from linear to constant harvesting rate may turn a stable stationary state to a periodic or chaotic oscillatory mode from a mathematical perspective.However, when deciding the constant level of harvesting, the instability of the constant harvest strategy calls for great care.In Ref. [29], Beddington et al. introduced a more realistic smooth harvesting function in which the fishing effort is limited upwards because the constant harvest cannot be achieved for small populations of the fish.The adjustment of a harvesting function can prevent extinction and increase the stability to some extent.Moreover, the dramatic increase of the predator or prey challenges the normal ecological balance and capacity of harvest.Thus, it is interesting to construct a new kind of harvesting rate that combines the advantages from both linear and constant harvesting rates.
Motivated by these ideas, in this paper, we consider a predator-prey model with a novel harvesting rate.Our ideas to develop the harvesting rate are derived from the capacity of treatments of diseases that had been well studied in the dynamical epidemic models [30][31][32].
Using model Equation ( 1) as our baseline model, we assume that harvesting takes place, but only the predator population is under harvesting and introduce harvesting function H(P) of the predator to prey-predator model Equation (1) for discussing its dynamical features.The interactive dynamics are governed by the following system Following the methods in [14][15][16][17]30], we investigate the existence and stability of multiple equilibria, bifurcations, and limit cycles, and study the effects of switched harvest on the dynamics of the predator-prey model.
This paper is organized as follow.Sections 2 and 3.1 represent the boundedness of model Equation (3) and existence of multiple equilibria.In Section 3.2, we study the stability of equilibria, bifurcations, and the existence and stability of a limit cycle.In Section 4, we give numerical simulations to verify our results.Brief discussions are presented finally in Section 5.

Model Formulation
Now, we consider model Equation (3).Firstly, we describe harvesting function H(P) of the predators in model Equation (3), which has the following form We assume that the harvesting rate is proportional to the predator population size until it reaches a threshold value due to limited facilities of harvesting or resource protection.The harvesting rate will then be kept as a constant.Denote the harvesting threshold value as h = mP 0 .
When 0 ≤ P ≤ P 0 , model Equation ( 3) is When P > P 0 , model Equation (3) becomes It is straightforward to verify that solutions of Equation (3) with positive initial conditions are all positive for t ≥ 0 and ultimately bounded.Thus the following set is positive invariant for system Equation (3).

Existence of Equilibria
In this section, we explore the existence of all nonnegative equilibria.First, the origin (0, 0) is still a trivial equilibrium and the predator-free equilibrium (K, 0) exists.Moreover, it is easy to see that there exists no positive equilibrium in region D if n 0 ≤ 1.We thus assume n 0 > 1 hereafter, and present our results of the existence of positive equilibria as follows.
Theorem 1. System Equation (3) has a positive coexistence equilibrium E * (N * , P * ), in the subregion of D with 0 < P ≤ P 0 , if and only if P 0 ≥ r a and n > 1, or P 0 < r a and 1 < n ≤ r r − aP 0 , where r a represents the maximum predator density for which the prey population can establish itself from a small initial population and n is the net reproductive number of the predator under harvesting defined by Proof.In the subregion of D with 0 < P ≤ P 0 , a positive equilibrium of Equation (3) satisfies Then it follows that and then In the mean time, it follows from that there exists a positive equilibrium in this subregion if To investigate the existence of positive equilibria in the subregion of D with P > P 0 , we first give the relation of roots and coefficients for a quadratic equation.
Quadratic equation (x − A)(B − x) = C, with constants A, B, and C positive, has two positive roots respectively.
The results for the existence of positive equilibria of Equation ( 3) in the subregion of D with P > P 0 are provided as follows.
Theorem 2. We assume P 0 < r a and define System Equation (3), in the subregion of D with P > P 0 , has where and In the subregion of D with P > P 0 , a positive equilibrium satisfies It follows from Equation (14a) that Substituting it into Equation (14b) yields It follows from the relation of roots and coefficients for the quadratic equation To have P i > P 0 , i = 1, 2, we need that is, Thus, if P 0 ≥ r/a, there is no positive equilibrium of Equation (15) in the subregion of D with P > P 0 .We assume P 0 < r/a.
Suppose there are two positive solutions, N 1 < N 2 , to Equation ( 15) and hence P 2 < P 1 .Then P 2 > P 0 if and only if Substituting N 2 in Equation ( 13) into Equation ( 17) yields Then if P 0 ≥ P, we have P 2 ≤ P 0 .Assume P 0 < P. Squaring both sides of Equation ( 18) yields It follows from the definition of h 1 that If n < 1, we have P 2 > P 0 .Otherwise, Equation ( 20) is equivalent to Thus, if h 2 < h < h 1 , P 2 > P 0 , and if h ≤ h 2 , P 2 < P 0 .We now consider P 1 > P 0 , that is, Substituting N 1 in Equation ( 13) into Equation ( 21), we have If P 0 ≤ P, inequality Equation ( 21) is satisfied which implies P 1 > P 0 .Suppose P 0 > P. Similarly as above, inequality Equation ( 22) is equivalent to Thus, if h < h 2 , P 1 > P 0 , and if h ≥ h 2 , P 1 < P 0 .By putting all together, the proof is completed.
Remark 1. System Equation (1) with linear predator harvest strategy leads to the predator extinct if the net reproductive number n < 1.However, by Theorem 2 (c), we find that for system Equation (1) with switched predator harvest strategy; that is, when the density of predator is below harvest level P 0 , the linear harvesting rate is applied to the system, whereas when the density of predators is higher than harvest level, the system adopts nonzero constant harvesting rate, even if n < 1, the prey and predator may coexist.

Stability of Equilibria
In this section, we discuss the stability of equilibria of model Equation (3).
Proof.It is easy to obtain that the characteristic roots to the linearized equation of system Equation ( 3) Thus, E 0 is locally asymptotically stable if n < 1 and unstable if n > 1.
Next, note that (0, 0) is always unstable.If n < 1, by Theorems 1 and 2, there is not any other equilibrium of system Equation ( 3) than E 0 in D\{(0, 0)}.Since D is the invariant set of system Equation ( 3) and E 0 is locally asymptotically stable, it follows from the Bendixson Theorem that every solution of system Equation (3) in D approaches E 0 when t tends to positive infinity.
Since E 0 is unstable as n > 1, following the similar arguments to Cantrell and Cosner ([33] Theorem 3.1) (see also ([34] Theorems 3,4)), which is based on the uniform persistence theory introduced by Hale and Waltman [35], we are able to conclude that system Equation ( 3) is permanent if n > 1.
The proof is complete.
Remark 2. It follows from Cantrell and Cosner ([33] Theorem 3.1) that system Equation ( 5) is permanent if and only if n > 1.Thus Theorem 3 suggests that if system Equation ( 5) is permanent, then so is system Equation (3).
For the three positive equilibria E * , E 1 , E 2 , we have the following results.Firstly, we consider the stability of E * .Theorem 4. The positive equilibrium E * (N * , P * ) of system Equation ( 3) is globally asymptotically stable if P 0 < r a and 1 < n ≤ r r − aP 0 .
Proof.According to Theorem 1, positive equilibrium E * (N * , P * ) exists if and only if P 0 < r a and . The Jacobian matrix of system Equation ( 3) at E * (N * , P * ) is Because N * and P * satisfy Equation (9), by means of Equation ( 9), the trace and determinant of J * are simplified as Therefore, all eigenvalues of matrix J * have negative real parts when P 0 < r a and 1 < n ≤ r r − aP 0 .It follows that E * (N * , P * ) is locally asymptotically stable.
Then, it suffices for us to prove the global attractiveness of E * (N * , P * ).Inspired by the work of McCluskey [36], we define a Lyapunov function where N * = K n , and N * − 1 and P(t) P * − ln P(t) P * − 1 ≥ 0 for all N(t), P(t) > 0. (The equality holds if and only if N(t) = N * , P(t) = P * .)From the definition of M(t), we know that M(t) is well-defined and M(t) ≥ 0. The equality holds if and only if N(t) = N * and P(t) = P * .
Differentiating M(t) along the solutions of system Equation (3), we obtain It follows that M(t) is bounded and non-increasing.Thus lim t→∞ M(t) exists.Note that dM dt = 0 if and only if N = N * .Substituting N = N * into the first equation of Equation ( 3), one can directly get P = P * .Therefore, the maximal compact invariant set in dM dt = 0 is the singleton E * .By the LaSalle invariance principle (see, for example, Theorem 5.3.1 in Hale and Verduyn Lunel [37]), positive equilibrium E * is globally attracting.Further, E * is globally asymptotically stable.Now we concentrate on the stability of coexistence equilibria E i (N i , P i ), i = 1, 2. The Jacobian matrix of system Equation (3), at E i (N i , P i ), is Theorem 5.If the coexistence equilibrium E 2 (N 2 , P 2 ) of system Equation (3) exists, it is unstable.
For the coexistence equilibrium E 1 (N 1 , P 1 ), similarly, the determinant of J 1 is but since we obtain det(J 1 ) > 0, and E 1 (N 1 , P 1 ) may be node, focus or center.By Equations ( 13)-( 15), the sign of the trace of matrix J 1 is determined by According to Equation (28), we obtain Summarizing the above discussions, we have the following results on the stability of equilibrium E 1 (N 1 , P 1 ).Theorem 6. Suppose that the coexistence equilibrium E 1 (N 1 , P 1 ) exists and if one of following conditions is satisfied:    2, where the horizontal line denotes the curve of the predator with predator-extinct equilibrium E 0 .Two positive equilibria E 1 and E 2 arise simultaneously at h = h 1 when the parameter h decreases.Remark 3. In Figure 2, we consider how to set up the harvesting threshold value h.We find that h 1 is an important harvesting amount.If the harvesting threshold value exceeds h 1 , the system does not have a positive equilibrium; that is, the predator eventually tends to extinction.If the harvesting threshold value is less than h 1 , the system has two positive equilibria among which one is unstable and the other may be stable, i.e., the predator and prey may coexist.A bifurcation diagram is illustrated in Figure 3, where the horizontal blue line presents the curve of the predator with the positive equilibrium E * .It displays that there is a bifurcation at h = h 1 when the parameter h reduces, which produces three equilibria E * ,E 1 and E 2 .The existence of limit cycles plays an important role in determining the dynamical behavior of the system.For example, if there is no limit cycle in system Equation (3) and its positive equilibrium is unique and locally asymptotically stable, then it must be globally stable.Now, we explore the existence of limit cycles in system Equation (3).Theorem 7. Suppose n > 1 and h < min{h 1 , h 2 }.If ϕ > 0, then system Equation (3) has at least a stable limit cycle which encircles E 1 .
Proof.For n > 1 and h < h 2 , it is known from Theorem 2 a) that the equilibrium E * of system Equation (3) does not exist.Furthermore, because n > 1, h < h 1 , and h < h 2 , it follows from Theorem 2 b) that the equilibrium E 2 of system Equation (3) does not exist, but the equilibrium E 1 exists.
It follows from ϕ > 0 that E 1 is an unstable focus or node.It is easy to see that the unstable manifold at the saddle point E 0 (K, 0) is in the first quadrant.As the set D is positively invariant for system Equation (3), and system Equation (3) does not have any equilibrium in the interior of D \ {E 1 }.It follows from the Poincar é-Bendixson theorem that system Equation (3) has at least a stable limit cycle which encircles E 1 .
In general, Dulac functions are only applied to smooth vector fields in the study of nonexistence of limit cycles.Since the right-hand sides of Equation ( 3), denoted by f 1 and f 2 , are not smooth, following the similar arguments as in Wang ([31] Lemma 3.2), which is based on Green s Theorem, we are able to obtain sufficient conditions for the nonexistence of limit cycles in system Equation (3).Theorem 8. System Equation (3) does not have a limit cycle if caK < d + r. 3), it is easy to see that the positive solutions of Equation ( 3) eventually enter and remain in the region C = {(N, P) : N ≤ K} Thus, if a limit cycle exists, it must lie in the region C. Take a Dulac function F = 1

Proof. By the first equation of Equation (
N .Then we have Hence, system Equation (3) does not have a limit cycle.

Numerical Simulation
In this section, we present numerical examples for system Equation (3).Selecting P 0 = 0.25, we get h = 0.0075 > h 2 .The equilibrium E * (0.2, 0.24) exists, but E 1 and E 2 do not exist (Theorem 2 b)).Its phase portrait is given in Figure 4, which shows that the unique positive equilibrium E * is globally asymptotically stable.3) when E * is globally asymptotically stable and E 0 is unstable.
If we choose P 0 = 0.2, then h = 0.006 < h 2 .Equilibrium E 1 (0.16044, 0.27165) exists, but E 2 and E * do not exist (Theorem 2 b)).Its phase portrait is given in Figure 5.The unique positive equilibrium E 1 is globally asymptotically stable in D.  3) when E 1 is locally asymptotically stable and E 0 is unstable.
The equilibrium E * in Figure 4 corresponds to some point on the curve of P * in Figure 1 (h > h 2 ), and the equilibrium E 1 in Figure 5 corresponds to some point on the curve of P 1 in Figure 1 (h < h 2 ).If we choose P 0 = 0.1, then P 0 < P. Equilibria E 1 (0.2, 0.5) and E 2 (0.325, 0.1875) exist, but equilibrium E * does not exist.Its phase portrait is illustrated in Figure 6.It shows that equilibria E 0 and E 1 are asymptotically stable.
The equilibria E 1 and E 2 in Figure 6 correspond to some points on the curves of P 1 and P 2 , respectively, in Figure 2 (h < h 1 ).3) is shown in Figure 7. Equilibria E * and E 1 are asymptotically stable, and E 0 and E 2 are unstable.So model Equation ( 3) has bistable positive equilibria E * and E 1 .It follows from Figure 7 that the stable region Γ s of the saddle point E 2 divides the positive invariant set into two regions.The attractive basin for the stable equilibrium E 1 is the region above Γ s and the region below Γ s is the basin of attraction for the stable equilibrium E * .
The equilibria E 1 and E 2 in Figure 7 correspond to some points on the curves of P 1 and P 2 , respectively, and E * corresponds to some point on the curve of P * in Figure 3 (h 2 < h < h 1 ).

Example 4.3 (II).
We set parameter values r = 0.03, K = 0.25, a = 0.5, c = 0.4, d = 0.01, m = 0.03, and P 0 = 0.0133.Thus we have n = 1.25, uh 1 = 0.00048, h 2 = 0.00036, P = 0.024, and h = 0.000399.These parameter values satisfy conditions a) and c) in Theorem 2 and the condition c) in Theorem 6.Thus, equilibria E 0 (0.25, 0), E * (0.The model with the parameter values in Example 4.3 (II) has a similar bifurcation diagram as that in Example 4.3 (I), and the stability of equilibria in Figure 8 is identical with that in Figure 7. Choosing P 0 = 0.008, we have h = 2 * 10 −5 .Thus, h < h 2 < h 1 , P 0 < P, caK > r, h > h 3 , and η < 0. Equilibrium E 1 (0.04069, 0.03186) exists but is unstable, and E * and E 2 do not exist.The parameter values satisfy the conditions of Theorem 7. Its phase portrait is given in Figure 9, which shows that model Equation (3) has a stable limit cycle which encircles E 1 .3) is shown in Figure 10.In Figure 10, equilibrium E 1 is stable and two periodic orbits encircle E 1 .We can see that the outside periodic orbit is stable.However, the inside periodic orbit is unstable.A trajectory (the dotted line) between the outside and inside periodic orbits ultimately tends to the outside periodic orbit, but a trajectory (the thin black line) starting from within the unstable periodic orbit finally tends to equilibrium E 1 .So the initial state is important for ultimate trends of trajectories.

Conclusions
In this paper, we proposed and studied a new predator-prey model with non-smooth switched harvest on the predator.If the density of the predator is below a switched value, the harvest has a linear harvesting rate.Otherwise, the harvesting rate is constant.Our model exhibits new dynamical features compared to those with a linear harvesting rate or a constant harvesting rate.
According to the Kolmogorov Theorem [38], under certain assumptions, the model with a linear predator harvesting rate has either a stable equilibrium or a stable limit cycle, whereas the model with a constant harvesting rate on the predator has richer dynamics [16,23,23].For example, for a class of predator-prey systems, Brauer and Soudack [23,23] obtained different types of dynamics for which the harvesting was in prey or a predator; Xiao and Jennings [16] further studied a ratio-dependent predator-prey model with a constant harvest on prey.They proved that the model could exhibit complicated bifurcation phenomena, including the Bogdanov-Takens bifurcation of cusp type, the heteroclinic bifurcation, or a separatrix connecting a saddle-node and a saddle bifurcation.
For the model studied in this paper, we showed that (see Theorem 2 and Corollary 1) a backward bifurcation from the predator-prey coexistence equilibrium may occur, which shows that reducing the net reproductive number of the predator to less than unity is not enough to eradicate the predator.On the other hand, when the net reproductive number of the predator is greater than unity, we showed that the predator always coexists with the prey permanently (Theorem 3), and the model may exhibit the following dynamics: (i) a unique globally asymptotically stable coexistence equilibrium; (ii) the coexistence of positive saddle equilibria connecting with either a locally asymptotically stable positive equilibrium (biostable) or a limit cycle; (iii) two stable positive equilibria coexisting with a saddle point.By numerical examples, we also showed that the model could exhibit more new dynamical features: (a) a limit cycle encircling a unique positive equilibrium (see Figure 9); (b) two cycles surround an identical positive equilibrium, with one stable and one unstable (see Figure 10).
We would like to point out that we have assumed the simple functional response of the bilinear type in our current model Equation (3).We may also consider other types of functional responses.The dynamics may be richer and more complex.Further investigations are planned in our future studies.

Corollary 1 .Example 3 . 1 .
(a) caK ≤ r; (b) caK > r and h ≤ h 3 ; (c) caK > r and h > h 3 and η > 0, then E 1 (N 1 , P 1 ) is locally asymptotically stable.It is unstable if caK > r, h > h 3 , and η < 0.Combining the existence and stability results of equilibria of model Equation (3), by following similar arguments to ([31] Corollary 2.3), we present the following corollary to give conditions for bifurcation.If n < 1 and h < min{h 1 , m P}, then system Equation (3) has a backward bifurcation of positive equilibria.Next, we give examples to demonstrate the bifurcation of multiple equilibria.For various parameter values, model Equation (3) has a forward bifurcation from one positive equilibrium to another positive equilibrium (see Example 3.1.)and a backward bifurcation with a predator-extinct equilibrium and two positive equilibria (Example 3.2.).Note that the conditions in Theorem 2 a) and Theorem 2 c) guarantee the existence of three positive equilibria E * and E 1 , E 2 (Example 3.3.).Using the following parameter values r = 0.1, K = 0.5, a = 0.25, c = 0.8, d = 0.01, and m = 0.03, we obtain n = 2.5, and h 2 = 0.0072.When h ≤ h 2 , a bifurcation diagram is shown in Figure 1.When the parameter h decreases, the bifurcation at h = h 2 is forward, and model Equation (3) has a unique positive equilibrium for h > 0, which is similar to ([30] Example 2.4). 0

1 Figure 2 .
Figure 2. The backward bifurcation diagram of P 1 and P 2 versus h for Equation (3).The solid line with P 1 and dotted line with P 2 represent the curves of the predator with coexistent equilibrium E 1 and E 2 , respectively.

Figure 3 .
Figure 3.The bifurcation diagram with P 1 , P 2 and P * versus h.The lines with P 1 and P 2 indicate the curves of the predator with coexistent equilibria E 1 and E 2 , respectively, and the line with P * indicates the curve of the predator with coexistent equilibrium E * .

Figure 4 .
Figure 4.The phase portrait of model Equation (3) when E * is globally asymptotically stable and E 0 is unstable.

1 Figure 5 .
Figure 5.The phase portrait of model Equation (3) when E 1 is locally asymptotically stable and E 0 is unstable.

Figure 7 .
Figure 7.The phase portrait of model Equation (3) with bistable positive equilibria E * and E 1 , and unstable equilibria E 0 and E 2 .

2 Figure 10 .
Figure 10.The phase portrait of model Equation (3) when E 1 is stable.
Choosing r = 0.1, K = 0.4, a = 0.1, c = 0.8, d = 0.01, and m = 0.03, we obtain n = 0.8 < 1, h 1 = 0.00378, and P = 0.34367.A backward bifurcation diagram is given in Figure * P 2 Figure 1.The forward bifurcation diagram from P * to P 1 versus u for Equation (3).The line with P 1 indicates the curve of the predator with coexistent equilibrium E 1 , and the line with P * indicates the curve of the predator with coexistent equilibrium E * .Example 3.2.