Multiple limit cycles in a Leslie – Gower-type predator – prey model considering weak Allee effect on prey ∗

In this work, a modified Leslie–Gower-type predator–prey model is analyzed, considering now that the prey population is affected by a weak Allee effect, complementing results obtained in previous papers in which the consequences of strong Allee effect for the same model were established. In order to simplify the calculations, a diffeomorphism is constructed to obtain a topological equivalent system for which we establish the boundedness of solutions, the nature of equilibrium points, the existence of a separatrix curve dividing the behavior of trajectories. Also, the existence of two concentric limit cycles surrounding a unique positive equilibrium point (generalized Hopf or Bautin bifurcation) is shown. Although the equilibrium point associated to the weak Allee effect lies in the second quadrant, the model has a rich dynamics due to this phenomenon, such as it happens when a strong Allee effect is considered in prey population. The model here analyzed has some similar behaviors with the model considering strong Allee effect, having both two limit cycles; nevertheless, they differ in the amount of positive equilibrium points and the existence in our model of a non-infinitesimal limit cycle, which exists when the positive equilibrium is a repeller node. The main results obtained are reinforced by means of some numerical simulations.


Introduction
One of the main problems in Population Dynamics is the study of the dynamical behaviour of predator-prey interactions.Varied models have been proposed for this interaction, being one of them established in the seminal work by P.H. Leslie in 1948 [28] and discussed by Leslie and Gower in 1960 [29], which is not conforming to the Lotka-Volterra framework [38].
The principal characteristic of Leslie [32] or Leslie-Gower-type model [6] is that the equation describing the predator growth is based on a direct generalization of the logistic equation [38], assuming that the environmental carrying capacity of predators K y is directly proportional to the available prey quantity [38], that is, K y (x) = nx, where x = x(t) indicates the prey population size at time t 0 [2,3].
This model retains the prey equation from the Volterra model [38] and the functional response, the change in the density of prey attacked per unit time per predator, as the prey density changes [32], is expressed by the function h(x) = qx with q > 0. Implicitly, this assumption considers that the predator subsistence exclusively depends on prey population, i.e., the predator is specialist [38].
Meanwhile, the predator population growth rate depends on the predator population size ratio to the prey population size, a positive growth rate is predicted when the absolute prey and predator densities are significantly low [38].In spite of anomaly, these models have been employed to describe a real predator-prey interaction without considering this objection [25].
It is well known, the Leslie-Gower model has a unique equilibrium point positive (in the interior of the first quadrant), which is globally asymptotically stable [38].Nonetheless, it can assume that predators are generalist having an alternative food when its preferred prey is scarce.So, the environmental carrying capacity is given by K y (x) = nx+c with c > 0; it is said that the model is represented by a Leslie-Gower scheme, also known as modified Leslie-Gower model [4,6].
In this work, the mathematical results obtained in a previous paper [21] are complemented in which the consequences of strong Allee effect for the Leslie-Gower-type predator-prey model were established.Any mechanism leading to a positive relationship between a component of individual fitness and the number or density of conspecifics can be named as a Allee effect [37].
It describes a biological scenario characterized by a positive correlation between the population sizes and its per capita growth rate at low population sizes, increasing their likelihood of extinction [13]; it has been denominated in different ways in Population Dynamics [30] and depensation in Fisheries Sciences [11,17,30].
The main characterization of Allee effect is that the per capita growth rate is positive for low population sizes.This is a common phenomenon in some animal populations and various mechanisms have been proposed as potential sources of this biological phenomena (Table 1 in [8] or Table 2.1 in [12]).
Recents ecological research suggest that two or more Allee effects can lead to these mechanisms acting simultaneously on a single population (Table 2 in [8]); the combined influence of some of these phenomena is known as multiple (double) Allee effect [8,20].
This phenomenon has become crucial for population dynamics since in fact it has a surprising number of ramifications towards different branches of ecology [12], and the knowledge of this effect on simple models is essential to understand more complicated ones, specially when the weak Allee effect is considered [31].
To model the growth rate of a population affected by the Allee effect, the most simplest mathematical form will be employed, which is described by the cubic polynomial https://www.mii.vu.lt/NA differential equation [7,24,26] where −K < m K and r, K > 0. When m > 0, the population growth rate decreases if the population size is below the threshold level m and the population goes to extinction.In this case, equation (1) describes the strong Allee effect [39][40][41].If m 0, it is said that the population is affected by a weak Allee effect [11,18,41].In Fisheries Sciences, the same phenomena are called critical and pure depensation, respectively [11,30,41].
The function f (x) = dx/(x dt), expresses the per capita growth rate; clearly, for −K < m K.This function grows as the population size increases from x 0 until it attains a maximum value; then, it decreases for larger population size, as the logistic growth function.Notice that if m < −K, the per capita growth rate f (x) is negative for all x 0, and the growth function is of compensatory type [11] as the logistic growth function and equation ( 1) ceases of to represent the weak Allee effect [23].But, mathematically the added factor (x − m) have influence in the dynamics of this equation [26].
Populations with weak Allee dynamic experimenting lower per capita growth rates at low densities, but never lead to negative per capita growth rates and therefore have no critical threshold to be exceed for population to survive [34].To our knowledge no other work studies the weak Allee effect considering m < 0 in equation (1).Other articles mentioning weak Allee effect in continuous time models, only consider the special case m = 0 [22] or take other expression to describe this effect [27,31,43].
Nevertheless, the study of weak Allee effect is usually neglected since the per capita population growth rate remains positive for any population sizes such as it happens with the logistic growth [23].
Moreover, there is relatively limited evidences of strong Allee effects in natural populations, and thus for a key part of the theoretical framework for Allee effects [34] and a weak Allee effect may appear over a given range of densities [41]; so, the influence on dynamics of model is especially important to study.
Specifically, we would like to know how the weak Allee effect affect the population dynamics of both prey and predator; particularly if the weak form of this phenomenon can help to both populations avoid the extintion.
Different mathematical forms have been proposed to model the Allee effect [8] and most of them are topologically equivalent [19].The consequences of this different mathematical formalizations has been studied for Gause-type predator-prey models in [20,22,35,43,44] and for Leslie-Gower in [21,33].Particularly, in a previous work, we have proved that by considering two different ways for describe the Allee effect, in a same predator-prey model, changes in the quantity of limit cycles surrounding a positive equilibrium point can occur [20].We consider, the determination of the number of limit cycles bifurcating from a weak focus is a good criterion for the classification of predatorprey models [23].
The rest of this paper is organized as follows: The proposed model is presented in the next section; the main properties are established in Section 3, and its consequences are shown in Section 4. In addition, simulations are given to verify the main results.

