DDNSDiscrete Dynamics in Nature and Society1607-887X1026-0226Hindawi Publishing Corporation61698210.1155/2009/616982616982Research ArticleGlobal Asymptotic Stability for a Fourth-Order Rational Difference EquationGülpinarMeseret Tuba1BayramMustafa2BravermanElena1Department of MathematicsFaculty of Arts and SciencesYıldız Technical University34210 EsenlerİstanbulTurkeyyildiz.edu.tr2Department of MathematicsFaculty of Arts and SciencesFatih University34500 BuyukcekmeceİstanbulTurkeyfatih.edu.tr20091407200920090805200904062009040620092009Copyright © 2009This 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.

Our aim is to investigate the global behavior of the following fourth-order rational difference equation: xn+1=(xnxn2xn3+xn+xn2+xn3+a)/(xnxn2+xnxn3+xn2xn3+1+a), n=0,1,2, where a[0,) and the initial values x3,x2,x1,x0(0,). To verify that the positive equilibrium point of the equation is globally asymptotically stable, we used the rule of the successive lengths of positive and negative semicycles of nontrivial solutions of the aforementioned equation.

1. Introduction

There has been a great interest in studying global behaviors of rational difference equations. One can easily see that it is hard to understand thoroughly the global behaviors of solutions of rational difference equations although they have simple forms. And there has not been any general method to identify the global behaviors of rational difference equations of order greater than one and so far .

Let us consider the following fourth-order difference equation: xn+1=xnxn-2xn-3+xn+xn-2+xn-3+axnxn-2+xnxn-3+xn-2xn-3+1+a  ,n=0,1,2,, where a[0,) and the initial values x-3,x-2,x-1,x0(0,) in this paper. By determining the rule for the positive and negative semicycles, we assigned the global behavior of the positive equilibrium point. The unique positive equilibrium point x̅ of (1.1) is obtained x̅=1 by solving x̅=x̅3+3x̅+a3x̅2+1+a.

Definition 1.1.

A solution {xn}n=-3 of (1.1) is said to be eventually trivial if xn is eventually equal to x̅=1; otherwise, the solution is said to be nontrivial .

Definition 1.2.

A solution {xn}n=-3 of (1.1) is said to be eventually positive if xn is eventually greater than x̅=1.

A solution {xn}n=-3 of (1.1) is said to be eventually negative if xn is eventually less than x̅=1 [1, 2, 46].

Definition 1.3.

A positive semicycle of a solution {xn}n=-3 of (1.1) consists of a “string” of terms {xl,xl+1,,xm}, all greater than or equal to the equilibrium point x̅, with l-3 and m< such that eitherl=-3orl>-3,xl-1<x̅,eitherm=orm<,xm+1<x̅. A negative semicycle of a solution {xn}n=-3 of (1.1) consists of a “string” of terms {xl,xl+1,,xm}, all greater than or equal to the equilibrium point x̅, with l-3 and m< such that eitherl=-3orl>-3,xl-1x̅,eitherm=orm<,xm+1x̅. And also the lengths of a semicycle is (m-l+1), the number of the terms contained in it. And we denote that the lengths of a positive semicycle are by (m-l+1)+ and the lengths of a negative semicycle are by (m-l+1)- [1, 2, 46].

Definition 1.4.

A solution {xn}n=-3 of (1.1) is called nonoscillatory about x̅, or simply nonoscillatory, if there exists N-3 such that either xnx̅,nN or xn<x̅,nN. Otherwise, the solution {xn}n=-3 is called oscillatory about x̅, or simply oscillatory [1, 2].

Lemma 1.5.

A positive solution {xn}n=-3 of (1.1) is eventually trivial if and only if (x-3-1)(x-2-1)(x-1-1)(x0-1)=0.

Proof.

To prove the lemma, first we assume that (x-3-1)(x-2-1)(x-1-1)(x0-1)0 and then we must show that xn1 for any n1.

Assume that for some N1, xN=1 and that xn1 for -3nN-1. So xN=xN-1xN-3xN-4+xN-1+xN-3+xN-4+axN-1xN-3+xN-1xN-4+xN-3xN-4+1+a=1, and we obtain (xN-4-1)(xN-3-1)(xN-1-1)=0; hence xN-4=1,xN-3=1,xN-1=1, when we solve the equation above. It is easy to see that xN-4=1, xN-3=1, or xN-1=1 contradicts with xn1 for -3nN-1.

If (1.7) holds, it is clear that the following conclusions hold:

if x-3=1, xn=1 for n1,

if x-2=1, xn=1 for n1,

if x-1=1, xn=1 for n2,

if x0=1, xn=1 for n1.

It is obvious that if the initial conditions do not satisfy (1.7), then the positive solution {xn}n=-3 of (1.1) is eventually nontrivial.

Lemma 1.6.

