Javascript required
Skip to content Skip to sidebar Skip to footer

Find a Closed Form Solution for the Difference Equation

  • Research
  • Open Access
  • Published:

Closed-form solution of a general three-term recurrence relation

Advances in Difference Equations volume 2014, Article number:196 (2014) Cite this article

  • 5675 Accesses

  • 2 Citations

  • Metrics details

Abstract

We present a closed-form solution for n th term of a general three-term recurrence relation with arbitrary given n-dependent coefficients. The derivation and corresponding proof are based on two approaches, which we develop and describe in detail. First, the recursive-sum theory, which gives the exact solution in a compact finite form using a recursive indexing. Second, the discrete dimensional-convolution procedure, which transforms the solution to the non-recursive expression of n, including a finite number of elementary operations and functions.

1 Introduction

A general three-term recurrence relation is usually defined by the following expression:

W n + 1 = A n W n + B n W n 1 ,

(1)

where n 1 ( n N ), W n is unknown function of n, A n and B n are arbitrary given functions of n, and W 0 = C 0 , W 1 = C 1 are initial conditions (we assume not to have the trivial case when C 1 = C 0 = 0 and W n 0 ).

This well-known relation has a large number of applications and plays an important role in many areas of mathematics and physics. We recall here a few of them to underline the importance of the proposed statement.

First, the recurrence relation Eq. (1) corresponds to the finite difference equation for the general second order differential equation with unknown function f ( x ) and arbitrary given U ( x ) , namely: f U ( x ) f = 0 . Therefore, it is widely used for the analytical and numerical analysis (and approximations) in corresponding physical and mathematical applications; see for example [1].

Second, three-term recurrence relations appear naturally when one uses the Frobenius method for solving some linear differential equations and studying some special functions; see [2].

Third, Eq. (1) corresponds to a continued fraction with arbitrary given coefficients. Basically, the three-term recurrence relation corresponds to the expressions for the numerators and denominators of this continued fraction (see Chapter IV in [3]), which were derived by Euler.

In this manuscript we obtain a closed-form solution of a canonical three-term recurrence relation, which is equivalent to Eq. (1) in the case of ( A n 0 , n). Our goal is to obtain the solution consisting of a finite number of terms, rather than a variety of methods with infinite series [4]. In that way we develop two approaches for the expressions with a finite number of terms. We start with the recursive-sum theory, which is presented in Section 3. It allows one to obtain the exact solution of the three-term recurrence relation in a compact form using recursive indexing. Next, we develop the discrete dimensional-convolution procedure, which allows one to eliminate recursive indexing and to represent the solution as a closed-form expression. It means that the final closed-form expression depends on arbitrary given coefficients and includes a finite number of elementary operations, such as: '+', '−', '×', '÷'; and elementary functions, such as the Heaviside step function (or unit step function) and the floor function for integer division.

2 Preliminaries

Let us start from the simplification of Eq. (1). For this we assume that A n 0 , n (in the opposite case, the recurrence relation near numbers of zero A-coefficients has other fundamental properties and it should be considered separately). Thus, we can use the following substitution:

W n + 1 = a n + 1 i = 1 n A i .

(2)

According to this substitution and Eq. (1), a n fulfills the following three-term recurrence relation:

a n + 1 = a n + B n A n A n 1 a n 1 .

(3)

Since Eq. (1), Eq. (3) are linear and n, A n , B n are arbitrary, we can assume without loss of generality the following initial conditions: a 0 = 1 , a 1 = 1 . Finally, we specify d n = B n ( A n A n 1 ) 1 as an arbitrary given function of n, and write below the expression, which we further call the canonical three-term recurrence relation:

a n + 1 = a n + d n a n 1 , a 0 = 1 , a 1 = 1 .

(4)

We note that the solution of Eq. (4) consists of a sum of d i -products ( i N ) with different powers. Here the product power p is equal to the total number of d i in the product. The maximum power of a n + 1 , i.e. the largest power among all the products in a n + 1 , is equal to ( n + 1 ) / 2 , since

( n + 1 ) / 2 = max ( n / 2 , 1 + ( n 1 ) / 2 ) , n ,

(5)

where we use Gauss notation for the integer division: x = x = max ( m Z m x ) . Because of this fundamental characteristic it is natural to express a n + 1 , i.e. the solution of Eq. (4), in the following form:

