RESEARCH Open Access Dynamic behavior of a nonlinear rational difference equation and generalization Qihong Shi * , Qian Xiao, Guoqiang Yuan and Xiaojun Liu * Correspondence: shiqh03@163. com Department of Basic Courses, Hebei Finance University, Baoding 071000, PR China Abstract This paper is concerned about the dynamic behavior for the following high order nonlinear difference equation x n =(x n-k + x n-m + x n-l )/(x n-k x n-m + x n-m x n-l +1) with the initial data {x −l , x −l+1 , , x −1 }∈R l + and 1 ≤ k ≤ m ≤ l. The convergence of solution to this equation is investigated by introducing a new sequence, which extends and includes corresponding results obtained in the references (Li in J Math Anal Appl 312:103-111, 2005; Berenhaut et al. Appl. Math. Lett. 20:54-58, 2007; Papaschinopoulos and Schinas J Math Anal Appl 294:614-620, 2004) to a large extent. In addition, some propositions for generalized equations are reported. Keywords: Nonlinear; Difference equation, Global stability, Positive solution 1 Introduction Our aim in this paper is to study the dynamical behavior of the following equation x n = x n−k + x n−m + x n−l x n − k x n − m + x n − m x n − l +1 , n = 0,1,2, . (1:1) where the initial data {x −l , x −l+1 , , x −1 }∈R l + and 1 ≤ k ≤ m ≤ l The study of properties of similar difference equations has been an area of intense interest in recent years [1-3]. There have been a lot of work concerning the behavior of the solution. In particular, Çinar [4] studied the properties of positive solution to x n+1 = x n−1 1+x n x n −1 , n =0,1, . (1:2) Yang et al. [5] investigated the qualitative behavior of the recursive sequence x n+1 = ax n−1 + bx n−2 c + dx n −1 x n −2 , n =0,1, , (1:3) Li et al. [6] studied the global as ymptotic of the following nonlinear difference equa- tion x n+1 = x n−1 x n−2 x n−3 + x n−1 + x n−2 + x n−3 + a 1+x n −1 x n −2 + x n −1 x n − 3 + x n −2 x n − 3 + a , n =0,1, , (1:4) with a ≥ 0. For more similar work, one can refer to [7-9] and references therein. Investigation of the equation (1.1) is motivated by the above studies. However, due to the special non- linear relation, the methods mentio ned in the references [4,5,7] do not always work for Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 © 2011 Shi et al; licensee Springer. This is an Open Access article distributed under the terms of th e Creative Commons Attribution License (http://creativecom mons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. the equation (1.1). In fact, equation (1.1) has lost the perfect symmetry. To this end, we introduce a simple transformed sequence to construct a contraction to prove the convergence of solutions, and apply this way solving a class of general equation. The rest of this paper proceeds as follows. In Sect. 2, we introduce some definitions and preliminary lemmas. Section 3 contains the main results and their proofs. In Sect. 4, we prove the stability for generalized rational difference equations and present our conjectures for similar equations. 2 Preliminaries In this section, we introdu ce some basic but important preliminary lemmas and nota- tion. For any x i Î ℝ + , we define a new sequence as x ∗ i =max{x i ,1/x i } . With the help of the transformed sequence {x ∗ i } , we can deduce the following conclusion. Lemma 1. Suppose the function f is defined by f (x, y, z)= x + y + z x y + y z +1 , (2:1) then f is decreasing in x and z if and onl y if y > 1 and increasing in x and z if and only if y < 1. Similarly, f is decreasing in y if and only if x + z > 1, conversely, it is increasing in y. Proof. This conclusion follows directly from the fact ∂ ∂x f (x, y, z)= 1 − y 2 ( xy + yz +1 ) 2 , (2:2) and ∂ ∂y f (x, y, z)= 1 − (x + z) 2 ( xy + yz +1 ) 2 . (2:3) Since x and z is symmetrical, then the proposition is obvious. □ Moreover, we can also prove the following contraction lemma which is useful in showing convergence of solutions in the transformed space mentioned in first para- graph of this section. Lemma 2. Suppose x n satisfying the equation (1.1), for any n ≥ l and (x n−k , x n−m , x n−l ) ∈ R 3 + , we have 1 ≤ x ∗ n ≤ x ∗ n − m . (2:4) Proof. Noticed that x n−k + x n−m + x n−l − ( 1+x n−k x n−m + x n−m x n−l ) = − ( x n−m − 1 )( x n−k + x n−l − 1 ) (2:5) and hence from (1.1), x n ≤ 1wheneverx n -m -1andx n-k + x n- l -1areofthesame signs, otherwise, x n ≥ 1. Le t x n- k = u, x n-m = v, x n- l = w. The RHS of (2.4) is obvious. Next we prove the LHS part. Indeed we have eight cases to consider. when (1 - v)(u + w -1)≥ 0, then x ∗ n = x n = u + v + w uv + vw +1 . Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 Page 2 of 8 Case (1) (u ≤ 1, v ≤ 1, w ≥ 1, u + w ≥ 1).Here,bylemma1,notethatv* ≥ 1, we have x ∗ n = 1 u ∗ + 1 v ∗ + w ∗ 1+ 1 u ∗ 1 v ∗ + 1 v ∗ w ∗ = u ∗ + v ∗ + u ∗ v ∗ w ∗ 1+u ∗ v ∗ + u ∗ w ∗ ≤ v ∗ . (2:6) Case (2) (u ≥ 1, v ≤ 1, w ≤ 1, u + w ≥ 1). Here, since v* ≥ 1, w* ≥ 1, we have x ∗ n = u ∗ + 1 v ∗ + 1 w ∗ u ∗ 1 v ∗ + 1 v ∗ 1 w ∗ +1 = w ∗ + v ∗ + u ∗ v ∗ w ∗ 1+w ∗ u ∗ + v ∗ w ∗ ≤ v ∗ . Case (3) (u ≥ 1, v ≤ 1, w ≥ 1, u + w ≥ 1). Similarly, x ∗ n = u ∗ + 1 v ∗ + w ∗ 1+u ∗ 1 v ∗ + 1 v ∗ w ∗ = u ∗ v ∗ +1+v ∗ w ∗ u ∗ + v ∗ + w ∗ ≤ v ∗ . Case (4) (u ≤ 1, v ≥ 1, w ≤ 1, u + w ≤ 1). Here, x ∗ n = 1 u ∗ + v ∗ + 1 w ∗ 1+ 1 u ∗ v ∗ + v ∗ 1 w ∗ = u ∗ + w ∗ + u ∗ v ∗ w ∗ u ∗ v ∗ + u ∗ w ∗ + v ∗ w ∗ ≤ v ∗ . (2:7) Oppositely, if (1 - v)(u + w -1)≤ 0, from the definition of x*, it is obvious that x ∗ n = 1 x n = uv + vw +1 u + v + w . Case (5) (u ≤ 1, v ≤ 1, w ≤ 1, u + w ≤ 1). By definition of x ∗ n , Lemma 1 and the fact v* ≥ 1, we have x ∗ n = 1+ 1 u ∗ 1 v ∗ + 1 v ∗ 1 w ∗ 1 u ∗ + 1 v ∗ + 1 w ∗ = u ∗ + w ∗ + u ∗ v ∗ w ∗ v ∗ w ∗ + u ∗ v ∗ + u ∗ w ∗ ≤ v ∗ . (2:8) Case (6) (u ≥ 1, v ≥ 1, w ≤ 1, u + w ≥ 1). Here, we have x ∗ n = 1+u ∗ v ∗ + v ∗ 1 w ∗ u ∗ + v ∗ + 1 w ∗ = w ∗ + v ∗ + u ∗ v ∗ w ∗ 1+w ∗ v ∗ + u ∗ w ∗ ≤ v ∗ . (2:9) Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 Page 3 of 8 Case (7) (u ≤ 1, v ≥ 1, w ≥ 1, u + w ≥ 1). Similarly, we have x ∗ n = 1+ 1 u ∗ v ∗ + v ∗ w ∗ w ∗ + v ∗ + 1 u ∗ = u ∗ + v ∗ + u ∗ v ∗ w ∗ 1+u ∗ v ∗ + u ∗ w ∗ ≤ v ∗ . (2:10) Case (8) (u ≥ 1, v ≥ 1, w ≥ 1, u + w ≥ 1). By the same way, we have x ∗ n = 1+u ∗ v ∗ + v ∗ w ∗ u ∗ + v ∗ + w ∗ ≤ v ∗ . (2:11) Inequalities (2.6)-(2.13) suggest our claim. □ Remark. In fact, by Lemma 1 and in view of u* ≥ 1andw* ≥ 1, the result x ∗ n ≤ u ∗ , x ∗ n ≤ w ∗ can also be derived from the argument for front eight different cases. Now let X n =max n−l ≤ i ≤ n−1 {x ∗ i } for all n ≥ l.ByLemma2,wecandeducethefollowing consequence. Lemma 3. The sequence {X i } is monotonically non-increasing in i which is much greater than l. Since X i ≥ 1fori ≥ l, Lemma 3 implies that as i tends to i nfinity, the sequence {X i } convergence to some limit, denote X, where X ≥ 1. 3 Convergence of solutions In what follows, we state and prove our main result in the sequence space. Theorem 1. Suppose the initial data of equation (1.1) (x −l , x −l+1 , , x −1 ) ∈ R l + . Then the solution sequence {x i } converges to the unique positive equilibrium ¯ x = 1 . Proof. Note that it suffices to show that the transformed sequence {x ∗ i } converges to 1. By the definition of X i , the values of X i are taken o n by entries in the sequence {x ∗ i } , and as well, by Lemma 2, {x ∗ i }∈[1, X i ] for i ≥ m.SupposeX >1,thenforanyε Î (0, X), we can find an N such that {x ∗ N }∈[X , X + ε ] , and for i ≥ N - l, {x ∗ i }∈[1, X + ε ] . Next we consider the eight possible cases again, and show that X = 1. From the defi- nitions of x ∗ i , X i and X, the result follows. Case (1) (u ≤ 1, v ≤ 1, w ≥ 1, u + w ≥ 1). Here, by lemma 1, we have X ≤ x ∗ n = 1 u ∗ + 1 v ∗ + w ∗ 1+ 1 u ∗ 1 v ∗ + 1 v ∗ w ∗ ≤ 1+ 1 X+ε + X + ε 2+ 1 X+ ε . (3:1) Hence 2X 2 +2Xε + X ≤ (1 + X + ε)(X + ε)+1 , ⇒ X 2 ≤ 1+ε + ε 2 . (3:2) Case (2) (u ≥ 1, v ≤ 1, w ≤ 1, u + w ≥ 1). The argument is identical to that in Case (1). Case (3) (u ≥ 1, v ≤ 1, w ≥ 1, u + w ≥ 1). Here, X ≤ x ∗ n = u ∗ + 1 v ∗ + w ∗ 1+u ∗ 1 v ∗ + 1 v ∗ w ∗ ≤ 2(X + ε)+ 1 X+ε 3 . (3:3) Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 Page 4 of 8 Therefore 3X 2 +3Xε ≤ 2(X + ε) 2 +1, ⇒ (X − ε 2 ) 2 ≤ 1+ 9 4 ε 2 , ⇒ X ≤ 1+ 9 4 ε 2 + ε 2 . (3:4) Case (4) (u ≤ 1, v ≤ 1, w ≤ 1, u + v ≤ 1). Here, X ≤ x ∗ n = 1 u ∗ + v ∗ + 1 w ∗ 1+ 1 u ∗ v ∗ + v ∗ 1 w ∗ ≤ 2 X+ε + X + ε 3 . (3:5) From this, we have ⇒ (X + ε) 2 ≤ 1+ 9 1 6 ε 2 . (3:6) Namely, X ≤ 1+ 9 16 ε 2 − ε 4 . (3:7) Case (5) (u ≤ 1, v ≤ 1, w ≤ 1, u + w ≤ 1). We have X ≤ x ∗ n = 1+ 1 u ∗ 1 v ∗ + 1 v ∗ 1 w ∗ 1 u ∗ + 1 v ∗ + 1 w ∗ ≤ 1+( 1 X+ε ) 2 +( 1 X+ε ) 2 3 X+ ε (3:8) which also implies X ≤ 1+ 9 16 ε 2 − ε 4 . (3:9) Case (6) (u ≥ 1, v ≥ 1, w ≤ 1, u + w ≥ 1). Here, X ≤ x ∗ n = 1+u ∗ v ∗ + v ∗ 1 w ∗ u ∗ + v ∗ + 1 w ∗ ≤ 1+(X + ε) 2 +(X + ε) 1+2(X + ε) . (3:10) We have 2X(X + ε)+X ≤ 1+(X + ε) 2 +(X + ε) , ⇒ X 2 ≤ 1+ε + ε 2 (3:11) Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 Page 5 of 8 Case (7) (u ≤ 1, v ≥ 1, w ≥ 1, u + w ≥ 1). Here, it follows X ≤ x ∗ n = 1+ 1 u ∗ v ∗ + v ∗ w ∗ w ∗ + v ∗ + 1 u ∗ . (3:12) By the same argument with Case(6), we have X 2 ≤ 1+ε + ε 2 . (3:13) Case (8) (u ≥ 1, v ≥ 1, w ≥ 1, u + w ≥ 1). It here derives X ≤ x ∗ n = 1+u ∗ v ∗ + v ∗ w ∗ u ∗ + v ∗ + w ∗ ≤ 1+2(X + ε) 2 3 ( X + ε ) . (3:14) Hence X ≤ 1+ 9 4 ε 2 + ε 2 . (3:15) Collecting all above inequalities which imply X = 1 since ε > 0 is arbitrary, we com- plete the proof. □ 4 Generalization As mentioned above, the global asymptotic stability of positive solutions to the various equation listed above suggests that the same potentially holds for similar rational equa- tions. We can deduce the following natural generalization of (1.1) and (1.4). Corollary.Lets Î N + and Z s denote the set Z s = {1, 2, , s}. Suppose that {x i } satis- fies the form x n = s i=1 x n−k i s i =1 i = j , j ∈ Z s x n−k i x n−k j +1 , n =0,1, . (4:1) with initial value x -k , x -k+1 , , x -1 Î ℝ + , here k =max 1 ≤ i ≤ s {k i } . Then the sequence {x i } con- verges to the unique equilibrium 1. Remark. If we consider the equation which is added a constant a onto numerator and denominator of (4.1), the result is still viable. Indeed this corollary covers the results in [6]. Moreover, consulting the results of article [6,7,10], by the similar way to Lemma 2, we have the following generalization. Theorem 2.Suppose f (x n−k 1 , , x n−k r ) ∈ C(R r + , R + ), g(x n−m 1 , , x n−m s ) ∈ C(R s + , R + ) and h(x n−l 1 , , x n−l t ) ∈ C(R t + , R + ) satisfying [g(x n−m 1 , , x n−m s )] ∗ ≤ x ∗ n−m 1 . Then the equation x n+1 = f + g + h fg + g h +1 , n =0,1,2, . (4:2) Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 Page 6 of 8 with the corresponding positive initial data has a unique positive equilibrium ¯ x = 1 , and every solution of (4.2) converges to this point. Proof.Let {x n } ∞ n=− p be a solution sequence of equation (4.2) with initial data x -p , x -p+1 , x 0 Î ℝ + , where p =max{k r , m s , l t }. By the definition of x ∗ n , From the equation (4.2), the arguments in Lemma 2 and the hypothesis, it follows that for any n ≥ 0, 1 ≤ x ∗ n+1 = f + g + h fg + gh +1 ∗ ≤ [g(x n−m 1 , , x n−m s )] ∗ ≤ x ∗ n−m 1 , (4:3) from which we get that for any n ≥ 0and0≤ i ≤ m 1 , 1 ≤ x ∗ i+ ( n+1 )( m 1 +1 ) ≤ x ∗ i+n ( m 1 +1 ) . Hence the sequence {x ∗ i+n ( m 1 +1 ) } ∞ n= 0 with 0 ≤ i ≤ m 1 is convergent. Denote the limit as lim n →∞ x ∗ i+n(m 1 +1) = A i ,thenA i ≥ 1. Write M =max{A 0 , A 1 , , A m 1 } and A i+n ( m 1 +1 ) = A i for any integer n. Then there exists some 0 ≤ j ≤ m 1 such that lim n →∞ x ∗ j+n(m 1 +1) = M . From (4.3), it suggests M = g(M, A j −1−m 2 , , A j −1−m s )= M . Combining the facts 1 + ab ≥ a + b and ab+1+bc a + b + c ≤ ab+1+bc+abc a + b + c + ac , where a ≥ 1, b ≥ 1 and c ≥ 1, for the different situation in Theorem 1, we have x ∗ n+1 ≤ ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ f ∗ + g ∗ + f ∗ g ∗ h ∗ 1+f ∗ g ∗ + f ∗ h ∗ ≤ f ∗ + g ∗ + f ∗ g ∗ h ∗ f ∗ + g ∗ + f ∗ h ∗ , forCase (1,2,6,7), (4.4 ) f ∗ g ∗ +1+g ∗ h ∗ f ∗ + g ∗ + h ∗ ≤ f ∗ g ∗ +1+g ∗ h ∗ + f ∗ g ∗ h ∗ f ∗ + g ∗ + h ∗ + f ∗ h ∗ , for Case (3, 8), (4.5 ) f ∗ + h ∗ + f ∗ g ∗ h ∗ g ∗ h ∗ + f ∗ g ∗ + f ∗ h ∗ ≤ f ∗ + h ∗ + f ∗ g ∗ h ∗ f ∗ + h ∗ + f ∗ h ∗ , for Case (4, 5). (4.6 ) Therefore 1 ≤ M ≤ ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ f ∗ + M + f ∗ Mh ∗ f ∗ + M + f ∗ h ∗ , forCase(1,2,6,7), (4.7 ) f ∗ M +1+Mh ∗ + f ∗ Mh ∗ f ∗ + M + h ∗ + f ∗ h ∗ , for Case(3, 8), (4.8 ) f ∗ + h ∗ + f ∗ Mh ∗ f ∗ + h ∗ + f ∗ h ∗ , for Case(4, 5), (4.9 ) from which it follows M = 1. This implies A i =1for0≤ i ≤ m 1 and lim n →∞ x ∗ n = 1 . Since 1/x ∗ n ≤ x n ≤ x ∗ n , we obtain lim n →∞ x n = 1 . □ Remark. The stability of solut ion to equation (4.2) is ever p roposed to consider as a conjecture by K.S.Berenhaut etc. in [7]. Indeed, Theorem 1 proved the conjecture partially. In addition, gathering lots of relevant work listed in reference, we put forward the following conjecture. Conjecture. Let s Î N + , Z s = {1, 2, , s} and l ij ≥ 0. Suppose that {x i } satisfies x n = s j=1 i∈Z s x l ij n−k ij s−1 j=1 i∈Z s x l ij n−k i j +1 , n =0,1,2, . (4:10) with x -k , x -k+1 , , x -1 Î ℝ + , k =max i, j ∈Z s {k i,j } , then the sequence {x i } ∞ i = 0 converges to the unique equilibrium 1. Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 Page 7 of 8 Acknowledgements The authors would like to thank the reviewers and the editors for their valuable suggestions and comments; The authors wish to express their deep gratitude to Professor C.Y. Wang for his valuable advice and constant encouragement for this work supported in part by Natural Science Foundation for Colleges and Universities in Hebei Province(Z2011111, Z2011162) and Human Resources and Social Security Subject of Hebei Province(JRS-2011-1042). Authors’ contributions QS completed the main part of this paper, QX and GY corrected the main theorems. XL participated in the design and coordination. All authors read and approved the final manuscript. Competing interests The authors declare that they have no competing interests. Received: 17 May 2011 Accepted: 27 September 2011 Published: 27 September 2011 References 1. Li, W., Sun, H.: Global attractivity in a rational recursive sequence. Dyn. Syst. Appl. 11, 339–346 (2002) 2. Kocic, VL, Ladas, G: Global Behavior of Nonlinear Difference Equations of Higher Order with Application. Kiuwer Academic Publishers, Dordrecht. (1993) 3. Kulenovic, MR., Ladas, G: Dynamics of Second Order Rational Difference Equations with Open Problems and Conjectures. Chapman Hall/CRC, Boca Raton. (2001) 4. Çinar, C: On the positive solutions of the difference equation x n+1 = x n−1 1+x n x n−1 . Appl. Math. Comput. 150,21–24 (2004). doi:10.1016/S0096-3003(03)00194-2 5. Yang, X, Su, W, Chen, B, Megson, GM, Evans, DJ: On the recursive sequence. Appl. Math. Comput. 162, 1485–1497 (2005). doi:10.1016/j.amc.2004.03.023 6. Li, X: Global behavior for a fourth-order rational di_erenceequation. J. Math. Anal. Appl. 312, 103–111 (2005) 7. Berenhaut, KS, Foley, JD, Stevic, S: The global attractivity of the rational difference equation y n =(y n-k + y n-m )>/(1 + y n- k y n-m ). Appl. Math. Lett. 20,54–58 (2007). doi:10.1016/j.aml.2006.02.022 8. Papaschinopoulos, G, Schinas, CJ: Global asymptotic stability and oscillation of a family of difference equations. J. Math. Anal. Appl. 294, 614–620 (2004). doi:10.1016/j.jmaa.2004.02.039 9. Sun, T., Xi, H.: Global attractivity for a family of nonlinear difference equations. Appl. Math. Lett. 20, 741–745 (2007). doi:10.1016/j.aml.2006.08.024 10. Li, Z., Zhu, D.: Global asymptotic stability of a higher order nonlinear difference equation. Appl. Math. Lett. 19, 926–930 (2006). doi:10.1016/j.aml.2005.10.014 doi:10.1186/1687-1847-2011-36 Cite this article as: Shi et al.: Dynamic behavior of a nonlinear rational difference equation and generalization. Advances in Difference Equations 2011 2011:36. Submit your manuscript to a journal and benefi t from: 7 Convenient online submission 7 Rigorous peer review 7 Immediate publication on acceptance 7 Open access: articles freely available online 7 High visibility within the fi eld 7 Retaining the copyright to your article Submit your next manuscript at 7 springeropen.com Shi et al. Advances in Difference Equations 2011, 2011:36 http://www.advancesindifferenceequations.com/content/2011/1/36 Page 8 of 8 . RESEARCH Open Access Dynamic behavior of a nonlinear rational difference equation and generalization Qihong Shi * , Qian Xiao, Guoqiang Yuan and Xiaojun Liu * Correspondence:. Schinas, CJ: Global asymptotic stability and oscillation of a family of difference equations. J. Math. Anal. Appl. 294, 614–620 (2004). doi:10.1016/j.jmaa.2004.02.039 9. Sun, T., Xi, H.: Global attractivity. the main part of this paper, QX and GY corrected the main theorems. XL participated in the design and coordination. All authors read and approved the final manuscript. Competing interests The authors