Next Article in Journal
A Design for Genetically Oriented Rules-Based Incremental Granular Models and Its Application
Next Article in Special Issue
Supersymmetric Higher Spin Models in Three Dimensional Spaces
Previous Article in Journal
Graph Cellular Automata with Relation-Based Neighbourhoods of Cells for Complex Systems Modelling: A Case of Traffic Simulation
Previous Article in Special Issue
Evaporation and Antievaporation Instabilities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boundary Value Problems for Some Important Classes of Recurrent Relations with Two Independent Variables

1
Mathematical Institute of the Serbian Academy of Sciences, Knez Mihailova 36/III, 11000 Beograd, Serbia
2
Faculty of Electrical Engineering, Belgrade University, Bulevar Kralja Aleksandra 73, 11000 Beograd, Serbia
3
Faculty of Electrical Engineering and Communication, Department of Mathematics, Brno University of Technology, Technicka 3058/10, CZ-616 00 Brno, Czech Republic
*
Author to whom correspondence should be addressed.
Symmetry 2017, 9(12), 323; https://doi.org/10.3390/sym9120323
Submission received: 3 December 2017 / Revised: 17 December 2017 / Accepted: 18 December 2017 / Published: 20 December 2017
(This article belongs to the Special Issue Symmetry: Feature Papers 2017)

Abstract

:
It is shown that complex-valued boundary value problems for several classes of recurrent relations with two independent variables, of some considerable interest, are solvable on the following domain: C = { ( n , k ) : 0 k n , k N 0 , n N } , the so called combinatorial domain. The recurrent relations include some of the most important combinatorial ones, which, among other things, serve as a motivation for the investigation. The methods for solving the boundary value problems are presented and explained in detail.

1. Introduction

