Qualitative analysis for two fractional difference equations

Mohammed B. Almatrafi; Marwa M. Alzubaidi

doi:10.1515/nleng-2020-0014

Article Open Access

Qualitative analysis for two fractional difference equations

Mohammed B. Almatrafi and Marwa M. Alzubaidi

Published/Copyright: May 13, 2020

Published by

Become an author with De Gruyter Brill

Submit Manuscript Author Information

From the journal Nonlinear Engineering Volume 9 Issue 1

Abstract

Some difference equations are generally studied by investigating their long behaviours rather than their exact solutions. The proposed equations cannot be solved analytically. Hence, this article discusses the main qualitative behaviours of two rational difference equations. Some appropriate hypotheses are examined and given to show the local and global attractivity. Special cases from the considered equations are solved analytically. The periodicity is also proved in this work. We also illustrate the achieved results in some 2D figures.

Keywords: attractivity; boundedness; difference equations; equilibrium; periodicity; stability

1 Introduction

Difference equations are widely utilized in describing some phenomena occurred in nonlinear sciences such as biology, physics, chemical reactions, economy, probability theory, population growth, genetics, computer science and so on. They are sometimes used to discretise partial and ordinary derivatives appeared in partial and ordinary differential equations, respectively. Such equations have gained their popularity in recent years due to their use in nonlinear problems. It is well known fact that the analytical solutions of most fractional recursive relations cannot be often constructed due to the lack of the relevant mathematical methods. This has prevented some experts from searching the future pattern of most natural situations described by difference equations. Therefore, the majority of researchers have made great efforts to study these phenomena by means of some auxiliary elements such as local stability, global character, periodicity, boundedness and several others. For instance, El-Dessoky et al. [1] investigated the dynamics and periodicity of the fifth order difference equation

yn+1=αyn+βynyn−3Ayn−4+Byn−3.

Almatrafi and Alzubaidi [2] analysed the asymptotic attractivity, periodic nature and boundedness of the eighth order difference equation

xn+1=c1xn−3+c2xn−3c3xn−3−c4xn−7.

Moreover, Kalabušić et al. [3] studied the global character of the second order equation

xn+1=xnxn−1+αxn+βxn−1axnxn−1+bxn−1.

The author in [4] highlighted the semi-cycle behaviour, stability character and solutions for the second order difference equation

xn+1=xn−1a−xnxn−1.

El-Owaidy et al. [5] explored the global stability of the difference equation

xn+1=αxn−1β+yxn+1p.

In [6], the author utilized Fibonacci sequence to construct the solutions of special type of Riccati difference equations given by

xn+1=11+xn,yn+1=1yn−1.

Almatrafi [7] presented forms of exact solutions to the fractional system

xn+1=xn−1yn−3yn−1−1−xn−1yn−3,yn+1=yn−1xn−3xn−1±1±yn−1xn−3,n=0,1,…,

In order to obtain more qualitative results on difference equations, see refs.[8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28].

The basic task of this work is to point out the local and global stability, periodicity and forms of solutions to the following two recursive equations:

um+1=aum−1+bum−1um−4cum−4−dum−6,m=0,1,…,um+1=aum−1−bum−1um−4cum−4−dum−6,m=0,1,…,

where the coefficients a, b, c and d are supposed to be positive real numbers and the initial data u_ifor alli = −6, −5, …, 0, are arbitrary non-zero real numbers.

2 Analysis of first equation

This section is devoted to study the difference equation given by

um+1=aum−1+bum−1um−4cum−4−dum−6,m=0,1,…,(1)

where the constants a, b, c and d are assumed to be positive real numbers and the initial conditions u_ifor alli = −6, −5, …, 0, are arbitrary non-zero real numbers.

2.1 Equilibrium and local stability

The discussion in this part concentrates on the local behaviour around the equilibrium point. The equilibrium point is obtained as follows:

u¯=au¯+bu¯2cu¯−du¯,

which leads to a unique equilibrium point given by u = 0, if b ≠ (1 − a)(c − d). In order to linearise about the fixed point, we define a function g : (0, ∞)³ ⟶ (0, ∞) by

g(x,y,z)=ax+bxycy−dz.(2)

Hence,

∂g(x,y,z)∂x=a+by(cy−dz),(3)

