# Global Stability of a Discrete Mutualism Model

Authors: Kun Yang; Xiangdong Xie; Fengde Chen
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

A discrete mutualism model is studied in this paper. By using the linear
approximation method, the local stability of the interior equilibrium of the system
is investigated. By using the iterative method and the comparison principle of
difference equations, sufficient conditions which ensure the global asymptotical stability of the interior equilibrium of the system are obtained. The conditions which ensure the local stability of the positive equilibrium is enough to ensure the global attractivity are proved.


---

## Body

## 1. Introduction

There are many examples where the interaction of two or more species is to the advantage of all; we call such a situation the mutualism. For example, cellulose of white ants’ gut provides nutrients for flagellates, while flagellates provide nutrients for white ants through the decomposition of cellulose to glucose. As was pointed out by Chen et al. [1] “the mutual advantage of mutualism or symbiosis can be very important. As a topic of theoretical ecology, even for two species, this area has not been as widely studied as the others even though its importance is comparable to that of predator-prey and competition interactions.” Thus, it seems interesting to study some relevant topics on the symbiosis system.The following model was proposed by Chen et al. [1] to describe the mutualism mechanism:
(1)
d
N
1
d
t
=
r
1
N
1
(
t
)
[
K
1
+
α
1
N
2
(
t
)
1
+
N
2
(
t
)
-
N
1
(
t
)
]
,
d
N
2
d
t
=
r
2
N
2
(
t
)
[
K
2
+
α
2
N
1
(
t
)
1
+
N
1
(
t
)
-
N
2
(
t
)
]
,
where r
i refers to the intrinsic rate of population N
i and α
i
>
K
i
(
i
=
1,2
). In the absence of other species, the carrying capacity of the species N
i is K
i. Thanks to the cooperation of the other species, the carrying capacity of the species N
i becomes (
K
i
+
α
i
N
3
-
i
)
/
(
1
+
N
3
-
i
).Li [2] argued that the nonautonomous one is more appropriate, and he proposed the following two-species cooperative model:
(2)
d
N
1
d
t
=
r
1
(
t
)
N
1
(
t
)
×
[
K
1
(
t
)
+
α
1
(
t
)
N
2
(
t
-
τ
2
(
t
)
)
1
+
N
2
(
t
-
τ
2
(
t
)
)
-
N
1
(
t
-
σ
1
(
t
)
)
]
,
d
N
2
d
t
=
r
2
(
t
)
N
2
(
t
)
×
[
K
2
(
t
)
+
α
2
(
t
)
N
1
(
t
-
τ
1
(
t
)
)
1
+
N
1
(
t
-
τ
1
(
t
)
)
-
N
2
(
t
-
σ
2
(
t
)
)
]
,
where r
i
,
K
i
,
α
i
,
τ
i
,
σ
i
∈
C
(
R
,
R
+
), α
i
>
K
i, r
i
,
K
i
,
α
i
,
τ
i
,
σ
i
(
i
=
1,2
) are periodic functions of period ω
>
0. Here the author incorporates the time delays to the model, which means that the cooperation effect needs to spend some time to realize, but not immediately realize. By applying the coincidence degree theory, Li showed that the system has at least one positive periodic solution. For more works related to the system (1) and (2), one could refer to [1–9] and the references cited therein.It is well known that the discrete time models governed by difference equations are more appropriate than the continuous ones when the populations have nonoverlapping generations, and discrete time models can also provide efficient computational models of continuous models for numerical simulations. Corresponding to system (2), Li [10] proposed the following delayed discrete model of mutualism:
(3)
x
1
(
k
+
1
)
=
x
1
(
k
)
exp
⁡
{
r
1
(
k
)
[
K
1
(
k
)
+
α
1
(
k
)
x
2
(
k
-
τ
2
(
k
)
)
1
+
x
2
(
k
-
τ
2
(
k
)
)
-
x
1
(
k
-
σ
1
(
k
)
)
K
1
(
k
)
+
α
1
(
k
)
x
2
(
k
-
τ
2
(
k
)
)
1
+
x
2
(
k
-
τ
2
(
k
)
)
]
}
,
x
1
(
k
+
1
)
=
x
1
(
k
)
exp
⁡
{
r
1
(
k
)
[
K
2
(
k
)
+
α
2
(
k
)
x
1
(
k
-
τ
1
(
k
)
)
1
+
x
1
(
k
-
τ
1
(
k
)
)
-
x
2
(
k
-
σ
2
(
k
)
)
K
2
(
k
)
+
α
2
(
k
)
x
1
(
k
-
τ
1
(
k
)
)
1
+
x
1
(
k
-
τ
1
(
k
)
)
]
}
.
Under the assumption that r
i
,
K
i
,
α
i
,
τ
i
,
σ
i
(
i
=
1
,
2
) are positive periodic sequences with a common cycle ω, and α
i
>
K
i holds, by applying coincidence degree theory, he showed that system (3) has at least one positive ω-periodic solution, where ω is a positive integer. Chen [11] argued that the general nonautonomous nonperiodic system is more appropriate, and he showed that the system (3) is permanent. For more work about cooperative system, we can refer to [12–19].It brings to our attention that neither Li [10] nor Chen [11] investigated the stability property of the system (3), which is one of the most important topics on the study of population dynamics. We mention here that, with σ
i
(
k
)
≠
0, system (3) is a pure-delay system, and it is not an easy thing to investigate the stability property of the system. This motivated us to discuss a simple system, that is, the following autonomous cooperative system:
(4)
x
1
(
k
+
1
)
=
x
1
(
k
)
exp
⁡
{
r
1
[
K
1
+
α
1
x
2
(
k
)
1
+
x
2
(
k
)
-
x
1
(
k
)
]
}
,
x
2
(
k
+
1
)
=
x
2
(
k
)
exp
⁡
{
r
2
[
K
2
+
α
2
x
1
(
k
)
1
+
x
1
(
k
)
-
x
2
(
k
)
]
}
,
where x
i
(
k
)
(
i
=
1,2
) are the population density of the ith species at k-generation. r
i
,
K
i, i
=
1
,
2, are all positive constants.Throughout this paper, we assume that the coefficients of system (4) satisfy(
H
1
)
r
i
,
K
i
,
α
i
(
i
=
1
,
2
) are all positive constants and α
i
>
K
i
(
i
=
1,2
).The aim of this paper is, by further developing the analysis technique of [15], to obtain a set of sufficient conditions to ensure the global asymptotical stability of the interior equilibrium of system (4). More precisely, we will prove the following result.Theorem 1.
In addition to(
H
1
), further assume that(
H
2
)
r
i
α
i
≤
1
,
(
i
=
1,2
)holds; then the unique positive equilibrium (
x
1
*
,
x
2
*
) of the system (4) is globally asymptotically stable.The rest of the paper is arranged as follows. In Section2 we will introduce a useful lemma and investigate the local stability property of the positive equilibrium. With the help of several useful lemmas, the global attractivity of positive equilibrium of the system (4) is investigated in Section 3. An example together with its numeric simulation is presented in Section 4 to show the feasibility of our results. We end this paper by a brief discussion.


