Global Dynamics of Leslie-Gower Competitive Systems in the Global Dynamics of Leslie-Gower Competitive Systems in the Plane Plane

: Global dynamic results are obtained for families of competitive systems of difference equations of the form x n + 1 =


Introduction
The system of difference equations, where the parameters b 1 , b 2 are positive real numbers, and α 1 , α 2 , c 1 , and c 2 and the initial conditions x 0 and y 0 are arbitrary non-negative numbers, is considered as a major discrete model that describes the competition of two species, see [1][2][3][4][5][6][7].In this paper we consider the effect of terms α 1 , α 2 on the global dynamics of system (1).The global dynamics of (1) was considered in the case where the parameters α 1 , α 2 are positive in [1][2][3]5] and the complete description of the dynamics was given in [5] where the following result was obtained: Assuming, without loss of generality, that α 1 = α 2 = 1, it has been shown in [1] that under the condition b 1 > 1 and b 2 > 1, the points are equilibria of Equation (1), and that for some values of the parameters there exists an additional equilibrium point E 3 located in the open positive quadrant, given by Important subsets of parameter space are described in Table 1, together with corresponding behavior of equilibrium solutions established in [1].
Table 1.Global behavior of solutions to Equation (1) when b 1 > 1 and b 2 > 1. Equality relations are not represented for the sake of a simpler description.
E 0 is a repeller E x is a saddle E y is a saddle E is an interior local attractor.
Case 2. E 0 is a repeller E x is a saddle E y attractor on [0, ∞) × (0, ∞) No interior fixed point exists.Theorem 4.1 in [5] applies when parameters vary from Case 1 to Case 4 of Table 1.Set and define T a to be the map of Equation ( 1) restricted to R = [0, ∞) × (0, ∞), that is, Therefore, Theorem 4.1 in [5] gives global behavior of solutions to system (1) on R = [0, ∞) × (0, ∞) for a ∈ A. In particular, a bifurcation occurs when the equilibrium E changes its local character from a locally stable equilibrium to a saddle point.This happens when the parameters cross the critical surface It is also shown in [1] that the open, positive semiaxis (0, ∞) × {0} is attracted to E x , and that the open, positive semiaxis {0} × (0, ∞) is attracted to E y .The following two results describe the global dynamics of system (1) in all cases.The first result gives the global dynamics in the hyperbolic case and the second result in the non-hyperbolic case.
is the graph of a continuous, increasing function of the first coordinate.Furthermore, a solution {x n } converges to E x whenever x 0 is above W s (E) in South-east ordering, and {x n } converges to E y whenever x 0 is below W s (E) in South-east ordering.(ii) Suppose that c , and E y attracts all points on the open semiaxis {0} × (0, ∞).
See Figure 1 for graphical interpretation.
. The non-hyperbolic case when was not considered in [1].When (3) holds, a direct calculation gives that the fixed points of T a are E 0 (0, 0) and all points on the segment E := {E t : 0 ≤ t ≤ 1 }, where The eigenvalues of the Jacobian of T a at E t are and corresponding eigenvectors are It is shown in [8] that, for system (1), the hypotheses of Theorem 5 in [6] are satisfied and that all solutions fall inside an invariant rectangular region.Therefore, every solution of (1) converges to an equilibrium point.A direct calculation shows that the origin is a repeller.We conclude that every nonzero solution converges to a point (x, y) ∈ E .Also, with an argument similar to the one used in [9], one has that the equilibrium depends continuously on the initial condition.That is, if T * (x, y) := lim T n a (x, y), then T * is continuous.These observations, together with an application of Theorem 1 in [6] lead to the following result.Every nonzero solution to system (1) converges to an equlibrium (x, y) ∈ E .(ii) For every (x, y) ∈ E with x = 0 and y = 0, the stable set W s (x,y) is an unbounded, increasing curve C with endpoint (0, 0).(iii) The limiting equilibrium varies continuously with the initial condition.
See Figure 1 for graphical interpretation.
Statement (ii) excludes equilibria of the form (0, y) and (x, 0) since the hypotheses of Theorem 1 in [6] are not satisfied at these points.
In this paper, we consider two related systems, namely and where all present coefficients are positive and the initial conditions are non-negative and such that x 0 + y 0 > 0. We derive the global dynamics of both systems (4) and ( 5), which explains the effect of the coefficients α 1 , α 2 on the global dynamics.Related systems are considered in [2,3,[10][11][12].The paper is organized as follows.Section 2 presents some basic preliminary results about competitive systems, which is our main tool in proving the results.Section 3 contains the global dynamics of system (4) and Section 4 gives the global dynamics of system (5).Actually we show that the global dynamics of all three systems is determined by their linearized dynamics.We use different techniques to prove the results for systems (4) and (5).In particular, we will use the fact that system (4) is homogeneous to obtain the exact equation of the global stable manifold of positive equilibrium solution of system (4).We will show that the presence of more parameters will create more dynamic scenarios in such a way that system (1) exhibits one additional dynamic scenario compared to system (5) which in turn will have one additional dynamic scenario more than system (4).This is in contrast with the global dynamics of second order difference equation with all non-negative parameters and initial conditions such that A + Bx n + Cx n−1 > 0, n = 0, 1, . .., where the most complicated dynamics occurs in the special case of Lyness' equation when γ = A = B = 0 [13].We will also show that systems (1) and ( 4) are solvable for the special values of parameters while there is not a formula for the exact solution of system (5) at this time for any choice of parameters.

Preliminaries
In this section we provide some basic facts about competitive maps and systems of difference equations in the plane.Definition 1.Let R be a subset of R 2 with nonempty interior, and let T : R → R be a map (i.e., a continuous function).Set T(x, y) = ( f (x, y), g(x, y)).The map T is competitive if f (x, y) is non-decreasing in x and non-increasing in y, and g(x, y) is non-increasing in x and non-decreasing in y.If both f and g are non-decreasing in x and y, we say that T is cooperative.If T is competitive (cooperative), the associated system of difference equations is said to be competitive (cooperative).The map T and associated difference equations system are said to be strongly competitive (strongly cooperative) if the adjectives non-decreasing and non-increasing are replaced by increasing and decreasing.
If T is differentiable, a sufficient condition for T to be strongly competitive is that the Jacobian matrix of T at any x ∈ R has the sign configuration Competitive systems of the form (6)  Denote with se the South-East partial order in the plane whose non-negative cone is the standard fourth quadrant {(x, y) : x ≥ 0, y ≤ 0}, that is, (x 1 , y 1 ) (x 2 , y 2 ) if and only if x 1 ≤ x 2 and y 1 ≥ y 2 .The North-East partial order ne is defined analogously with the non-negative cone given by the standard first quadrant {(x, y) : x ≥ 0, y ≥ 0}.
Competitive maps T in the plane preserve the South-East ordering: T(u) se T(v) whenever u se v. Similarly, cooperative maps in the plane preserve the North-East ordering.In fact, the concepts of competitive and cooperative (for maps) may be defined in terms of the order preserving properties of the maps.Thus the theory of competitive maps is a special case of the theory of order preserving maps (or monotone operators).
Order preserving maps in R n , and in particular competitive maps in R 2 , may have chaotic dynamics.Smale [19] showed that any continuous time vector field on the standard (n − 1)-simplex in R n can be embedded on a smooth, competitive vector field in R n for which the simplex is an attractor.In the case of a planar system (6), this means that any first order difference equation, including chaotic, can be embedded into a competitive system (6) in the plane.An effective method to do this is provided by Smith in [7].
Let be a partial order on R n with non-negative cone P. For x, y ∈ R n the order interval x, y is the set of all z such that x z y.We say x ≺ y if x y and x = y, and x y if y − x ∈ P • , where P • is the interior of P. A map T on a subset of R n is order preserving if T(x) T(y) whenever x ≺ y, strictly order preserving if T(x) ≺ T(y) whenever x ≺ y, and strongly order preserving if T(x) T(y) whenever x ≺ y.
Let T : R → R be a map with a fixed point x and let R be an invariant subset of R that contains x.We say that x is stable (asymptotically stable) relative to R if x is a stable (asymptotically stable) fixed point of the restriction of T to R .Definition 2. Let S be a nonempty subset of R 2 .A competitive map T : S → S is said to satisfy condition (O+) if for every x, y in S, T(x) ne T(y) implies x ne y, and T is said to satisfy condition (O−) if for every x, y in S, T(x) ne T(y) implies y ne x.
The following theorem was proved by deMottoni-Schiaffino for the Poincaré map of a periodic competitive Lotka-Volterra system of differential equations.Smith generalized the proof to competitive and cooperative maps [20,21].Theorem 3. Let S be a nonempty subset of R 2 .If T is a competitive map for which (O+) holds then for all x ∈ S, {T n (x)} is eventually componentwise monotone.If the orbit of x has compact closure, then it converges to a fixed point of T. If instead (O−) holds, then for all x ∈ S, {T 2n } is eventually componentwise monotone.If the orbit of x has compact closure in S, then its omega limit set is either a period-two orbit or a fixed point.

System (5) is an example of what might happen if the orbit of x has no compact closure in S.
The next two results are stated for order-preserving maps on R n .These results are known but are given here for completeness.See Theorem 2.1 in [5] and Corollary 1 in [5].See [15] for a more general version that is valid in ordered Banach spaces.There exists a fixed point of T in a, b .ii.
If T is strongly order preserving, then there exists a fixed point in a, b which is stable relative to a, b .iii.
If there is only one fixed point in a, b , then it is a global attractor in a, b and therefore asymptotically stable relative to a, b .
Corollary 1.If the non-negative cone of is a generalized quadrant in R n , and if T has no fixed points in u 1 , u 2 other than u 1 and u 2 , then the interior of u 1 , u 2 is either a subset of the basin of attraction of u 1 or a subset of the basin of attraction of u 2 .
Our main tool will be results from [4][5][6] regarding the existence of the global stable and unstable manifolds of competitive maps in the plane.
The non-hyperbolic equilibrium solution of system ( 6) is said to be of stable (resp.unstable) type if the second eigenvalue of the Jacobian matrix evaluated at the equilibrium solution is by absolute value less than 1 (resp.bigger than 1).

Global Dynamics of System (4)
First we give some basic results about the global behavior of system (4).Denote by the map associated with system (4).System (4) is homogeneous and was partially investigated in [22].

Lemma 1.
(a) Every solution of system (4) det J T (x, y) = 0 for every (x, y), where J T denotes the Jacobian matrix of the map T.
then the map T has an invariant line Proof.The Jacobian matrix J T of the map T has the form System (4) always has two equilibrium solutions on the axes, E x (b 1 , 0), E y (0, b 2 ).It can also have either exactly one interior equilibrium solution E or an infinite number of equilibrium solutions E t .Since the interior equilibrium solution E is an intersection of two equilibrium curves (isoclines) x-intercept of C 1 smaller than x-intercept of C 2 and y-intercept of C 1 bigger than y-intercept of C 2 ).These two geometrical conditions can be unified as condition (7).
Condition (7) implies that c 1 c 2 = 1, in which case the interior equilibrium E( x, ȳ) is given as: 2 summarizing the equilibrium points of system (4).

