Dynamics of a two-dimensional competitive system of rational difference equations with quadratic terms

We investigate global dynamics of the following systems of difference equations: {xn+1=b1xn2A1+yn2,yn+1=a2+c2yn2xn2,n=0,1,2,…, where the parameters b1, a2, A1, c2 are positive numbers and the initial condition y0 is an arbitrary nonnegative number and x0 is a positive number. We show that this system has rich dynamics which depends on the part of a parametric space. We find precisely the basins of attraction of all attractors including the points at ∞. MSC:39A10, 39A30, 37E99, 37D10.


Introduction
In this paper we study the global dynamics of the following rational system of difference equations: where the parameters b  , a  , A  , c  are positive numbers and the initial condition y  is an arbitrary nonnegative number and x  is a positive number.
The related system of difference equations where the parameters A  , β  , α  and γ  are positive numbers and the initial conditions x  > , y  ≥ , was considered in [], where it was shown that this system has simple dynamics.Precisely, it was shown that system () has no equilibrium points if β  ≤ A  and that it has a unique equilibrium point if β  > A  , in which case this equilibrium point is a saddle point.Furthermore, the following result describes the global dynamics of system ().http://www.advancesindifferenceequations.com/content/2014/1/301Theorem  Consider system ().
() Assume that β  > A  and γ  A  = α  .Then there exist a set C ⊂ R which is invariant and a subset of the basin of attraction of E. The set C is a graph of a strictly increasing continuous function of the first variable on an interval (and so is a manifold) and separates R into two connected and invariant components, namely W -:= {x ∈ R\C : ∃y ∈ C with x se y} and W + := {x ∈ R\C : ∃y ∈ C with y se x}, which satisfy: lim n→∞ (x n , y n ) = (, ∞) for every (x  , y  ) ∈ W -, and lim n→∞ (x n , y n ) = (∞, ) for every (x  , y  ) ∈ W + .
() Assume that γ  A  = α  .Then system () can be decoupled as follows: and every solution of this system (depending of the choice of the initial condition (x  , y  )) is either bounded and converges to an equilibrium point, or increases monotonically to infinity.
Thus every solution of system () either converges to the unique equilibrium point or is asymptotic to one of the points at infinity, precisely to either (, ∞) or to (∞, ).In all cases, either solution is eventually monotonic or the subsequences of even indexed and odd indexed terms are eventually monotonic.Introduction of quadratic terms into the system will substantially change the dynamics by introducing new equilibrium points (up to three) with different local character and minimal period-two solutions (up to ).Again, most of the solutions of system () will be asymptotic to (∞, ) or (, ∞), but the separatrix between the two basins of attraction may consist of several global stable manifolds of either saddle point equilibrium points or non-hyperbolic equilibrium points or minimal periodtwo solutions.In one case, when there exists a unique non-hyperbolic equilibrium point, it is possible that this point will have a basin of attraction of positive Lebesgue measure.
System () is a competitive system, and our results are based on recent results about competitive systems in the plane, see [, ].System () can be used as a mathematical model for competition in population dynamics.The first systematic study of a specific competitive system with quadratic terms was performed in [] where system of the form where the parameters a, b >  and the initial conditions x  , y  ≥ , was considered.It was shown that the dynamics of system () is very similar to the dynamics of the corresponding linear fractional system , n = , , . . ., with the same conditions on parameters and initial conditions.Both systems have nine parametric regions with different dynamical behavior.
As noted, the introduction of quadratic terms in system () dramatically changes the dynamics.The techniques used to study system () were straightforward calculations, while the techniques which will be used to study system () are a combination of techniques for studying real algebraic curves and implicit function theorem as neither equilibrium points nor period-two solutions are explicitly computable.Some of our calculations are performed by using Mathematica and outputs are included in the Appendix.
The paper is organized as follows.Section  contains some necessary results on competitive systems in the plane.Section  provides some basic facts about the equilibrium points and injectivity of the map associated with system ().Section  contains local stability analysis of both equilibrium solutions and minimal period-two solutions.Section  gives global dynamics in different cases.

Preliminaries
A first-order system of difference equations x n+ = f (x n , y n ), y n+ = g(x n , y n ), n = , , , . . ., () where S ⊂ R  , (f , g) : S → S, f , g are continuous functions, 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 nondecreasing in y.If both f and g are non-decreasing in x and y, system () is cooperative.
Competitive and cooperative maps are defined similarly.Strongly competitive systems of difference equations or strongly competitive maps are those for which the functions f and g are coordinate-wise strictly monotone.Competitive and cooperative systems have been investigated by many authors, see [, , -].Special attention to discrete competitive and cooperative systems in the plane was given in [, , , , , , -].One of the reasons for paying special attention to two-dimensional discrete competitive and cooperative systems is their applicability to mathematical models in biology and economics, the former involves competition or cooperation between two species.Another reason is that the theory of two-dimensional discrete competitive and cooperative systems is very well developed, unlike such theory for three and higher dimensional systems.Part of the reason for this situation is de Mottoni-Schiaffino theorem given below, which provides relatively simple scenarios for possible behavior of many two-dimensional discrete competitive and cooperative systems.However, this does not mean that one cannot encounter chaos in such systems as has been shown by Smith, see [].
y ≥ v}, and so on.Define the south-east partial order se on R  by (x, y) se (s, t) if and only if x ≤ s and y ≥ t.
Similarly, we define the north-east partial order ne on R  by (x, y) ne (s, t) if and only if x ≤ s and y ≤ t.For A ⊂ R  and x ∈ R  , define the distance from x to A as dist(x, A) := inf{ x -y : y ∈ A}.By int A we denote the interior of a set A.
It is easy to show that a map F is competitive if it is non-decreasing with respect to the south-east partial order, that is, if the following holds: For standard definitions of attracting fixed point, saddle point, stable manifold, and related notions, see [, ].
We now state three results for competitive maps in the plane.The following definition is from [].
Definition  Let S be a nonempty subset of R  .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 de Mottoni and Schiaffino [] for the Poincaré map of a periodic competitive Lotka-Volterra system of differential equations.Smith generalized the proof to competitive and cooperative maps [].
Theorem  Let S be a nonempty subset of R  .If T is a competitive map for which (O+) holds, then, for all x ∈ S, {T n (x)} is eventually component-wise 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 n (x)} is eventually component-wise 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.
The following result is from [], with the domain of the map specialized to be the cartesian product of intervals of real numbers.It gives a sufficient condition for conditions (O+) and (O-).
The following result is a direct consequence of the trichotomy theorem of Dancer and Hess, see [] and [], and is helpful for determining the basins of attraction of the equilibrium points.
Corollary  If the nonnegative cone of is a generalized quadrant in R n , and if T has no fixed points in Ju  , u  K other than u  and u  , then the interior of Ju  , u  K is either a subset of the basin of attraction of u  or a subset of the basin of attraction of u  .
The next result is a well-known global attractivity result which holds in partially ordered Banach spaces as well, see [].http://www.advancesindifferenceequations.com/content/2014/1/301Theorem  Let T be a monotone map on a closed and bounded rectangular region R ⊂ R  .Suppose that T has a unique fixed point ē in R. Then ē is a global attractor of T on R.
The following theorems were proved by Kulenović and Merino [] for competitive systems in the plane, when one of the eigenvalues of the linearized system at an equilibrium (hyperbolic or non-hyperbolic) is by an absolute value smaller than , while the other has an arbitrary value.These results are useful for determining basins of attraction of fixed points of competitive maps.

