ISRN International Scholarly Research Notices 2356-7872 Hindawi Publishing Corporation 10.1155/2014/760502 760502 Research Article On the System of High Order Rational Difference Equations Zhang Qianhong 1 Zhang Wenzhuan 1 Shao Yuanfu 2 Liu Jingzhong 3 Zhang Chunrui 1 Guizhou Key Laboratory of Economics System Simulation Guizhou University of Finance and Economics, Guiyang, Guizhou 550004 China gufe.edu.cn 2 School of Science, Guilin University of Technology Guilin, Guangxi 541000 China gutcie.cn 3 Department of Mathematics and Physics, Hunan Institute of Technology Hengyang, Hunan 421002 China hnit.edu.cn 2014 29 10 2014 2014 15 03 2014 02 07 2014 22 07 2014 29 10 2014 2014 Copyright © 2014 Qianhong Zhang et al. 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.

This paper is concerned with the boundedness, persistence, and global asymptotic behavior of positive solution for a system of two rational difference equations xn+1=A+(xn/i=1kyn-i),  yn+1=B+(yn/i=1kxn-i),n=0,1,,k{1,2,}, where A,B0,,x-i0,, and y-i0,,i=0,1,2,,k.

1. Introduction

In this paper, we study the global behavior of solutions of the following system: (1)xn+1=A+xni=1kyn-i,yn+1=B+yni=1kxn-i,n=0,1,, where A,B are positive constants and initial conditions x-i,y-i0,,  i=0,1,2,,k.

A pair of sequences of positive real numbers {(xn,yn)} that satisfies (1) is a positive solution of (1). If a positive solution of (1) is a pair of positive constants (x,y), then the solution is the equilibrium solution.

A positive solution {(xn,yn)} of (1) is bounded and persists, if there exist positive constants M,N such that (2)Mxn,ynN,n=-2,-1,.

In 1998, DeVault et al.  proved that every positive solution of the difference equation (3)xn+1=A+xnxn-1,n=0,1,, where A(0,), oscillates about the positive equilibrium c=1+A of (3). Moreover, every positive solution of (3) is bounded away from zero and infinity. Also the positive equilibrium of (3) is globally asymptotically stable.

In 2003, Abu-Saris and DeVault  studied the following recursive difference equation: (4)xn+1=A+xnxn-k,n=0,1,, where A1,+,  x-k,x-k+1,,x0 are positive real numbers.

Papaschinopoulos and Schinas  investigated the global behavior for a system of the following two nonlinear difference equations: (5)xn+1=A+ynxn-p,yn+1=A+xnyn-q,n=0,1,, where A is a positive real number, p,q are positive integers, and x-p,,x0,  y-q,,y0 are positive real numbers.

In 2012, Zhang et al.  investigated the global behavior for a system of the following third-order nonlinear difference equations: (6)xn+1=A+xnyn-1+yn-2,yn+1=B+ynxn-1+xn-2, where A,B(0,), and the initial values x-i,y-i(0,),  i=0,1. For other related results, the reader can refer to .

Motivated by the discussion above, we study the global asymptotic behavior of solutions for system (1). More precisely, we prove the following: if A>1/k,B>1/k then every positive solution {(xn,yn)} of (1) is persistent and bounded. Moreover, we prove that every positive solution {(xn,yn)} of (1) converges the unique positive equilibrium (x,y) as n.

2. Main Results

In the following lemma, we show boundedness and persistence of the positive solutions of (1).

Lemma 1.

Consider (1). Suppose that (7)A>1k,B>1k are satisfied. Then, every positive solution (xn,yn) of (1) is satisfied, for n=k+1,k+2,(8)Axn1kBn-kxk-kABkB-1+kABkB-1Byn1kAn-kyk-kABkA-1+kABkA-1.

Proof.

Let {(xn,yn)} be a positive solution of (1). Since xn>0 and yn>0 for all n1, (1) implies that (9)xnA,ynB,n=1,2,3,

Moreover, using (1) and (9), we have (10)xnA+1kBxn-1,ynB+1kAyn-1,n=k+1,k+2,. Let vn,wn be the solution of the system, respectively, (11)vn=A+1kBvn-1,wn=B+1kAwn-1,nk+1, such that (12)vi=xi,wi=yi,i=1,2,,k. We prove by induction that (13)xnvn,ynwn,nk+1. Suppose that (13) is true for n=mk+1. Then, from (10), we get (14)xm+1A+1kBxmA+1kBvm=vm+1,  ym+1B+1kAymB+1kBwm=wm+1. Therefore, (13) is true. From (11) and (12), we have (15)vn=1kBn-kxk-kABkB-1+kABkB-1,wn=1kAn-kyk-kABkA-1+kABkA-1,nk. Then, from (9), (13), and (15), the proof of the relation (8) follows immediately.

Theorem 2.

Consider the system of difference equation (1). If relation (7) is satisfied, then the following statements are true.