∂g(x,y,z)∂y=−bdxz(cy−dz)2,(4)

∂g(x,y,z)∂z=bdxy(cy−dz)2.(5)

Evaluating Eqs. (3), (4) and Eq. (5) at ū gives us

∂g(u¯,u¯,u¯)∂x=a+bu¯(cu¯−du¯)=a+bc−d:=−q1,∂g(u¯,u¯,u¯)∂y=−bdu¯2(cu¯−du¯)2=−bd(c−d)2:=−q2,∂g(u¯,u¯,u¯)∂z=bdu¯2(cu¯−du¯)2=bd(c−d)2:=−q3.

Hence, the linearised equation of Eq. (1) around ū is shown as follows:

Un+1+q1Un−1+q2Un−4+q3Un−6=0.

Theorem 1

Assume that

a+bc−d<1−2bd(c−d)2.(6)

Then, the equilibrium point of Eq.(1) is locally asymptotically stable.

Proof

The proof is established according to the hypotheses of Theorem A in [13]. Therefore, the local stability occurs if

q1+q2+q3<1,(7)

which implies that

−(a+bc−d)+|bd(c−d)2|+−(bd(c−d)2)<1.(8)

Inequality (8) can be simply reduced to

a+bc−d<1−2bd(c−d)2,

which is the required condition.

2.2 Global attractivity

The global behaviour in the neighbourhood of the equilibrium point is explained in this subsection according to the conditions of Theorem C in [29].

Theorem 2

Leta+bc−d>0.Then, the equilibrium point of Eq.(1) is a global attractor if a ≠ 1.

Proof

Let r₁, r₂ ∈ ℝ, assume that g : [r₁, r₂]³ ⟶ [r₁, r₂] is a function defined by Eq. (2), and assume that a+bc−d>0. Then, it can be easily observed from Eqs. (3), (4) and Eq. (5) that g is increasing in x and in z and decreasing in y. Suppose that (ζ, η) is a solution to the following system:

ζ=g(ζ,η,ζ),η=g(η,ζ,η).

Or,

ζ=g(ζ,η,ζ)=aζ+bζηcη−dζ,η=g(η,ζ,η)=aη+bηζcζ−dη.

Simplifying this gives

cζη−dζ2=acζη−adζ2+bζη,(9)

cζη−dη2=acζη−adη2+bζη.(10)

Subtracting Eq.(10) from Eq.(9) yields

d(η2−ζ2)=ad(η2−ζ2),

Hence, if a ≠ 1, then ζ = η. Consequently, Theorem C [29] concludes that every solution of Eq. (1) converges to ū.

Theorem 3

Assume thata+bc−d<0.Then, the equilibrium point of Eq.(1) is a global attractor ifd ≠ ac + b.

Proof

The proof is similar to the previous proof. Thus, it is omitted.

2.3 Exact solution of um+1=um−1+um−1um−4um−4−um−6

We now begin with presenting the exact solution of the equation

um+1=um−1+um−1um−4um−4−um−6.(11)

The coefficients a, b, c and d are positive real numbers and the initial data u_ifor alli = −6, −5, …, 0, are arbitrary non-zero real numbers.

Theorem 4

Let{um}m=−6∞be a solution to Eq. (11)and suppose thatu₋₆ = a, u₋₅ = b, u₋₄ = c, u₋₃ = d, u₋₂ = e, u₋₁ = r, u₀ = k. Then, form = 0, 1, …, we have

u10m−6=(F2m+1d+F2m−1b)(F2m+1r+F2m−1d)(F2mk+F2m−2e)(F2mc+F2m−2a)(F2me+F2m−2c)ec(d−b)(r−d),u10m−5=(F2m+1k+F2m−1e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)(F2md+F2m−2b)(F2mr+F2m−2d)d(c−a)(e−c)(k−e),u10m−4=(F2m+1d+F2m−1b)(F2m+1r+F2m−1d)(F2mk+F2m−2e)(F2me+F2m−2c)(F2m+2c+F2ma)ec(d−b)(r−d),u10m−3=(F2m+2d+F2mb)(F2mr+F2m−2d)(F2m+1k+F2m−1e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)d(c−a)(e−c)(k−e),u10m−2=(F2m+2e+F2mc)(F2m+2c+F2ma)(F2m+1d+F2m−1b)(F2m+1r+F2m−1d)(F2mk+F2m−2e)ec(d−b)(r−d),u10m−1=(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m−1e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)d(c−a)(e−c)(k−e),u10m=(F2m+2k+F2me)(F2m+2e+F2mc)(F2m+2c+F2ma)(F2m+1d+F2m−1b)(F2m+1r+F2m−1d)ec(d−b)(r−d),u10m+1=(F2m+3c+F2m+1a)(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m−1e)(F2m+1e+F2m−1c)d(c−a)(e−c)(k−e),u10m+2=(F2m+3d+F2m+1b)(F2m+1r+F2m−1d)(F2m+2k+F2me)(F2m+2e+F2mc)(F2m+2c+F2ma)ec(d−b)(r−d),u10m+3=(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m−1e)(F2m+3e+F2m+1c)(F2m+3c+F2m+1a)d(c−a)(e−c)(k−e).

Here, {Fm}m=0∞ = {1, 1, 2, 3, 5, 8, … }, i.e, F_m = F_m−1 + F_m−2, F₋₁ = 0 and F₋₂ = −1, is called Fibonacci sequence.

Proof

The relations are true for m = 0. Now, assume that m > 0 and that our assumption holds for m − 1. That is,

u10m−16=(F2m−1d+F2m−3b)(F2m−1r+F2m−3d)(F2m−2k+F2m−4e)(F2m−2c+F2m−4a)(F2m−2e+F2m−4c)ec(d−b)(r−d),u10m−15=(F2m−1k+F2m−3e)(F2m−1e+F2m−3c)(F2m−1c+F2m−3a)(F2m−2d+F2m−4b)(F2m−2r+F2m−4d)d(c−a)(e−c)(k−e),u10m−14=(F2m−1d+F2m−3b)(F2m−1r+F2m−3d)(F2m−2k+F2m−4e)(F2m−2e+F2m−4c)(F2m−2c+F2m−4a)ec(d−b)(r−d),u10m−13=(F2md+F2m−2b)(F2m−2r+F2m−4d)(F2m−1k+F2m−3e)(F2m−1e+F2m−3c)(F2m−1c+F2m−3a)d(c−a)(e−c)(k−e),u10m−12=(F2me+F2m−2c)(F2mc+F2m−2a)(F2m−1d+F2m−3b)(F2m−1r+F2m−3d)(F2m−2k+F2m−4e)ec(d−b)(r−d),u10m−11=(F2md+F2m−2b)(F2mr+F2m−2d)(F2m−1k+F2m−3e)(F2m−1e+F2m−3c)(F2m−1c+F2m−3a)d(c−a)(e−c)(k−e),u10m−10=(F2mk+F2m−2e)(F2me+F2m−2c)(F2mc+F2m−2a)(F2m−1d+F2m−3b)(F2m−1r+F2m−3d)ec(d−b)(r−d),u10m−9=(F2m+1c+F2m−1a)(F2md+F2m−2b)(F2mr+F2m−2d)(F2m−1k+F2m−3e)(F2m−1e+F2m−3c)d(c−a)(e−c)(k−e),u10m−8=(F2m+1d+F2m−1b)(F2m−1r+F2m−3d)(F2mk+F2m−2e)(F2me+F2m−2c)(F2mc+F2m−2a)ec(d−b)(r−d),u10m−7=(F2md+F2m−2b)(F2mr+F2m−2d)(F2m−1k+F2m−3e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)d(c−a)(e−c)(k−e).

Two random relations will be now selected and proved. It can be easily observed from Eq. (11) that

u10m−5=u10m−7+y10m−7y10m−10y10m−10−y10m−12=u10m−71+y10m−10y10m−10−y10m−12=u10m−71+F2mk−F2m−2eF2mk−F1m−2e−F2m−2k+F2m−4e=u10m−71+F2mk−F2m−2eF2m−1k−F2m−3e=u10m−7F2m+1k−F2m−1eF2m−1k−F2m−3e=(F2md+F2m−2b)(F2mr+F2m−2d)(F2m−1k+F2m−3e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)d(c−a)(e−c)(k−e)F2m+1k−F2m−1eF2m−1k−F2m−3e=(F2m+1k+F2m−1e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)(F2md+F2m−2b)(F2mr+F2m−2d)d(c−a)(e−c)(k−e).

