# Stability Analysis of a System of Exponential Difference Equations

Authors: Q. Din; K. A. Khan; A. Nosheen
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

We study the boundedness character and persistence, existence and uniqueness of positive equilibrium, local and global behavior, and rate of convergence of positive solutions of the following system of exponential difference equations:x
n
+
1
=
(
α
1
+
β
1
e
-
x
n
+
γ
1
e
-
x
n
-
1
)
/
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
), y
n
+
1
=
(
α
2
+
β
2
e
-
y
n
+
γ
2
e
-
y
n
-
1
)
/
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
), where the parameters α
i
,
β
i
,
γ
i
,
a
i
,
b
i, and c
i for i
∈
{
1,2
} and initial conditions x
0
,
x
-
1
,
y
0, and y
-
1 are positive real numbers. Furthermore, by constructing a discrete Lyapunov function, we obtain the global asymptotic stability of the positive equilibrium. Some numerical examples are given to verify our theoretical results.


---

## Body

## 1. Introduction

Many population models are governed by exponential difference equations. We refer to [1–6] and the references therein. Systems of nonlinear difference equations of higher-order are of paramount importance in applications. Such equations also appear naturally as discrete analogues and as numerical solutions of systems differential and delay differential equations which model diverse phenomena in biology, ecology, physiology, physics, engineering, and economics. For applications and basic theory of rational difference equations we refer to [7–9]. In [10–17] applications of difference equations in mathematical biology are given. It is very interesting to investigate the behavior of solutions of a system of nonlinear difference equations and to discuss the local asymptotic stability of their equilibrium points.El-Metwally et al. [1] investigated boundedness character, asymptotic behavior, periodicity nature of the positive solutions, and stability of equilibrium point of the following population model:
(1)
x
n
+
1
=
α
+
β
x
n
-
1
e
-
x
n
.
Papaschinopoulos et al. [2] studied the boundedness, the persistence, and the asymptotic behavior of positive solutions of the following two directional interactive and invasive species models:
(2)
x
n
+
1
=
a
+
b
x
n
-
1
e
-
y
n
,
y
n
+
1
=
c
+
d
y
n
-
1
e
-
x
n
.
Papaschinopoulos et al. [3] investigated the asymptotic behavior of the solutions of the following three systems of difference equations of exponential form:
(3)
x
n
+
1
=
α
+
β
e
-
y
n
γ
+
y
n
-
1
,
y
n
+
1
=
δ
+
ϵ
e
-
x
n
ζ
+
x
n
-
1
,
x
n
+
1
=
α
+
β
e
-
y
n
γ
+
x
n
-
1
,
y
n
+
1
=
δ
+
ϵ
e
-
x
n
ζ
+
y
n
-
1
,
x
n
+
1
=
α
+
β
e
-
x
n
γ
+
y
n
-
1
,
y
n
+
1
=
δ
+
ϵ
e
-
y
n
ζ
+
x
n
-
1
.
Recently, Papaschinopoulos and Schinas [4] studied the asymptotic behavior of the positive solutions of the systems of the two difference equations:
(4)
x
n
+
1
=
a
+
b
y
n
-
1
e
-
y
n
,
y
n
+
1
=
c
+
d
x
n
-
1
e
-
x
n
,
x
n
+
1
=
a
+
b
y
n
-
1
e
-
x
n
,
y
n
+
1
=
c
+
d
x
n
-
1
e
-
y
n
.
Motivated by the above study, our aim in this paper is to investigate the qualitative behavior of positive solutions of the following system of exponential difference equations:
(5)
x
n
+
1
=
α
1
+
β
1
e
-
x
n
+
γ
1
e
-
x
n
-
1
a
1
+
b
1
y
n
+
c
1
y
n
-
1
,
y
n
+
1
=
α
2
+
β
2
e
-
y
n
+
γ
2
e
-
y
n
-
1
a
2
+
b
2
x
n
+
c
2
x
n
-
1
,
where the parameters α
i, β
i, γ
i, a
i, b
i, and c
i for i
∈
{
1,2
} and initial conditions x
0
,
x
-
1, y
0, and y
-
1 are positive real numbers.More precisely, we investigate the boundedness character and persistence, existence and uniqueness of positive steady state, local asymptotic stability and global behavior of unique positive equilibrium point, and rate of convergence of positive solutions of system (5) which converge to its unique positive equilibrium point. Some special cases of system (5) can be treated as population models of two species [3].


## 2. Main Results

The following theorem shows that every solution of (5) is bounded and persists.Theorem 1.
Every positive solution{
(
x
n
,
y
n
)
} of system (5) is bounded and persists.Proof.
For any positive solution{
(
x
n
,
y
n
)
} of system (5), one has
(6)
x
n
+
1
≤
α
1
+
β
1
+
γ
1
a
1
=
U
1
,
y
n
+
1
≤
α
2
+
β
2
+
γ
2
a
2
=
U
2
,
p
p
p
p
p
p
p
p
p
n
=
0,1
,
2
,
…
.
Furthermore, from systems (5) and (6) we obtain that
(7)
x
n
+
1
≥
α
1
+
β
1
e
-
U
1
+
γ
1
e
-
U
1
a
1
+
b
1
U
2
+
c
1
U
2
=
L
1
,
y
n
+
1
≥
α
2
+
β
2
e
-
U
2
+
γ
2
e
-
U
2
a
2
+
b
2
U
1
+
c
2
U
1
=
L
2
,
p
p
p
p
p
p
p
p
p
p
p
p
p
p
p
p
i
n
=
2,3
,
…
.
From (6) and (7), it follows that
(8)
L
1
≤
x
n
≤
U
1
,
L
2
≤
y
n
≤
U
2
,
n
=
3,4
,
…
.
Hence, the theorem is proved.Lemma 2.
Let{
(
x
n
,
y
n
)
} be a positive solution of system (5). Then, [
L
1
,
U
1
]
×
[
L
2
,
U
2
] is invariant set for system (5).Proof.
The proof follows by induction.


### 2.1. Stability Analysis