The model
When both the Allee effect phenomenon on prey and the functional response are modelled in the most simple way, the Leslie-Gower predator-prey model is given by the following bidimensional differential equation system: where x = x(t) and y = y(t) represent the prey and predators population size, respectively, for t 0 (measured as biomass, number of individuals or density by unit of area or volume) with µ = (r, K, q, s, n, m) ∈ R 5 + ×]−K, K[ and the parameters having different biological meanings [21].
In this work, for the analysis of system (2), we will be consider the prey population is affected by a weak Allee effect assuming that −K < m < 0, i.e., m ∈ ] − K, 0].The particular case of weak Allee effect when m = 0 was studied in [21]; thus, the present work complements the results obtained in our previous paper [21].Notice that if predator is generalist (considering c > 0 in the predator carrying capacity K y ), the model would have a new equilibrium point (0, c) over the y-axis; the implications of the added parameter will be not studied here.
For the strong Allee effect, conditions such that an equilibrium point corresponding to a positive steady state may undergo saddle-node, Hopf and Bogdanov-Takens bifurcations are given in [1] for Leslie-Gower-type predator-prey model.There the Allee effect is modelled by means of a general term that multiplies the classic logistic growth rate in the prey equation (multiplicative Allee effect) as equation (1).Nevertheless, we orientate our analysis to describe the behavior of system (2) by a classification in the parameter space (diagram bifurcation), mainly interested in to establish the number of limit cycles encircling a positive equilibrium point when m < 0.
It is clear that model ( 1) is a Kolmogorov-type system [15] not defined at the y-axis, particularly at the point (0, 0), but both isoclines of system (1) pass through this point; however, (0, 0) is a point of particular interest having a strong incidence on behavior of system (3) or vector field X µ , which is defined on the set: The equilibrium points of system (2) in Ω are (K, 0) and (x e , nx e ), where x e satisfies the quadratic equation Then, the equation can have positive roots if and only if Kr + mr − Knq > 0 and δ = (Kr + mr − Knq) 2 − 4Kmr 2 0, or none if Kr + mr − Knq 0. We note that (m, 0) lies in the second quadrant.
In order to simplify the calculations, we make a change of variables and a time rescaling [20,24,33,36] described in the following.Proposition 1. System (2) is topologically equivalent to the polynomial system given by Z η (u, v): Proof.Let x = Ku and y = nKv the new dimensionless variables.Substituting into system (1) and simplifying, we obtain
Remark 1. 1.We have constructed the diffeomorphism [10,14] The Jacobian matrix of ϕ is and we have that det Dϕ(u, v, τ ) = nKu/r > 0.Then, the diffeomorphism ϕ is a smooth change of variables with a rescaling of the time preserving the time orientation; thus, the vector field X µ (x, y), is topologically equivalent to the vector field )∂/∂v and the associated differential equation system is given by the polynomial system (3) of fourth degree and of Kolmogorov type.
2. This parameterization and time rescaling permits us to obtain a representative system, with the least amount of parameters possible, which describes the dynamical behaviors of all those systems that are topologically equivalent to system (2).Then, more important than knowing the influence of a particular parameter in the dynamical behavior of system (2), it is more useful the knowledge of the relationships between some of them for a better view of their implications in the properties of system.The singularities of system (3) in Ω are (0, 0), (1, 0) and (u e , u e ), where u e satisfies the equation The point (M, 0) lies in the second quadrant, and it has no ecological interest.

Main results
For the vector field Z η or system (3) with M < 0, we have Lemma 1.
The solutions are bounded.
Proof. 1. System (3) is of Kolmogorov type, then the axis are invariant sets.Considering u = 1, we have that where du/dτ < 0, anything the sign of dv/dτ , all trajectories point to Γ .
2. We define the new function w(τ ) = u(τ ) + v(τ )/S; the time derivative of w along a solution of system (3) is ) with α a parameter for determine.
Collecting the expression in the right of the last equality respect to u and after some simple algebraic manipulations, we obtain Then, considering 0 u 1 and as v satisfy the second equation of system (3), it has 0 v 1; thus, ).Thus, there exists δ > 0 such as being a first-order linear inequality.
Applying the theorem on differential inequality [9], we obtain Clearly, when τ → ∞, 0 w(u, v) Qδ/M and the solutions are bounded.
The result above assures that there exists a set such that B is the region in which all the solutions of system (3) with initial conditions in Ω are confined.This property assures that the model is well-posed, i.e., it is no possible the existence of many predators and a little prey population size, which would be a counterintuitive situation, not sustainable ecologically.We recall that if M 0, system (3) can have until two positive equilibrium points [21].The number of positives singularities of system (3) differs when M < 0 or M 0 as is demonstrated in the following lemma.Lemma 2. There exists a unique positive equilibrium point.From each case above, system (3) has only one positive equilibrium point.
The Jacobian matrix of system (3) necessary to obtain the nature of equilibrium points is 3.1 Nature of equilibrium points Lemma 3.For all parameter values, the point (1, 0) is a hyperbolic saddle point.
Proof.Evaluating the Jacobian matrix, we obtain Lemma 4. The point (0, 0) is a non hyperbolic equilibrium point having a hyperbolic and a parabolic sector determined for the line v = ((S + M )/S)u, i.e., there exists a separatrix curve Σ in the phase plane that divides the behavior of trajectories.
Proof.The Jacobian matrix is the null matrix and the proof follows as that given for the case M > 0 [21].
Notice that the equilibrium (0, 0) of the vector field Z η is a non hyperbolic attractor point.Moreover, the straight line v = ((S + M )/S)u is tangent to the separatrix curve in the phase plane dividing the behavior of the trajectories; thus, the point (0, 0) is an attractor point for certain trajectories and a saddle point for others, but the trajectories above this separatrix have different α-limit sets.
As M < 0, the slope of this line is positive for S > −M > 0 and negative for 0 < S < −M ; in the last case, lies in the second quadrant.

Lemma 5.
For an open subset of parameter values, there exists a heteroclinic cycle γ h in the first quadrant containing the equilibria (0, 0) and (1, 0).
Clearly, the α-limit of W s (0, 0) and the ω-limit of W u (1, 0) are bounded in the direction of the v-axis.Neither the ω-limit of W u (1, 0) is on the u-axis.Let u * be such that 0 < u * < 1.Then, there are points (u * , v s ) ∈ W s (0, 0) and (u * , v u ) ∈ W u (1, 0) https://www.mii.vu.lt/NA with v s and v u depending on the parameter values such that v s = s(Q, S, M ) and It is easy to see that, if 0 < u * 1, then v s < v u and the stable manifold W s (M, 0) is below the unstable manifold W u (1, 0).If M u * < 1, then v s > v u and the stable manifold W s (M, 0) is above the unstable manifold W u (1, 0).Since the vector field Y η is continuous respect to the parameters, the stable manifold W s (0, 0) intersect the unstable manifold W u (1, 0) for some parameter values.
Hence, there exists Moreover, by uniqueness of solutions of system (3), this intersection must occur along a whole trajectory γ 1M .Therefore, the equation s(Q, S, M ) = u(Q, S, M ) defines a codimension-one submanifold in the parameters space for which the heteroclinic curve γ 10 exists in R 2 + , connecting the points (1, 0) and (0, 0).Then, γ 10 ⊂ W s (0, 0)∩ W u (1, 0) and it lies entirely on a segment of the u-axis and exists for any parameter value such that 0 < H < 1.It follows that a heteroclinic cycle γ h exists for certain parameter values on the same submanifold.More precisely, γ h = (1, 0) ∪ γ 10 ∪ (M, 0) ∪ γ 01 .Let (u s , v s ) ∈ Σ with Σ the separatrix curve determined by the stable manifold of the non hyperbolic equilibrium point (0, 0).Now, we define the set i.e., Λ is determined by the separatrix curve Σ, the u-axis and the straight line u = 1.Thus, the region Λ is compact.This implies that Poincaré-Bendixon theorem is applicable there, which will determine the nature of the unique equilibrium point (H, H) with H = (A + W )/2; the Jacobian matrix is The sign of tr DZ η (H, H) depends on the factor Moreover, to determine if the positive equilibrium point is node or focus, it is necessary to known the sign of the quantity P defined by Nonlinear Anal.Model.Control, 22(3):347-365 Figure 1.For 0 < H < 1, the two possible relative positions between the stable manifold W s (0, 0) of singularity (0, 0) and the unstable manifold W u (1, 0) of saddle point (1, 0) are shown.
In Fig. 1 we show the relative positions of the stable manifold W s (0, 0) and the unstable manifold W u (1, 0) described in the above lemma.2. Supposing that v s < v u ; then (H, H) is a repeller focus or node and the equilibrium point (0, 0) or else, a non-infinitesimal stable limit cycle, can be ω-limit of all trajectories of the system.

(i) (H, H) is an attractor equilibrium point if and only if S > H(1 + M − 2H).
(ii) (H, H) is a repeller equilibrium point if and only if S < H(1 + M − 2H).As the region Λ is compact, by Hopf Bifurcation Theorem, a stable limit cycle is generated at this point.The sign of P determines if (H, H) is a focus or a node.
2. Supposing that v s < v u .In this case, the stable manifold W s (0, 0) lies under the unstable manifold W u (1, 0).Then, T > 0 and the point (H, H) is a repeller node.By the Existence and Uniqueness Theorem, all trajectories with initial conditions near of the point (H, H) can not intersect W u (1, 0) nor W s (0, 0).Since (1, 0) is a saddle point and (0, 0) is a non hyperbolic saddle point, the trajectories have the origin (0, 0) as their ω-limit, and this point is globally asymptotically stable.
Moreover, if T > 0, the limit cycles that appear by Hopf bifurcations grows until disappearing when the heteroclinic cycle γ 10 joining points (0, 0) and (1, 0) is attained; later, this heteroclinic is broken; then, there exists a subset on the parameters space for which (0, 0) is globally asymptotically stable.
So, a trajectory exists, which is originated on (H, H) and ending on (0, 0), forming a new heteroclinic curve.
The above analysis considers all cases, so that the statement is proved.
1.At least one infinitesimal limit cycle bifurcates when T changes of sign (see Fig. 2).Moreover, as v s decreases towards v u , this periodic orbit gets larger.At the limit when v s = v u , this periodic orbit becomes the heteroclinic cycle γ h from Lemma 3 and then disappears as 2)), the point (H, H) is an attractor point since T > 0 and tr DZ η (H, H) < 0 and there are no limit cycles.
1.There exist conditions on the parameter values for which an homoclinic curve is determined by the stable and unstable manifold non-hyperbolic equilibrium point (0, 0). 2. There exists a non-infinitesimal limit cycle [10,16] that bifurcates from the homoclinic surrounding the equilibrium point (H, H).
Proof. 1.Let W u (0, 0) and W s (0, 0) be the unstable and the stable manifold of the nonhyperbolic equilibrium point (0, 0).The trajectory determined by W u (0, 0) cannot cross the line u = 1 towards the right since Γ is an invariant region.The α-limit of stable manifold W s (0, 0) must be the point (H, H) or lies at infinity over u-axis.
Then, the ω-limit of stable manifold W u (0, 0) must be: (i) The point (H, H) when this point is an attractor or else a stable (attractor) limit cycle.(ii) The non-hyperbolic equilibrium point (0, 0) when the point (H, H) is a repeller.(iii) The trajectory determined by the unstable manifold W u (1, 0) of the hyperbolic equilibrium point (1, 0).
Therefore, there exists a set of parameter values for which the trajectories determine by W u (0, 0) and W s (0, 0) intersect forming a homoclinic curve.
2. Is straightforward by introduction of a single cross section, transverse to the homoclinic solution, that if this homoclinic is broken, then there exists a subset on the parameters space for which a stable non-infinitesimal limit cycle appears.A attracting singleround periodic orbit bifurcates inside the invariant region determined by the homoclinic loop.

Multiple limit cycles
The problem of determining the quantity of limit cycles, which can be born throughout the bifurcation of a weak (fine) focus [10], is not an easy task, and it is related to solving the well-known Hilbert 16th Problem (proposed by the mathematician David Hilbert in 1900) on the maximum number and relative position of limit cycles [16].Yet this question remains unsolved in the case of a two-dimensional polynomial differential equation system, which degree must be less than or equal to p ∈ N and p 2 [16], existing an example with four limit cycles for p = 2 [16].
To determine the weakness of the focus (H, H), i.e., the number of limit cycles bifurcating of a weak focus [10,14], in the next theorem, we will use the method for calculations of the Lyapunov numbers (or quantities) [10].
Theorem 3. The singularity (H, H) of vector field Z η is at least a third-order weak focus if and only if Proof.System (3) can be expressed by where then the new system translated to origin of coordinates system is Zν : and the Jacobian matrix of system (5) at the point (0, 0) is The Jordan matrix associated to vector field Z ν is [5] Then, using the matrix change of basis [5] B = 0 . https://www.mii.vu.lt/NAAfter a large algebraic calculations, we obtain Zν : By means of a time rescaling given by T = Lτ , we derive to the normal form [10] given by Zν : Using the Mathematica package [42], we obtain that the second Lyapunov quantity [10] is Clearly, f 2 (M, H) > 0 for H > 1/6 and η 2 < 0; then, (0, 0) is a one-order weak focus of vector field Zν ; thus, a unique limit cycle exists in vector field Z η for η 1 = tr DZ η (H, H) > 0 (the point (H, H) is a repeller) and η 2 < 0.
We note that f 2 (M, H) can change of sign for H < 1/6; then, the equilibrium (0, 0) of vector field Zν is surrounded by at least two limit cycles.If (0, 0) is a third-order weak focus of vector field Zν , then η 3 must be the first non-zero Lyapunov quantity.We obtain where The sign of f 3 (M, H) could change, and therefore, η 4 must be calculated to determine the existence of a third-limit cycle.However, such task will not be tackled in this work.When f 3 (M, H) < 0, η 3 is positive.If tr D Zν (0, 0) < 0, f 2 (M, H) > 0 and η 3 is negative, the point (0, 0) of the vector field Zν is a third-order weak focus.Then, the equilibrium point (H, H) is at least a third-order weak focus of the vector field Z η or system (3).
Corollary 1.There exists at least two limit cycles surrounding the singularity (H, H) of vector field Z η , the innermost unstable and the outermost stable if and only if S < H(1 Proof.We note that the first condition implies the negativity of tr DZ η (H, H), implying the point (H, H) is a local attractor; the second inequality determines the positivity of η 2 .Thus, when η 3 is negative, there exist two limit cycles.Thus, the Lyapunov quantities have different sign, and the two limit cycles are the innermost unstable and the outermost stable.
Remark 3.Although the calculus of weakness of the focus (H, H) requires two stringent relationships between model parameters, a tiny deviation in one of them causes the condition for the existence of one or more limit cycles.So, these relations determine a subset of measure non-zero for the existence of a unique or two limit cycles.This is important biologically, inasmuch as in reality none of the equalities given in Theorem 11 will be possible maintain by long time, due that any little change in any of the involved parameters will mean inequality rather than equality, and hence, we are working here with a sort of structurally unstable system.
Corollary 12 implies a significant difference with the model in which the strong Allee effect is considered since this system has two limit cycles in opposite stability, i.e., the innermost stable and the outermost unstable [21].

Some simulations
In the following, some simulations are given to reinforces the results obtained for system (3).

Discussion
In this work, we have analyzed a simple Leslie-Gower predator-prey model, considering the phenomenon known as weak Allee effect showing a significant impact in the dynamics of model.The main results have been expressed in funcion of the parameter S = s/(rK), indicating as the relation among r, s and K is important in the dynamical behavior of system (2).
It is well known that Leslie-Gower model is structurally stable as it has a unique equilibrium point at the interior of the first quadrant that is globally asymptotically stable.Ecologically speaking, this means the prey and predator species coexist for any parameter values.
We have modified the Leslie-Gower model by introducing a new factor, always positive, in the prey growth rate to describe the weak Allee effect and a compensatory curve of growth is obtained [11] as logistic growth function.However, the new model here studied has a rich and varied dynamics, where two of them are not topologically equivalent, unlike the original Leslie-Gower model.
We have shown the importance of the singularity (0, 0) in the dynamics of the modified Leslie-Gower model, demonstrating the existence of a bistability phenomenon since (0, 0) is always an attractor in addition to either a stable positive equilibrium point or a stable limit cycle for determined parameter values.
Furthermore, the inclusion of weak Allee effect has a significant consequence in the dynamics of system since until two limit cycles can appear surrounding the unique positive equilibrium point.Limit cycles do not exist in the original Leslie-Gower model since the unique positive equilibrium point is globally asymptotically stable for a wide set of parameter values.
For system (3), we proved the existence of a separatrix curve determined by the stable manifold of point (0, 0) dividing the behavior of trajectories.This implies that the model is highly sensitive to the initial conditions.Moreover, for a wide set of the parameter values, the model can have two different attractors, and there exist solutions having the origin as its ω-limit.This implies in system (2) the possibility of depletion for both species is high, yet the ratio x/y is small.
We have shown that the weak Allee effect may be a destabilizing force in predatorprey models since the equilibrium point of the Leslie-Gower model changes from stable to unstable or otherwise.The analyzed model involves a difference the special case of weak Allee effect when m = 0 [21], respect to the quantity of limit cycles bifurcating of a fine (weak) focus.System (3) has a unique positive equilibrium point when m < 0, which is surrounded by two limit cycles, for certain parameter values, the innermost unstable and the outermost stable.Meanwhile, system (3) has a unique limit cycle, when m = 0 [21], which can disappear when the other parameters vary; however, the trajectories of system (3) can also have different ω-limit sets strongly depending on the initial conditions.
The Leslie-Gower-type predator-prey model, considering strong Allee effect studied in [21] and the weak Allee effect here analyzed, have an similar behavior with respect to the non-uniqueness of the limit cycle, but differs in the amount of positive equilibrium points.Assuming M > 0, there may be up to two positive equilibrium points with the possibility of homoclinic orbits [21].
Moreover, for the model with weak Allee, we show by simulations, that a non-infinitesimal limit cycle [10] surrounding the unique positive equilibrium point exits; the proof of the existence of this limit cycle is left as an open problem.https://www.mii.vu.lt/NAFurthermore, it can see the varied and rich dynamic that has the modified Leslie-Gower model here studied in comparison with the Gause-type predator-prey model using the same function to model the Allee effect and the linear functional response, which is described in the Kot's book [26].
Taking the Leslie-Gower-type predator-prey model [32,38], interesting works could be analyzed in future in which: i) a generalist predators be considered, or ii) other expressions to the Allee effect be considered in the prey growth function.
In short, the consideration of weak Allee effect on prey population has a high impact on the dynamics of the classical Leslie-Gower model, resulting in a new system with various different mathematical behaviors, which have adequate biological meanings.Thus, considering the obtained results, we affirm that the number of limit cycles bifurcating from a weak focus, it is a good mathematical criterion for the classification of predatorprey models.