Next Article in Journal
On the Bounds for a Two-Dimensional Birth-Death Process with Catastrophes
Next Article in Special Issue
An Implicit Hybrid Method for Solving Fractional Bagley-Torvik Boundary Value Problem
Previous Article in Journal
On Some Sufficiency-Type Stability and Linear State-Feedback Stabilization Conditions for a Class of Multirate Discrete-Time Systems
Previous Article in Special Issue
A Numerical Method for Solving a Class of Nonlinear Second Order Fractional Volterra Integro-Differntial Type of Singularly Perturbed Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Behavior of Certain Nonautonomous Linearizable Three Term Difference Equations

by
E. J. Janowski
and
M. R. S. Kulenović
*
Department of Mathematics, University of Rhode Island, Kingston, RI 02881-0816, USA
*
Author to whom correspondence should be addressed.
Mathematics 2018, 6(5), 79; https://doi.org/10.3390/math6050079
Submission received: 10 April 2018 / Revised: 1 May 2018 / Accepted: 7 May 2018 / Published: 9 May 2018
(This article belongs to the Special Issue Advances in Differential and Difference Equations with Applications)

Abstract

:
We investigate the nonautonomous difference equation with real initial conditions and coefficients g i , i = 0 , 1 which are in general functions of n and/or the state variables x n , x n 1 , , and satisfy g 0 + g 1 = 1 . We also obtain some global results about the behavior of solutions of the nonautonomous non-homogeneous difference equation where g i , i = 0 , 1 , 2 are functions of n and/or the state variables x n , x n 1 , , with g 0 + g 1 = 1 . Our results are based on the explicit formulas for solutions. We illustrate our results by numerous examples.
JEL Classification:
AMS 2010 Mathematics Subject Classification: 39A10; 39A20; 39A30

1. Introduction

We investigate the global behavior of the second order difference equation
x n + 1 = g 0 x n + g 1 x n 1 , n = 0 , 1 , ,
where the functions g i : R k R , i = 0 , 1 ; k = 1 , 2 , , are such that
g 0 + g 1 = 1 , n = 0 , 1 , .
Here, the functions g i , i = 0 , 1 are in general functions of n and/or the state variables x n , x n 1 , . So Equation (1) is indeed a nonautonomous equation. It can also be an equation of order higher than two. As we have shown in [1] if | g 0 | + | g 1 | C < 1 , then the zero equilibrium of Equation (1) is globally asymptotically stable, and if | g 0 | + | g 1 | 1 , then the zero equilbrium of Equation (1) is stable. Furthermore, as we have shown in [1], if g 0 + g 1 r > 1 , where r is a constant independent of n and the state variables, then Equation (1) has unbounded solutions. Thus, the most interesting case is when condition (2) is satisfied. This case was addressed in [1] for the general k-th order difference equation, and some attractivity results of the equilibrium solutions were given there. However, a restriction that g i 0 , i = 0 , 1 was imposed. Similar results were obtained in [2] where attractivity results of the periodic solutions were obtained.
When Equation (1) is a linear equation with constant coefficients, condition (2) implies that we have the non-hyperbolic case since λ = 1 would be a characteristic root. In general, difference equations can be very sensitive to any kind of perturbations in the non-hyperbolic case as we can see from the following two examples:
Example 1.
Consider the difference equation
x n + 1 = 1 n x n + n 1 n x n 1 , n = 1 , 2 , .
Starting with the two initial conditions x 0 , x 1 we obtain the corresponding solution to be x n = x 1 , n 1 . The limiting case of Equation (3), as n , is the equation
x n + 1 = x n 1 , n = 1 , 2 , .
where every non-constant solution is the period-two solution { x 0 , x 1 } .
Example 2.
Consider the difference equation
x n + 1 = 1 2 + 1 2 n x n + 1 2 1 2 n x n 1 , n = 1 , 2 , .
Starting with two initial conditions x 0 , x 1 we obtain the corresponding solution to be x n = x 1 , n 1 . The limiting case of Equation (4), as n , is the equation
x n + 1 = 1 2 x n + 1 2 x n 1 , n = 1 , 2 , .
where every non-constant solution is given as:
x n = 1 3 x 0 + 2 x 1 + 2 ( x 0 x 1 ) 1 2 n , n = 0 , 1 ,
Thus the limiting equation has all solutions asymptotic to x 0 + 2 x 1 3 .
Example 3.
The following two difference equations
x n + 1 = x n 2 a + x n + a x n 1 a + x n , n = 1 , 2 ,
and
x n + 1 = b x n b + x n 1 + x n 1 2 b + x n 1 , n = 1 , 2 ,
where a , b > 0 , x 0 , x 1 0 , which have respective linearizations
x n + 1 = x n a + x n x n + a a + x n x n 1 , n = 1 , 2 ,
and
x n + 1 = b a + x n 1 x n + x n 1 b + x n 1 x n 1 , n = 1 , 2 , ,
with the condition (2) satisfied, possible models in population dynamics as all transition functions are Holling type functions, see [3] for further references.
There are many ways that a general second order difference equation of the form
x n + 1 = f ( n , x n , x n 1 ) , n = 0 , 1 ,
where f is a real function from R 3 to R , can be transformed into the linearized form (1). See [1,2,4] for linearization methods based on algebraic manipulations and especially [5] for symmetries based methods.
The linear Equation (1) can also be used to obtain generalized identities which can be used to study the behavior of the solutions of the difference equation, see [1,2]. Identities are the building blocks of semicycle analysis, method of invariants and rate of convergence. See [4,6,7,8,9,10,11,12,13,14].
We also extend our investigation to the second order nonautonomous non-homogeneous difference equation
x n + 1 = g 0 x n + g 1 x n 1 + g 2 , n = 0 , 1 , ,
where g 0 , g 1 and g 2 are real valued functions which depend on n and satisfy condition (2).
The technique which will be used to obtain our results is based on the explicit formulas for the solutions of Equations (1) and (7). These formulas are of independent interest and can be used in a variety of situations.
We illustrate our results with several examples.

2. Invariants

