Startseite Qualitative analysis for two fractional difference equations
Artikel Open Access

Qualitative analysis for two fractional difference equations

  • Mohammed B. Almatrafi EMAIL logo und Marwa M. Alzubaidi
Veröffentlicht/Copyright: 13. Mai 2020
Veröffentlichen auch Sie bei De Gruyter Brill

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.

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+βynyn3Ayn4+Byn3.

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

xn+1=c1xn3+c2xn3c3xn3c4xn7.

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

xn+1=xnxn1+αxn+βxn1axnxn1+bxn1.

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

xn+1=xn1axnxn1.

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

xn+1=αxn1β+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=1yn1.

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

xn+1=xn1yn3yn11xn1yn3,yn+1=yn1xn3xn1±1±yn1xn3,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=aum1+bum1um4cum4dum6,m=0,1,,um+1=aum1bum1um4cum4dum6,m=0,1,,

where the coefficients a, b, c and d are supposed to be positive real numbers and the initial data uifor 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=aum1+bum1um4cum4dum6,m=0,1,,(1)

where the constants a, b, c and d are assumed to be positive real numbers and the initial conditions uifor 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)(cd). In order to linearise about the fixed point, we define a function g : (0, ∞)3 ⟶ (0, ∞) by

g(x,y,z)=ax+bxycydz.(2)

Hence,

g(x,y,z)x=a+by(cydz),(3)
g(x,y,z)y=bdxz(cydz)2,(4)
g(x,y,z)z=bdxy(cydz)2.(5)

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

g(u¯,u¯,u¯)x=a+bu¯(cu¯du¯)=a+bcd:=q1,g(u¯,u¯,u¯)y=bdu¯2(cu¯du¯)2=bd(cd)2:=q2,g(u¯,u¯,u¯)z=bdu¯2(cu¯du¯)2=bd(cd)2:=q3.

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

Un+1+q1Un1+q2Un4+q3Un6=0.

Theorem 1

Assume that

a+bcd<12bd(cd)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+bcd)+|bd(cd)2|+(bd(cd)2)<1.(8)

Inequality (8) can be simply reduced to

a+bcd<12bd(cd)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+bcd>0.Then, the equilibrium point of Eq.(1) is a global attractor if a ≠ 1.

Proof

Let r1, r2 ∈ ℝ, assume that g : [r1, r2]3 ⟶ [r1, r2] is a function defined by Eq. (2), and assume that a+bcd>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+bcd<0.Then, the equilibrium point of Eq.(1) is a global attractor ifdac + b.

Proof

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

2.3 Exact solution of um+1=um1+um1um4um4um6

We now begin with presenting the exact solution of the equation

um+1=um1+um1um4um4um6.(11)

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

Theorem 4

Let{um}m=6be a solution to Eq. (11)and suppose thatu−6 = a, u−5 = b, u−4 = c, u−3 = d, u−2 = e, u−1 = r, u0 = k. Then, form = 0, 1, …, we have

u10m6=(F2m+1d+F2m1b)(F2m+1r+F2m1d)(F2mk+F2m2e)(F2mc+F2m2a)(F2me+F2m2c)ec(db)(rd),u10m5=(F2m+1k+F2m1e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)(F2md+F2m2b)(F2mr+F2m2d)d(ca)(ec)(ke),u10m4=(F2m+1d+F2m1b)(F2m+1r+F2m1d)(F2mk+F2m2e)(F2me+F2m2c)(F2m+2c+F2ma)ec(db)(rd),u10m3=(F2m+2d+F2mb)(F2mr+F2m2d)(F2m+1k+F2m1e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)d(ca)(ec)(ke),u10m2=(F2m+2e+F2mc)(F2m+2c+F2ma)(F2m+1d+F2m1b)(F2m+1r+F2m1d)(F2mk+F2m2e)ec(db)(rd),u10m1=(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m1e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)d(ca)(ec)(ke),u10m=(F2m+2k+F2me)(F2m+2e+F2mc)(F2m+2c+F2ma)(F2m+1d+F2m1b)(F2m+1r+F2m1d)ec(db)(rd),u10m+1=(F2m+3c+F2m+1a)(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m1e)(F2m+1e+F2m1c)d(ca)(ec)(ke),u10m+2=(F2m+3d+F2m+1b)(F2m+1r+F2m1d)(F2m+2k+F2me)(F2m+2e+F2mc)(F2m+2c+F2ma)ec(db)(rd),u10m+3=(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m1e)(F2m+3e+F2m+1c)(F2m+3c+F2m+1a)d(ca)(ec)(ke).