a n + 1 = 1 + p = 1 ( n + 1 ) / 2 S ( n , p ) ,

(6)

where S ( n , p ) includes all terms of the power p. In the next sections we demonstrate that each S ( n , p ) is equal to the corresponding recursive-sum, and we rigorously prove all the propositions.

3 Recursive-sum theory

In this section we start from the definition of a general recursive-sum (R-sum). Next, we present and prove some of its properties, which determine the R-sum algebra. It could be useful for solving different recurrence relations, but, particularly for the solution of the canonical three-term recurrence relation it is sufficient to use a particular case of the R-sum, namely the reduced R-sum, which is described in the next section.

Definition 1 A general R-sum (R) is defined for the ordered sequence ( d i , i Z ) by the following expression:

R = R ( N , k , Δ , Δ 0 ) = m = 1 k [ i m = i m 1 + Δ N + m Δ d i m ] ,

(7)

where ( N , k , Δ , Δ 0 Z ), i 0 Δ 0 , N is the order of the R-sum, k 1 is the power of the R-sum, Δ is the recursive shift of the R-sum, Δ 0 is the initial shift of the R-sum. For the correct definition it is also assumed that N Δ 0 and d i is defined for all indices used in Eq. (7). The product in Eq. (7) corresponds to the standard left-to-right order of multiplication, so the definition can be written also in the following way:

R ( N , 1 , Δ , Δ 0 ) = i 1 = Δ 0 + Δ N + Δ d i 1 , R ( N , 2 , Δ , Δ 0 ) = i 1 = Δ 0 + Δ N + Δ d i 1 i 2 = i 1 + Δ N + 2 Δ d i 2 , R ( N , 3 , Δ , Δ 0 ) = i 1 = Δ 0 + Δ N + Δ d i 1 i 2 = i 1 + Δ N + 2 Δ d i 2 i 3 = i 2 + Δ N + 3 Δ d i 3 , R ( N , k , Δ , Δ 0 ) = i 1 = Δ 0 + Δ N + Δ d i 1 i 2 = i 1 + Δ N + 2 Δ d i 2 i k 1 = i k 2 + Δ N + ( k 1 ) Δ d i k 1 i k = i k 1 + Δ N + k Δ d i k .

(8)

Next, we present some of the elementary properties of the R-sum, P.(1)-P.(3), which follow directly from the definition.

P.(1):

R ( N , k , Δ , N ) = j = 1 k d N + j Δ .

(9)

P.(2):

R ( N , k + 1 , Δ , Δ 0 ) = j = Δ 0 + Δ N + Δ d j R ( N + Δ , k , Δ , j ) .

(10)

P.(3):

R ( N , k , Δ , Δ 0 + 1 ) = R ( N , k , Δ , Δ 0 ) d Δ 0 + Δ R ( N + Δ , k 1 , Δ , Δ 0 + Δ ) .

(11)

Finally, we consider the R-sum key lemma: If N Δ 0 + 1 , then

R ( N , k + 1 , Δ , Δ 0 ) = R ( N 1 , k + 1 , Δ , Δ 0 ) + d N + ( k + 1 ) Δ R ( N , k , Δ , Δ 0 ) .

(12)

Proof of the key lemma First, we consider the difference D = [ R ( N , k + 1 , Δ , Δ 0 ) d N + ( k + 1 ) Δ R ( N , k , Δ , Δ 0 ) ] and try to eliminate the index numbers which correspond to zero terms in the original expression. For the correct procedure (when an upper bound in each sum is never less than the lower one) we will eliminate the zero terms and delete the corresponding index numbers step by step from i 1 to i k + 1 . We have

D = i 1 = Δ 0 + Δ N + Δ d i 1 i 2 = i 1 + Δ N + 2 Δ d i 2 i k = i k 1 + Δ N + k Δ d i k [ i k + 1 = i k + Δ N + ( k + 1 ) Δ d i k + 1 d N + ( k + 1 ) Δ ] .

(13)

At the next step we note that for i 1 = N + Δ according to the definition we have only one value, i 2 = N + 2 Δ , and one term, d N + 2 Δ , in the second sum in Eq. (13). According to the definition and P.(1) the same is true for the other indices: i 3 = N + 3 Δ , , i k = N + k Δ . The last term in Eq. (13) gives [ i k + 1 = N + k Δ + Δ N + ( k + 1 ) Δ d i k + 1 d N + ( k + 1 ) Δ ] 0 . This means that all the terms in Eq. (13) which correspond to i 1 = N + Δ are equal to zero. Since the condition ( N Δ 0 + 1 ), we can then write