## 2. Local Stability

In view of the actual ecological implications of system (4), we assume that the initial value x
i
(
0
)
>
0
(
i
=
1
,
2
). Obviously, any solution of system (4) with positive initial condition is well defined on Z
+, where Z
+
=
{
0,1
,
2
,
…
} and remains positive for all n
≥
0.We determine the positive equilibrium of the system (4) through solving the following equations:
(5)
x
1
=
x
1
exp
⁡
[
r
1
(
K
1
+
α
1
x
2
1
+
x
2
-
x
1
)
]
,
x
2
=
x
2
exp
⁡
[
r
2
(
K
2
+
α
2
x
1
1
+
x
1
-
x
2
)
]
,
which is equivalent to
(6)
x
1
=
K
1
+
α
1
x
2
1
+
x
2
,
x
2
=
K
2
+
α
2
x
1
1
+
x
1
,
and so
(7)
(
α
1
+
1
)
x
2
2
+
(
1
-
α
1
α
2
+
K
1
-
K
2
)
x
2
-
(
K
2
+
α
2
K
1
)
=
0
.
Since α
1
+
1
>
0, K
2
+
α
2
K
1
>
0, (7) admits a unique positive solution x
¯
2. From the first equation of system (6), one could obtain x
¯
1; therefore, system (4) admits a unique positive equilibrium E
+
(
x
¯
1
,
x
¯
2
).Following we will discuss the local stability of equilibriumE
+
(
x
¯
1
,
x
¯
2
). The Jacobian matrix of system (4) at E
+
(
x
¯
1
,
x
¯
2
) is as follows:
(8)
J
(
E
+
)
=
(
1
-
r
1
x
¯
1
r
2
x
¯
2
α
2
-
K
2
(
1
+
x
¯
1
)
2
r
1
x
¯
1
α
1
-
K
1
(
1
+
x
¯
2
)
2
1
-
r
2
x
¯
2
)
.
The characteristic equation of J
(
E
+
) is
(9)
F
(
λ
)
=
λ
2
+
B
λ
+
C
=
0
,
where
(10)
B
=
-
(
2
-
r
1
x
¯
1
-
r
2
x
¯
2
)
,
C
=
(
1
-
r
1
x
¯
1
)
(
1
-
r
2
x
¯
2
)
-
r
1
r
2
x
¯
1
x
¯
2
(
α
1
-
K
1
)
(
α
2
-
K
2
)
(
1
+
x
¯
1
)
2
(
1
+
x
¯
2
)
2
.Lemma 2 (see [20]).
LetF
(
λ
)
=
λ
2
+
B
λ
+
C
=
0, where B and C are constants. Suppose F
(
1
)
>
0 and λ
1
,
λ
2 are two roots of F
(
λ
)
=
0. Then |
λ
1
|
<
1 and |
λ
2
|
<
1 if and only if F
(
-
1
)
>
0 and C
<
1.Letλ
1 and λ
2 be the two roots of (9), which are called eigenvalues of equilibrium E
+
(
x
¯
1
,
x
¯
2
). From [20] we know that if |
λ
1
|
<
1 and |
λ
2
|
<
1, then E
+
(
x
¯
1
,
x
¯
2
) is locally asymptotically stable.Theorem 3.
Assume that(
H
1
) and (
H
2
) hold; then E
+
(
x
¯
1
,
x
¯
2
) is locally asymptotically stable.Proof.
Since (7) has two real solutions, it follows that its discriminant is positive; that is,
(11)
Δ
=
(
1
-
α
1
α
2
+
K
1
-
K
2
)
2
+
4
(
α
1
+
1
)
(
K
2
+
α
2
K
1
)
=
(
1
+
α
1
α
2
+
K
1
+
K
2
)
2
-
4
(
α
1
-
K
1
)
(
α
2
-
K
2
)
>
0
,
and so
(12)
(
1
+
α
1
α
2
+
K
1
+
K
2
)
2
>
4
(
α
1
-
K
1
)
(
α
2
-
K
2
)
.
Equation (6) combined with the above inequality implies that
(13)
(
1
+
x
¯
1
)
2
(
1
+
x
¯
2
)
2
=
(
1
+
K
1
+
α
1
x
¯
2
1
+
x
¯
2
)
2
(
1
+
x
¯
2
)
2
=
(
(
α
1
+
1
)
x
¯
2
+
(
K
1
+
1
)
)
2
=
(
1
+
α
1
α
2
+
K
1
+
K
2
+
Δ
2
)
2
>
(
1
+
α
1
α
2
+
K
1
+
K
2
)
2
4
>
4
(
α
1
-
K
1
)
(
α
2
-
K
2
)
4
=
(
α
1
-
K
1
)
(
α
2
-
K
2
)
.
The above inequality leads to
(14)
1
-
(
α
1
-
K
1
)
(
α
2
-
K
2
)
(
1
+
x
¯
1
)
2
(
1
+
x
¯
2
)
2
>
0
.
From (9) and (14) it follows that
(15)
F
(
1
)
=
1
+
B
+
C
=
r
1
r
2
x
¯
1
x
¯
2
(
1
-
(
α
1
-
K
1
)
(
α
2
-
K
2
)
(
1
+
x
¯
1
)
2
(
1
+
x
¯
2
)
2
)
>
0
.
Also, from (6) and (
H
2
), one has
(16)
x
¯
i
<
α
i
,
r
i
x
¯
i
<
1
(
i
=
1
,
2
)
.
According to inequality (14) and (16), we have
(17)
F
(
-
1
)
=
1
-
B
+
C
=
r
1
r
2
x
¯
1
x
¯
2
(
1
-
(
α
1
-
K
1
)
(
α
2
-
K
2
)
(
1
+
x
¯
1
)
2
(
1
+
x
¯
2
)
2
)
+
(
4
-
2
r
1
x
¯
1
-
2
r
2
x
¯
2
)
>
0
,
C
-
1
=
-
r
1
x
¯
1
-
r
2
x
¯
2
+
r
1
r
2
x
¯
1
x
¯
2
(
1
-
(
α
1
-
K
1
)
(
α
2
-
K
2
)
(
1
+
x
¯
1
)
2
(
1
+
x
¯
2
)
2
)
<
-
r
1
x
¯
1
(
1
-
r
2
x
¯
2
)
-
r
2
x
¯
2
<
0
.
From Lemma 2 we can obtain that E
+
(
x
¯
1
,
x
¯
2
) is locally asymptotically stable. This completes the proof of Theorem 3.


## 3. Global Stability

We will give a strict proof of Theorem1 in this section. To achieve this objective, we introduce several useful lemmas.Lemma 4 (see [20]).
Letf
(
u
)
=
u
exp
⁡
(
α
-
β
u
), where α and β are positive constants; then f
(
u
) is nondecreasing for u
∈
(
0
,
(
1
/
β
)
].Lemma 5 (see [20]).
Assume that sequence{
u
(
k
)
} satisfies
(18)
u
(
k
+
1
)
=
u
(
k
)
exp
⁡
(
α
-
β
u
(
k
)
)
,
k
=
1
,
2
,
…
,
where α and β are positive constants and u
(
0
)
>
0. Then(i)
ifα
<
2, then lim
⁡
k
→
+
∞
u
(
k
)
=
α
/
β;
(ii)
ifα
<
1, then u
(
k
)
≤
(
1
/
β
), k
=
2,3
,
….Lemma 6 (see [21]).
Suppose that functionsf
,
g
:
Z
+
×
[
0
,
∞
)
→
[
0
,
∞
) satisfy f
(
k
,
x
)
≤
g
(
k
,
x
)
(
f
(
k
,
x
)
≥
g
(
k
,
x
)
) for k
∈
Z
+ and x
∈
[
0
,
∞
) and g
(
k
,
x
) is nondecreasing with respect to x. If {
x
(
k
)
} and {
u
(
k
)
} are the nonnegative solutions of the difference equations
(19)
x
(
k
+
1
)
=
f
(
k
,
x
(
k
)
)
,
u
(
k
+
1
)
=
g
(
k
,
u
(
k
)
)
,
respectively, and x
(
0
)
≤
u
(
0
)
(
x
(
0
)
≥
u
(
0
)
), then
(20)
x
(
k
)
≤
u
(
k
)
(
x
(
k
)
≥
u
(
k
)
)
∀
k
≥
0
.Proof of Theorem1.
Let(
x
1
(
k
)
,
x
2
(
k
)
) be arbitrary solution of system (4) with x
1
(
0
)
>
0 and x
2
(
0
)
>
0. Denote
(21)
U
1
=
limsup
⁡
k
→
+
∞
x
1
(
k
)
,
V
1
=
liminf
⁡
k
→
+
∞
x
1
(
k
)
,
U
2
=
limsup
⁡
k
→
+
∞
x
2
(
k
)
,
V
2
=
liminf
⁡
k
→
+
∞
x
2
(
k
)
.
We claim that U
1
=
V
1
=
x
¯
1 and U
2
=
V
2
=
x
¯
2.
From the first equation of system (4), we obtain
(22)
x
1
(
k
+
1
)
≤
x
1
(
k
)
exp
⁡
{
r
1
α
1
-
r
1
x
1
(
k
)
}
,
k
=
0
,
1
,
2
,
…
.
Consider the auxiliary equation as follows:
(23)
u
(
k
+
1
)
=
u
(
k
)
exp
⁡
{
r
1
α
1
-
r
1
u
(
k
)
}
,
k
=
0
,
1
,
2
,
…
.
Because of 0
<
r
1
α
1
≤
1, according to (ii) of Lemma 5, we can obtain u
(
k
)
≤
(
1
/
r
1
) for all k
≥
2, where u
(
k
) is arbitrary positive solution of (23) with initial value u
(
0
)
>
0. From Lemma 4, f
(
u
)
=
u
exp
⁡
(
r
1
α
1
-
r
1
u
) is nondecreasing for u
∈
(
0
,
(
1
/
r
1
)
]. According to Lemma 6 we can obtain x
(
k
)
≤
u
(
k
) for all k
≥
2, where u
(
k
) is the solution of (23) with the initial value u
(
2
)
=
x
(
2
). According to (i) of Lemma 5, we can obtain
(24)
U
1
=
limsup
⁡
k
→
+
∞
x
1
(
k
)
≤
lim
⁡
k
→
+
∞
u
(
k
)
=
α
1
=
def
M
1
x
1
.
From the second equation of system (4), we obtain
(25)
x
2
(
k
+
1
)
≤
x
2
(
k
)
exp
⁡
{
r
2
α
2
-
r
2
x
2
(
k
)
}
,
k
=
0
,
1
,
2
,
…
.
Similar to the above analysis, we have
(26)
U
2
=
limsup
⁡
k
→
+
∞
x
2
(
k
)
≤
lim
⁡
k
→
+
∞
u
(
k
)
=
α
2
=
def
M
1
x
2
.
Then, for sufficiently small constant ε
>
0, there is an integer k
1
>
2 such that
(27)
x
1
(
k
)
≤
M
1
x
1
+
ε
,
x
2
(
k
)
≤
M
1
x
2
+
ε
∀
k
>
k
1
.
According to the first equation of system (4) we can obtain
(28)
x
1
(
k
+
1
)
≥
x
1
(
k
)
exp
⁡
{
r
1
K
1
-
r
1
x
1
(
k
)
}
.
Consider the auxiliary equation as follows:
(29)
u
(
k
+
1
)
=
u
(
k
)
exp
⁡
{
r
1
K
1
-
r
1
u
(
k
)
}
.
Since 0
<
r
1
K
1
≤
1, according to (ii) of Lemma 5, we can obtain u
(
k
)
≤
(
1
/
K
1
) for all k
≥
2, where u
(
k
) is arbitrary positive solution of (23) with initial value u
(
0
)
>
0. From Lemma 4, f
(
u
)
=
u
exp
⁡
(
r
1
K
1
-
r
1
u
) is nondecreasing for u
∈
(
0
,
(
1
/
K
1
)
]. According to Lemma 6 we can obtain x
(
k
)
≥
u
(
k
) for all k
≥
2, where u
(
k
) is the solution of (23) with the initial value u
(
k
1
)
=
x
(
k
1
). According to (i) of Lemma 5, we have
(30)
V
1
=
liminf
⁡
k
→
+
∞
x
1
(
k
)
≥
lim
⁡
k
→
+
∞
u
(
k
)
=
K
1
=
def
N
1
x
1
.
From the second equation of system (4), we obtain
(31)
x
2
(
k
+
1
)
≥
x
2
(
k
)
exp
⁡
{
r
2
K
2
-
r
2
x
2
(
k
)
}
.
Similar to the above analysis, we have
(32)
V
2
=
liminf
⁡
k
→
+
∞
x
2
(
k
)
≥
lim
⁡
k
→
+
∞
u
(
k
)
=
K
2
=
def
N
1
x
2
.
Then, for sufficiently small constant ε
>
0, there is an integer k
2
>
k
1 such that
(33)
x
1
(
k
)
≥
N
1
x
1
-
ε
,
x
2
(
k
)
≥
N
1
x
2
-
ε
∀
k
>
k
2
.
Equation (24) combined with the first equation of system (4) leads to
(34)
x
1
(
k
+
1
)
≤
x
1
(
k
)
exp
⁡
{
r
1
[
K
1
+
α
1
(
M
1
x
2
+
ε
)
1
+
(
M
1
x
2
+
ε
)
-
x
1
(
k
)
]
}
,
k
>
k
2
.
Similar to the analysis of (23) and (24), we have
(35)
U
1
=
limsup
⁡
k
→
+
∞
x
1
(
k
)
≤
K
1
+
α
1
(
M
1
x
2
+
ε
)
1
+
(
M
1
x
2
+
ε
)
.
Because of arbitrariness of ε
>
0, we have U
1
≤
M
2
x
1, where
(36)
M
2
x
1
=
K
1
+
α
1
M
1
x
2
1
+
M
1
x
2
<
M
1
x
1
.
Equation (24) combined with the second equation of system (4) leads to
(37)
x
2
(
k
+
1
)
≤
x
2
(
k
)
exp
⁡
{
r
2
[
K
2
+
α
2
(
M
1
x
1
+
ε
)
1
+
(
M
1
x
1
+
ε
)
-
x
2
(
k
)
]
}
,
k
>
k
2
.
Similar to the analysis of (23) and (24), we can obtain
(38)
U
2
=
limsup
⁡
k
→
+
∞
x
2
(
k
)
≤
K
2
+
α
2
(
M
1
x
1
+
ε
)
1
+
(
M
1
x
1
+
ε
)
.
Because of arbitrariness of ε
>
0, we have U
2
≤
M
2
x
2, where
(39)
M
2
x
2
=
K
2
+
α
2
M
1
x
1
1
+
M
1
x
1
<
M
1
x
2
.
Then, for sufficiently small constant ε
>
0, there is an integer k
3
>
k
2 such that
(40)
x
1
(
k
)
≤
M
2
x
1
-
ε
,
x
2
(
k
)
≤
M
2
x
2
-
ε
,
∀
k
>
k
3
.
Equation (27) combined with the first equation of system (4) leads to
(41)
x
1
(
k
+
1
)
≥
x
1
(
k
)
exp
⁡
{
r
1
[
K
1
+
α
1
(
N
1
x
2
-
ε
)
1
+
(
N
1
x
2
-
ε
)
-
x
1
(
k
)
]
}
,
k
>
k
3
.
From this, we can finally obtain
(42)
V
1
=
liminf
⁡
k
→
+
∞
x
1
(
k
)
≥
K
1
+
α
1
(
N
1
x
2
-
ε
)
1
+
(
N
1
x
2
-
ε
)
.
Because of the arbitrariness of ε, we have V
1
≥
N
2
x
1, where
(43)
N
2
x
1
=
K
1
+
α
1
N
1
x
2
1
+
N
1
x
2
>
K
1
=
N
1
x
1
.
Equation (27) combined with the first equation of system (4) leads to
(44)
x
2
(
k
+
1
)
≥
x
2
(
k
)
exp
⁡
{
r
1
[
K
2
+
α
2
(
N
1
x
1
-
ε
)
1
+
(
N
1
x
1
-
ε
)
-
x
1
(
k
)
]
}
,
k
>
k
3
.
From the above inequality we can obtain
(45)
V
2
=
liminf
⁡
k
→
+
∞
x
2
(
k
)
≥
K
2
+
α
2
(
N
1
x
1
-
ε
)
1
+
(
N
1
x
1
-
ε
)
.
Because of the arbitrariness of ε, we have V
2
≥
N
2
x
2, where
(46)
N
2
x
2
=
K
2
+
α
2
N
1
x
1
1
+
N
1
x
1
>
K
2
=
N
1
x
2
.
Then, for sufficiently small constant ε
>
0, there is an integer k
4
>
k
3 such that
(47)
x
1
(
k
)
≥
N
2
x
1
-
ε
,
x
2
(
k
)
≥
N
2
x
2
-
ε
∀
k
>
k
4
.
Continuing the above steps, we can get four sequences {
M
k
x
1
},  {
M
k
x
2
},  {
N
k
x
1
}, and {
N
k
x
1
} such that
(48)
M
k
x
1
=
K
1
+
α
1
M
k
-
1
x
2
1
+
M
k
-
1
x
2
,
M
k
x
2
=
K
2
+
α
2
M
k
-
1
x
1
1
+
M
k
-
1
x
1
;
(49)
N
k
x
1
=
K
1
+
α
1
N
k
-
1
x
2
1
+
N
k
-
1
x
2
,
N
k
x
2
=
K
2
+
α
2
N
k
-
1
x
1
1
+
N
k
-
1
x
1
.
Clearly, we have
(50)
N
k
x
1
≤
V
1
≤
U
1
≤
M
k
x
1
,
N
k
x
2
≤
V
2
≤
U
2
≤
M
k
x
2
,
k
=
0
,
1
,
2
,
…
.
Now, we will prove{
M
k
x
i
}
(
i
=
1,2
) is monotonically decreasing and {
N
k
x
i
}
(
i
=
1,2
) is monotonically increasing by means of inductive method.
First of all, it is clear thatM
2
x
i
≤
M
1
x
i, N
2
x
i
≥
N
1
x
i
(
i
=
1
,
2
). For i
≥
2, we assume that M
i
x
1
≤
M
i
-
1
x
1 and N
i
x
1
≥
N
i
-
1
x
1 hold; then
(51)
M
i
+
1
x
2
=
K
2
+
α
2
M
i
x
1
1
+
M
i
x
1
≤
K
2
+
α
2
M
i
-
1
x
1
1
+
M
i
-
1
x
1
=
M
i
x
2
,
(52)
N
i
+
1
x
2
=
K
2
+
α
2
N
i
x
1
1
+
N
i
x
1
≥
K
2
+
α
2
N
i
-
1
x
1
1
+
N
i
-
1
x
1
=
N
i
x
2
,
(53)
M
i
+
1
x
1
=
K
1
+
α
1
M
i
x
2
1
+
M
i
x
2
≤
K
1
+
α
1
M
i
-
1
x
2
1
+
M
i
-
1
x
2
=
M
i
x
1
,
(54)
N
i
+
1
x
1
=
K
1
+
α
1
N
i
x
2
1
+
N
i
x
2
≥
K
1
+
α
1
N
i
-
1
x
2
1
+
N
i
-
1
x
1
=
N
i
x
1
.
Equations (51)–(54) show that {
M
k
x
1
} and {
M
k
x
2
} are monotonically decreasing and {
N
k
x
1
} and {
N
k
x
2
} are monotonically increasing. Consequently, lim
⁡
k
→
+
∞
{
M
k
x
i
} and lim
⁡
k
→
+
∞
{
N
k
x
i
}
(
i
=
1,2
) both exist. Let
(55)
lim
⁡
k
→
+
∞
M
k
x
i
=
X
i
*
,
lim
⁡
k
→
+
∞
N
k
x
i
=
x
i
*
,
i
=
1
,
2
.
From (48) and (55), we have
(56)
X
1
*
=
K
1
+
α
1
X
2
*
1
+
X
2
*
,
X
2
*
=
K
2
+
α
2
X
1
*
1
+
X
1
*
.
From (49) and (55), we get
(57)
x
1
*
=
K
1
+
α
1
x
2
*
1
+
x
2
*
,
x
2
*
=
K
2
+
α
2
x
1
*
1
+
x
1
*
.
Equations (56) and (57) show that (
X
1
*
,
X
2
*
) and (
x
1
*
,
x
2
*
) are all solutions of system (6). However, system (6) has unique positive solution (
x
¯
1
,
x
¯
2
). Therefore
(58)
U
i
=
V
i
=
lim
⁡
k
→
+
∞
x
i
(
k
)
=
x
¯
i
,
i
=
1
,
2
.
That is, E
+
(
x
¯
1
,
x
¯
2
) is globally attractive.
From Theorem3, we get that equilibrium E
+
(
x
¯
1
,
x
¯
2
) is locally asymptotically stable. And so, E
+
(
x
¯
1
,
x
¯
2
) is globally asymptotically stable. This ends the proof of Theorem 1.


## 4. Example

In this section, we will give an example to illustrate the feasibility of the main result.Example. Consider the following cooperative system:
(59)
x
1
(
k
+
1
)
=
x
1
(
k
)
exp
⁡
{
2
[
0.3
+
0.4
x
2
(
k
)
1
+
x
2
(
k
)
-
x
1
(
k
)
]
}
,
x
2
(
k
+
1
)
=
x
2
(
k
)
exp
⁡
{
4
[
0.2
+
0.25
x
1
(
k
)
1
+
x
1
(
k
)
-
x
2
(
k
)
]
}
.
By calculating, we have that positive equilibrium E
+
(
x
¯
1
,
x
¯
2
)
=
(
0.317495,0.212049
), r
1
α
1
=
0.8
<
1, r
2
α
2
=
1, K
i
<
α
i
(
i
=
1,2
) and the coefficients of system (59) satisfy (
H
1
) and (
H
2
). From Theorem 1, positive equilibrium E
+
(
x
¯
1
,
x
¯
2
) is globally asymptotically stable. Numeric simulation also supports our finding (see Figure 1).Figure 1
Dynamic behaviors of the solution(
x
1
(
t
)
,
x
2
(
t
)
) of system (59), with the initial conditions (
x
1
(
0
)
,
x
2
(
0
)
)
=
(
0.1,0.25
)
,
(
0.28,0.3
)
,
(
0.15,0.14
), and (
0.4,0.1
), respectively.


## 5. Discussion

It is well known [6] that, for autonomous two-species Lotka-Volterra mutualism model, the conditions which ensure the existence of positive equilibrium are enough to ensure that the equilibrium is globally stable. However, for the two-species discrete Lotka-Volterra mutualism model, Lu and Wang [22] proved that a cooperative system cannot be permanent. That is, the dynamic behaviors of discrete Lotka-Volterra mutualism model are very different to the continuous ones.Recently, by using the iterative method, Xie et al. [8] showed that, for a mutualism model with infinite delay, conditions which ensure the permanence of the system are enough to ensure the global stability of the system. As a corollary of their result, one could draw the conclusion that system (1) admits a unique positive equilibrium, which is globally stable. One interesting issue is proposed. For the discrete type mutualism model (4), is there any relationship between the existence of positive equilibrium and the stability property of the positive equilibrium?In this paper, by using the linear approximation, comparison principle of difference equations, and method of iteration scheme, we showed that the conditions which ensure the local stability property of the positive equilibrium ((
H
2
)
0
<
r
i
α
i
≤
1
(
i
=
1,2
)) are also enough to guarantee the global stability of the positive equilibrium E
+
(
x
¯
1
,
x
¯
2
).At the end of this paper, we would like to mention here that, for the Lotka-Volterra type mutualism system with time delay, delay is one of the most important factors to influence the dynamic behaviors of the system [23–25]. It seems interesting to incorporate the time delay to the system (4) and investigate the dynamic behaviors of the system; we leave this for future study.



---

*Source: 709124-2014-06-23.xml*