Here, {Fm}m=0 = {1, 1, 2, 3, 5, 8, … }, i.e, Fm = Fm−1 + Fm−2, F−1 = 0 and F−2 = −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,

u10m16=(F2m1d+F2m3b)(F2m1r+F2m3d)(F2m2k+F2m4e)(F2m2c+F2m4a)(F2m2e+F2m4c)ec(db)(rd),u10m15=(F2m1k+F2m3e)(F2m1e+F2m3c)(F2m1c+F2m3a)(F2m2d+F2m4b)(F2m2r+F2m4d)d(ca)(ec)(ke),u10m14=(F2m1d+F2m3b)(F2m1r+F2m3d)(F2m2k+F2m4e)(F2m2e+F2m4c)(F2m2c+F2m4a)ec(db)(rd),u10m13=(F2md+F2m2b)(F2m2r+F2m4d)(F2m1k+F2m3e)(F2m1e+F2m3c)(F2m1c+F2m3a)d(ca)(ec)(ke),u10m12=(F2me+F2m2c)(F2mc+F2m2a)(F2m1d+F2m3b)(F2m1r+F2m3d)(F2m2k+F2m4e)ec(db)(rd),u10m11=(F2md+F2m2b)(F2mr+F2m2d)(F2m1k+F2m3e)(F2m1e+F2m3c)(F2m1c+F2m3a)d(ca)(ec)(ke),u10m10=(F2mk+F2m2e)(F2me+F2m2c)(F2mc+F2m2a)(F2m1d+F2m3b)(F2m1r+F2m3d)ec(db)(rd),u10m9=(F2m+1c+F2m1a)(F2md+F2m2b)(F2mr+F2m2d)(F2m1k+F2m3e)(F2m1e+F2m3c)d(ca)(ec)(ke),u10m8=(F2m+1d+F2m1b)(F2m1r+F2m3d)(F2mk+F2m2e)(F2me+F2m2c)(F2mc+F2m2a)ec(db)(rd),u10m7=(F2md+F2m2b)(F2mr+F2m2d)(F2m1k+F2m3e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)d(ca)(ec)(ke).

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

u10m5=u10m7+y10m7y10m10y10m10y10m12=u10m71+y10m10y10m10y10m12=u10m71+F2mkF2m2eF2mkF1m2eF2m2k+F2m4e=u10m71+F2mkF2m2eF2m1kF2m3e=u10m7F2m+1kF2m1eF2m1kF2m3e=(F2md+F2m2b)(F2mr+F2m2d)(F2m1k+F2m3e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)d(ca)(ec)(ke)F2m+1kF2m1eF2m1kF2m3e=(F2m+1k+F2m1e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)(F2md+F2m2b)(F2mr+F2m2d)d(ca)(ec)(ke).

We can also see from Eq. (11) that

u10m1=u10m3+u10m3u10m6u10m6u10m8=u10m31+u10m6u10m6u10m8=u10m31+F2m+1rF2m1dF2m+1rF1m1dF2m1r+F2m3d=u10m31+F2m+1rF2m1dF2mrF2m2d=u10m3F2m+2rF2mdF2mrF2m2d=(F2m+2d+F2mb)(F2mr+F2m2d)(F2m+1k+F2m1e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)d(ca)(ec)(ke)F2m+2rF2mdF2mrF2m2d=(F2m+2d+F2mb)(F2m+2r+F2md)(F2m+1k+F2m1e)(F2m+1e+F2m1c)(F2m+1c+F2m1a)d(ca)(ec)(ke).

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=um1um1um4um4um6,m=0,1,,(12)

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

Theorem 5

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

a,b,c,d,e,r,k,arac,bkbd,acr(ac)(ce),bdk(bd)(dr),acer(ac)(ce)(ek),bdk(ac)c(bd)(dr),acer(bd)d(ac)(ce)(ek),bdk(ac)(ce)ce(bd)(dr),ace(bd)(dr)d(ac)(ce)(ek),bd(ac)(ce)(ek)ce(bd)(dr),c^2e(bd)(dr)d(ac)(ce)(ek),d^2(ac)(ce)(ek)ce(bd)(dr),ce^2(bd)(dr)d(ac)(ce)(ek),dr(ac)(ce)(ek)ce(bd)(dr),cek(bd)(dr)d(ac)(ce)(ek),adr(ce)(ek)ce(bd)(dr),bcek(dr)d(ac)(ce)(ek),adr(ek)e(bd)(dr),bcek(ac)(ce)(ek),adr(bd)(dr),bek(ce)(ek),ardr,bkek,a,b,c,d,e,r,k,arac,bkbd,acr(ac)(ce),bdk(bd)(dr),acer(ac)(ce)(ek),bdk(ac)c(bd)(dr),acer(bd)d(ac)(ce)(ek),bdk(ac)(ce)ce(bd)(dr),ace(bd)(dr)d(ac)(ce)(ek),bd(ac)(ce)(ek)ce(bd)(dr),c^2e(bd)(dr)d(ac)(ce)(ek),d^2(ac)(ce)(ek)ce(bd)(dr),ce^2(bd)(dr)d(ac)(ce)(ek),dr(ac)(ce)(ek)ce(bd)(dr),cek(bd)(dr)d(ac)(ce)(ek),adr(ce)(ek)ce(bd)(dr),bcek(dr)d(ac)(ce)(ek),adr(ek)e(bd)(dr),bcek(ac)(ce)(ek),adr(bd)(dr),bek(ce)(ek),ardr,bkek,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

u60m66=a,u60m65=b,u60m64=c,u60m63=d,u60m62=e,u60m61=r,u60m60=k,u60m59=arac,u60m58=bkbd,u60m57=acr(ac)(ce),u60m56=bdk(bd)(dr)u60m55=acer(ac)(ce)(ek),u60m54=bdk(ac)c(bd)(dr),u60m53=acer(bd)d(ac)(ce)(ek),u60m52=bdk(ac)(ce)ce(bd)(dr),u60m51=ace(bd)(dr)d(ac)(ce)(ek),u60m50=bd(ac)(ce)(ek)ce(bd)(dr),u60m49=c2e(bd)(dr)d(ac)(ce)(ek),u60m48=d2(ac)(ce)(ek)ce(bd)(dr),u60m47=ce2(bd)(dr)d(ac)(ce)(ek),u60m46=dr(ac)(ce)(ek)ce(bd)(dr),u60m45=cek(bd)(dr)d(ac)(ce)(ek),u60m44=adr(ce)(ek)ce(bd)(dr),u60m43=bcek(dr)d(ac)(ce)(ek),u60m42=adr(ek)e(bd)(dr),u60m41=bcek(ac)(ce)(ek),u60m40=adr(bd)(dr),u60m39=bek(ce)(ek),u60m38=ardr,u60m37=bkek,u60m36=a,u60m35=b,u60m34=c,u60m33=d,u60m32=e,u60m31=r,u60m30=k,u60m29=arac,u60m28=bkbd,u60m27=acr(ac)(ce),u60m26=bdk(bd)(dr),u60m25=acer(ac)(ce)(ek),u60m24=bdk(ac)c(bd)(dr),u60m23=acer(bd)d(ac)(ce)(ek),u60m22=bdk(ac)(ce)ce(bd)(dr),u60m21=ace(bd)(dr)d(ac)(ce)(ek),u60m20=bd(ac)(ce)(ek)ce(bd)(dr),u60m19=c2e(bd)(dr)d(ac)(ce)(ek),u60m18=d2(ac)(ce)(ek)ce(bd)(dr),u60m17=ce2(bd)(dr)d(ac)(ce)(ek),u60m16=dr(ac)(ce)(ek)ce(bd)(dr),u60m15=cek(bd)(dr)d(ac)(ce)(ek),u60m14=adr(ce)(ek)ce(bd)(dr),u60m13=bcek(dr)d(ac)(ce)(ek),u60m12=adr(ek)e(bd)(dr),u60m11=bcek(ac)(ce)(ek),u60m10=adr(bd)(dr),u60m9=bek(ce)(ek),u60m8=ardr,u60m7=bkek,u60m6=a,u60m5=b,u60m4=c,u60m3=d,u60m2=e,u60m1=r,u60m=k.

Next, some of these relations are verified.

u60m6=u60m8u60m8u60m11u60m11u60m13=ardrardrbcek(ac)(ce)(ek)bcek(ac)(ce)(ek)+bcek(dr)d(ac)(ce)(ek)=ardraddr=a.

Similarly,

u60m5=u60m7u60m7u60m10u60m10u60m12=bkekbkekadr(bd)(dr)adr(bd)(dr)adr(ek)e(bd)(dr)=bkekbeek=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−6 = 0.01, u−5 = −0.02, u−4 = 0.03, u−3 = −0.2, u−2 = 0.3, u−1 = −0.3, u0 = 0.4.

Figure 1 Local Stability of the Equilibrium.
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−6 = 0.1, u−5 = 1, u−4 = −3, u−3 = 2, u−2 = −4, u−1 = 5, u0 = −6.

Figure 2 Global Behaviour about the Equilibrium.
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−6 = 2.8, u−5 = 3.5, u−4 = 1.6, u−3 = 2.5, u−2 = 1, u−1 = 3, u0 = 2.

Figure 3 The Periodicity of Eq. (12).
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−6 = 0.2, u−5 = 2, u−4 = 1.8, u−3 = 2.5, u−2 = 0.5, u−1 = 1, u0 = 0.2.

Figure 4 Stability of the Equilibrium.
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/9129354Suche 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.0010Suche 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/math6010010Suche 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.024Suche in Google Scholar

[5] El-Owaidy H.M., Ahmed A.M., Youssef A.M., The dynamics of the recursive sequence xn+1=(αxn1)/(β+yxn+1p), Appl. Math. Let., 2005, 18(9), 1013–1018.10.1016/j.aml.2003.09.014Suche 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-174Suche 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.0032Suche 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.454999Suche 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.Suche 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.Suche in Google Scholar

[11] Amleh A.M., Grove E.A., Georgiou D.A., Ladas G., On the recursive sequence xn+1 = α + xn−1/xn, J. Math. Anal. Appl., 1999, 233, 790–798.10.1006/jmaa.1999.6346Suche 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.369Suche in Google Scholar

[13] Elabbasy E.M., El-Metawally H., Elsayed E.M., On the difference equation xn+1=(axn2+bxn1xnk)/(cxn2+dxn1xnk), Sarajevo J. Math., 2008, 4(17), 1–10.Suche 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/030202Suche 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.483Suche 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.Suche in Google Scholar

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

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

[19] Khaliq A., Hassan Sk.S., Dynamics of a rational difference equation xn+1 = axn+(α+βxnk)/(A + Bxnk), Int. J. Adv. Math., 2018, 1(2018), 159–179.Suche 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.Suche 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-8Suche 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.Suche in Google Scholar

[23] Liu K., Li P., Han F., Zhong W., Global dynamics of nonlinear difference equation xn+1 = A + xn/xn−1xn−2, J. Comput. Anal. Appl., 2018, 24(6), 1125–1132.Suche 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.Suche 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/math7070585Suche in Google Scholar

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

[27] Saleh M., Alkoumi N., Farhat A., On the dynamic of a rational difference equation xn+1 = α+βxn + yxnk/Bxn + Cxnk, Chaos, Solitons and Fractals, 2017, 96, 76–84.10.1016/j.chaos.2016.12.026Suche in Google Scholar

[28] Simsek D., Cinar C., Yalcinkaya I., On the recursive sequence xn+1=xn31+xn1, Int. J. Contemp. Math. Sci., 2006, 1(10), 475–480.10.12988/ijcms.2006.06052Suche 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/9781420035384Suche in Google Scholar

Received: 2018-09-04
Accepted: 2019-12-18
Published Online: 2020-05-13

© 2020 M. B. Almatrafi and Marwa M. Alzubaidi, published by De Gruyter

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

Artikel in diesem Heft

  1. Comparison of the method of variation of parameters to semi-analytical methods for solving nonlinear boundary value problems in engineering
  2. Nonlinear H-infinity control for switched reluctance machines
  3. Energy flow of a 2018 FIA F1 racing car and proposed changes to the powertrain rules
  4. Risk index to monitor an anaerobic digester using a dynamic model based on dilution rate, temperature, and pH
  5. MHD Peristaltic flow of a nanofluid in a constricted artery for different shapes of nanosized particles
  6. Comparative study of homotopy perturbation transformation with homotopy perturbation Elzaki transform method for solving nonlinear fractional PDE
  7. Approximate method for solving strongly fractional nonlinear problems using fuzzy transform
  8. Numerical approach to MHD flow of power-law fluid on a stretching sheet with non-uniform heat source
  9. Entropy generation in an inclined porous channel with suction/injection
  10. Heat transfer from convecting-radiating fin through optimized Chebyshev polynomials with interior point algorithm
  11. Two dimensional simulation of laminar flow by three- jet in a semi-confined space
  12. Influence of temperature-dependent properties on a gravity-driven thin film along inclined plate
  13. Simulation and time-frequency analysis of the longitudinal train dynamics coupled with a nonlinear friction draft gear
  14. Study of differential transform technique for transient hydromagnetic Jeffrey fluid flow from a stretching sheet
  15. Generalized second-order slip for unsteady convective flow of a nanofluid: a utilization of Buongiorno’s two-component nonhomogeneous equilibrium model
  16. Numerical treatment for the solution of singularly perturbed pseudo-parabolic problem on an equidistributed grid
  17. Relative sea-level rise and land subsidence in Oceania from tide gauge and satellite GPS
  18. On finite series solutions of conformable time-fractional Cahn-Allen equation
  19. A generalized perspective of Fourier and Fick’s laws: Magnetized effects of Cattaneo-Christov models on transient nanofluid flow between two parallel plates with Brownian motion and thermophoresis
  20. MHD natural convection flow of Casson fluid in an annular microchannel containing porous medium with heat generation/absorption
  21. Numerical simulation of variable thermal conductivity on 3D flow of nanofluid over a stretching sheet
  22. Two meshless methods for solving nonlinear ordinary differential equations in engineering and applied sciences
  23. Thermoelastic analysis of FGM hollow cylinder for variable parameters and temperature distributions using FEM
  24. Qualitative analysis for two fractional difference equations
  25. MHD fractionalized Jeffrey fluid over an accelerated slipping porous plate
  26. Nonlinear analysis of high accuracy and reliability in traffic flow prediction
  27. Numerical solution of time-dependent Emden-Fowler equations using bivariate spectral collocation method on overlapping grids
  28. A reliable analytical technique for fractional Caudrey-Dodd-Gibbon equation with Mittag-Leffler kernel
  29. Accelerated HPSTM: An efficient semi-analytical technique for the solution of nonlinear PDE’s
  30. Effect of magnetized variable thermal conductivity on flow and heat transfer characteristics of unsteady Williamson fluid
  31. Couple stress fluid flow due to slow steady oscillations of a permeable sphere
  32. State-of-the-art of MW-level capacity oceanic current turbines
  33. Approximate solution for fractional attractor one-dimensional Keller-Segel equations using homotopy perturbation sumudu transform method
  34. Nonlinear absolute sea-level patterns in the long-term-trend tide gauges of the West Coast of North America
  35. Insight into the dynamics of non-Newtonian Casson fluid over a rotating non-uniform surface subject to Coriolis force
  36. Mixed convection flow in a vertical channel in the presence of wall conduction, variable thermal conductivity and viscosity
  37. A new structure formulations for cubic B-spline collocation method in three and four-dimensions
  38. Mathematical and numerical optimality of non-singular fractional approaches on free and forced linear oscillator
  39. MHD mixed convection on an inclined stretching plate in Darcy porous medium with Soret effect and variable surface conditions
  40. Comparative study of two techniques on some nonlinear problems based ussing conformable derivative
Heruntergeladen am 25.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/nleng-2020-0014/html
Button zum nach oben scrollen