We can also see from Eq. (11) that

u10m−1=u10m−3+u10m−3u10m−6u10m−6−u10m−8=u10m−31+u10m−6u10m−6−u10m−8=u10m−31+F2m+1r−F2m−1dF2m+1r−F1m−1d−F2m−1r+F2m−3d=u10m−31+F2m+1r−F2m−1dF2mr−F2m−2d=u10m−3F2m+2r−F2mdF2mr−F2m−2d=(F2m+2d+F2mb)(F2mr+F2m−2d)(F2m+1k+F2m−1e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)d(c−a)(e−c)(k−e)F2m+2r−F2mdF2mr−F2m−2d=(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m−1e)(F2m+1e+F2m−1c)(F2m+1c+F2m−1a)d(c−a)(e−c)(k−e).

Similarly, other relations can be proved.

3 Analysis of second equation

In this section, we are interested in analysing the periodicity of the recursive equation

um+1=um−1−um−1um−4um−4−um−6,m=0,1,…,(12)

where the constants a, b, c and d are positive real numbers and the initial data u_ifor alli = −6, −5, …, 0, are arbitrary non-zero real numbers.

Theorem 5

Let{um}m=−6∞be a solution to Eq. (12)whereu_k ≠ u_k−2for allk = −4, −3, −2, … . Then, every solution of Eq. (12)is periodic with period 60. Furthermore, {um}m=−6∞takes the form

a,b,c,d,e,r,k,ara−c,bkb−d,acr(a−c)(c−e),bdk(b−d)(d−r),acer(a−c)(c−e)(e−k),−bdk(a−c)c(b−d)(d−r),−acer(b−d)d(a−c)(c−e)(e−k),bdk(a−c)(c−e)ce(b−d)(d−r),ace(b−d)(d−r)d(a−c)(c−e)(e−k),−bd(a−c)(c−e)(e−k)ce(b−d)(d−r),c^2e(b−d)(d−r)d(a−c)(c−e)(e−k),−d^2(a−c)(c−e)(e−k)ce(b−d)(d−r),ce^2(b−d)(d−r)d(a−c)(c−e)(e−k),−dr(a−c)(c−e)(e−k)ce(b−d)(d−r),cek(b−d)(d−r)d(a−c)(c−e)(e−k),−adr(c−e)(e−k)ce(b−d)(d−r),bcek(d−r)d(a−c)(c−e)(e−k),−adr(e−k)e(b−d)(d−r),bcek(a−c)(c−e)(e−k),−adr(b−d)(d−r),−bek(c−e)(e−k),ard−r,bke−k,−a,−b,−c,−d,−e,−r,−k,−ara−c,−bkb−d,−acr(a−c)(c−e),−bdk(b−d)(d−r),−acer(a−c)(c−e)(e−k),bdk(a−c)c(b−d)(d−r),acer(b−d)d(a−c)(c−e)(e−k),−bdk(a−c)(c−e)ce(b−d)(d−r),−ace(b−d)(d−r)d(a−c)(c−e)(e−k),bd(a−c)(c−e)(e−k)ce(b−d)(d−r),−c^2e(b−d)(d−r)d(a−c)(c−e)(e−k),d^2(a−c)(c−e)(e−k)ce(b−d)(d−r),−ce^2(b−d)(d−r)d(a−c)(c−e)(e−k),dr(a−c)(c−e)(e−k)ce(b−d)(d−r),−cek(b−d)(d−r)d(a−c)(c−e)(e−k),adr(c−e)(e−k)ce(b−d)(d−r),−bcek(d−r)d(a−c)(c−e)(e−k),adr(e−k)e(b−d)(d−r),−bcek(a−c)(c−e)(e−k),adr(b−d)(d−r),bek(c−e)(e−k),−ard−r,−bke−k,a,b,c,d,e,r,k,…

Proof

The results are true for m = 0. Assume that m > 0 and the results hold for m − 1. That is

