# Global Character of a Six-Dimensional Nonlinear System of Difference Equations

Authors: Mehmet Gümüş; Yüksel Soykan
Journal: Discrete Dynamics in Nature and Society
 (2016)

---

## Abstract

The aim of this paper is to study the dynamical behavior of positive solutions for a system of rational difference equations of the following form:u
n
+
1
=
α
u
n
-
1
/
β
+
γ
v
n
-
2
p, v
n
+
1
=
α
1
v
n
-
1
/
β
1
+
γ
1
u
n
-
2
p, n
=
0,1
,
…, where the parameters α
,
β
,
γ
,
α
1
,
β
1
,
γ
1
,
p and the initial values u
-
i
,
v
-
i for i
=
0,1
,
2 are positive real numbers.


---

## Body

## 1. Introduction

Higher-order fractional difference equations are of great importance in applications. Such equations also seem surely as discrete analogues and as numerical solutions of differential and delay differential equations which model various diverse phenomena in biology, ecology, physiology, physics, engineering, economics, and so on. For example, see [1, 2]. The theory of difference equations is the focus point in applicable analysis. That is, the importance of the theory of difference equations in mathematics as a whole will go on. Hence, it is very valuable to investigate the behavior of solutions of a system of fractional difference equations and to present the stability character of equilibrium points.In the recent times, the behavior of solutions of various systems of rational difference equations has been one of the main topics in the theory of difference equations (see [3–7] and the references cited therein).According to us, it is of great importance to investigate not only rational nonlinear difference equations and their systems, but also those equations and systems which contain powers of arbitrary positive numbers (see [8–12]).In [13], Kurbanlı et al. studied the behavior of positive solutions of the following system of difference equations:(1)
x
n
+
1
=
x
n
-
1
y
n
x
n
-
1
+
1
,
y
n
+
1
=
y
n
-
1
x
n
y
n
-
1
+
1
,
n
∈
N
,where the initial conditions are arbitrary nonnegative real numbers.In [14], Papaschinopoulos and Schinas studied the system of two nonlinear difference equations:(2)
x
n
+
1
=
A
+
y
n
x
n
-
p
,
y
n
+
1
=
A
+
x
n
y
n
-
q
,
n
∈
N
,where A
∈
(
0
,
∞
) and p
,
q are positive integers.In [15], Zhang et al. studied the dynamical behavior of positive solutions for a system for third-order rational difference equations:(3)
x
n
+
1
=
x
n
-
2
B
+
y
n
y
n
-
1
y
n
-
2
,
y
n
+
1
=
y
n
-
2
A
+
x
n
x
n
-
1
x
n
-
2
,
n
∈
N
.In [16], Din et al. studied the dynamics of a system of fourth-order rational difference equations:(4)
x
n
+
1
=
α
x
n
-
3
β
+
γ
y
n
y
n
-
1
y
n
-
2
y
n
-
3
,
y
n
+
1
=
α
1
y
n
-
3
β
1
+
γ
1
x
n
x
n
-
1
x
n
-
2
x
n
-
3
,
n
∈
N
.In [17], Touafek and Elsayed investigated the behavior of solutions of systems of difference equations:(5)
x
n
+
1
=
x
n
-
3
±
1
±
x
n
-
3
y
n
-
1
,
y
n
+
1
=
y
n
-
3
±
1
±
y
n
-
3
x
n
-
1
,
n
∈
N
,with a nonzero real number’s initial conditions.In [18], El-Owaidy et al. investigated the global behavior of the following difference equation:(6)
x
n
+
1
=
α
x
n
-
1
β
+
γ
x
n
-
2
p
,
n
∈
N
,with nonnegative parameters and nonnegative initial values.Motivated by all above-mentioned works, in the paper we investigate the equilibrium points, the local asymptotic stability of these points, the global behavior of positive solutions, the existence unbounded solutions, and the existence of the prime two-periodic solutions of the following system:(7)
u
n
+
1
=
α
u
n
-
1
β
+
γ
v
n
-
2
p
,
v
n
+
1
=
α
1
v
n
-
1
β
1
+
γ
1
u
n
-
2
p
,
n
=
0,1
,
…
,where the parameters α
,
β
,
γ
,
α
1
,
β
1
,
γ
1
,
p and the initial values u
-
i
,
v
-
i for i
=
0,1
,
2 are positive real numbers. Our results extend and complement some results in the literature.Note that system (7) can be reduced to the following system of difference equations:(8)
x
n
+
1
=
r
x
n
-
1
1
+
y
n
-
2
p
,
y
n
+
1
=
s
y
n
-
1
1
+
x
n
-
2
p
,
n
=
0,1
,
…by the change of variables u
n
=
β
1
/
γ
1
1
/
p
x
n and v
n
=
β
/
γ
1
/
p
y
n with r
=
α
/
β and s
=
α
1
/
β
1. So, in order to study system (7), we investigate system (8).By using the induction and the equations in (8) we see that if x
-
i
,
y
-
i are positive real numbers for i
∈
{
0,1
,
2
}, then(9)
min
⁡
x
n
,
y
n
>
0
,
n
≥
-
2
,which means that positive initial values generate positive solutions of system (8).As far as we examine, there is no paper dealing with system (7). Therefore, in this paper we focus on system (7) in order to fill in the gap.