As usual, N stands for the set of all positive integers { 1 , 2 , 3 , } , whereas N 0 = N { 0 } . If m and n are two numbers from N 0 such that m n , then for the set of all i N 0 such that m i n , we use the following notation i = m , n ¯ . Throughout the paper, as usual, we adopt the convention that:
j = k l a k = 0 ,
if k , l N 0 are such that l < k and that:
j = k k 1 a k = 1 ,
for each k N 0 .
Let:
H n = j = 1 n 1 j , n N ,
(the harmonic numbers).
By C k n , 0 k n , are denoted the Binomial coefficients, which can be introduced, for example, as the coefficients of the polynomial P n ( x ) = ( 1 + x ) n , n N 0 , that is, we have:
P n ( x ) = C 0 n + C 1 n x + + C n n x n .
Recall that, for the coefficients, the following formula holds:
C k n = n ! k ! ( n k ) ! ,
where:
m ! = j = 1 m j and 0 ! : = 1 .
Note that:
C 0 n = C n n = 1 , n N 0 .
One of the basic relations involving the coefficients is the following two-dimensional recurrent relation:
C k n = C k n 1 + C k 1 n 1 ,
for k , n N , such that 1 k < n , which is easily obtained from the following obvious identity ( 1 + x ) n 1 ( 1 + x ) = ( 1 + x ) n . Many other equalities and relations involving the coefficients, can be found, for example, in the following books: [1,2,3,4,5,6]. More specifically, some basic relations can be found in [1,2], while some more complex ones can be found in [3]. Book [4] is of an encyclopedic character and presents a large list of various relations. More complex methods, among others, for investigating relations with the coefficients can be found in [5]. Use of the coefficients in combinatorics, as well as their combinatorial meaning can be found in many books, for example, in [2,6].
Formula (1), the equalities in (2), as well as the recurrent relation (3) were the main motivations in [7], in which the relation (3), together with (2), was considered in a new way.
One of the main things in [7], connected to the recurrent relation (3), is the solvability, which is an old topic (see the classics [1,3,5,8,9,10] for well-known methods for solving difference equations). The topic has reattracted some interest recently, especially after 2004 when Stevo Stević used a transformation for showing the solvability of the nonlinear difference equation:
x n + 1 = x n 1 a + b x n 1 x n , n N 0 ,
(for some extensions of the results, see [11,12], while the corresponding results for some related systems of difference equations can be found in [13]). Since that time, related transformations have been frequently used on difference equations ([14,15,16]), as well as on close to symmetric systems (see [15,17,18] and numerous references therein), an area essentially initiated by Papaschinopoulos and Schinas (see [19,20,21,22,23,24,25]). Somewhat more complex methods can be found in [26]. For some related topics, such as finding invariants, special types of solutions and applications of solvable difference equations (see, for example, [21,22,23,25,27,28,29] and the references therein).
Quite frequently, the solvability was essentially shown by using some special cases of the following difference equation:
x n + 1 = a n x n + b n , n N ,
where coefficients ( a n ) n N , ( b n ) n N , and the initial value x 0 are real or complex (see [11,12,13,14,15,17,18,26,30]), which shows how useful the equation is (for how Equation (4) is solved, consult, for example, [3,8]; the book [3] has a nice explanation of three methods for solving it). Recently, we have studied several classes of product-type equations and systems (see [31,32,33,34] and the references therein), which cannot be directly solved by Equation (4), but behind their solvability is hidden the equation. Namely, these papers on product-type equations and systems use the solvability of some special cases of the following product-type analog of Equation (4):
z n = b n z n 1 a n , n N 0
(note that, if a n , b n , z n are positive, by taking the logarithm, the equation is transformed to Equation (4), which also shows the importance of the equation).
It has been noticed in [7] that Equation (4) can be useful not only in solving some ordinary difference equations, but also in solving some important equations with two independent variables, which has motivated us to conduct further investigations in this direction ([35,36]). In fact, Refs. [7,35,36] are devoted to the solvability of some boundary value problems. This paper is a continuation of our investigation of the solvability of difference equations, as well as the solvability of some boundary value problems for difference equations with two independent variables. It should be pointed out that not all of the coefficients of the equations studied here will be constant, which makes the study more complex. Note also that, for the special case a n = 1 of Equation (4), which is used in [7], is a simpler case, but the behavior of the solutions to Equation (4) can be quite complex (see [37,38,39] for the case of metric spaces). Let x ( 0 ) : = 1 and:
x ( n ) : = j = 0 n 1 ( x j ) ,
for n N , the so called falling factorial.
It is easy to see that x ( n ) is a polynomial of degree n, so that the following equality holds
x ( n ) = k = 0 n S k n x k ,
for some integer numbers S k n , k , n N 0 . These numbers S k n are called (signed) Stirling’s numbers of the first kind (see, for example, [3,9]).
Now, we will present some basic facts on the Stirling numbers. First, from Equality (5), we immediately see that:
S 0 n = 0 , n N ,
S n n = 1 , n N 0 .
From the definition of x ( n ) , Equality (5) and by using Viète formulas, it is easy to see that:
S k n = ( 1 ) n k 1 i 1 < i 2 < < i n k n 1 i 1 i 2 i n k ,
where i j { 1 , 2 , , n 1 } , j = 1 , n k ¯ , and 1 k n 1 , from which it follows that:
S 1 n = ( 1 ) n 1 ( n 1 ) ! , n N ,
S 2 n = ( 1 ) n 2 ( n 1 ) ! j = 1 n 1 1 j = ( 1 ) n 2 ( n 1 ) ! H n 1 , n 2 ,
S n 1 n = 1 j n 1 j = C 2 n , n 2 .
From (5), since:
x ( n + 1 ) = x ( n ) ( x n ) ,
and by some calculations, we also have that:
x ( n + 1 ) = k = 0 n + 1 S k n + 1 x k = ( x n ) k = 0 n S k n x k = k = 0 n S k n x k + 1 k = 0 n n S k n x k = S n n x n + 1 + k = 1 n ( S k 1 n n S k n ) x k n S 0 n ,
from which by comparing the coefficients of x k are obtained the equations:
S 0 n + 1 = n S 0 n , n N ,
S k n + 1 = S k 1 n n S k n , 1 k n , n N ,
S n + 1 n + 1 = S n n , n N .
From (12)–(14), Formulas (6), (7), (9)–(11) can also be obtained. Namely, from (12) and since S 0 1 = 0 , we again obtain (6), whereas from (14) and since S 1 1 = 1 , we again obtain (7). From (13) with k = 1 and by using (6),
S 1 n = ( n 1 ) S 1 n 1 ,
is obtained, from which, along with S 1 1 = 1 , it follows again that (9) holds. From (13) with k = n and by using (7), is obtained:
S n 1 n = S n 2 n 1 ( n 1 ) ,
from which by summing up this equality from 1 to n, and by using the fact S 0 1 = 0 , it again easily follows that (11) holds. Formula (10) can also be found, essentially reducing recurrent relation (13) with k = 2 , to a special case of Equation (4). Namely, using (9) into (13), we have that:
S 2 i = ( i 1 ) S 2 i 1 + S 1 i 1 = ( i 1 ) S 2 i 1 + ( 1 ) i 2 ( i 2 ) ! ,
for 3 i n .
Multiplying (15) by j = i n 1 ( j ) , summing up such obtained equalities for i = 3 , n ¯ , and using the fact that S 2 2 = 1 , we obtain:
S 2 n = S 2 2 j = 2 n 1 ( j ) + i = 3 n 1 ( i 1 ) j = 1 n 1 ( j ) = ( 1 ) n 2 ( n 1 ) ! H n 1 .
The formula can be found in [9], but the author did not notice its connection with Equation (4).
This process can be continued for other values of k. It is a natural problem to try to find a general formula for solutions to Equation (13), including the (special) solution S k n . As far as we know, there are no formulas of this type in the literature and our task will be to present some.
If we use the change of variables:
d n , k = ( 1 ) n k S k n ,
then the recurrent relation that is obtained from (13) when n n 1 , is transformed to the following one:
d n , k = d n 1 , k 1 + ( n 1 ) d n 1 , k ,
where 1 k < n , k , n N , and, from (8) and (16), it is clear that sequence | S k n | is a solution to Equation (17), which is a partial difference equation (for some classical results on the equations, mostly connected to the their solvability, see, for example, [9,10], for some results up to 2003, see [40]; see also an interesting paper [41]). The numbers | S k n | , 0 k n , are called the unsigned Stirling’s numbers of the first kind. To reduce working with expressions with negative signs to a great extent, from now on, we will consider Equation (17) instead of the equation obtained from (13) when n n 1 . Further details on the Stirling numbers can be found, for example, in [3,5,6,9].
From (17), we see that the Stirling numbers are a solution to a two-dimensional recurrent relation with given boundary conditions S n n and S 0 n , n N . Since there is a closed form formula for it, that is, a formula for the solution to Equation (17) with the conditions:
d n , 0 = 0 , d n , n = 1 ,
n N , it is expected that the equation can be also solved on the following domain:
C = { ( n , k ) : 0 k n , k N 0 , n N } ,
the so called combinatorial domain ([7]), for any values of d n , 0 and d n , n , which is the case with the so called Binomial partial difference equation (see [7]). Our analysis has shown that the methods presented in [9,10] do not seem suitable for solving Equation (17) on domain C , although we do not exclude a possibility that some modification of the methods could be suitable.
Using the idea and method in [7], we show, among other things, that there are closed form formulas for solutions to Equation (17) on domain C in terms of given boundary values d n , 0 and d n , n , n N .