If {xn}n=-3 is a nontrivial positive solution of (1.1), then the following conclusions are satisfied:

(xn+1-1)(xn-1)(xn-2-1)(xn-3-1)>0,

(xn+1-xn)(xn-1)<0,

(xn+1-xn-2)(xn-2-1)<0,

(xn+1-xn-3)(xn-3-1)<0.

Proof.

(i) The proof of the inequality (i) is obtained by subtracting 1 from (1.1) xn+1-1=xnxn-2xn-3+xn+xn-2+xn-3+axnxn-2+xnxn-3+xn-2xn-3+1+a-1=(xn-1)(xn-2-1)(xn-3-1)xnxn-2+xnxn-3+xn-2xn-3+1+a. The dominator of this fraction is positive so (xn+1-1)(xn-1)(xn-2-1)(xn-3-1)>0 also.

(ii) If we subtract xn from (1.1), we obtain xn+1-xn=xnxn-2xn-3+xn+xn-2+xn-3+axnxn-2+xnxn-3+xn-2xn-3+1+a-xn=-(xn-1)[xn-2(xn+1)+xn-3(xn+1)+a]xnxn-2+xnxn-3+xn-2xn-3+1+a. The expression [xn-2(xn+1)+xn-3(xn+1)+a]/(xnxn-2+xnxn-3+xn-2xn-3+1+a) is positive, and so we get (xn+1-xn)(xn-1)<0. The proofs for inequalities (iii) and (iv) are similar to the one for (ii).

2. Main Results and Their Proofs

The trajectory of (1.1) and global asymptotic stability of the positive solution are considered in this part of the paper.

Theorem 2.1.

Let {xn}n=-3 be a strictly oscillatory solution of (1.1). Then the positive and negative semicycles of (1.1) are ,3+,3-,3+,3-,3+,3-,3+,3-, or ,2+,1-,2+,1-,2+,1-,2+,1-, or ,2-,1+,2-,1+,2-,1+,2-,1+, or ,1+,1-,1+,1-,1+,1-,1+,1-,.

Proof.

Assume that {xn}n=-3 is a strictly oscillatory solution of (1.1), then the initial values must satisfy one of the following four cases:

xp-3>1, xp-2<1, xp-1>1, xp>1,

xp-3>1, xp-2<1, xp-1>1, xp<1,

xp-3>1, xp-2<1, xp-1<1, xp>1,

xp-3>1, xp-2<1, xp-1<1, xp<1.

If (i) occurs, it follows from Lemma 1.6(i) that xp-3>1,xp-2<1,xp-1>1,xp>1,xp+1<1,  xp+2>1,xp+3>1,xp+4<1,  xp+5>1,xp+6>1,  xp+7<1,xp+8>1,  xp+9>1,xp+10<1,  xp+11>1,  xp+12>1,xp+13<1,xp+14>1,  xp+15>1,xp+16<1,  xp+17>1,  xp+18>1,xp+19<1,xp+20>1,xp+21>1,. It means that the rule for the lengths of positive and negative semicycles of the solution of (1.1) occurs successively as ,2+,1-,2+,1-,2+,1-,2+,1-,.

If (ii) happens, the positive and negative semicycles are xp-3>1,xp-2<1,  xp-1>1,  xp<1,xp+1>1,xp+2<1,  xp+3>1,xp+4<1,xp+5>1,  xp+6<1,xp+7>1,xp+8<1,xp+9>1,xp+10<1,xp+11>1,xp+12<1,xp+13>1,xp+14<1,xp+15>1,  xp+16<1,xp+17>1,xp+18<1,xp+19>1,  xp+20<1,xp+21>1,xp+22<1,xp+23>1,. The regulation for the lengths of positive and negative semicycles which occur successively is ,1+,1-,1+,1-,1+,1-,1+,1-,.

The other cases can be shown similarly.

Theorem 2.2.

The positive equilibrium point of (1.1) is globally asymptotically stable.

Proof.

Let us show that the positive equilibrium point x̅=1 of (1.1) is locally asymptotically stable and also globally attractive, then it is globally asymptotically stable. The linearized equation of (1.1) about the positive equilibrium x̅=1 is yn+1=q0yn+q1yn-1+q2yn-2+q3yn-3, where qi=(F/ui)(x̅,x̅,x̅,x̅) and F(u0,u1,u2,u3)=(u0u2u3+u0+u2+u3+a)/(u0u2+u0u3+u2u3+1+a).

And we obtain yn+1=q0yn+q1yn-1+q2yn-2+q3yn-3=0·yn+0·yn-1+0·yn-2+0·yn-3=0, thereby x̅=1 is locally asymptotically stable.

Now we must show that limnxn=x̅=1. The proof is as follows.

If initial values of (1.1) satisfy (1.7), then xn=1 according to Lemma 1.5, so limnxn=x̅=1.

If the initial values of (1.1) do not satisfy (1.7), then for any solution of (1.1), xn1 for n-3.

If the solution is nonoscillatory about the positive equilibrium point of (1.1), then {xn} is monotonic and bounded because of Lemma 1.6. So, the limit limnxn=L exists and is finite. If we take limits on both sides of (1.1), then we obtain L=(L3+3L+a)/(3L2+1+a) and thereby L=1. So limnxn=L=1.

If the solution is strictly oscillatory, then trajectory structure of nontrivial solutions of (1.1) is   ,3+,3-,3+,3-,3+,3-,3+,3-, or ,2+,1-,2+,1-,2+,1-,2+,1-, or ,2-,1+,2-,1+,2-,1+,2-,1+, or ,1+,1-,1+,1-,1+,1-,1+,1-,.

First, we investigate the case where the rule of the trajectory structure is ,3+,3-,3+,3-,3+,3-,3+,3-, in a period. We denote positive semicycles by {xp,xp+1,xp+2}+ and negative semicycles by {xp+3,xp+4,xp+5}-. The rule for the positive and negative semicycles can be “periodically’’ expressed as follows: {xp+6n,xp+6n+1,xp+6n+2}+,{xp+6n+3,xp+6n+4,xp+6n+5}-,n=0,1,. By using Lemma 1.6 we obtain

xp+6n+2<xp+6n+1<xp+6n; xp+6n+3<xp+6n+4<xp+6n+5,

xp+6n+6<xp+6n+2; xp+6n+9>xp+6n+5.

This relations give rise to xp+6n+6<xp+6n+2<xp+6n+1<xp+6n,xp+6n+3<xp+6n+4<xp+6n+5<xp+6n+9. It is easy to see that {xp+6n}n=0 is decreasing with its lower bound 1 because of the inequality xp+6n+6<xp+6n+2<xp+6n+1<xp+6n. Hence the limit exists and is finite. The limit is then limnxp+6n=limnxp+6n+1=limnxp+6n+2=L. Similarly {xp+6n+3}n=0 is increasing and its upper bound 1 because of the inequality xp+6n+3<xp+6n+4<xp+6n+5<xp+6n+9. So, the limit limnxp+6n+3=M exists and is finite. And we derive limnxp+6n+4=limnxp+6n+5=M.

Now, we must show that L=M=1. To do this, let us take xp+6n+6=xp+6n+5xp+6n+3xp+6n+2+xp+6n+5+xp+6n+3+xp+6n+2+axp+6n+5xp+6n+3+xp+6n+5xp+6n+2+xp+6n+3xp+6n+2+1+a    , and taking the limit on both sides of the above equality, we obtain L=M·M·L+M+M+L+aM·M+M·L+M·L+1+a. By solving this equation we have (L-1)[2M(L+1)+a]=0, and then L=1. Similarly to obtain M by using 1xp+6n+2<xp+6n+1 we are taking xp+6n+5=xp+6n+4xp+6n+2xp+6n+1+xp+6n+4+xp+6n+2+xp+6n+1+axp+6n+4xp+6n+2+xp+6n+4xp+6n+1+xp+6n+2xp+6n+1+1+a. And taking the limit on both sides of the above equality we get M=M·L·L+M+L+L+aM·L+M·L+L·L+1+a. By solving this equation we have (M-1)[2L(M+1)+a]=0 and M=1. Hence, we derive L=M=1, so limnxn=1. It can be shown that limnxn=1 for the other rules of the trajectory structures with the same manner. Therefore, the positive equilibrium point is globally asymptotically stable.

CamouzisE.LadasG.Dynamics of Third-Order Rational Difference Equations with Open Problems and Conjectures20085Boca Raton, Fla, USAChapman & Hall/CRCxxii+554Advances in Discrete Mathematics and ApplicationsMR2363297KulenovićM. R. S.LadasG.Dynamics of Second Order Rational Difference Equations with Open Problems and Conjectures2002Boca Raton, Fla, USAChapman & Hall/CRCxii+218MR1935074ÇınarC.StevićS.Yalçınkayaİ.A note on global asymptotic stability of a family of rational equationsRostocker Mathematisches Kolloquium2005594149MR2169497ZBL1083.39003LiX.Qualitative properties for a fourth-order rational difference equationJournal of Mathematical Analysis and Applications2005311110311110.1016/j.jmaa.2005.02.063MR2165465ZBL1082.39004LiX.Global behavior for a fourth-order rational difference equationJournal of Mathematical Analysis and Applications2005312255556310.1016/j.jmaa.2005.03.097MR2179097ZBL1083.39007LiD.lds1010@sina.comLiP.lipingping@ujs.edu.cnLiX.Dynamical properties for a class of fourth-order nonlinear difference equationsAdvances in Difference Equations2008200813678402EID2-s2.0-4194910797710.1155/2008/678402ZBL1146.39012