Our first result gives a connection between Equation (1) under the condition (2) and the existence of an invariant.
Lemma 1.
Consider Equation (1) under the condition (2) where g 0 , g 1 : R k + 1 R . Then Equation (1) has an invariant of the form I n = x n + g 1 x n 1 if and only if g 1 is a constant.
Proof. 
First, assume condition (2), where g 1 = A R for n 0 , is satisfied. Hence Equation (1) becomes
x n + 1 = ( 1 A ) x n + A x n 1 , n 0 .
Thus
x n + 1 + A x n = x n + A x n 1 , n 0
and so
I 0 = x 0 + A x 1 = x n + A x n 1 , n 0
is an invariant for Equation (1).
Second, assume that Equation (1) has an invariant of the form I n = x n + g 1 x n 1 . Then the following relation must be satisfied
x n + 1 + g 1 x n = x n + g 1 x n 1 , n 0 .
If g 1 is not a constant for all n 0 , then g 1 = A n . Thus
x n + 1 + A n + 1 x n = I n + 1 = I n = x n + A n x n 1 , n 0 .
Since by condition (2) g 0 = 1 A n for all n 0 then Equation (1) becomes
x n + 1 = ( 1 A n ) x n + A n x n 1 , n 0 ,
which by substitution in Equation (9) gives
( 1 A n ) x n + A n x n 1 + A n + 1 x n = x n + A n x n 1
and so x n ( A n + 1 A n ) = 0 for all n 0 .
Note that the invariant must hold for every solution of Equation (1) including the infinite number of constant solutions x n = x 0 , n 0 . Thus A n + 1 = A n for all n 0 , so g 1 is a constant for n 0 . ☐
Remark 1.
The restriction that g 0 and g 1 are constants is a severe restriction, yet the obtained invariant can be useful. This invariant will be used to reduce the order of difference Equation (1) resulting in first order linear difference equation which is solvable. See Proposition 1 and Corollary 2.
Lemma 1 can be extended to any three term difference equation of the form
x n + l = g i x n i + g k x n k , n = 0 , 1 , ,
where l { 1 , 2 , } and i , k Z such that l < i < k with
g i + g k = 1 , n = 0 , 1 , ,
Lemma 2.
Let l { 1 , 2 , } and i , k Z such that l < i < k . Consider Equation (10) where g i , g k : R k + l R and condition (11) is satisfied. Then Equation (10) has an invariant of the form
I n = x n + l 1 + + x n i + g k x n 1 i + + g k x n k , n = 0 , 1 ,
if and only if g k is a constant.
Proof. 
First assume condition (11) and that g k = A is a constant. Then g i = 1 A for n 0 and so Equation (10) becomes
x n + l = x n i A x n i + A x n k , n = 0 , 1 , .
Then Equation (12) is an invariant. Indeed
I n + 1 = x n + l + x n + l 1 + + x n + 1 i + A x n i + + A x n + 1 k = ( x n i A x n i + A x n k ) + x n + l 1 + + x n + 1 i + A x n i + + A x n + 1 k = x n + l 1 + + A n x n k , n = 0 , 1 , .
Second, assume that Equation (12) is an invariant of Equation (10). If g k is not a constant for all n 0 , then g k is a function of n that is, g k = A n , n = 0 , 1 , . Then Equation (10) with condition (11) becomes
x n + l = ( 1 A n ) x n i + A x n k , n = 0 , 1 ,
and I n + 1 = I n implies that for n 0
( x n i + x n 1 i + + x n + 1 k ) ( A n + 1 A n ) = 0 .
Therefore A n + 1 = A n for all n 0 , which completes the proof. ☐

3. Representation of Solutions

Next we give some basic results on the representation of the solutions of Equation (13). These results have importance of their own and have straight-forward extensions to the difference equations in linear normed spaces.
Lemma 3.
Consider the difference equation
x n + 1 = x n + F n ( x n x n 1 ) + D n , n = 0 , 1 , ,
where F n and D n are real valued functions for n = 0 , 1 , , . Then for n = 0 , 1 ,
(a) 
x n + 1 = x n + ( x 0 x 1 ) j = 0 n F j + k = 0 n D k i = k + 1 n F i ,
where by definition i = k + 1 n F i = 1 when k + 1 > n ;
(b) 
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + i = 0 n D i 1 + k = 1 n i j = i + 1 k + i F j ,
where by definition k = 1 n i j = i + 1 k + i F j = 0 when n < i + 1 .
Proof. 
Set u n = x n x n 1 , n = 0 , 1 , , in Equation (13). We obtain:
u n + 1 = F n u n + D n , n = 0 , 1 ,
which is a first order linear nonautonomous non-homogeneous difference equation whose solution is
u n = k = 0 n 1 F k u 0 + k = 0 n 1 D k j = k + 1 n 1 F j , n = 1 , 2 , .
Thus x n satisfies a first order linear nonautonomous non-homogeneous difference equation
x n = x n 1 + u n , n = 0 , 1 , ,
whose solution is given by Equation (14). By iterating Equation (14) we obtain Equation (15). ☐
From Lemma 3 we can obtain an explicit formula for the solution of the difference Equation (7) where g 0 , g 1 and g 2 are real valued functions which depend on n and satisfy condition (2). For convenience of notation set g 0 = A n , g 1 = B n and g 2 = D n at n = 0 , 1 , . Taking into account that A n + B n = 1 Equation (7) becomes
x n + 1 = x n B n ( x n x n 1 ) + D n , n = 0 , 1 , .
By Lemma 3 part (b) with F n = B n , n 0 we get that for n = 0 , 1 ,
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n ( 1 ) i + 1 j = 0 i B j + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i B j
and so for n = 0 , 1 ,
x n = 1 + j = 0 n 1 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 n 1 ( 1 ) j i = 0 j B i x 1 + i = 0 n 1 D i 1 + k = 1 n i 1 ( 1 ) k j = 1 + i k + i B j .

4. Main Results

In this section we present our main results. Consider Equation (7) where the function g 2 = D n , n = 0 , 1 , , in which case Equation (7) becomes
x n + 1 = g 0 x n + g 1 x n 1 + D n , n = 0 , 1 , .
Remark 2.
Consider Equation (18) where the functions g 0 and g 1 satisfies condition (2).
(1) 
If for some N = 0 , 1 , D N 0 , then Equation (18) does not have a constant solution. Indeed, assume that Equation (18) has a constant solution x ¯ . Then D N = x ¯ g 0 x ¯ g 1 x ¯ = ( 1 g 0 g 1 ) x ¯ = 0 , which is a contradiction.
(2) 
If x 0 = x 1 and D n = 0 for all n 0 , then for n 0 , x n + 1 = x 0 , which follows from Lemma 3 part (b).
Theorem 1.
Consider Equation (18) where the functions g 0 and g 1 satisfies condition (2).
Let r ( 0 , ) and assume that x n R . Assume that either x 0 x 1 or D n 0 for some n 0 .
(1) 
Suppose that g 1 1 for n 0 .
(a) 
If x 0 x 1 and D n 0 for n 0 , then lim n x n = ;
(b) 
If x 0 x 1 and D n 0 for n 0 , then lim n x n = ;
(2) 
Suppose that 0 < | g 1 | r < 1 for n 0 and i = 0 D i converges.
(a) 
If either D n = 0 for n 0 or D n 0 for n 0 , then every solution of Equation (18) converges to a limit point;
(b) 
If D n 0 for n 0 , then every solution of Equation (18) is bounded;
(3) 
Suppose that g 1 = 0 for all n 0 , and let a R .
(a) 
If D n = 0 for n 0 , then x n = x 0 ;
(b) 
If i = 0 D i converges, then every solution of Equation (18) converges to a limit point;
(c) 
If i = 0 D i = ± , then lim n x n = ± ;
(d) 
If D 2 i = a and D 2 i + 1 = a , i = 0 , 1 , , then every solution of Equation (18) is a period-two sequence of the form , x 0 + a , x 0 , .
(4) 
Suppose that g 1 = 0 for finite number of n 0
(a) 
If D n = 0 for n 0 , then every solution of Equation (18) is eventually equal to a constant;
(b) 
If D n 0 for some n 0 , then the solution of Equation (18) may have all kinds of dynamic scenarios which depend on the properties of the nonhomogeneous part.
(5) 
Suppose that g 1 = 1 for n 0 .
(a) 
If D n = 0 for n 0 , then every solution of Equation (18) is a period-two sequence of the form , x 1 , x 0 , .;
(b) 
If either i = 0 D 2 i = ± or i = 0 D 2 i + 1 = ± , then every solution of Equation (18) is unbounded;
(c) 
If i = 0 D 2 i and i = 0 D 2 i + 1 converge, then every solution of Equation (18) converges to (not necessarily minimal) period-two sequence;
(6) 
Suppose that 1 < r g 1 for n 0 . Assume that for n 0 either
x 0 x 1 , D 2 n 0 , D 2 n + 1 0
or
x 0 x 1 , D 2 n 0 , D 2 n + 1 0 .
Then every solution of Equation (18) is unbounded.
Proof. 
When we apply condition (2) to Equation (18) we get that
x n + 1 = x n g 1 ( x n x n 1 ) + D n , n = 0 , 1 ,
For convenience of notation set F n = g 1 and H n = g 1 at each n = 0 , 1 , . For n 0 let R n be the sum of positive terms and let S n be the sum of negative terms.
(1)
Suppose that g 1 1 for n 0 . Then for n 0 , F n 1 . By Lemma 3 part (b) we get that for n 0
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + i = 0 n D i 1 + k = 1 n i j = 1 + i k + i F j .
If D n 0 for some n 0 , then let N be the first n 0 such that D n 0 . Then for n 0
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + D N 1 + k = 1 n N j = 1 + N k + N F j + i = N + 1 n D i 1 + k = 1 n i j = 1 + i k + i F j .
Observe that for i = 0 , 1 , 1 j = 0 i F j and for k j , 1 j = 1 + i k + i F j . Thus for n 0 n < i = 0 n j = 0 i F j and for n N n N < 1 + k = 1 n N j = 1 + N k + N F j .
(a)
Assume that x 0 x 1 and D n 0 , n = 0 , 1 , . Then for n N
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + D N 1 + k = 1 n N j = 1 + N k + N F j + R n > x 0 + ( x 0 x 1 ) n + ( n N ) D N .
Thus x n as n .
(b)
Assume that x 0 x 1 and D n 0 , n = 0 , 1 , . Then for n N
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + D N 1 + k = 1 n N j = 1 + N k + N F j + S n < x 0 + ( x 0 x 1 ) n + ( n N ) D N .
Thus x n as n .
(2)
Suppose that 0 < | g 1 | r < 1 for n 0 and i = 0 D i converges. Then for n 0 0 < | H n | r < 1 . By Lemma 3 (b) with F n = | H n | we have that for n 0
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n ( 1 ) i + 1 j = 0 i | H j | + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | .
Observe that for i = 0 , 1 , , 0 < | ( 1 ) i + 1 j = 0 i H j | r i + 1 . Hence the series i = 0 n ( 1 ) i + 1 j = 0 i H j converges by the comparison test for series to the limiting value.
(a)
If D n = 0 for n 0 , then x n converges to a limit as n . Now assume that D n 0 for n 0 . Observe that for n 0 and for any i = 0 , 1 , , such that n i + 1 we have the following:
When n i is odd, then
0 1 | H i + 1 | + ( 1 | H i + 3 | ) | H i + 1 | | H i + 2 | + + ( 1 | H n | ) j = 1 + i n 1 | H j | = 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | < 1 + k = 1 n i j = 1 + i k + i | H j | < k = 0 r k = 1 1 r .
When n i is even, then
0 1 | H i + 1 | + ( 1 | H i + 3 | ) | H i + 1 | | H i + 2 | + + ( 1 | H n 1 | ) j = 1 + i n 2 | H j | + j = 1 + i n | H j | = 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | < 1 + k = 1 n i j = 1 + i k + i | H j | < k = 0 r k = 1 1 r .
Hence for n 0 and for any i = 0 , 1 , , such that n i + 1 and D i > 0
0 D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | < D i 1 r
and so for n 0
0 i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | < 1 1 r i = 0 n D i .
Thus the series
i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j |
converges by comparison test and so x n converges to a finite limit.
(b)
Assume that D n 0 for n 0 . Then there exists L ( , 0 ) such that i = 0 D i = L . From part (a) we have that for n i + 1 , i = 0 , 1 , ,
0 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | < 1 1 r .
Then for n i + 1 , i = 0 , 1 , , and D i < 0
0 D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | > D i 1 r
and so
0 i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | > 1 1 r i = 0 n D i > 1 1 r i = 0 D i = L 1 r .
If x 0 = x 1 , then
x 0 x 0 + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | > x 0 + L 1 r .
Thus x 0 x n + 1 > x 0 + L 1 r for n 0 .
Now assume that x 0 x 1 . From part (a) (with i = 0 ) we have that for n 1
0 1 + k = 1 n ( 1 ) k j = 1 k | H j | < 1 1 r .
If | H 0 | ( x 0 x 1 ) < 0 , then from Equation (19) we get that for n 0
x 0 x 0 + ( x 0 x 1 ) k = 0 n ( 1 ) k + 1 j = 0 k | H j | x 0 | H 0 | ( x 0 x 1 ) 1 r
and so for n 0
x 0 x 0 + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | x 0 + ( x 0 x 1 ) i = 0 n ( 1 ) i + 1 j = 0 i | H j | + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | x 0 | H 0 | ( x 0 x 1 ) 1 r + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | x 0 | H 0 | ( x 0 x 1 ) 1 r + L 1 r .
Thus x 0 x n + 1 x 0 + L | H 0 | ( x 0 x 1 ) 1 r for n 0 .
If | H 0 | ( x 0 x 1 ) < 0 , then from Equation (19) we get that for n 0
x 0 x 0 + ( x 0 x 1 ) i = 0 n ( 1 ) i + 1 j = 0 i | H j | x 0 | H 0 | ( x 0 x 1 ) 1 r
and so
x 0 + L 1 r x 0 + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | x 0 + ( x 0 x 1 ) i = 0 n ( 1 ) i + 1 j = 0 i | H j | + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | x 0 | H 0 | ( x 0 x 1 ) 1 r + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i | H j | x 0 | H 0 | ( x 0 x 1 ) 1 r .
Thus for n 0
x 0 + L 1 r x n + 1 x 0 | H 0 | ( x 0 x 1 ) 1 r
and so every solution is bounded.
(3)
Suppose that g 1 = 0 for all n 0 . Then Equation (18) becomes x n + 1 = x n + D n , n = 0 , 1 , , and so x n + 1 = x 0 + i = 0 n D i , n = 0 , 1 , , from which (a), (b) and (c) follow.
(d)
Let a R \ { 0 } and assume that for i = 0 , 1 , , D 2 i = a and D 2 i + 1 = a . Equation (18) gives for n = 0 , 1 , ,
x 2 n + 1 = x 0 + i = 0 2 n D i = x 0 + a , x 2 n + 2 = x 0 + i = 0 2 n + 1 D i = x 0 .
(4)
Suppose that g 1 = 0 for finite number of n 0 . By Lemma 3 part (b) with F n = H n we get that
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n ( 1 ) i + 1 j = 0 i H j + i = 0 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i H j .
Let N 0 be the first n such that g 1 = 0 and let M = 0 , 1 , , be the last n 0 such that g 1 = 0 . Then H N = H M = 0 .
(a)
Assume that D n = 0 for all n 0 . Then by Equation (20) for n N 1
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 N 1 ( 1 ) i + 1 j = 0 i H j .
(b)
Assume that D n 0 for some n 0 . Then by Equation (20) for n M
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 N 1 ( 1 ) i + 1 j = 0 i H j + i = 0 M 1 D i 1 + k = 1 M i 1 ( 1 ) k j = 1 + i k + i H j + i = M n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i H j = c o n s t a n t + i = M n D i + i = M n D i k = 1 n i ( 1 ) k j = 1 + i k + i H j .
(5)
Suppose that g 1 = 1 for all n 0 . Then Equation (18) becomes x n + 1 = x n 1 + D n , n = 0 , 1 , . Thus for n 0
x 2 n + 1 = x 1 + i = 0 n D 2 i , x 2 n + 2 = x 0 + i = 0 n D 2 i + 1 .
(a)
If D n = 0 for all n 0 , then x 2 n + 1 = x 1 and x 2 n + 2 = x 0 for n = 0 , 1 , .
(b)
If i = 0 D 2 i = ± , then x 2 n + 1 ± as n . If i = 0 D 2 i + 1 = ± , then x 2 n + 2 ± as n .
(c)
If i = 0 D 2 i converges and i = 0 D 2 i + 1 converges, then x 2 n + 1 and x 2 n + 2 converge as n .
(6)
Suppose that 1 < r g 1 for all n 0 . Then 1 < r H n for all n 0 . From Lemma 3 part (b) with F n = H n we get that for n 0
x 2 n + 2 = x 0 + ( x 0 x 1 ) i = 0 2 n + 1 ( 1 ) i + 1 j = 0 i H j + i = 0 2 n + 1 D i 1 + k = 1 2 n + 1 i ( 1 ) k j = 1 + i k + i H j = x 0 + ( x 0 x 1 ) i = 0 n ( H 2 i + 1 1 ) j = 0 2 i H j i = 0 n D 2 i k = 0 n i ( H 2 k + 2 i + 1 1 ) j = 1 + 2 i 2 k + 2 i H j + i = 0 n D 2 i + 1 1 + k = 1 n i ( H 2 k + 2 i + 1 1 ) j = 2 + 2 i 2 k + 2 i H j .
Observe that for n 0 , ( r 1 ) r 2 n < ( H 2 n + 1 1 ) j = 0 2 n H j and lim n r 2 n N = for any fixed N. Let N be the first n 0 such that D n 0 . Suppose that N is even. Then for n N / 2 we have that
x 2 n + 2 = x 0 + ( x 0 x 1 ) i = 0 n 1 ( H 2 i + 1 1 ) j = 0 2 i H j + ( x 0 x 1 ) ( H 2 n + 1 1 ) j = 0 2 n H j D N k = 0 n N / 2 ( H 2 k + N + 1 1 ) j = 1 2 k H j + N i = N + 2 2 n D 2 i k = 0 n i ( H 2 k + 2 i + 1 1 ) j = 1 + 2 i 2 k + 2 i H j + i = N / 2 n D 2 i + 1 1 + k = 1 n i ( H 2 k + 2 i + 1 1 ) j = 2 + 2 i 2 k + 2 i H j .
For n > N / 2
( r 1 ) r 2 n N < ( r 1 ) + ( r 1 ) r 2 + + ( r 1 ) r 2 n N ( H N + 1 1 ) + ( H N + 3 1 ) j = 1 2 H j + N + + ( H 2 n + 1 1 ) j = 1 2 n N H j + N = k = 0 n N / 2 ( H 2 k + N + 1 1 ) j = 1 2 k H j + N .
(a)
Assume that x 0 x 1 , D 2 n 0 , D 2 n + 1 0 , n = 0 , 1 , . Then for n N / 2
x 2 n + 2 = x 0 + ( x 0 x 1 ) ( H 2 n + 1 1 ) j = 0 2 n H j D N k = 0 n N / 2 ( H 2 k + N + 1 1 ) j = 0 2 k H j + N + R n > x 0 + ( x 0 x 1 ) ( r 1 ) r 2 n D N ( r 1 ) r 2 n N .
Thus x 2 n + 2 as n .
(b)
Assume that x 0 x 1 , D 2 n 0 , D 2 n + 1 0 , n = 0 , 1 , . Then for n > N / 2
x 2 n + 2 = x 0 + ( x 0 x 1 ) ( H 2 n + 1 1 ) j = 0 2 n H j D N k = 0 n N / 2 ( H 2 k + N + 1 1 ) j = 0 2 k H j + N + S n < x 0 + ( x 0 x 1 ) ( r 1 ) r 2 n D N ( r 1 ) r 2 n N .
Thus x 2 n + 2 as n .
Now suppose that N is odd. Then for n N we have that
x 2 n + 2 = x 0 + ( x 0 x 1 ) i = 0 n 1 ( H 2 i + 1 1 ) j = 0 2 i H j + ( x 0 x 1 ) ( H 2 n + 1 1 ) j = 0 2 n H j + D N 1 + k = 1 n N 1 2 ( H 2 k + N 1 ) j = 1 2 k 1 H j + N + i = N + 1 2 n D 2 i + 1 1 + k = 1 n i ( H 2 k + 2 i + 1 1 ) j = 2 + 2 i 2 k + 2 i H j i = N + 1 2 n D 2 i k = 0 n i ( H 2 k + 2 i + 1 1 ) j = 1 + 2 i 2 k + 2 i H j .
For n N
( r 1 ) r 2 n N < 1 + ( r 1 ) + ( r 1 ) r 2 + + ( r 1 ) r 2 n N 1 + ( H N + 2 1 ) j = 1 1 H j + N + + ( H 2 n + 1 1 ) j = 1 2 n N H j + N = 1 + k = 1 n N 1 2 ( H 2 k + N 1 ) j = 1 2 k 1 H j + N .
(a)
Assume that x 0 x 1 , D 2 n 0 , D 2 n + 1 0 , n = 0 , 1 , . Then for n N
x 2 n + 2 = x 0 + ( x 0 x 1 ) ( H 2 n + 1 1 ) j = 0 2 n H j + D N 1 + k = 1 n N 1 2 ( H 2 k + N 1 ) j = 0 2 k 1 H j + N + R n > x 0 + ( x 0 x 1 ) ( r 1 ) r 2 n + D N ( r 1 ) r 2 n N .
Thus x 2 n + 2 as n .
(b)
Assume that x 0 x 1 , D 2 n 0 , D 2 n + 1 0 , n = 0 , 1 , . Then for n N
x 2 n + 2 = x 0 + ( x 0 x 1 ) ( H 2 n + 1 1 ) j = 0 2 n H j + D N 1 + k = 1 n N 1 2 ( H 2 k + N 1 ) j = 1 2 k 1 H j + N + S n < x 0 + ( x 0 x 1 ) ( r 1 ) j = 0 2 n H j + D N ( r 1 ) r 2 n N .
Thus x 2 n + 2 as n . ☐
Example 4.
The difference equation
x n + 1 = x n e x n + ( a + 1 ) x n x n 1 a + e x n , n = 0 , 1 , ,
where a 0 , when rewritten as
x n + 1 = e x n + ( a + 1 ) a + e x n x n + 1 a + e x n x n 1 = g 0 x n + g 1 x n 1 , n = 0 , 1 , ,
satisfies all conditions of Theorem 1 part (2a) when a > 1 , that is g 0 + g 1 = 1 and | g 1 | = 1 a + e x n 1 a and so every solution of Equation (21), when a > 1 converges to a limit point.
In addition, Equation (21) implies
x n + 1 x n = x n x n 1 a + e x n , n = 0 , 1 , ,
which shows that when x 0 > x 1 (resp. x 0 < x 1 ) the solution { x n } is increasing (resp. decreasing) and so it has to have a finite or infinite limit.
Suppose that 0 a 1 . Assume that x 0 < x 1 and x 0 < 0 . Then e x 0 > 1 and so there exists c ( 0 , ) such that e x n e x 0 c > 1 for n 0 . Thus for n 0
0 > 1 a + e x n 1 a + c > 1 .
Therefore by Theorem 1 part (2a) lim n x n = L ( , ) .
Now suppose that a [ 0 , 1 ) . Assume that x 0 > x 1 such that x 0 ln ( 1 a ) . Then e x n e x 0 1 a and so 1 a + e x n 1 for n 0 . Thus by Theorem 1 part (1a) lim n x n = .
In some very special cases of the terms D n in Equation (18) one can obtain results not included in Theorem 1. We illustrate this with the following result.
Theorem 2.
Consider the difference equation
x n + 1 = A n x n + B n x n 1 + D n , n = 0 , 1 ,
where A n , B n R \ { 0 } such that A n + B n = 1 for all n 0 . Let r ( 0 , ) . Then every solution of Equation (24) is unbounded provided one of the following holds:
(1) 
Assume that 1 < r B n < 0 and either i = 0 D i = , where D n 0 or i = 0 D i = , where D n 0 for n 0 ;
(2) 
Assume that 0 < B n r < 1 for n 0 .
(a) 
Suppose that for all n 0 either
D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i = or D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i = ;
(b) 
Suppose that for all n 0 either
D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i + 1 = or D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i + 1 = .
Proof. 
First we rewrite Equation (24) in the form
x n + 1 = x n B n ( x n x n 1 ) + D n , n = 0 , 1 , .
Let R n (resp. S n ) be the sum of positive (resp. negative) terms.
(1)
Assume that 1 < r B n < 0 for n 0 . Setting F n = B n Lemma 3 part (b) implies for n 0
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + i = 0 n D i k = 0 n i j = 1 + i k + i F j + i = 0 n D i .
Since 0 < j = 0 i F j r i for i = 0 , 1 , , then the series i = 0 n j = 0 i F j converges as n .
If for n 0 , D n 0 and i = 0 D i = , then
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + i = 0 n D i + R n .
Thus x n as n .
If for n 0 , D n 0 and i = 0 D i = , then
x n + 1 = x 0 + ( x 0 x 1 ) i = 0 n j = 0 i F j + i = 0 n D i + S n .
Thus x n as n .
(2)
Assume that 0 < B n r < 1 for n 0 .
(a)
By Lemma 3 part (b) with F n = B n we get that for n 0
x 2 n + 1 = x 0 + ( x 0 x 1 ) i = 0 2 n ( 1 ) i + 1 j = 0 i B j + i = 0 2 n D i 1 + k = 1 n i ( 1 ) k j = 1 + i k + i B j = x 0 + ( x 0 x 1 ) i = 0 2 n ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i 1 k = 1 n i ( 1 B 2 k + 2 i ) j = 1 + 2 i 2 k + 2 i 1 B j + i = 0 n D 2 i + 1 k = 1 n i ( 1 B 2 k + 2 i ) j = 2 + 2 i 2 k + 2 i 1 B j = x 0 + ( x 0 x 1 ) i = 0 2 n ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i + i = 0 n ( D 2 i + 1 B 2 i + 1 D 2 i ) k = 1 n i ( 1 B 2 k + 2 i ) j = 2 + 2 i 2 k + 2 i 1 B j .
Observe that the series i = 0 ( 1 ) i + 1 j = 0 i B j converges by the alternating series test.
If D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i = , then for n 0
x 2 n + 1 = x 0 + ( x 0 x 1 ) i = 0 2 n ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i + R n .
Thus x 2 n + 1 as n .
If D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i = , then for n 0
x 2 n + 1 = x 0 + ( x 0 x 1 ) i = 0 2 n ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i + S n .
Thus x 2 n + 1 as n .
(b)
By Lemma 3 part (b) with F n = B n we get that for n 0
x 2 n + 2 = x 0 + ( x 0 x 1 ) i = 0 2 n + 1 ( 1 ) i + 1 j = 0 i B j + i = 0 2 n + 1 D i 1 + k = 1 2 n + 1 i ( 1 ) k j = 1 + i k + i B j = x 0 + ( x 0 x 1 ) i = 0 2 n + 1 ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i + 1 1 k = 1 n i ( 1 B 2 k + 2 i + 1 ) j = 2 + 2 i 2 k + 2 i B j + i = 0 n D 2 i 1 B 2 i + 1 + k = 1 n i ( 1 B 2 k + 2 i + 1 ) j = 1 + 2 i 2 k + 2 i B j = x 0 + ( x 0 x 1 ) i = 0 2 n + 1 ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i ( 1 B 2 i + 1 ) + i = 0 n D 2 i + 1 + i = 0 n ( B 2 i + 1 D 2 i D 2 i + 1 ) k = 1 n i ( 1 B 2 k + 2 i + 1 ) j = 2 + 2 i 2 k + 2 i B j .
If D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i + 1 = , then for n 0
x 2 n + 2 = x 0 + ( x 0 x 1 ) i = 0 2 n + 1 ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i + 1 + R n .
Thus x 2 n + 2 as n .
If D 2 n + 1 B 2 n + 1 D 2 n , D n 0 and i = 0 D 2 i + 1 = , then for n 0
x 2 n + 2 = x 0 + ( x 0 x 1 ) i = 0 2 n + 1 ( 1 ) i + 1 j = 0 i B j + i = 0 n D 2 i + 1 + S n .
Thus x 2 n + 2 as n . ☐
Remark 3.
In Theorem 1 cases (2) and (6) the existence of r is a sufficient condition but not a necessary condition. To show this for case (2) consider Equation (3) where n = 2 , 3 , . Since g 1 = n 1 n 1 as n , there is no r such that g 1 r < 1 . Rewriting Equation (3) as
x n + 1 = x n n 1 n ( x n x n 1 ) , n = 2 , 3 , ,
we find its solution as
x n + 1 = x 2 ( x 2 x 1 ) i = 2 n ( 1 ) i 1 i , n = 2 , 3 , ,
which converges to a constant.
Now for the case (5), consider the difference equation
x n + 1 = 1 n x n + n + 1 n x n 1 , n = 1 , 2 , ,
Clearly g 0 + g 1 = 1 , g 1 1 as n and there is no r such that g 1 r > 1 .
Rewriting Equation (25) as
x n + 1 = x n n + 1 n ( x n x n 1 ) , n = 1 , 2 , ,
we find its solution as
x n + 1 = x 1 + ( x 1 x 0 ) i = 1 n ( 1 ) i ( i + 1 ) ,
which is unbounded.
Example 5.
Consider the following equation
x n + 1 = 1 2 + a n x n + 1 2 a n x n 1 , n = 0 , 1 , ,
where a n is a real sequence, which depends on n and on one or more state variables x n . An immediate application of Theorem 1 with D n = 0 , n = 0 , 1 , , leads to the following result for Equation (26). Here r is a positive constant.
Corollary 1.
The global behavior of Equation (26) is as follows:
(a) 
If 3 2 a n for n 0 , then every non-constant solution of Equation (26) is unbounded.
(b) 
If 0 < | 1 2 a n | r < 1 for n 0 , then every non-constant solution of Equation (26) converges to a limit.
(c) 
If a n = 1 2 for all n 0 , then x n = x 0 .
(d) 
If a n = 1 2 for a finite number of n 0 , then every solution of Equation (26) is eventually equal to a constant.
(e) 
If a n = 1 2 for n 0 , then every non-constant solution of Equation (26) is a period-two solution.
(f) 
If a n 1 2 r < 1 2 , r > 1 for n 0 , then every non-constant solution of Equation (26) is unbounded.
In the case when Equation (1) has periodic coefficients it can be reduced to an associated system of second order equations with constant coefficients which can be solved. We illustrate this case by the Example 6.
Example 6.
Consider the following second order equation
x n + 1 = A n x n + B n x n 1 , n = 0 , 1 , ,
where B n R \ { 0 } and A n + B n = 1 for all n 0 . Assume that A n and B n are period-two sequences, that is
A n = A i , B n = B i , n = 2 j + i , i = 0 , 1 ; j = 0 , 1 , .
A straight-forward simplification leads to the following autonomous system for even and odd indexed terms:
x 2 n + 3 = ( 1 + B 0 B 1 ) x 2 n + 1 B 0 B 1 x 2 n 1 x 2 n + 4 = ( 1 + B 0 B 1 ) x 2 n + 2 B 0 B 1 x 2 n , n = 0 , 1 , .
Let y n α = x 2 n + 1 + α , α = 0 , 1 , n = 0 , 1 , . Then Equation (27) is equivalent with the autonomous equation
y n + 1 α = ( 1 + B 0 B 1 ) y n α B 0 B 1 y n 1 α , α = 0 , 1 , n = 0 , 1 , .
In view of Lemma 1 and Theorem 1, we have the following result.
Proposition 1.
Consider Equation (27) with period-two coefficients where B n R \ { 0 } and A n + B n = 1 for all n 0 . Then
(1) 
If | B 0 B 1 | < 1 , then every nonconstant solution { x n } of Equation (27) converges to the limit
lim n x n = ( 1 B 0 ) x 0 + B 0 ( 1 B 1 ) x 1 1 B 0 B 1 ;
(2) 
If B 0 B 1 = 1 , then every nonconstant solution { x n } of Equation (27) is of period 4;
(3) 
If either B 0 B 1 1 or B 0 B 1 < 1 , then every nonconstant solution { x n } of Equation (27) is unbounded.
Proof. 
By Theorem 1 with D n = 0 , n = 0 , 1 , , we have
(1)
If | B 0 B 1 | < 1 , then every nonconstant solution { y n α } of Equation (28) converges to the limit L α , that is
lim n y n α = L α , α = 0 , 1 .
(2)
If B 0 B 1 = 1 , then every nonconstant solution { y n α } of Equation (28) is of period 2;
(3)
If either B 0 B 1 1 or B 0 B 1 < 1 , then every nonconstant solution { y n α } of Equation (28) is unbounded.
Now the results follow from (1)–(3). In particular, by Lemma 1 Equation (28) has the invariant I n α = y n α B 0 B 1 y n 1 α , where α = 0 , 1 and n = 0 , 1 , . If | B 0 B 1 | < 1 , then every nonconstant solution { y n α } of Equation (28) converges to the limit L α , which implies
lim n I n α = L α ( 1 B 0 B 1 ) = I 0 α = y 0 α B 0 B 1 y 1 α .
Now by using Equation (17), Lemma 1 and Theorem 1 we can obtain explicitly the behavior of the solutions of Equation (1) when g 0 and g 1 are periodic and condition (2) is satisfied.
Corollary 2.
Let p { 1 , 2 , } and let B n R \ { 0 } , n 0 . Consider Equation (27) where A n + B n = 1 for all n 0 and that
A n = A k , B n = B k , k = 0 , 1 , , p 1 , n = p j + k ; j = 0 , 1 , .
(a) 
If either ( 1 ) p 1 i = 0 p 1 B i 1 or ( 1 ) p 1 i = 0 p 1 B i > 1 , then every non-constant solution { x n } of Equation (27) is unbounded.
(b) 
If ( 1 ) p 1 i = 0 p 1 B i = 1 , then every non-constant solution { x n } of Equation (27) is periodic of period 2 p . In particular, the periodic solution is { x 1 , x γ } where
x γ = x 0 + ( x 0 x 1 ) i = 0 γ 1 ( 1 ) i + 1 j = 0 i B j ,
and γ = 0 , 1 , , 2 p 2 .
(c) 
If | i = 0 p 1 B i | < 1 , then every non-constant solution { x n } of Equation (27) converges to a limit, that is
lim n x n = 1 + j = 0 p 2 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 p 1 ( 1 ) j i = 0 j B i x 1 1 + ( 1 ) p 1 i = 0 p 1 B i .
Proof. 
Equation (27) implies that for β = 0 , 1 ,
x p n + β + 1 = ( 1 B β ) x p n + β + B β x p n + β 1 , n = 0 , 1 ,
and so
( 1 B β ) x p n + β = x p n + β + 1 B β x p n + β 1 .
Then for β = 0 , 1 , ; n 0
1 + j = 0 β ( 1 ) j + 1 i = 0 j B i x p n + β = 1 + j = 0 β 1 ( 1 ) j + 1 i = 0 j B i x p n + β + 1 + ( 1 ) β + 1 i = 0 β B i x p n 1 .
Indeed for β = 0 , 1 , , we have that for n 0
( 1 B β + 1 ) ( 1 + j = 0 β ( 1 ) j + 1 i = 0 j B i ) x p n + β + 1 = 1 + j = 0 β ( 1 ) j + 1 i = 0 j B i x p n + β + 2
1 + j = 0 β ( 1 ) j + 1 i = 0 j B i B p + 1 x p n + β .
Multiplying Equation (30) by B β + 1 and substituting into Equation (31) we obtain
1 + j = 0 β + 1 ( 1 ) j + 1 i = 0 j B i x p n + β + 1 = 1 + j = 0 β ( 1 ) j + 1 i = 0 j B i x p n + β + 2 + ( 1 ) β + 2 i = 0 β + 1 B i x p n 1 .
Now from Equation (29) when β = p 1 we have that for n 0
x p n + p = ( 1 B p 1 ) x p n + p 1 + B p 1 x p n + p 2 .
By using Equation (29) with β = p and Equation (32) we obtain
x p n + p + 1 = ( 1 B p 1 + B p 1 B p ) x p n + p 1 + B p 1 ( 1 B p ) x p n + p 2 .
Now continuing this process we get that for δ = 0 , 1 ,
x p n + p + δ = 1 + j = 0 δ ( 1 ) j + 1 i = 0 j B i + p 1 x p n + p 1 + B p 1 1 + j = 0 δ 1 ( 1 ) j + 1 i = 0 j B i + p x p n + p 2 .
Observe that since B i = B i + p for i = 0 , 1 , , then
j = 0 p 1 ( 1 ) j + 1 i = 0 j B i + p 1 = j = 0 p 2 ( 1 ) j + 1 i = 0 j B i + p 1 + ( 1 ) p i = 0 p 1 B i + p 1 = B p 1 1 + j = 0 p 3 ( 1 ) j + 1 i = 0 j B i + ( 1 ) p i = 0 p 1 B i .
Now multiplying Equation (30) where β = p 2 with B p 1 and substituting into Equation (33) where δ = p 1 we obtain
x p n + 2 p 1 = 1 + ( 1 ) p i = 0 p 1 B i x p n + p 1 + ( 1 ) p 1 i = 0 p 1 B i x p n 1 .
Similarly, we can show that for α = 1 , 2 , , p 1
x p n + 2 p 1 + α = 1 + ( 1 ) p i = 0 p 1 B i x p n + p 1 + α + ( 1 ) p 1 i = 0 p 1 B i x p n 1 + α .
For α = 0 , , p 1 let y n α = x p n + p 1 + α , n = 0 , 1 , . Then for α = 0 , , p 1 , y n + 1 α = x p n + 2 p 1 + α and y n 1 α = x p n 1 + α . Thus for α = 0 , , p 1 Equation (34) becomes
y n + 1 α = 1 + ( 1 ) p i = 0 p 1 B i y n α + ( 1 ) p 1 i = 0 p 1 B i y n 1 α .
(a)
Suppose that either ( 1 ) p 1 i = 0 p 1 B i 1 or ( 1 ) p 1 i = 0 p 1 B i > 1 . By Theorem 1 parts (1) and (6), where g 1 = ( 1 ) p 1 i = 0 p 1 B i , and D n = 0 , n = 0 , 1 , , every non-constant solution { y n α } , α = 0 , 1 , , p 1 of Equation (35) is unbounded. Thus every non-constant subsequence { x p n + p 1 + α } of Equation (27) is unbounded.
(b)
Suppose that ( 1 ) p 1 i = 0 p 1 B i = 1 . Then for α = 0 , 1 , , p 1 Equation (35) becomes y n + 1 α = y n 1 α , n = 0 , 1 , . Thus for α = 0 , 1 , , p 1 the solution { y n α } is a period-two sequence. Since there are p subsequences, then every non-constant solution { x n } of Equation (27) is periodic of period 2 p .
(c)
Suppose that | i = 0 p 1 B i | < 1 . Then by Theorem 1 part (2), where g 1 = ( 1 ) p 1 i = 0 p 1 B i and D n = 0 , n = 0 , 1 , , every non-constant solution { y n α } , α = 0 , 1 , , p 1 of Equation (35) converges to a limit point. Let L α = lim n y n α for α = 0 , , p 1 . Then by Lemma 1
L α 1 + ( 1 ) p 1 i = 0 p 1 B i = y 0 α + ( ( 1 ) p 1 i = 0 p 1 B i ) y 1 α .
Thus for α = 0 , , p 1
lim n y n α = y 0 α + ( 1 ) p 1 i = 0 p 1 B i y 1 α 1 + ( 1 ) p 1 i = 0 p 1 B i
and so for α = 0 , , p 1
lim n x p n + p 1 + α = x p 1 + α + ( 1 ) p 1 i = 0 p 1 B i x 1 + α 1 + ( 1 ) p 1 i = 0 p 1 B i .
In view of Equation (17) with D n = 0 , n = 0 , 1 , , we have that for α = 2 , 3 , , p 1
x p 1 + α = 1 + j = 0 p 2 + α ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 p 2 + α ( 1 ) j i = 0 j B i x 1
and
x α 1 = 1 + j = 0 α 2 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 α 2 ( 1 ) j i = 0 j B i x 1 .
Then for α = 2 , 3 , , p 1 and after simplifying terms we obtain
x p 1 + α + ( 1 ) p 1 i = 0 p 1 B i x α 1 = 1 + j = 0 p 2 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 p 1 ( 1 ) j i = 0 j B i x 1 .
It should be noticed that by using Equation (37) and the initial conditions we also get Equation (38) for the cases α = 0 , 1 . Thus for α = 0 , 1 , , p 1
lim n x p n + p 1 + α = 1 + j = 0 p 2 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 p 1 ( 1 ) j i = 0 j B i x 1 1 + ( 1 ) p 1 i = 0 p 1 B i .
Since all subsequences have the same limit, we finally obtain the formula for the limit. ☐
Example 7.
Consider the difference equation
x n + 1 = ( 1 y n + 1 ) x n + y n + 1 x n 1 , n = 0 , 1 ,
where
y n + 1 = 1 + y n y n 1 , n = 0 , 1 , .
Assume that y 1 , y 0 { 1 , 0 , 1 ± 5 2 } and y 1 + y 0 + 1 0 . Then the solution { y n } of Equation (40) is periodic of period 5, that is the solution is
, y 1 , y 0 , 1 + y 0 y 1 , 1 + y 0 + y 1 y 0 y 1 , 1 + y 1 y 0 , y 1 , y 0 ,
Equation (40), see [12] also has the invariant
I n = 1 + 1 y n 1 1 + 1 y n 1 + y n + y n 1 , n = 0 , 1 ,
and so for n 0
I n = I 0 = 1 + 1 y 1 1 + 1 y 0 1 + y 0 + y 1 = y 1 y 2 y 3 y 1 y 0 .
Thus from Corollary 2 we obtain the following:
(1) 
If I 0 1 or I 0 > 1 , then every nonconstant solution of Equation (39) is unbounded.
(2) 
If I 0 = 1 , then every nonconstant solution of Equation (39) is of period 10.
(3) 
If | I 0 | < 1 , then every nonconstant solution of Equation (39) converges to a limit.
Suppose that y 1 = 1 / 2 , y 0 = 1 / 4 , then the corresponding solution of Equation (40) is , 3 / 2 , 2 , 2 , 1 / 2 , 1 / 4 , , and so I 0 = 3 / 4 < 1 . Thus every nonconstant solution of Equation (39) converges to a limit even though g 1 = y n + 1 < 1 for some n 0 and for other n, g 1 is strictly greater than 1.
Remark 4.
The same idea applies to Equation (15) with D n = 0 , n = 0 , 1 , . For example, if there exists a constant r such that r < 1 and 0 < i = 0 j F i r j , j = 0 , 1 , , then the solutions of Equation (15) will converge to a limit even though there may be some F i > 1 for some i = 0 , , j . Since Equation (15) is used extensively to obtain our results, the behavior of the solutions of Equations (1), (17), (18) and (24) will depend on the product of the coefficients rather than the individual coefficients.
If for some n , A n = 1 , B n = 0 , then the solution of Equation (27) eventually becomes constant as the following result shows.
Theorem 3.
Suppose that A n + B n = 1 for all n 0 and A k = 1 , B k = 0 for some k 0 . Let N be the first k such that B n = 0 . Then the solution { x n } of Equation (27) eventually becomes a constant. More precisely for n N
x n = 1 + j = 0 N 1 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 N 1 ( 1 ) j i = 0 j B i x 1 .
Proof. 
Let N 0 be the first n such that B n = 0 , that is B N 1 0 , B N = 0 . Then Equation (17) with D n = 0 , n 0 becomes for N 0
x N = 1 + j = 0 N 1 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 N 1 ( 1 ) j i = 0 j B i x 1 .
Note that if N = 0 , then B 0 = 0 and Equation (42) is satisfied. From Equation (42) we get that for N 0
x N + 1 = 1 + j = 0 N ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 N ( 1 ) j i = 0 j B i x 1 = 1 + j = 0 N 1 ( 1 ) j + 1 i = 0 j B i x 0 + j = 0 N 1 ( 1 ) j i = 0 j B i x 1 = x N .
Similarly we can show that x n = x N , n N . ☐

