Abstract and Applied Analysis

Volume 2013 (2013), Article ID 179423, 10 pages

http://dx.doi.org/10.1155/2013/179423

## On the Period-Two Cycles of

^{1}School of Mathematical Sciences, The National of University of Malaysia, 43600 Bangi, Selangor, Malaysia^{2}Department of Basic Sciences, King Saud bin Abdulaziz University for Health Sciences, Riyadh 11426, Saudi Arabia

Received 2 January 2013; Revised 13 April 2013; Accepted 13 April 2013

Academic Editor: Douglas Anderson

Copyright © 2013 S. Atawna 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.

#### Abstract

We consider the higher order nonlinear rational difference equation , where the parameters are positive real numbers and the initial conditions are nonnegative real numbers, . We give a necessary and sufficient condition for the equation to have a prime period-two solution. We show that the period-two solution of the equation is locally asymptotically stable.

#### 1. Introduction

Recently, dynamics of nonnegative solutions of higher order rational difference equation has been an area of intense interest. Related to this subject, researches are done by Dehghan et al. [1–4], Zayed [5–7], Huang and Knopf [8, 9], Karatas [10, 11], and others. For the general theory of difference equations, one can refer to the monographes of Kocic and Ladas [12], Elaydi [13], Agarwal [14], Kulenovic and Ladas [15], and Camouzis and Ladas [16]. Other related results can be found in [17–24].

Our aim in this paper is to study the higher order nonlinear rational difference equation where the parameters are positive real numbers and the initial conditions are nonnegative real numbers, . Our concentration is on the periodic character of all positive solutions of (1).

The periodic character of positive solutions of (1) for has been investigated by the authors in [25]. They showed that the period-two solution of (1) for is locally asymptotically stable if it exists.

Motivated by the above results, our interest is now to study and generalize the previous results to the general case depicted in (1).

The change of variable reduces (1) to where are positive real numbers and the initial conditions are nonnegative real numbers.

This paper is organized besides this introduction in three sections. In Section 2, we present some preliminaries and some results which can be mainly deduced from the general situation studied in [12–16, 26]. Our main results are presented in Section 3; we give a necessary and sufficient condition for the equation to have a prime period-two solution, in addition to providing a necessary and sufficient conditions for the prime period-two solution of the equation to be locally asymptotically stable. In order to illustrate the results of the previous section and to support our theoretical discussion, we consider several numerical examples in Section 4; we use MATLAB to see how the behaviors of (1) look like. Finally, we conclude in Section 5 with suggestions for future research.

#### 2. Preliminaries

For the sake of self-containment and convenience, we recall the following definitions and results from [16].

Let be some interval of real numbers and let be a continuously differentiable function. Then for every set of initial conditions , the difference equation has a unique solution .

A solution of (6) that is constant for all is called an equilibrium solution of (6). If is an equilibrium solution of (6), then is called an equilibrium point or simply an equilibrium of (6).

*Definition 1. *(i) The equilibrium point of (6) is called locally stable if for every , there exists such that if , and , we have

(ii) The equilibrium point of (6) is called locally asymptotically stable if it is locally stable, and if there exists , if , and , we have

(iii) The equilibrium point of (6) is called a global attractor if for every , we have

(iv) The equilibrium point of (6) is called globally asymptotically stable if it is locally stable and a global attractor.

(v) The equilibrium point of (6) is called unstable if it is not stable.

*Definition 2. *(i) A solution of (6) is said to be periodic with period if

(ii) A solution of (6) is said to be periodic with prime period if is the least positive integer for which (11) holds.

*Definition 3. *Let
denote the partial derivatives of evaluated at the equilibrium of (6). Then the equation
is called the linearized equation associated with (6) about the equilibrium point . Its characteristic equation is

Theorem 4 (linearized stability). *(a) If all roots of (14) lie in the open unit disk , then the equilibrium of (6) is locally asymptotically stable. **(b) If at least one of the roots of (14) has absolute value greater than one, then is unstable. *

The following result from [26] will become handy in the sequel.

Lemma 5. *If
**
then satisfies the following recursive formula:
*

The aforementioned lemma leads to the following conclusion.

Corollary 6. *If
**
then
*

#### 3. Main Result

In this section, we give a necessary and sufficient condition for (1) to have a prime period-two solution. We show that the period-two solution of (1) is locally asymptotically stable.

Equation (3) has a unique positive equilibrium given by The linearized equation associated with (3) about the equilibrium is given by and its characteristic equation is

Theorem 7. *(a) If
**
?then (3) has no nonnegative prime period-two solution. **(b) If
**
?then (3) has prime period-two solution
**
?if and only if is odd and
**
?where the values of and are the positive and distinct solutions of the quadratic equation
*

*Proof. *Assume that there exist distinct nonnegative real numbers and , such that
is a prime period-two solution of (3); there are two cases to be considered.*Case??1??(k is even).* In this case and satisfy
Furthermore,
Subtracting (30) from (29), we have
so
This contradicts the hypothesis that and are distinct nonnegative real numbers. Also, contradicts the hypothesis that and are positive real numbers.*Case??2??(k is odd).* (a) If
then in this case and satisfy
Furthermore,
Subtracting (35) from (36), we have
But, , this implies that which contradicts the hypothesis that are distinct positive real numbers.

(b) If
then in this case and satisfy
Moreover,
Subtracting (41) from (42), we have
Furthermore, one adds (41) to (42), makes use of (43), and then does some elementary algebraic manipulation; we have
Equation (44) leads to the following conclusion:
that follows from the facts that
Notice that when then adding (41) to (42) gives , which is impossible.

Construct the quadratic equation
So and are the positive and distinct solutions of the above quadratic equation, that is,

Theorem 8. *Suppose (3) has a prime period-two solution. Then, the period-two solution is locally asymptotically stable. *

*Proof. *To investigate the local stability of the two cycles
we first vectorize (3) by introducing the following change of variables:
and write (3) in the equivalent form:
where
Now and generate a period-two solution of (3) only if
is a fixed point of , the second iterate of . Furthermore,
where
The prime period-two solution of (3) is asymptotically stable if the eigenvalues of the Jacobian matrix , evaluated at lie inside the unit disk.

We have

Now let be the characteristic polynomial of . Then, by the Laplace expansion in the row,

However; by Lemma 5, Corollary 6, and the fact that is odd,
Therefore,

Hence, the characteristic polynomial is given by
where

Assume that . Then, by (39),
Hence,
Similarly, we observe that
Furthermore, since , (43) implies the sum of , is less than and, a fortiori, each is less than 1. Indeed, we have

With that in mind, it is clear that
In addition, with understanding that
and the fact that
we have to establish
First, we will establish inequality (69). To this end, observe that inequality (69) is equivalent to
which is true if and only if
which is true if and only if
which is true if and only if
which is true if and only if
which is true if and only if
which is clearly satisfied.

Next we will establish inequality (70). Observe that inequality (70) is equivalent to
which is true if and only if
which is true if and only if
which is true if and only if
which is true if and only if
which is true if and only if
Now observe that the righthand side of (82) is
The lefthand side of (82) is
Hence, inequality (70) is true if and only if
or equivalently
which is clearly satisfied (condition (25)).

Now, by applying Theorem 4 we shall show that the zeros of in (60) lie in the open unit disk . To do so, suppose to the contrary that has a zero such that . Then, by the triangle inequality,
Thus,

However, by the Descartes’ Rule of Signs has either two or no positive zeros. Furthermore,
and so, by the Intermediate Value Theorem, has two positive zeros in the open interval . Moreover, since , we conclude that for all which contradicts inequality (88).

The proof is complete.

*Remark 9. *The characteristic equation of the linearized equation at the equilibrium solution is given by
Since the magnitude of the constant term is less than , the equation has at least one root inside the unit disk. As such, by the Stable Manifold Theorem, there is a manifold of solutions, of dimension bigger than or equal to 1, that converge to the equilibrium solution. Hence, the period-two solution cannot be globally asymptotically stable.

#### 4. Numerical Examples

In order to illustrate the results of the previous section and to support our theoretical discussion, we consider several numerical examples generated by MATLAB.

*Case??1??(k is even). *For this case we consider the following example:
The dynamics of (91) is shown in Figure 1, no prime period-two solution.

*Case??2??(k is odd)*. There are two cases to be considered.

*Subcase??2.1?? **.* For this case we consider the following example:
The dynamics of (92) is shown in Figure 2, no prime period-two solution.

*Subcase??2.2?? ** and **.* For this case we consider the following example:
The dynamics of (93) is shown in Figure 3; it has prime period-two solution.

#### 5. Conclusion

In this paper, we showed that the period-two solution of the higher order nonlinear rational difference equation where the parameters are positive real numbers and the initial conditions are nonnegative real numbers, , is locally asymptotically stable if it exists.

We consider the aforementioned result as a step forward in investigating bigger classes of difference equations which afford the ELAS property; that is, the existence of a periodic solution implies its local asymptotic stability.

#### Acknowledgment

The authors are very grateful to the anonymous referees for carefully reading the paper and for their comments and valuable suggestions that lead to an improvement in the paper.

#### References

- M. Dehghan, M. J. Douraki, and M. Razzaghi, “Global stability of a higher order rational recursive sequence,”
*Applied Mathematics and Computation*, vol. 179, no. 1, pp. 161–174, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Dehghan and R. Mazrooei-Sebdani, “The characteristics of a higher-order rational difference equation,”
*Applied Mathematics and Computation*, vol. 182, no. 1, pp. 521–528, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Dehghan and R. Mazrooei-Sebdani, “Dynamics of a higher-order rational difference equation,”
*Applied Mathematics and Computation*, vol. 178, no. 2, pp. 345–354, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Dehghan and N. Rastegar, “On the global behavior of a high-order rational difference equation,”
*Computer Physics Communications*, vol. 180, no. 6, pp. 873–878, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. M. E. Zayed, “Dynamics of the rational difference equation ${x}_{n+1}=(a{x}_{n}+b{x}_{n-k}+a{x}_{n}+b{x}_{n-k})/(c+d{x}_{n}{x}_{n-k})$,”
*Communications on Applied Nonlinear Analysis*, vol. 18, no. 4, pp. 39–56, 2011. View at Google Scholar · View at MathSciNet - E. M. E. Zayed and M. A. El-Moneam, “On the rational recursive sequence ${x}_{n+1}=a{x}_{n}+(\beta {x}_{n}+\gamma {x}_{n-k})/(b{x}_{n}{x}_{n+k})$,”
*Communications on Applied Nonlinear Analysis*, vol. 16, no. 3, pp. 91–106, 2009. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. M. E. Zayed and M. A. El-Moneam, “On the rational recursive sequence ${x}_{n+1}=(a{x}_{n}+b{x}_{n-k}+\beta {x}_{n}+\gamma {x}_{n-k})/(C{x}_{n-k}+d{x}_{n-k})$,”
*Acta Applicandae Mathematicae*, vol. 111, no. 3, pp. 287–301, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. S. and P. M. Knopf, “Boundedness of positive solutions of second-order rational difference equations,”
*Journal of Difference Equations and Applications*, vol. 10, no. 11, pp. 935–940, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. S. Huang and P. M. Knopf, “Boundedness and some convergence properties of the difference equation ${x}_{n+1}=(\gamma {x}_{n-1}+\delta {x}_{n-2})/(B{x}_{n}+D{x}_{n2})$,”
*Journal of Difference Equations and Applications*, vol. 18, no. 1, pp. 27–55, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Karatas, “Global behavior of a higher order difference equation,”
*Computers & Mathematics with Applications*, vol. 60, no. 3, pp. 830–839, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Karatas and A. Gelişken, “Qualitative behavior of a rational difference equation,”
*Ars Combinatoria*, vol. 100, pp. 321–326, 2011. View at Google Scholar · View at MathSciNet - V. L. Kocić and G. Ladas,
*Global Behavior of Nonlinear Difference Equations of Higher Order with Applications*, vol. 256, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1993. View at MathSciNet - S. N. Elaydi,
*An Introduction to Difference Equations*, Springer, New York, NY, USA, 2nd edition, 1999. View at MathSciNet - R. P. Agarwal,
*Difference Equations and Inequalities: Theory, Methods, and Applications*, vol. 228, Marcel Dekker, New York, NY, USA, 2nd edition, 2000. View at MathSciNet - M. R. S. Kulenović and G. Ladas,
*Dynamics of second-Order Rational Difference Equations with Open Problem and Conjectures*, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2002. View at MathSciNet - E. Camouzis and G. Ladas,
*Dynamics of Third-Order Rational Difference Equations with Open Problem and Conjectures*, vol. 5, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2008. View at MathSciNet - R. Abu-Saris, C. Çinar, and I. Yalçinkaya, “On the asymptotic stability of ${x}_{n+1}=(a+{x}_{n}{x}_{n-k})/({x}_{n}+{x}_{n-k})$,”
*Computers & Mathematics with Applications*, vol. 56, no. 5, pp. 1172–1175, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. S. Berenhaut and S. Stević, “The global attractivity of a higher order rational difference equation,”
*Journal of Mathematical Analysis and Applications*, vol. 326, no. 2, pp. 940–944, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y.-H. Su and W.-T. Li, “Global attractivity of a higher order nonlinear difference equation,”
*Journal of Difference Equations and Applications*, vol. 11, no. 10, pp. 947–958, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Sun and H. Xi, “Global asymptotic stability of a higher order rational difference equation,”
*Journal of Mathematical Analysis and Applications*, vol. 330, no. 1, pp. 462–466, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Wang, “On the dynamics of ${x}_{n+1}=(\beta {x}_{n}+\gamma {x}_{n-k})/(b{x}_{n}+C{x}_{n-k}+\alpha )$,”
*Journal of Difference Equations and Applications*, vol. 15, no. 10, pp. 949–961, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Xi and T. Sun, “Global behavior of a higher-order rational difference equation,”
*Advances in Difference Equations*, Article ID 27637, 7 pages, 2006. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X.-X. Yan, W.-T. Li, and Z. Zhao, “Global asymptotic stability for a higher order nonlinear rational difference equations,”
*Applied Mathematics and Computation*, vol. 182, no. 2, pp. 1819–1831, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Yang and F. B. Lv, “Dynamics of a higher order rational difference equation,”
*Advanced Materials Research*, vol. 216, pp. 50–55, 2011. View at Publisher · View at Google Scholar - S. Atawna, R. Abu-Saris, and I. Hashim, “Local stability of period two cycles of second order rational difference equation,”
*Discrete Dynamics in Nature and Society*, vol. 2012, Article ID 969813, 11 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. W. Hager,
*Applied Numerical Linear Algebra*, Prentice-Hall International Editions, 1988.