Theorem  Let T be a competitive map on a rectangular region
) is nonempty (i.e., x is not the NW or SE vertex of R), and T is strongly competitive on .Suppose that the following statements are true.
(a) The map T has a C  extension to a neighborhood of x.
(b) The Jacobian J T (x) of T at x has real eigenvalues λ, μ such that  < |λ| < μ, where |λ| < , and the eigenspace E λ associated with λ is not a coordinate axis.
Then there exists a curve C ⊂ R through x that is invariant and a subset of the basin of attraction of x such that C is tangential to the eigenspace E λ at x, and C is the graph of a strictly increasing continuous function of the first coordinate on an interval.Any endpoints of C in the interior of R are either fixed points or minimal period-two points.In the latter case, the set of endpoints of C is a minimal period-two orbit of T .
The situation where the endpoints of C are boundary points of R is of interest.The following result gives a sufficient condition for this case.
Theorem  For the curve C of Theorem  to have endpoints in ∂R, it is sufficient that at least one of the following conditions is satisfied.
(i) The map T has no fixed points nor periodic points of minimal period two in .
(ii) The map T has no fixed points in , det J T (x) > , and T(x) = x has no solutions x ∈ .(iii) The map T has no points of minimal period two in , det J T (x) < , and T(x) = x has no solutions x ∈ .
The next result is useful for determining basins of attraction of fixed points of competitive maps.(i) W -is invariant, and dist(T n (x), (iii) For every x ∈ W -, there exists (iv) For every x ∈ W + , there exists If T is a map on a set R and if x is a fixed point of T, the stable set W s (x) of x is the set {x ∈ R : T n (x) → x} and unstable set W u (x) of x is the set When T is non-invertible, the set W s (x) may not be connected and made up of infinitely many curves, or W u (x) may not be a manifold.The following result gives a description of the stable and unstable sets of a saddle point of a competitive map.If the map is a diffeomorphism on R, the sets W s (x) and W u (x) are the stable and unstable manifolds of x.
Theorem  In addition to the hypotheses of part (B) of Theorem , suppose that μ >  and that the eigenspace E μ associated with μ is not a coordinate axis.If the curve C of Theorem  has endpoints in ∂R, then C is the stable set W s (x) of x, and the unstable set W u (x) of x is a curve in R that is tangential to E μ at x and such that it is the graph of a strictly decreasing function of the first coordinate on an interval.Any endpoints of W u (x) in R are fixed points of T.
The following result gives information on local dynamics near a fixed point of a map when there exists a characteristic vector whose coordinates have negative product and such that the associated eigenvalue is hyperbolic.This is a well-known result, valid in a much more general setting: we include it here for completeness.A point (x, y) is a subsolution if T(x, y) se (x, y), and (x, y) is a supersolution if (x, y) se T(x, y).An order interval J(a, b), (c, d)K is the cartesian product of the two compact intervals [a, c] and [b, d].
Theorem  Let T be a competitive map on a rectangular set R ⊂ R  with an isolated ) , v () ) ∈ R  be an eigenvector of the Jacobian of T at x, with associated eigenvalue μ ∈ R. If v () v () < , then there exists an order interval I which is also a relative neighborhood of x such that for every relative neighborhood U ⊂ I of x the following statements are true.