Let us consider four-dimensional discrete dynamical system of the form:(9)
x
n
+
1
=
f
(
x
n
,
x
n
-
1
,
y
n
,
y
n
-
1
)
y
n
+
1
=
g
(
x
n
,
x
n
-
1
,
y
n
,
y
n
-
1
)
,
n
=
0,1
,
…
,
where f
:
I
2
×
J
2
→
I and g
:
I
2
×
J
2
→
J are continuously differentiable functions and I
,
J are some intervals of real numbers. Furthermore, a solution {
(
x
n
,
y
n
)
}
n
=
-
1
∞ of system (9) is uniquely determined by initial conditions (
x
i
,
y
i
)
∈
I
×
J for i
∈
{
-
1,0
}. Along with system (9) we consider the corresponding vector map F
=
(
f
,
x
n
,
g
,
y
n
). An equilibrium point of (9) is a point (
x
¯
,
y
¯
) that satisfies
(10)
x
¯
=
f
(
x
¯
,
x
¯
,
y
¯
,
y
¯
)
,
y
¯
=
g
(
x
¯
,
x
¯
,
y
¯
,
y
¯
)
.
The point (
x
¯
,
y
¯
) is also called a fixed point of the vector map F.Definition 3.
Let(
x
¯
,
y
¯
) be an equilibrium point of the system (9).(i)
An equilibrium point(
x
¯
,
y
¯
) is said to be stable if for every ɛ
>
0 there exists δ
>
0 such that for every initial condition (
x
i
,
y
i
), i
∈
{
-
1,0
}, if ∥
∑
i
=
-
1
0
‍
(
x
i
,
y
i
)
-
(
x
¯
,
y
¯
)
∥
<
δ implies ∥
(
x
n
,
y
n
)
-
(
x
¯
,
y
¯
)
∥
<
ɛ for all n
>
0, where ∥
·
∥ is usual Euclidian norm in R
2.
(ii)
An equilibrium point(
x
¯
,
y
¯
) is said to be unstable if it is not stable.
(iii)
An equilibrium point(
x
¯
,
y
¯
) is said to be asymptotically stable if there exists η
>
0 such that ∥
∑
i
=
-
1
0
‍
(
x
i
,
y
i
)
-
(
x
¯
,
y
¯
)
∥
<
η and (
x
n
,
y
n
)
→
(
x
¯
,
y
¯
) as n
→
∞.
(iv)
An equilibrium point(
x
¯
,
y
¯
) is called global attractor if (
x
n
,
y
n
)
→
(
x
¯
,
y
¯
) as n
→
∞.
(v)
An equilibrium point(
x
¯
,
y
¯
) is called asymptotic global attractor if it is a global attractor and stable.Definition 4.
Let(
x
¯
,
y
¯
) be an equilibrium point of a map F
=
(
f
,
x
n
,
g
,
y
n
), where f and g are continuously differentiable functions at (
x
¯
,
y
¯
). The linearized system of (9) about the equilibrium point (
x
¯
,
y
¯
) is
(11)
X
n
+
1
=
F
(
X
n
)
=
F
J
X
n
,
where X
n
=
(
x
n
y
n
x
n
-
1
y
n
-
1
) and F
J is Jacobian matrix of system (9) about the equilibrium point (
x
¯
,
y
¯
).To construct corresponding linearized form of system (5) we consider the following transformation:
(12)
(
x
n
,
y
n
,
x
n
-
1
,
y
n
-
1
)
⟼
(
f
,
g
,
f
1
,
g
1
)
,
where f
=
x
n
+
1, g
=
y
n
+
1, f
1
=
x
n, and g
1
=
y
n. The linearized system of (5) about (
x
¯
,
y
¯
) is given by
(13)
Z
n
+
1
=
F
J
(
x
¯
,
y
¯
)
Z
n
,
where Z
n
=
(
x
n
y
n
x
n
-
1
y
n
-
1
) and the Jacobian matrix about the fixed point (
x
¯
,
y
¯
) under transformation (12) is given by
(14)
F
J
(
x
¯
,
y
¯
)
=
(
A
1
A
2
A
3
A
4
B
1
B
2
B
3
B
4
1
0
0
0
0
1
0
0
)
,
where
(15)
A
1
=
-
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
,
A
2
=
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
,
A
3
=
-
γ
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
,
A
4
=
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
,
B
1
=
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
,
B
2
=
-
β
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
,
B
3
=
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
,
B
4
=
-
γ
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
.Lemma 5 (see [9]).
Assume thatX
n
+
1
=
F
(
X
n
), n
=
0,1
,
…, is a system of difference equations such that X
¯ is a fixed point of F. If all eigenvalues of the Jacobian matrix J
F about X
¯ lie inside the open unit disk |
λ
|
<
1, then X
¯ is locally asymptotically stable. If one of them has a modulus greater than one, then X
¯ is unstable.Theorem 6.
System (5) has a unique positive equilibrium point (
x
¯
,
y
¯
)
∈
[
L
1
,
U
1
]
×
[
L
2
,
U
2
], if the following condition is satisfied:
(16)
a
2
+
L
1
(
b
2
+
c
2
)
<
α
2
+
(
β
2
+
γ
2
)
e
-
K
K
,
where
(17)
K
=
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
L
1
(
b
1
+
c
1
)
.Proof.
Consider the following system of equations:(18)
x
=
α
1
+
(
β
1
+
γ
1
)
e
-
x
a
1
+
(
b
1
+
c
1
)
y
,
y
=
α
2
+
(
β
2
+
γ
2
)
e
-
y
a
2
+
(
b
2
+
c
2
)
x
.
From (18), it follows that
(19)
y
=
α
1
+
(
β
1
+
γ
1
)
e
-
x
-
a
1
x
(
b
1
+
c
1
)
x
,
x
=
α
2
+
(
β
2
+
γ
2
)
e
-
y
-
a
2
y
(
b
2
+
c
2
)
y
.
Set
(20)
F
(
x
)
=
α
2
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
-
a
2
f
(
x
)
(
b
2
+
c
2
)
f
(
x
)
-
x
,
where f
(
x
)
=
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
-
a
1
x
)
/
(
(
b
1
+
c
1
)
x
) and x
∈
[
L
1
,
U
1
]. Then, it follows that
(21)
F
(
L
1
)
=
L
1
(
b
1
+
c
1
)
×
(
a
2
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
L
1
(
b
1
+
c
1
)
α
2
+
(
β
2
+
γ
2
)
e
-
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
/
L
1
(
b
1
+
c
1
)
-
a
2
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
L
1
(
b
1
+
c
1
)
)
×
(
(
b
2
+
c
2
)
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
)
-
1
-
L
1
.
Furthermore, it is easy to see that
(22)
(
a
1
+
(
b
1
+
c
1
)
U
2
)
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
)
>
a
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
U
1
)
.
From (22) it follows that
(23)
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
>
0
.
Let K
=
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
/
(
L
1
(
b
1
+
c
1
)
)
>
0. Then F
(
L
1
) can be expressed as
(24)
F
(
L
1
)
=
α
2
+
(
β
2
+
γ
2
)
e
-
K
-
a
2
K
(
b
2
+
c
2
)
K
-
L
1
.
Suppose that a
2
+
L
1
(
b
2
+
c
2
)
<
(
α
2
+
(
β
2
+
γ
2
)
e
-
K
)
/
K; then, it follows that F
(
L
1
)
>
0. Furthermore, we have
(25)
F
(
U
1
)
=
U
1
(
b
1
+
c
1
)
(
a
2
(
β
1
+
γ
1
)
(
1
-
e
-
U
1
)
U
1
(
b
1
+
c
1
)
(
β
2
+
γ
2
)
e
(
β
1
+
γ
1
)
(
1
-
e
-
U
1
)
/
U
1
(
b
1
+
c
1
)
+
α
2
+
a
2
(
β
1
+
γ
1
)
(
1
-
e
-
U
1
)
U
1
(
b
1
+
c
1
)
)
×
(
(
b
2
+
c
2
)
(
α
1
+
e
-
U
1
(
β
1
+
γ
1
)
-
a
1
U
1
)
)
-
1
-
U
1
.
Then it is easy to see that
(26)
α
1
+
e
-
U
1
(
β
1
+
γ
1
)
-
a
1
U
1
=
(
β
1
+
γ
1
)
(
e
-
U
1
-
1
)
<
0
,
which gives that F
(
U
1
)
<
0. Hence, F
(
x
)
=
0 has at least one positive solution in [
L
1
,
U
1
]. Moreover, we obtain that
(27)
F
′
(
x
)
=
-
1
-
f
′
(
x
)
(
b
2
+
c
2
)
f
(
x
)
×
[
a
2
+
(
b
2
+
c
2
)
x
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
]
,
where
(28)
f
′
(
x
)
=
-
α
1
+
(
β
1
+
γ
1
)
(
1
+
x
)
e
-
x
(
b
1
+
c
1
)
x
2
.
Then, from (28) it follows that f
(
U
1
)
<
f
(
x
)
<
f
(
L
1
) and using (26) we obtain
(29)
F
′
(
x
)
=
-
1
+
α
1
+
(
β
1
+
γ
1
)
(
1
+
x
)
e
-
x
(
b
1
+
c
1
)
(
b
2
+
c
2
)
x
2
f
(
x
)
×
[
a
2
+
(
b
2
+
c
2
)
x
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
]
≤
-
1
+
α
1
+
(
β
1
+
γ
1
)
(
1
+
x
)
e
-
x
(
b
1
+
c
1
)
(
b
2
+
c
2
)
x
2
f
(
U
1
)
×
[
a
2
+
(
b
2
+
c
2
)
x
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
]
<
0
.
Hence, F
(
x
)
=
0 has a unique positive solution in [
L
1
,
U
1
]. The proof is therefore completed.Theorem 7.
The unique positive equilibrium point of system (5) is locally asymptotically stable under the following condition:
(30)
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
+
(
1
+
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
×
(
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
+
(
β
2
+
γ
2
)
e
-
L
2
a
2
+
(
b
2
+
c
2
)
L
1
)
+
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
×
(
(
b
1
+
c
1
)
(
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
)
(
a
1
+
(
b
1
+
c
1
)
L
2
)
2
+
γ
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
<
1
.Proof.
The characteristic polynomial of Jacobian matrixF
J
(
x
¯
,
y
¯
) about (
x
¯
,
y
¯
) is given by
(31)
P
(
λ
)
=
λ
4
-
(
A
1
+
B
2
)
λ
3
-
(
A
2
B
1
-
A
1
B
2
+
B
3
+
B
4
)
λ
2
-
(
(
A
3
+
A
4
)
B
1
-
(
B
3
+
B
4
)
A
1
)
λ
.
Clearly, one root of P
(
λ
) is 0. To check the behavior of the other three roots of P
(
λ
), we let Φ
(
λ
)
=
λ
3 and Ψ
(
λ
)
=
(
A
1
+
B
2
)
λ
2
+
(
A
2
B
1
-
A
1
B
2
+
B
3
+
B
4
)
λ
+
(
A
3
+
A
4
)
B
1
-
(
B
3
+
B
4
)
A
1. Assume that (30) holds and |
λ
|
=
1; then, one has
(32)
|
Ψ
(
λ
)
|
<
|
A
1
+
B
2
|
+
|
A
2
B
1
-
A
1
B
2
+
B
3
+
B
4
|
+
|
(
A
3
+
A
4
)
B
1
-
(
B
3
+
B
4
)
A
1
|
<
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
+
(
1
+
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
)
×
(
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
+
(
β
2
+
γ
2
)
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
)
+
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
×
(
(
b
1
+
c
1
)
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
+
γ
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
)
<
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
+
(
1
+
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
×
(
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
+
(
β
2
+
γ
2
)
e
-
L
2
a
2
+
(
b
2
+
c
2
)
L
1
)
+
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
×
(
(
b
1
+
c
1
)
(
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
)
(
a
1
+
(
b
1
+
c
1
)
L
2
)
2
+
γ
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
<
1
.
Then, by Rouche’s Theorem, Φ
(
λ
) and Φ
(
λ
)
-
Ψ
(
λ
) have the same number of zeroes in an open unit disk |
λ
|
<
1. Hence, all the roots of (31) satisfy |
λ
|
<
1, and it follows from Lemma 5 that the unique positive equilibrium point (
x
¯
,
y
¯
) of the system (5) is locally asymptotically stable.Theorem 8.
The unique positive equilibrium point(
x
¯
,
y
¯
) of system (5) is globally asymptotically stable, if the following condition is satisfied:
(33)
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
<
x
¯
(
a
1
+
(
b
1
+
c
1
)
L
2
)
,
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
<
y
¯
(
a
2
+
(
b
2
+
c
2
)
L
1
)
.Proof.
Arguing as in [18], we consider the following discrete time analogue of Lyapunov function:
(34)
V
n
=
x
¯
(
x
n
x
¯
-
1
-
ln
⁡
x
n
x
¯
)
+
y
¯
(
y
n
y
¯
-
1
-
ln
⁡
y
n
y
¯
)
.
Then nonnegativity of V
n follows from the following inequality:
(35)
x
-
1
-
ln
⁡
⁡
x
≥
0
∀
x
>
0
.
Furthermore, we have
(36)
-
ln
⁡
⁡
(
x
n
+
1
x
n
)
=
ln
⁡
⁡
(
1
-
(
1
-
x
n
x
n
+
1
)
)
≤
-
x
n
+
1
-
x
n
x
n
+
1
,
-
ln
⁡
⁡
(
y
n
+
1
y
n
)
=
ln
⁡
⁡
(
1
-
(
1
-
y
n
y
n
+
1
)
)
≤
-
y
n
+
1
-
y
n
y
n
+
1
.
Assume that (33) holds true; then, it follows that
(37)
V
n
+
1
-
V
n
=
x
¯
(
x
n
+
1
x
¯
-
1
-
ln
⁡
x
n
+
1
x
¯
)
+
y
¯
(
y
n
+
1
y
¯
-
1
-
ln
⁡
y
n
+
1
y
¯
)
-
x
¯
(
x
n
x
¯
-
1
-
ln
⁡
x
n
x
¯
)
-
y
¯
(
y
n
y
¯
-
1
-
ln
⁡
y
n
y
¯
)
≤
(
x
n
+
1
-
x
n
)
+
(
y
n
+
1
-
y
n
)
-
x
¯
x
n
+
1
(
x
n
+
1
-
x
n
)
-
y
¯
y
n
+
1
(
y
n
+
1
-
y
n
)
=
(
x
n
+
1
-
x
n
)
(
1
-
x
¯
x
n
+
1
)
+
(
y
n
+
1
-
y
n
)
(
1
-
y
¯
y
n
+
1
)
≤
(
U
1
-
L
1
)
×
(
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
-
x
¯
(
a
1
+
(
b
1
+
c
1
)
L
2
)
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
)
+
(
U
2
-
L
2
)
×
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
-
y
¯
(
a
2
+
(
b
2
+
c
2
)
L
1
)
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
≤
0
for all n
≥
0 so that V
n
≥
0 is monotonically decreasing sequence. It follows that lim
⁡
n
→
∞
V
n
≥
0. Hence, we obtain that
(38)
lim
⁡
n
→
∞
(
V
n
+
1
-
V
n
)
=
0
.
Then it follows that lim
⁡
n
→
∞
x
n
+
1
=
x
¯ and lim
⁡
n
→
∞
y
n
+
1
=
y
¯. Furthermore, V
n
≤
V
0 for all n
≥
0, which gives that (
x
¯
,
y
¯
)
∈
[
L
1
,
U
1
]
×
[
L
2
,
U
2
] is uniformly stable. Hence, unique positive equilibrium point (
x
¯
,
y
¯
)
∈
[
L
1
,
U
1
]
×
[
L
2
,
U
2
] of system (5) is globally asymptotically stable.


### 2.2. Rate of Convergence

In this section we will determine the rate of convergence of a solution that converges to the unique positive equilibrium point of the system (5).The following result gives the rate of convergence of solutions of a system of difference equations:(39)
X
n
+
1
=
(
A
+
B
(
n
)
)
X
n
,
where X
n is an m-dimensional vector, A
∈
C
m
×
m is a constant matrix, and B
:
Z
+
→
C
m
×
m is a matrix function satisfying
(40)
∥
B
(
n
)
∥
⟶
0
as n
→
∞, where ∥
·
∥ denotes any matrix norm which is associated with the vector norm
(41)
∥
(
x
,
y
)
∥
=
x
2
+
y
2
.Proposition 9 ((Perron’s Theorem) [19]).
Suppose that condition (40) holds. If X
n is a solution of (39), then either X
n
=
0 for all large n or
(42)
ρ
=
lim
⁡
n
→
∞
(
∥
X
n
∥
)
1
/
n
exists and is equal to the modulus of one of the eigenvalues of matrix A.Proposition 10 (see [19]).
Suppose that condition (40) holds. If X
n is a solution of (39), then either X
n
=
0 for all large n or
(43)
ρ
=
lim
⁡
n
→
∞
∥
X
n
+
1
∥
∥
X
n
∥
exists and is equal to the modulus of one of the eigenvalues of matrix A.Let{
(
x
n
,
y
n
)
} be any solution of the system (5) such that lim
⁡
n
→
∞
x
n
=
x
¯ and lim
⁡
n
→
∞
y
n
=
y
¯, where x
¯
∈
[
L
1
,
U
1
] and y
¯
∈
[
L
2
,
U
2
]. To find the error terms, one has from the system (5)
(44)
x
n
+
1
-
x
¯
=
α
1
+
β
1
e
-
x
n
+
γ
1
e
-
x
n
-
1
a
1
+
b
1
y
n
+
c
1
y
n
-
1
-
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
=
β
1
(
e
-
x
n
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
x
¯
)
(
x
n
-
x
¯
)
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
y
n
-
y
¯
)
+
γ
1
(
e
-
x
n
-
1
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
1
-
x
¯
)
(
x
n
-
1
-
x
¯
)
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
y
n
-
1
-
y
¯
)
,
y
n
+
1
-
y
¯
=
α
2
+
β
2
e
-
y
n
+
γ
2
e
-
y
n
-
1
a
2
+
b
2
x
n
+
c
2
x
n
-
1
-
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
=
β
2
(
e
-
y
n
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
x
¯
)
(
y
n
-
y
¯
)
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
x
n
-
x
¯
)
+
γ
2
(
e
-
y
n
-
1
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
1
-
y
¯
)
(
y
n
-
1
-
y
¯
)
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
x
n
-
1
-
x
¯
)
.
Let e
n
1
=
x
n
-
x
¯ and e
n
2
=
y
n
-
y
¯; then, one has
(45)
e
n
+
1
1
=
a
n
e
n
1
+
b
n
e
n
2
+
c
n
e
n
-
1
1
+
d
n
e
n
-
1
2
,
e
n
+
1
2
=
f
n
e
n
1
+
g
n
e
n
2
+
h
n
e
n
-
1
1
+
k
n
e
n
-
1
2
,
where
(46)
a
n
=
β
1
(
e
-
x
n
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
x
¯
)
,
b
n
=
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
,
c
n
=
γ
1
(
e
-
x
n
-
1
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
1
-
x
¯
)
,
d
n
=
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
,
f
n
=
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
,
g
n
=
β
2
(
e
-
y
n
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
x
¯
)
,
h
n
=
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
,
k
n
=
γ
2
(
e
-
y
n
-
1
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
1
-
y
¯
)
.
Moreover,
(47)
lim
⁡
n
→
∞
a
n
=
-
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
=
A
1
,
lim
⁡
n
→
∞
b
n
=
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
=
A
2
,
lim
⁡
n
→
∞
c
n
=
-
γ
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
=
A
3
,
lim
⁡
n
→
∞
d
n
=
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
=
A
4
,
lim
⁡
n
→
∞
f
n
=
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
=
B
1
,
lim
⁡
n
→
∞
g
n
=
-
β
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
=
B
2
,
lim
⁡
n
→
∞
h
n
=
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
=
B
3
,
lim
⁡
n
→
∞
k
n
=
-
γ
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
=
B
4
.
Now the limiting system of error terms can be written as
(48)
[
e
n
+
1
1
e
n
+
1
2
e
n
1
e
n
2
]
=
[
A
1
A
2
A
3
A
4
B
1
B
2
B
3
B
4
1
0
0
0
0
1
0
0
]
[
e
n
1
e
n
2
e
n
-
1
1
e
n
-
1
2
]
,
which is similar to linearized system of (5) about the equilibrium point (
x
¯
,
y
¯
).Using Proposition9, one has the following result.Theorem 11.
Assume that{
(
x
n
,
y
n
)
} is a positive solution of the system (5) such that lim
⁡
n
→
∞
x
n
=
x
¯, and lim
⁡
n
→
∞
y
n
=
y
¯, where x
¯
∈
[
L
1
,
U
1
] and y
¯
∈
[
L
2
,
U
2
]. Then, the error vector e
n
=
(
e
n
1
e
n
2
e
n
-
1
1
e
n
-
1
2
) of every solution of (5) satisfies both of the following asymptotic relations:
(49)
lim
⁡
n
→
∞
(
∥
e
n
∥
)
1
/
n
=
|
λ
1,2
,
3,4
F
J
(
x
¯
,
y
¯
)
|
,
lim
⁡
n
→
∞
∥
e
n
+
1
∥
∥
e
n
∥
=
|
λ
1,2
,
3,4
F
J
(
x
¯
,
y
¯
)
|
,
where λ
1,2
,
3,4
F
J
(
x
¯
,
y
¯
) are the characteristic roots of Jacobian matrix F
J
(
x
¯
,
y
¯
).


