Discrete Dynamics in Nature and Society
Volume 2016 (2016), Article ID 6842521, 7 pages
http://dx.doi.org/10.1155/2016/6842521

Research Article

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

Mehmet Gümüş and Yüksel Soykan

Department of Mathematics, Faculty of Science and Arts, Bülent Ecevit University, 67100 Zonguldak, Turkey

Received 25 November 2015; Accepted 2 June 2016

Academic Editor: Zhan Zhou

Copyright © 2016 Mehmet Gümüş and Yüksel Soykan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 [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:

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 , then

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:

, 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

where

and 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 [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) 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 that

Suppose 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

Then

from 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 have

such 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.

References

Q. Din and E. M. Elsayed, “Stability analysis of a discrete ecological model,” Computational Ecology and Software, vol. 4, no. 2, pp. 89–103, 2014.

A. Q. Khan, Q. Din, M. N. Qureshi, and T. F. Ibrahim, “Global behavior of an anti-competitive system of fourth-order rational difference equations,” Computational Ecology and Software, vol. 4, no. 1, pp. 35–46, 2014.

E. M. Elsayed, “Solutions of rational difference systems of order two,” Mathematical and Computer Modelling, vol. 55, no. 3-4, pp. 378–384, 2012.

A. Q. Khan, M. N. Qureshi, and Q. Din, “Global dynamics of some systems of higher-order rational difference equations,” Advances in Difference Equations, vol. 2013, article 354, 2013.

A. Y. Özban, “On the system of rational difference equations ,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 833–837, 2007.

I. Yalcinkaya, “On the global asymptotic stability of a second-order system of difference equations,” Discrete Dynamics in Nature and Society, vol. 2008, Article ID 860152, 12 pages, 2008.

I. Yalcinkaya, C. Cinar, and D. Simsek, “Global asymptotic stability of a system of difference equations,” Applicable Analysis, vol. 87, no. 6, pp. 677–687, 2008.

R. Abo-Zeid, “Global attractivity of a higher-order difference equation,” Discrete Dynamics in Nature and Society, vol. 2012, Article ID 930410, 11 pages, 2012.

D. Chen, X. Li, and Y. Wang, “Dynamics for nonlinear difference equation ,” Advances in Difference Equations, vol. 2009, Article ID 235691, 2009.

H. M. El-Owaidy, A. M. Youssef, and A. M. Ahmed, “On the Dynamics of ,” Rostocker Mathematisches Kolloquium, vol. 59, pp. 11–18, 2005.

M. Gümüş, “The periodicity of positive solutions of the non-linear difference equation ,” Discrete Dynamics in Nature and Society, vol. 2013, Article ID 742912, 3 pages, 2013.

Ö. Öcalan, H. Ogünmez, and M. Gümüş, “Global behavior test for a nonlinear difference equation with a period-two coefficient,” Dynamics of Continuous, Discrete and Impulsive Systems Series A: Mathematical Analysis, vol. 21, no. 3-4, pp. 307–316, 2014.

A. S. Kurbanlı, C. Çinar, and I. Yalçinkaya, “On the behavior of positive solutions of the system of rational difference equations , ,” Mathematical and Computer Modelling, vol. 53, no. 5, pp. 1261–1267, 2011.

G. Papaschinopoulos and C. J. Schinas, “On a system of two nonlinear difference equations,” Journal of Mathematical Analysis and Applications, vol. 219, no. 2, pp. 415–426, 1998.

Q. Zhang, L. Yang, and J. Liu, “Dynamics of a system of rational third-order difference equation,” Advances in Difference Equations, vol. 2012, article 136, 2012.

Q. Din, M. N. Qureshi, and A. Q. Khan, “Dynamics of a fourth-order system of rational difference equations,” Advances in Difference Equations, vol. 2012, no. 1, article 215, 15 pages, 2012.

N. Touafek and E. M. Elsayed, “On the solutions of systems of rational difference equations,” Mathematical and Computer Modelling, vol. 55, no. 7-8, pp. 1987–1997, 2012.

H. M. El-Owaidy, A. M. Ahmed, and A. M. Youssef, “The dynamics of the recursive sequence ,” Applied Mathematics Letters, vol. 18, no. 9, pp. 1013–1018, 2005.

E. Camouzis and G. Ladas, Dynamics of Third-Order Rational Difference Equations with Open Problems and Conjectures, vol. 5, CRC Press, 2007.

V. L. Kocic and G. Ladas, Global Behavior of Nonlinear Difference Equations of Higher Order with Applications, Kluwer Academic, Dordrecht, The Netherlands, 1993.

M. R. S. Kulenović and G. Ladas, Dynamics of Second Order Rational Difference Equations, Chapman & Hall/CRC Press, 2001.

H. Sedaghat, Nonlinear Difference Equations Theory with Applications to Social Science Models, vol. 15, Springer, 2003.