Condition
Equilibrium Points The following result describes the local stability character of all equilibrium solutions.

(a)
The equilibrium solution E x is locally asymptotically stable In each case, the eigenvectors associated with the eigenvalues The equilibrium solution E y is locally asymptotically stable In each case, the eigenvectors associated with the eigenvalues Proof.
(a) In view of (9), we have which implies that the eigenvalues of the Jacobian matrix are The corresponding eigenvectors are as stated.(b) In view of (9), we have which implies that the eigenvalues of the Jacobian matrix are The corresponding eigenvectors are as stated.(c) The eigenvalues of the Jacobian matrix evaluated at the equilibrium E, λ 1 = 0 and λ 2 = Tr(J T (E)), correspond to the roots of the characteristic polynomial p(t) = t 2 − Tr(J T (E))t.
Note that λ 2 > 0 by (9).Furthermore , then the equilibrium solution E x (resp.E y ) is locally asymptotically stable and E y (resp.E x ) is a saddle point.The basin of attraction of E x (resp.E y ) is the first quadrant of initial conditions without the positive part of the y-axis (resp.x-axis), which is attracted by E y (resp.E   Proof.

(a)
First we show that T does not have any period-two solutions.Our condition implies c 1 c 2 > 1.
By direct calculation one can show that a period-two solution satisfies the equation Please note that which means that both terms of such solution are complex conjugate and so there is no period-two solution in the first quadrant.
Taking into account that the Jacobian matrix evaluated at E has all non-zero entries, Theorem 5 of [4] implies the existence and uniqueness of both global stable and unstable manifolds W s (E) and W u (E) and so W s (E) = S .Furthermore, Theorem 5 of [4] implies that every (x 0 , y 0 ) below ] for some n ≥ N. In view of Corollary 1, T n ((x 0 , y 0 )) → E x .In a similar way, we can treat the case when (x 0 , y 0 ) is above S .(b) In view of Lemma 2 part (a), the eigenvectors which correspond to E x and E y point to the interior of the fourth and the second quadrant, which means that the local unstable manifolds W u loc (E x ) and W u loc (E y ) exist and point strictly toward E. Thus there exist points u, v in the interior of E y , E x , arbitrarily close to E y and E x such that u se T(u) se E se T(v) se v. Now, statement iii of Theorem 3 implies that E is a global attractor in [[u, v]], which completes the proof.
implies that E x is locally asymptotically stable and E y is a saddle point.In view of Lemma 2 part (b), the eigenvector which corresponds to E y points to the interior of the fourth quadrant, which means that the local unstable manifold W u loc (E y ) exists and points strictly toward E x .Thus there exists a point u in the interior of E y , E x , arbitrarily close to E y such that u se T(u).However, then this shows that the map T has a lower solution in every neighborhood of E y , which in view of Theorem 6 in [6] implies that the interior of E y , E x is a subset of the basin of attraction of E x .The result follows.
The proof when b 2 > b 1 c 2 , b 1 < b 2 c 1 is similar and will be omitted.