### 2.3. Examples

In order to verify our theoretical results and to support our theoretical discussions, we consider several interesting numerical examples in this section. These examples represent different types of qualitative behavior of solutions to the system of nonlinear difference equations (5). The first and last examples show that positive equilibrium of system (5) is unstable with suitable parametric choices. Moreover, from the remaining examples it is clear that unique positive equilibrium point of system (5) is globally asymptotically stable with different parametric values.Example 1.
Letα
1
=
6.5, β
1
=
8.5, γ
1
=
12.3, a
1
=
1.6, b
1
=
5.02, c
1
=
1.2, α
2
=
3.5, β
2
=
5.5, γ
2
=
1.3, a
2
=
0.1, b
2
=
0.003, and c
2
=
2.2. Then, system (5) can be written as
(50)
x
n
+
1
=
6.5
+
12.3
e
-
x
n
+
8.5
e
-
x
n
-
1
1.6
+
1.2
y
n
+
5.02
y
n
-
1
,
y
n
+
1
=
3.5
+
1.3
e
-
y
n
+
5.5
e
-
y
n
-
1
0.1
+
2.2
x
n
+
0.003
x
n
-
1
,
with initial conditions x
-
1
=
1.2
,
x
0
=
1.1
,
y
-
1
=
2.1
, and y
0
=
2.2.
In this case the positive equilibrium point of the system (50) is unstable. Moreover, in Figure 1 the plot of x
n is shown in Figure 1(a), the plot of y
n is shown in Figure 1(b), and a phase portrait of the system (50) is shown in Figure 1(c).Plots for the system (50).
(a)
Plot ofx
n for the system  (50)
(b)
Plot ofy
n for the system  (50)
(c)
Phase portrait of system  (50)Example 2.
Letα
1
=
1.5, β
1
=
7.5, γ
1
=
13.5, a
1
=
2.1, b
1
=
7.92, c
1
=
0.95, α
2
=
4.1, β
2
=
7.2, γ
2
=
16.7, a
2
=
0.01, b
2
=
0.2, and c
2
=
2.3. Then, system (5) can be written as
(51)
x
n
+
1
=
1.5
+
7.5
e
-
x
n
+
13.5
e
-
x
n
-
1
2.1
+
7.92
y
n
+
0.95
y
n
-
1
,
y
n
+
1
=
4.1
+
7.2
e
-
y
n
+
16.7
e
-
y
n
-
1
0.01
+
0.2
x
n
+
2.3
x
n
-
1
,
with initial conditions x
-
1
=
0.37
,
x
0
=
0.36
,
y
-
1
=
4.5, and y
0
=
4.6.
In this case the unique positive equilibrium point of the system (51) is given by (
x
¯
,
y
¯
)
=
(
0.370864,4.62495
). Moreover, in Figure 2 the plot of x
n is shown in Figure 2(a), the plot of y
n is shown in Figure 2(b), and an attractor of the system (51) is shown in Figure 2(c).Plots for the system (51).
(a)
Plot ofx
n for the system  (51)
(b)
Plot ofy
n for the system  (51)
(c)
An attractor of the system  (51)Example 3.
Letα
1
=
8.98, β
1
=
75, γ
1
=
135, a
1
=
21, b
1
=
79, c
1
=
5, α
2
=
41, β
2
=
71.9, γ
2
=
16, a
2
=
6, b
2
=
2, andc
2
=
30. Then, system (5) can be written as
(52)
x
n
+
1
=
8.98
+
75
e
-
x
n
+
135
e
-
x
n
-
1
21
+
79
y
n
+
5
y
n
-
1
,
y
n
+
1
=
41
+
71.9
e
-
y
n
+
16
e
-
y
n
-
1
6
+
2
x
n
+
30
x
n
-
1
,
with initial conditions x
-
1
=
0.634, x
0
=
0.63, y
-
1
=
1.99, and y
0
=
2.
In this case the unique positive equilibrium point of the system (52) is given by (
x
¯
,
y
¯
)
=
(
0.634497,2.00755
). Moreover, in Figure 3 the plot of x
n is shown in Figure 3(a), the plot of y
n is shown in Figure 3(b), and an attractor of the system (52) is shown in Figure 3(c).Plots for the system (52).
(a)
Plot ofx
n for the system  (52)
(b)
Plot ofy
n for the system  (52)
(c)
An attractor of the system  (52)Example 4.
Letα
1
=
23.8, β
1
=
350, γ
1
=
66, a
1
=
55, b
1
=
14, c
1
=
111, α
2
=
66, β
2
=
35, γ
2
=
87, a
2
=
1.7, b
2
=
80, and c
2
=
1.4. Then, system (5) can be written as
(53)
x
n
+
1
=
23.8
+
350
e
-
x
n
+
66
e
-
x
n
-
1
55
+
14
y
n
+
111
y
n
-
1
,
y
n
+
1
=
66
+
35
e
-
y
n
+
87
e
-
y
n
-
1
1.7
+
80
x
n
+
1.4
x
n
-
1
,
with initial conditions x
-
1
=
0.82
,
x
0
=
0.81
,
y
-
1
=
1.2, and y
0
=
1.3.
In this case the unique positive equilibrium point of the system (53) is given by (
x
¯
,
y
¯
)
=
(
0.883148,1.33372
). Moreover, in Figure 4 the plot of x
n is shown in Figure 4(a), the plot of y
n is shown in Figure 4(b), and an attractor of system (53) is shown in Figure 4(c).Plots for the system (53).
(a)
Plot ofx
n for the system  (53)
(b)
Plot ofy
n for the system  (53)
(c)
An attractor of the system  (53)Example 5.
Letα
1
=
8, β
1
=
3, γ
1
=
6, a
1
=
0.05, b
1
=
4, c
1
=
16, α
2
=
9, β
2
=
17, γ
2
=
7, a
2
=
0.07, b
2
=
22, and c
2
=
0.2. Then, system (5) can be written as
(54)
x
n
+
1
=
8
+
3
e
-
x
n
+
6
e
-
x
n
-
1
0.05
+
4
y
n
+
16
y
n
-
1
,
y
n
+
1
=
9
+
17
e
-
y
n
+
7
e
-
y
n
-
1
0.07
+
22
x
n
+
0.2
x
n
-
1
,
with initial conditions x
-
1
=
0.5, x
0
=
0.4, y
-
1
=
1.2, and y
0
=
1.3.
In this case the positive equilibrium point of the system (54) is unstable. Moreover, in Figure 5 the plot of x
n is shown in Figure 5(a), the plot of y
n is shown in Figure 5(b), and a phase portrait of system (54) is shown in Figure 5(c).Plots for the system (54).
(a)
Plot ofx
n for the system  (54)
(b)
Plot ofy
n for the system  (54)
(c)
Phase portrait of system  (54)


