Next Article in Journal
A Refinement of Schwarz–Pick Lemma for Higher Derivatives
Next Article in Special Issue
Stability and Bifurcation Analysis on a Predator–Prey System with the Weak Allee Effect
Previous Article in Journal
A Study of an EOQ Model under Lock Fuzzy Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Dynamics of Leslie-Gower Competitive Systems in the Plane

by
Mustafa R. S. Kulenović
* and
David T. McArdle
Department of Mathematics, University of Rhode Island, Kingston, RI 02881-0816, USA
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(1), 76; https://doi.org/10.3390/math7010076
Submission received: 5 December 2018 / Revised: 4 January 2019 / Accepted: 9 January 2019 / Published: 12 January 2019

Abstract

:
Global dynamic results are obtained for families of competitive systems of difference equations of the form x n + 1 = b 1 x n α 1 + x n + c 1 y n , y n + 1 = b 2 y n α 2 + c 2 x n + y n n = 0 , 1 , , where the parameters b 1 , b 2 are positive numbers, and α 1 , α 2 , c 1 , and c 2 and the initial conditions x 0 and y 0 are arbitrary non-negative numbers, when one or both of α i , i = 1 , 2 equalls 0. We assume that the denominators of both equations are always positive. We will show that the presence of more parameters will create more dynamic scenarios.

1. Introduction

The system of difference equations,
x n + 1 = b 1 x n α 1 + x n + c 1 y n , y n + 1 = b 2 y n α 2 + c 2 x n + y n n = 0 , 1 , ,
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
E 0 ( 0 , 0 ) , E x ( b 1 1 , 0 ) , E y ( 0 , b 2 1 )
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
E b 2 1 c 1 c 2 1 c 1 b 1 1 b 2 1 , b 1 1 c 1 c 2 1 c 2 b 2 1 b 1 1 .
Important subsets of parameter space are described in Table 1, together with corresponding behavior of equilibrium solutions established in [1].
Theorem 4.1 in [5] applies when parameters vary from Case 1 to Case 4 of Table 1. Set
A = { a R + 4 : a = ( b 1 , b 2 , c 1 , c 2 ) and c 1 ( b 2 1 ) > b 1 1 > 0 }
and define T a to be the map of Equation (1) restricted to R = [ 0 , ) × ( 0 , ) , that is,
T a ( x , y ) = b 1 x 1 + x + c 1 y , b 2 y 1 + c 2 x + y .
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 Γ ( b 1 , b 2 , c 1 , c 2 ) = c 2 ( b 1 1 ) b 2 + 1 = 0 .
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.
Theorem 1.
Consider system (1).
(i) 
Suppose that c 1 ( b 2 1 ) > b 1 1 > 0 . If c 2 ( b 1 1 ) > b 2 1 , then E y is globally asymptotically stable on [ 0 , ) × ( 0 , ) , and E x attracts all points on the open semiaxis ( 0 , ) × { 0 } . If c 2 ( b 1 1 ) < b 2 1 , then the stable manifold W s ( E ) in [ 0 , ) × [ 0 , ) 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 1 ( b 2 1 ) < b 1 1 > 0 . If c 2 ( b 1 1 ) > b 2 1 , then E x is globally asymptotically stable on [ 0 , ) × ( 0 , ) , and E y attracts all points on the open semiaxis { 0 } × ( 0 , ) . If c 2 ( b 1 1 ) < b 2 1 , then E is globally asymptotically stable on ( 0 , ) × ( 0 , ) , E y attracts all points on the open semiaxis { 0 } × ( 0 , ) , and E x attracts all points on the open semiaxis ( 0 , ) × { 0 } .
See Figure 1 for graphical interpretation.
The non-hyperbolic case when
c 1 ( b 2 1 ) = b 1 1 and c 2 ( b 1 1 ) = b 2 1
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
E t : = ( b 1 1 ) ( 1 t ) , ( b 2 1 ) t , 0 t 1 .
The eigenvalues of the Jacobian of T a at E t are
λ 1 = 1 and λ 2 = ( 1 t ) 1 b 1 + t 1 b 2 , 0 t 1 ,
and corresponding eigenvectors are
e 1 = 1 b 1 1 b 2 , 1 and e 2 = b 2 ( 1 b 1 ) 2 ( 1 t ) , b 1 ( 1 b 2 ) 2 t , 0 t 1 .
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 a n ( x , y ) , then T * is continuous. These observations, together with an application of Theorem 1 in [6] lead to the following result.
Theorem 2.
Assume (3) holds. Then,
(i) 
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 ( x ¯ , y ¯ ) s 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
x n + 1 = b 1 x n x n + c 1 y n , y n + 1 = b 2 y n c 2 x n + y n n = 0 , 1 ,
and
x n + 1 = b 1 x n α 1 + x n + c 1 y n , y n + 1 = b 2 y n c 2 x n + y n n = 0 , 1 , ,
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
x n + 1 = α + β x n + γ x n 1 A + B x n + C x n 1 . n = 0 , 1 , ,
with all non-negative parameters and initial conditions such that A + B x n + C x 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.

