Global Analysis of Solutions of a New Class of Rational Difference Equation
Abstract: The study suggests asymptotic behavior of the solution to a new class of difference equations: . where a, bi, α and β are positive real numbers for i = 0, 1, · · · , k , and the initial conditions ψ-j, ψ-j+1, · · ·, ψ0 are randomly positive real numbers where j = 2k + 1. Accordingly, we consider the stability, boundedness and periodicity of the solutions of this recursive sequence. Indeed, we give some interesting counter examples in order to verify our strong results.

1. Introduction

The global asymptotic behavior of the solutions and oscillation of solution are two such qualitative properties which are very important for applications in many areas such as control theory, mathematical biology, neural networks, etc. It is impossible to use computer based (numerical) techniques to study the oscillation or the asymptotic behavior of all solutions of a given equation due to the global nature of these properties. Therefore, these properties have received the attention of several mathematicians and engineers.

Currently, much attention has given to study the properties of the solutions of the recursive sequences from scientists in various disciplines. Specifically, the topics dealt with include the following:

- Finding equilibrium points for the recursive sequences;

- Investigating the local stability of the solutions of the recursive sequences;

- Finding conditions which insure that the solutions of the recursive sequences are bounded;

- Investigating the global asymptotic stability of the solutions of the recursive sequences;

- Finding conditions which insure that the solutions of equation are periodic with positive prime period two or more;

- Finding conditions for oscillation of solutions.

Closely related global convergence results were well-gained from these articles  -  . Khuong in  studied the dynamics the recursive sequences

${\psi }_{\eta +1}=a+{\left(\frac{{\psi }_{\eta -k}}{{\psi }_{\eta -r}}\right)}^{p}.$

For further related and special cases of this difference equations see    ,    .

Elsayed  studied the periodicity, the boundedness of the positive solution of the recursive sequences

${\psi }_{\eta +1}=a+\frac{b{\psi }_{\eta -l}+c{\psi }_{\eta -k}}{c{\psi }_{\eta -l}+d{\psi }_{\eta -k}}.$

Abdelrahman  considered analytical investigation of the solution of the recursive sequence

${\psi }_{\eta +1}=a{\psi }_{\eta -k}+\delta {\left(\frac{b{\psi }_{\eta -l}}{c{\psi }_{\eta -k}+d{\psi }_{\eta -l}}\right)}^{\alpha }.$

By new method, Elsayed  investigated the periodic solution of the equation

${\psi }_{\eta +1}=\alpha +\beta \frac{{\psi }_{\eta }}{{\psi }_{\eta -1}}+\gamma \frac{{\psi }_{\eta -1}}{{\psi }_{\eta }}.$

Also, Moaaz  completed the results of  .

In this work, we deal with some qualitative behaviour of the solutions of the recursive sequence

