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: , , , where the parameters and the initial values for are positive real numbers.

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 [37] 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 [812]).

In [13], Kurbanlı et al. studied the behavior of positive solutions of the following system of difference equations:where the initial conditions are arbitrary nonnegative real numbers.

In [14], Papaschinopoulos and Schinas studied the system of two nonlinear difference equations:where and are positive integers.

In [15], Zhang et al. studied the dynamical behavior of positive solutions for a system for third-order rational difference equations:

In [16], Din et al. studied the dynamics of a system of fourth-order rational difference equations:

In [17], Touafek and Elsayed investigated the behavior of solutions of systems of difference equations:with a nonzero real number’s initial conditions.

In [18], El-Owaidy et al. investigated the global behavior of the following difference equation: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:where the parameters and the initial values for 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:by the change of variables and with and . So, in order to study system (7), we investigate system (8).

By using the induction and the equations in (8) we see that if are positive real numbers for , thenwhich 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:, where and are continuously differentiable functions and are some intervals of real numbers. Also, a solution of system (10) is uniquely determined by initial values for .

Definition 1. An equilibrium point of system (10) is a point that satisfies Together with system (10), if we consider the associated vector map , then the point is also called a fixed point of the vector map

Definition 2. Let be an equilibrium point of system (10).(i)An equilibrium point is said to be stable if, for every , there exists such that for every initial value , with , , implying , for .(ii)If an equilibrium point is not stable, then it is said to be unstable.(iii)If an equilibrium point is stable and there exists such that , and as , then it is said to be asymptotically stable.(iv)If as , then an equilibrium point is said to be a global attractor.(v)If an equilibrium point is both global attractor and stable, then it is said to be globally asymptotically stable.

Definition 3. Let be an equilibrium point of the map where and are continuously differentiable functions at The linearized system of (10) about the equilibrium point is whereand is a Jacobian matrix of system (10) about the equilibrium point

Definition 4. For the system , of difference equations such that is a fixed point of , if no eigenvalues of the Jacobian matrix about have absolute value equal to one, then is called hyperbolic. If there exists an eigenvalue of the Jacobian matrix about with absolute value equal to one, then 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 of system (10).

Theorem 5. For the system , of difference equations such that is a fixed point of , if all eigenvalues of the Jacobian matrix about lie inside the open unit disk , then is locally asymptotically stable. If one of them has a modulus greater than one, then is unstable.

For other basic knowledge about difference equations and their systems, the reader is referred to books [1922].

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) is always the equilibrium point of system (8).(ii)If and , then system (8) has the equilibrium point .(iii)If and , then system (8) has the equilibrium point .(iv)If and , then system (8) has the equilibrium point .(v)If , , and is an even positive integer, then system (8) has the equilibrium point .(vi)If , , and is an even positive integer, then system (8) has the equilibrium point .(vii)If and is an even positive integer, then it also has the positive equilibrium point .

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: where , , , , , and . The Jacobian matrix about the fixed point under the above transformation is as follows: where

Theorem 7. (i) If and , then the zero equilibrium point is locally asymptotically stable.
(ii) If or , then the zero equilibrium point is locally unstable.

Proof. (i) The linearized system of (8) about the equilibrium point is given by where and The characteristic equation of is as follows: The roots of are Since all eigenvalues of the Jacobian matrix about lie inside the open unit disk , the zero equilibrium point is locally asymptotically stable.
(ii) It is easy to see that if or , then there exists at least one root of such that . Hence by Theorem 5, if or , then is unstable. Thus, the proof is complete.

Theorem 8. (i) If and , then the positive equilibrium point is locally unstable.
(ii) If , , and is an even positive integer, then the positive equilibrium point is locally unstable.

Proof. (i) The linearized system of (8) about the equilibrium point is given by where and The characteristic equation of is as follows: It is clear that has a root in the interval since This completes the proof.
(ii) The proof is similar to the proof of (i), so it will be omitted.

Theorem 9. (i) If and , then the equilibrium point is nonhyperbolic point.
(ii) If and , then the equilibrium point is nonhyperbolic point.
(iii) If , , and is an even positive integer, then the equilibrium point is nonhyperbolic point.
(iv) If , , and is an even positive integer, then the equilibrium point is nonhyperbolic point.

Proof. (i) The linearized system of (8) about the equilibrium point is given by where and The characteristic equation of is as follows: The roots of are Thus, the equilibrium point 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. If and , then the zero equilibrium point is globally asymptotically stable.

Proof. We know by Theorem 7 that the equilibrium point of system (8) is locally asymptotically stable. Hence, it suffices to show for any solution of system (8) that Since we have by induction Thus, for and , we obtain This completes the proof.

Theorem 11. Consider system (8) and suppose that Then, we obtain the following invariant intervals, for : (i) for (ii) for

Proof. (i) Let for From system (8), we have We prove by induction thatSuppose that (34) is true for Then from system (8), we have Therefore, (34) is true for all . 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 Then, we obtain the following invariant intervals, for : (i) for .(ii) for .

Theorem 13. Assume that ; then there exist unbounded solutions of system (8).

Proof. From Theorem 11, we can assume without loss of generality that the solution of system (8) is such that Thenfrom which it follows that This completes the proof.

Corollary 14. Assume that and ; then there exist unbounded solutions of system (8).

Theorem 15. If , then system (8) possesses the prime period two solution which is one of the following forms: with .

Proof. Let be a period two solution of system (8). Then, we havesuch that and . Firstly, we consider the case , . Then, we obtain from (42) that which is a contradiction. Similarly, the case , is impossible with (43). Thus, one of them must be equal to zero. Now we assume that , , , . In this sense, we obtain from (42) and (43) that . Therefore, is a two-periodic solution of system (8) with . 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.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.