(d)
By Theorem 1 of [6], for each E t there exists the set W s (E t ) passing through E t and asymptotic to (0, 0), which is the graph of a continuous, non-decreasing function, which is exactly the basin of attraction of E t .The continuity of the limiting equilibrium solution as a function of initial conditions follows as in [9].(e) The proof is similar to the proof of part (c) and will be omitted.
Remark 1.In the special case c 1 = c 2 = 1 system (4) has an explicit solution.Indeed in this case we have that and so the solution of system (4) is x 0 +y 0 ), n = 0, 1, . .., which is simplified version of possible competitive scenarios, where the first two are competitive exclusion and the third is competitive coexistence.
In the special case c 1 = c 2 = 1, α 1 = α 2 > 0 system (1) is solvable, although the solution formula is complicated.Using this formula one can similarly explained competitive coexistence and competitive exclusion scenarios.(5) In this section, we present the global behavior of system (5).Denote by T(x, y) = b 1 x

Global Dynamics of System
the map associated with system (5).
Please note that a direct consequence of the above form of J T (x, y) is that T is strongly competitive.This will play a pivotal role in the global behavior of the system.System (5) always has an equilibrium solution on the y-axis, E y (0, b 2 ).Provided that α 1 < b 1 , there exists an equilibrium solution on the x-axis, E x (b 1 − α 1 , 0).Depending on the values of the parameters α 1 , b 1 , b 2 , c 1 and c 2 , there is also the possibility of either exactly one interior equilibrium solution E or an infinite number of interior equilibrium solutions E t .The interior equilibrium solution is an intersection of two equilibrium curves C 1 : > b 2 (i.e., the x-intercept of C 1 is smaller than the x-intercept of C 2 and the y-intercept of C 1 is bigger than the y-intercept of C 2 ).As in (7), these two geometrical conditions can be unified as where . Condition (12) implies that c 1 c 2 = 1 and the interior equilibrium E( x, ȳ) is given as: Please note that if ∆ 1 ∆ 2 < 0, ∆ 1 = 0 and ∆ 2 = 0, or ∆ 2 = 0 and ∆ 1 = 0 then there does not exist an interior equilibrium solution.Since α 1 ≥ b 1 implies that ∆ 1 ∆ 2 ≤ 0, we must have α 1 < b 1 in order for an interior equilibrium point to exist.Furthermore, if ∆ 1 = ∆ 2 = 0 then the two equilibrium curves C 1 and C 2 coincide and every point of the segment c The equilibrium points for system (5) are summarized in Table 3.