${\psi }_{\eta +1}=a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{{\psi }_{\eta -\left(2i+1\right)}}{\alpha {\psi }_{\eta -2i}+\beta {\psi }_{\eta -\left(2i+1\right)}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\eta =0,1,2,\cdots$ (1.1)

where $a,{b}_{i},\alpha$ and $\beta$ are positive real numbers for $i=0,1,\cdots ,k$ , and the the initial conditions ${\psi }_{-j},{\psi }_{-j+1},\cdots ,{\psi }_{0}$ are arbitrary positive real numbers where $j=2k+1$ .

In the next, we will and to many of the basic concepts. Before anything, the concept of equilibrium point is essential in the study of the dynamics of any physical system. A point $\stackrel{¯}{\psi }$ in the domain of the function $\Phi$ is called an equilibrium point of the equation

${\psi }_{\eta +1}=\Phi \left({\psi }_{\eta },{\psi }_{\eta -1},\cdots ,{\psi }_{\eta -k}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\eta =0,1,2,\cdots$ (1.2)

if $\stackrel{¯}{\psi }$ is a fixed point of $\Phi$ [ $\Phi \left(\stackrel{¯}{\psi },\stackrel{¯}{\psi },\cdots ,\stackrel{¯}{\psi }\right)=\stackrel{¯}{\psi }$ ]. For a stability of equilibrium point, equilibrium point $\stackrel{¯}{\psi }$ of equation (2) is said to be locally stable if for all $\epsilon >0$ there exists $\delta >0$ such that, if ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ with ${\sum }_{i=0}^{k}|{\psi }_{-i}-\stackrel{¯}{\psi }|<\delta$ . As well, $\stackrel{¯}{\psi }$ is said to be locally asymptotically stable if it is locally stable and there exists $\gamma >0$ such that, if ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ with ${\sum }_{i=0}^{k}|{\psi }_{-i}-\stackrel{¯}{\psi }|<\gamma$ , then ${\mathrm{lim}}_{\eta \to \infty }{\psi }_{\eta }=\stackrel{¯}{\psi }$ . Also, $\stackrel{¯}{\psi }$ is said to be a global attractor if used for every ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ , we have ${\mathrm{lim}}_{\eta \to \infty }{\psi }_{\eta }=\stackrel{¯}{\psi }$ . On the other hand, $\stackrel{¯}{\psi }$ is said to be unstable if it is not locally stable.

Finally, Equation (1.2) is called permanent and bounded if there exists numbers r and R with $0 such that for any initial conditions ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ there exists a positive integer N which depends on these initial conditions such as $r<{\psi }_{\eta } for all $\eta \ge N$ .

The linearized equation of Equation (1.1) about the equilibrium point $\stackrel{^}{\psi }$ is

${y}_{\eta +1}=\underset{i=0}{\overset{k}{\sum }}{p}_{i}{y}_{\eta -i}$ (1.3)

where

${p}_{i}=\frac{\partial F}{\partial {\psi }_{\eta -i}}\left(\stackrel{^}{\psi },\stackrel{^}{\psi },\cdots ,\stackrel{^}{\psi }\right).$

Theorem 1.1.  Assume that ${p}_{i}\in ℝ$ for $i=0,1,\cdots ,k$ . The $+ve$ equilibrium of (1.1) is locally asymptotically stable if

$|{p}_{0}|+|{p}_{1}|+\cdots +|{p}_{k}|<1.$ (1.4)

2. Local Stability of Equation (1.1)

The $+ve$ equilibrium point of Equation (1.1) is

$\stackrel{^}{\psi }=a+\underset{i=0}{\overset{k}{\sum }}\frac{{b}_{i}\stackrel{^}{\psi }}{\alpha \stackrel{^}{\psi }+\beta \stackrel{^}{\psi }},$

and so,

$\stackrel{^}{\psi }=a+\frac{B}{\alpha +\beta },$

where

$B=\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}.$

Let $f\in C\left({\left(0,\infty \right)}^{2k+2},\left(0,\infty \right)\right)$ defined by

$f\left({u}_{0},{u}_{1},\cdots ,{u}_{2k+1}\right)=a+\underset{i=0}{\overset{k}{\sum }}\frac{{b}_{i}{u}_{2i+1}}{\alpha {u}_{2i}+\beta {u}_{2i+1}}.$ (2.1)

Therefore it follows that

$\frac{\partial f}{\partial {u}_{2r}}=-\frac{{b}_{r}\alpha {u}_{2r+1}}{{\left(\alpha {u}_{2r}+\beta {u}_{2r+1}\right)}^{2}}$ (2.2)

and

$\frac{\partial f}{\partial {u}_{2r+1}}=\frac{{b}_{r}\alpha {u}_{2r}}{{\left(\alpha {u}_{2r}+\beta {u}_{2r+1}\right)}^{2}},$ (2.3)

for $r=0,1,\cdots ,k$ .

Theorem 2.1. Let $\stackrel{^}{\psi }$ be $+ve$ equilibrium of Equation (1.1). If

$\left(\alpha -\beta \right)B

than $\stackrel{^}{\psi }$ is locally stable.

Proof. From (2.2) to (2.3), we obtain

$\frac{\partial f}{\partial {u}_{2r}}\left(\stackrel{^}{\psi },\cdots ,\stackrel{^}{\psi }\right)=-\frac{{b}_{r}\alpha }{\left(\alpha +\beta \right)\left(a\left(\alpha +\beta \right)+B\right)}={p}_{2r},$

and

$\frac{\partial f}{\partial {u}_{2r+1}}\left(\stackrel{^}{\psi },\cdots ,\stackrel{^}{\psi }\right)=\frac{{b}_{r}\alpha }{\left(\alpha +\beta \right)\left(a\left(\alpha +\beta \right)+B\right)}={p}_{2r+1},$

for $r=0,1,\cdots ,k$ . Thus, the linearized equation of (1.1) is

${y}_{\eta +1}={p}_{0}{y}_{\eta }+{p}_{1}{y}_{\eta -1}+\cdots +{p}_{2k+1}{y}_{\eta -\left(2k+1\right)}.$

It follows by Theorem 1.1 that Equation (1.1) is locally stable if

$|\frac{{b}_{0}\alpha }{\left(\alpha +\beta \right)}|+|\frac{{b}_{0}\alpha }{\left(\alpha +\beta \right)\rho }|+\cdots +|\frac{{b}_{k}\alpha }{\left(\alpha +\beta \right)\rho }|+|\frac{{b}_{k}\alpha }{\left(\alpha +\beta \right)\rho }|<1,$

where $\rho =\left(a\left(\alpha +\beta \right)+B\right)$ , and hence,

$\frac{2\alpha }{\left(\alpha +\beta \right)\rho }B<1.$

Thus, we find

$2\alpha B<\left(\alpha +\beta \right)\rho ,$

and so,

$\left(\alpha -\beta \right)B

Hence, the proof is complete.

In order to verify and support our theoretical outcomes and discussions, in this concern, we investigate several interesting numerical examples.

Example 2.1. By Theorem 2.1, the $+ve$ equilibrium Equation (1.1) with $a=2$ , $k=2$ , ${b}_{i}=0.3$ , $\beta =0.1$ and $\alpha =1$ , is locally stable (see Figure 1).

3. Global Stability of Equation (1.1)

In the following theorem, we check into the global stability of the recursive sequence (1.1).

Theorem 3.1. The $+ve$ equilibrium $\stackrel{^}{\psi }$ of Equation (1.1) is global attractor if

$B=a\left(\alpha -\beta \right).$

Figure 1. The stable solution corresponding to difference Equation (1.1).

Proof. We consider the function as follow:

$f\left({u}_{0},{u}_{1},\cdots ,{u}_{2k+1}\right)=a+\underset{i=0}{\overset{k}{\sum }}\frac{{b}_{i}{u}_{2i+1}}{\alpha {u}_{2i}+\beta {u}_{2i+1}}.$

From (2.2) and (2.3), we note that $f$ is increasing in ${u}_{2r+1}$ and decreasing in ${u}_{2r}$ for all $r=0,1,\cdots ,k$ . Suppose that $\left(\lambda ,\mu \right)$ is a solution of the system

$\lambda =f\left(\mu ,\lambda ,\mu ,\lambda ,\cdots ,\lambda \right)$

$\mu =f\left(\lambda ,\mu ,\lambda ,\mu ,\cdots ,\mu \right).$

Then, we find

$\lambda =a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{\lambda }{\alpha \mu +\beta \lambda },$

and

$\mu =a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{\mu }{\alpha \lambda +\beta \mu }.$

Hence, we get

$\lambda =a+\frac{B\lambda }{\alpha \mu +\beta \lambda },$ (3.1)

and

$\mu =a+\frac{B\mu }{\alpha \lambda +\beta \mu }.$ (3.2)

By (3.1) and (3.2), we obtain

$\beta \left({\lambda }^{2}-{\mu }^{2}\right)-\left(a\beta +B-a\alpha \right)\left(\lambda -\mu \right)=0.$

Thus,

$\left(\lambda -\mu \right)\left(\beta \left(\lambda +\mu \right)-\left(a\beta +B-a\alpha \right)\right)=0.$

Since $B=a\left(\alpha -\beta \right)$ , we have that $\mu =\lambda$ . Hence, the proof of Theorem 3.1 is complete. □

4. Periodic Solutions

In this section, we enumerate some basic facts concerning the existence of two period solutions.

Theorem 4.1. Equation (1.1) has prime period-two solutions if

$\left(a\beta +B-a\alpha \right)\left(\alpha -\beta \right)>4a\alpha \beta .$ (4.1)

Proof. Assume that Equation (1.1) has a prime period-two solution

$\cdots ,\rho ,\sigma ,\rho ,\sigma ,\rho ,\sigma ,\cdots$

We shall prove that condition (4.1) holds. From Equation (1.1), we see that

${\psi }_{\eta }={\psi }_{\eta -2}=\cdots ={\psi }_{\eta -2k}=\sigma ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\psi }_{\eta +1}={\psi }_{\eta -1}=\cdots ={\psi }_{\eta -\left(2k+1\right)}=\rho ,$