## 2.1. Stability Analysis

Let us consider four-dimensional discrete dynamical system of the form:(9)
x
n
+
1
=
f
(
x
n
,
x
n
-
1
,
y
n
,
y
n
-
1
)
y
n
+
1
=
g
(
x
n
,
x
n
-
1
,
y
n
,
y
n
-
1
)
,
n
=
0,1
,
…
,
where f
:
I
2
×
J
2
→
I and g
:
I
2
×
J
2
→
J are continuously differentiable functions and I
,
J are some intervals of real numbers. Furthermore, a solution {
(
x
n
,
y
n
)
}
n
=
-
1
∞ of system (9) is uniquely determined by initial conditions (
x
i
,
y
i
)
∈
I
×
J for i
∈
{
-
1,0
}. Along with system (9) we consider the corresponding vector map F
=
(
f
,
x
n
,
g
,
y
n
). An equilibrium point of (9) is a point (
x
¯
,
y
¯
) that satisfies
(10)
x
¯
=
f
(
x
¯
,
x
¯
,
y
¯
,
y
¯
)
,
y
¯
=
g
(
x
¯
,
x
¯
,
y
¯
,
y
¯
)
.
The point (
x
¯
,
y
¯
) is also called a fixed point of the vector map F.Definition 3.
Let(
x
¯
,
y
¯
) be an equilibrium point of the system (9).(i)
An equilibrium point(
x
¯
,
y
¯
) is said to be stable if for every ɛ
>
0 there exists δ
>
0 such that for every initial condition (
x
i
,
y
i
), i
∈
{
-
1,0
}, if ∥
∑
i
=
-
1
0
‍
(
x
i
,
y
i
)
-
(
x
¯
,
y
¯
)
∥
<
δ implies ∥
(
x
n
,
y
n
)
-
(
x
¯
,
y
¯
)
∥
<
ɛ for all n
>
0, where ∥
·
∥ is usual Euclidian norm in R
2.
(ii)
An equilibrium point(
x
¯
,
y
¯
) is said to be unstable if it is not stable.
(iii)
An equilibrium point(
x
¯
,
y
¯
) is said to be asymptotically stable if there exists η
>
0 such that ∥
∑
i
=
-
1
0
‍
(
x
i
,
y
i
)
-
(
x
¯
,
y
¯
)
∥
<
η and (
x
n
,
y
n
)
→
(
x
¯
,
y
¯
) as n
→
∞.
(iv)
An equilibrium point(
x
¯
,
y
¯
) is called global attractor if (
x
n
,
y
n
)
→
(
x
¯
,
y
¯
) as n
→
∞.
(v)
An equilibrium point(
x
¯
,
y
¯
) is called asymptotic global attractor if it is a global attractor and stable.Definition 4.
Let(
x
¯
,
y
¯
) be an equilibrium point of a map F
=
(
f
,
x
n
,
g
,
y
n
), where f and g are continuously differentiable functions at (
x
¯
,
y
¯
). The linearized system of (9) about the equilibrium point (
x
¯
,
y
¯
) is
(11)
X
n
+
1
=
F
(
X
n
)
=
F
J
X
n
,
where X
n
=
(
x
n
y
n
x
n
-
1
y
n
-
1
) and F
J is Jacobian matrix of system (9) about the equilibrium point (
x
¯
,
y
¯
).To construct corresponding linearized form of system (5) we consider the following transformation:
(12)
(
x
n
,
y
n
,
x
n
-
1
,
y
n
-
1
)
⟼
(
f
,
g
,
f
1
,
g
1
)
,
where f
=
x
n
+
1, g
=
y
n
+
1, f
1
=
x
n, and g
1
=
y
n. The linearized system of (5) about (
x
¯
,
y
¯
) is given by
(13)
Z
n
+
1
=
F
J
(
x
¯
,
y
¯
)
Z
n
,
where Z
n
=
(
x
n
y
n
x
n
-
1
y
n
-
1
) and the Jacobian matrix about the fixed point (
x
¯
,
y
¯
) under transformation (12) is given by
(14)
F
J
(
x
¯
,
y
¯
)
=
(
A
1
A
2
A
3
A
4
B
1
B
2
B
3
B
4
1
0
0
0
0
1
0
0
)
,
where
(15)
A
1
=
-
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
,
A
2
=
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
,
A
3
=
-
γ
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
,
A
4
=
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
,
B
1
=
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
,
B
2
=
-
β
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
,
B
3
=
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
,
B
4
=
-
γ
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
.Lemma 5 (see [9]).
Assume thatX
n
+
1
=
F
(
X
n
), n
=
0,1
,
…, is a system of difference equations such that X
¯ is a fixed point of F. If all eigenvalues of the Jacobian matrix J
F about X
¯ lie inside the open unit disk |
λ
|
<
1, then X
¯ is locally asymptotically stable. If one of them has a modulus greater than one, then X
¯ is unstable.Theorem 6.
System (5) has a unique positive equilibrium point (
x
¯
,
y
¯
)
∈
[
L
1
,
U
1
]
×
[
L
2
,
U
2
], if the following condition is satisfied:
(16)
a
2
+
L
1
(
b
2
+
c
2
)
<
α
2
+
(
β
2
+
γ
2
)
e
-
K
K
,
where
(17)
K
=
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
L
1
(
b
1
+
c
1
)
.Proof.
Consider the following system of equations:(18)
x
=
α
1
+
(
β
1
+
γ
1
)
e
-
x
a
1
+
(
b
1
+
c
1
)
y
,
y
=
α
2
+
(
β
2
+
γ
2
)
e
-
y
a
2
+
(
b
2
+
c
2
)
x
.
From (18), it follows that
(19)
y
=
α
1
+
(
β
1
+
γ
1
)
e
-
x
-
a
1
x
(
b
1
+
c
1
)
x
,
x
=
α
2
+
(
β
2
+
γ
2
)
e
-
y
-
a
2
y
(
b
2
+
c
2
)
y
.
Set
(20)
F
(
x
)
=
α
2
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
-
a
2
f
(
x
)
(
b
2
+
c
2
)
f
(
x
)
-
x
,
where f
(
x
)
=
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
-
a
1
x
)
/
(
(
b
1
+
c
1
)
x
) and x
∈
[
L
1
,
U
1
]. Then, it follows that
(21)
F
(
L
1
)
=
L
1
(
b
1
+
c
1
)
×
(
a
2
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
L
1
(
b
1
+
c
1
)
α
2
+
(
β
2
+
γ
2
)
e
-
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
/
L
1
(
b
1
+
c
1
)
-
a
2
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
L
1
(
b
1
+
c
1
)
)
×
(
(
b
2
+
c
2
)
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
)
-
1
-
L
1
.
Furthermore, it is easy to see that
(22)
(
a
1
+
(
b
1
+
c
1
)
U
2
)
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
)
>
a
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
U
1
)
.
From (22) it follows that
(23)
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
>
0
.
Let K
=
(
α
1
+
e
-
L
1
(
β
1
+
γ
1
)
-
a
1
L
1
)
/
(
L
1
(
b
1
+
c
1
)
)
>
0. Then F
(
L
1
) can be expressed as
(24)
F
(
L
1
)
=
α
2
+
(
β
2
+
γ
2
)
e
-
K
-
a
2
K
(
b
2
+
c
2
)
K
-
L
1
.
Suppose that a
2
+
L
1
(
b
2
+
c
2
)
<
(
α
2
+
(
β
2
+
γ
2
)
e
-
K
)
/
K; then, it follows that F
(
L
1
)
>
0. Furthermore, we have
(25)
F
(
U
1
)
=
U
1
(
b
1
+
c
1
)
(
a
2
(
β
1
+
γ
1
)
(
1
-
e
-
U
1
)
U
1
(
b
1
+
c
1
)
(
β
2
+
γ
2
)
e
(
β
1
+
γ
1
)
(
1
-
e
-
U
1
)
/
U
1
(
b
1
+
c
1
)
+
α
2
+
a
2
(
β
1
+
γ
1
)
(
1
-
e
-
U
1
)
U
1
(
b
1
+
c
1
)
)
×
(
(
b
2
+
c
2
)
(
α
1
+
e
-
U
1
(
β
1
+
γ
1
)
-
a
1
U
1
)
)
-
1
-
U
1
.
Then it is easy to see that
(26)
α
1
+
e
-
U
1
(
β
1
+
γ
1
)
-
a
1
U
1
=
(
β
1
+
γ
1
)
(
e
-
U
1
-
1
)
<
0
,
which gives that F
(
U
1
)
<
0. Hence, F
(
x
)
=
0 has at least one positive solution in [
L
1
,
U
1
]. Moreover, we obtain that
(27)
F
′
(
x
)
=
-
1
-
f
′
(
x
)
(
b
2
+
c
2
)
f
(
x
)
×
[
a
2
+
(
b
2
+
c
2
)
x
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
]
,
where
(28)
f
′
(
x
)
=
-
α
1
+
(
β
1
+
γ
1
)
(
1
+
x
)
e
-
x
(
b
1
+
c
1
)
x
2
.
Then, from (28) it follows that f
(
U
1
)
<
f
(
x
)
<
f
(
L
1
) and using (26) we obtain
(29)
F
′
(
x
)
=
-
1
+
α
1
+
(
β
1
+
γ
1
)
(
1
+
x
)
e
-
x
(
b
1
+
c
1
)
(
b
2
+
c
2
)
x
2
f
(
x
)
×
[
a
2
+
(
b
2
+
c
2
)
x
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
]
≤
-
1
+
α
1
+
(
β
1
+
γ
1
)
(
1
+
x
)
e
-
x
(
b
1
+
c
1
)
(
b
2
+
c
2
)
x
2
f
(
U
1
)
×
[
a
2
+
(
b
2
+
c
2
)
x
+
(
β
2
+
γ
2
)
e
-
f
(
x
)
]
<
0
.
Hence, F
(
x
)
=
0 has a unique positive solution in [
L
1
,
U
1
]. The proof is therefore completed.Theorem 7.
The unique positive equilibrium point of system (5) is locally asymptotically stable under the following condition:
(30)
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
+
(
1
+
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
×
(
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
+
(
β
2
+
γ
2
)
e
-
L
2
a
2
+
(
b
2
+
c
2
)
L
1
)
+
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
×
(
(
b
1
+
c
1
)
(
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
)
(
a
1
+
(
b
1
+
c
1
)
L
2
)
2
+
γ
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
<
1
.Proof.
The characteristic polynomial of Jacobian matrixF
J
(
x
¯
,
y
¯
) about (
x
¯
,
y
¯
) is given by
(31)
P
(
λ
)
=
λ
4
-
(
A
1
+
B
2
)
λ
3
-
(
A
2
B
1
-
A
1
B
2
+
B
3
+
B
4
)
λ
2
-
(
(
A
3
+
A
4
)
B
1
-
(
B
3
+
B
4
)
A
1
)
λ
.
Clearly, one root of P
(
λ
) is 0. To check the behavior of the other three roots of P
(
λ
), we let Φ
(
λ
)
=
λ
3 and Ψ
(
λ
)
=
(
A
1
+
B
2
)
λ
2
+
(
A
2
B
1
-
A
1
B
2
+
B
3
+
B
4
)
λ
+
(
A
3
+
A
4
)
B
1
-
(
B
3
+
B
4
)
A
1. Assume that (30) holds and |
λ
|
=
1; then, one has
(32)
|
Ψ
(
λ
)
|
<
|
A
1
+
B
2
|
+
|
A
2
B
1
-
A
1
B
2
+
B
3
+
B
4
|
+
|
(
A
3
+
A
4
)
B
1
-
(
B
3
+
B
4
)
A
1
|
<
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
+
(
1
+
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
)
×
(
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
+
(
β
2
+
γ
2
)
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
)
+
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
×
(
(
b
1
+
c
1
)
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
+
γ
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
)
<
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
+
(
1
+
β
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
×
(
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
+
(
β
2
+
γ
2
)
e
-
L
2
a
2
+
(
b
2
+
c
2
)
L
1
)
+
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
(
a
2
+
(
b
2
+
c
2
)
L
1
)
2
×
(
(
b
1
+
c
1
)
(
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
)
(
a
1
+
(
b
1
+
c
1
)
L
2
)
2
+
γ
1
e
-
L
1
a
1
+
(
b
1
+
c
1
)
L
2
)
<
1
.
Then, by Rouche’s Theorem, Φ
(
λ
) and Φ
(
λ
)
-
Ψ
(
λ
) have the same number of zeroes in an open unit disk |
λ
|
<
1. Hence, all the roots of (31) satisfy |
λ
|
<
1, and it follows from Lemma 5 that the unique positive equilibrium point (
x
¯
,
y
¯
) of the system (5) is locally asymptotically stable.Theorem 8.
The unique positive equilibrium point(
x
¯
,
y
¯
) of system (5) is globally asymptotically stable, if the following condition is satisfied:
(33)
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
<
x
¯
(
a
1
+
(
b
1
+
c
1
)
L
2
)
,
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
<
y
¯
(
a
2
+
(
b
2
+
c
2
)
L
1
)
.Proof.
Arguing as in [18], we consider the following discrete time analogue of Lyapunov function:
(34)
V
n
=
x
¯
(
x
n
x
¯
-
1
-
ln
⁡
x
n
x
¯
)
+
y
¯
(
y
n
y
¯
-
1
-
ln
⁡
y
n
y
¯
)
.
Then nonnegativity of V
n follows from the following inequality:
(35)
x
-
1
-
ln
⁡
⁡
x
≥
0
∀
x
>
0
.
Furthermore, we have
(36)
-
ln
⁡
⁡
(
x
n
+
1
x
n
)
=
ln
⁡
⁡
(
1
-
(
1
-
x
n
x
n
+
1
)
)
≤
-
x
n
+
1
-
x
n
x
n
+
1
,
-
ln
⁡
⁡
(
y
n
+
1
y
n
)
=
ln
⁡
⁡
(
1
-
(
1
-
y
n
y
n
+
1
)
)
≤
-
y
n
+
1
-
y
n
y
n
+
1
.
Assume that (33) holds true; then, it follows that
(37)
V
n
+
1
-
V
n
=
x
¯
(
x
n
+
1
x
¯
-
1
-
ln
⁡
x
n
+
1
x
¯
)
+
y
¯
(
y
n
+
1
y
¯
-
1
-
ln
⁡
y
n
+
1
y
¯
)
-
x
¯
(
x
n
x
¯
-
1
-
ln
⁡
x
n
x
¯
)
-
y
¯
(
y
n
y
¯
-
1
-
ln
⁡
y
n
y
¯
)
≤
(
x
n
+
1
-
x
n
)
+
(
y
n
+
1
-
y
n
)
-
x
¯
x
n
+
1
(
x
n
+
1
-
x
n
)
-
y
¯
y
n
+
1
(
y
n
+
1
-
y
n
)
=
(
x
n
+
1
-
x
n
)
(
1
-
x
¯
x
n
+
1
)
+
(
y
n
+
1
-
y
n
)
(
1
-
y
¯
y
n
+
1
)
≤
(
U
1
-
L
1
)
×
(
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
-
x
¯
(
a
1
+
(
b
1
+
c
1
)
L
2
)
α
1
+
(
β
1
+
γ
1
)
e
-
L
1
)
+
(
U
2
-
L
2
)
×
(
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
-
y
¯
(
a
2
+
(
b
2
+
c
2
)
L
1
)
α
2
+
(
β
2
+
γ
2
)
e
-
L
2
)
≤
0
for all n
≥
0 so that V
n
≥
0 is monotonically decreasing sequence. It follows that lim
⁡
n
→
∞
V
n
≥
0. Hence, we obtain that
(38)
lim
⁡
n
→
∞
(
V
n
+
1
-
V
n
)
=
0
.
Then it follows that lim
⁡
n
→
∞
x
n
+
1
=
x
¯ and lim
⁡
n
→
∞
y
n
+
1
=
y
¯. Furthermore, V
n
≤
V
0 for all n
≥
0, which gives that (
x
¯
,
y
¯
)
∈
[
L
1
,
U
1
]
×
[
L
2
,
U
2
] is uniformly stable. Hence, unique positive equilibrium point (
x
¯
,
y
¯
)
∈
[
L
1
,
U
1
]
×
[
L
2
,
U
2
] of system (5) is globally asymptotically stable.


## 2.2. Rate of Convergence

In this section we will determine the rate of convergence of a solution that converges to the unique positive equilibrium point of the system (5).The following result gives the rate of convergence of solutions of a system of difference equations:(39)
X
n
+
1
=
(
A
+
B
(
n
)
)
X
n
,
where X
n is an m-dimensional vector, A
∈
C
m
×
m is a constant matrix, and B
:
Z
+
→
C
m
×
m is a matrix function satisfying
(40)
∥
B
(
n
)
∥
⟶
0
as n
→
∞, where ∥
·
∥ denotes any matrix norm which is associated with the vector norm
(41)
∥
(
x
,
y
)
∥
=
x
2
+
y
2
.Proposition 9 ((Perron’s Theorem) [19]).
Suppose that condition (40) holds. If X
n is a solution of (39), then either X
n
=
0 for all large n or
(42)
ρ
=
lim
⁡
n
→
∞
(
∥
X
n
∥
)
1
/
n
exists and is equal to the modulus of one of the eigenvalues of matrix A.Proposition 10 (see [19]).
Suppose that condition (40) holds. If X
n is a solution of (39), then either X
n
=
0 for all large n or
(43)
ρ
=
lim
⁡
n
→
∞
∥
X
n
+
1
∥
∥
X
n
∥
exists and is equal to the modulus of one of the eigenvalues of matrix A.Let{
(
x
n
,
y
n
)
} be any solution of the system (5) such that lim
⁡
n
→
∞
x
n
=
x
¯ and lim
⁡
n
→
∞
y
n
=
y
¯, where x
¯
∈
[
L
1
,
U
1
] and y
¯
∈
[
L
2
,
U
2
]. To find the error terms, one has from the system (5)
(44)
x
n
+
1
-
x
¯
=
α
1
+
β
1
e
-
x
n
+
γ
1
e
-
x
n
-
1
a
1
+
b
1
y
n
+
c
1
y
n
-
1
-
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
=
β
1
(
e
-
x
n
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
x
¯
)
(
x
n
-
x
¯
)
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
y
n
-
y
¯
)
+
γ
1
(
e
-
x
n
-
1
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
1
-
x
¯
)
(
x
n
-
1
-
x
¯
)
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
y
n
-
1
-
y
¯
)
,
y
n
+
1
-
y
¯
=
α
2
+
β
2
e
-
y
n
+
γ
2
e
-
y
n
-
1
a
2
+
b
2
x
n
+
c
2
x
n
-
1
-
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
=
β
2
(
e
-
y
n
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
x
¯
)
(
y
n
-
y
¯
)
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
x
n
-
x
¯
)
+
γ
2
(
e
-
y
n
-
1
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
1
-
y
¯
)
(
y
n
-
1
-
y
¯
)
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
x
n
-
1
-
x
¯
)
.
Let e
n
1
=
x
n
-
x
¯ and e
n
2
=
y
n
-
y
¯; then, one has
(45)
e
n
+
1
1
=
a
n
e
n
1
+
b
n
e
n
2
+
c
n
e
n
-
1
1
+
d
n
e
n
-
1
2
,
e
n
+
1
2
=
f
n
e
n
1
+
g
n
e
n
2
+
h
n
e
n
-
1
1
+
k
n
e
n
-
1
2
,
where
(46)
a
n
=
β
1
(
e
-
x
n
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
x
¯
)
,
b
n
=
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
,
c
n
=
γ
1
(
e
-
x
n
-
1
-
e
-
x
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
(
x
n
-
1
-
x
¯
)
,
d
n
=
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
(
a
1
+
b
1
y
n
+
c
1
y
n
-
1
)
,
f
n
=
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
,
g
n
=
β
2
(
e
-
y
n
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
x
¯
)
,
h
n
=
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
,
k
n
=
γ
2
(
e
-
y
n
-
1
-
e
-
y
¯
)
(
a
2
+
b
2
x
n
+
c
2
x
n
-
1
)
(
y
n
-
1
-
y
¯
)
.
Moreover,
(47)
lim
⁡
n
→
∞
a
n
=
-
β
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
=
A
1
,
lim
⁡
n
→
∞
b
n
=
-
b
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
=
A
2
,
lim
⁡
n
→
∞
c
n
=
-
γ
1
e
-
x
¯
a
1
+
(
b
1
+
c
1
)
y
¯
=
A
3
,
lim
⁡
n
→
∞
d
n
=
-
c
1
(
α
1
+
(
β
1
+
γ
1
)
e
-
x
¯
)
(
a
1
+
(
b
1
+
c
1
)
y
¯
)
2
=
A
4
,
lim
⁡
n
→
∞
f
n
=
-
b
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
=
B
1
,
lim
⁡
n
→
∞
g
n
=
-
β
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
=
B
2
,
lim
⁡
n
→
∞
h
n
=
-
c
2
(
α
2
+
(
β
2
+
γ
2
)
e
-
y
¯
)
(
a
2
+
(
b
2
+
c
2
)
x
¯
)
2
=
B
3
,
lim
⁡
n
→
∞
k
n
=
-
γ
2
e
-
y
¯
a
2
+
(
b
2
+
c
2
)
x
¯
=
B
4
.
Now the limiting system of error terms can be written as
(48)
[
e
n
+
1
1
e
n
+
1
2
e
n
1
e
n
2
]
=
[
A
1
A
2
A
3
A
4
B
1
B
2
B
3
B
4
1
0
0
0
0
1
0
0
]
[
e
n
1
e
n
2
e
n
-
1
1
e
n
-
1
2
]
,
which is similar to linearized system of (5) about the equilibrium point (
x
¯
,
y
¯
).Using Proposition9, one has the following result.Theorem 11.
Assume that{
(
x
n
,
y
n
)
} is a positive solution of the system (5) such that lim
⁡
n
→
∞
x
n
=
x
¯, and lim
⁡
n
→
∞
y
n
=
y
¯, where x
¯
∈
[
L
1
,
U
1
] and y
¯
∈
[
L
2
,
U
2
]. Then, the error vector e
n
=
(
e
n
1
e
n
2
e
n
-
1
1
e
n
-
1
2
) of every solution of (5) satisfies both of the following asymptotic relations:
(49)
lim
⁡
n
→
∞
(
∥
e
n
∥
)
1
/
n
=
|
λ
1,2
,
3,4
F
J
(
x
¯
,
y
¯
)
|
,
lim
⁡
n
→
∞
∥
e
n
+
1
∥
∥
e
n
∥
=
|
λ
1,2
,
3,4
F
J
(
x
¯
,
y
¯
)
|
,
where λ
1,2
,
3,4
F
J
(
x
¯
,
y
¯
) are the characteristic roots of Jacobian matrix F
J
(
x
¯
,
y
¯
).


## 2.3. Examples

In order to verify our theoretical results and to support our theoretical discussions, we consider several interesting numerical examples in this section. These examples represent different types of qualitative behavior of solutions to the system of nonlinear difference equations (5). The first and last examples show that positive equilibrium of system (5) is unstable with suitable parametric choices. Moreover, from the remaining examples it is clear that unique positive equilibrium point of system (5) is globally asymptotically stable with different parametric values.Example 1.
Letα
1
=
6.5, β
1
=
8.5, γ
1
=
12.3, a
1
=
1.6, b
1
=
5.02, c
1
=
1.2, α
2
=
3.5, β
2
=
5.5, γ
2
=
1.3, a
2
=
0.1, b
2
=
0.003, and c
2
=
2.2. Then, system (5) can be written as
(50)
x
n
+
1
=
6.5
+
12.3
e
-
x
n
+
8.5
e
-
x
n
-
1
1.6
+
1.2
y
n
+
5.02
y
n
-
1
,
y
n
+
1
=
3.5
+
1.3
e
-
y
n
+
5.5
e
-
y
n
-
1
0.1
+
2.2
x
n
+
0.003
x
n
-
1
,
with initial conditions x
-
1
=
1.2
,
x
0
=
1.1
,
y
-
1
=
2.1
, and y
0
=
2.2.
In this case the positive equilibrium point of the system (50) is unstable. Moreover, in Figure 1 the plot of x
n is shown in Figure 1(a), the plot of y
n is shown in Figure 1(b), and a phase portrait of the system (50) is shown in Figure 1(c).Plots for the system (50).
(a)
Plot ofx
n for the system  (50)
(b)
Plot ofy
n for the system  (50)
(c)
Phase portrait of system  (50)Example 2.
Letα
1
=
1.5, β
1
=
7.5, γ
1
=
13.5, a
1
=
2.1, b
1
=
7.92, c
1
=
0.95, α
2
=
4.1, β
2
=
7.2, γ
2
=
16.7, a
2
=
0.01, b
2
=
0.2, and c
2
=
2.3. Then, system (5) can be written as
(51)
x
n
+
1
=
1.5
+
7.5
e
-
x
n
+
13.5
e
-
x
n
-
1
2.1
+
7.92
y
n
+
0.95
y
n
-
1
,
y
n
+
1
=
4.1
+
7.2
e
-
y
n
+
16.7
e
-
y
n
-
1
0.01
+
0.2
x
n
+
2.3
x
n
-
1
,
with initial conditions x
-
1
=
0.37
,
x
0
=
0.36
,
y
-
1
=
4.5, and y
0
=
4.6.
In this case the unique positive equilibrium point of the system (51) is given by (
x
¯
,
y
¯
)
=
(
0.370864,4.62495
). Moreover, in Figure 2 the plot of x
n is shown in Figure 2(a), the plot of y
n is shown in Figure 2(b), and an attractor of the system (51) is shown in Figure 2(c).Plots for the system (51).
(a)
Plot ofx
n for the system  (51)
(b)
Plot ofy
n for the system  (51)
(c)
An attractor of the system  (51)Example 3.
Letα
1
=
8.98, β
1
=
75, γ
1
=
135, a
1
=
21, b
1
=
79, c
1
=
5, α
2
=
41, β
2
=
71.9, γ
2
=
16, a
2
=
6, b
2
=
2, andc
2
=
30. Then, system (5) can be written as
(52)
x
n
+
1
=
8.98
+
75
e
-
x
n
+
135
e
-
x
n
-
1
21
+
79
y
n
+
5
y
n
-
1
,
y
n
+
1
=
41
+
71.9
e
-
y
n
+
16
e
-
y
n
-
1
6
+
2
x
n
+
30
x
n
-
1
,
with initial conditions x
-
1
=
0.634, x
0
=
0.63, y
-
1
=
1.99, and y
0
=
2.
In this case the unique positive equilibrium point of the system (52) is given by (
x
¯
,
y
¯
)
=
(
0.634497,2.00755
). Moreover, in Figure 3 the plot of x
n is shown in Figure 3(a), the plot of y
n is shown in Figure 3(b), and an attractor of the system (52) is shown in Figure 3(c).Plots for the system (52).
(a)
Plot ofx
n for the system  (52)
(b)
Plot ofy
n for the system  (52)
(c)
An attractor of the system  (52)Example 4.
Letα
1
=
23.8, β
1
=
350, γ
1
=
66, a
1
=
55, b
1
=
14, c
1
=
111, α
2
=
66, β
2
=
35, γ
2
=
87, a
2
=
1.7, b
2
=
80, and c
2
=
1.4. Then, system (5) can be written as
(53)
x
n
+
1
=
23.8
+
350
e
-
x
n
+
66
e
-
x
n
-
1
55
+
14
y
n
+
111
y
n
-
1
,
y
n
+
1
=
66
+
35
e
-
y
n
+
87
e
-
y
n
-
1
1.7
+
80
x
n
+
1.4
x
n
-
1
,
with initial conditions x
-
1
=
0.82
,
x
0
=
0.81
,
y
-
1
=
1.2, and y
0
=
1.3.
In this case the unique positive equilibrium point of the system (53) is given by (
x
¯
,
y
¯
)
=
(
0.883148,1.33372
). Moreover, in Figure 4 the plot of x
n is shown in Figure 4(a), the plot of y
n is shown in Figure 4(b), and an attractor of system (53) is shown in Figure 4(c).Plots for the system (53).
(a)
Plot ofx
n for the system  (53)
(b)
Plot ofy
n for the system  (53)
(c)
An attractor of the system  (53)Example 5.
Letα
1
=
8, β
1
=
3, γ
1
=
6, a
1
=
0.05, b
1
=
4, c
1
=
16, α
2
=
9, β
2
=
17, γ
2
=
7, a
2
=
0.07, b
2
=
22, and c
2
=
0.2. Then, system (5) can be written as
(54)
x
n
+
1
=
8
+
3
e
-
x
n
+
6
e
-
x
n
-
1
0.05
+
4
y
n
+
16
y
n
-
1
,
y
n
+
1
=
9
+
17
e
-
y
n
+
7
e
-
y
n
-
1
0.07
+
22
x
n
+
0.2
x
n
-
1
,
with initial conditions x
-
1
=
0.5, x
0
=
0.4, y
-
1
=
1.2, and y
0
=
1.3.
In this case the positive equilibrium point of the system (54) is unstable. Moreover, in Figure 5 the plot of x
n is shown in Figure 5(a), the plot of y
n is shown in Figure 5(b), and a phase portrait of system (54) is shown in Figure 5(c).Plots for the system (54).
(a)
Plot ofx
n for the system  (54)
(b)
Plot ofy
n for the system  (54)
(c)
Phase portrait of system  (54)


## 3. Concluding Remarks

In literature several articles are related to qualitative behavior of exponential systems of rational difference equations. It is a very interesting mathematical problem to study the dynamics of such systems because these are closely related to models in population dynamics and biological sciences. This work is related to qualitative behavior of an exponential system of second-order rational difference equations. We have investigated the existence and uniqueness of positive steady state of system (5). Under certain parametric conditions the boundedness and persistence of positive solutions are proved. Moreover, we have shown that unique positive equilibrium point of system (5) is locally as well as globally asymptotically stable. The main objective of dynamical systems theory is to predict the global behavior of a system based on the knowledge of its present state. An approach to this problem consists of determining the possible global behaviors of the system and determining which parametric conditions lead to these long-term behaviors. By constructing a discrete Lyapunov function, we have obtained the global asymptotic stability of the positive equilibrium of (5). Furthermore, rate of convergence of positive solutions of (5) which converge to its unique positive equilibrium point is demonstrated. Finally, some illustrative examples are provided to support our theoretical discussion.



---

*Source: 375890-2014-07-24.xml*