1. Introduction
The system of difference equations,
where the parameters
are positive real numbers, and
and
and the initial conditions
and
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
on the global dynamics of system (
1). The global dynamics of (
1) was considered in the case where the parameters
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
, it has been shown in [
1] that under the condition
and
, the points
are equilibria of Equation (
1), and that for some values of the parameters there exists an additional equilibrium point
located in the open positive quadrant, given by
Important subsets of parameter space are described in
Table 1, together with corresponding behavior of equilibrium solutions established in [
1].
Theorem 4.1 in [
5] applies when parameters vary from Case 1 to Case 4 of
Table 1. Set
and define
to be the map of Equation (
1) restricted to
, that is,
Therefore, Theorem 4.1 in [
5] gives global behavior of solutions to system (
1) on
for
. In particular, a bifurcation occurs when the equilibrium
E changes its local character from a locally stable equilibrium to a saddle point. This happens when the parameters cross the critical surface
.
It is also shown in [
1] that the open, positive semiaxis
is attracted to
, and that the open, positive semiaxis
is attracted to
. 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. - (i)
Suppose that . If , then is globally asymptotically stable on , and attracts all points on the open semiaxis . If , then the stable manifold in is the graph of a continuous, increasing function of the first coordinate. Furthermore, a solution converges to whenever is above in South-east ordering, and converges to whenever is below in South-east ordering.
- (ii)
Suppose that . If , then is globally asymptotically stable on , and attracts all points on the open semiaxis . If , then E is globally asymptotically stable on , attracts all points on the open semiaxis , and attracts all points on the open semiaxis .
See Figure 1 for graphical interpretation. The non-hyperbolic case when
was not considered in [
1]. When (
3) holds, a direct calculation gives that the fixed points of
are
and all points on the segment
, where
The eigenvalues of the Jacobian of
at
are
and corresponding eigenvectors are
It is shown in [
8] that, for system (
1), the hypotheses of Theorem 5 in [
6] are satisfied and that all solutions fall inside an invariant rectangular region. Therefore, every solution of (
1) converges to an equilibrium point. A direct calculation shows that the origin is a repeller. We conclude that every nonzero solution converges to a point
. 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
, then
is continuous. These observations, together with an application of Theorem 1 in [
6] lead to the following result.
Theorem 2. - (i)
Every nonzero solution to system (1) converges to an equlibrium . - (ii)
For every with and , the stable set is an unbounded, increasing curve with endpoint .
- (iii)
The limiting equilibrium varies continuously with the initial condition.
See Figure 1 for graphical interpretation. Statement (ii) excludes equilibria of the form
and
since the hypotheses of Theorem 1 in [
6] are not satisfied at these points.
In this paper, we consider two related systems, namely
and
where all present coefficients are positive and the initial conditions are non-negative and such that
. We derive the global dynamics of both systems (
4) and (
5), which explains the effect of the coefficients
on the global dynamics. Related systems are considered in [
2,
3,
10,
11,
12].
The paper is organized as follows.
Section 2 presents some basic preliminary results about competitive systems, which is our main tool in proving the results.
Section 3 contains the global dynamics of system (
4) and
Section 4 gives the global dynamics of system (
5). Actually we show that the global dynamics of all three systems is determined by their linearized dynamics. We use different techniques to prove the results for systems (
4) and (
5). In particular, we will use the fact that system (
4) is homogeneous to obtain the exact equation of the global stable manifold of positive equilibrium solution of system (
4). We will show that the presence of more parameters will create more dynamic scenarios in such a way that system (
1) exhibits one additional dynamic scenario compared to system (
5) which in turn will have one additional dynamic scenario more than system (
4). This is in contrast with the global dynamics of second order difference equation
with all non-negative parameters and initial conditions such that
, where the most complicated dynamics occurs in the special case of Lyness’ equation when
[
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 with nonempty interior, and let be a map (i.e., a continuous function). Set . The map T is competitive
if is non-decreasing in x and non-increasing in y, and 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 equationsis 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
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 the South-East partial order in the plane whose non-negative cone is the standard fourth quadrant , that is, if and only if and . The North-East partial order is defined analogously with the non-negative cone given by the standard first quadrant .
Competitive maps T in the plane preserve the South-East ordering: whenever . 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
, and in particular competitive maps in
, may have chaotic dynamics. Smale [
19] showed that any continuous time vector field on the standard
-simplex in
can be embedded on a smooth, competitive vector field in
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 with non-negative cone P. For the order interval is the set of all z such that . We say if and , and if , where is the interior of P. A map T on a subset of is order preserving if whenever , strictly order preserving if whenever , and strongly order preserving if whenever .
Let be a map with a fixed point and let be an invariant subset of R that contains . We say that is stable (asymptotically stable) relative to if is a stable (asymptotically stable) fixed point of the restriction of T to .
Definition 2. Let be a nonempty subset of . A competitive map is said to satisfy condition () if for every x, y in , implies , and T is said to satisfy condition () if for every x, y in , implies .
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 be a nonempty subset of . If T is a competitive map for which () holds then for all , is eventually componentwise monotone. If the orbit of x has compact closure, then it converges to a fixed point of T. If instead () holds, then for all , is eventually componentwise monotone. If the orbit of x has compact closure in , 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
. The next two results are stated for order-preserving maps on
. 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 and ⪯ a partial order on , let be an order preserving map, and let be such that and . If and , then is invariant and
- i.
There exists a fixed point of T in .
- ii.
If T is strongly order preserving, then there exists a fixed point in which is stable relative to .
- iii.
If there is only one fixed point in , then it is a global attractor in and therefore asymptotically stable relative to .
Corollary 1. If the non-negative cone of ⪯ is a generalized quadrant in , and if T has no fixed points in other than and , then the interior of is either a subset of the basin of attraction of or a subset of the basin of attraction of .
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
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 . - (b)
for every , where denotes the Jacobian matrix of the map T.
- (c)
, for every .
- (d)
Every solution of system (4) satisfies the difference equationwhere . - (e)
Ifthen the map T has an invariant line
Proof. The Jacobian matrix
of the map
T has the form
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
is exactly an equilibrium of the equation in (d). □
System (
4) always has two equilibrium solutions on the axes,
,
. It can also have either exactly one interior equilibrium solution
E or an infinite number of equilibrium solutions
. Since the interior equilibrium solution
E is an intersection of two equilibrium curves (isoclines)
it will exist if either
(
x-intercept of
smaller than
x-intercept of
and
y-intercept of
bigger than
y-intercept of
) or
(
x-intercept of
smaller than
x-intercept of
and
y-intercept of
bigger than
y-intercept of
). These two geometrical conditions can be unified as condition (
7).
Condition (
7) implies that
, in which case the interior equilibrium
is given as:
Notice that
,
implies
and
,
implies
. If either
,
and
, or
and
then there are no interior equilibrium points. Furthermore, if
and
, the two equilibrium curves
coincide and every point on the segment
is an equilibrium solution
. See
Table 2 summarizing the equilibrium points of system (
4).
The following result describes the local stability character of all equilibrium solutions.
Lemma 2. - (a)
The equilibrium solution is locally asymptotically stable if , non-hyperbolic of stable type if and a saddle point if . In each case, the eigenvectors associated with the eigenvalues and are and .
- (b)
The equilibrium solution is locally asymptotically stable if , non-hyperbolic of stable type if and a saddle point if . In each case, the eigenvectors associated with the eigenvalues and are and .
- (c)
The interior equilibrium solution E is a saddle point when and and is locally asymptotically stable when and .
- (d)
The interior equilibrium solutions are non-hyperbolic of the stable type and the eigenvector which corresponds to is given as .
Proof. - (a)
In view of (
9), we have
which implies that the eigenvalues of the Jacobian matrix are
. The corresponding eigenvectors are as stated.
- (b)
In view of (
9), we have
which implies that the eigenvalues of the Jacobian matrix are
. The corresponding eigenvectors are as stated.
- (c)
The eigenvalues of the Jacobian matrix evaluated at the equilibrium
E,
and
, correspond to the roots of the characteristic polynomial
. Note that
by (
9). Furthermore
Consequently, if
and
then
and
and if
and
then
and
. It follows that
E is a saddle point when
and
and is locally asymptotically stable when
and
.
- (d)
In this case, the eigenvalues of the Jacobian matrix evaluated at the equilibrium are . The eigenvector that corresponds to is , where satisfies and points towards the first quadrant.
□
Now, global behavior of system (
4) is described by the following result:
Theorem 5. - (a)
If , then the equilibrium solutions are locally asymptotically stable and the interior equilibrium E is a saddle point. The separatrix , which is a graph of a continuous, non-decreasing curve, is the basin of attraction of E and the region below (resp. above) is the basin of attraction of (resp. ).
- (b)
If , then the equilibrium solutions 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 (resp. ).
- (c)
If (resp. ), then the equilibrium solution (resp. ) is locally asymptotically stable and (resp. ) is a saddle point. The basin of attraction of (resp. ) is the first quadrant of initial conditions without the positive part of the y-axis (resp. x-axis), which is attracted by (resp. ).
- (d)
If and , then there is an infinite family of equilibrium solutions for which there exists the global stable manifold , which is the graph of a continuous, non-decreasing function asymptotic to and is exactly the basin of attraction of . The limiting equilibrium varies continuously with the initial condition.
- (e)
If (resp. ) is non-hyperbolic and (resp. ) is locally asymptotically stable then (resp. ) attracts the first quadrant of initial conditions except the positive part of x-axis (resp. y-axis) which is attracted by (resp. ). If (resp. ) is non-hyperbolic and (resp. ) is a saddle point then (resp. ) attracts the first quadrant of initial conditions except the positive part of y-axis (resp. x-axis) which is attracted by (resp. ).
See Figure 2 for graphical interpretation. Proof. - (a)
First we show that
T does not have any period-two solutions. Our condition implies
. By direct calculation one can show that a period-two solution satisfies the equation
Please note that
which means that both terms of such solution are complex conjugate and so there is no period-two solution in the first quadrant.
Taking into account that the Jacobian matrix evaluated at
E has all non-zero entries, Theorem 5 of [
4] implies the existence and uniqueness of both global stable and unstable manifolds
and
and so
. Furthermore, Theorem 5 of [
4] implies that every
below
will satisfy
for some
. In view of Corollary 1,
. In a similar way, we can treat the case when
is above
.
- (b)
In view of Lemma 2 part (a), the eigenvectors which correspond to and point to the interior of the fourth and the second quadrant, which means that the local unstable manifolds and exist and point strictly toward E. Thus there exist points in the interior of , arbitrarily close to and such that . Now, statement iii of Theorem 3 implies that E is a global attractor in , which completes the proof.
- (c)
Assume that
which implies that
is locally asymptotically stable and
is a saddle point. In view of Lemma 2 part (b), the eigenvector which corresponds to
points to the interior of the fourth quadrant, which means that the local unstable manifold
exists and points strictly toward
. Thus there exists a point
u in the interior of
, arbitrarily close to
such that
. However, then this shows that the map
T has a lower solution in every neighborhood of
, which in view of Theorem 6 in [
6] implies that the interior of
is a subset of the basin of attraction of
. The result follows.
The proof when is similar and will be omitted.
- (d)
By Theorem 1 of [
6], for each
there exists the set
passing through
and asymptotic to
, which is the graph of a continuous, non-decreasing function, which is exactly the basin of attraction of
. 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 system (4) has an explicit solution. Indeed in this case we have that , where and so the solution of system (4) isThus implies and implies while implies , 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 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
the map associated with system (
5).
Lemma 3. - (a)
Every solution of system (5) satisfies ; - (b)
satisfies () condition on and so T has no period-two points;
- (c)
For every , ;
- (d)
For every , .
Proof. Part (a) and (c) follow by immediate checking and part (d) follows by solving the resulting Beverton-Holt one-dimensional equation. Notice
This is equivalent to
This implies that
and
and thus
, that is
. □
Please note that a direct consequence of the above form of is that 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,
. Provided that
, there exists an equilibrium solution on the
x-axis,
. Depending on the values of the parameters
and
, there is also the possibility of either exactly one interior equilibrium solution
E or an infinite number of interior equilibrium solutions
. The interior equilibrium solution is an intersection of two equilibrium curves
and
. This solution will exist if either
and
(i.e., the
x-intercept of
is smaller than the
x-intercept of
and the
y-intercept of
is bigger than the
y-intercept of
) or
and
(i.e., the
x-intercept of
is smaller than the
x-intercept of
and the
y-intercept of
is bigger than the
y-intercept of
). As in (
7), these two geometrical conditions can be unified as
where
and
. Condition (
12) implies that
and the interior equilibrium
is given as:
Please note that if
,
and
, or
and
then there does not exist an interior equilibrium solution. Since
implies that
, we must have
in order for an interior equilibrium point to exist. Furthermore, if
then the two equilibrium curves
and
coincide and every point of the segment
,
is an equilibrium solution
. The equilibrium points for system (
5) are summarized in
Table 3.
The local stability character of , , E and are presented in Lemma 4. The proof requires Proposition 1.
Proposition 1. The eigenvalues λ and μ of are positive.
Proof. In view of (
14) we have
which implies
Please note that the equilibrium point
E exists under the hypothesis
, which means that either
and
or
and
. In either case, we have
, and consequently,
. Since system (
5) is strongly competitive, by the Perron-Frobenius Theorem [
15,
23], the largest eigenvalue of
is positive, which completes the proof. □
Lemma 4. - (a)
The equilibrium solution exists if . It is locally asymptotically stable if , non-hyperbolic of stable type if and a saddle point if . In each case, the eigenvectors associated with the eigenvalues and are and , respectively.
- (b)
The equilibrium solution always exists and it is locally asymptotically stable if , non-hyperbolic of stable type if and a saddle point if . In each case, the eigenvectors associated with the eigenvalues and are and , respectively.
- (c)
The interior equilibrium solution E exists if and it is locally asymptotically stable if and and a saddle point if and .
- (d)
The interior equilibrium solutions exist if , and . They are non-hyperbolic of the stable type and the eigenvector associated with where is .
Proof. - (a)
The Jacobian matrix associated with the map
has the form
In view of (
14), we have
which implies that the eigenvalues of the Jacobian matrix are
. The corresponding eigenvectors are as stated.
- (b)
In view of (
14), we have
which implies that the eigenvalues of the Jacobian matrix are
. The corresponding eigenvectors are as stated.
- (c)
Denote the eigenvalues of
by
and
, which represent the roots of the characteristic polynomial,
. By Proposition 1,
and
are real and positive. Notice
If
and
then
and by (
15),
. Combining this with the fact that
, it follows that
E is a saddle point. If
and
then
and from (
15) we have
and
. Combining this with
, we conclude that
E is locally asymptotically stable.
- (d)
For , the eigenvalues of are and . Since we clearly have that are non-hyperbolic equilibrium points of the stable type. It follows by immediate checking that the eigenvector associated with is , which points towards the first quadrant for .
□
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. - (a)
If then 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 and every solution which starts on the positive x-axis converges to the singular point . - (b)
For , if , (resp. , ) then system (5) has equilibrium solutions and where (resp. ) is locally asymptotically stable and (resp. ) is a saddle point. The basin of attraction of (resp. ) is the first quadrant of initial conditions without the positive part of the y-axis (resp. x-axis), which is attracted by (resp. ). - (c)
If , , and then system (5) has equilibrium solutions , and E. The equilibrium solutions and 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 (resp. ). - (d)
If , , and then system (5) has equilibrium solutions , and E. The equilibrium solutions and are locally asymptotically stable and the interior equilibrium E is a saddle point. There exists the global stable manifold and the global unstable manifold , where is the graph of a continuous, non-decreasing function and 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 is the basin of attraction of (resp. ) and the curve is the basin of attraction of E. - (e)
If , , and then there is an infinite family of equilibrium solutions for which there exists the global stable manifold for all , which is the graph of a continuous, non-decreasing function asymptotic to and is exactly the basin of attraction of . The limiting equilibrium varies continuously with the initial condition.
- (f)
If (resp. ) is non-hyperbolic and (resp. ) is locally asymptotically stable then (resp. ) attracts the first quadrant of initial conditions except the positive part of the x-axis (resp. y-axis), which is attracted by (resp. ). If (resp. ) is non-hyperbolic and (resp. ) is a saddle point then (resp. ) attracts the first quadrant of initial conditions except the positive part of the y-axis (resp. x-axis), which is attracted by (resp. ).
See Figure 3 for graphical interpretation. Proof. - (a)
Let
. Lemma 3(c) and (d) guarantee that for initial conditions on the positive
y-axis,
and for initial conditions on the positive
x-axis,
. To treat the dynamics in the interior of
, consider
. By Theorem 2 of [
4],
is invariant. The region
also attracts the interior of
. To verify this, suppose that
with
. In this case
and
It follows that there exists such that for all , and thus is attracting. To conclude the proof, suppose with . In this case , and as a consequence of the invariance of , is a decreasing sequence while is a non-decreasing sequence. Therefore, . The above arguments prove that the basins of attraction for and the singular point are and .
- (b)
Let
,
, and
. Lemma 3(c) and (d) guarantee that for all initial conditions on the positive
y-axis,
and for all initial conditions on the positive
x-axis,
. To treat the interior of
, consider
shown in
Figure 4.
Please note that
is an invariant region by Theorem 2 of [
4]. Consider
with
and notice
As a consequence of the invariance of
,
is a non-decreasing sequence and
is a non-increasing sequence. Therefore, using basic properties of sequences and the fact that
is strongly competitive,
. Finally, suppose
with
. By Lemma 3(a),
for all
. Choose
with
such that
. Since
is strongly competitive, notice
Therefore, . We have arrived at the desired result that the basins of attraction for and are and .
The proof for the case when , is similar and will be omitted.
- (c)
Let
,
, and
. As in part (b), Lemma 3(c) and (d) guarantee that the positive part of the
y-axis is a subset of
and the positive part of the
x-axis is a subset of
. To treat the interior of
, consider the region
shown in
Figure 5.
Please note that
is invariant by Theorem 2 of [
4]. Provided that
with
, monotonicity properties (similar to part (a) and (b)) along with Lemma 3(b) can be used to prove that
. Suppose
with
. By Lemma 3(a) we know that
for all
. Moreover, since
then
for all
. Consequently, there must exist an
such that
. Now, choose
such that
. Since
is strongly competitive we have
Therefore, for all . We have reached the desired result that the basins of attraction for and are , and .
- (d)
Let
,
, and
. In light of Lemma 4(c), Theorems 1 and 5 of [
6] guarantees that there exist the global stable and unstable manifolds for
E,
and
respectively, with the above mentioned properties. An immediate checking shows that
and that the interior of the ordered interval
is a subset of
, while the interior of the ordered interval
is a subset of
. Now, take any point
such that
(i.e., above
). Then
, where
. By Lemma 3(c) and the monotonicity of
, for
,
Since
, (
18) implies that
enters the ordered interval
and so converges to
. In a similar way, one can show that the ordered interval
attracts all points below
, and so all such points converge to
.
- (e)
By Theorem 1 of [
6], for each
there exists the set
passing through
and asymptotic to
, which is the graph of a continuous, non-decreasing function, which is exactly the basin of attraction of
. 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
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
which is possible because of the fact that
. In this case,
is a globally asymptotically stable equilibrium for system (
1), while the basin of attraction of the singular point
is an empty set for system (
4), and the basin of attraction of the singular point
is the non-negative part of the
x-axis for system (
5). The major difference between systems (
4) and (
5) is the case
which is possible only for system (
5) when the singular point
has non-empty basin of attraction. Thus the presence of
will introduce new dynamic scenario while the presence of both
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.