Equation (1) has a unique positive equilibrium (x,y) given by (16)x=k2AB-1kkB-1,y=k2AB-1kkA-1.

Every positive solution (xn,yn) of system (1) tends to the positive equilibrium (x,y) of (1) as n.

Proof.

(i) Let x and y be positive numbers such that (17)x=A+xky,y=B+ykx. Then, from (7) and (17), we have that the positive solution (x,y) is given by (16). This completes the proof of Part (i).

(ii) From (1) and (8), we have (18)limnsupxn=L1,limninfxn=l1,limnsupyn=L2,limninfyn=l2, where li,Li(0,),  i=1,2. Then, from (1) and (18), we get (19)L1A+L1kl2,l1A+l1kL2,  L2B+L2kl1,l2B+l2kL1, from which we have (20)L1kB-1l2kA-1,L2kA-1l1kB-1. Then, relations (7) and (20) imply that L1L2l1l2, from which it follows that (21)L1L2=l1l2. We claim that (22)L1=l1,L2=l2. Suppose on the contrary that l1<L1. Then, from (21), we have L1L2=l1l2<L1l2 and so L2<l2 which is a contradiction. So L1=l1. Similarly, we can prove that L2=l2. Therefore, (22) is true. Hence, from (1) and (22), there exist the limxn and limyn, as n and (23)limnxn=x,limnyn=y, where (x,y) is the unique positive equilibrium of (1). This completes the proof of Part (ii). The proof of Theorem 2 is completed.

Theorem 3.

Consider the system of difference equation (1). If relation (7) is satisfied and assuming that (24)k2AB-1kA-1+k2AB-1kB-1<1, then the unique positive equilibrium (x,y) is locally asymptotically stable.

Proof.

From Theorem 2, the system of difference equation (1) has a unique equilibrium (x,y). The linearized equation of system (1) about the equilibrium point (x,y) is (25)Ψn+1=BΨn, where Ψn=(xn,,xn-k,yn,,yn-k)T, and(26)B(2k+2)×(2k+2)=1ky0000-xk2y2-xk2y2-xk2y2100000000001000000-yk2x2-yk2x2-yk2x21kx0000000100000000010. Let λ1,λ2,,λ2k+2 denote the eigenvalues of matrix B and let D=diag(d1,d2,,d2k+2) be a diagonal matrix, where d1=dk+2=1,  di=dk+1+i=1-iɛ(i=2,,k+1), and (27)0<ɛ<min1k+11-x+yky2,1k+11-x+ykx2. Clearly, D is invertible. Computing matrix DBD-1, we obtain that(28)DBD-1=1ky0000-xk2y2d1dk+3-1-xk2y2d1d2k+1-1-xk2y2d1d2k+2-1d2d1-1000000000dk+1dk-1000000-yk2x2dk+2d2-1-yk2x2dk+2dk-1-yk2x2dk+2dk+1-11kx0000000dk+3dk+2-1000000000d2k+2d2k+1-10.From d1>d2>>dk+1>0 and dk+2>dk+3>>d2k+2>0, imply that (29)d2d1-1<1,d3d2-1<1,,dk+1dk-1<1,  dk+3dk+2-1<1,,d2k+2d2k+1-1<1. Furthermore, noting (7), (24), and (27), we have (30)1ky+xk2y2d1dk+3-1++xk2y2d1d2k+2-1=1ky+xk2y211-2ɛ++11-(k+1)ɛ<1ky+xky211-k+1ɛ<1,1kx+yk2x2dk+2d2-1++yk2x2dk+2dk+1-1=1kx+yk2x211-2ɛ++11-(k+1)ɛ<1kx+ykx211-k+1ɛ<1. It is well known that B has the same eigenvalues as DBD-1; we have that (31)max1i2k+2λiDBD-1=max+yk2x211-2ɛ++11-(k+1)ɛd2d1-1,,dk+1dk-1,dk+3dk+2-1,,d2k+2d2k+1-1,1ky+xk2y211-2ɛ++11-k+1ɛ,1kx+yk2x211-2ɛ++11-(k+1)ɛ<1. This implies that the equilibrium (x,y) of (1) is locally asymptotically stable.

Combining Theorem 2 with Theorem 3, we obtain the following theorem.

Theorem 4.

Consider the system of difference equation (1). If relations (7) and (24) are satisfied, then the unique positive equilibrium (x,y) is globally asymptotically stable.

3. Some Numerical Examples

In order to illustrate the results of the previous sections 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 nonlinear difference equations and system of nonlinear difference equations.

Example 1.

Consider the following difference equations: (32)xn+1=0.8+xnyn-1+yn-2+yn-3,  yn+1=0.6+ynxn-1+xn-2+xn-3, with the initial values x-i=y-i=0.5(i=1,2,3). Then, the solution (xn,yn) of system (32) is bounded and persists and the system has a unique equilibrium (x,y)=(1.3833,0.7905) which is globally asymptotically stable (see Figure 1).

The dynamics of system (32).

Example 2.