Some basic facts
In this section we give some basic facts which will be used later.The map T associated to system () is given by

Equilibrium points
The equilibrium points (x, ȳ) of system () satisfy the equations By eliminating x =  from (), we get Similarly, we can eliminate variable ȳ from system () to obtain Then the following statements hold: Proof The proof of (a) follows from Descartes' rule of signs.Let The following matrix, called the discrimination matrix of f (y) and f (y) in [], is actually the Sylvester matrix of f (y) and f (y) with some permuted rows.
Let D k denote the determinant of the submatrix of Discr( f ), formed by the first k row and the first k columns, for k = , . . ., m.So, by a straightforward calculation, one can see that Assume that D  > .The sign list of the sequence

Injectivity, (O+) and (O-)
Lemma  Assume that (x, ȳ) is an equilibrium of the map T .Then the following hold: is invariant under the map T. Furthermore, the following holds: By substituting this into equation (), we obtain Since a  = A  c  , equations () and () become . By substituting this into () we get from which the proof follows.
() The Jacobian matrix of the map T has the form The determinant of () at any point is equal to The proof of () and () follows from Theorem .

Linearized stability analysis
The determinant of () at the equilibrium point is given by The trace of () at the equilibrium point is given by The characteristic equation has the form can be given explicitly as functions of y: Note that x g -(x) is always negative.Let x g (y) denote x g + (y).We consider only x f (y) and x g (y).Let x(y) = x f (y)x g (y).http://www.advancesindifferenceequations.com/content/2014/1/301 Lemma  Let T = (f , g) be the map defined by ().Then f x (x, y) > , and the following is true: Proof The first derivative of x f (y) is given by Since f y (x, y) < , we get f x (x, y) > .Further, Now, the proof follows from Lemma  Let T be the map defined by (), and let be the Jacobian matrix of T at a fixed point (x, ȳ).Then the Jacobian matrix () has real and distinct eigenvalues λ  and λ  such that |λ  | < λ  and λ  > .Furthermore, the following holds: Proof Implicit differentiation of the equations defining C f and C g at (x, ȳ) gives The characteristic equation associated with the Jacobian matrix of T at (x, ȳ) is given by Since the map T is competitive, then the eigenvalues of the Jacobian matrix of the map T at the equilibrium (x, ȳ) are real and distinct and |λ  | < λ  .By (), we have Theorem  Assume that  >  and λ  and λ  are eigenvalues of J T (x, ȳ).Then there exists the unique equilibrium point E = (x, ȳ) and the following hold: (b) If (ȳ) > , then E is a saddle point.Furthermore, the following hold: Proof In view of () and Lemma , we have that the function has one zero ȳ of multiplicity one.In view of Lemma , the map T has a unique equilibrium point.Since f () = -a  b   <  and lim y→+∞ f (y) = +∞, we have f (y) <  for y < ȳ and f (y) >  for y > ȳ.By Lemmas  and  from [], the equilibrium curves C f and C g intersect transversally at (x, ȳ), i.e., x (ȳ) = .In view of Lemma  and by the continuity of function x(y), there exists a neighborhood U ȳ of ȳ such that x (y) >  for y ∈ U ȳ, which implies () From () and Lemma  we obtain λ  <  and where p(λ) = (λλ  )(λλ  ).In view of () and p(-) = (λ  + )(λ  + ), we obtain statement (b) of the theorem.
Lemma  Suppose that all the assumptions of Theorem  are satisfied.Let Then the following statements are true.
(a) (ȳ) >  if and only if one of the following inequalities holds:

and only if the following hold:
Proof The function f (y) has one simple zero ȳ, which implies f (y) <  for  ≤ y < ȳ and f (y) >  for y > ȳ.Then . Now the proof follows from the fact that f (y) =  has real roots  and λ ()  are the eigenvalues of Proof In view of Lemma , equation () has three positive roots of multiplicity one.Since then by (a) of Lemma  we obtain that the map T has three equilibrium points that we denote by E  , E  and E  .Given points lie on the increasing curve which implies that the points are in the north-east ordering.Descartes' rule of signs and () imply that det J T (x, y) >  when a  < A  c  .In view of () and Lemma , we have that the polynomial has three zeros ȳi , i = , , , of multiplicity one.Since f () = -a  b   <  and lim y→+∞ f (y) = +∞, we have f (y) <  for y ∈ (, ȳ ) ∪ (ȳ  , ȳ ) and f (y) >  for y ∈ (ȳ  , ȳ ) ∪ (ȳ  , +∞).
By Lemmas  and  from [], the equilibrium curves C f and C g intersect transversally at E  , E  and E  , i.e., x (ȳ i ) = , i = , , .By this and Lemma  and by the continuity of function x(y), there exists a neighborhood U (i)  ȳi of ȳi such that x (y) >  for y ∈ U (i) ȳi for i = ,  and x (y) <  for y ∈ U ()  ȳi .Using this we get In view of (), we have det  .This completes the proof.
Theorem  Assume that  =  and  =  Then there exist two distinct equilibrium points in the positive quadrant E  = (x  , ȳ ) and E  = (x  , ȳ ) such that E  ne E  .Let λ (i)  and λ (i)   be the eigenvalues of J T (E i ), i = , .Then the following hold: (a) Exactly one of the roots ȳ or ȳ of () has multiplicity two.http://www.advancesindifferenceequations.com/content/2014/1/301(b) If ȳ is a root of () of multiplicity two, then the equilibrium point E  is non-hyperbolic and E  is a saddle point.Furthermore, Proof In view of Lemma , equation () has two positive zeros, one of multiplicity one and another one of multiplicity two, which implies statement (a).Since x = (A  + ȳ )/b  > , we obtain that the map T has two equilibrium points that we denote by E  and E  .Descartes' rule of signs and () imply that a  < A  c  ⇒ det J T (x, y) > .Now, we prove statement (b).
Similarly as in the proof of Theorem , one can see that E  is a saddle point.In view of Lemmas  and , from [] we have that x (ȳ  ) = , since ȳ is the root of () of multiplicity two.By Lemma  we obtain λ ()  = , λ ()  > .The proof of statement (c) is similar and we will skip it.
Theorem  Assume that  =  and  = .Then there exists one equilibrium point in the positive quadrant Proof In view of Lemma , ȳ is zero of () of multiplicity three.In view of Lemmas  and , from [] we have that x (ȳ  ) = .The rest of the proof is similar to that in Theorem  and we skip it.

Period-two solution
Let where and Period-two solution {( , ), T( , )} satisfies the system () http://www.advancesindifferenceequations.com/content/2014/1/301 The Jacobian matrix of the map T  at (x, y) has the form The determinant of () at (x, y) is given by The trace of () at (x, y) is given by = y} be the period-two curves, that is, the curves the intersection of which is a period-two solution.Then, for all y > , there exist exactly one x F (y) >  and exactly one x G (y) >  such that F(x F (y), y) = x and G(x G (y), y) = y.Furthermore, x F (y) and x G (y) are continuous functions and x F (y) > .
Proof Since F(x, y) = x and G(x, y) = y if and only if respectively, in view of Descartes' rule of signs, we have that for all y >  there exist exactly one x F (y) >  and exactly one x G (y) >  such that F(x F (y), y) = x and G(x G (y), y) = y.Taking derivatives of F(x, y) = x with respect to y, we get Since F y (x, y) < , we get x F (y) > .
Theorem  If a  ≤ A  c  , then T has no minimal period-two solution.If a  > A  c  and T has a minimal period-two solution {( , ), T( , )}, then {( , ), T( , )} is unstable.If μ  and μ  (μ  < μ  ) are the eigenvalues of J T  ( , ), then μ  >  and μ  > .All period-two solutions are ordered with respect to the north-east ordering.http://www.advancesindifferenceequations.com/content/2014/1/301Proof If a  ≤ A  c  , the statement follows from Lemma .If a  > A  c  , then from the first equation of () we have that By substituting this into () we obtain The rest of the proof follows from the fact that tr are algebraic curves.By using software Mathematica, one can see that the resultant of the polynomials F(x, y) and G(x, y) in variable x is given by Suppose that {(  ,  ), T(  ,  )} is a non-hyperbolic minimal period-two solution.This as in Lemma , one can see that where p  (μ) is the characteristic equation of the matrix From Theorem  we have that  < μ  < μ  and μ  > .Since {(  ,  ), T(  ,  )} is nonhyperbolic, we obtain that μ  = , from which it follows that x F (  )x G (  ) = .Since R( F, G) ≡ , we have that C F and C G have no common component.By Lemmas  and , from [], the curves C F and C G intersect transversally at (  ,  ) (i.e., y F (  )y G(  ) = ) if and only if  is zero of p(y) of multiplicity one.By Theorem . [], p(y) has zeros of multiplicity greater than one if and only if the discriminant D(p) of the polynomial p(y) is equal to zero, which proves the first statement of the lemma.
Assume that {(  ,  ), T(  ,  )} and {(  ,  ), T(  ,  )} are two minimal period-two solutions such that T has no other minimal period-two solutions in J(  ,  ), (  ,  )K = {(x, y) : (  ,  ) ne (x, y) ne (  ,  )} and D(p) = .From the previous discussion we ) < .Indeed assume, for example, that x F (  )x G (  ) <  and x F (  )x G (  ) < .Then there exists >  such that x F (y)x G (y) <  for y ∈ (  ,  + ) and x F (y)x G (y) >  for y ∈ (  -,  ).Since x F (y)x G (y) is a continuous function, this implies that there exists ∈ (  ,  ) such that x F ( )x G ( ) = , which is a contradiction.The rest of the proof follows from the fact that e i >  and g i < , i = , .
Proof By Theorem  we have that T has one equilibrium point E(x, ȳ), which is a repeller.This and Lemma  imply that T n (x  , y  ) is asymptotic to either (, ∞) or (∞, ), or a minimal period-two solution, for all (x  , y  ) ∈ R. Let B(, ∞) be the basin of attraction of (, ∞), and let B(∞, ) be the basin of attraction of (∞, ).By using Theorem  one can prove that int(Q Let S  denote the boundary of B(∞, ) considered as a subset of Q  (E), and let S  denote the boundary of B(∞, ) considered as a subset of Q  (E).It is easy to see that E ∈ S  , E ∈ S  and T(R) ⊂ int(R).Now we prove the following claim.
Claim  Let S  and S  be the sets defined as above.Then (a) Proof (a) The statement follows from T n (x  , y  ) se T n (x  , y  ) se (∞, ) and The claim (b) follows from the observation that there exists a ball centered at (x  , y  ) with the property that all its points (x, y) satisfy (x, y) se (x  , y  ).But one of these points necessarily lies in B(∞, ), so by (a) there exists (x  , y  ) ∈ B(∞, ).Furthermore, there exists a ball centered at (x  , y  ) with the property that all its points (x, y) satisfy (x, y) ∈ B(∞, ), which implies (x  , y  ) ∈ int(B)(∞, ).
In view of Claim , we have that (S  ∪ S  , ne ) is a totally ordered set which is invariant under T. If (x  , y  ) ∈ S  ∪ S  , then {T (n) (x  , y  )} is eventually component-wise monotone.Then there exists a minimal period-two solution {( , ), T( , )} ∈ S  ∪ S  ⊂ Q  (E) ∪ Q  (E) such that T (n) (x  , y  ) → ( , ) as n → ∞.By Theorem , {( , ), T( , )} is a non-hyperbolic or a saddle point.Assume that T has no minimal period-two solutions which are non-hyperbolic points and, for example, that ( , ), T( , ) ∈ S  such that E ne ( , ) ne T( , ).Since {( , ), T( , )} is a saddle point, in view of Theorems ,  and , we have that the global stable manifolds W s ({( , ), T( , )}) are the union of two curves W s ( , ) and W s (T( , )) whose endpoints are repeller points such that T(W s ( , )) = W s (T( , )) and E ne W s ( , ) ne W s (T( , )).If P  and P  (P  ne P  ) are endpoints of W s ( , ), then T(P  ) and T(P  ) are endpoints of T(W s ( , )) and either T(P  ) ne T(P  ) or T(P  ) ne T(P  ).Assume, for example, that P  ne W s ( , ) ne P  ne T(P  ) ne W s (T( , )) ne T(P  ).By Theorem  between two repellers P  and T(P  ), there exists a saddle point S  where its stable manifold is the union of two invariant curves W s (S  ) and W s (T(S  )) whose endpoints are repellers such that P  ne W s (S  ) ne T(P  ).Continuing in this way, we obtain that T has infinitely many minimal period-two solutions {P i , T(P i )}, which is in contradiction with the fact that T has at most eleven minimal period-two solutions.Hence ( , ) ne E ne T( , ).
Corollary  Assume that a  > A  c  , (ȳ) <  and D(p) = .Then there exists one equilibrium point E which is a repeller.Further, the set int(Q  (E)) ∪ int(Q  (E)) contains an odd number of minimal period-two solutions {( i , i ), ( ˜ i , ˜ i )}, i = , . . ., n + , such that ( i+ , i+ ) ne ( i , i ) ne E and E ne ( ˜ i , ˜ i ) ne ( ˜ i+ , ˜ i+ ), where ( ˜ i , ˜ i ) = T( i , i ).Furthermore, odd indexed period-two points are saddles and even indexed period-two points are repellers.
Proof By Theorem  we have that T has one equilibrium point E(x, ȳ), which is a repeller.By Theorem  all minimal period-two solutions are hyperbolic and the number of minimal period-two solutions is finite.In view of Claim , T has at least one minimal period-two solution which is a saddle point.Let {P  , T(P  )} be a minimal periodtwo solution which is a saddle point such that P  ne E ne T(P  ) and T has no minimal period-two solutions in JE, T(P  )K and JP  , EK.Such a minimal period-two solution exists in view of Theorem .The map T  satisfies all conditions of Theorems ,  and , which yields the existence of the global stable manifolds W s ({P  , P }) which are the union of two curves W s (P  ) and W s ( P ) that have a common endpoint E. If T has minimal period-two solutions in int(Q  (T(P  ))) ∪ int(Q  (P  )), let {P  , P } (P  ne P ) denote minimal periodtwo solutions such that T has no other minimal period-two solutions in JT(P  ), T(P  )K and JP  , P  K. Then W s (P  ) has the second endpoint at P  and W s (T(P  )) has the second endpoint at T(P  ) and P  ne P  ne E ne T(P  ) ne T(P  ).Furthermore, a minimal period-two solution {P  , T(P  )} is a repeller.Similarly as in Theorem , one can prove that int(Q  (T(P  ))) ∪ int(Q  (P  )) contains at least one minimal period-two solution which is a saddle point.Since the number of minimal period-two solutions is finite, continuing in this way, we will end with a minimal period-two solution which is a saddle point, from which the proof follows.

Corollary  If a
Furthermore, even indexed period-two points are saddles and odd indexed period-two points are repellers.
Proof The proof is similar as the proof of Corollary  and it will be omitted.
Based on a series of numerical simulations, we propose the following conjecture.
Conjecture  System () has at most one minimal period-two solution.

Global behavior
In this section we present global dynamics of system () in different parametric regions.We have five parametric regions with different dynamics which will be characterized by the following five theorems.

The case where the equilibrium points and period-two solutions are hyperbolic points ( 1 = 0 and D(p) = 0)
Theorem  Assume that  < .Then system () has three equilibrium solutions E  ne E  ne E  , where E  and E  are saddle points and E  is a repeller.In this case there exist four invariant continuous curves have end points at E  , and are graphs of increasing functions.The curves W u (E  ), W u (E  ) are the graphs of decreasing functions.Every solution {(x n , y n )} which starts below W s (E  ) ∪ W s (E  ) in the south-east ordering is asymptotic to (, ∞), and every solution {(x n , y n )} which starts above W s (E  ) ∪ W s (E  ) in the south-east ordering is asymptotic to (∞, ).The first quadrant of initial condition Q  = {(x  , y  ) : x  > , y  ≥ } is the union of five disjoint basins of attraction, i.e., where Proof Theorem  implies that there exist three equilibrium points, namely E  , E  and E  , such that E  ne E  ne E  .In this case, E  and E  are saddle points and E  is a repeller.In view of (), the map T is competitive on R and strongly competitive on int(R).It follows from the Perron-Frobenius theorem and a change of variables [] that at each point the Jacobian matrix of a strongly competitive map has two real and distinct eigenvalues, the larger one in absolute value being positive, and that corresponding eigenvectors may be http://www.advancesindifferenceequations.com/content/2014/1/301chosen to point in the direction of the second and first quadrant, respectively.Also, one can show that if the map is strongly competitive, then no eigenvector is aligned with a coordinate axis.Since  <  implies a  < A  c  , we have that det J T (E i ) > , i = , .Hence, all conditions of Theorems ,  and  are satisfied, which yields the existence of the global stable manifolds W s (E  ), W s (E  ) and the global unstable manifolds W u (E  ), W u (E  ), where W s (E  ), W s (E  ) are passing through the point E  and are graphs of increasing functions.The curves W u (E  ), W u (E  ) are the graphs of decreasing functions.Let In view of Theorem , since E  and E  are saddle points, we obtain that for all (x  , The sequence {u n } is monotone decreasing and bounded from below and {v n } is monotone increasing.Since T has no other equilibrium point except E  , E  and E  , this implies that the sequence {u n } has finite limit and R, the proof is similar and we skip it. Another way of completing the proof is by using () of Lemma .
 and (ȳ) >  and T has no minimal period-two solution, then system () has one equilibrium solution E(x, ȳ), which is a saddle point.There exists the global stable manifold W s (E) which is the graph of a continuous increasing function and the global unstable manifold W u (E), which is the graph of a continuous decreasing function.Every solution {(x n , y n )} which starts below W s (E) in the south-east ordering is asymptotic to (, ∞), and every solution {(x n , y n )} which starts above W s (E) in the south-east ordering is asymptotic to (∞, ).The first quadrant of initial condition Proof The conditions of this theorem and a  = A  c  imply all assumptions of Theorems ,  and  for R. The proof of this theorem in this case is similar to the proof of Theorem  and will be skipped.Now, we assume that a  = A  c  .In view of Lemma , the set is invariant and contains equilibrium E, and T(x, y) = E for (x, y) ∈ I.In view of the uniqueness of the global stable manifold, we conclude that W s (E) = I.Let W -= (x, y)|(x, y) se (x  , ỹ ) for some (x  , ỹ ) ∈ W s (E) ; Take any point (x  , y  ) ∈ W + (E).Then there exists the point (x l , y l ) ∈ I such that (x l , y l ) se (x  , y  ).Since the map T is strongly competitive, then , then there exists the point (x r , y r ) ∈ I such that (x  , y  ) se (x r , y r ).Since the map T is strongly competitive, then Theorem  If a  > A  c  and (ȳ) < , then there exists one equilibrium point E which is a repeller Furthermore, the odd indexed period-two points are saddles and the even indexed period-two points are repellers.The http://www.advancesindifferenceequations.com/content/2014/1/301global stable manifolds are given by where W s ( k+ , k+ ) and W s ( ˜ k+ , ˜ k+ ) are the graphs of a continuous strictly increasing function such that W s ( ˜ k+ , ˜ k+ ) = T(W s ( k+ , k+ )) with endpoints at ( ˜ k , ˜ k ), ( ˜ k+ , ˜ k+ ) and ( k+ , k+ ), ( k , k ), k = , . . ., n -, respectively, where separates R into two components W -and W + , which are basins of attraction of (, ∞) and (∞, ), respectively, where Further, for k = , . . ., n, we have B ( k+ , k+ ), ( ˜ k+ , ˜ k+ ) = W s ( k+ , k+ ) ∪ W s ( ˜ k+ , ˜ k+ ); B ( k , k ), ( ˜ k , ˜ k ) = ( k , k ), ( ˜ k , ˜ k ) .

The global unstable manifolds are given by
where W u ( k+ , k+ ) and W u ( ˜ k+ , ˜ k+ ) are the graphs of continuous strictly decreasing functions such that W u ( ˜ k+ , ˜ k+ ) = T(W u ( k+ , k+ )) with endpoints at (, ∞) and (∞, ).
Furthermore, the odd indexed period-two points are saddles and the even indexed period-two points are repellers.The map T  satisfies all conditions of Theorems ,  and , which yields the existence of the global stable and unstable manifolds with the above properties.In view of Theorem , for (x  , y  ) ∈ W -∩ R, there exists n  >  such that and for (x  , y  ) ∈ W + ∩ R, there exists n  >  such that Similarly as in the proof of Theorem , one can prove that  Theorem  Assume that  = ,  = .Then there exist two equilibrium points E  (x  , ȳ ) and E  (x  , ȳ ) such that E  ne E  , and the following hold: (a) If E  is non-hyperbolic, then there exist two invariant curves C +  and C +  that are contained in Q  (E  ) with endpoint in ∂Q  (E  ), which are the graphs of continuous strictly increasing functions.Further, there exists the global stable manifold W s (E  ) with endpoint at E  , which is a graph of a continuous increasing function and the global unstable manifold W u (E  ), which is a graph of a continuous decreasing function.Every solution {(x n , y n )} which starts below W s (E  ) ∪ C +  in the south-east ordering is asymptotic to (, ∞), and every solution {(x n , y n )} which starts above W s (E  ) ∪ C +  in the south-east ordering is asymptotic to (∞, ).The first quadrant of initial condition Q  = {(x  , y  ) : (b) If E  is non-hyperbolic, then there exist two invariant curves C -  and C -  that are contained in Q  (E  ) with endpoint in ∂Q  (E  ), which are the graphs of continuous strictly increasing functions.Further, there exists the global stable manifold W s (E  ) with endpoint at E  , which is a graph of a continuous increasing function, and the global unstable manifold W u (E  ), which is a graph of a continuous decreasing function.Every solution {(x n , y n )} which starts below W s (E  ) ∪ C -  in the south-east ordering is asymptotic to (, ∞), and every solution {(x n , y n )} which starts above W s (E  ) ∪ C -  in the south-east ordering is asymptotic to (∞, ).The first quadrant of initial condition Q  = {(x  , y  ) : Proof We will prove statement (a).The proof of statement (b) is similar.By Theorem  we have that T has two equilibrium points E  and E  such that E  ne E  .Then either E  is non-hyperbolic and E  is a saddle point or E  is non-hyperbolic and E  is a saddle point.
Assume that E  is non-hyperbolic.From Descartes' rule of signs applied to () we have that a  < A  c  .Lemma  implies that T n (x  , y  ) is asymptotic to either (, ∞) or (∞, ) or to an equilibrium point for all (x  , y  ) ∈ R. Let B(, ∞) be the basin of attraction of (, ∞), and let B(∞, ) be the basin of attraction of (∞, ).By using Theorem  one can prove that  and E ∈ C +  .Now Theorems ,  and  yield the existence of the global stable manifold W s (E  ) and the global unstable manifold W u (E  ), where W s (E  ) has an endpoint at E  and it is a graph of an increasing function.Similarly as in the proof of Theorem , one can see that Similarly as in Claim , one can prove the following claim.http://www.advancesindifferenceequations.com/content/2014/1/301Claim  Let C +  and C +  be the sets defined as above.Then we have the following: )) for all (x  , y  ) se (x  , y  ) and (x  , y  ) / ∈ B(∞, ) for all (x  , y  ) se (x  , y  ); The corresponding assertion for C +  is proved in a similar fashion.Now, we prove that C +  and C +  are the graphs of continuous strictly increasing functions.Let J be the projection of C +  onto the first coordinate.From (g), no two distinct elements of C +  can be weakly related nor can they have the same projection onto the first coordinate.Further, J contains [x  , x  ] whenever x  ≤ x  belong to J. Indeed, if x  ≤ x  belong to J, then there exist y  and y  such that (x  , y  ), (x  , y  ) ∈ C +  .From (c), we have that int(Q  (x  , y  )) ⊂ int(B(, ∞)) and int(Q  (x  , y  )) ⊂ int(B(∞, )).By using this it is easy to see that (x, ŷ) ∈ C +  if and only if ŷ = inf{y : y ∈ [y  , y  ] and (x, y) ∈ int(B(, ∞))}, where x ∈ [x  , x  ].This implies that J contains [x  , x  ].By (e), C +  ∩ int(Q  (E  )) = .From this it follows that J is an interval such that int(J) = ∅ and C +  is a connected set.Since points on C +  are non-comparable, C +  is the graph of a strictly increasing function f (x) of x ∈ J.If there is a jump discontinuity at x  ∈ J, let y -and y + respectively be the left and right (distinct) limits of f (x) as x approaches x  , respectively.The points (x  , y -) and (x  , y + ) are comparable in se -ordering, which is in contradiction with (g).Thus f (x) is a continuous function.Since T has no equilibrium points in Q  (E  ), similarly as in the proof of Theorem  [], one can prove that C +  has its endpoints in ∂Q  (E  ).The proof that considers C +  is similar and will be skipped.Take Based on a series of numerical simulations, we propose the following conjecture.See Figure  for visual illustration of Theorem .
Conjecture  Suppose that all assumptions of Theorem  are satisfied, then the following hold: Theorem  Assume that  = ,  =  and a  < A  c  .Then there exist the unique nonhyperbolic equilibrium point E(x, ȳ) and two invariant curves C -and C + with endpoint in http://www.advancesindifferenceequations.com/content/2014/1/301∂R, which are the graphs of continuous strictly increasing functions.Every solution {(x n , y n )} which starts below C  in the south-east ordering is asymptotic to (, ∞) and every solution {(x n , y n )} which starts above C  in the south-east ordering is asymptotic to (∞, ).The first Proof The proof is similar to the proof of Theorem  and will be skipped.
Based on a series of numerical simulations, we propose the following conjecture.Conjecture  Suppose that all assumptions of Theorem  are satisfied, then the following holds: Appendix: Values of coefficients p i for i = 0, . . ., 22  http://www.advancesindifferenceequations.com/content/2014/1/301 Theorem  (A) Assume the hypotheses of Theorem , and let C be the curve whose existence is guaranteed by Theorem .If the endpoints of C belong to ∂R, then C separates R into two connected components, namely W -:= {x ∈ R \ C : ∃y ∈ C with x se y} and W + := {x ∈ R \ C : ∃y ∈ C with y se x}, such that the following statements are true.
(a) Consider equation ().Then all its real roots are positive numbers.Furthermore, equation () has one, two, or three real roots.(b) If  > , then equation () has one real root and two pairs of distinct conjugate imaginary roots.(c) If  < , then equation () has three distinct real roots and one pair of conjugate imaginary roots.(d) If  =  and  = , then equation () has one pair of conjugate imaginary roots and two real roots, one real root of multiplicity one and other one of multiplicity two.(e) If  =  and  = , then equation () has one pair of conjugate imaginary roots and one real root of multiplicity three.http://www.advancesindifferenceequations.com/content/2014/1/301 ) from which it follows that the number of sign changes of the revised sign list of list () is two.Now, statement (b) follows in view of Theorem  [].Assume that D  < .If D  ≥ , then we obtain that f (y) has three pairs of conjugate imaginary roots, which is a contradiction.Hence, D  < .The sign list of the sequence {D  , D  , D  , D  , D  } is given by [, -, -, -, -], () which implies that the number of sign changes of the revised sign list of () is one.Now, statement (c) follows in view of Theorem  [].Similarly, one can prove statements (d) and (e).http://www.advancesindifferenceequations.com/content/2014/1/301 and Lemma .Theorem  If the map T has a minimal period-two solution{(  ,  ), T(  ,  )}, which is non-hyperbolic, then D(p) = , where D(p) is the discriminant of the polynomial p(y) := p  y  + p  y  + • • • + p  y + p  ,where the coefficients p i , i = , . . ., , are in the Appendix.If {(  ,  ), T(  ,  )} and {(  ,  ), T(  ,  )} are two minimal period-two solutions such that T has no other minimal period-two solutions in J(  ,  ), (  ,  )K = {(x, y) : (  ,  ) ne (x, y) ne (  ,  )} and D(p) = , then one of them is a saddle point and the other is a repeller.Proof Period-two solution curves C F = {(x, y) ∈ R : F(x, y) = } and C G = {(x, y) ∈ R : G(x, y) = }, where implies that F(  ,  ) = , G(  ,  ) = .By Theorem . [], F and G have a common non-constant factor if and only R( F, G) = , which implies that system F(x, y) = , G(x, y) =  has a solution if and only if R( F, G) = .Since f (  ) = , it must be p(  ) = .Similarly http://www.advancesindifferenceequations.com/content/2014/1/301
from which the proof follows in this case.See Figure  for visual illustration of Theorems  and .

Figure 2
Figure 2 Visual illustration of Theorem 19.
which completes the proof.See Figure  for visual illustration of Theorem .

5. 2
The case where at least one of the equilibrium points is non-hyperbolic and all minimal period-two solutions are hyperbolic ( 1 = 0 and D(p) = 0)

Figure 3
Figure 3 Visual illustration of statements (a) and (b) of Theorem 20.

Figure 4
Figure 4 Visual illustration of Theorem 21.
If, in addition to the hypotheses of part (A), x is an interior point of R and T is C  and strongly competitive in a neighborhood of x, then T has no periodic points in the boundary of Q  (x) ∪ Q  (x) except for x, and the following statements are true.