We study the global stability, the periodic character, and the boundedness character of the positive solutions of the difference equation xn+1=(α−βxn)/(γ−δxn−xn−k),n=0,1,2,…,k∈{1,2,…}, in the two cases: (i) δ≥0,α>0,γ>β>0; (ii) δ≥0,α=0,γ,β>0, where the coefficients α,β,γ,andδ, and the initial conditions x−k,x−k+1,…,x−1,x0 are real numbers. We show that the positive equilibrium of this equation is a global attractor with a basin that depends on certain conditions posed on the coefficients of this equation.
1. Introduction
The asymptotic
stability of the rational recursive sequence, xn+1=α+βxnγ+∑i=0kγixn−i,n=0,1,2,…, was investigated
when the coefficients α,β,γ,andγi are nonnegative real numbers (see [1–3]). Studying the asymptotic behavior
of the rational sequence (1.1) when some of the coefficients are negative was
suggested in [3]. Recently, Aboutaleb et al. [4] studied the rational recursive
sequence, xn+1=α−βxnγ+xn−1,n=0,1,2,…,where α,β,andγ are nonnegative real numbers and obtained
sufficient conditions for the global attractivity of the positive equilibria.
Yan et al. [5] studied recently the rational recursive
sequence, xn+1=α−βxnγ−xn−k,n=0,1,2,…, where α≥0,γ,β>0 are real numbers while k≥1 is an integer number, and the initial
conditions x−k,x−k+1,…,x−1,x0 are arbitrary real numbers. They proved that the positive equilibrium x¯ of (1.3) is a global attractor with a basin
that depends on certain conditions of the coefficients. He et al. [6]
studied recently the rational recursive sequence, xn+1=a−bxn−kA+xn,n=0,1,2,…, where a≥0,A,b>0 are real numbers while k≥1 is an integer number and the initial
conditions x−k,x−k+1,…,x−1,x0 are arbitrary real numbers. They proved the global attractivity and
periodic character of the positive solution of (1.4). Stević [7] studied
recently the rational recursive sequence, xn+1=α+βxnγ−xn−k,n=0,1,2,…, where the parameters α,β,andγ are nonnegative real numbers and k≥1 is an integer number while the initial
conditions x−k,x−k+1,…,x−1,x0 are arbitrary real numbers. Other related results can be found in
[8–19].
Our aim in this paper is to study the global
attractivity, the periodicity, and the boundedness of the positive solution of
the following rational recursive sequence: xn+1=α−βxnγ−δxn−xn−k,n=0,1,2,…, in the two cases (i) δ≥0,α>0,γ>β>0, (ii) δ≥0,α=0,γ,β>0, where the coefficients α,β,γ,andδ are real numbers and k≥1 is an integer number, while the initial
conditions x−k,x−k+1,…,x−1,x0 are arbitrary real numbers. We will prove that the positive
equilibrium x¯ of (1.6) is a global attractor with a basin
that depends on certain conditions of these coefficients.
2. Local Stability and Permanence
We first recall some results which will be useful in
the sequel. Let I be some real interval and let F be a
continuous function defined on Ik+1. Then, for initial conditions x−k,x−k+1,…,x−1,x0∈I, it is easy to see that the difference equation,xn+1=F(xn,xn−1,…,xn−k),n=0,1,2,…,k≥1, has a unique solution {xn}.
Definition 2.1.
A point x¯ is called an equilibrium of (2.1), if x¯=F(x¯,…,x¯). That is, xn=x¯ for n≥0 is a solution of (2.1), or equivalently, is a
fixed point of F.
Definition 2.2.
An interval J⊂I is called an invariant interval of (2.1) if
the initial conditions x−k,x−k+1,…,x−1,x0∈J imply that the solution xn∈J for n>0. That is, every solution of (2.1) with initial conditions in J remains in J.
Definition 2.3.
The difference equation (2.1) is said to be
permanent if there exist numbers P and Q with 0<P≤Q<∞ such that for any initial conditions x−k,x−k+1,…,x−1,x0 there exists a positive integer N which depends on the initial
conditions such that P≤xn≤Q, for all n≥N.
The linearized equation associated with (2.1) about
the equilibrium x¯ is yn+1=∑i=0k∂F(x¯,…,x¯)∂uiyn−i,n=0,1,2,…. Its characteristic equation
is λn+1=∑i=0k∂F(x¯,…,x¯)∂uiλn−i,n=0,1,2,….
Theorem 2.4 (see [3]).
Assume that F is a C1-function and let x¯ be an equilibrium of (2.1). Then, the
following statements are true:
if all the
roots of (2.3) lie in the open unit disk |λ|<1, then the equilibrium x¯ of (2.1) is locally asymptotically stable;
if at least
one root of (2.3) has absolute value greater than one, then the equilibrium x¯ of (2.1) is unstable.
Theorem 2.5 (see [3, 8]).
Assume that p,q∈ℝ,
and k∈{1,2,…}. Then, |p|+|q|<1 is a sufficient condition for the
asymptotic stability of the difference equation xn+1−pxn+qxn−k=0,n=0,1,2,…. Suppose in addition that one of
the following two cases holds: (i) k is odd and q<0, or (ii)
k is even and pq<0. Then, (2.4) is also a necessary condition for
the asymptotic stability of (2.5) (see [6]).
First, we study the rational recursive
sequence xn+1=α−βxnγ−δxn−xn−k,n=0,1,2,…, together with the
conditions δ≥0,α>0,γ>β>0,k∈{1,2,…}.
The unique positive equilibrium point x¯ of (2.6) is the solution of the
equation x¯=α−βx¯γ−(δ+1)x¯, which is given by x¯=(γ+β)±T2(δ+1), where T=(γ+β)2−4α(δ+1). If (2.7) holds and α=(γ+β)2/4(δ+1), then (2.6) has a unique positive equilibrium x¯0=(γ+β)/2(δ+1). If (2.7) holds and α<(γ+β)2/4(δ+1) then (2.6) has two positive equilibria x¯1,2 given by (2.9).
The linearized equation of (2.6) about the equilibrium x¯i(i=0,1,2) is given by yn+1+β−δx¯i[γ−(δ+1)x¯i]yn−x¯i[γ−(δ+1)x¯i]yn−k=0. The characteristic equation
associated with (2.6) about x¯0 is λk+1+[2βγ−β−δ(γ+β)(δ+1)(γ−β)]λk−γ+β(δ+1)(γ−β)=0. Now, we have the following
results:
if 0≤δ<2β/(γ−β),
then (γ+β)/(δ+1)(γ−β)>1 and hence the equilibrium x¯0 of (2.6) is unstable (see Figure 1);
if δ>2β/(γ−β), then |2βγ−β−δ(γ+β)(δ+1)(γ−β)|+|γ+β(δ+1)(γ−β)|=1.
Thus, the linearized stability
analysis fails. On the other hand, the characteristic equation associated with
(2.6) about x¯1 is λk+1+[2βγ−β−T−δ(γ+β+T)(δ+1)(γ−β−T)]λk−γ+β+T(δ+1)(γ−β−T)=0. Now, we have the following results:
if 0≤δ<2β/(γ−β),
then it is obvious that |γ+β+T(δ+1)(γ−β−T)|≥(γ−βγ+β)[γ+β+T][γ−β−T]>1, hence the equilibrium x¯1 of (2.6) is unstable;
if δ≥2β/(γ−β),
then it is easy to see that 2β(δ+1)<δ[γ+β+T], and consequently, we have |2βγ−β−T−δ[γ+β+T](δ+1)[γ−β−T]|+|γ+β+T(δ+1)[γ−β−T]|=γ−β+Tγ−β−T>1, and hence the equilibrium x¯1 of (2.6) is unstable.
For the positive equilibrium x¯2, in view of conditions (2.7) and α<(γ+β)2/4(δ+1),
we have x¯2=γ+β−T2(δ+1)<γ+β2(δ+1)<γδ+1. Hence, if 0<α≤β(γ−β)δ+1, then T≥(γ+β)2−4β(γ−β)>(γ+β)2−(γ+3β)(γ−β)=(γ+β)2−(γ+β)2+4β2=2β. Consequently, we
have |β−δx¯2γ−(δ+1)x¯2|+|x¯2γ−(δ+1)x¯2|<β+(δ+1)x¯2γ−(δ+1)x¯2=3β+γ−Tγ−β+T<3β+γ−2βγ−β+2β=1, which by Theorem 2.5 implies that x¯2 is locally asymptotically stable (see Figure 2).
Lemma 2.6.
Let f(u,v)=(α−βu)/(γ−δu−v) and assume that conditions (2.7) and
(2.18) hold. Then, the following statements are true:
0<x¯2<α/β,α/β<x¯1<∞;
f(x,x) is a
strictly decreasing function in (−∞,α/β);
let u,v∈(−∞,α/β), then the function f(u,v) is a strictly
decreasing function in u and a strictly increasing function in v.
Proof.
We prove (a) only. The proofs of
(b) and (c) are omitted here. In view of (2.7) and (2.18), we
have x¯2=γ+β−T2(δ+1)<γ+β2(δ+1)<γδ+1. From (2.8) and (2.21), we have α−βx¯2>0 and so x¯2<α/β. Also, in view of (2.7) and (2.18), we have 0<α−βx¯1γ−(δ+1)x¯1=x¯1=γ+β+T2(δ+1)≥γ+β+(γ+β)2−4β(γ−β)2(δ+1)=γ+β+(γ−β)2+4β22(δ+1)>γ+β+(γ−β)22(δ+1)=γδ+1, and so γ−x¯1(δ+1)<0. Consequently, α−βx¯1<0 which implies that x¯1>α/β. The proof is completed.
Theorem 2.7.
Assume that the conditions (2.7)
and (2.18) hold. Let {xn} be any solution of (2.6). If xi∈(−∞,α/β],for i=−k,−k+1,…,−1 and if x0∈[0,α/β], then 0≤xn≤αβ,n=1,2,…. That is the solution {xn} is bounded.
Proof.
By part (c) of Lemma 2.6, we
have 0=α−β⋅(α/β)γ−δx0−x−k≤x1=α−βx0γ−δx0−x−k≤α−β⋅0γ−δ⋅0−α/β=αγ−α/β=βαγβ−α. From (2.18), we deduce that γβ−α>β2,
and then we have 0≤x1≤αβ. Also, we have 0=α−β⋅(α/β)γ−δx1−x−k+1≤x2=α−βx1γ−δx1−x−k+1≤α−β⋅0γ−δ⋅0−α/β<αβ. Thus, 0≤x2≤αβ. The result (2.23) now follows by
induction. The proof is completed.
3. Global Attractivity
In this section, we will study the global attractivity
of positive solutions of (2.6). We show that the positive equilibrium x¯ of (2.6) is a global attractor with a basin
that depends on certain conditions imposed on the coefficients.
Theorem 3.1.
Assume
that the conditions (2.7) and (2.18) hold. Then, the equilibrium point x¯2 of (2.6) is globally asymptotically stable.
Proof.
In Section 2, we have shown under the assumptions (2.7) and (2.18) that the
equilibrium x¯2 is locally asymptotically stable. It remains to prove that the
equilibrium x¯2 is a global attractor. To this end, set I=limn→∞infxn and S=limn→∞supxn which by Theorem 2.7 exist and are positive numbers. Then, from (2.6) we
deduce that S≤α−βSγ−(δ+1)I,I≥α−βIγ−(δ+1)S. Consequently, we have −α+(γ+β)S≤(δ+1)IS≤−α+(γ+β)I, from which it follows that I=S. Thus, the the proof of Theorem 3.1 is completed.
Lemma 3.2 (see [8]).
Consider the difference
equation xn+1=f(xn,xn−k),k≥1,n=0,1,2,…. Let [a,b] be some interval of real numbers, and assume
that f:[a,b]×[a,b]→[a,b] is a continuous function satisfying the
following properties:
f(u,v) is a nonincreasing function in u, and a
nondecreasing function in v;
if (m,M)∈[a,b]×[a,b] is a solution of the system m=f(M,m),M=f(m,M), then, m=M.
Then, (3.3) has
a unique equilibrium point x¯ and every solution of (3.3) converges to x¯.
Theorem 3.3.
Assume
that the conditions (2.7) and (2.18) hold. Then, the positive equilibrium x¯ of (2.6) is a global attractor with a basin S*=[0,α/β]k+1.
Proof.
For u,v∈[0,α/β], set f(u,v)=α−βuγ−δu−v. We claim that f:[0,α/β]×[0,α/β]→[0,α/β]. In fact, if we set a=0,b=α/β,
then f(b,a)=α−βbγ−δb−a=α−αγ−δ(α/β)=0=a, and in view of the condition
(2.18), we have f(a,b)=α−βaγ−δa−b=αγ−α/β=βαγβ−α<αβ=b. Since f(u,v) is decreasing in u and increasing in v,
it follows that a≤f(u,v)≤b, for all u,v∈[a,b], which implies that our assertion is true. On
the other hand, conditions (a) and (b) of Lemma 3.2 are clearly true. Let {xn} be a solution of (2.6) with the initial conditions (x−k,x−k+1,…,x−1,x0)∈S. By Lemma 3.2, we have limn→∞xn=x¯. The proof is completed.
Theorem 3.4.
Assume that conditions (2.7)
and (2.18) hold. Then, the positive equilibrium x¯ of (2.6) is a global attractor with a basin S∗=(−∞,α/β]k×[0,α/β].
Proof.
Let {xn} be a solution of (2.6) with the initial conditions x−k,x−k+1,…,x−1,x0∈S*. Then, by Theorem 2.7, we have xn∈[0,αβ],n=1,2,…. By Theorem 3.3, we have limn→∞xn+k=x¯ and so limn→∞xn=x¯. The proof is completed.
Theorem 3.5.
Assume that conditions (2.7)
hold with 0≤δ<1.
Also, assume that k is an odd positive integer. Then, the necessary and
sufficient condition for (2.6) to have positive solutions of prime period two
is that β(γ−β)<α<(γ−β)4[γ+3β−δ(γ−β)].
Proof.
First, suppose that there exist distinctive
positive solutions of prime period two, …,P,Q,P,Q,…, of the difference equation
(2.6).
If k is odd, then xn+1=xn−k. It follows from the difference equation (2.6) that P=α−βQγ−δQ−P,Q=α−βPγ−δP−Q. Consequently, we
obtain P+Q=γ−β,PQ=α−β(γ−β)1−δ. Thus, we deduce
that α>β(γ−β),0≤δ<1. Now it is clear that P,Q are two
positive distinct real roots of the quadratic equation t2−(P+Q)t+PQ=0. Therefore, we
have (γ−β)2>4[α−β(γ−β)]1−δ. From (3.13) and (3.15) we obtain
condition (3.9). Conversely, suppose that the condition (3.9) is valid.
Then, we deduce that (3.13) and (3.15) hold. Consequently, there exists two
positive distinct real numbers P and Q such that P=γ−β2−K2,Q=γ−β2+K2, where K>0 is given by K=(γ−β)2−4[α−β(γ−β)1−δ],0≤δ<1. Thus, P and Q given by (3.16) and (3.17) represent two
positive distinct real roots of the quadratic equation (3.14). Now, we are going
to prove that P and Q given by (3.16) and (3.17) are positive
solutions of prime period two of the difference equation (2.6). To this end, we
assume that x−k=P,x−k+1=Q,…,x−1=P,x0=Q. We wish to prove that x1=P and x2=Q.
It follows from the difference equation (2.6) and the
formulas (3.16) and (3.17) that x1=α−βx0γ−δx0−x−k=α−βQγ−δQ−P=2α−β[γ−β+K]2γ−δ[γ−β+K]−[γ−β−K]=β[2α/β−(γ−β)−K]2γ−(1+δ)(γ−β)+(1−δ)K=(β1−δ){2α/β−(γ−β)−K}{2γ/(1−δ)−((1+δ)/(1−δ))(γ−β)−K}{2γ/(1−δ)−((1+δ)/(1−δ))(γ−β)+K}{2γ/(1−δ)−((1+δ)/(1−δ))(γ−β)−K}. After some reduction, we deduce
that x1=(β1−δ)[α(1−δ)+β2][γ−β−K]/β(1−δ)2[α(1−δ)+β2]/(1−δ)2=γ−β−K2=P. Similarly, we can show
that, x2=α−βx1γ−δx1−x−k+1=α−βPγ−δP−Q=Q. By using the induction, we
have xn=P,xn+1=Q,∀n≥−k. Thus, the difference equation
(2.6) has positive solutions of prime period two. Hence, the proof of Theorem 3.5
is completed.
Theorem 3.6.
Assume that conditions (2.7)
hold. If k is even, then (2.6) has no positive solutions of prime period
two.
Proof.
Suppose that there exists distinctive positive solutions of prime
period two, …,P,Q,P,Q,…, of the difference equation
(2.6).
If k is even, then xn=xn−k. It follows from the difference equation (2.6) that P=α−βQγ−(δ+1)Q,Q=α−βPγ−(δ+1)P. From which we have (γ−β)(P−Q)=0 and by using (2.7), we deduce that P=Q. This is a contradiction. Thus, the proof of
Theorem 3.6 is completed.
4. The Case α=0
Secondly, we study the rational recursive
sequence xn+1=−βxnγ−δxn−xn−k,n=0,1,2,…, where δ≥0,γ,β>0 are real numbers and k∈{1,2,…}. By putting xn=βyn, (4.1) yields yn+1=−ynA−δyn−yn−k,n=0,1,2,…, where A=γ/β.
Equation (4.2) has two equilibrium points y¯1=0,y¯2=1+A1+δ. The linearized equation
associated with (4.2) about the equilibria y¯i,(i=1,2) is zn+1+1−δy¯iA−(δ+1)y¯izn−y¯iA−(δ+1)y¯izn−k=0. The characteristic equation of
(4.4) about the equilibrium y¯2=(1+A)/(1+δ) is λk+1+[δA−1δ+1]λk+A+1δ+1=0. Now, we deduce from (4.5) the
following results:
if δ=0, and since A+1>1, then the equilibrium y¯2 is unstable (see [7]);
if A>δ>0, and since (A+1)/(δ+1)>1, then the equilibrium y¯2 is unstable;
if A=δ, then |δA−1δ+1|+|δ+1δ+1|=|δ−1|+1.
Now, we have the following
results from case (c): (i) if A=δ>1,
then the equilibrium y¯2 is unstable; (ii) if 0<A=δ<1,
then the equilibrium y¯2 is unstable; (iii) if A=δ=1,
then the linearized stability analysis fails;
if 1<A<δ, |δA−1δ+1|+|A+1δ+1|=δA−1δ+1+A+1δ+1=A(δ+1)δ+1=A>1, and hence the equilibrium y¯2 is unstable;
if A<δ≤1,|δA−1δ+1|+|A+1δ+1|=1−δA+1+Aδ+1≥1+A2(1−δ)≥1, and hence the equilibrium y¯2 is unstable.
The
characteristic equation of (4.4) about the equilibrium y¯1=0 is λk+1+1Aλk=0. This equation has two
roots λ1=0,λ2=−1A. Now, we deduce from (4.10) the
following results:
if A>1, then the equilibrium y¯1=0 is locally asymptotically stable (see Figure 3);
if 0<A<1, then the equilibrium y¯1=0 is unstable (see Figure 4);
if A=1, then the linearized stability analysis fails.
In the following results, we assume that A≥δ+2, where δ≥0.
Lemma 4.1.
Assume that the initial conditions y−i∈[−1,1], for i=1,2,…,k and y0∈[−1,0]. Then, {y2n−1} is nonnegative and monotonically decreasing to zero, while {y2n} is nonpositive and monotonically increasing to zero.
Proof.
Suppose
that y−i∈[−1,1], for i=1,2,…,k and y0∈[−1,0]. Clearly, 0≤y1≤1 and −1≤y2≤0. By induction, we can see that 0≤y2n−1≤1 and −1≤y2n≤0 for n≥1.
If A≥δ+2,δ≥0 we have y2n−1y2n+1=(A−δy2n−y2n−k)(A−δy2n−1−y2n−k−1)>1, and hence y2n−1>y2n+1,n=1,2,…. Similarly, we can show that y2n<y2n+2,n=1,2,…. The proof of Lemma 4.1 is completed.
On using arguments similar to that used in Lemma 4.1, we
can easily prove the following lemma.
Lemma 4.2.
Assume that the initial conditions y−i∈[−1,1], for i=1,2,…,k and y0∈[0,1]. Then, {y2n−1} is nonpositive and monotonically increasing to zero, while {y2n} is nonnegative and monotonically decreasing to zero.
Corollary 4.3.
The equilibrium point y¯1=0 of (4.1) is a global attractor with a basin S*=[−1,1]k+1.
Theorem 4.4.
The
equilibrium point y¯1=0 of (4.1) is a global attractor with a basin S*=(−∞,1]k×[(−A+1)/(δ+1),(A−1)/(δ+1)], where δ≥0.
Proof.
Assuming that the initial conditions y−k,y−k+1,…,y−1,y0∈S*. If A≥δ+2, with δ≥0, then we deduce that −1≤(1−A)/(1+δ)A−δy0−y−k≤y1=−y0A−δy0−y−k≤(A−1)/(δ+1)(A−1)/(δ+1)=1,−1≤−1A−δy1−y−k+1≤y2=−y1A−δy1−y−k+1≤1A−δ−1≤1. By induction, it follows that yi∈[−1,1] for i≥1. Thus, the proof of Theorem 4.4 follows from Corollary 4.3.
Theorem 4.5.
If A>1, then the equilibrium point y¯1=0 of (4.2) is globally asymptotically stable.
Finally, on using arguments similar to that used in
Theorems 3.5 and 3.6, we can prove easily the following results.
Theorem 4.6.
Assume
that δandA>1. If k is an odd positive integer, then the
necessary and sufficient condition for (4.2) to have positive solutions of
prime period two is that (see Figure 5) (A−1)δ>A+3.
Theorem 4.7.
If k is an even positive integer,
then (4.2) has no positive solutions of prime period two (see Figure 6).
KocićV. L.LadasG.Global attractivity in a second-order nonlinear difference equation19931801144150MR125028510.1006/jmaa.1993.1390ZBL0802.39001KocićV. L.LadasG.RodriguesI. W.On rational recursive sequences19931731127157MR120591410.1006/jmaa.1993.1057ZBL0777.39002KocićV. L.LadasG.1993256Dordrecht, The NetherlandsKluwer Academic Publishersxii+228Mathematics and Its ApplicationsMR1247956ZBL0787.39001AboutalebM. T.El-SayedM. A.HamzaA. E.Stability of the recursive sequence xn+1=(α−βxn)/(γ+xn−1)20012611126133MR185096110.1006/jmaa.2001.7481ZBL0990.39009YanX.-X.LiW. T.SunH.-R.Global attractivity in a higher order nonlinear difference equation200225158MR1979410ZBL1004.39010HeW.-S.LiW.-T.Attractivity in a nonlinear delay difference equation200444853MR2077780ZBL1071.39007StevićS.On the recursive sequence xn+1=(α+βxn)/(γ−xn−k)20043216170MR2037745DeVaultR.KosmalaW.LadasG.SchultzS. W.Global behavior of yn+1=(p+yn−k)/(qyn+yn−k)200147747434751MR1975867ZBL1042.3952310.1016/S0362-546X(01)00586-7El-OwaidyH. M.El-AfifiM. M.A note on the periodic cycle of xn+2=(1+xn+1)/(xn)20001092-3301306MR173819510.1016/S0096-3003(99)00031-4ZBL1023.39010GibbonsC. H.KulenovićM. R. S.LadasG.VoulovH. D.On the trichotomy character of xn+1=(α+βxn+γxn−1)/(A+xn)2002817592MR188459310.1080/10236190211940ZBL1005.39017GroveE. A.LadasG.PredescuM.RadinM.On the global character of yn+1=(pyn−1+yn−2)/(q+yn−2)2001572539MR1883898ZBL0994.39012KulenovićM. R. S.LadasG.ProkupN. R.On the recursive sequence xn+1=(αxn+βxn−1)/(A+xn)200065563576MR180244710.1080/10236190008808246ZBL0966.39003KulenovićM. R. S.LadasG.ProkupN. R.A rational difference equation2001415-6671678MR1822594ZBL0985.39017KuruklisS. A.The asymptotic stability of xn+1−axn+bxn−k=019941883719731MR130548010.1006/jmaa.1994.1457ZBL0842.39004KulenovićM. R. S.LadasG.MartinsL. F.RodriguesI. W.The dynamics of xn+1=(α+βxn)/(A+Bxn+Cxn−1): facts and conjectures2003456–910871099MR2000580ZBL1077.39004LadasG.On the recursive sequence xn+1=(α+βxn+γxn−1)/(A+Bxn+Cxn−1)199513317321MR135044710.1080/10236199508808030ZBL0855.39006CamouzisE.LadasG.VoulovH. D.On the dynamics of xn+1=(α+γxn−1+δxn−2)/(A+xn−2)200398731738MR199290610.1080/1023619021000042153ZBL1050.39003StevićS.On the recursive sequence xn+1=(−1/xn)+(A/xn−1)200127116MR188205410.1155/S0161171201010614ZBL1005.39016StevićS.The recursive sequence xn+1=g(xn,xn−1)/(A+xn)2002153305308MR1891551ZBL1029.3900710.1016/S0893-9659(01)00135-5