Condition
Equilibrium Points The local stability character of E x , E y , E and E t are presented in Lemma 4. The proof requires Proposition 1.
Proposition 1.The eigenvalues λ and µ of J T (E) are positive.
Proof.In view of ( 14) we have Please note that the equilibrium point E exists under the hypothesis In either case, we have det(J T (E)) > 0, and consequently, λ • µ > 0. Since system ( 5) is strongly competitive, by the Perron-Frobenius Theorem [15,23], the largest eigenvalue of J T (E) is positive, which completes the proof.

(a)
The equilibrium solution E x exists if . In each case, the eigenvectors associated with the eigenvalues The equilibrium solution E y always exists and it is locally asymptotically stable if b In each case, the eigenvectors associated with the eigenvalues are e 1 = (0, 1) and The interior equilibrium solution E exists if They are non-hyperbolic of the stable type and the eigenvector associated with , 1). Proof.

(a)
The Jacobian matrix associated with the map T has the form In view of ( 14), we have , which implies that the eigenvalues of the Jacobian matrix are . The corresponding eigenvectors are as stated.(b) In view of ( 14), we have which implies that the eigenvalues of the Jacobian matrix are The corresponding eigenvectors are as stated.(c) Denote the eigenvalues of J T (E) by λ 1 and λ 2 , which represent the roots of the characteristic polynomial, p(t) = t 2 − TrJ T (E)t + det J T (E).By Proposition 1, λ 1 and λ 2 are real and positive.Notice 2 and by (15), p(1) < 0. Combining this with the fact that p(0 2 and from (15) we have p(1) > 0 and p (1) > 0. Combining this with p(0) > 0, we conclude that E is locally asymptotically stable.
and λ 2 = 1.Since α 1 < b 1 we clearly have that E t are non-hyperbolic equilibrium points of the stable type.It follows by immediate checking that the eigenvector associated with , 1), which points towards the first quadrant for t ∈ [0, b 2 ).
The global behavior of system ( 5) is described by the following result.Please note that the proofs presented for Theorem 6 differ from those of Theorem 5 in order to depict an alternative approach.Theorem 6.Consider system (5).