D = i 1 = Δ 0 + Δ ( N 1 ) + Δ d i 1 i 2 = i 1 + Δ N + 2 Δ d i 2 i k = i k 1 + Δ N + k Δ d i k [ i k + 1 = i k + Δ N + ( k + 1 ) Δ d i k + 1 d N + ( k + 1 ) Δ ] .

(14)

The same procedure can be repeated for any successive indices. Consider the m th sum ( m k ) in Eq. (13). For each i m = N + m Δ we have only one term in the next sum, corresponding to i m + 1 = N + ( m + 1 ) Δ , and so on until the last term: [ i k + 1 = N + k Δ + Δ N + ( k + 1 ) Δ d i k + 1 d N + ( k + 1 ) Δ ] 0 . So we can delete all the corresponding index numbers and write

D = i 1 = Δ 0 + Δ ( N 1 ) + Δ d i 1 i 2 = i 1 + Δ ( N 1 ) + 2 Δ d i 2 i k = i k 1 + Δ ( N 1 ) + k Δ d i k [ i k + 1 = i k + Δ N + ( k + 1 ) Δ d i k + 1 d N + ( k + 1 ) Δ ] .

(15)

Finally, we note that [ i k + 1 = i k + Δ N + ( k + 1 ) Δ d i k + 1 d N + ( k + 1 ) Δ ] = i k + 1 = i k + Δ ( N 1 ) + ( k + 1 ) Δ d i k + 1 . Thus we obtain

D = i 1 = Δ 0 + Δ ( N 1 ) + Δ d i 1 i 2 = i 1 + Δ ( N 1 ) + 2 Δ d i 2 i k = i k 1 + Δ ( N 1 ) + k Δ d i k i k + 1 = i k + Δ ( N 1 ) + ( k + 1 ) Δ d i k + 1 = R ( N 1 , k + 1 , Δ , Δ 0 ) ,

(16)

which proves the key lemma. □

4 Exact solution of a canonical three-term recurrence relation via finite R-sum expansion

In this section we start with a definition of a reduced R-sum, a particular case of the general R-sum. Then we construct the exact solution of the canonical three-term recurrence relation, by using the R-sum key lemma.

Definition 2 A reduced R-sum ( R ˜ ) is a particular case of the general R-sum (see Eqs. (7)-(8)), defined by the following expression:

R ˜ = R ˜ ( N , k ) = R ( N , k , 2 , 1 ) = i 1 = 1 N + 2 d i 1 i 2 = i 1 + 2 N + 4 d i 2 i k 1 = i k 2 + 2 N + 2 ( k 1 ) d i k 1 i k = i k 1 + 2 N + 2 k d i k .

(17)

The key lemma Eq. (12) gives the following for the reduced R-sum:

R ˜ ( N , k + 1 ) = R ˜ ( N 1 , k + 1 ) + d N + 2 ( k + 1 ) R ˜ ( N , k ) .

(18)

Since N and k are independent and arbitrary, let us consider them as functions of new numbers n and p: N = N ( n , p ) = n 2 p ; k = k ( n , p ) = p . In that way we consider a new subsidiary function:

S ( n , p ) = R ˜ ( n 2 p , p ) .

(19)

According to the previous notations we have N ( n , p ) = n 2 p = n 2 k = N ( n + 2 , p + 1 ) = N ( n + 1 , p + 1 ) + 1 . Thus we can write the analog of the key lemma for the S function, by using Eqs. (18) and (19):

S ( n + 2 , p + 1 ) = S ( n + 1 , p + 1 ) + d n + 2 S ( n , p ) .

(20)

Now, we will rigorously prove a theorem about the exact solution of a canonical three-term recurrence relation Eq. (4).

Theorem 1 An exact solution of Eq. (4) is the following:

a n + 1 = 1 + p = 1 ( n + 1 ) / 2 S ( n , p ) .

(21)

Proof Our proof is based on mathematical induction. Base cases:

n = 1 : a 2 = 1 + d 1 = 1 + p = 1 1 S ( 1 , p ) , n = 2 : a 3 = 1 + d 1 + d 2 = 1 + p = 1 1 S ( 2 , p ) , n = 3 : a 4 = 1 + d 1 + d 2 + d 3 + d 1 d 3 = 1 + p = 1 2 S ( 3 , p ) , n = 4 : a 5 = 1 + d 1 + d 2 + d 3 + d 4 + d 1 d 3 + d 1 d 4 + d 2 d 4 = 1 + p = 1 2 S ( 4 , p ) .