Consider the following difference equations: (33)xn+1=0.8+xnyn-1+yn-2+yn-3+yn-4,  yn+1=0.6+ynxn-1+xn-2+xn-3+xn-4, with the initial values x-i=y-i=1.5(i=1,2,3,4). Then, the solution (xn,yn) of system (33) is bounded and persists and the system has a unique equilibrium (x,y)=(1.1929,0.7591) which is globally asymptotically stable (see Figure 2).

The dynamics of system (33).

4. Conclusion

In this paper, we study the dynamics of a system of high order difference equation (34)xn+1=A+xni=1kyn-i,yn+1=B+yni=1kxn-i,    n=0,1,,k1,2,. It concluded that, under condition A>1/k,B>1/k, the positive solution (xn,yn) of this system is bounded and persists; moreover, if k2AB-1/kA-1+k2AB-1/kB-1<1, it converges asymptotically the unique equilibrium (x,y).

We conclude the paper by presenting the following open problem.

Open Problem. Consider the system of difference equation (1) with A1/k and B1/k. Find the set of all initial conditions that generate bounded solutions. In addition, investigate global behavior of these solutions.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (Grant no. 11361012), the Scientific Research Foundation of Guizhou Provincial Science and Technology Department (J2083, J2061), and the Natural Science Foundation of Guizhou Provincial Educational Department (no. 2008040).

DeVault R. Ladas G. Schultz S. W. On the recursive sequence xn+1=A/xn+1/xn-2 Proceedings of the American Mathematical Society 1998 126 11 3257 3261 10.1090/S0002-9939-98-04626-7 Abu-Saris R. M. DeVault R. Global stability of yn+1=A+yn/yn-k Applied Mathematics Letters 2003 16 173 178 10.1016/S0893-9659(03)80028-9 Papaschinopoulos G. Schinas C. J. On a system of two nonlinear difference equations Journal of Mathematical Analysis and Applications 1998 219 2 415 426 10.1006/jmaa.1997.5829 MR1606350 2-s2.0-0000475245 Zhang Q. Yang L. Liao D. Global asymptotic behavior of positive solutions to the system of rational difference equations Journal of Southwest University 2012 7 12 15 Zhang Q. Yang L. Liu J. Dynamics of a system of rational third-order difference equation Advances in Difference Equations 2012 2012, article 136 10.1186/1687-1847-2012-136 MR2992065 2-s2.0-84871311577 Amleh A. M. Grove E. A. Ladas G. Georgiou D. On the recursive sequence xn+1=A+(xn-1/xn) Journal of Mathematical Analysis and Applications 1999 233 2 790 798 10.1006/jmaa.1999.6346 MR1689579 He W. Li W. Yan X. Global attractivity of the difference equation xn+1=a+xn-k/xn Applied Mathematics and Computation 2004 151 3 879 885 10.1016/S0096-3003(03)00528-9 MR2052466 DeVault R. Ladas G. Schultz S. W. Necessary and sufficient conditions the boundedness of xn+1=A/xnp+B/xn-1q Journal of Difference Equations and Applications 1998 3 259 266 Elaydi S. N. An Introduction to Difference Equations 1996 New York, NY, USA Springer 10.1007/978-1-4757-9168-6 MR1410259 Agarwal R. P. Li W. Pang P. Y. H. Asymptotic behavior of a class of nonlinear delay difference equations Journal of Difference Equations and Applications 2002 8 8 719 728 10.1080/1023619021000000735 MR1914599 ZBL1010.39003 2-s2.0-0036699887 Kocic V. L. Ladas G. Global Behavior of Nonlinear Difference Equations of Higher Order with Applications 1993 Dordrecht, Netherlands Kluwer Academic 10.1007/978-94-017-1703-8 MR1247956 Kulenonvic M. R. S. Ladas G. Dynamics of Second Order Rational Difference Equations with Open Problems and Conjectures 2002 Boca Raton, Fla, USA Chapaman & Hall/CRC Li W. Sun H. Dynamics of a rational difference equation Applied Mathematics and Computation 2005 163 2 577 591 10.1016/j.amc.2004.04.002 MR2121812 2-s2.0-13544252478 Su Y. Li W. Global attractivity of a higher order nonlinear difference equation Journal of Difference Equations and Applications 2005 11 10 947 958 10.1080/10236190500273333 MR2173724 2-s2.0-27344435108 Hu L. Li W. Global stability of a rational difference equation Applied Mathematics and Computation 2007 190 2 1322 1327 10.1016/j.amc.2007.02.101 MR2339725 2-s2.0-34250647876 Ibrahim T. F. Two-dimensional fractional system of nonlinear difference equations in the modelling competitive populations International Journal of Basic & Applied Sciences 2012 12 5 103 121 Ibrahim T. F. Closed form solution of a symmetric competitive system of rational difference equations Studies in Mathematical Sciences 2012 5 1 49 57 Ibrahim T. F. Zhang Q. Stability of an anti-competitive system of rational difference equations Archives des Sciences 2013 66 44 58