(a)
If α 1 ≥ b 1 then E y is the unique equilibrium solution of system (5) and it is locally asymptotically stable.Every solution in the first quadrant which starts off of the x-axis converges to E y and every solution which starts on the positive x-axis converges to the singular point (0, 0).

(b)
For  ) then there is an infinite family of equilibrium solutions E t for which there exists the global stable manifold W s (E t ) for all t ∈ [0, b 2 ], which is the graph of a continuous, non-decreasing function asymptotic to (0, 0) and is exactly the basin of attraction of E t .
The limiting equilibrium varies continuously with the initial condition.(f) If E x (resp.E y ) is non-hyperbolic and E y (resp.E x ) is locally asymptotically stable then E y (resp.E x ) attracts the first quadrant of initial conditions except the positive part of the x-axis (resp.y-axis), which is attracted by E x (resp.E y ).If E x (resp.E y ) is non-hyperbolic and E y (resp.E x ) is a saddle point then E x (resp.E y ) attracts the first quadrant of initial conditions except the positive part of the y-axis (resp. x-axis), which is attracted by E y (resp.E x ).
See Figure 3 for graphical interpretation. Proof.
To treat the dynamics in the interior of R 2 + , consider R a := {(x, y) : x, y ≥ 0 and c 2 x + y ≤ b 2 }.By Theorem 2 of [4], R a is invariant.The region R a also attracts the interior of R 2 + .To verify this, suppose that (x 0 , y 0 ) ∈ R a with x 0 , y 0 > 0. In this case c 2 x 0 + y 0 > b 2 and It follows that there exists N > 0 such that for all n ≥ N, (x n , y n ) ∈ R a and thus R a is attracting.
To conclude the proof, suppose (x 0 , y 0 ) ∈ R a with x 0 , y 0 > 0. In this case x 1 < x 0 , y 1 ≥ y 0 and as a consequence of the invariance of R a , {x n } is a decreasing sequence while {y n } is a non-decreasing sequence.Therefore, lim n→∞ Tn (x 0 , y 0 ) = E y .The above arguments prove that the basins of attraction for E y and the singular point (0, 0) are B(E y ) = [0, ∞) × (0, ∞) and B(0 . Lemma 3(c) and (d) guarantee that for all initial conditions on the positive y-axis, T(x 0 , y 0 ) = E y and for all initial conditions on the positive x-axis, lim n→∞ Tn (x 0 , y 0 ) = E x .To treat the interior of R 2 + , consider R b := {(x, y) : x, y ≥ 0, c 2 x + y ≥ b 2 and α 1 + x + c 1 y ≤ b 1 } shown in Figure 4.
Please note that R b is an invariant region by Theorem 2 of [4].Consider (x 0 , y 0 ) ∈ R b with x 0 , y 0 > 0 and notice As a consequence of the invariance of R b , {x n } is a non-decreasing sequence and {y n } is a non-increasing sequence.Therefore, using basic properties of sequences and the fact that T is strongly competitive, lim n→∞ Tn (x 0 , y 0 ) = E x .Finally, suppose (x 0 , y 0 ) ∈ R b with x 0 , y 0 > 0.
By Lemma 3(a), (x n , Therefore, lim n→∞ Tn (x 0 , y 0 ) = E x .We have arrived at the desired result that the basins of attraction for E x and E y are B(E The proof for the case when b ) is similar and will be omitted.
. As in part (b), Lemma 3(c) and (d) guarantee that the positive part of the y-axis is a subset of B(E y ) and the positive part of the x-axis is a subset of B(E x ).To treat the interior of R 2 + , consider the region R c shown in Figure 5. Please note that R c is invariant by Theorem 2 of [4].Provided that (x 0 , y 0 ) ∈ R c with x 0 , y 0 > 0, monotonicity properties (similar to part (a) and (b)) along with Lemma 3(b) can be used to prove that lim n→∞ Tn (x 0 , y 0 ) = E. Suppose (x 0 , y 0 ) ∈ R c with x 0 , y 0 > 0. By Lemma 3(a) we know that (x n , y ] for all n ≥ 1.Moreover, since (x 0 , y 0 ) se (x 0 , 0) then Tn (x 0 , y 0 ) se Tn (x 0 , 0) for all n ≥ 1.Consequently, there must exist an N ≥ 1 such that Therefore, lim n→∞ Tn (x 0 , y 0 ) = E for all (x 0 , y 0 ) ∈ R c .We have reached the desired result that the basins of attraction for E, E x and E y are B . In light of Lemma 4(c), Theorems 1 and 5 of [6] guarantees that there exist the global stable and unstable manifolds for E, W s (E) and W u (E) respectively, with the above mentioned properties.An immediate checking shows that E y se E se E x and that the interior of the ordered interval E y , E is a subset of B(E y ), while the interior of the ordered interval E, E x is a subset of B(E x ).Now, take any point (x 0 , y 0 ) ∈ R 2 + such that (x 0 , y 0 ) ≺ se W s (E) (i.e., above W s (E)).Then (0, y 0 ) se (x 0 , y 0 ) ≺ se (x W s (E) , y 0 ), where (x W s (E) , y 0 ) ∈ W s (E).By Lemma 3(c) and the monotonicity of T, for n ≥ 1, E y = Tn ((0, y 0 )) se Tn ((x 0 , y 0 )) se Tn ((x W s (E) , y 0 )).
Since lim n→∞ Tn ((x W s (E) , y 0 )) = E, (18) implies that Tn ((x 0 , y 0 )) enters the ordered interval E y , E and so converges to E y .In a similar way, one can show that the ordered interval E, E x attracts all points below W s (E), and so all such points converge to E x .(e) By Theorem 1 of [6], for each E t there exists the set W s (E t ) passing through E t and asymptotic to (0, 0), which is the graph of a continuous, non-decreasing function, which is exactly the basin of attraction of E t .The continuity of the limiting equilibrium solution as a function of initial conditions follows as in [9].(f) The proof is similar to the proof of part (b) and will be omitted here.
Based on Figures 2 and 3, the global dynamics of systems ( 4) and ( 5) are similar.However, the techniques of the proofs are different since the determinant of the map corresponding to system (4) is identically zero in the first quadrant while the determinant of the map corresponding to ( 5) is positive and the map satisfies (O+) condition.This condition greatly simplifies the proof for system (5).The qualitative difference between system (1) and systems (4) and ( 5) is in the case when b 1 ≤ 1, b 2 ≤ 1 which is possible because of the fact that α 1 , α 2 > 0. In this case, E 0 (0, 0) is a globally asymptotically stable equilibrium for system (1), while the basin of attraction of the singular point E 0 (0, 0) is an empty set for system (4), and the basin of attraction of the singular point E 0 (0, 0) is the non-negative part of the x-axis for system (5).The major difference between systems (4) and ( 5) is the case α 1 ≥ b 1 which is possible only for system (5) when the singular point (0, 0) has non-empty basin of attraction.Thus the presence of α 1 will introduce new dynamic scenario while the presence of both α 1 , α 2 will introduce one additional dynamic scenario.As we mentioned in Remark 1 systems (1) and ( 4) are solvable for the special values of parameters while there is not a formula for the exact solution of system (5) at this time for any choice of parameters.

Conclusions
The presented results will serve as motivation for deriving the general results for monotone (competitive and cooperative) systems in the plane such as in [24], which will lead to some standard global dynamic scenarios results.The global dynamic scenarios results require detailed local stability analysis which can be fairly complicated and tedious as soon as at least one quadratic term is present, see for instance [25].The presence of quadratic terms may cause the emergence of new dynamic scenarios such as the Allee effect, periodic solutions, Naimark-Sacker bifurcation etc.Furthermore, we are interested in extending the presented results to higher dimensions to so called multispecies Leslie-Gower model which will consist of k ≥ 3 equation of type (1).Some results in this direction were obtained in [11,12].We would like to use the full strength of monotone systems theory which at this moment is not available in dimensions higher than 2. The k ≥ 3 dimensional version of System (4) will be an example of linear fractional homogeneous system which dynamics can be reduced to dynamics of k − 1 linear fractional homogeneous system, which means that the dynamics 3-dimensional version of System (4) can be reduced to dynamics of this system.In some special cases analogue to the ones in Remark 1 one can find the explicit solution of such system and so obtain the dynamics in such case.

Theorem 4 .
For a nonempty set R ⊂ R n and a partial order on R n , let T : R → R be an order preserving map, and let a, b ∈ R be such that a ≺ b and a, b ⊂ R. If a T(a) and T(b) b, then a, b is invariant and i.
which implies (b).Parts (a) and (c) follow by immediate checking.Part (d) follows by dividing equations of system (4).Part (e) follows from (d) since b 2 −b 1 c 2 b 1 −b 2 c 1 is exactly an equilibrium of the equation in (d).

1 ) and e 2 = ( b 1 c 1 c 2 b 2 ,
−1).(c) The interior equilibrium solution E is a saddle point when b 2 < b 1 c 2 and b 1 < b 2 c 1 and is locally asymptotically stable when b 2 > b 1 c 2 and b 1 > b 2 c 1 .(d) The interior equilibrium solutions E t are non-hyperbolic of the stable type and the eigenvector which corresponds to λ 1 = 0 is given as e 1 = (1, y).

Theorem 5 .
) Consequently, if b 2 < b 1 c 2 and b 1 < b 2 c 1 then p(1) < 0 and p (1) > 0 and if b 2 > b 1 c 2 and b 1 > b 2 c 1 then p(1) > 0 and p (1) > 0. It follows that E is a saddle point when b 2 < b 1 c 2 and b 1 < b 2 c 1 and is locally asymptotically stable when b 2 > b 1 c 2 and b 1 > b 2 c 1 .(d)In this case, the eigenvalues of the Jacobian matrix evaluated at the equilibrium E t are λ 1 = 0, λ 2 = 1.The eigenvector that corresponds to λ 1 = 0 is e 1 = (1, y), where y > 0 satisfies x + c 1 y = b 1 and points towards the first quadrant.Now, global behavior of system (4) is described by the following result: Consider system (4).

2 < b 1 c 2 ,
b 1 < b 2 c 1 , then the equilibrium solutions E x , E y are locally asymptotically stable and the interior equilibrium E is a saddle point.The separatrix S , which is a graph of a continuous, non-decreasing curve, is the basin of attraction of E and the region below (resp.above) S is the basin of attraction of E x (resp.E y ).(b)If b 2 > b 1 c 2 , b 1 > b 2 c 1, then the equilibrium solutions E x , E y are saddle points and the interior equilibrium E is locally asymptotically stable.Every solution in the first quadrant which starts off the coordinate axes converges to E. Every solution which starts on the positive part of the x-axis (resp.y-axis) is attracted by E x (resp.E y ). b

c 1 (
i.e., the x-intercept of C 2 is smaller than the x-intercept of C 1 and the y-intercept of C 2 is bigger than the y-intercept of C 1 ) or b 1 − α 1 < b 2 c 2 and b 1 −α 1 c 1 and b 2 < c 2 (b 1 − α 1 ) then system (5) has equilibrium solutions E x , E y and E. The equilibrium solutions E x and E y are locally asymptotically stable and the interior equilibrium E is a saddle point.There exists the global stable manifold W s (E) and the global unstable manifold W u (E), where W s (E) is the graph of a continuous, non-decreasing function and W u (E) is the graph of a continuous, non-increasing function which connects all three equilibrium solutions.The region in the first quadrant above (resp.below) the curve W s (E) is the basin of attraction of E y (resp.E x ) and the curve W s (E) \ {(0, 0)} is the basin of attraction of E. (e) If α 1 < b 1 , b 1 = α 1 + b 2 c 1 , and b 2 = c 2 (b 1 − α 1
a saddle No interior fixed point exists.
then there is an infinite family of equilibrium solutions E t , t ∈ [0, b 2 ] for which there exists the global stable manifold W s (E t ), which is the graph of a continuous, non-decreasing function asymptotic to (0, 0) and is exactly the basin of attraction of E t .The limiting equilibrium varies continuously with the initial condition.(e)IfE x (resp.E y ) is non-hyperbolic and E y (resp.E x ) is locally asymptotically stable then E y (resp.E x ) attracts the first quadrant of initial conditions except the positive part of x-axis (resp.y-axis) which is attracted by E x (resp.E y ).If E x (resp.E y ) is non-hyperbolic and E y (resp.E x ) is a saddle point then E x (resp.E y ) attracts the first quadrant of initial conditions except the positive part of y-axis (resp.x-axis) which is attracted by E y (resp.E x ).
5) has equilibrium solutions E x and E y where E x (resp.E y ) is locally asymptotically stable and E y (resp.E x ) is a saddle point.The basin of attraction of E x (resp.E y ) is the first quadrant of initial conditions without the positive part of the y-axis (resp.x-axis),which is attracted by E y (resp.E x ).(c)Ifα 1 < b 1 , b 1 > α 1 + b 2 c 1 , and b 2 > c 2 (b 1 − α 1 )then system (5) has equilibrium solutions E x , E y and E. The equilibrium solutions E x and E y are saddle points and E is locally asymptotically stable.Every solution in the first quadrant which starts off of the coordinate axes converges to E and every solution which starts on the positive x-axis (resp.y-axis) converges to E x (resp.E y ).