## 2. Preliminaries

For the completeness in the paper, we find it useful to remember some basic concepts of the difference equations theory as follows.Let us introduce the six-dimensional discrete dynamical system:(10)
x
n
+
1
=
f
1
x
n
,
x
n
-
1
,
x
n
-
2
,
y
n
,
y
n
-
1
,
y
n
-
2
,
y
n
+
1
=
f
2
x
n
,
x
n
-
1
,
x
n
-
2
,
y
n
,
y
n
-
1
,
y
n
-
2
,
n
∈
N, where f
1
:
I
1
3
×
I
2
3
→
I
1 and f
2
:
I
1
3
×
I
2
3
→
I
2 are continuously differentiable functions and I
1
,
I
2 are some intervals of real numbers. Also, a solution {
(
x
n
,
y
n
)
}
n
=
-
2
∞ of system (10) is uniquely determined by initial values (
x
-
i
,
y
-
i
)
∈
I
1
×
I
2 for i
∈
{
0,1
,
2
}.Definition 1.
An equilibrium point of system (10) is a point (
x
¯
,
y
¯
) that satisfies (11)
x
¯
=
f
1
x
¯
,
x
¯
,
x
¯
,
y
¯
,
y
¯
,
y
¯
,
y
¯
=
f
2
x
¯
,
x
¯
,
x
¯
,
y
¯
,
y
¯
,
y
¯
.Together with system (10), if we consider the associated vector map F
=
(
f
1
,
x
n
,
x
n
-
1
,
x
n
-
2
,
f
2
,
y
n
,
y
n
-
1
,
y
n
-
2
), then the point (
x
¯
,
y
¯
) is also called a fixed point of the vector map F
.Definition 2.
Let(
x
¯
,
y
¯
) be an equilibrium point of system (10).(i)
An equilibrium point(
x
¯
,
y
¯
) is said to be stable if, for every ε
>
0, there exists δ
>
0 such that for every initial value (
x
-
i
,
y
-
i
)
∈
I
1
×
I
2, with ∑
i
=
-
2
0
x
i
-
x
¯
<
δ, ∑
i
=
-
2
0
y
i
-
y
¯
<
δ, implying x
n
-
x
¯
<
ε, y
n
-
y
¯
<
ε for n
∈
N.
(ii)
If an equilibrium point(
x
¯
,
y
¯
) is not stable, then it is said to be unstable.
(iii)
If an equilibrium point(
x
¯
,
y
¯
) is stable and there exists γ
>
0 such that ∑
i
=
-
2
0
x
i
-
x
¯
<
γ, ∑
i
=
-
2
0
y
i
-
y
¯
<
γ and (
x
n
,
y
n
)
→
(
x
¯
,
y
¯
) as n
→
∞, then it is said to be asymptotically stable.
(iv)
If(
x
n
,
y
n
)
→
(
x
¯
,
y
¯
) as n
→
∞, then an equilibrium point (
x
¯
,
y
¯
) is said to be a global attractor.
(v)
If an equilibrium point(
x
¯
,
y
¯
) is both global attractor and stable, then it is said to be globally asymptotically stable.Definition 3.
Let(
x
¯
,
y
¯
) be an equilibrium point of the map F where f
1 and f
2 are continuously differentiable functions at (
x
¯
,
y
¯
)
. The linearized system of (10) about the equilibrium point (
x
¯
,
y
¯
) is (12)
X
n
+
1
=
F
X
n
=
B
X
n
,where(13)
X
n
=
x
n
x
n
-
1
x
n
-
2
y
n
y
n
-
1
y
n
-
2and B is a Jacobian matrix of system (10) about the equilibrium point (
x
¯
,
y
¯
)
.Definition 4.
For the systemX
n
+
1
=
F
(
X
n
)
,
n
=
0,1
,
…, of difference equations such that X
¯ is a fixed point of F, if no eigenvalues of the Jacobian matrix B about X
¯ have absolute value equal to one, then X
¯ is called hyperbolic. If there exists an eigenvalue of the Jacobian matrix B about X
¯ with absolute value equal to one, then X
¯ is called nonhyperbolic.The following result, known as the Linearized Stability Theorem, is very practical in confirming the local stability character of the equilibrium point(
x
¯
,
y
¯
) of system (10).Theorem 5.
For the systemX
n
+
1
=
F
(
X
n
)
,
n
=
0,1
,
…, of difference equations such that X
¯ is a fixed point of F, if all eigenvalues of the Jacobian matrix B 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.For other basic knowledge about difference equations and their systems, the reader is referred to books [19–22].