2. Main Results

In this section, we show, among other things, how general solutions to the boundary-value problems for Equation (17) on domain C can be found by using the main idea and method in [7], which has been also recently employed in [35] for the case of another partial difference equation. Namely, domain C is divided into naturally chosen half-lines, and Equation (17) is regarded on each of these lines as an equation of type (4), that is, as a one-dimensional difference equation. Such one-dimensional equations are “solved" on the half lines, and, based on the obtained formulas, the general solution to the boundary-value problems for Equation (17) is obtained. In order to solve the equation, the method requires certain modifications, which are done/explained in the procedure that follows. There are also two natural ways of how to slice the domain C by lines. Based on the choice of slicing, some different looking solutions of the same problem are obtained. Besides Equation (17), an extension of it will be considered by using a different slicing, as well as another equation with a coefficient appearing in a different place.

2.1. Solving Equation (17)

To solve Equation (17), we will use the slicing of domain C by the lines:
k = c ,
where c N .
To be more suggestive, we will present the first three steps in detail and then give the general solution of the boundary-value problem on the domain. Let k = 1 , and then (17) becomes:
d n , 1 = ( n 1 ) d n 1 , 1 + d n 1 , 0 ,
for n 2 .
Multiplying the relation:
d i , 1 = ( i 1 ) d i 1 , 1 + d i 1 , 0 ,
by j = i n 1 j , summing up such obtained equalities for 2 i n and canceling the equal terms appearing on both sides of the summed equality, we get:
d n , 1 = ( n 1 ) ! d 1 , 1 + i = 1 n 1 d i , 0 j = i + 1 n 1 j = ( n 1 ) ! d 1 , 1 + i = 1 n 1 d i , 0 i ! ,
for n 2 . In fact, by using above-mentioned conventions, we see that (19) also holds for n = 1 .
Letting k = 2 , then (17) becomes:
d n , 2 = ( n 1 ) d n 1 , 2 + d n 1 , 1 ,
for n 3 .
Multiplying the relation:
d i , 2 = ( i 1 ) d i 1 , 2 + d i 1 , 1 ,
by j = i n 1 j , summing up such obtained equalities for 3 i n and canceling the equal terms appearing on both sides of the summed equality, we get:
d n , 2 = ( n 1 ) ! 1 ! d 2 , 2 + i = 2 n 1 d i , 1 j = i + 1 n 1 j = ( n 1 ) ! d 2 , 2 1 ! + i = 2 n 1 d i , 1 i ! ,
for n 3 . In fact, (21) also holds for n = 2 .
By using (19) into (21), we get:
d n , 2 = ( n 1 ) ! d 2 , 2 1 ! + i = 2 n 1 d i , 1 i ! = ( n 1 ) ! d 2 , 2 1 ! + i = 2 n 1 1 i ! ( i 1 ) ! d 1 , 1 + j = 1 i 1 d j , 0 j ! = ( n 1 ) ! d 2 , 2 1 ! + d 1 , 1 i = 2 n 1 1 i + i = 2 n 1 1 i j = 1 i 1 d j , 0 j ! ,
for n 3 .
Letting k = 3 , then (17) becomes:
d n , 3 = ( n 1 ) d n 1 , 3 + d n 1 , 2 ,
for n 4 .
Multiplying the relation:
d i , 3 = ( i 1 ) d i 1 , 3 + d i 1 , 2 ,
by j = i n 1 j , summing up such obtained equalities for all 4 i n and canceling the equal terms appearing on both sides of such obtained equality, it follows that:
d n , 3 = ( n 1 ) ! 2 ! d 3 , 3 + i = 3 n 1 d i , 2 j = i + 1 n 1 j = ( n 1 ) ! d 3 , 3 2 ! + i = 3 n 1 d i , 2 i ! ,
for every n 4 (in fact, it is immediately seen that the equality (24) holds also for n = 3 ).
By using (22) into (24), we get:
d n , 3 = ( n 1 ) ! d 3 , 3 2 ! + i = 3 n 1 d i , 2 i ! = ( n 1 ) ! d 3 , 3 2 ! + i = 3 n 1 1 i ! ( i 1 ) ! d 2 , 2 1 ! + d 1 , 1 j = 2 i 1 1 j + j = 2 i 1 1 j l = 1 j 1 d l , 0 l ! = ( n 1 ) ! d 3 , 3 2 ! + d 2 , 2 1 ! i = 3 n 1 1 i + d 1 , 1 i = 3 n 1 1 i j = 2 i 1 1 j + i = 3 n 1 1 i j = 2 i 1 1 j l = 1 j 1 d l , 0 l ! ,
for n 3 .
Assume that, for a fixed k such that k n , we have proved that:
d n , k = ( n 1 ) ! ( d k , k ( k 1 ) ! + d k 1 , k 1 ( k 2 ) ! i k = k n 1 1 i k + + d 1 , 1 i k = k n 1 1 i k i 3 = 3 i 4 1 1 i 3 i 2 = 2 i 3 1 1 i 2 + i k = k n 1 1 i k i 2 = 2 i 3 1 1 i 2 i 1 = 1 i 2 1 d i 1 , 0 i 1 ! ) ,
for n k .
If, in (17), we replace k by k + 1 , we get:
d n , k + 1 = ( n 1 ) d n 1 , k + 1 + d n 1 , k ,
for n k + 2 .
Multiplying the relation:
d i , k + 1 = ( i 1 ) d i 1 , k + 1 + d i 1 , k ,
by j = i n 1 j , summing up such obtained equalities for all k + 2 i n and canceling the equal terms appearing on both sides of the such obtained equality, it follows that:
d n , k + 1 = ( n 1 ) ! k ! d k + 1 , k + 1 + i k + 1 = k + 1 n 1 d i k + 1 , k j = i k + 1 + 1 n 1 j = ( n 1 ) ! d k + 1 , k + 1 k ! + i k + 1 = k + 1 n 1 d i k + 1 , k i k + 1 ! ,
for n k + 2 (in fact, for n k + 1 ).
Using (26) into (28), we get:
d n , k + 1 = ( n 1 ) ! d k + 1 , k + 1 k ! + i k + 1 = k + 1 n 1 d i k + 1 , k i k + 1 ! = ( n 1 ) ! ( d k + 1 , k + 1 k ! + i k + 1 = k + 1 n 1 1 i k + 1 ( d k , k ( k 1 ) ! + d k 1 , k 1 ( k 2 ) ! i k = k i k + 1 1 1 i k + + d 1 , 1 i k = k i k + 1 1 1 i k i 2 = 2 i 3 1 1 i 2 + i k = k i k + 1 1 1 i k i 2 = 2 i 3 1 1 i 2 i 1 = 1 i 2 1 d i 1 , 0 i 1 ! ) ) = ( n 1 ) ! ( d k + 1 , k + 1 k ! + d k , k ( k 1 ) ! i k + 1 = k + 1 n 1 1 i k + 1 + d k 1 , k 1 ( k 2 ) ! i k + 1 = k + 1 n 1 1 i k + 1 i k = k i k + 1 1 1 i k + + d 1 , 1 i k + 1 = k + 1 n 1 1 i k + 1 i k = k i k + 1 1 1 i k i 2 = 2 i 3 1 1 i 2 + i k + 1 = k + 1 n 1 1 i k + 1 i k = k i k + 1 1 1 i k i 2 = 2 i 3 1 1 i 2 i 1 = 1 i 2 1 d i 1 , 0 i 1 ! ) ,
for n k + 2 . Note that (29) also holds for n = k + 1 .
From (19), (29) and the method of induction, it follows that Formula (26) holds for n k .
From the above conducted detailed analysis and consideration, we see that the following result holds.
Theorem 1.
If ( a k ) k N , ( b k ) k N are given sequences of complex numbers. Then, the solution to the partial difference Equation (17) on the domain C , with the boundary value conditions,
d k , 0 = a k a n d d k , k = b k , k N ,
is given by:
d n , k = ( n 1 ) ! j = 1 k b j ( j 1 ) ! i k = k n 1 1 i k i j + 1 = 2 i j + 2 1 1 i j + 1 + i k = k n 1 1 i k i 2 = 2 i 3 1 1 i 2 i 1 = 1 i 2 1 a i 1 i 1 ! .
Proof. 
Using all the conditions given in (30) in Formula (26), Formula (31) is obtained. ☐
From Theorem 1 and (16), the following corollary follows.
Corollary 1.
If ( a k ) k N , ( b k ) k N are given sequences of real numbers, then the solution to the partial difference equation:
s n , k = s n 1 , k 1 ( n 1 ) s n 1 , k ,
on the domain C , with the boundary value conditions:
s k , 0 = a k a n d s k , k = b k , k N ,
is given by:
s n , k = ( 1 ) n k ( n 1 ) ! j = 1 k b j ( j 1 ) ! i k = k n 1 1 i k i j + 1 = 2 i j + 2 1 1 i j + 1 + i k = k n 1 1 i k i 2 = 2 i 3 1 1 i 2 i 1 = 1 i 2 1 a i 1 i 1 ! .
Remark 1.
Note that, in the first three steps above, as well as in the main part of the inductive argument, we essentially used a method for solving Equation (4) (see, how we "solve" Equations (18), (20), (23) and (27)).

2.2. Partial Difference Equations with an Interchanged Coefficient

In this section, we consider the following partial difference equation:
d n , k = ( n 1 ) d n 1 , k 1 + d n 1 , k ,
on domain C .
We will solve the equation also by slicing domain C by the lines k = c , c N and use the above procedure for solving the corresponding one-dimensional difference equations on such obtained half-lines that will lead to the general solution to the boundary-value problem for Equation (35) on the domain.
Assuming first that k = 1 , then (35) becomes:
d n , 1 = ( n 1 ) d n 1 , 0 + d n 1 , 1 ,
for n 2 .
Summing up the following relations:
d i , 1 = ( i 1 ) d i 1 , 0 + d i 1 , 1 ,
for 2 i n , we get:
d n , 1 = d 1 , 1 + i = 1 n 1 i d i , 0 ,
for n 2 . Note that (37) also holds for n = 1 .
Letting k = 2 , then (35) becomes:
d n , 2 = ( n 1 ) d n 1 , 1 + d n 1 , 2 ,
for n 3 .
Summing up the relations:
d i , 2 = ( i 1 ) d i 1 , 1 + d i 1 , 2 ,
for 3 i n , canceling the equal terms and using Formula (37), we get:
d n , 2 = d 2 , 2 + j = 2 n 1 j d j , 1 = d 2 , 2 + j = 2 n 1 j d 1 , 1 + i = 1 j 1 i d i , 0 = d 2 , 2 + d 1 , 1 j = 2 n 1 j + j = 2 n 1 j i = 1 j 1 i d i , 0 ,
for n 3 . Note that (39) also holds for n = 2 .
Assume that the following formula has been proved for a k N and n k ,
d n , k = d k , k + d k 1 , k 1 i k = k n 1 i k + d k 2 , k 2 i k = k n 1 i k i k 1 = k 1 i k 1 i k 1 + + d 1 , 1 i k = k n 1 i k i 3 = 3 i 4 1 i 3 i 2 = 2 i 3 1 i 2 + i k = k n 1 i k i 3 = 3 i 4 1 i 3 i 2 = 2 i 3 1 i 2 i 1 = 1 i 2 1 i 1 d i 1 , 0 .
Then, by using (35) with k k + 1 , we get:
d n , k + 1 = ( n 1 ) d n 1 , k + d n 1 , k + 1 ,
for n k + 2 .
Summing up the relations:
d i , k + 1 = ( i 1 ) d i 1 , k + d i 1 , k + 1 ,
for k + 2 i n , and using (40), it follows that:
d n , k + 1 = d k + 1 , k + 1 + i k + 1 = k + 1 n 1 i k + 1 d i k + 1 , k = d k + 1 , k + 1 + i k + 1 = k + 1 n 1 i k + 1 ( d k , k + d k 1 , k 1 i k = k i k + 1 1 i k + d k 2 , k 2 i k = k i k + 1 1 i k i k 1 = k 1 i k 1 i k 1 + + d 1 , 1 i k = k i k + 1 1 i k i 3 = 3 i 4 1 i 3 i 2 = 2 i 3 1 i 2 + i k = k i k + 1 1 i k i 3 = 3 i 4 1 i 3 i 2 = 2 i 3 1 i 2 i 1 = 1 i 2 1 i 1 d i 1 , 0 ) = d k + 1 , k + 1 + d k , k i k + 1 = k + 1 n 1 i k + 1 + d k 1 , k 1 i k + 1 = k + 1 n 1 i k + 1 i k = k i k + 1 1 i k + + d 1 , 1 i k + 1 = k + 1 n 1 i k + 1 i k = k i k + 1 1 i k i 3 = 3 i 4 1 i 3 i 2 = 2 i 3 1 i 2 + i k + 1 = k + 1 n 1 i k + 1 i k = k i k + 1 1 i k i 3 = 3 i 4 1 i 3 i 2 = 2 i 3 1 i 2 i 1 = 1 i 2 1 i 1 d i 1 , 0 ,
for n k + 2 ((42) obviously holds for n = k + 1 ).
From (37), (42) and the method of mathematical induction, it follows that Formula (40) really holds for n k .
Now, note that the procedure described above can be applied to every partial difference equation in the following form:
d n , k = f ( n 1 ) d n 1 , k 1 + d n 1 , k ,
on domain C , where ( f ( m ) ) m N is a complex sequence. The following theorem holds (the proof is omitted since it is similar to the one given in the case f ( m ) = m ).
Theorem 2.
Let ( u j ) j N , ( v j ) j N be given sequences of complex numbers. Then, the solution to Equation (43) on the domain C , with the following conditions
d j , 0 = u j a n d d j , j = v j , j N ,
is given by:
d n , k = v k + v k 1 i k = k n 1 f ( i k ) + v k 2 i k = k n 1 f ( i k ) i k 1 = k 1 i k 1 f ( i k 1 ) + + v 1 i k = k n 1 f ( i k ) i 3 = 3 i 4 1 f ( i 3 ) i 2 = 2 i 3 1 f ( i 2 ) + i k = k n 1 f ( i k ) i 3 = 3 i 4 1 f ( i 3 ) i 2 = 2 i 3 1 f ( i 2 ) i 1 = 1 i 2 1 f ( i 1 ) u i 1 ,
for n k .
Theorem 2 shows that the calculation of iterated sums plays an important role in solving Equation (43).

2.3. An Extension to Equation (17)

Here, we consider the following class of partial difference equations:
d n , k = d n 1 , k 1 + f ( n , k ) d n 1 , k ,
where ( n , k ) C , which is an extension to Equation (17), and whose special cases appear in the literature.
In contrast to the previous two equations, this time, we will find a solution to the equation by slicing domain C in a different way. Namely, the domain will be sliced now by the lines:
y l ( k ) = k + l , k N ,
where l N .
For n = k + 1 , Equation (45) becomes:
d k + 1 , k = d k , k 1 + f ( k + 1 , k ) d k , k ,
for k N .
Summing up equalities (46) from 1 to k,
d k + 1 , k = d 1 , 0 + j = 1 k f ( j + 1 , j ) d j , j ,
is obtained for k N . Note that (47) holds also for k = 0 .
For n = k + 2 , Equation (45) becomes:
d k + 2 , k = d k + 1 , k 1 + f ( k + 2 , k ) d k + 1 , k ,
for k N .
Summing up equalities (48) from 1 to k, and using (47) in such obtained equality, we get:
d k + 2 , k = d 2 , 0 + j = 1 k f ( j + 2 , j ) d j + 1 , j = d 2 , 0 + j = 1 k f ( j + 2 , j ) d 1 , 0 + i = 1 j f ( i + 1 , i ) d i , i = d 2 , 0 + d 1 , 0 j = 1 k f ( j + 2 , j ) + j = 1 k f ( j + 2 , j ) i = 1 j f ( i + 1 , i ) d i , i ,
for k N . Note that (49) holds also for k = 0 .
Assume that, for an l N ,
d k + l , k = d l , 0 + d l 1 , 0 i l = 1 k f ( i l + l , i l ) + d l 2 , 0 i l = 1 k f ( i l + l , i l ) i l 1 = 1 i l f ( i l 1 + l 1 , i l 1 ) + + d 1 , 0 i l = 1 k f ( i l + l , i l ) i l 1 = 1 i l f ( i l 1 + l 1 , i l 1 ) i 2 = 1 i 3 f ( i 2 + 2 , i 2 ) + i l = 1 k f ( i l + l , i l ) i l 1 = 1 i l f ( i l 1 + l 1 , i l 1 ) i 1 = 1 i 2 f ( i 1 + 1 , i 1 ) d i 1 , i 1 ,
for k N 0 .
For n = k + l + 1 , Equation (45) becomes:
d k + l + 1 , k = d k + l , k 1 + f ( k + l + 1 , k ) d k + l , k ,
for k N .
Summing up equalities (51) from 1 to k, and using (50) in such obtained equalities, we get:
d k + l + 1 , k = d l + 1 , 0 + i l + 1 = 1 k f ( i l + 1 + l + 1 , i l + 1 ) d i l + 1 + l , i l + 1 = d l + 1 , 0 + i l + 1 = 1 k f ( i l + 1 + l + 1 , i l + 1 ) ( d l , 0 + d l 1 , 0 i l = 1 i l + 1 f ( i l + l , i l ) + + d 1 , 0 i l = 1 i l + 1 f ( i l + l , i l ) i l 1 = 1 i l f ( i l 1 + l 1 , i l 1 ) i 2 = 1 i 3 f ( i 2 + 2 , i 2 ) + i l = 1 i l + 1 f ( i l + l , i l ) i l 1 = 1 i l f ( i l 1 + l 1 , i l 1 ) i 1 = 1 i 2 f ( i 1 + 1 , i 1 ) d i 1 , i 1 ) = d l + 1 , 0 + d l , 0 i l + 1 = 1 k f ( i l + 1 + l + 1 , i l + 1 ) + d l 1 , 0 i l + 1 = 1 k f ( i l + 1 + l + 1 , i l + 1 ) i l = 1 i l + 1 f ( i l + l , i l ) + + d 1 , 0 i l + 1 = 1 k f ( i l + 1 + l + 1 , i l + 1 ) i l = 1 i l + 1 f ( i l + l , i l ) i 2 = 1 i 3 f ( i 2 + 2 , i 2 ) + i l + 1 = 1 k f ( i l + 1 + l + 1 , i l + 1 ) i l = 1 i l + 1 f ( i l + l , i l ) i 1 = 1 i 2 f ( i 1 + 1 , i 1 ) d i 1 , i 1 ,
for k N . Since (52) obviously holds for k = 0 , it holds for every k N 0 .
From (47), (52) and the method of induction, we see that (50) holds for every k N 0 and l N .
From the above considerations, we have the following results.
Theorem 3.
Let ( u j ) j N , ( v j ) j N be given sequences of complex numbers. Then, the solution to Equation (45) on the domain C , with the following boundary-value conditions:
d j , 0 = u j a n d d j , j = v j , j N ,
is given by:
d n , k = u n k + u n k 1 i n k = 1 k f ( i n k + n k , i n k ) + u n k 2 i n k = 1 k f ( i n k + n k , i n k ) i n k 1 = 1 i n k f ( i n k 1 + n k 1 , i n k 1 ) + + u 1 i n k = 1 k f ( i n k + n k , i n k ) i 3 = 1 i 4 f ( i 3 + 3 , i 3 ) i 2 = 1 i 3 f ( i 2 + 2 , i 2 ) + i n k = 1 k f ( i n k + n k , i n k ) i 2 = 1 i 3 f ( i 2 + 2 , i 2 ) i 1 = 1 i 2 f ( i 1 + 1 , i 1 ) v i 1 ,
for k , n N such that n k + 1 .
Proof. 
If we apply the change of variables l = n k in Formula (50) and then use the boundary value conditions (53), we obtain (54). ☐
If we apply Theorem 3 for the case f ( n , k ) = n 1 , we obtain another formula for the solution to the boundary value problem in Theorem 1. Namely, the following result holds.
Corollary 2.
Let ( u j ) j N , ( v j ) j N be given sequences of complex numbers. Then the solution to Equation (17) on the domain C , with the following boundary-value conditions:
d j , 0 = u j a n d d j , j = v j , j N ,
is given by:
d n , k = u n k + u n k 1 i n k = 1 k ( i n k + n k 1 ) + u n k 2 i n k = 1 k ( i n k + n k 1 ) i n k 1 = 1 i n k ( i n k + n k 2 ) + + u 1 i n k = 1 k ( i n k + n k 1 ) i 3 = 1 i 4 ( i 3 + 2 ) i 2 = 1 i 3 ( i 2 + 1 ) + i n k = 1 k ( i n k + n k 1 ) i 2 = 1 i 3 ( i 2 + 1 ) i 1 = 1 i 2 i 1 v i 1 ,
for k , n N such that n k + 1 .
If, in Corollary 2, we choose d j , 0 = 0 and d j , j = 1 , for j N , and use (16), we get the following formula for the Stirling numbers of the first kind:
S k n = ( 1 ) n k i n k = 1 k ( i n k + n k 1 ) i 2 = 1 i 3 ( i 2 + 1 ) i 1 = 1 i 2 i 1 ,
or by changing the order of the summation:
S k n = ( 1 ) n k i 1 = 1 k i 1 i 2 = i 1 k ( i 2 + 1 ) i n k = i n k 1 k ( i n k + n k 1 ) ,
for k , n N , such that n k + 1 .

3. Conclusions

Here, we have presented another influence and application of linear first-order difference equations, this time on the solvability of some boundary value problems for some important classes of recurrent relations with two independent variables on the so-called combinatorial domain. Our method of half-lines is used and suitably modified to show the solvability of the problems. To a certain extent, we have been motivated by some recurrent relations in combinatorics, particularly by the recurrent relation generating Stirling’s numbers of the first kind. Using this fact, as well as the fact that the recurrent relation for the binomial coefficients can be solved by the method of half-lines, it is highly expected that some suitable modifications of the method can be successfully used in considering the other related recurrent relations that stem from combinatorics, as well as the other areas of mathematics and science, which will be one of our topics for further investigation.

Acknowledgments

The work of Bratislav Iričanin was supported by the Serbian Ministry of Education and Science projects III 41025 and OI 171007. The work of Zdeněk Šmarda was supported by the project FEKT-S-17-4225 of Brno University of Technology.

Author Contributions

The paper is based on some ideas by Stevo Stević. All of the authors have contributed to the writing of this paper. They read and approved the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Krechmar, V.A. A Problem Book in Algebra; Mir Publishers: Moscow, Russia, 1974. [Google Scholar]
  2. Mitrinović, D.S. Mathematical Induction, Binomial Formula, Combinatorics; Gradjevinska Knjiga: Beograd, Serbia, 1980. (In Serbian) [Google Scholar]
  3. Mitrinović, D.S.; Kečkić, J.D. Methods for Calculating Finite Sums; Naučna Knjiga: Beograd, Serbia, 1984. (In Serbian) [Google Scholar]
  4. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series; Nauka: Moscow, Russian, 1981. (In Russian) [Google Scholar]
  5. Riordan, J. Combinatorial Identities; John Wiley & Sons Inc.: New York, NY, USA; London, UK; Sydney, Australia, 1968. [Google Scholar]
  6. Yablonskiy, S.V. Introduction to Discrete Mathematics; Mir Publishers: Moscow, Russia, 1989. [Google Scholar]
  7. Stević, S. Note on the binomial partial difference equation. Electron. J. Qual. Theory Differ. Equ. 2015, 2015. Article 96. [Google Scholar]
  8. Brand, L. Differential and Difference Equations; John Wiley & Sons, Inc.: New York, NY, USA, 1966. [Google Scholar]
  9. Jordan, C. Calculus of Finite Differences; Chelsea Publishing Company: New York, NY, USA, 1956. [Google Scholar]
  10. Levy, H.; Lessman, F. Finite Difference Equations; Dover Publications: New York, NY, USA, 1992. [Google Scholar]
  11. Stević, S. On the difference equation xn=xnk/(b+cxn−1xnk). Appl. Math. Comput. 2012, 218, 6291–6296. [Google Scholar]
  12. Stević, S.; Diblik, J.; Iričanin, B.; Šmarda, Z. On the difference equation xn=anxnk/(bn+cnxn−1xnk). Abstr. Appl. Anal. 2012, 2012. Article 409237. [Google Scholar]
  13. Stević, S.; Diblik, J.; Iričanin, B.; Šmarda, Z. On a third-order system of difference equations with variable coefficients. Abstr. Appl. Anal. 2012, 2012. Article 508523. [Google Scholar]
  14. Papaschinopoulos, G.; Stefanidou, G. Asymptotic behavior of the solutions of a class of rational difference equations. Int. J. Differ. Equ. 2010, 5, 233–249. [Google Scholar]
  15. Stević, S.; Diblik, J.; Iričanin, B.; Šmarda, Z. On some solvable difference equations and systems of difference equations. Abstr. Appl. Anal. 2012, 2012. Article 541761. [Google Scholar]
  16. Stević, S.; Diblik, J.; Iričanin, B.; Šmarda, Z. Solvability of nonlinear difference equations of fourth order. Electron. J. Differ. Equ. 2014, 2014. Article 264. [Google Scholar]
  17. Berg, L.; Stević, S. On some systems of difference equations. Appl. Math. Comput. 2011, 218, 1713–1718. [Google Scholar] [CrossRef]
  18. Stević, S.; Iričanin, B.; Šmarda, Z. On a close to symmetric system of difference equations of second order. Adv. Differ. Equ. 2015, 2015. Article 264. [Google Scholar]
  19. Fotiades, N.; Papaschinopoulos, G. On a system of difference equations with maximum. Appl. Math. Comput. 2013, 221, 684–690. [Google Scholar]
  20. Papaschinopoulos, G.; Schinas, C.J. On a system of two nonlinear difference equations. J. Math. Anal. Appl. 1998, 219, 415–426. [Google Scholar] [CrossRef]
  21. Papaschinopoulos, G.; Schinas, C.J. On the behavior of the solutions of a system of two nonlinear difference equations. Commun. Appl. Nonlinear Anal. 1998, 5, 47–59. [Google Scholar]
  22. Papaschinopoulos, G.; Schinas, C.J. Invariants for systems of two nonlinear difference equations. Differ. Equ. Dyn. Syst. 1999, 7, 181–196. [Google Scholar] [CrossRef]
  23. Papaschinopoulos, G.; Schinas, C.J. Invariants and oscillation for systems of two nonlinear difference equations. Nonlinear Anal. Theory Methods Appl. 2001, 46, 967–978. [Google Scholar] [CrossRef]
  24. Papaschinopoulos, G.; Schinas, C.J. On the dynamics of two exponential type systems of difference equations. Comput. Math. Appl. 2012, 64, 2326–2334. [Google Scholar] [CrossRef]
  25. Papaschinopoulos, G.; Schinas, C.J.; Stefanidou, G. On a k-order system of Lyness-type difference equations. Adv. Differ. Equ. 2007, 2007. Article 31272. [Google Scholar] [CrossRef]
  26. Stević, S.; Diblik, J.; Iričanin, B.; Šmarda, Z. On a solvable system of rational difference equations. J. Differ. Equ. Appl. 2014, 20, 811–825. [Google Scholar] [CrossRef]
  27. Berezansky, L.; Braverman, E. On impulsive Beverton-Holt difference equations and their applications. J. Differ. Equ. Appl. 2004, 10, 851–868. [Google Scholar] [CrossRef]
  28. Diblik, J.; Halfarová, H. Explicit general solution of planar linear discrete systems with constant coefficients and weak delays. Adv. Differ. Equ. 2013, 2013. Article 50. [Google Scholar] [CrossRef]
  29. Iričanin, B.; Stević, S. Eventually constant solutions of a rational difference equation. Appl. Math. Comput. 2009, 215, 854–856. [Google Scholar] [CrossRef]
  30. Stević, S.; Diblik, J.; Iričanin, B.; Šmarda, Z. On the difference equation xn+1=xnxnk/(xnk+1(a+bxnxnk)). Abstr. Appl. Anal. 2012, 2012. Article 108047. [Google Scholar]
  31. Stević, S. First-order product-type systems of difference equations solvable in closed form. Electron. J. Differ. Equ. 2015, 2015. Article 308. [Google Scholar]
  32. Stević, S.; Iričanin, B.; Šmarda, Z. On a product-type system of difference equations of second order solvable in closed form. J. Inequal. Appl. 2015, 2015. Article 327. [Google Scholar]
  33. Stević, S.; Iričanin, B.; Šmarda, Z. Solvability of a close to symmetric system of difference equations. Electron. J. Differ. Equ. 2016, 2016. Article 159. [Google Scholar]
  34. Stević, S.; Iričanin, B.; Šmarda, Z. Two-dimensional product-type system of difference equations solvable in closed form. Adv. Differ. Equ. 2016, 2016. Article 253. [Google Scholar]
  35. Stević, S. Solvability of boundary value problems for a class of partial difference equations on the combinatorial domain. Adv. Differ. Equ. 2016, 2016. Article 262. [Google Scholar]
  36. Stević, S. Solvability of boundary-value problems for a linear partial difference equation. Electron. J. Differ. Equ. 2017, 2017. Article 17. [Google Scholar]
  37. Mitrinović, D.S.; Adamović, D.D. Sequences and Series; Naučna Knjiga: Beograd, Serbia, 1980. (In Serbian) [Google Scholar]
  38. Polya, G.; Szegö, G. Aufgaben Und Lehrsätze Aus Der Analysis; Verlag von Julius Springer: Berlin, Germany, 1925. [Google Scholar]
  39. Ašić, M.D.; Adamović, D.D. Limit points of sequences in metric spaces. Am. Math. Mon. 1970, 77, 613–616. [Google Scholar] [CrossRef]
  40. Cheng, S.S. Partial Difference Equations; Taylor & Francis: London, UK; New York, NY, USA, 2003. [Google Scholar]
  41. Musielak, A.; Popenda, J. On the hyperbolic partial difference equations and their oscillatory properties. Glas. Mat. 1998, 33, 209–221. [Google Scholar]

Share and Cite

MDPI and ACS Style

Stević, S.; Iričanin, B.; Šmarda, Z. Boundary Value Problems for Some Important Classes of Recurrent Relations with Two Independent Variables. Symmetry 2017, 9, 323. https://doi.org/10.3390/sym9120323

AMA Style

Stević S, Iričanin B, Šmarda Z. Boundary Value Problems for Some Important Classes of Recurrent Relations with Two Independent Variables. Symmetry. 2017; 9(12):323. https://doi.org/10.3390/sym9120323

Chicago/Turabian Style

Stević, Stevo, Bratislav Iričanin, and Zdeněk Šmarda. 2017. "Boundary Value Problems for Some Important Classes of Recurrent Relations with Two Independent Variables" Symmetry 9, no. 12: 323. https://doi.org/10.3390/sym9120323

APA Style

Stević, S., Iričanin, B., & Šmarda, Z. (2017). Boundary Value Problems for Some Important Classes of Recurrent Relations with Two Independent Variables. Symmetry, 9(12), 323. https://doi.org/10.3390/sym9120323

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