5. Conclusions

In this paper we consider the general nonautonomous second order difference Equation (1) where the functions g i : R k R , i = 0 , 1 ; k = 1 , 2 , satisfy condition (2). The condition (2) implies the existence of an invariant which then reduces Equation (1) to the first order nonautonomous difference equation such as Equation (13) which has closed form solution (14). Using the closed form solution we were able to obtain the global behavior of solutions of some nonautonomous second order difference equations, such as Equation (26). The different types of asymptotic dynamics that the general nonautonomous second order difference Equation (1) exhibits are either convergence to a constant or a to a period-two sequence, divergence to ± , or eventual equality to a constant or a period-two sequence. Some results are given for the case when the coefficients of Equation (1) are periodic sequences when one can reduce the dynamics to the autonomous case, but our technique here was based on explicit solutions formula, see Example 6. The presented results have potential to be extended to systems of difference equations.

Author Contributions

Both authors contributed equally in the writing of this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Janowski, E.; Kulenović, M.R.S. Attractivity and global stability for linearizable difference equations. Comput. Math. Appl. 2009, 57, 1592–1607. [Google Scholar] [CrossRef]
  2. Janowski, E.; Kulenović, M.R.S.; Silić, E. Periodic Solutions of Linearizable Difference Equations. Int. J. Differ. Equ. 2011, 6, 113–125. [Google Scholar]
  3. Bilgin, A.; Kulenović, M.R.S. Global Asymptotic Stability for Discrete Single Species Population Models. Discret. Dyn. Nat. Soc. 2017, 2017, 5963594. [Google Scholar] [CrossRef]
  4. Sedaghat, H. Nonlinear Difference Equations, Theory with Applications to Social Science Models; Mathematical Modelling: Theory and Applications; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003; p. 15. [Google Scholar]
  5. Sedaghat, H. Form Symmetries and Reduction of Order in Difference Equations; Advances in Discrete Mathematics and Applications; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  6. Elaydi, S. An Introduction to Difference Equations, 3rd ed.; Undergraduate Texts in Mathematics; Springer: New York, NY, USA, 2005; 539p. [Google Scholar]
  7. Feuer, J.; Janowski, E.; Ladas, G. Invariants for some rational recursive sequences with periodic coefficients. J. Differ. Equ. Appl. 1996, 2, 167–174. [Google Scholar] [CrossRef]
  8. Gibbons, C.; Kulenović, M.R.S.; Ladas, G.; Voulov, H.D. On the trichotomy character of xn+1 = (α + βxn + γxn−1)/(A + xn). J. Differ. Equ. Appl. 2002, 8, 75–92. [Google Scholar] [CrossRef]
  9. Kalabušić, S.; Kulenović, M.R.S. Rate of convergence of solutions of rational difference equation of second order. Adv. Differ. Equ. 2004, 1, 121–140. [Google Scholar] [CrossRef]
  10. Kocic, V.L.; Ladas, G. Global Behavior of Nonlinear Difference Equations of Higher Order with Applications; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1993. [Google Scholar]
  11. Kulenović, M.R.S. Invariants and related Liapunov functions for difference equations. Appl. Math. Lett. 2000, 13, 1–8. [Google Scholar] [CrossRef]
  12. Kulenović, M.R.S.; Ladas, G. Dynamics of Second Order Rational Difference Equations, with Open Problems and Conjectures; Chapman & Hall/CRC Press: Boca Raton, FL, USA, 2001. [Google Scholar]
  13. Kulenović, M.R.S.; Merino, O. Discrete Dynamical Systems and Difference Equations with Mathematica; Chapman and Hall/CRC: Boca Raton, FL, USA; London, UK, 2002. [Google Scholar]
  14. Nesemann, T. Invariants and Liapunov functions for nonautonomous systems. Advances in difference equations, III. Comput. Math. Appl. 2001, 42, 385–392. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Janowski, E.J.; Kulenović, M.R.S. Global Behavior of Certain Nonautonomous Linearizable Three Term Difference Equations. Mathematics 2018, 6, 79. https://doi.org/10.3390/math6050079

AMA Style

Janowski EJ, Kulenović MRS. Global Behavior of Certain Nonautonomous Linearizable Three Term Difference Equations. Mathematics. 2018; 6(5):79. https://doi.org/10.3390/math6050079

Chicago/Turabian Style

Janowski, E. J., and M. R. S. Kulenović. 2018. "Global Behavior of Certain Nonautonomous Linearizable Three Term Difference Equations" Mathematics 6, no. 5: 79. https://doi.org/10.3390/math6050079

APA Style

Janowski, E. J., & Kulenović, M. R. S. (2018). Global Behavior of Certain Nonautonomous Linearizable Three Term Difference Equations. Mathematics, 6(5), 79. https://doi.org/10.3390/math6050079

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop