Next Article in Journal
Visual Inspection Method for Metal Rolls Based on Multi-Scale Spatial Location Feature
Next Article in Special Issue
Numerical Integration of Stiff Differential Systems Using Non-Fixed Step-Size Strategy
Previous Article in Journal
A Characterization for the Neighbor-Distinguishing Index of Planar Graphs
Previous Article in Special Issue
On Nonlinear Biharmonic Problems on the Heisenberg Group
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symmetries, Reductions and Exact Solutions of a Class of (2k + 2)th-Order Difference Equations with Variable Coefficients

by
Mensah Folly-Gbetoula
School of Mathematics, University of the Witwatersrand, Johannesburg 2000, South Africa
Symmetry 2022, 14(7), 1290; https://doi.org/10.3390/sym14071290
Submission received: 29 May 2022 / Revised: 13 June 2022 / Accepted: 17 June 2022 / Published: 21 June 2022

Abstract

:
We perform a Lie analysis of ( 2 k + 2 ) th-order difference equations and obtain k + 1 non-trivial symmetries. We utilize these symmetries to obtain their exact solutions. Sufficient conditions for convergence of solutions are provided for some specific cases. We exemplify our theoretical analysis with some numerical examples. The results in this paper extend to some work in the recent literature.

1. Introduction

Rational ordinary difference equations have been studied in the literature by many researchers [1,2,3,4,5,6,7,8]. Recently, the well-known symmetry methods for differential equations have been extended to difference equations [9,10,11]. The concept is similar to the one developed for differential equations (see [12]), and it consists of solving the difference equations using the group of transformations that leave the equation invariant. Hydon [10] developed a systematic method that enables one to find the group of transformations for difference equations. Though his algorithm is valid for any given difference equation, he mainly applied it to second-order difference equations. This could be due to the fact that for higher-order difference equations, the calculations are cumbersome. For more references on recurrence equations via the symmetry approach, see [13,14,15].
In this paper, we investigate the solutions of difference equations of the form
x n + 1 = x n ( 2 k + 1 ) a n + b n x n k x n ( 2 k + 1 ) ,
for some arbitrary sequences a n and b n using a symmetry-based method.
Special cases of the above equation exist in the recent literature. In [1], Cinar studied the positive solutions of the special case where a n = 1 , b n = 1 , and k = 0 :
x n + 1 = x n 1 1 + x n 1 x n .
In [2], the author studied the special case where a n = 1 , b n = 1 , and k = 0 :
x n + 1 = x n 1 1 + x n 1 x n .
In [5], the author studied the special case where a n = 1 , b n = a , and k = 0 :
x n + 1 = x n 1 1 + a x n 1 x n .
In [4], the author studied the special case where a n = 1 , b n = a , and k = 0 :
x n + 1 = x n 1 1 + a x n 1 x n .
In [3], the author studied the special case where a n = 1 / a , b n = b / a , and k = 0 :
x n + 1 = a x n 1 1 + b x n 1 x n .
In [15], the author obtained the solutions of Equation (6) with less restrictions on the initial conditions.
In [6], the author studied the special case where a n = a , b n = 1 , and k = 0 :
x n + 1 = x n 1 a x n 1 x n .
In [7], the case where a n = ± 1 , b n = ± 1 , and k = 1 such that
x n + 1 = x n 3 ± 1 ± x n 1 x n 3 ,
was studied by Elsayed, and exact solutions were obtained.
In [8], the author studied and obtained exact solutions of the special case where a n = 1 , b n = 1 , and k = 2 :
x n + 1 = x n 5 1 + x n 3 x n 5 .
In [16], the author studied the positive solutions and attractivenesss of the special case where a n = 1 and b n = 1 :
x n + 1 = x n ( 2 k + 1 ) 1 + x n k x n ( 2 k + 1 ) .
In [17], the author studied the solutions of the special case where a n = 1 and b n = 1 / a :
x n + 1 = a x n ( 2 k + 1 ) a + x n k x n ( 2 k + 1 ) .
For more work on difference equations, please see [18].
For definiteness, we study the difference equation
u n + 2 k + 2 = u n A n + B n u n + k + 1 u n ,
instead of Equation (1). To derive the solutions for Equation (12) using a symmetry-based method, we will first find the Lie group of point transformations of Equation (12). Then, we will reduce the order via the invariants and construct the solutions. Furthermore, we will explain how the solutions for Equation (1) can be obtained from the solutions for Equation (12). Finally, we will show how one can obtain the results in the literature using our results.

Preliminaries

In this section, we provide background to Lie symmetry analysis of difference equations. In this paper, we adopt the same notation as that in [10].
Definition 1.
Let G be a local group of transformations acting on a manifold M. A subset S M is called G-invariant, and G is called a symmetry group of S . If whenever x S , and g G is such that g · x is defined, then g · x S [12].
Definition 2.
Let G be a connected group of transformations acting on a manifold M. A smooth real-valued function ζ : M R is an invariant function for G if and only if
X ( ζ ) = 0 f o r a l l x M ,
and every infinitesimal generator X of G [12].
Definition 3.
A parameterized set of point transformations
Γ ε : x x ^ ( x ; ε ) ,
where x = x i , i = 1 , , p are continuous variables is a one-parameter local Lie group of transformations if the following conditions are satisfied [10]:
  • Γ 0 is the identity map if x ^ = x when ε = 0 ;
  • Γ a Γ b = Γ a + b for every a and b sufficiently close to 0;
  • Each x i ^ can be represented as a Taylor series (in a neighborhood of ε = 0 that is determined by x), and therefore
x i ^ ( x : ε ) = x i + ε ξ i ( x ) + O ( ε 2 ) , i = 1 , , p .
Consider the ordinary difference equation
u n + 2 k + 2 = ω ( n , u n , u n + 1 , , u n + 2 k + 1 ) , n D
for some smooth function ω and a regular domain D Z . To find a symmetry group of Equation (13), we consider the group of point transformations given by
G ε : ( n , u n ) ( n , u n + ε Q ( n , u n ) ) ,
where ε is the parameter and Q is a continuous function which we shall refer to as a characteristic. Let
X = Q ( n , u n ) u n + Q ( n + 1 , u n + 1 ) u n + 1 + Q ( n + 2 , u n + 2 ) u n + 2 + + Q ( n + 2 k + 1 , u n + 2 k + 1 ) u n + 2 k + 1
be the corresponding infinitesimal of G ε . The substitution of the new variable (obtained using G ε given in Equation (14)) in Equation (13) yields the linearized symmetry condition
S ( 2 k + 2 ) Q ( n , u n ) X ω = 0
whenever Equation (13) is true. The shift operator S acts on n as follows: S : n n + 1 . Once the characteristic Q = Q ( n , u n ) is known, the invariant V n may be obtained by solving the characteristic system [12]
d u n Q = d u n + 1 S Q = = d u n + 2 k + 1 S n + 2 k + 1 Q = d V n 0
or by introducing the canonical coordinate [19]
s n = d u n Q ( n , u n ) .
We use the standard conventions
j = s l θ j = 1 when s > l and j = s l θ j = 0 when s > l .

2. Symmetries

Consider the ( 2 k + 2 ) th-order difference Equation (12), which is
u n + 2 k + 2 = u n A n + B n u n + k + 1 u n .
To obtain the symmetries, we impose the infinitesimal criterion of invariance (Equation (15)) to obtain
Q ( n + 2 k + 2 , u n + 2 k + 2 ) + B n u n 2 ( B n u n u n + k + 1 + A n ) 2 Q ( n + k + 1 , u n + k + 1 ) A n ( B n u n u n + k + 1 + A n ) 2 Q ( n , u n ) = 0 .
The latter is a functional equation for the characteristic Q, making Equation (16) difficult to solve. To eliminate the first undesirable argument u n + 2 k + 2 , we differentiate implicitly with respect to u n (keeping u n + 2 k + 2 fixed and regarding u n + k + 1 as a function of u n and u n + 2 k + 2 ). This leads to
1 ( B n u n u n + k + 1 + A n ) 2 u n + k + 1 Q ( n + k + 1 , u n + k + 1 ) 1 ( B n u n u n + k + 1 + A n ) 2 u n Q ( n , u n ) + 2 u n ( B n u n u n + k + 1 + A n ) 2 Q ( n , u n ) = 0
After simplification, we obtain
u n + k + 1 Q ( n + k + 1 , u n + k + 1 ) u n Q ( n , u n ) + 2 u n Q ( n , u n ) = 0 .
To eliminate the second undesirable variable u n + k + 1 , we differentiate with respect to u n , and we obtain the following second-order differential equation involving only the argument u n :
d 2 d u n 2 Q ( n , u n ) 2 u n d d u n Q ( n , u n ) + 2 u n 2 Q ( n , u n ) = 0 .
The general solution of Equaiton (17) is then given by
Q = f ( n ) u n + g ( n ) u n 2 ,
where f and g are arbitrary functions of n. To eliminate the dependency between these two functions, we substitute Equation (18) into Equation (16) and simplify the resulting equation to obtain
B n g ( n + k + 1 ) u n + k + 1 2 u n 2 + B n f ( n + k + 1 ) u n + k + 1 u n 2 + B n f ( n + 2 k + 2 ) u n + k + 1 u n 2 A n g ( n ) u n 2 A n f ( n ) u n + A n f ( n + 2 k + 2 ) u n + g ( n + 2 k + 2 ) u n 2 = 0 .
We separate by powers of shifts of u n , and we readily reduce the resulting over-determining system to obtain
f ( n ) + f ( n + k + 1 ) = 0 , g ( n ) = 0 .
The solutions for Equation (19) are
exp ± i ( 2 s + 1 ) n π k + 1 , 0 s ( k 1 ) / 2 if k odd ( 1 ) n , exp ± i ( 2 s + 1 ) n π k + 1 , 0 s ( k 2 ) / 2 if k even .
This means that we have k + 1 non-trivial characteristics
Q s = β s n u n , Q ^ s = β ¯ s n u n , 0 s k 1 2 if k odd
and
Q 1 = ( 1 ) n u n , Q s = β s n u n , Q ^ s = β ¯ s n u n , 0 s k 2 2 if k even ,
where β s = exp i ( 2 s + 1 ) π k + 1 . We will refer to X i as the corresponding symmetry generator of Q i .

3. Exact Solutions

Here, we introduce the canonical coordinate (i.e., the variable S n such that X S n = 1 ). We use the well-known choice [19]
S n = d u n Q ( n , u n ) .
Using any one of the symmetry generators, such as X s , we have that
S n = d u n Q s ( n , u n ) = 1 ( β s ) n ln | u n | .
Using Equation (19), we have proven that
X s ( β s ) n + k + 1 S n + k + 1 + ( β s ) n S n = 0
Therefore, the following is an invariant function of X s :
r n = ( β s ) n + k + 1 S n + k + 1 + ( β s ) n S n
Since the equation under study is a rational difference equation, it is favorable to use
| r ˜ n | = exp ( r n ) ,
In other words, r ˜ n = ± 1 / u n u n + k + 1 . It turns out (using the plus sign together with Equation (12)) that
r ˜ n + k + 1 = A n r ˜ n + B n .
It is easy to verify that the solution for Equation (22) in closed form is given by
r ˜ ( k + 1 ) n + j = r ˜ j k 1 = 0 n 1 A ( k + 1 ) k 1 + j + l = 0 n 1 B ( k + 1 ) l + j k 2 = l + 1 n 1 A ( k + 1 ) k 2 + j ,
for j = 0 , 1 , , k . By going up the hierarchy created by the change of variables, we obtain the solution for Equation (12) in a unified manner as follows:
  • For k being odd:
    | u n | = exp { s = 0 k 1 2 c s β s n + c ˜ s β ¯ s n + s = 0 k 1 2 k s = 0 n 1 1 k + 1 β s n β ¯ s k s ln r ˜ k s + s = 0 k 1 2 k s = 0 n 1 1 k + 1 β s ¯ n β s k s ln r ˜ k s } | u n | = Γ n exp s = 0 k 1 2 k s = 0 n 1 2 k + 1 Re γ s ( n , k s ) ln r ˜ k s
  • For k being even:
    | u n | = θ n exp { s = 0 k 2 2 k s = 0 n 1 2 k + 1 Re γ s ( n , k s ) ln r ˜ k s + ( 1 ) n j = 0 n 1 1 k + 1 ( 1 ) j ln r ˜ ( k + 1 ) j }
where
Γ n = exp s = 0 k 1 2 c s β s n + c ˜ s β ¯ s n , θ n = exp c ( 1 ) n + s = 0 k 2 2 c s β s n + c ˜ s β ¯ s n , γ s ( n , k s ) = β s n β ¯ s k s .
The properties of Γ n , θ n and γ ( n , k ) are as follows:
Γ ( 2 k + 2 ) n + j = Γ j = | u j | , 0 j k Γ ( 2 k + 2 ) n + k + 1 + j = 1 Γ j = 1 | u j | , 0 j k
θ ( 2 k + 2 ) n + j = θ j = | u j | , 0 j k θ ( 2 k + 2 ) n + k + 1 + j = 1 θ j = 1 | u j | , 0 j k
γ s ( ( 2 k + 2 ) n , j ) = γ s ( 0 , j ) , γ s ( n + k + 1 , j ) = γ s ( n , j + k + 1 ) = γ s ( n , j ) , γ s ( j , j ) = 1 .
In regard to these properties, one can further simplify the solution through Equations (24) and (25). We have that
| u ( 2 k + 2 ) n | = Γ ( 2 k + 2 ) n exp s = 0 k 1 2 k s = 0 ( 2 k + 2 ) n 1 2 k + 1 Re γ s ( ( 2 k + 2 ) n , k s ) ln r ˜ k s = u 0 exp s = 0 k 1 2 k s = 0 ( 2 k + 2 ) n 1 2 k + 1 Re γ s ( 0 , k s ) ln r ˜ k s = u 0 exp { 2 k + 1 k 0 = 0 ( 2 k + 2 ) n 1 Re [ exp k 0 k + 1 i π + exp 3 k 0 k + 1 i π + + exp k k 0 k + 1 i π ] ln r ˜ k 0 } = u 0 exp { 2 k + 1 k 0 = 0 2 n 1 Re [ exp k 0 i π + exp 3 k 0 i π + + exp k k 0 i π ] ln r ˜ k 0 } = u 0 exp k 0 = 0 2 n 1 ( 1 ) k 0 ln r ˜ k 0 ( k + 1 ) = u 0 k 0 = 0 n 1 r ˜ 2 k 0 ( k + 1 ) r ˜ ( 2 k 0 + 1 ) ( k + 1 ) .
Note: It can be shown, using the expression of r ˜ n given in Equation (21), that one can write Equation (26) without using the symbol of the absolute value.
On this note, and in the same way, we have proven that
u ( 2 k + 2 ) n + j = u j k 0 = 0 n 1 r ˜ 2 k 0 ( k + 1 ) + j r ˜ ( 2 k 0 + 1 ) ( k + 1 ) + j , 0 j k ,
u ( 2 k + 2 ) n + k + 1 + j = u k + 1 + j k 0 = 0 n 1 r ˜ ( 2 k 0 + 1 ) ( k + 1 ) + j r ˜ ( 2 k 0 + 2 ) ( k + 1 ) + j , 0 j k .
We therefore obtain, by combining Equations (23), (27) and (28) and by replacing r ˜ i with 1 / ( u i u i + k + 1 ) , the solution for Equation (12) in closed form as
u ( 2 k + 2 ) n + j = u j k 0 = 0 n 1 k 1 = 0 2 k 0 1 A ( k + 1 ) k 1 + j + u j u j + k + 1 l = 0 2 k 0 1 B ( k + 1 ) l + j k 2 = l + 1 2 k 0 1 A ( k + 1 ) k 2 + j k 1 = 0 2 k 0 A ( k + 1 ) k 1 + j + u j u j + k + 1 l = 0 2 k 0 B ( k + 1 ) l + j k 2 = l + 1 2 k 0 A ( k + 1 ) k 2 + j ,
u ( 2 k + 2 ) n + k + 1 + j = u k + 1 + j k 0 = 0 n 1 k 1 = 0 2 k 0 A ( k + 1 ) k 1 + j + u j u j + k + 1 l = 0 2 k 0 B ( k + 1 ) l + j k 2 = l + 1 2 k 0 A ( k + 1 ) k 2 + j k 1 = 0 2 k 0 + 1 A ( k + 1 ) k 1 + j + u j u j + k + 1 l = 0 2 k 0 + 1 B ( k + 1 ) l + j k 2 = l + 1 2 k 0 + 1 A ( k + 1 ) k 2 + j ,
where 0 j k and provided that
u j u j + k + 1 l = 0 m B ( k + 1 ) l + j k 2 = l + 1 m A ( k + 1 ) k 2 + j k 1 = 0 m A ( k + 1 ) k 1 + j ,
where m < 2 n . Recall that we enacted the shift operator 2 k + 1 times in Equation (1) to obtain Equation (12), and therefore we can write the solutions for Equation (1) as follows
x ( 2 k + 2 ) n ( 2 k + 1 ) + j = x j 2 k 1 k 0 = 0 n 1 k 1 = 0 2 k 0 1 a ( k + 1 ) k 1 + j + x j 2 k 1 x j k l = 0 2 k 0 1 b ( k + 1 ) l + j k 2 = l + 1 2 k 0 1 a ( k + 1 ) k 2 + j k 1 = 0 2 k 0 a ( k + 1 ) k 1 + j + x j 2 k 1 x j k l = 0 2 k 0 b ( k + 1 ) l + j k 2 = l + 1 2 k 0 a ( k + 1 ) k 2 + j ,
x ( 2 k + 2 ) n k + j = x j k k 0 = 0 n 1 k 1 = 0 2 k 0 a ( k + 1 ) k 1 + j + x j 2 k 1 x j k l = 0 2 k 0 b ( k + 1 ) l + j k 2 = l + 1 2 k 0 a ( k + 1 ) k 2 + j k 1 = 0 2 k 0 + 1 a ( k + 1 ) k 1 + j + x j 2 k 1 x j k l = 0 2 k 0 + 1 b ( k + 1 ) l + j k 2 = l + 1 2 k 0 + 1 a ( k + 1 ) k 2 + j ,
provided that
x j 2 k 1 x j k l = 0 m b ( k + 1 ) l + j k 2 = l + 1 m a ( k + 1 ) k 2 + j k 1 = 0 m a ( k + 1 ) k 1 + j , m < 2 n .

3.1. Case with ( K + 1 ) -Periodic Sequences ( A n ) and ( B n )

Let ( a n ) = ( a 0 , a 1 , , a k , a 0 , a 1 , ) , ( b n ) = ( b 0 , b 1 , , b k , b 0 , b 1 , ) and Φ j = x j 2 k 1 x j k b j . The solutions given in Equations (32), (33) and (34) simplify to
x ( 2 k + 2 ) n + j ( 2 k + 1 ) = x j 2 k 1 k 0 = 0 n 1 a j 2 k 0 + Φ j l = 0 2 k 0 1 a j l a j 2 k 0 + 1 + Φ j l = 0 2 k 0 ( a j ) l ,
x ( 2 k + 2 ) n k + j = x j k k 0 = 0 n 1 a j 2 k 0 + 1 + Φ j l = 0 2 k 0 a j l a j 2 k 0 + 2 + Φ j l = 0 2 k 0 + 1 a j l ,
with
Φ j l = 0 2 k 0 ( a j ) l a j 2 k 0 + 1 , k 0 < n , Φ j l = 0 2 k 0 + 1 ( a j ) l a j 2 k 0 + 2 , k 0 < n , j k .

Case with A j 1

The solutions given in (35), (36) and (37) become
x ( 2 k + 2 ) n + j ( 2 k + 1 ) = x j 2 k 1 k 0 = 0 n 1 a j 2 k 0 + Φ j 1 ( a j ) 2 k 0 1 a j a j 2 k 0 + 1 + Φ j 1 ( a j ) 2 k 0 + 1 1 a j
x ( 2 k + 2 ) n k + j = x j k k 0 = 0 n 1 a j 2 k 0 + 1 + Φ j 1 ( a j ) 2 k 0 + 1 1 a j a j 2 k 0 + 2 + Φ j 1 ( a j ) 2 k 0 + 2 1 a j ,
with
Φ j 1 ( a j ) 2 k 0 + 1 ( 1 a j ) a j 2 k 0 + 1 ,
Φ j 1 ( a j ) 2 k 0 + 2 ( 1 a j ) a j 2 k 0 + 2 , k 0 < n , j k .

3.2. Case with One-Periodic Sequences ( A n ) and ( B n )

Let ( a n ) = ( a 0 , a 0 , ) , ( b n ) = ( b 0 , b 0 , ) and Φ j = x j 2 k 1 x j k b 0 . The solutions given in Equations (35), (36) and (37) simplify to
x ( 2 k + 2 ) n + j ( 2 k + 1 ) = x j 2 k 1 k 0 = 0 n 1 a 0 2 k 0 + Φ j l = 0 2 k 0 1 a 0 l a 0 2 k 0 + 1 + Φ j l = 0 2 k 0 ( a 0 ) l
x ( 2 k + 2 ) n k + j = x j k k 0 = 0 n 1 a 0 2 k 0 + 1 + Φ j l = 0 2 k 0 a 0 l a 0 2 k 0 + 2 + Φ j l = 0 2 k 0 + 1 ( a 0 ) l ,
with
Φ j l = 0 2 k 0 ( a 0 ) l a 0 2 k 0 + 1 , Φ j l = 0 2 k 0 + 1 ( a 0 ) l a 0 2 k 0 + 2 , k 0 < n , j k .

3.2.1. Case with A 0 1

The solutions given in Equations (42), (43) and (44) become
x ( 2 k + 2 ) n + j ( 2 k + 1 ) = x j 2 k 1 k 0 = 0 n 1 a 0 2 k 0 + Φ j 1 ( a 0 ) 2 k 0 1 a 0 a 0 2 k 0 + 1 + Φ j 1 ( a 0 ) 2 k 0 + 1 1 a 0
x ( 2 k + 2 ) n k + j = x j k k 0 = 0 n 1 a 0 2 k 0 + 1 + Φ j 1 ( a 0 ) 2 k 0 + 1 1 a 0 a 0 2 k 0 + 2 + Φ j 1 ( a 0 ) 2 k 0 + 2 1 a 0 ,
with
Φ j 1 ( a 0 ) 2 k 0 + 1 ( 1 a 0 ) a 0 2 k 0 + 1 , Φ j 1 ( a 0 ) 2 k 0 + 2 ( 1 a 0 ) a 0 2 k 0 + 2 ,
k 0 < n , j k .

3.2.2. Case with A 0 = 1

The solutions given in Equations (42), (43) and (44) become
x ( 2 k + 2 ) n + j ( 2 k + 1 ) = x j 2 k 1 k 0 = 0 n 1 1 + ( 2 k 0 ) Φ j 1 + ( 2 k 0 + 1 ) Φ j
x ( 2 k + 2 ) n k + j = x j k k 0 = 0 n 1 1 + ( 2 k 0 + 1 ) Φ j 1 + ( 2 k 0 + 2 ) Φ j ,
with
( 2 k 0 + 1 ) Φ j 1 , ( 2 k 0 + 2 ) Φ j 1 , k 0 < n , j k .

4. Results

In this section, we verify the results in [1,2,3,4,5,6,7,8,16,17] by utilizing different combinations of values of a n and b n in Equation (1):
  • If we set k = 0 , b 0 = 1 in Equations (42)–(44), we obtain the result in [1] (see Theorem 2.1) for Equation (2). However, the restriction ( x 1 and x 0 are positive real numbers) in [1] is a special case of our restrictions ( ( 2 k 0 + 1 ) x 1 x 0 1 and ( 2 k 0 + 1 ) x 0 x 1 1 , k 0 < n );
  • If we set k = 0 , b 0 = 1 in Equations (45)–(47), we obtain the result in [2] for Equation (3) (see Theorem 2.1). The restriction ( x 1 x 0 1 ) in [2] coincides with our restriction;
  • If we set k = 0 , b 0 = a in Equations (48)–(50), we obtain the result in [5] for Equation (4) (see Theorem 2.1). However, the restriction ( x 1 and x 0 are nonnegative real numbers) in [5] is a special case of our restrictions ( ( 2 k 0 + 1 ) a x 1 x 0 1 and ( 2 k 0 + 2 ) a x 1 x 0 1 , k 0 < n );
  • If we set k = 0 and b 0 = a in Equations (45)–(47), we obtain the result in [4] for Equation (5) (see Theorem 2.1), and the restriction ( a x 1 x 0 1 ) in [4] coincides with our restriction;
  • If we set k = 0 , a 0 = ( 1 / a ) , and b 0 = ( b / a ) in Equations (42)–(44) we obtain the result in [3] for Equation (6) (see Theorem 2.1). However, the restriction ( a , b , x 1 , and x 0 are nonnegative real numbers) in [3] is a special case of our restrictions ( b x 1 x 0 l = 0 2 k 0 a l 1 and b x 1 x 0 l = 0 2 k 0 + 1 a l 1 , k 0 < n );
  • If we set k = 0 , a 0 = a , and b 0 = 1 in Equations (42)–(44), we obtain the result (for the case a 1 ) in [6] for Equation (7) (see Theorem 2), and the restriction ( x 1 x 0 a j ( 1 a ) / ( 1 a j ) ) in [6] coincides with our restrictions ( x 1 x 0 1 a 2 k 0 + 1 ( 1 a ) a 2 k 0 + 1 , x 1 x 0 1 a 2 k 0 + 2 ( 1 a ) a 2 k 0 + 2 , and k 0 < n ). Additionally, the solution for the case a = 1 (see Theorem 5 in [6]) corresponds to our solution with the same restrictions on the initial conditions;
  • If we set k = 1 , a 0 = 1 , and b 0 = 1 (resp b 0 = 1 ) in Equations (48)–(50), we obtain the result in [7] for Equation (8) (see Theorem 1 (resp. Theorem 4)). However, the restriction ( x 1 and x 0 are nonzero positive real numbers (resp. j x 0 x 2 1 , j x 1 x 3 1 for j = 1 , 2 , 3 , )) in [7] is a special case of our restrictions ( ( 2 k 0 + 1 ) x 1 x 3 1 , ( 2 k 0 + 2 ) x 0 x 2 1 (resp. [ 2 k 0 + 1 ] x 1 x 3 1 , and [ 2 k 0 + 2 ] x 0 x 2 1 )). On the other hand, if a 0 = 1 and b 0 = ± 1 , the results are the same as in [7] (see Theorems 6 and 9) with the same restrictions on the initial conditions ( x 0 x 2 1 and x 1 x 3 1 );
  • If we set k = 2 , a 0 = 1 , and b 0 = 1 in Equations (45)–(47), we obtain the result in [8] for Equation (9) (see Theorem 1), and the restrictions ( x 0 x 3 1 , x 1 x 4 1 , and x 5 x 2 1 ) in [8] coincide with our restriction;
  • If we set a n = 1 and b n = 1 in Equations (32)–(34), we obtain the result in [16] for Equation (10) (see Theorem 1). The author in [16] restricted himself to positive solutions. However, our result is valid for negative solutions as well, provided that Equation (34) is satisfied.
  • If we set a 0 = 1 and b 0 = 1 / a in Equations (45)–(47), we obtain the result in [17] for Equation (11) (see Theorem 1). However, the restriction ( x 0 x ( k + 1 ) a ; x 1 x ( k + 2 ) a ; x 2 2 x ( k + 3 ) a , , x k x ( 2 k + 1 ) a ) in [17] is a special case of our restriction ( x j 2 k 1 x j k a , j k ).

5. Numerical Examples

Below are some graphs that show the behavior of the solutions.
Example 1.
For the case with A n = 1 , B n = B and k = 0 , where u n + 2 = u n 1 + B u n u n + 1 , setting A n = 1 , B n = B , and k = 0 in Equations (29) and (30) yields
u 2 n = u 0 k 0 = 0 n 1 1 + 2 k 0 B u 0 u 1 1 + ( 2 k 0 + 1 ) B u 0 u 1
and
u 2 n + 1 = u 1 k 0 = 0 n 1 1 + ( 2 k 0 + 1 ) B u 0 u 1 1 + ( 2 k 0 + 2 ) B u 0 u 1 .
A sufficient condition for this solution to converge is given by 1 + 2 k 0 B u 0 u 1 < 1 + ( 2 k 0 + 1 ) B u 0 u 1 and 1 + ( 2 k 0 + 1 ) B u 0 u 1 < 1 + ( 2 k 0 + 2 ) B u 0 u 1 for 0 k 0 n 1 . This implies that B u 0 u 1 > 0 . It is then clear that if B u 0 u 1 > 0 , then the solution converges. On the other hand, a sufficient condition for this solution to diverge is given by 1 + 2 k 0 B u 0 u 1 > 1 + ( 2 k 0 + 1 ) B u 0 u 1 or 1 + ( 2 k 0 + 1 ) B u 0 u 1 > 1 + ( 2 k 0 + 2 ) B u 0 u 1 for 0 k 0 n 1 . This implies that B u 0 u 1 < 0 . It is then clear that if B u 0 u 1 < 0 , then the solution diverges.
Figure 1 shows the graph of u n + 2 = u n 1 + u n u n + 1 with u 0 = 2 and u 1 = 5 . The solution converges as expected, since B u 0 u 1 > 0 .
Figure 2 shows the graph of u n + 2 = u n 1 0.28 u n u n + 1 with u 0 = 0.08 and u 1 = 0.67 . The solution diverges as expected, since B u 0 u 1 < 0 .
Example 2.
In the case with A n = A 1 , B n = B , and k = 1 such that u n + 4 = u n A + B u n u n + 2 , setting A n = A , B n = B , and k = 1 in Equations (29) and (30) yields
u 4 n = u 0 A + B u 0 u 2 k 0 = 1 n 1 A 2 k 0 + B u 0 u 2 l = 0 2 k 0 1 A l A 2 k 0 + 1 + B u 0 u 2 l = 0 2 k 0 A l ,
u 4 n + 1 = u 1 A + B u 1 u 3 k 0 = 1 n 1 A 2 k 0 + B u 1 u 3 l = 0 2 k 0 1 A l A 2 k 0 + 1 + B u 1 u 3 l = 0 2 k 0 A l ,
u 4 n + 2 = u 2 k 0 = 0 n 1 A 2 k 0 + 1 + B u 0 u 2 l = 0 2 k 0 A l A 2 k 0 + 2 + B u 0 u 2 l = 0 2 k 0 + 1 A l ,
u 4 n + 3 = u 3 k 0 = 0 n 1 A 2 k 0 + 1 + B u 1 u 3 l = 0 2 k 0 A l A 2 k 0 + 2 + B u 1 u 3 l = 0 2 k 0 + 1 A l .
Similarly, a sufficient condition for this solution to converge is for the numerators in Equations (51)–(54) to be less than their corresponding denominators. This implies that B u 0 u 2 > 1 A and B u 1 u 3 > 1 A , where A > 0 . It is then clear that if B u 0 u 2 > 1 A and B u 1 u 3 > 1 A , then the solution converges. On the other hand, a sufficient condition for this solution to diverge is for at least one of the numerators in Equations (51)–(54) to be greater than its corresponding denominator. This implies that B u 0 u 2 < 1 A or B u 1 u 3 < 1 A , where A > 0 . It is then clear that if B u 0 u 2 < 1 A or B u 1 u 3 < 1 A , then the solution diverges.
Figure 3 shows the graph of u n + 4 = u n 2 2 u n u n + 2 with u 0 = 0.63 , u 1 = 0.33 , u 2 = 0.72 , and u 3 = 0 . The solution converges as expected, since B u 0 u 2 > 1 A and B u 1 u 3 > 1 A .
Figure 4 shows the graph of u n + 4 = u n 0.78 + 0.18 u n u n + 2 with u 0 = 0.36 , u 1 = 0.63 , u 2 = 0.01 , and u 3 = 0 . The solution diverges as expected, since B u 0 u 2 < 1 A .

6. Conclusions

We obtained all the Lie point symmetries of the difference equation (Equation (12)). We used a symmetry-based method to derive its solutions. The solutions were given in ‘single’ form and then separated into ( 2 k + 2 ) categories. We explained how the solutions of Equation (1) can be obtained using those of Equation (12) to show that the results in the recent literature can be obtained as special cases of our generalized solutions.

Funding

This research was funded by the National Research Foundation (NRF) of South Africa, grant number: 132108.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Cinar, C. On the positive solutions of the difference equation xn+1 = xn−1/(1 + xnxn−1). Appl. Math. Comput. 2004, 150, 21–24. [Google Scholar]
  2. Cinar, C. On the difference equation xn+1 = xn−1/(−1 + xnxn−1). Appl. Math. Comput. 2004, 158, 813–816. [Google Scholar]
  3. Cinar, C. On the positive solutions of the difference equation xn+1 = axn−1/(1 + bxnxn−1). Appl. Math. Comput. 2004, 156, 587–590. [Google Scholar]
  4. Cinar, C. On the solutions of the difference equation xn+1 = xn−1/(−1 + axnxn−1). Appl. Math. Comput. 2004, 158, 793–797. [Google Scholar]
  5. Cinar, C. On the positive solutions of the difference equation xn+1 = xn−1/(1 + axnxn−1). Appl. Math. Comput. 2004, 158, 809–812. [Google Scholar]
  6. Aloqeili, M. Dynamics of a rational difference equation. Appl. Math. Comput. 2006, 176, 768–774. [Google Scholar] [CrossRef]
  7. Elsayed, E.M. On the solution of some difference equation. Eur. J. Pure Appl. Math. 2011, 4, 287–303. [Google Scholar]
  8. Elsayed, E.M. On the difference equation. Int. J. Contemp. Math. Sci. 2008, 3, 1657–1664. [Google Scholar]
  9. Quispel, G.R.W.; Sahadevan, R. Lie symmetries and the integration of difference equations. Phys. Lett. A 1993, 184, 64–70. [Google Scholar] [CrossRef]
  10. Hydon, P.E. Difference Equations by Differential Equation Methods; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  11. Levi, D.; Vinet, L.; Winternitz, P. Lie group formalism for difference equations. J. Phys. A Math. Gen. 1997, 30, 633–649. [Google Scholar] [CrossRef]
  12. Olver, P. Applications of Lie Groups to Differential Equations, 2nd ed.; Springer: New York, NY, USA, 1993. [Google Scholar]
  13. Folly-Gbetoula, M.; Kara, A.H. Symmetries, conservation laws, and integrability of difference equations. Adv. Differ. Eq. 2014, 2014, 1–14. [Google Scholar] [CrossRef] [Green Version]
  14. Hydon, P.E. Symmetries and first integrals of ordinary difference equations. Proc. R. Soc. Lond. A 2000, 456, 2835–2855. [Google Scholar] [CrossRef] [Green Version]
  15. Folly-Gbetoula, M. Symmetry, reductions and exact solutions of the difference equation un+2 = aun/(1 + bunun+1). J. Diff. Eq. Appl. 2017, 23, 1017–1024. [Google Scholar] [CrossRef]
  16. Yalcinkaya, I. On the global attractivity of positive solutions of a rational difference equation. Selcuk J. Appl. Math. 2008, 9, 3–8. [Google Scholar]
  17. Karatas, R. On the solutions of the recursive sequence x n + 1 = a x n ( 2 k + 1 ) a + x n k x n ( 2 k + 1 ) . Fasc. Math. 2010, 45, 37–45. [Google Scholar]
  18. Abdelrahman, M.A.E.; Chatzarakis, G.E.; Li, T.; Moaaz, O. On the difference equation xn+1 = axn−l + bxnk + f(xn−l,xnk). Adv. Differ. Eq. 2018, 431, 2018. [Google Scholar]
  19. Joshi, N.; Vassiliou, P. The existence of Lie Symmetries for First-Order Analytic Discrete Dynamical Systems. J. Math. Anal. Appl. 1995, 195, 872–887. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The graph of u n + 2 = u n 1 + u n u n + 1 .
Figure 1. The graph of u n + 2 = u n 1 + u n u n + 1 .
Symmetry 14 01290 g001
Figure 2. The graph of u n + 2 = u n 1 0.28 u n u n + 1 .
Figure 2. The graph of u n + 2 = u n 1 0.28 u n u n + 1 .
Symmetry 14 01290 g002
Figure 3. The graph of u n + 4 = u n 2 2 u n u n + 2 .
Figure 3. The graph of u n + 4 = u n 2 2 u n u n + 2 .
Symmetry 14 01290 g003
Figure 4. The graph of u n + 4 = u n 0.78 + 0.18 u n u n + 2 .
Figure 4. The graph of u n + 4 = u n 0.78 + 0.18 u n u n + 2 .
Symmetry 14 01290 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Folly-Gbetoula, M. Symmetries, Reductions and Exact Solutions of a Class of (2k + 2)th-Order Difference Equations with Variable Coefficients. Symmetry 2022, 14, 1290. https://doi.org/10.3390/sym14071290

AMA Style

Folly-Gbetoula M. Symmetries, Reductions and Exact Solutions of a Class of (2k + 2)th-Order Difference Equations with Variable Coefficients. Symmetry. 2022; 14(7):1290. https://doi.org/10.3390/sym14071290

Chicago/Turabian Style

Folly-Gbetoula, Mensah. 2022. "Symmetries, Reductions and Exact Solutions of a Class of (2k + 2)th-Order Difference Equations with Variable Coefficients" Symmetry 14, no. 7: 1290. https://doi.org/10.3390/sym14071290

APA Style

Folly-Gbetoula, M. (2022). Symmetries, Reductions and Exact Solutions of a Class of (2k + 2)th-Order Difference Equations with Variable Coefficients. Symmetry, 14(7), 1290. https://doi.org/10.3390/sym14071290

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