Stability analysis for a rational difference equation

Abstract Most qualitative behaviours of difference equations are rapidly investigated nowadays. This can be attributed to the fact that it is often sophisticated to construct the exact solutions of most difference equations. This article is written to analyse the local stability, global attractor and the boundedness of the solution of the seventh order difference equation given by xn+1=c1xn−1−c2xn−1xn−4c3xn−4+c4xn−6,n=0,1,…, where the coefficients ci,foralli=1,…,4, are supposed to be positive real numbers and the initial conditions xiforalli=−6,−5,…,0, are arbitrary non-zero real numbers. Under some suitable conditions, the stability, boundedness and a special case equation from the considered equation are presented in 2 D figures.


Introduction
The theory of difference equations has attained its popularity in recent decades due to its use in modelling most real life phenomena.It is worth mentioning that the difference equations are utilized in modelling some problems such as biological, physical, economical problems and several others.In order to understand the long behaviour of such models, one can study the qualitative and quantitative behaviours of these problem.This can be easily achieved by exploring the stability, periodicity and boundedness of the solutions.
Most researchers are mainly interested in studying difference equations of low order.As a result, simple analysis is required to investigate such equations.However, considering difference equations of high order leads to a very interesting, powerful and rigorous analysis.In this work, we present local and global stability for seventh order difference equations which can be generalized for other equations with higher orders.Indeed, such equations appear naturally as discrete analogues for differential equations which model a massive number of natural phenomena (Agarwal, 1992;Fisher, 1984;Kocic & Ladas, 1993).Thus, we are excited to study such proposed equations.

Local stability of the equilibrium point
The local stability in the neighbourhood of the fixed point is extensively highlighted in this section.We provide special conditions under which the equilibrium point is locally asymptotically stable.The equilibrium point is given by then the unique equilibrium point is given by x ¼ 0: We now introduce a new function to obtain the stability of Equation (2).Define a function g : ð0, 1Þ 3 !ð0, 1Þ by gð/, w, hÞ Hence, @gð/, w, hÞ @gð/, w, hÞ @gð/, w, hÞ The following step is embodied in evaluating Equations ( 4), ( 5) and Equation ( 6) at x: That is Thus, the linearized equation of Equation ( 2) around x can be presented as follows: Then, the equilibrium point of Equation ( 2) is locally asymptotically stable.
Proof.According to Theorem A in (Elabbasy et al., 2008), the local stability of the equilibrium point occurs if Substituting which can be reduced to Then, the equilibrium point of Equation ( 2) is locally asymptotically stable.
Proof.The proof is similar to the proof of Theorem (1).Hence, it omitted.

Analysis of global stability
This section is included to illustrate the global attractivity of the equilibrium point.The investigation of this part is established by employing Theorem B in (Elsayed, 2010).
Theorem 3. Assume that c 1 > c 2 w c 3 wþc 4 h : Then, the equilibrium point of Equation ( 2) is a global attractor if c 1 <1: Proof.Let r, s 2 R and assume that g : ½r, s 3 !½r, s is a function defined by Equation (3).Then, it can be apparently noted from Equations ( 4), ( 5) and Equation ( 6) that g is increasing in / and in h and decreasing in w: Theorem B in (Elsayed, 2010) requires a solution (say ðf, gÞ) to the following system: which can be expanded as follows: Now, we subtract Equation (10) from Equation ( 9) to obtain which can be rearranged as Therefore, if c 1 <1, then f ¼ g: As a result, Theorem B (Elsayed, 2010) concludes that the equilibrium point is a global attractor.
Theorem 4. Assume that c 1 < c 2 w c 3 wþc 4 h : Then, the equilibrium point of Equation ( 2 Proof.The proof is similar to the previous proof.Thus, it is omitted.

Exact solution to
This section provides sophisticated forms of solutions to the following rational recursive relation: where the initial data as shown previously.
Theorem 5. Let fx n g 1 n¼À6 be a solution to Equation (11) and assume that x Proof.It is simply shown that the solutions are true for n ¼ 0: Suppose that nS> 0 and that our assumption holds fornÀ1: That is, We now turn to prove some formulae.It can be easily obtained from Equation (11) that it can be concluded that  x 10nÀ6 ¼ a Another formula can be proved as follows.
The rest of the relations can be similarly confirmed.

Confirmation of the results
Our main goal in this part is to verify the appropriate and vital results obtained in this article.

Conclusion
In brief, we have shown the stability of Equation ( 2) under some powerful and effective hypotheses.In particular, the local stability of the fixed point occurs if c 2 ðc 4 Àc 3 Þ<ð1Àc 1 Þðc 3 þ c 4 Þ 2 , and the global behaviour of the equilibrium point takes place if c 1 <1, as can be observed in Figure 1 and Figure 2, respectively.Other conditions for the local and global stability are discovered and illustrated.In Section 4, the exact solution of Equation ( 11) has been perfectly determined using Fibonacci sequence.Simsek, D., Cinar, C., & Yalcinkaya, I. (2006).On the recursive sequence x nþ1 ¼ xnÀ3 1þxnÀ1 : International Journal of Contemporary Mathematical Sciences, 1, 475-480. doi:10. 12988/ijcms.2006.06052

Figure 2 .
Figure 2. Global nature about the equilibrium.