(22)

Inductive step: we will prove that the statement Eq. (21) for a n and a n 1 , n, gives the following:

a n + d n a n 1 = a n + 1 .

(23)

To prove it for all n, we consider below two cases. First case: n is an arbitrary even number, n = 2 m , m N . Second case: n is an arbitrary odd number, n = 2 m + 1 , m N .

First case, n = 2 m .

a n + d n a n 1 = a 2 m + d 2 m a 2 m 1 = 1 + d 2 m + l = 1 m S ( 2 m 1 , l ) + d 2 m p = 1 m 1 S ( 2 m 2 , p ) .

(24)

Note that l = 1 m S ( 2 m 1 , l ) = S ( 2 m 1 , 1 ) + l = 2 m S ( 2 m 1 , l ) = S ( 2 m 1 , 1 ) + p = 1 m 1 S ( 2 m 1 , p + 1 ) . After the substitution we have

a n + d n a n 1 = 1 + d 2 m + S ( 2 m 1 , 1 ) + p = 1 m 1 S ( 2 m 1 , p + 1 ) + d 2 m p = 1 m 1 S ( 2 m 2 , p ) = 1 + d 2 m + S ( 2 m 1 , 1 ) + p = 1 m 1 [ S ( 2 m 1 , p + 1 ) + d 2 m S ( 2 m 2 , p ) ] .

(25)

According to Eq. (20) we can write S ( 2 m 1 , p + 1 ) + d 2 m S ( 2 m 2 , p ) = S ( 2 m , p + 1 ) . Also, according to the definitions Eq. (17) and Eq. (19) we can write d 2 m + S ( 2 m 1 , 1 ) = d 2 m + j = 1 2 m 1 d j = j = 1 2 m d j = S ( 2 m , 1 ) . Thus, we can obtain

a n + d n a n 1 = 1 + S ( 2 m , 1 ) + p = 1 m 1 S ( 2 m , p + 1 ) = 1 + S ( 2 m , 1 ) + p = 2 m S ( 2 m , p ) = 1 + p = 1 m S ( 2 m , p ) = a n + 1 ,

(26)

which proves the first case of the theorem.

Second case, n = 2 m + 1 .

We have

a n + d n a n 1 = a 2 m + 1 + d 2 m + 1 a 2 m = 1 + d 2 m + 1 + l = 1 m S ( 2 m , l ) + d 2 m + 1 p = 1 m S ( 2 m 1 , p ) .

(27)

Similar to the previous case, we write l = 1 m S ( 2 m , l ) = S ( 2 m , 1 ) + p = 1 m 1 S ( 2 m , p + 1 ) . Also we note that p = 1 m S ( 2 m 1 , p ) = p = 1 m 1 S ( 2 m 1 , p ) + S ( 2 m 1 , m ) . Thus, we can write

a n + d n a n 1 = 1 + d 2 m + 1 + S ( 2 m , 1 ) + d 2 m + 1 S ( 2 m 1 , m ) + p = 1 m 1 [ S ( 2 m , p + 1 ) + d 2 m + 1 S ( 2 m 1 , p ) ] .

(28)

According to Eq. (20) we again can rewrite S ( 2 m , p + 1 ) + d 2 m + 1 S ( 2 m 1 , p ) = S ( 2 m + 1 , p + 1 ) . In addition we note that p = 1 m 1 S ( 2 m + 1 , p + 1 ) = p = 2 m S ( 2 m + 1 , p ) . Summarizing, we obtain

a n + d n a n 1 = 1 + d 2 m + 1 + S ( 2 m , 1 ) + d 2 m + 1 S ( 2 m 1 , m ) + p = 2 m S ( 2 m + 1 , p ) .

(29)

Finally, we note that d 2 m + 1 S ( 2 m 1 , m ) = d 2 m + 1 R ˜ ( 1 , m ) = d 2 m + 1 j = 1 m d 2 j 1 = R ˜ ( 1 , m + 1 ) = S ( 2 m + 1 , m + 1 ) . In addition, we can write d 2 m + 1 + S ( 2 m , 1 ) = d 2 m + 1 + j = 1 2 m d j = j = 1 2 m + 1 d j = S ( 2 m + 1 , 1 ) . This gives us the final result:

a n + d n a n 1 = 1 + S ( 2 m + 1 , 1 ) + S ( 2 m + 1 , m + 1 ) + p = 2 m S ( 2 m + 1 , p ) = 1 + p = 1 m + 1 S ( 2 m + 1 , p ) = a n + 1 ,

(30)

which proves the second case and Theorem 1. □

5 Discrete dimensional-convolution procedure

In this section we develop a procedure which could be used for the transformation of the exact solution Eq. (21) to the expression without recursive indexing.

Let us start from a simple case to demonstrate the idea of the discrete dimensional-convolution procedure. We consider a simple expression, which can be associated with a particular case of R-sum:

E 2 = i 1 = 1 N 1 i 2 = i 1 + δ N 2 f ( i 1 , i 2 ) ,

(31)

where N 1 , N 2 N ( N 1 N 2 ) are arbitrary given natural numbers; δ N ( 0 δ [ N 2 N 1 ] ) is an arbitrary given natural number (shift); and f ( i 1 , i 2 ) is an arbitrary given function of natural index numbers i 1 and i 2 . The previous expression E 2 can be rewritten as follows:

E 2 = i 1 = 1 N 1 i 2 = 1 N 2 F ( i 1 , i 2 ) ,

(32)

where F ( i 1 , i 2 ) = f ( i 1 , i 2 ) H ( i 2 i 1 δ ) , and H ( x ) is Heaviside step function (or unit step function):

H ( x ) = { 1 , x 0 , 0 , x < 0 .

(33)

Let us consider now Eq. (32). The indices i 1 and i 2 can be associated with two dimensions by the following way. Let us consider a two dimensional plot (array with elements F ( i 1 , i 2 ) ) with horizontal numbering related to the i 2 index and vertical numbering related to the i 1 index;

[ F ( 1 , 1 ) F ( 1 , 2 ) F ( 1 , N 2 ) F ( 2 , 1 ) F ( 2 , 2 ) F ( 2 , N 2 ) F ( N 1 , 1 ) F ( N 1 , 2 ) F ( N 1 , N 2 ) ] .

For numbering of this array of elements it is possible to use a traversal rule with one global index q [ 1 , , N 1 N 2 ] instead of i 1 [ 1 , , N 1 ] and i 2 [ 1 , , N 2 ] . Consider, for example, the following rule: q = i 1 + N 1 ( i 2 1 ) , with the corresponding relations:a i 1 = i 1 ( q ) = 1 + N 1 { ( q 1 ) / N 1 } = q N 1 ( q 1 ) / N 1 , i 2 = i 2 ( q ) = 1 + ( q 1 ) / N 1 . Since, for the presented rule, a unique combination of ( i 1 , i 2 ) corresponds to a certain unique number q and vice versa (i.e. we have one-to-one mapping), we can rewrite Eq. (32) in the following form:

E 2 = q = 1 N 1 N 2 F ( 1 + N 1 { ( q 1 ) / N 1 } , 1 + ( q 1 ) / N 1 ) .

(34)

The presented procedure, i.e. the reduction from two dimensions to one dimension, is a particular case of discrete dimensional-convolution procedure. Let us now generalize it to a multidimensional case. In that way we consider a multidimensional analog of Eq. (32):

E k = i 1 = 1 N 1 i 2 = 1 N 2 i k = 1 N k F ( i 1 , i 2 , , i k ) ,

(35)

where ( k N ), k 3 , is an arbitrary given natural number, ( N 1 , N 2 , , N k N ) is an array of arbitrary given natural numbers, and F ( i 1 , i 2 , , i k ) is an arbitrary given function of i 1 , i 2 , , i k .

Let us now define a global index q as follows:

q = i 1 + N 1 ( i 2 1 ) + N 1 N 2 ( i 3 1 ) + + [ j = 1 k 1 N j ] ( i k 1 ) .

(36)

Note that the minimum value q min = 1 corresponds to the case that every index number is equal to unity. The maximum value q max = N 1 N 2 N k corresponds to the case that every index number reaches its maximum, since we have the relation

N 1 + N 1 ( N 2 1 ) + N 1 N 2 ( N 3 1 ) + + [ j = 1 k 1 N j ] ( N k 1 ) = [ j = 1 k N j ] .

(37)

Now we construct one-to-one mapping between the global index q and the index numbers ( i 1 , , i k ) . For that we need to solve Eq. (36), i.e. to express any certain index number as a function of q (not depending on any other index numbers). We do that separately for i 1 , i k , and all others, i r ( 2 r [ k 1 ] , r N ).

Obtaining of i 1 : Since ( i 1 1 ) < N 1 , the expression { ( q 1 ) / N 1 } does not depend on any index numbers except i 1 . Thus, we obtain

i 1 = i 1 ( q ) = 1 + N 1 { ( q 1 ) / N 1 } .

(38)

Obtaining of i k : Since i 1 + N 1 ( i 2 1 ) + N 1 N 2 ( i 3 1 ) + + [ j = 1 k 2 N j ] ( i k 1 1 ) 1 < [ j = 1 k 1 N j ] (see Eq. (36)), the expression ( q 1 ) / j = 1 k 1 N j does not depend on any index numbers except i k . Thus we obtain

i k = i k ( q ) = 1 + ( q 1 ) / j = 1 k 1 N j .

(39)

Obtaining of i r ( 2 r [ k 1 ] ): Here we make two steps, similar to the previous cases, in order to eliminate the dependence of higher and lower index numbers separately. Consider the following expression:

h r = h r ( i 1 , , i r ) = i 1 + + [ j = 1 r 1 N j ] ( i r 1 ) .

(40)

According to Eq. (37) we have h r 1 < [ j = 1 r N j ] . In addition, we have the simple relations { u / v } v u for u < v , and { u / v } v 0 for ( u / v ) N . Thus we obtain

h r = 1 + { ( q 1 ) / j = 1 r N j } j = 1 r N j .

(41)

Next, we perform a step similar to the one we performed when obtaining i k (see also Eq. (39)):

i r = 1 + ( h r 1 ) / j = 1 r 1 N j .

(42)

Finally, we can write

i r = i r ( q ) = 1 + N r { ( q 1 ) / j = 1 r N j } .

(43)

As can be seen, Eq. (38) corresponds to Eq. (43) if we put formally r = 1 . In addition, we can write ( q 1 ) { ( q 1 ) / j = 1 k N j } j = 1 k N j . This means that Eq. (39) also corresponds to Eq. (43) if we put formally r = k . Thus we can write the described substitution for the discrete dimensional-convolution procedure in the following form:

q = i 1 + N 1 ( i 2 1 ) + N 1 N 2 ( i 3 1 ) + + [ j = 1 k 1 N j ] ( i k 1 ) , q [ 1 , , j = 1 k N j ] ; i r = 1 + N r { ( q 1 ) / j = 1 r N j } , r [ 1 , , k ] .

(44)

According to the previous equations every certain unique index number combination corresponds to the unique value of q. In addition, the total number of different index number combinations is equal to the total number of different values of q. Thus, we can conclude that Eq. (44) describes the one-to-one mapping for Eq. (35), and we can write

E k = q = 1 N 1 N 2 N k F ( i 1 ( q ) , , i k ( q ) ) .

(45)

Summarizing, the proposed procedure allows one to calculate a reduced R-sum without recursive indexing.

6 Closed-form solution of the canonical three-term recurrence relation

In this section we apply the previous results, in order to obtain the closed-form solution of Eq. (4) without recursive indexing. By using the exact solution via R-sum expansion Eq. (21), the definition of the reduced R-sum Eq. (17), and the procedure described in the previous section, we can write

a n + 1 = 1 + p = 1 ( n + 1 ) / 2 R ˜ ( n 2 p , p ) , R ˜ ( n 2 p , p ) = i 1 = 1 N 1 i 2 = 1 N 2 i p = 1 N p ( m = 1 p [ d i m H ( i m i m 1 2 ) ] ) R ˜ ( n 2 p , p ) = i 1 = 1 N 1 i 2 = 1 N 2 i p = 1 N p F ˜ ( i 1 , i 2 , , i p ) ,

(46)

where i 0 1 , N j = n 2 p + 2 j ( 1 j p , j N ), and H ( x ) is the unit step function (see Eq. (33)).

Now, we apply the results from the previous section Eqs. (44)-(45) in order to replace the recursive indexing by one global index q. For that we introduce new integer-valued function g = g ( m , n , p , q ) which is associated with a recursive index i m ( 1 m p , m N ):

g ( m , n , p , q ) = 1 + ( n 2 p + 2 m ) { ( q 1 ) / j = 1 m ( n 2 p + 2 j ) } .

(47)

By using this function and Eqs. (45)-(46) we can obtain finally the closed-form solutionb

a n + 1 = 1 + p = 1 ( n + 1 ) / 2 q = 1 M ( n , p ) m = 1 p [ d g ( m , n , p , q ) H [ g ( m , n , p , q ) g ( m 1 , n , p , q ) 2 ] ] , where M ( n , p ) = j = 1 p N j = j = 1 p ( n 2 p + 2 j ) = n ! ! ( n 2 p + 2 ) ( n 2 p + 2 ) ! ! .

(48)

Remark The integer-valued function g ( m , n , p , q ) [ 1 , , n ] determines the index number of arbitrary given d (index number) in the canonical three-term recurrence relation; see Eq. (4). The first term in the d-products in Eq. (48) is d g ( 1 , n , p , q ) according to Eq. (46) and the corresponding condition i 0 1 , so for the brevity in our notation we imply: g ( 0 , n , p , q ) 1 .

7 Fibonacci numbers

An important particular case of the canonical three-term recurrence relation is the well-known Fibonacci sequence, which corresponds to the case d n 1 in Eq. (4). Thus, the Fibonacci numbers ( F n ) satisfy the following relation:

F n + 1 = F n + F n 1 , F 0 = 1 , F 1 = 1 .

(49)

The Fibonacci number ( F n + 1 ) corresponds to the total number of terms (d-products) in the solution a n + 1 of the canonical three-term recurrence relation since the exact solution consists of a sum of unique combinations of d i products; see Eqs. (4)-(5) and Eqs. (46)-(48).

The closed-form solution of Fibonacci sequence has well been known at least since the discovery of Abraham de Moivre, however, we can obtain it also by using our results, namely, Eq. (21) or Eq. (48):

F n + 1 = 1 + p = 1 ( n + 1 ) / 2 ( i 1 = 1 n 2 p + 2 i 2 = i 1 + 2 n 2 p + 4 i p 1 = i p 2 + 2 n 2 i p = i p 1 + 2 n 1 ) = 1 + p = 1 ( n + 1 ) / 2 q = 1 M ( n , p ) m = 1 p H [ g ( m , n , p , q ) g ( m 1 , n , p , q ) 2 ] ,

(50)

where g ( m , n , p , q ) and M ( n , p ) are introduced in Eqs. (47)-(48), and H ( x ) is the unit step function; see Eq. (33).

This expression gives the exact solution of the Fibonacci sequence without involving irrational numbers; the solution includes only certain combinations of natural numbers.

8 Conclusions

In summary, we obtain the closed-form solution of a canonical three-term recurrence relation Eq. (4) with an arbitrary given n-dependent coefficient d n . The final non-recursive expression Eq. (48) includes a finite number of elementary operations and functions.

Possible applications of the developed approaches, namely the R-sum theory and the discrete dimensional-convolution procedure, are not limited by the considered statement. Due to its universality, they could be used for solving other recursive problems, in particular many-term recurrence relations.

An interesting and open question for the author is how the solution Eq. (48) could be efficiently used for approximations and solving differential equations.

Endnotes

We use the following notations: { x } = x x , x = x = max ( m Z m x ) .

We define double factorial for an arbitrary number m N as follows: m ! ! = m ( m 2 ) ( m 4 ) 1 .

References

Acknowledgements

IG would like to thank his wife Alla for her help and understanding during the preparation of the manuscript.

Author information

Affiliations

  1. Department of Physics, Umeå University, Umeå, 90187, Sweden

    Ivan Gonoskov

Corresponding author

Correspondence to Ivan Gonoskov.

Additional information

Competing interests

The author declares that he has no competing interests.

Authors' contributions

IG has made all substantive intellectual contributions to the published study.

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gonoskov, I. Closed-form solution of a general three-term recurrence relation. Adv Differ Equ 2014, 196 (2014). https://doi.org/10.1186/1687-1847-2014-196

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI : https://doi.org/10.1186/1687-1847-2014-196

Keywords

  • Exact Solution
  • Recurrence Relation
  • Index Number
  • Global Index
  • Elementary Operation

Find a Closed Form Solution for the Difference Equation

Source: https://advancesindifferenceequations.springeropen.com/articles/10.1186/1687-1847-2014-196