u60m−66=a,u60m−65=b,u60m−64=c,u60m−63=d,u60m−62=e,u60m−61=r,u60m−60=k,u60m−59=ara−c,u60m−58=bkb−d,u60m−57=acr(a−c)(c−e),u60m−56=bdk(b−d)(d−r)u60m−55=acer(a−c)(c−e)(e−k),u60m−54=−bdk(a−c)c(b−d)(d−r),u60m−53=−acer(b−d)d(a−c)(c−e)(e−k),u60m−52=bdk(a−c)(c−e)ce(b−d)(d−r),u60m−51=ace(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−50=−bd(a−c)(c−e)(e−k)ce(b−d)(d−r),u60m−49=c2e(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−48=−d2(a−c)(c−e)(e−k)ce(b−d)(d−r),u60m−47=ce2(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−46=−dr(a−c)(c−e)(e−k)ce(b−d)(d−r),u60m−45=cek(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−44=−adr(c−e)(e−k)ce(b−d)(d−r),u60m−43=bcek(d−r)d(a−c)(c−e)(e−k),u60m−42=−adr(e−k)e(b−d)(d−r),u60m−41=bcek(a−c)(c−e)(e−k),u60m−40=−adr(b−d)(d−r),u60m−39=−bek(c−e)(e−k),u60m−38=ard−r,u60m−37=bke−k,u60m−36=−a,u60m−35=−b,u60m−34=−c,u60m−33=−d,u60m−32=−e,u60m−31=−r,u60m−30=−k,u60m−29=−ara−c,u60m−28=−bkb−d,u60m−27=−acr(a−c)(c−e),u60m−26=−bdk(b−d)(d−r),u60m−25=−acer(a−c)(c−e)(e−k),u60m−24=bdk(a−c)c(b−d)(d−r),u60m−23=acer(b−d)d(a−c)(c−e)(e−k),u60m−22=−bdk(a−c)(c−e)ce(b−d)(d−r),u60m−21=−ace(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−20=bd(a−c)(c−e)(e−k)ce(b−d)(d−r),u60m−19=−c2e(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−18=d2(a−c)(c−e)(e−k)ce(b−d)(d−r),u60m−17=−ce2(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−16=dr(a−c)(c−e)(e−k)ce(b−d)(d−r),u60m−15=−cek(b−d)(d−r)d(a−c)(c−e)(e−k),u60m−14=adr(c−e)(e−k)ce(b−d)(d−r),u60m−13=−bcek(d−r)d(a−c)(c−e)(e−k),u60m−12=adr(e−k)e(b−d)(d−r),u60m−11=−bcek(a−c)(c−e)(e−k),u60m−10=adr(b−d)(d−r),u60m−9=bek(c−e)(e−k),u60m−8=−ard−r,u60m−7=−bke−k,u60m−6=a,u60m−5=b,u60m−4=c,u60m−3=d,u60m−2=e,u60m−1=r,u60m=k.

Next, some of these relations are verified.

u60m−6=u60m−8−u60m−8u60m−11u60m−11−u60m−13=−ard−r−ard−rbcek(a−c)(c−e)(e−k)−bcek(a−c)(c−e)(e−k)+bcek(d−r)d(a−c)(c−e)(e−k)=−ard−r−add−r=a.

Similarly,

u60m−5=u60m−7−u60m−7u60m−10u60m−10−u60m−12=−bke−k−−bke−kadr(b−d)(d−r)adr(b−d)(d−r)−adr(e−k)e(b−d)(d−r)=−bke−k−bee−k=b.

Thus, the rest of the relations can be similarly proved.

4 Numerical results

In this section, we are mainly interested in confirming the obtained theoretical results by providing some numerical examples.

Example 1

In Figure 1, we confirm that the local stability of the equilibrium occurs if condition (6) is satisfied. Here, the constants and initial values have been taken by a = 0.1, b = 2, c = 0.8, d = 4, u₋₆ = 0.01, u₋₅ = −0.02, u₋₄ = 0.03, u₋₃ = −0.2, u₋₂ = 0.3, u₋₁ = −0.3, u₀ = 0.4.

Figure 1

Local Stability of the Equilibrium.

Example 2

The global stability of the equilibrium point of Eq. (1) is shown in Figure 2 under the values a = 0.3, b = 0.7, c = 4, d = 0.9, u₋₆ = 0.1, u₋₅ = 1, u₋₄ = −3, u₋₃ = 2, u₋₂ = −4, u₋₁ = 5, u₀ = −6.

Figure 2

Global Behaviour about the Equilibrium.

Example 3

Figure 3 illustrates the periodicity of the solution of Eq. (12) when we assume that u₋₆ = 2.8, u₋₅ = 3.5, u₋₄ = 1.6, u₋₃ = 2.5, u₋₂ = 1, u₋₁ = 3, u₀ = 2.

Figure 3

The Periodicity of Eq. (12).

Example 1

The stability of Eq. (12) is depicted in Figure 4 under the values a = 0.2, b = 0.1, c = 0.3, d = 0.5, u₋₆ = 0.2, u₋₅ = 2, u₋₄ = 1.8, u₋₃ = 2.5, u₋₂ = 0.5, u₋₁ = 1, u₀ = 0.2.

Figure 4

Stability of the Equilibrium.

5 Conclusion

This paper has introduced some new results for tow different equations. The local stability of the obtained equilibrium points are shown. Moreover, we obtained that every solution of the second proposed equation is periodic with period 60. The stability and the periodicity are depicted in some of the presented figures. The used technique in determining the exact solution of the first proposed equation can be straightforwardly applied for equations with high order.

References

[1] El-Dessoky M.M., Elabbasy E.M., Asiri A., Dynamics and solutions of a fifth-Order nonlinear difference equation, Disc. Dyn. Natu. Soci., 2018, 2018, 1–21.10.1155/2018/9129354Search in Google Scholar

[2] Almatrafi M.B., Alzubaidi M.M., Analysis of the qualitative behaviour of an eighth-order fractional difference equation, Open J. Disc. Appl. Math., 2019, 2(1), 41–47.10.30538/psrp-odam2019.0010Search in Google Scholar

[3] Kalabušić S., Nurkanović M., Nurkanović Z., Global dynamics of certain mix monotone difference equation, Mathematics, 2018, 6(10), 1–13.10.3390/math6010010Search in Google Scholar

[4] Aloqeili M., Dynamics of a rational difference equation, Appl. Math. Comput., 2006, 176(2), 768–774.10.1016/j.amc.2005.10.024Search in Google Scholar

[5] El-Owaidy H.M., Ahmed A.M., Youssef A.M., The dynamics of the recursive sequence xn+1=(αxn−1)/(β+yxn+1p), Appl. Math. Let., 2005, 18(9), 1013–1018.10.1016/j.aml.2003.09.014Search in Google Scholar

[6] Tollu D.T., Yazlik Y., Taskara N., On the solutions of two special types of riccati difference equation via Fibonacci numbers, Adv. Differ. Equ., 2013, 174, 1–7.10.1186/1687-1847-2013-174Search in Google Scholar

[7] Almatrafi M.B., Solutions structures for some systems of fractional difference equations, Open J. Math. Anal., 2019, 3(1), 51–61.10.30538/psrp-oma2019.0032Search in Google Scholar

[8] Almatrafi M.B., Elsayed E.M., Alzahrani F., Qualitative behavior of two rational difference equations, Fundam. J. Math. Appl., 2018, 1(2), 194–204.10.33401/fujma.454999Search in Google Scholar

[9] Almatrafi M.B., Elsayed E.M., Solutions and formulae for some systems of difference equations, MathLAB J., 2018, 1(3), 356–369.Search in Google Scholar

[10] Almatrafi M.B., Elsayed E.M., Alzahrani F., Qualitative behavior of a quadratic second-order rational difference equation, Int. J. Adv. Math., 2019, 2019(1), 1–14.Search in Google Scholar

[11] Amleh A.M., Grove E.A., Georgiou D.A., Ladas G., On the recursive sequence x_n+1 = α + x_n−1/x_n, J. Math. Anal. Appl., 1999, 233, 790–798.10.1006/jmaa.1999.6346Search in Google Scholar

[12] Belhannache F., Touafek N., Abo-zeid R., On a higher-order rational difference equation, J. Appl. Math. & Informatics, 2016, 34(5–6), 369–382.10.14317/jami.2016.369Search in Google Scholar

[13] Elabbasy E.M., El-Metawally H., Elsayed E.M., On the difference equation xn+1=(axn2+bxn−1xn−k)/(cxn2+dxn−1xn−k), Sarajevo J. Math., 2008, 4(17), 1–10.Search in Google Scholar

[14] El-Moneam M.A., Zayed E.M.E., Dynamics of the rational difference equation, Inf. Sci. Let., 2014, 3(2), 45–53.10.12785/isl/030202Search in Google Scholar

[15] Elsayed E.M., Dynamics of recursive sequence of order two, Kyungpook Math. J., 2010, 50, 483–497.10.5666/KMJ.2010.50.4.483Search in Google Scholar

[16] Garić-Demirović M., Nurkanović M., Nurkanović Z., Stability, periodicity and Neimark-Sacker bifurcation of certain homogeneous fractional difference equations, Int. J. Differ. Equ., 2017, 12(1), 27–53.Search in Google Scholar

[17] Ghazel M., Elsayed E.M., Matouk A.E., Mousallam A.M., Investigating dynamical behaviors of the difference equation x_n+1 = Cx_n−5/(A + Bx_n−2x_n−5), J. Nonlinear Sci. Appl., 2017, 10, 4662–4679.10.22436/jnsa.010.09.09Search in Google Scholar

[18] Hamza A.E., On the recursive sequence x_n+1 = α + x_n−1/x_n, J. Math. Anal. Appl., 2006, 322, 668–674.10.1016/j.jmaa.2005.09.029Search in Google Scholar

[19] Khaliq A., Hassan Sk.S., Dynamics of a rational difference equation x_n+1 = ax_n+(α+βx_n−k)/(A + Bx_n−k), Int. J. Adv. Math., 2018, 1(2018), 159–179.Search in Google Scholar

[20] Khyat T., Kulenović M.R.S., The invariant curve caused by Neimark-Sacker bifurcation of a perturbed Beverton-Holt difference equation, Int. J. Differ. Equ., 2017, 12(2), 267–280.Search in Google Scholar

[21] Kocic V.L., Ladas G., Global behaviour of nonlinear difference equations of higher order with applications, Kluwer Academic Publishers, Dordrecht, 1993.10.1007/978-94-017-1703-8Search in Google Scholar

[22] Kostrov Y., Kudlak Z., On a second-order rational difference equation with a quadratic term, Int. J. differ. Equ., 2016, 11(2), 179–202.Search in Google Scholar

[23] Liu K., Li P., Han F., Zhong W., Global dynamics of nonlinear difference equation x_n+1 = A + x_n/x_n−1x_n−2, J. Comput. Anal. Appl., 2018, 24(6), 1125–1132.Search in Google Scholar

[24] Moranjkić S., Nurkanović Z., Local and global dynamics of certain second-order rational difference equations containing quadratic terms, Adv. Dyn. Sys. Appl., 2017, 12(2), 123–157.Search in Google Scholar

[25] Moaaz O., Chalishajar D., Bazighifan O., Some qualitative behavior of solutions of general class of difference equations, Mathematics, 2019, 7, 1–12.10.3390/math7070585Search in Google Scholar

[26] Saleh M., Farhat A., Global asymptotic stability of the higher order equation xn+1=axn+bxn−kA+Bxn−k, J. Appl. Math. Comput., 2017, 55(1), 135–148.10.1007/s12190-016-1029-4Search in Google Scholar

[27] Saleh M., Alkoumi N., Farhat A., On the dynamic of a rational difference equation x_n+1 = α+βx_n + yx_n−k/Bx_n + Cx_n−k, Chaos, Solitons and Fractals, 2017, 96, 76–84.10.1016/j.chaos.2016.12.026Search in Google Scholar

[28] Simsek D., Cinar C., Yalcinkaya I., On the recursive sequence xn+1=xn−31+xn−1, Int. J. Contemp. Math. Sci., 2006, 1(10), 475–480.10.12988/ijcms.2006.06052Search in Google Scholar

[29] Kulenovic M.R.S., Ladas G., Dynamics of second order rational difference equations: with open problems and conjectures, Chapman & Hall/CRC Press, Boca Raton, 2001.10.1201/9781420035384Search in Google Scholar

Received: 2018-09-04

Accepted: 2019-12-18

Published Online: 2020-05-13

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

https://doi.org/10.1515/nleng-2020-0014

Keywords for this article

attractivity; boundedness; difference equations; equilibrium; periodicity; stability

Creative Commons

BY 4.0