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
Almatrafi and Alzubaidi [2] analysed the asymptotic attractivity, periodic nature and boundedness of the eighth order difference equation
Moreover, Kalabušić et al. [3] studied the global character of the second order equation
The author in [4] highlighted the semi-cycle behaviour, stability character and solutions for the second order difference equation
El-Owaidy et al. [5] explored the global stability of the difference equation
In [6], the author utilized Fibonacci sequence to construct the solutions of special type of Riccati difference equations given by
Almatrafi [7] presented forms of exact solutions to the fractional system
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:
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
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:
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, ∞)3 ⟶ (0, ∞) by
Hence,
Evaluating Eqs. (3), (4) and Eq. (5) at ū gives us
Hence, the linearised equation of Eq. (1) around ū is shown as follows:
Theorem 1
Assume that
Then, the equilibrium point of Eq.(1) is locally asymptotically stable.
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
Let
Proof
Let r1, r2 ∈ ℝ, assume that g : [r1, r2]3 ⟶ [r1, r2] is a function defined by Eq. (2), and assume that
Or,
Simplifying this gives
Subtracting Eq.(10) from Eq.(9) yields
Hence, if a ≠ 1, then ζ = η. Consequently, Theorem C [29] concludes that every solution of Eq. (1) converges to ū.
Theorem 3
Assume that
Proof
The proof is similar to the previous proof. Thus, it is omitted.
2.3 Exact solution of u m + 1 = u m − 1 + u m − 1 u m − 4 u m − 4 − u m − 6
We now begin with presenting the exact solution of the equation
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
Here,
Proof
The relations are true for m = 0. Now, assume that m > 0 and that our assumption holds for m − 1. That is,
Two random relations will be now selected and proved. It can be easily observed from Eq. (11) that
We can also see from Eq. (11) that
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
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
Proof
The results are true for m = 0. Assume that m > 0 and the results hold for m − 1. That is
Next, some of these relations are verified.
Similarly,
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.

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.

Global Behaviour about 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
[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
[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+(α+βxn−k)/(A + Bxn−k), 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
[27] Saleh M., Alkoumi N., Farhat A., On the dynamic of a rational difference equation xn+1 = α+βxn + yxn−k/Bxn + Cxn−k, 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
[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
© 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
- Comparison of the method of variation of parameters to semi-analytical methods for solving nonlinear boundary value problems in engineering
- Nonlinear H-infinity control for switched reluctance machines
- Energy flow of a 2018 FIA F1 racing car and proposed changes to the powertrain rules
- Risk index to monitor an anaerobic digester using a dynamic model based on dilution rate, temperature, and pH
- MHD Peristaltic flow of a nanofluid in a constricted artery for different shapes of nanosized particles
- Comparative study of homotopy perturbation transformation with homotopy perturbation Elzaki transform method for solving nonlinear fractional PDE
- Approximate method for solving strongly fractional nonlinear problems using fuzzy transform
- Numerical approach to MHD flow of power-law fluid on a stretching sheet with non-uniform heat source
- Entropy generation in an inclined porous channel with suction/injection
- Heat transfer from convecting-radiating fin through optimized Chebyshev polynomials with interior point algorithm
- Two dimensional simulation of laminar flow by three- jet in a semi-confined space
- Influence of temperature-dependent properties on a gravity-driven thin film along inclined plate
- Simulation and time-frequency analysis of the longitudinal train dynamics coupled with a nonlinear friction draft gear
- Study of differential transform technique for transient hydromagnetic Jeffrey fluid flow from a stretching sheet
- Generalized second-order slip for unsteady convective flow of a nanofluid: a utilization of Buongiorno’s two-component nonhomogeneous equilibrium model
- Numerical treatment for the solution of singularly perturbed pseudo-parabolic problem on an equidistributed grid
- Relative sea-level rise and land subsidence in Oceania from tide gauge and satellite GPS
- On finite series solutions of conformable time-fractional Cahn-Allen equation
- 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
- MHD natural convection flow of Casson fluid in an annular microchannel containing porous medium with heat generation/absorption
- Numerical simulation of variable thermal conductivity on 3D flow of nanofluid over a stretching sheet
- Two meshless methods for solving nonlinear ordinary differential equations in engineering and applied sciences
- Thermoelastic analysis of FGM hollow cylinder for variable parameters and temperature distributions using FEM
- Qualitative analysis for two fractional difference equations
- MHD fractionalized Jeffrey fluid over an accelerated slipping porous plate
- Nonlinear analysis of high accuracy and reliability in traffic flow prediction
- Numerical solution of time-dependent Emden-Fowler equations using bivariate spectral collocation method on overlapping grids
- A reliable analytical technique for fractional Caudrey-Dodd-Gibbon equation with Mittag-Leffler kernel
- Accelerated HPSTM: An efficient semi-analytical technique for the solution of nonlinear PDE’s
- Effect of magnetized variable thermal conductivity on flow and heat transfer characteristics of unsteady Williamson fluid
- Couple stress fluid flow due to slow steady oscillations of a permeable sphere
- State-of-the-art of MW-level capacity oceanic current turbines
- Approximate solution for fractional attractor one-dimensional Keller-Segel equations using homotopy perturbation sumudu transform method
- Nonlinear absolute sea-level patterns in the long-term-trend tide gauges of the West Coast of North America
- Insight into the dynamics of non-Newtonian Casson fluid over a rotating non-uniform surface subject to Coriolis force
- Mixed convection flow in a vertical channel in the presence of wall conduction, variable thermal conductivity and viscosity
- A new structure formulations for cubic B-spline collocation method in three and four-dimensions
- Mathematical and numerical optimality of non-singular fractional approaches on free and forced linear oscillator
- MHD mixed convection on an inclined stretching plate in Darcy porous medium with Soret effect and variable surface conditions
- Comparative study of two techniques on some nonlinear problems based ussing conformable derivative
Artikel in diesem Heft
- Comparison of the method of variation of parameters to semi-analytical methods for solving nonlinear boundary value problems in engineering
- Nonlinear H-infinity control for switched reluctance machines
- Energy flow of a 2018 FIA F1 racing car and proposed changes to the powertrain rules
- Risk index to monitor an anaerobic digester using a dynamic model based on dilution rate, temperature, and pH
- MHD Peristaltic flow of a nanofluid in a constricted artery for different shapes of nanosized particles
- Comparative study of homotopy perturbation transformation with homotopy perturbation Elzaki transform method for solving nonlinear fractional PDE
- Approximate method for solving strongly fractional nonlinear problems using fuzzy transform
- Numerical approach to MHD flow of power-law fluid on a stretching sheet with non-uniform heat source
- Entropy generation in an inclined porous channel with suction/injection
- Heat transfer from convecting-radiating fin through optimized Chebyshev polynomials with interior point algorithm
- Two dimensional simulation of laminar flow by three- jet in a semi-confined space
- Influence of temperature-dependent properties on a gravity-driven thin film along inclined plate
- Simulation and time-frequency analysis of the longitudinal train dynamics coupled with a nonlinear friction draft gear
- Study of differential transform technique for transient hydromagnetic Jeffrey fluid flow from a stretching sheet
- Generalized second-order slip for unsteady convective flow of a nanofluid: a utilization of Buongiorno’s two-component nonhomogeneous equilibrium model
- Numerical treatment for the solution of singularly perturbed pseudo-parabolic problem on an equidistributed grid
- Relative sea-level rise and land subsidence in Oceania from tide gauge and satellite GPS
- On finite series solutions of conformable time-fractional Cahn-Allen equation
- 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
- MHD natural convection flow of Casson fluid in an annular microchannel containing porous medium with heat generation/absorption
- Numerical simulation of variable thermal conductivity on 3D flow of nanofluid over a stretching sheet
- Two meshless methods for solving nonlinear ordinary differential equations in engineering and applied sciences
- Thermoelastic analysis of FGM hollow cylinder for variable parameters and temperature distributions using FEM
- Qualitative analysis for two fractional difference equations
- MHD fractionalized Jeffrey fluid over an accelerated slipping porous plate
- Nonlinear analysis of high accuracy and reliability in traffic flow prediction
- Numerical solution of time-dependent Emden-Fowler equations using bivariate spectral collocation method on overlapping grids
- A reliable analytical technique for fractional Caudrey-Dodd-Gibbon equation with Mittag-Leffler kernel
- Accelerated HPSTM: An efficient semi-analytical technique for the solution of nonlinear PDE’s
- Effect of magnetized variable thermal conductivity on flow and heat transfer characteristics of unsteady Williamson fluid
- Couple stress fluid flow due to slow steady oscillations of a permeable sphere
- State-of-the-art of MW-level capacity oceanic current turbines
- Approximate solution for fractional attractor one-dimensional Keller-Segel equations using homotopy perturbation sumudu transform method
- Nonlinear absolute sea-level patterns in the long-term-trend tide gauges of the West Coast of North America
- Insight into the dynamics of non-Newtonian Casson fluid over a rotating non-uniform surface subject to Coriolis force
- Mixed convection flow in a vertical channel in the presence of wall conduction, variable thermal conductivity and viscosity
- A new structure formulations for cubic B-spline collocation method in three and four-dimensions
- Mathematical and numerical optimality of non-singular fractional approaches on free and forced linear oscillator
- MHD mixed convection on an inclined stretching plate in Darcy porous medium with Soret effect and variable surface conditions
- Comparative study of two techniques on some nonlinear problems based ussing conformable derivative