## 3. Main Results

In this section we prove our main results.Theorem 6.
One has the following cases for the equilibrium points of (8):(i)
(
x
¯
0
,
y
¯
0
)
=
(
0,0
) is always the equilibrium point of system (8).
(ii)
Ifr
>
1 and s
>
1, then system (8) has the equilibrium point (
x
¯
1
,
y
¯
1
)
=
(
(
s
-
1
)
1
/
p
,
(
r
-
1
)
1
/
p
).
(iii)
Ifr
>
1 and s
=
1, then system (8) has the equilibrium point (
x
¯
2
,
y
¯
2
)
=
(
0
,
(
r
-
1
)
1
/
p
).
(iv)
Ifr
=
1 and s
>
1, then system (8) has the equilibrium point (
x
¯
3
,
y
¯
3
)
=
(
(
s
-
1
)
1
/
p
,
0
).
(v)
Ifr
∈
(
0,1
), s
=
1, and 1
/
p is an even positive integer, then system (8) has the equilibrium point (
x
¯
4
,
y
¯
4
)
=
(
0
,
(
r
-
1
)
1
/
p
).
(vi)
Ifr
=
1, s
∈
(
0,1
), and 1
/
p is an even positive integer, then system (8) has the equilibrium point (
x
¯
5
,
y
¯
5
)
=
(
(
s
-
1
)
1
/
p
,
0
).
(vii)
Ifr
,
s
∈
(
0,1
) and 1
/
p is an even positive integer, then it also has the positive equilibrium point (
x
¯
6
,
y
¯
6
)
=
(
(
s
-
1
)
1
/
p
,
(
r
-
1
)
1
/
p
).Proof.
The proof is easily obtained from the definition of equilibrium point.Before we give the following theorems about the local asymptotic stability of the aforementioned equilibrium points, we build the corresponding linearized form of system (8) and consider the following transformation: (14)
x
n
,
x
n
-
1
,
x
n
-
2
,
y
n
,
y
n
-
1
,
y
n
-
2
⟶
f
,
f
1
,
f
2
,
g
,
g
1
,
g
2
,where f
=
r
x
n
-
1
/
1
+
y
n
-
2
p, f
1
=
x
n, f
2
=
x
n
-
1, g
=
s
y
n
-
1
/
1
+
x
n
-
2
p, g
1
=
y
n, and g
2
=
y
n
-
1. The Jacobian matrix about the fixed point (
x
¯
,
y
¯
) under the above transformation is as follows: (15)
B
x
¯
,
y
¯
=
0
r
1
+
y
¯
p
0
0
0
-
r
p
x
¯
y
¯
p
-
1
1
+
y
¯
p
2
1
0
0
0
0
0
0
1
0
0
0
0
0
0
-
s
p
y
¯
x
¯
p
-
1
1
+
x
¯
p
2
0
s
1
+
x
¯
p
0
0
0
0
1
0
0
0
0
0
0
1
0
,where r
,
s
,
p
∈
(
0
,
∞
)
.Theorem 7.
(i) Ifr
<
1 and s
<
1, then the zero equilibrium point (
x
¯
0
,
y
¯
0
) is locally asymptotically stable.
(ii) Ifr
>
1 or s
>
1, then the zero equilibrium point (
x
¯
0
,
y
¯
0
) is locally unstable.Proof.
(i) The linearized system of (8) about the equilibrium point (
x
¯
0
,
y
¯
0
) is given by (16)
X
n
+
1
=
B
x
¯
0
,
y
¯
0
X
n
, where X
n
=
(
x
n
,
x
n
-
1
,
x
n
-
2
,
y
n
,
y
n
-
1
,
y
n
-
2
)
T and (17)
B
x
¯
0
,
y
¯
0
=
0
r
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
s
0
0
0
0
1
0
0
0
0
0
0
1
0
. The characteristic equation of B
(
x
¯
0
,
y
¯
0
) is as follows: (18)
P
λ
=
λ
6
-
r
+
s
λ
4
+
r
s
λ
2
=
0
. The roots of P
(
λ
) are (19)
λ
1
=
s
,
λ
2
=
-
s
,
λ
3,4
=
0
,
λ
5
=
r
,
λ
6
=
-
r
. Since all eigenvalues of the Jacobian matrix B about (
x
¯
0
,
y
¯
0
) lie inside the open unit disk λ
<
1, the zero equilibrium point is locally asymptotically stable.
(ii) It is easy to see that ifr
>
1 or s
>
1, then there exists at least one root λ of P
(
λ
) such that λ
>
1. Hence by Theorem 5, if r
>
1 or s
>
1, then (
x
¯
0
,
y
¯
0
) is unstable. Thus, the proof is complete.Theorem 8.
(i) Ifr
>
1 and s
>
1, then the positive equilibrium point (
x
¯
1
,
y
¯
1
) is locally unstable.
(ii) Ifr
<
1, s
<
1, and 1
/
p is an even positive integer, then the positive equilibrium point (
x
¯
6
,
y
¯
6
) is locally unstable.Proof.
(i) The linearized system of (8) about the equilibrium point (
x
¯
1
,
y
¯
1
)
=
(
(
s
-
1
)
1
/
p
,
(
r
-
1
)
1
/
p
) is given by (20)
X
n
+
1
=
B
x
¯
1
,
y
¯
1
X
n
, where X
n
=
(
x
n
,
x
n
-
1
,
x
n
-
2
,
y
n
,
y
n
-
1
,
y
n
-
2
)
T and (21)
B
x
¯
1
,
y
¯
1
=
0
1
0
0
0
-
p
s
-
1
1
/
p
r
-
1
p
-
1
/
p
r
1
0
0
0
0
0
0
1
0
0
0
0
0
0
-
p
r
-
1
1
/
p
s
-
1
p
-
1
/
p
s
0
1
0
0
0
0
1
0
0
0
0
0
0
1
0
. The characteristic equation of B
(
x
¯
1
,
y
¯
1
) is as follows: (22)
P
λ
=
λ
6
-
2
λ
4
+
λ
2
-
p
2
r
s
r
-
1
s
-
1
=
0
. It is clear that P
(
λ
) has a root in the interval (
1
,
∞
) since (23)
P
1
=
-
p
2
r
s
r
-
1
s
-
1
<
0
,
lim
λ
→
∞
⁡
P
λ
=
∞
. This completes the proof.
(ii) The proof is similar to the proof of (i), so it will be omitted.Theorem 9.
(i) Ifr
>
1 and s
=
1, then the equilibrium point (
x
¯
2
,
y
¯
2
) is nonhyperbolic point.
(ii) Ifr
=
1 and s
>
1, then the equilibrium point (
x
¯
3
,
y
¯
3
) is nonhyperbolic point.
(iii) Ifr
<
1, s
=
1, and 1
/
p is an even positive integer, then the equilibrium point (
x
¯
4
,
y
¯
4
) is nonhyperbolic point.
(iv) Ifr
=
1, s
<
1, and 1
/
p is an even positive integer, then the equilibrium point (
x
¯
5
,
y
¯
5
) is nonhyperbolic point.Proof.
(i) The linearized system of (8) about the equilibrium point (
x
¯
2
,
y
¯
2
)
=
(
0
,
(
r
-
1
)
1
/
p
) is given by (24)
X
n
+
1
=
B
x
¯
2
,
y
¯
2
X
n
, where X
n
=
(
x
n
,
x
n
-
1
,
x
n
-
2
,
y
n
,
y
n
-
1
,
y
n
-
2
)
T and (25)
B
x
¯
2
,
y
¯
2
=
0
1
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
s
0
0
0
0
1
0
0
0
0
0
0
1
0
. The characteristic equation of B
(
x
¯
2
,
y
¯
2
) is as follows: (26)
P
λ
=
λ
6
-
s
+
1
λ
4
+
s
λ
2
=
0
. The roots of P
(
λ
) are (27)
λ
1
=
s
,
λ
2
=
-
s
,
λ
3,4
=
0
,
λ
5
=
1
,
λ
6
=
-
1
. Thus, the equilibrium point (
x
¯
2
,
y
¯
2
)
=
(
0
,
(
r
-
1
)
1
/
p
) is nonhyperbolic point.
(ii) The proof is similar to the proof of (i), so it will be omitted.
(iii) The proof is easily seen from the proof of (i).
(iv) The proof is easily seen from the proof of (i).Now, we will study the global behavior of zero equilibrium point.Theorem 10.
Ifr
<
1 and s
<
1, then the zero equilibrium point is globally asymptotically stable.Proof.
We know by Theorem7 that the equilibrium point (
x
¯
0
,
y
¯
0
) of system (8) is locally asymptotically stable. Hence, it suffices to show for any x
n
,
y
n solution of system (8) that (28)
lim
n
→
∞
⁡
x
n
,
y
n
=
x
¯
0
,
y
¯
0
. Since (29)
0
≤
x
n
+
1
=
r
x
n
-
1
1
+
y
n
-
2
p
<
r
x
n
-
1
,
0
≤
y
n
+
1
=
s
y
n
-
1
1
+
x
n
-
2
p
<
s
y
n
-
1
, we have by induction (30)
x
2
n
-
1
<
r
n
x
-
1
,
x
2
n
<
r
n
x
0
,
y
2
n
-
1
<
s
n
y
-
1
,
y
2
n
<
s
n
y
0
. Thus, for r
<
1 and s
<
1, we obtain (31)
lim
n
→
∞
⁡
x
n
,
y
n
=
0,0
.This completes the proof.Theorem 11.
Consider system (8) and suppose that (32)
r
>
1
,
s
>
1
.Then, we obtain the following invariant intervals, for i
∈
{
0,1
,
2
}: (i)
(
x
-
i
,
y
-
i
)
∈
(
0
,
(
s
-
1
)
1
/
p
)
×
(
(
r
-
1
)
1
/
p
,
∞
)
⇒
(
x
n
,
y
n
)
∈
(
0
,
(
s
-
1
)
1
/
p
)
×
(
(
r
-
1
)
1
/
p
,
∞
) for n
≥
1
.
(ii)
(
x
-
i
,
y
-
i
)
∈
(
(
s
-
1
)
1
/
p
,
∞
)
×
(
0
,
(
r
-
1
)
1
/
p
)
⇒
(
x
n
,
y
n
)
∈
(
(
s
-
1
)
1
/
p
,
∞
)
×
(
0
,
(
r
-
1
)
1
/
p
) for n
≥
1
.Proof.
(i) Let(
x
-
i
,
y
-
i
)
∈
(
0
,
(
s
-
1
)
1
/
p
)
×
(
(
r
-
1
)
1
/
p
,
∞
) for i
∈
{
0,1
,
2
}
. From system (8), we have (33)
x
1
=
r
x
-
1
1
+
y
-
2
p
<
r
x
¯
1
1
+
y
¯
1
p
=
x
¯
1
=
s
-
1
1
/
p
,
y
1
=
s
y
-
1
1
+
x
-
2
p
>
s
y
¯
1
1
+
x
¯
1
p
=
y
¯
1
=
r
-
1
1
/
p
. We prove by induction that(34)
x
n
,
y
n
∈
0
,
s
-
1
1
/
p
×
r
-
1
1
/
p
,
∞
,
∀
n
>
1
.Suppose that (34) is true for n
=
k
>
1
. Then from system (8), we have (35)
x
k
+
1
=
r
x
k
-
1
1
+
y
k
-
2
p
<
r
x
¯
1
1
+
y
¯
1
p
=
x
¯
1
=
s
-
1
1
/
p
,
y
k
+
1
=
s
y
k
-
1
1
+
x
k
-
2
p
>
s
y
¯
1
1
+
x
¯
1
p
=
y
¯
1
=
r
-
1
1
/
p
.Therefore, (34) is true for all n
≥
-
2. This completes the proof of (i). Similarly, we can obtain the proof of (ii) which will be omitted.Corollary 12.
Consider system (8) and suppose that (36)
r
<
1
,
s
<
1
,
1
p
∈
2
Z
+
. Then, we obtain the following invariant intervals, for i
∈
{
0,1
,
2
}: (i)
(
x
-
i
,
y
-
i
)
∈
(
0
,
(
s
-
1
)
1
/
p
)
×
(
(
r
-
1
)
1
/
p
,
∞
)
⇒
(
x
n
,
y
n
)
∈
(
0
,
(
s
-
1
)
1
/
p
)
×
(
(
r
-
1
)
1
/
p
,
∞
) for n
≥
1.
(ii)
(
x
-
i
,
y
-
i
)
∈
(
(
s
-
1
)
1
/
p
,
∞
)
×
(
0
,
(
r
-
1
)
1
/
p
)
⇒
(
x
n
,
y
n
)
∈
(
(
s
-
1
)
1
/
p
,
∞
)
×
(
0
,
(
r
-
1
)
1
/
p
) for n
≥
1.Theorem 13.
Assume thatr
,
s
∈
(
1
,
∞
); then there exist unbounded solutions of system (8).Proof.
From Theorem11, we can assume without loss of generality that the solution x
n
,
y
n of system (8) is such that (37)
x
n
<
x
¯
1
=
s
-
1
1
/
p
,
y
n
>
y
¯
1
=
r
-
1
1
/
p
,
for
n
≥
-
2
.Then(38)
x
n
+
1
=
r
x
n
-
1
1
+
y
n
-
2
p
<
r
x
n
-
1
1
+
r
-
1
=
x
n
-
1
,
y
n
+
1
=
s
y
n
-
1
1
+
x
n
-
2
p
>
s
y
n
-
1
1
+
s
-
1
=
y
n
-
1
,from which it follows that (39)
lim
n
→
∞
⁡
x
n
=
0
,
lim
n
→
∞
⁡
y
n
=
∞
.This completes the proof.Corollary 14.
Assume thatr
,
s
∈
(
0,1
) and 1
/
p
∈
2
Z
+; then there exist unbounded solutions of system (8).Theorem 15.
Ifr
=
s
=
1, then system (8) possesses the prime period two solution which is one of the following forms: (40)
…
,
0
,
b
,
0
,
d
,
0
,
b
,
0
,
d
,
…
…
,
0
,
b
,
c
,
0
,
0
,
b
,
c
,
0
,
…
…
,
a
,
0
,
c
,
0
,
a
,
0
,
c
,
0
,
…
…
,
a
,
0
,
0
,
d
,
a
,
0
,
0
,
d
,
…with a
,
b
,
c
,
d
>
0.Proof.
Let(41)
…
,
a
,
b
,
c
,
d
,
a
,
b
,
c
,
d
,
… be a period two solution of system (8). Then, we have(42)
a
=
r
a
1
+
d
p
,
b
=
r
b
1
+
c
p
,
(43)
c
=
s
c
1
+
b
p
,
d
=
s
d
1
+
a
psuch that a
≠
b and c
≠
d. Firstly, we consider the case a
≠
0, b
≠
0. Then, we obtain from (42) that c
=
d
=
(
r
-
1
)
1
/
p which is a contradiction. Similarly, the case c
≠
0, d
≠
0 is impossible with (43). Thus, one of them must be equal to zero. Now we assume that a
=
0, b
≠
0, c
=
0, d
≠
0. In this sense, we obtain from (42) and (43) that r
=
s
=
1. Therefore, (44)
…
,
0
,
b
,
0
,
d
,
0
,
b
,
0
,
d
,
… is a two-periodic solution of system (8) with b
,
d
>
0. The other cases are similar and the proofs of them are omitted.


## 4. Conclusion and Future Work

System (7) is the extension of the third-order equation in [18]. Here, we studied some dynamics of a six-dimensional discrete system. We investigated the local asymptotic behavior of solutions of system (7) using linearization method. We also obtained the global asymptotic stability result for the zero equilibrium point by using analytical techniques. At the end of the paper, the existence of prime two-periodic solutions and the existence of unbounded solutions were proved.Fractional difference equations and their systems which contain powers of arbitrary positive degrees will next our aim to study.



---

*Source: 6842521-2016-06-28.xml*