and hence,

$\rho =a+\frac{B\rho }{\alpha \sigma +\beta \rho },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sigma =a+\frac{B\sigma }{\alpha \rho +\beta \sigma }.$

Thus, we get

$\alpha \rho \sigma +\beta {\rho }^{2}=a\alpha \sigma +a\beta \rho +B\rho ,$ (4.2)

and

$\alpha \rho \sigma +\beta {\rho }^{2}=a\alpha \sigma +a\beta \sigma +B\sigma .$ (4.3)

From (4.3) and (4.2), we have

$\beta \left({\rho }^{2}-{\sigma }^{2}\right)=a\alpha \left(\sigma -\rho \right)+a\beta \left(\rho -\sigma \right)+B\left(\rho -\sigma \right).$

Dividing $\left(\rho -\sigma \right)$ , then we find

$\left(\rho +\sigma \right)=\frac{a\beta +B-a\alpha }{\beta }$ (4.4)

By combining (4.2) and (4.3), we obtain

$2\alpha \rho \sigma +\beta \left({\rho }^{2}+{\sigma }^{2}\right)=\left(a\alpha +a\beta +B\right)\left(\rho +\sigma \right).$

Since ${\rho }^{2}+{\sigma }^{2}={\left(\rho +\sigma \right)}^{2}-2\rho \sigma$ , we get

$\rho \sigma =\frac{a\alpha \left(a\beta +B-a\alpha \right)}{\beta \left(\alpha -\beta \right)}.$ (4.5)

Now, evident is that (4.4) and (4.5) that $\rho$ and $\sigma$ are both two positive distinct roots of the quadratic equation

${u}^{2}+\left(\rho +\sigma \right)u+\rho \sigma =0.$ (4.6)

Hence, we obtain

$\frac{\left(a\beta +B-a\alpha \right)}{\beta }>\frac{4a\alpha }{\alpha -\beta },$

which has the same extent as

$\left(a\beta +B-a\alpha \right)\left(\alpha -\beta \right)>4a\alpha \beta .$

Hence, the proof is complete. □

The next numerical example is mimicry to enhance our results.

Example 4.1. By Theorem 4.1, Equation (1.1) with $\alpha =0.02$ , $\beta =0.01$ , $a=500$ , ${b}_{0}=2$ , ${b}_{1}=200$ and ${b}_{2}=20$ , has prime period two solution (see Figure 2)

5. Boundedness

Theorem 5.1. Every solution of Equation (1.1) is bounded and persists.

Proof. Let ${\left\{{\psi }_{\eta }\right\}}_{\eta =-k}^{\infty }$ be a Solution (1.1), we can conclude from (1.1) that

${\psi }_{\eta +1}=a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{{\psi }_{\eta -\left(2i+1\right)}}{\alpha {\psi }_{\eta -2i}+\beta {\psi }_{\eta -\left(2i+1\right)}}>a.$

Then

${\psi }_{\eta }>a\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{all}\text{\hspace{0.17em}}\eta >-\left(2k+1\right).$

Also, from Equation (1.1), we see that

$\begin{array}{c}{\psi }_{\eta +1}=a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{{\psi }_{\eta -\left(2i+1\right)}}{\alpha {\psi }_{\eta -2i}+\beta {\psi }_{\eta -\left(2i+1\right)}}\\

then,

$a<{\psi }_{\eta }\le a+\frac{B}{\beta }\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{all}\text{\hspace{0.17em}}\eta >-\left(2k+1\right).$

Thus, the solution is bounded and persists and the proof is complete. □

Conclusion 1. In this paper, we study a asymptotic behavior of solutions of a general class of difference Equation (1.1). Our results extend and generalize to the earlier ones. Moreover, we obtain the next results:

- The $+ve$ equilibrium point $\stackrel{^}{\psi }$ of Equation (1.1) is local stable if $\left(\alpha -\beta \right)B . Also, if $B=a\left(\alpha -\beta \right)$ , then $\stackrel{^}{\psi }$ is global attractor.

Figure 2. Prime period two solution of Equation (1.1).

- Equation (1.1) has a prime period-two solutions if $\left(a\beta +B-a\alpha \right)\left(\alpha -\beta \right)>4a\alpha \beta$ .

- Every solution of (1.1) is bounded and persists.

Acknowledgements

The author is very grateful to the reviewers for their valuable suggestions and useful comments on this paper.

Cite this paper: Moaaz, O. , Elabbasy, E. and Alsaeed, S. (2017) Global Analysis of Solutions of a New Class of Rational Difference Equation. American Journal of Computational Mathematics, 7, 495-503. doi: 10.4236/ajcm.2017.74036.
References

   Abdelrahman, M.A.E. and Moaaz, O. (2017) Investigation of the New Class of the Nonlinear Rational Difference Equations. Fundamental Research and Development International, 7, 59-72.

   Abdelrahman, M.A.E. and Moaaz, O. (2018) On the New Class of the Nonlinear Difference Equations. Electronic Journal of Mathematical Analysis and Applications, 6, 117-125.

   Agarwal, R. (1992) Difference Equations and Inequalities. Tteory, Methods and Applications, Marcel Dekker Inc., New York.

   Amleh, A.M., Grove, E.A., Georgiou, A. and Ladas, G. (1999) On the Recursive Sequence ωη+1 = α + ωη-1η. The Journal of Mathematical Analysis and Applications, 233, 790-798.
https://doi.org/10.1006/jmaa.1999.6346

   Berenhaut, K.S. and Stevic, S. (2006) The Behaviour of the Positive Solutions of the Difference Equation ψη = a + (ψη-2η-1)σ. Journal of Difference Equations and Applications, 12, 909-918.
https://doi.org/10.1080/10236190600836377

   Camouzis, E., DeVault, R. and Ladas, G. (2001) On the Recursive Sequence ωη+1 = -1 + ωη-1η. Journal of Difference Equations and Applications, 7, 477-482.
https://doi.org/10.1080/10236190108808283

   Elabbasy, E.M., EL-Metwally, H.A. and Elsayed, E.M. (2006) On the Difference Equation ωη+1 = aωη - bωη/(cωη - dωη-1). Advances in Difference Equations, 2006, 1-10.

   Elaydi, S.N. (1996) An Introduction to Difference Equations, Undergraduate Texts in Mathematics. Springer, New York.
https://doi.org/10.1007/978-1-4757-9168-6

   Elsayed, E.M. and El-Metwally, H.A. (2012) Qualitative Studies of Scalars and Systems of Difference Equations. Publisher LAP Lambert Academic, Saarbrücken.

   Elsayed, E.M. (2015) New Method to Obtain Periodic Solutions of Period Two and Three of a Rational Difference Equation. Nonlinear Dynamics, 79, 241-250.
https://doi.org/10.1007/s11071-014-1660-2

   Elsayed, E.M. (2016) Dynamics and Behavior of a Higher Order Rational Difference Equation. The Journal of Nonlinear Sciences and Applications, 9, 1463-1474.

   Elsayed, E.M. (2010) On the Global Attractivity and the Periodic Character of a Recursive sequence. Opuscula Mathematica, 30, 431-446.
https://doi.org/10.7494/OpMath.2010.30.4.431

   Grove, E.A. and Ladas, G. (2005) Periodicities in Nonlinear Difference Equations, Vol. 4. Chapman & Hall/CRC, London.

   Khuong, V.V. (2008) On the Positive Nonoscillatory Solution of the Difference Equations ωη+1 = α + (ωη-kη-r)p. Applied Mathematics-A Journal of Chinese Universities, 24, 45-48.

   Kulenovic, M.R.S. and Ladas, G. (2001) Dynamics of Second Order Rational Difference Equations with Open Problems and Conjectures. Chapman & Hall/CRC, Florida.
https://doi.org/10.1201/9781420035384

   El-Metwally, H., Ladas, G., Grove, E.A. and Voulov, H.D. (2001) On the Global Attractivity and the Periodic Character of Some Difference Equations. Journal of Difference Equations and Applications, 7, 837-850.
https://doi.org/10.1080/10236190108808306

   Metwally, H.E., Grove, E.A. and Ladas, G. (2000) A Global Convergence Result with Applications to Periodic Solutions. Journal of Mathematical Analysis and Applications, 245, 161-170.
https://doi.org/10.1006/jmaa.2000.6747

   Moaaz, O. (2016) Comment on New Method to Obtain Periodic Solutions of Period Two and Three of a Rational Difference Equation. Nonlinear Dynamics, 79, 241-250.
https://doi.org/10.1007/s11071-016-3293-0

   Moaaz, O. and Abdelrahman, M.A.E. (2016) Behaviour of the New Class of the Rational Difference Equations. Journal of Mathematical Analysis and Applications, 4, 129-138.

   Moaaz, O. and Abdelrahman, M.A.E. (2017) On the Class of the Rational Difference Equations. International Journal of Advances in Mathematics, 2017, 46-55.

   Ocalan, O. (2014) Dynamics of the Difference Equation ωη+1 = Pη + ωη-kη with a Period-Two Coefficient. Applied Mathematics and Computation, 228, 31-37.
https://doi.org/10.1016/j.amc.2013.11.020

   Saleh, M. and Aloqeili, M. (2005) On the Difference Equation ψη+1 = A + ψηη-k. Applied Mathematics and Computation, 171, 862-869.
https://doi.org/10.1016/j.amc.2005.01.094

   Stevic, S. (2005) On the Recursive Sequence ψη+1 = a + βψη-k/f(ψη,···, ψη-k+1). Taiwanese Journal of Mathematics, 9, 583-593.

   Sun, T. and Xi, H. (2007) On Convergence of the Solutions of the Difference Equation ωη+1 = 1 + ψηη-k. Journal of Mathematical Analysis and Applications, 325, 1491-1494.
https://doi.org/10.1016/j.jmaa.2006.02.080

   Zhang, L., Zhang, G. and Liu, H. (2005) Periodicity and Attractivity for a Rational Recursive Sequence. Journal of Applied Mathematics and Computing, 19, 191-201.
https://doi.org/10.1007/BF02935798

Top