2. 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
x n + 1 = f ( x n , y n ) y n + 1 = g ( x n , y n ) , n = 0 , 1 , 2 , , ( x 1 , x 0 ) R
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) have been studied by many authors such as Clark, Hess, Hirsch and Smith, Kulenović, Merino, Nurkanović, Leonard and May, Selgrade, Smale, Smith [4,6,7,14,15,16,17,18,19] and others. In [1,2] the authors gave an interesting applications of this theory to some basic models in population dynamics. Many open problems and conjectures about the dynamics of competitive systems of linear fractional difference equations are gien in [8].
Denote with s e 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 n e 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 ) s e T ( v ) whenever u s e 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 ) n e T ( y ) implies x n e y , and T is said to satisfy condition ( O ) if for every x, y in S , T ( x ) n e T ( y ) implies y n e 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 2 n } 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.
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. 
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).

3. Global Dynamics of System (4)

First we give some basic results about the global behavior of system (4). Denote by
T ( x , y ) = b 1 x x + c 1 y , b 2 y c 2 x + y ,
the map associated with system (4). System (4) is homogeneous and was partially investigated in [22].
Lemma 1.
(a) 
Every solution of system (4) satisfies x n b 1 , y n b 2 , n 1 .
(b) 
det J T ( x , y ) = 0 for every ( x , y ) , where J T denotes the Jacobian matrix of the map T.
(c) 
T ( x 0 , 0 ) = E x ( b 1 , 0 ) , T ( 0 , y 0 ) = E y ( 0 , b 2 ) for every x 0 > 0 , y 0 > 0 .
(d) 
Every solution of system (4) satisfies the difference equation
r n + 1 = b 2 b 1 r n 1 + c 1 r n c 2 + r n , n 0 ,
where r n = y n / x n .
(e) 
If
( b 2 b 1 c 2 ) ( b 1 b 2 c 1 ) > 0 ,
then the map T has an invariant line
S : y = b 2 b 1 c 2 b 1 b 2 c 1 x .
Proof. 
The Jacobian matrix J T of the map T has the form
J T ( x , y ) = b 1 c 1 y ( x + c 1 y ) 2 b 1 c 1 x ( x + c 1 y ) 2 b 2 c 2 y ( c 2 x + y ) 2 b 2 c 2 x ( c 2 x + y ) 2 ,
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). □
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) C 1 : x + c 1 y = b 1 , C 2 : c 2 x + y = b 2 it will exist if either b 2 c 2 < b 1 , b 1 c 1 < b 2 (x-intercept of C 2 smaller than x-intercept of C 1 and y-intercept of C 2 bigger than y-intercept of C 1 ) or b 2 c 2 > b 1 , b 1 c 1 > b 2 (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 ¯ , y ¯ ) is given as:
x ¯ = b 1 b 2 c 1 1 c 1 c 2 , y ¯ = b 2 b 1 c 2 1 c 1 c 2 .
Notice that b 2 < b 1 c 2 , b 1 < b 2 c 1 implies c 1 c 2 > 1 and b 2 > b 1 c 2 , b 1 > b 2 c 1 implies c 1 c 2 < 1 . If either ( b 2 b 1 c 2 ) ( b 1 b 2 c 1 ) < 0 , b 2 = b 1 c 2 and b 1 b 2 c 1 , or b 1 = b 2 c 1 and b 2 b 1 c 2 then there are no interior equilibrium points. Furthermore, if b 2 = b 1 c 2 and b 1 = b 2 c 1 , the two equilibrium curves C 1 , C 2 coincide and every point on the segment x + c 1 y = b 1 , x , y 0 is an equilibrium solution E t ( b 1 c 1 t , t ) , t [ 0 , b 2 ] . See Table 2 summarizing the equilibrium points of system (4).
The following result describes the local stability character of all equilibrium solutions.
Lemma 2.
Consider system (4).
(a) 
The equilibrium solution E x is locally asymptotically stable if b 2 < b 1 c 2 , non-hyperbolic of stable type if b 2 = b 1 c 2 and a saddle point if b 2 > b 1 c 2 . In each case, the eigenvectors associated with the eigenvalues λ 1 = 0 and λ 2 = b 2 b 1 c 2 are e 1 = ( 1 , 0 ) and e 2 = ( b 1 c 1 c 2 b 2 , 1 ) .
(b) 
The equilibrium solution E y is locally asymptotically stable if b 1 < b 2 c 1 , non-hyperbolic of stable type if b 1 = b 2 c 1 and a saddle point if b 1 > b 2 c 1 . In each case, the eigenvectors associated with the eigenvalues λ 1 = 0 and λ 2 = b 1 b 2 c 1 are e 1 = ( 0 , 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 ) .
Proof. 
(a)
In view of (9), we have
J T ( E x ) = 0 c 1 0 b 2 b 1 c 2 ,
which implies that the eigenvalues of the Jacobian matrix are λ 1 = 0 , λ 2 = b 2 b 1 c 2 . The corresponding eigenvectors are as stated.
(b)
In view of (9), we have
J T ( E y ) = b 1 b 2 c 1 0 c 2 0 ,
which implies that the eigenvalues of the Jacobian matrix are λ 1 = 0 , λ 2 = b 1 b 2 c 1 . 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
p ( 1 ) = ( b 1 b 2 c 1 ) ( b 2 b 1 c 2 ) b 1 b 2 ( 1 c 1 c 2 ) and p ( 1 ) = b 2 ( b 1 b 2 c 1 ) + b 1 ( b 2 b 1 c 2 ) b 1 b 2 ( 1 c 1 c 2 ) .
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:
Theorem 5.
Consider system (4).
(a) 
If b 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 ).
(c) 
If b 2 < b 1 c 2 , b 1 > b 2 c 1 (resp. b 2 > b 1 c 2 , b 1 < b 2 c 1 ), 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 x ).
(d) 
If b 1 = b 2 c 1 and b 2 = b 1 c 2 , 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) 
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 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 ).
See Figure 2 for graphical interpretation.
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
b 1 c 2 ( c 1 c 2 1 ) x 2 + ( b 1 + b 2 c 1 ) ( b 1 c 2 b 2 ) x + b 1 b 2 ( b 1 + b 2 c 1 ) = 0 .
Please note that
D = ( b 1 + b 2 c 1 ) ( b 2 + b 1 c 2 ) ( b 1 c 2 ( b 1 b 2 c 1 ) + b 2 c 1 ( b 2 b 1 c 2 ) b 1 b 2 ( c 1 c 2 1 ) ) < 0 ,
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 S will satisfy T n ( ( x 0 , y 0 ) ) [ [ E , E x ] ] 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 l o c u ( E x ) and W l o c u ( 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 s e T ( u ) s e E s e T ( v ) s e v . Now, statement iii of Theorem 3 implies that E is a global attractor in [ [ u , v ] ] , which completes the proof.
(c)
Assume that b 2 < b 1 c 2 , b 1 > b 2 c 1 which 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 l o c u ( 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 s e 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 y n + 1 x n + 1 = B y n x n , where B = b 2 b 1 and so the solution of system (4) is
x n = b 1 x 0 x 0 + y 0 B n 1 , y n = b 1 y 0 B n x 0 + y 0 B n 1 , n = 0 , 1 , .
Thus B < 1 implies ( x n , y n ) ( b 1 , 0 ) , n and B > 1 implies ( x n , y n ) ( 0 , b 2 ) , n while B = 1 implies ( x n , y n ) = ( b 1 x 0 x 0 + y 0 , b 1 y 0 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.

4. Global Dynamics of System (5)

In this section, we present the global behavior of system (5). Denote by
T ˜ ( x , y ) = b 1 x α 1 + x + c 1 y , b 2 y c 2 x + y ,
the map associated with system (5).
Lemma 3.
(a) 
Every solution of system (5) satisfies x n < b 1 , y n b 2 , n 1 ;
(b) 
T ˜ ( x , y ) satisfies ( O + ) condition on ( 0 , ) 2 and so T has no period-two points;
(c) 
For every y 0 > 0 , T ˜ ( 0 , y 0 ) = ( 0 , b 2 ) ;
(d) 
For every x 0 > 0 , T ˜ n ( x 0 , 0 ) = b 1 α 1 1 ( α 1 b 1 ) n , 0 .
Proof. 
Part (a) and (c) follow by immediate checking and part (d) follows by solving the resulting Beverton-Holt one-dimensional equation. Notice
T ˜ x 1 y 1 n e T ˜ x 2 y 2 b 1 x 1 α 1 + x 1 + c 1 y 1 b 2 y 1 c 2 x 1 + y 1 n e b 1 x 2 α 1 + x 2 + c 1 y 2 b 2 y 2 c 2 x 2 + y 2 .
This is equivalent to
α 1 ( x 1 x 2 ) c 1 ( x 2 y 1 x 1 y 2 ) , x 2 y 1 x 1 y 2 0 .
This implies that x 1 x 2 and x 2 y 1 x 1 y 2 and thus x 1 x 2 , y 1 y 2 , that is ( x 1 , y 1 ) n e ( x 2 , y 2 ) . □
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 : α 1 + x + c 1 y = b 1 and C 2 : c 2 x + y = b 2 . This solution will exist if either b 2 c 2 < b 1 α 1 and b 2 > b 1 α 1 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 > 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
Δ 1 Δ 2 > 0 ,
where Δ 1 = ( b 2 c 2 ( b 1 α 1 ) ) and Δ 2 = ( b 1 α 1 b 2 c 1 ) . Condition (12) implies that c 1 c 2 1 and the interior equilibrium E ( x ¯ , y ¯ ) is given as:
x ¯ = b 1 α 1 b 2 c 1 1 c 1 c 2 , y ¯ = b 2 c 2 ( b 1 α 1 ) 1 c 1 c 2 .
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 2 x + y = b 2 , x , y 0 is an equilibrium solution E t ( b 1 α 1 c 1 t , t ) , t [ 0 , b 2 ] . The equilibrium points for system (5) are summarized in Table 3.
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
J T ˜ ( E ) = b 2 c 1 b 1 c 1 c 2 + α 1 b 1 ( 1 c 1 c 2 ) c 1 ( b 1 b 2 c 1 α 1 ) b 1 ( 1 c 1 c 2 ) c 2 ( b 2 c 2 ( b 1 α 1 ) ) b 2 ( 1 c 1 c 2 ) c 2 ( b 1 b 2 c 1 α 1 ) b 2 ( 1 c 1 c 2 ) ,
which implies
det ( J T ˜ ( E ) ) = c 2 α 1 ( b 1 α 1 b 2 c 1 ) b 1 b 2 ( 1 c 1 c 2 ) .
Please note that the equilibrium point E exists under the hypothesis Δ 1 Δ 2 > 0 , which means that either b 2 c 2 < b 1 α 1 and b 2 > b 1 α 1 c 1 or b 1 α 1 < b 2 c 2 and b 1 α 1 c 1 > b 2 . 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. □
Lemma 4.
Consider system (5).
(a) 
The equilibrium solution E x exists if α 1 < b 1 . It is locally asymptotically stable if b 2 < c 2 ( b 1 α 1 ) , non-hyperbolic of stable type if b 2 = c 2 ( b 1 α 1 ) and a saddle point if b 2 > c 2 ( b 1 α 1 ) . In each case, the eigenvectors associated with the eigenvalues λ 1 = α 1 b 1 and λ 2 = b 2 c 2 ( b 1 α 1 ) are e 1 = ( 1 , 0 ) and e 2 = ( c 1 c 2 ( b 1 α 1 ) 2 c 2 α 1 ( b 1 α 1 ) b 1 b 2 , 1 ) , respectively.
(b) 
The equilibrium solution E y always exists and it is locally asymptotically stable if b 1 < α 1 + b 2 c 1 , non-hyperbolic of stable type if b 1 = α 1 + b 2 c 1 and a saddle point if b 1 > α 1 + b 2 c 1 . In each case, the eigenvectors associated with the eigenvalues λ 1 = 0 and λ 2 = b 1 b 2 c 1 + α 1 are e 1 = ( 0 , 1 ) and e 2 = ( b 1 c 2 ( α 1 + b 2 c 1 ) , 1 ) , respectively.
(c) 
The interior equilibrium solution E exists if Δ 1 Δ 2 > 0 and it is locally asymptotically stable if b 1 > α 1 + b 2 c 1 and b 2 > c 2 ( b 1 α 1 ) and a saddle point if b 1 < α 1 + b 2 c 1 and b 2 < c 2 ( b 1 α 1 ) .
(d) 
The interior equilibrium solutions E t exist if α 1 < b 1 , b 1 = α 1 + b 2 c 1 and b 2 = c 2 ( b 1 α 1 ) . They are non-hyperbolic of the stable type and the eigenvector associated with λ 1 where | λ 1 | < 1 is e 1 = ( ( b 1 α 1 ) ( b 2 t ) b 1 c 2 t , 1 ) .
Proof. 
(a)
The Jacobian matrix associated with the map T ˜ has the form
J T ˜ ( x , y ) = b 1 ( α 1 + c 1 y ) ( α 1 + x + c 1 y ) 2 b 1 c 1 x ( α 1 + x + c 1 y ) 2 b 2 c 2 y ( c 2 x + y ) 2 b 2 c 2 x ( c 2 x + y ) 2 .
In view of (14), we have
J T ˜ ( E x ) = α 1 b 1 c 1 ( b 1 α 1 ) b 1 0 b 2 c 2 ( b 1 α 1 ) ,
which implies that the eigenvalues of the Jacobian matrix are λ 1 = α 1 b 1 , λ 2 = b 2 c 2 ( b 1 α 1 ) . The corresponding eigenvectors are as stated.
(b)
In view of (14), we have
J T ˜ ( E y ) = b 1 α 1 + b 2 c 1 0 c 2 0 ,
which implies that the eigenvalues of the Jacobian matrix are λ 1 = 0 , λ 2 = b 1 α 1 + b 2 c 1 . 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 Tr J T ˜ ( E ) t + det J T ˜ ( E ) . By Proposition 1, λ 1 and λ 2 are real and positive. Notice
p ( 1 ) = ( b 1 α 1 b 2 c 1 ) ( b 2 c 2 ( b 1 α 1 ) ) b 1 b 2 ( 1 c 1 c 2 ) and p ( 1 ) = b 2 ( b 1 α 1 b 2 c 1 ) + b 1 ( b 2 c 2 ( b 1 α 1 ) ) b 1 b 2 ( 1 c 1 c 2 )
If b 1 < α 1 + b 2 c 1 and b 2 < c 2 ( b 1 α 1 ) then 1 < c 1 c 2 and by (15), p ( 1 ) < 0 . Combining this with the fact that p ( 0 ) = det J T ˜ ( E ) > 0 , it follows that E is a saddle point. If b 1 > α 1 + b 2 c 1 and b 2 > c 2 ( b 1 α 1 ) then 1 > c 1 c 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.
(d)
For t [ 0 , b 2 ) , the eigenvalues of J T ˜ ( E t ) are λ 1 = α 1 ( b 2 t ) b 1 b 2 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 is e 1 = ( ( b 1 α 1 ) ( b 2 t ) b 1 c 2 t , 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 α 1 < b 1 , if b 1 > α 1 + b 2 c 1 , b 2 < c 2 ( b 1 α 1 ) (resp. b 1 < α 1 + b 2 c 1 , b 2 > c 2 ( b 1 α 1 ) ) then system (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 ).
(d) 
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 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 ) 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. 
(a)
Let α 1 b 1 . Lemma 3(c) and (d) guarantee that for initial conditions on the positive y-axis, T ˜ ( x 0 , y 0 ) = E y and for initial conditions on the positive x-axis, lim n T ˜ n ( x 0 , y 0 ) = ( 0 , 0 ) . 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
x 1 = b 1 x 0 α 1 + x 0 + c 1 y 0 < x 0 and y 1 = b 2 y 0 c 2 x 0 + y 0 < y 0 .
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 T ˜ n ( 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 , 0 ) = [ 0 , ) × { 0 } .
(b)
Let α 1 < b 1 , b 1 > α 1 + b 2 c 1 , and b 2 < c 2 ( b 1 α 1 ) . 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 T ˜ n ( 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
x 1 = b 1 x 0 α 1 + x 0 + c 1 y 0 x 0 and y 1 = b 2 y 0 c 2 x 0 + y 0 y 0 .
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 T ˜ n ( 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 , y n ) [ 0 , b 1 ) × [ 0 , b 2 ] for all n 1 . Choose ( u , v ) R b with u , v > 0 such that ( u , v ) s e ( x 1 , y 1 ) s e ( x 1 , 0 ) . Since T ˜ is strongly competitive, notice
T ˜ n ( u , v ) s e T ˜ n ( x 1 , y 1 ) s e T ˜ n ( x 1 , 0 ) .
Therefore, lim n T ˜ n ( 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 x ) = ( 0 , ) × [ 0 , ) and B ( E y ) = { 0 } × ( 0 , ) .
The proof for the case when b 1 < α 1 + b 2 c 1 , b 2 > c 2 ( b 1 α 1 ) is similar and will be omitted.
(c)
Let α 1 < b 1 , b 1 > α 1 + b 2 c 1 , and b 2 > c 2 ( b 1 α 1 ) . 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 T ˜ n ( 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 n ) [ 0 , b 1 ) × [ 0 , b 2 ] for all n 1 . Moreover, since ( x 0 , y 0 ) s e ( x 0 , 0 ) then T ˜ n ( x 0 , y 0 ) s e T ˜ n ( x 0 , 0 ) for all n 1 . Consequently, there must exist an N 1 such that ( x N , y N ) [ 0 , b 2 c 2 ) × [ 0 , b 2 ] . Now, choose ( u , v ) , ( s , t ) R c such that ( u , v ) s e ( x N , y N ) s e ( s , t ) . Since T ˜ is strongly competitive we have
T ˜ n ( u , v ) s e T ˜ n ( x N , y N ) s e T ˜ n ( s , t ) .
Therefore, lim n T ˜ n ( 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 ( E ) = ( 0 , ) × ( 0 , ) , B ( E x ) = ( 0 , ) × { 0 } and B ( E y ) = { 0 } × ( 0 , ) .
(d)
Let α 1 < b 1 , b 1 < α 1 + b 2 c 1 , and b 2 < c 2 ( b 1 α 1 ) . 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 s e E s e 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 ) s e W s ( E ) (i.e., above W s ( E ) ). Then ( 0 , y 0 ) s e ( x 0 , y 0 ) s e ( 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 = T ˜ n ( ( 0 , y 0 ) ) s e T ˜ n ( ( x 0 , y 0 ) ) s e T ˜ n ( ( x W s ( E ) , y 0 ) ) .
Since lim n T ˜ n ( ( x W s ( E ) , y 0 ) ) = E , (18) implies that T ˜ n ( ( 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 Figure 2 and Figure 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.

5. 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.

Author Contributions

Both authors contributed equally in the writing of this article.

Funding

Mustafa R. S. Kulenović is supported in part by Maitland P. Simmons Foundation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cushing, J.M.; Levarge, S.; Chitnis, N.; Henson, S.M. Some discrete competition models and the competitive exclusion principle. J. Differ. Equ. Appl. 2004, 10, 1139–1152. [Google Scholar] [CrossRef]
  2. Cushing, J.M.; Levarge, S. Some discrete competition models and the principle of competitive exclusion. In Difference Equations and Discrete Dynamical Systems; World Scientific: Hackensack, NJ, USA, 2005; pp. 283–301. [Google Scholar]
  3. Cushing, J.M.; Henson, S.M.; Blackburn, C.C. Multiple mixed-type attractors in a competition model. J. Biol. Dyn. 2007, 1, 347–362. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Kulenović, M.R.S.; Merino, O. Competitive-Exclusion versus Competitive-Coexistence for Systems in the Plane. Discret. Contin. Dyn. Syst. Ser. B 2006, 6, 1141–1156. [Google Scholar]
  5. Kulenović, M.R.S.; Merino, O. Global Bifurcation for Competitive Systems in the Plane. Discret. Contin. Dyn. Syst. B 2009, 12, 133–149. [Google Scholar]
  6. Kulenović, M.R.S.; Merino, O. Invariant Manifolds for Competitive Discrete Systems in the Plane. Int. J. Bifurc. Chaos 2010, 20, 2471–2486. [Google Scholar] [CrossRef]
  7. Smith, H.L. Planar Competitive and Cooperative Difference Equations. J. Differ. Equ. Appl. 1998, 3, 335–357. [Google Scholar] [CrossRef]
  8. Camouzis, E.; Kulenović, M.R.S.; Ladas, G.; Merino, O. Rational Systems in the Plane-Open Problems and Conjectures. J. Differ. Equ. Appl. 2009, 15, 303–323. [Google Scholar] [CrossRef]
  9. Burgić, D.; Kalabušić, S.; Kulenović, M.R.S. Non-hyperbolic Dynamics for Competitive Systems in the Plane and Global Period-doubling Bifurcations. Adv. Dyn. Syst. Appl. 2008, 3, 229–249. [Google Scholar]
  10. Sacker, R. Global stability in a multi-species periodic Leslie-Gower model. J. Biol. Dyn. 2011, 5, 549–562. [Google Scholar] [CrossRef]
  11. Sacker, R. An invariance theorem for mappings. J. Differ. Equ. Appl. 2012, 18, 163–166. [Google Scholar] [CrossRef]
  12. Sacker, R. A note: An invariance theorem for mappings II. J. Dyn. Differ. Equ. 2012, 24, 595–599. [Google Scholar] [CrossRef]
  13. Kulenović, M.R.S.; Ladas, G. Dynamics of Second Order Rational Difference Equations with Open Problems and Conjectures; Chapman and Hall/CRC: Boca Raton, FL, USA; London, UK, 2001. [Google Scholar]
  14. Clark, D.; Kulenović, M.R.S.; Selgrade, J.F. Global Asymptotic Behavior of a Two Dimensional Difference Equation Modelling Competition. Nonlinear Anal. TMA 2003, 52, 1765–1776. [Google Scholar] [CrossRef]
  15. Hess, P. Periodic-Parabolic Boundary Value Problems and Positivity; Pitman Research Notes in Mathematics Series, 247; Longman Scientific & Technical: Harlow, UK, 1991; pp. 619–620. [Google Scholar]
  16. Kalabušić, S.; Kulenović, M.R.S.; Pilav, E. Dynamics of a two-dimensional system of rational difference equations of Leslie-Gower type. Adv. Differ. Equ. 2011, 229, 29. [Google Scholar] [CrossRef]
  17. Kulenović, M.R.S.; Nurkanović, M. Asymptotic Behavior of a Linear Fractional System of Difference Equations. J. Inequal. Appl. 2005, 2005, 127–143. [Google Scholar] [CrossRef]
  18. Leonard, W.J.; May, R. Nonlinear aspects of competition between species. SIAM J. Appl. Math. 1975, 29, 243–275. [Google Scholar]
  19. Smale, S. On the differential equations of species in competition. J. Math. Biol. 1976, 3, 5–7. [Google Scholar] [CrossRef] [PubMed]
  20. de Mottoni, P.; Schiaffino, A. Competition systems with periodic coefficients: A geometric approach. J. Math. Biol. 1981, 11, 319–335. [Google Scholar] [CrossRef]
  21. Smith, H.L. Periodic competitive differential equations and the discrete dynamics of competitive maps. J. Differ. Equ. 1986, 64, 165–194. [Google Scholar] [CrossRef] [Green Version]
  22. Huang, Y.S.; Knopf, P.M. Global convergence properties of first-order homogeneous systems of rational difference equations. J. Differ. Equ. Appl. 2012, 18, 1683–1707. [Google Scholar] [CrossRef]
  23. Hirsch, M.; Smith, H. Monotone dynamical systems. In Handbook of Differential Equations: Ordinary Differential Equations; Elsevier B. V.: Amsterdam, The Netherlands, 2005; Volume II, pp. 239–357. [Google Scholar]
  24. Bertrand, E.; Kulenović, M.R.S. Global Dynamic Scenarios for Competitive Maps in the Plane. Adv. Differ. Equ. 2018. [Google Scholar] [CrossRef]
  25. Brett, A.; Kulenović, M.R.S. Two Species Competitive Model with the Allee Effect. Adv. Differ. Equ. 2014, 2014, 307. [Google Scholar] [CrossRef]
Figure 1. Global dynamics of System (1).
Figure 1. Global dynamics of System (1).
Mathematics 07 00076 g001
Figure 2. Global dynamics of System (4).
Figure 2. Global dynamics of System (4).
Mathematics 07 00076 g002
Figure 3. Global dynamics of System (5).
Figure 3. Global dynamics of System (5).
Mathematics 07 00076 g003
Figure 4. R b : = { ( x , y ) : x , y 0 , c 2 x + y b 2 and α 1 + x + c 1 y b 1 } .
Figure 4. R b : = { ( x , y ) : x , y 0 , c 2 x + y b 2 and α 1 + x + c 1 y b 1 } .
Mathematics 07 00076 g004
Figure 5. R c : = { ( x , y ) : x , y 0 and ( x + c 1 y + α 1 b 1 ) ( c 2 x + y b 2 ) 0 } .
Figure 5. R c : = { ( x , y ) : x , y 0 and ( x + c 1 y + α 1 b 1 ) ( c 2 x + y b 2 ) 0 } .
Mathematics 07 00076 g005
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.
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.
Condition c 1 ( b 2 1 ) < b 1 1 c 1 ( b 2 1 ) > b 1 1
c 2 ( b 1 1 ) < b 2 1 Case 1.
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.
c 2 ( b 1 1 ) > b 2 1 Case 3.
E 0 is a repeller
E x is attractor on ( 0 , ) × [ 0 , )
E y is a saddle
No interior fixed point exists.
Case 4.
E 0 is a repeller
E x is a local attractor
E y is a local attractor
E is an interior saddle.
Table 2. The equilibrium points of system (4).
Table 2. The equilibrium points of system (4).
ConditionEquilibrium Points
( b 2 b 1 c 2 ) ( b 1 b 2 c 1 ) < 0 ,
b 2 = b 1 c 2   and   b 1 b 2 c 1 ,
  or   b 2 b 1 c 2 and b 1 = b 2 c 1
E x , E y
( b 2 b 1 c 2 ) ( b 1 b 2 c 1 ) > 0 E x , E y , E
b 2 = b 1 c 2   and   b 1 = b 2 c 1 E x , E y , E t
Table 3. The equilibrium points for system (5).
Table 3. The equilibrium points for system (5).
ConditionEquilibrium Points
α 1 b 1 E y
α 1 < b 1 Δ 1 Δ 2 < 0
Δ 1 = 0 , Δ 2 0 ,
Δ 1 0 , Δ 2 = 0
E x , E y
Δ 1 Δ 2 > 0 E x , E y , E
Δ 1 = Δ 2 = 0 E x , E y , E t

Share and Cite

MDPI and ACS Style

Kulenović, M.R.S.; McArdle, D.T. Global Dynamics of Leslie-Gower Competitive Systems in the Plane. Mathematics 2019, 7, 76. https://doi.org/10.3390/math7010076

AMA Style

Kulenović MRS, McArdle DT. Global Dynamics of Leslie-Gower Competitive Systems in the Plane. Mathematics. 2019; 7(1):76. https://doi.org/10.3390/math7010076

Chicago/Turabian Style

Kulenović, Mustafa R. S., and David T. McArdle. 2019. "Global Dynamics of Leslie-Gower Competitive Systems in the Plane" Mathematics 7, no. 1: 76. https://doi.org/10.3390/math7010076

APA Style

Kulenović, M. R. S., & McArdle, D. T. (2019). Global Dynamics of Leslie-Gower Competitive Systems in the Plane. Mathematics, 7(1), 76. https://doi.org/10.3390/math7010076

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop