Jca 21 08
Jca 21 08
Jca 21 08
Classical
A nalysis
Volume 21, Number 2 (2023), 87–152 doi:10.7153/jca-2023-21-08
D IEGO D OMINICI
Abstract. We study recurrence relations satisfied by the moments n (z) of discrete linear func-
tionals whose first moment satisfies a holonomic differential equation. We consider all cases
when the order of the ODE is less or equal than 3 .
1. Introduction
= zz , (1)
where z is the derivative operator z = z. We use the notation
(x + c) = (x + c1 ) · · · (x + cm ) , c ∈ Km , m ∈ N,
(x + 0)
/ = x.
c , Zagreb 87
Paper JCA-21-08
88 D. D OMINICI
and say that {n }n0 ⊂ K[x] is a monic basis of K[x] if n is monic and deg (n ) = n
for all n ∈ N0 . The Pochhammer symbol (c)n is defined by [49, 5.2.4]
n−1
(c)0 = 1, (c)n = (c + j), n ∈ N, c ∈ K,
j=0
In this article, we continue the work started in [21], where we studied the moments
n = L p,q [n ] of linear functionals L p,q : K [x] → F (acting on the variable x) defined
by
(a)x zx
L p,q [u] = u (x) , u ∈ K[x], (2)
x=0 (b + 1)x x!
with a ∈ K p and b ∈ Kq . It follows from (2) that the first moment 0 (z) = L [1] satisfies
a differential equation (in the variable z) with polynomial coefficients p,q [y] = 0,
where the differential operator p,q is defined by
p,q = ( + b) − z ( + a),
Functionals of the form (2) are called discrete semiclassical [24], [43], [47]. If we
define the polynomials
2. Previous results
In this section, we give a brief description of the results from our previous work
that will be needed in this paper. For more details and proofs, we refer the reader to
[21].
Since the operator defined in (1) satisfies
[zx f ] = zx ( + x) [ f ] , f ∈ F,
and therefore
L [u] = u ( ) L [1], u ∈ K [x] . (7)
We conclude from (7) that the moments n (z) of L on any monic basis {n }n0 are
completely determined by the first moment
n−1
n (x) = (x − k ) , n ∈ N, 0 (x) = 1,
k=0
where {k }k0 ⊂ K is a fixed sequence. These polynomials satisfy the 2 -term recur-
rence relation
xn (x) = n+1 (x) + n n (x) . (9)
Setting x = in (9), multiplying by 0 (z) and using (8), we obtain.
it follows that
and the moments of L p,q on the falling factorial polynomials basis modified moments
Since
zn zn [zx ] = zn n (x) zx−n = n (x) zx ,
we have n (z) = zn zn [0 ] and using the formula [49, 16.3.1]
a (a)n a+n
zn p Fq ;z = p Fq ;z ,
b+1 (b + 1)n b+n+1
R ECURRENCE RELATIONS FOR THE MOMENTS 91
and using (10) with n = 0 in (17) we have p,q [0 ] = 0, where the standard moments
recurrence operator p,q is defined by
where the elementary symmetric polynomials en (c) are defined by the generating func-
tion [41]
m
en (c)t n = (1 + tci) , c ∈ Km .
n=0 i=1
Using (19), we can rewrite p,q in extended form
q n p
n
p,q = eq−k (b) S n+k+1 − z e p− j (a) S k+ j ,
k=0 k=0 k j=0
and the equation p,q [0 ] = 0 gives a recurrence for the standard moments
q n p
n
eq−k (b) n+k+1 − z k e p− j (a) k+ j = 0
k=0 k=0 j=0
and the equation p,q [n ] = 0 gives a recurrence for the modified moments
q k p k
em−k (b + n + 1) Sk,i n+i+1 − z e p−k (a + n) Sk,i n+i = 0,
k=0 i=0 k=0 i=0
of minimal order s + 1, where s is the class of the functional L p,q defined in (4).
2.2. Transformations
There are 4 canonical transformations of the functional L p,q :
1) The Christoffel transformation at , which we define by
(a) zx
LCp,q [u] = (x − )u (x) (b + 1)
x
, ∈ K, u ∈ K[x]. (22)
x=0 x x!
and in particular
nC = n+1 + (n − ) n . (23)
Since from (15) it follows that
(− + 1)x
x − = − , (24)
(− )x
Thus, 0C is a solution of the ODE Cp,q [y] = 0, where the differential operator Cp,q is
defined by
Cp,q = ( ) ( − − 1) − z ( ) ( − + 1), (25)
and the polynomials , where defined in (3).
R ECURRENCE RELATIONS FOR THE MOMENTS 93
From (18) and (25) we see that Cp,q 0C = 0, where the recurrence operator Cp,q
is defined by
Cp,q = S n+1 (S + b)(S − − 1) − z (S + 1)n (S + a)(S − + 1).
Using (24), we see that the first moment has the hypergeometric representation
G −1 a, −
0 (z; ) = − p+1 Fq+1 ;z ,
b + 1, − + 1
and therefore is a solution of the ODE Gp,q [y] = 0, where the differential operator Gp,q
is defined by
Gp,q = ( ) ( − ) − z ( ) ( − ) = p,q ( − ). (27)
From (9) and (26), we see that
G
n+1 + (n − ) nG = LGp,q [(x − )n ] = L p,q [n ] = n ,
and in particular
G
n+1 + (n − ) nG = n . (28)
From (18) and (27) we see that Gp,q 0G = 0, where the recurrence operator Gp,q
is defined by
Gp,q = p,q (S − ) .
Similarly, (20) and (27) give Gp,q nG = 0, with
Gp,q = p,q (S + n − ).
3) The Uvarov transformation at , which we define by
LUp,q [u] = L p,q [u] + u ( ) z , ∈ K, u ∈ K [x] .
and therefore
Up,q [0 ] = 0, Up,q [z ] = 0.
Comparing (29) with (25) and (27), we see that
as some authors do. Thus, we will not consider the double transformation U,G
p,q , since
we can simplify the composition
U,G G
p,q = p,q .
U( )
p,q = (S − − 1) p,q ,
and
U( )
p,q = (S + n − ) p,q .
The second possibility for a reduced-Uvarov transformation happens when ( ) =
0. We now have
U( )
p,q = ( − ) p,q ,
U( )
p,q = (S − ) p,q ,
R ECURRENCE RELATIONS FOR THE MOMENTS 95
and
U( )
p,q = (S + n − − 1) p,q .
N
(a) zx
LTp,q [u] = u (x) (b + 1)
x
, N ∈ N0 , u ∈ K[x].
x=0 x x!
and therefore is a solution of the ODE Tp,q [y] = 0, where the differential operator Tp,q
is defined by
From (18) and (30) we see that Tp,q 0T = 0, where the recurrence operator Tp,q
is defined by
Tp,q = (S − N − 1) p,q .
3. Examples
We now illustrate the application of the formulas that we have derived in the previ-
ous section. We will consider all discrete semiclassical functionals of class s 2, and
also look at the subclasses obtained by applying one or more of the canonical transfor-
mations of the functional L p,q .
1,0;N = − z ( − N) .
Modified moments
1,0;N = S − z (S + n − N).
0,1 = ( + b) − z.
zn −
n (z) = 0 F1 ;z
(b + 1)n b+1+n
n−b √
= z 2 (b + 1)Ib+n 2 z ,
where Iv (z) denotes the modified Bessel function of the first kind [49, 10.25.2].
Modified moments recurrence operator
R EMARK 4. If we write
n = An pn ,
then the recurrence (34) becomes
(n + b) z
pn+1 + pn − 2 pn−1 = 0.
A A
Choosing
1 z
= −2x, − = 1,
A A2
we get
pn+1 − 2 (n + b)xpn + pn−1 = 0. (35)
The orthogonal polynomials satisfying the 3-term recurrence relation (35) with initial
conditions
p0 = 1, p1 = 2bx
n = (−1)n qn ,
where the monic polynomials qn (b) satisfy the 3 -term recurrence relation
For additional references on the generalized Charlier polynomials, see [16], [35],
[54], [58].
R ECURRENCE RELATIONS FOR THE MOMENTS 99
(a)n a+n
n (z) = zn 1 F1 ;z .
(b + 1)n b+1+n
R EMARK 5. If we define
n = (−1)n pn ,
then the monic polynomials pn (b) satisfy the 3 -term recurrence relation
For additional references on the generalized Meixner polynomials, see [10], [16],
[29].
C0,0 = ( − − 1) − z ( − + 1),
a = − + 1, b = − − 1.
C0,0 = S (S + n − − 1) − z (S + n − + 1)
(38)
= S 2 + (n − − z) S − z (n − + 1).
R EMARK 6. Using (37), we see that the modified moments satisfy the first order
C [ ] = 0, with
recurrence 0,0
C
0,0 = (n − + z) S − z (n − + 1 + z).
G G
0,0 = 0,0 ( − ) = ( − ) − z ( − ),
a = − , b = − .
R ECURRENCE RELATIONS FOR THE MOMENTS 101
G
0,0 = 0,0 (S − ) = S
n+1
− z (S + 1)n (S − ).
we obtain
nG (z; ) = (−1)n (−z) (n − , −z).
Since the function (a, z) satisfies the recurrence [49, 8.8.1]
(a + 1, z) = a (a, z) − za e−z ,
we see that
G
n+1 + (n − ) nG = zn ez = n ,
in agreement with (28).
Modified moments recurrence operator
G 2
0,0 = 0,0 (S + n − ) = S + (n + 1 − − z) S − z (n − ).
( ) = ,
a = 0, b = −1.
First moment
N
zx
0T (z) = x!
= N (z) ,
x=0
where N (z) denotes the truncated exponential series [49, 8.4.11].
First moment differential operator
T0,0 = ( − N − 1)0,0 = ( − N − 1) − z ( − N) ,
a = −N, b = −N − 1.
U (a, a, z) = ez (1 − a, z),
we conclude that
zn ez
nT (z) = (N − n + 1, z).
(N − n)!
Comparing with (31), we see that
(N − n + 1, z)
nT (z) = n (z) , 0 n N.
(N − n)!
2,0;N = − z ( − N) ( + a).
n − N, a + n
n (z) = zn (−N)n (a)n 2 F0 ;z .
−
(0,0)
Using the hypergeometric representation (32) of the Charlier polynomials Pn (x; z) ,
we can write
(0,0)
n (z) = zn (−N)n (a)n PN−n −a − n; −z−1 .
Modified moments recurrence operator
2,0;N = S − z (S + n − N)(S + n + a)
= −zS 2 + [1 − z (2n + 1 − N + a)] S + z (N − n)(n + a).
R EMARK 7. If we set z−1 = x, we see that the modified moments are a family of
monic orthogonal polynomials pn (x) , satisfying the 3-term recurrence relation
R EMARK 8. If we set b = −x, we see that the modified moments are a family of
orthogonal polynomials pn (x) , satisfying the 3-term recurrence relation
xpn = (1 − z) pn+1 + [n − z (2n − 1 + a1 + a2 )] pn
+ z (n − 1 + a1) (n − 1 + a2) pn−1 .
a1 , a2 (b − a2)−a1
2 F1 ;1 = , Re (b − a2) > Re (a1 ) .
b (b)−a1
Choosing a1 = −N, N ∈ N, a2 = a, we get
(−N)n (a)n (b + 1 − a)N−n
n (1) = ,
(b + 1)n (b + 1 + n)N−n
or
(b + 1 − a)N (−N)n (a)n
n (1) = (−1)n ,
(b + 1)N (a − b − N)n
which is a solution of (41) when z = 1, a1 = −N, and a2 = a
[(b − a + N − n)S − (n − N)(n + a)] [n (1)] = 0.
For additional references, see [20], [30].
R ECURRENCE RELATIONS FOR THE MOMENTS 105
R EMARK 9. From (42), we see that the modified moments satisfy the first order
C [ ] = 0, where
recurrence 1,0
C
1,0 = (1 − z)(n − + z + az)S − z (n + a)(n + 1 − + z + az).
This agrees with the second order recurrence (43), since
C
(S + n + 1 − ) 1,0 = (n + 1 − + z + az)C1,0 .
106 D. D OMINICI
G
1,0 = 1,0 ( − ) = ( − ) − z ( + a)( − ),
a1 = a, a2 = − , b = − .
G
1,0 = 1,0 (S − ) = S
n+1
− z (S + 1)n (S + a) (S − ).
zn a + n, n −
nG (z; ) = (a)n 2 F1 ;z .
n− n−+1
a, 1 − b
2 F1 ; z = az−a Bz (a, b),
a+1
we conclude that
nG (z; ) = (a)n z Bz (n − , 1 − a − n),
where Bz (a, b) is the incomplete beta function defined by [49, 8.17.1]
1
Bz (a, b) = z a
t a−1 (1 − zt)b−1 dt. (44)
0
we see that
G
n+1 + (n − ) nG = zn (a)n (1 − z)−a−n = n ,
in agreement with (28).
Modified moments recurrence operator
G
1,0 = 1,0 (S + n − )
= (1 − z)S 2 + [n − + 1 − z (1 + a − + 2n)] S − z (n − )(n + a).
R ECURRENCE RELATIONS FOR THE MOMENTS 107
i) = 0 Linear functional
zx
u (x) (a)x
U(0)
L1,0 [u] = + u (0) .
x=0 x!
First moment differential operator
U(0)
1,0 = ( − 1)1,0 = ( − 1) − z ( + a) ,
which is a special case of (39) with
a1 = a, a2 = 0, b = −1.
Standard moments recurrence operator
U(0)
1,0 = (S − 1)1,0 = S n+1 (S − 1) − z (S + 1)n (S + a)S .
Modified moments recurrence operator
U(0)
1,0 = (S + n)1,0 = (1 − z)S 2 + [n − z (2n + 1 + a)] S − zn (n + a).
For additional references, see [4], [8], [25].
Since the incomplete beta function (44) has the hypergeometric representation [49,
8.17.7]
za 1 − b, a
Bz (a, b) = 2 F1 ;z ,
a a+1
we conclude that
(a)N+1 n
nT (z) = z (1 − z)−a−n B1−z (a + n, N − n + 1),
(N − n)!
and comparing with (33), we see that
(a + n)N−n+1
nT (z) = B1−z (a + n, N − n + 1) n (z) .
(N − n)!
Modified moments recurrence operator
T1,0 = (S + n − N)1,0
= (1 − z)S 2 + [n − N − z (2n + 1 − N + a)] S − z (n − N)(n + a).
R ECURRENCE RELATIONS FOR THE MOMENTS 109
i) = 0 Linear functional
1 zx
u (x)
U(0)
L0,1 [u] = + u (0) .
x=0 (b + 1)x x!
a = 0, b1 = b, b2 = −1.
Standard moments recurrence operator
U(0)
0,1 = (S − 1)0,1 = S n+1 (S + b)(S − 1) − z (S + 1)n S .
Modified moments recurrence operator
U(0)
0,1 = (S + n)0,1 = S (S + n + b)(S + n − 1)− z (S + n).
a = −N, b1 = b, b2 = −N − 1.
zn (a1 )n (a2 )n a1 + n, a2 + n
n (z) = 2 F2 ;z .
(b1 + 1)n (b2 + 1)n b1 + 1 + n, b2 + 1 + n
U0,0 = ( − ) ( − − 1)0,0
= ( − ) ( − − 1) − z ( − + 1)( − ) ,
a1 = − + 1, a2 = − , b1 = − , b2 = − − 1.
Modified moments
2
nC = n+1
C
+ (n − 2) nC
= n+2 + (2n + 1 − 1 − 2 ) n+1 + (n − 1) (n − 2) n
= z2 + (2n + 1 − 1 − 2 ) z + (n − 2) (n − 1) zn ez .
114 D. D OMINICI
2 n − 1 + 1, n − 2 + 1
nC (z) = (n − 1) (n − 2) zn 2 F2 ;z .
n − 1 , n − 2
G,C C
0,0 = 0,0 ( − 2 ) = ( − 1 − 1)( − 2 ) − z ( − 1 + 1)( − 2 ) ,
a1 = −1 + 1, a 2 = − 2 , b1 = −1 − 1, b 2 = − 2 .
G,C C
0,0 = 0,0 (S − 2 )
= S n+1 (S − 1 − 1) − z (S + 1)n (S − 1 + 1) (S − 2) .
n − 1 n n − 1 + 1, n − 2
nG,C (z) = z 2 F2 ;z .
n − 2 n − 1 , n − 2 + 1
G,C C
0,0 = 0,0 (S + n − 2)
= [S (S + n − 1 − 1) − z (S + n − 1 + 1)] (S + n − 2) .
( ) = ( − − 1), ( ) = − + 1,
= 0, + 1, − 1.
R ECURRENCE RELATIONS FOR THE MOMENTS 115
i) = 0
Linear functional
zx
(x − )u (x)
U(0),C
L0,0 [u] = + u (0) .
x=0 x!
a1 = − + 1, a2 = 0, b1 = − − 1, b2 = −1.
ii) = + 1
Linear functional
U( +1),C zx
L0,0 [u] = (x − )u (x) x!
+ u ( + 1)z +1 .
x=0
a1 = − + 1, a2 = − − 1, b1 = − − 1, b2 = − − 2.
iii) = − 1
Linear functional
U( −1),C zx
L0,0 [u] = (x − )u (x) x!
+ u ( − 1)z −1 .
x=0
a 1 = − 1 , a 2 = − 2 , b 1 = − 1 , b 2 = − 2 .
2 zn n − 1 , n − 2
nG (z) = 2 F2 ;z .
(n − 1) (n − 2) n − 1 + 1, n − 2 + 1
( ) = ( − ), ( ) = − , ∈
/ N0 ,
a1 = 0, a2 = − , b1 = −1, b2 = − .
= S n+1 (S − 1) − z (S + 1)n S (S − ).
T,G G
0,0 = ( − N − 1)0,0 = ( − N − 1)0,0 ( − )
= ( − N − 1)( − ) − z ( − N) ( − ),
a1 = −N, a2 = − , b1 = −N − 1, b2 = − .
T,G G
0,0 = (S − N − 1)0,0 = (S − N − 1)0,0 (S − )
= S n+1 (S − N − 1) − z (S + 1)n (S − N) (S − ).
T,G G
0,0 = (S + n − N)0,0 = (S + n − N)0,0 (S + n − )
= [S (S + n − N − 1) − z (S + n − N)] (S + n − ).
Linear functional
U 2 (1,0) zx
L0,0 [u] = u (x) x!
+ 1 u (0) + 2 u (1) z.
x=0
U 2 (1,0) U(0)
0,0 = ( − 2)0,0 = ( − 2)( − 1)0,0
= ( − 1)( − 2) − z ( − 1),
U 2 (1,0) U(0)
0,0 = (S − 2)0,0 = (S − 2)(S − 1)0,0
= S n+1 (S − 2) − z (S + 1)n S (S − 1).
U 2 (1,0) U(0)
0,0 = (S + n − 1)0,0 = (S + n − 1)(S + n)0,0
(49)
= [S (S + n − 2)− z (S + n)] (S + n − 1).
( ) = ( − N − 1), ( ) = − N, N ∈ N0 ,
= 0, N + 1, N.
i) = 0
Linear functional
N
zx
u (x)
U(0),T
L0,0 [u] = + u (0) , N ∈ N0 .
x=0 x!
a1 = −N, a2 = 0, b1 = −N − 1, b2 = −1.
120 D. D OMINICI
ii) = N + 1
Linear functional
N
zx
u (x)
U(N+1),T
L0,0 [u] = + u (N + 1)zN+1 , N ∈ N0 .
x=0 x!
iii) = N
Linear functional
N
zx
u (x)
U(N),T
L0,0 [u] = + u (N) , N ∈ N0 .
x=0 x!
a1 = a, a2 = − + 1, b1 = b, b2 = − − 1.
G
1,1 = 1,1 ( − ) = ( + b)( − ) − z ( + a)( − ) ,
122 D. D OMINICI
a1 = a, a2 = − , b1 = b, b2 = − .
G
1,1 = 1,1 (S − ) = S
n+1
(S + b) − z (S + 1)n (S + a) (S − )
G
1,1 = 1,1 (S + n − ) = [S (S + n + b)− z (S + n + a)] (S + n − ).
( ) = ( + b), ( ) = + a,
= 0, −b, −a.
i) = 0
Linear functional
(a)x zx
u (x)
U(0)
L1,1 [u] = + u (0) .
x=0 (b + 1)x x!
a1 = a, a2 = 0, b1 = b, b2 = −1.
ii) = −a
Linear functional
(a)x zx
u (x)
U(−a)
L1,1 [u] = + u (−a)z−a .
x=0 (b + 1)x x!
U(−a)
1,1 = ( + a)1,1 = ( + b)( + a) − z ( + a)( + a + 1),
a1 = a, a2 = a + 1, b1 = b, b2 = a.
U(−a)
1,1 = (S + a)1,1 = S n+1 (S + b) − z (S + 1)n (S + a + 1) (S + a).
U(−a)
1,1 = (S + n + a + 1)1,1 = [S (S + n + b) − z (S + n + a + 1)] (S + n + a).
iii) = −b
Linear functional
(a)x zx
u (x)
U(−b)
L1,1 [u] = + u (−b)z−b .
x=0 (b + 1)x x!
U(−b)
1,1 = ( + b − 1)1,1 = ( + b)( + b − 1) − z ( + a)( + b),
a1 = a, a2 = b, b1 = b, b2 = b − 1.
U(−b)
1,1 = (S + b − 1)1,1 = S n+1 (S + b − 1)− z (S + 1)n (S + a) (S + b).
U(−b)
1,1 = (S + n + b)1,1 = [S (S + n + b − 1)− z (S + n + a)] (S + n + b).
124 D. D OMINICI
a1 = a, a2 = −N, b1 = b, b2 = −N − 1.
3,0;N = − z ( − N) ( + a1 ) ( + a2) .
Standard moments recurrence operator
n − N, a1 + n, a2 + n
n (z) = (−N)n (a1 )n (a2 )n zn 3 F0 ;z .
−
First moment
−N, a1 , a2
0 (z) = 3 F1 ;z .
b+1
First moment differential operator
G
2,0;N = 2,0;N ( − ) = ( − ) − z ( − N) ( + a)( − )
a1 = −N, a2 = a, a3 = − , b = − .
G
2,0;N = 2,0;N (S − ) = S
n+1
− z (S + 1)n (S − N) (S + a) (S − )
G
2,0;N = 2,0;N (S + n − ) = [S − z (S + n − N)(S + n + a)] (S + n − ).
( ) = , ( ) = ( − N) ( + a), N ∈ N0 ,
= 0, N, −a.
i) = 0
Linear functional
N
zx
u (x) (−N)x (a)x x! + u (0) ,
U(0)
L2,0;N [u] = N ∈ N0 .
x=0
a1 = −N, a2 = a, a3 = 0, b = −1.
R ECURRENCE RELATIONS FOR THE MOMENTS 127
ii) = N
Linear functional
N
zx
u (x) (−N)x (a)x x! + u (N) zN .
U(N)
L2,0;N [u] =
x=0
a1 = −N, a2 = a, a3 = −N + 1, b = −N.
iii) = −a
Linear functional
N
zx
u (x) (−N)x (a)x x! + u (−a)z−a .
U(−a)
L2,0;N [u] =
x=0
a1 = −N, a2 = a, a3 = a + 1, b = a.
128 D. D OMINICI
First moment
a1 , a2 , a3
0 (z) = 3 F2 ;z .
b1 + 1, b2 + 1
First moment differential operator
a1 = a, a2 = − + 1, a3 = − , b1 = − , b2 = − − 1.
G,C C
1,0 = 1,0 ( − 2 ) = ( − 1 − 1)( − 2 ) − z ( + a)( − 1 + 1)( − 2 ) ,
a1 = a, a2 = −1 + 1, a 3 = − 2 , b1 = −1 − 1, b 2 = − 2 .
G,C C
1,0 = 1,0 (S − 2 )
= S n+1 (S − 1 − 1) − z (S + 1)n (S + a)(S − 1 + 1) (S − 2 ) .
G,C C
1,0 = 1,0 (S + n − 2)
= [(S + n − 1 − 1)S − z (S + n + a)(S + n − 1 + 1)] (S + n − 2) .
130 D. D OMINICI
Modified moments
2
nC = n+1
C
+ (n − 2) nC
= n+2 + (2n + 1 − 1 − 2 ) n+1 + (n − 1) (n − 2) n
= 2 z2 + 1 z + 0 zn (a)n (1 − z)−a−n−2 ,
where
2 = (a + 1) (a + 2) ,
1 = a − (a + 2)1 − (a + 1)2 + n(1 + 2a + 1 + 2)
0 = (n − 1) (n − 2) .
2 a + n, n − 1 + 1, n − 2 + 1
nC (z) = (n − 1) (n − 2) (a)n zn 3 F2 ;z .
n − 1 , n − 2
= 0, + 1, −a, − 1.
i) = 0
Linear functional
zx
(x − )u (x) (a)x x! + u (0) .
U(0),C
L1,0 [u] =
x=0
a1 = a, a2 = − + 1, a3 = 0, b1 = − − 1, b2 = −1.
ii) = + 1
Linear functional
U( +1),C zx
L1,0 [u] = (x − )u (x) (a)x x! + u ( + 1)z +1 .
x=0
a1 = a, a2 = − + 1, a3 = − − 1, b1 = − − 1, b2 = − − 2.
132 D. D OMINICI
iii) = −a
Linear functional
zx
(x − )u (x) (a)x x! + u (−a)z−a .
U(−a),C
L0,0 [u] =
x=0
a1 = a, a2 = − + 1, a3 = a + 1, b1 = − − 1, b2 = a.
iv) = − 1
Linear functional
U( −1),C zx
L1,0 [u] = (x − )u (x) (a)x x! + u ( − 1)z −1 .
x=0
a1 = a, a2 = − + 1, a3 = − + 2, b1 = − − 1, b2 = − + 1.
T,C C
1,0 = ( − N − 1)1,0
= ( − − 1)( − N − 1) − z ( + a)( − + 1)( − N),
a1 = a, a2 = − + 1, a3 = −N, b1 = − − 1, b2 = −N − 1.
T,C C
1,0 = (S − N − 1)1,0 = S
n+1
(S − − 1)(S − N − 1)
− z (S + 1)n (S + a)(S − + 1)(S − N) .
T,C C
1,0 = (S + n − N)1,0 = S (S + n − − 1)(S + n − N − 1)
− z (S + n + a)(S + n − + 1)(S + n − N).
a1 = a, a 2 = − 1 , a 3 = − 2 , b 1 = − 1 , b 2 = − 2 .
2 (a)n zn a + n, n − 1, n − 2
nG (z) = 3 F2 ;z .
(n − 1) (n − 2) n − 1 + 1, n − 2 + 1
( ) = ( − ) , ( ) = ( + a)( − ), ∈
/ N0 ,
= 0, −a.
i) = 0
Linear functional
u (x) zx
x−
U(0),G
L1,0 [u] = (a)x + u (0) , ∈
/ N0 .
x=0 x!
a1 = a, a2 = − , a3 = 0, b1 = − , b2 = −1.
R ECURRENCE RELATIONS FOR THE MOMENTS 135
ii) = −a
Linear functional
u (x) zx
x−
U(−a),G
L1,0 [u] = (a)x + u (−a) z−a , ∈
/ N0 .
x=0 x!
a1 = a, a2 = − , a3 = a + 1, b1 = − , b2 = a.
T,G G
1,0 = ( − N − 1)1,0 = ( − N − 1)1,0 ( − )
= ( − )( − N − 1) − z ( + a)( − )( − N) ,
136 D. D OMINICI
a1 = a, a2 = − , a3 = −N b1 = − , b2 = −N − 1.
T,G G
1,0 = (S − N − 1)1,0 = (S − N − 1)1,0 (S − )
= S n+1 (S − N − 1) − z (S + 1)n (S + a)(S − N) (S − ).
T,G G
1,0 = (S + n − N)1,0 = (S + n − N)1,0 (S + n − )
= [S (S + n − N − 1) − z (S + n + a)(S + n − N)] (S + n − ).
= 1, −a, or = 0, −a − 1.
i) = 1
Linear functional
zx
u (x) (a)x
U(1,0)
L1,0 [u] = + 1 u (0) + 2 u (1) z.
x=0 x!
ii) = −a
Linear functional
zx
u (x) (a)x
U(−a,0)
L1,0 [u] = + 1 u (0) + 2 u (−a) z−a .
x=0 x!
a1 = a, a2 = 0, a3 = a + 1, b1 = −1, b2 = a.
iii) = −a − 1
Linear functional
zx
u (x) (a)x
U(−a−1,−a)
L1,0 [u] = + 1 u (−a) z−a + 2 u (−a − 1)z−a−1 .
x=0 x!
a1 = a, a2 = a + 1, a3 = a + 2, b1 = a, b2 = a + 1.
138 D. D OMINICI
( ) = ( − N − 1), ( ) = ( + a)( − N) , N ∈ N0 ,
= 0, N + 1, −a, N.
i) = 0
Linear functional
N
zx
u (x) (a)x x! + u (0) ,
U(0),T
L1,0 [u] = N ∈ N0 .
x=0
a1 = a, a2 = −N, a3 = 0, b1 = −N − 1, b2 = −1.
ii) = N + 1
Linear functional
N
zx
u (x) (a)x x! + u (N + 1)zN+1 ,
U(N+1),T
L1,0 [u] = N ∈ N0 .
x=0
a1 = a, a2 = −N, a3 = −N − 1, b1 = −N − 1, b2 = −N − 2.
iii) = −a
Linear functional
N
zx
u (x) (a)x x! + u (−a)z−a ,
U(−a),T
L1,0 [u] = N ∈ N0 .
x=0
a1 = a, a2 = −N, a3 = a + 1, b1 = −N − 1, b2 = a.
iv) = N
Linear functional
N
zx
u (x) (a)x x! + u (N) zN ,
U(N),T
L1,0 [u] = N ∈ N0 .
x=0
a1 = a, a2 = −N, a3 = −N + 1, b1 = −N − 1, b2 = −N.
a3 = − + 1, b1 = b, b2 = − − 1.
(a1 )n (a2 )n n a + n, a2 + n, n − + 1
nC (z) = (n − ) z 3 F2 1 ;z .
(b + 1)n b + 1 + n, n −
C2,1 = S (S + n + b)(S + n − − 1)
− z (S + n + a1) (S + n + a2) (S + n − + 1).
G
2,1 = 2,1 ( − ) = ( + b)( − ) − z ( + a1 ) ( + a2 ) ( − ) ,
a3 = − , b1 = b, b2 = − .
G
2,1 = 2,1 (S − ) = S
n+1
(S + b) − z (S + 1)n (S + a1 ) (S + a2) (S − )
G
2,1 = 2,1 (S + n − )
= [S (S + n + b) − z (S + n + a1) (S + n + a2)] (S + n − ).
( ) = ( + b), ( ) = ( + a1 ) ( + a2) ,
i) = 0
Linear functional
(a1 )x (a2 )x zx
u (x)
U(0)
L2,1 [u] = + u (0) .
x=0 (b + 1)x x!
a3 = 0, b1 = b, b2 = −1.
a3 = b, b1 = b, b2 = b − 1.
4. Summary
In this section, we list all the functionals and their transformations. Note that we
have
s=0 s=1 s=2
Charlier (0, 0) 4 cases (1, 1) 15 cases (2, 2)
Meixner (1, 0) 5 cases (2, 1) 20 cases (3, 2)
Generalized Charlier (0, 1) 5 cases (1, 2)
Generalized Meixner (1, 1) 6 cases (2, 2)
Generalized Krawtchouk (2, 0; N) 5 cases (3, 1; N)
Generalized Hahn (2, 1) 7 cases (3, 2)
Charlier
Differential operator:
00 = − z.
C0,0 = ( − − 1) − z ( − + 1).
G,C
0,0 = ( − 1 − 1)( − 2 ) − z ( − 1 + 1)( − 2 ) ,
U(0),C
0,0 = ( − − 1)( − 1) − z ( − + 1) ,
U( +1),C
0,0 = ( − − 1)( − − 2) − z ( − + 1)( − − 1),
U( −1),C
0,0 = ( − − 1)( − + 1) − z ( − + 1)( − + 2),
T,C
0,0 = ( − − 1)( − N − 1) − z ( − + 1)( − N) .
G
0,0 = ( − ) − z ( − ).
C,G G,C
0,0 = 0,0 ,
2
G
0,0 = ( − 1 ) ( − 2 ) − z ( − 1 ) ( − 2 ) ,
R ECURRENCE RELATIONS FOR THE MOMENTS 145
U(0),G
0,0 = ( − ) ( − 1) − z ( − ) ,
U( ),G
0,0 = G
0,0 ,
T,G
0,0 = ( − ) ( − N − 1) − z ( − )( − N) .
U 2 (1,0)
0,0 = ( − 1)( − 2) − z ( − 1),
T,U(0)
0,0 = ( − 1)( − N − 1) − z ( − N).
T0,0 = ( − N − 1) − z ( − N).
U(N+1),T
0,0 = ( − N − 1)( − N − 2) − z ( − N) ( − N − 1),
U(N),T
0,0 = ( − N − 1)( − N) − z ( − N) ( − N + 1).
Meixner
Differential operator:
1,0 = − z ( + a).
G,C
1,0 = ( − 1 − 1)( − 2 ) − z ( + a)( − 1 + 1)( − 2 ) ,
U(0),C
1,0 = ( − − 1)( − 1) − z ( + a)( − + 1) ,
U( +1),C
1,0 = ( − − 1)( − − 2) − z ( + a)( − + 1)( − − 1),
U(−a),C
1,0 = ( − − 1)( + a) − z ( + a)( − + 1)( + a + 1),
U( −1),C
1,0 = ( − − 1)( − + 1) − z ( + a)( − + 1)( − + 2),
T,C
1,0 = ( − − 1)( − N − 1) − z ( + a)( − + 1)( − N) .
G
1,0 = ( − ) − z ( + a)( − ) .
C,G G,C
1,0 = 1,0 ,
2
G
1,0 = ( − 1 ) ( − 2 ) − z ( + a)( − 1 ) ( − 2 ) ,
U(0),G
1,0 = ( − )( − 1) − z ( + a)( − ) ,
U(−a),G
1,0 = ( − )( + a) − z ( + a)( − )( + a + 1),
U( ),G
1,0 = G
1,0 ,
G,T
1,0 = ( − ) ( − N − 1) − z ( + a)( − ) ( − N) .
U 2 (1,0)
1,0 = ( − 1)( − 2) − z ( + a) ( − 1),
U 2 (−a,0)
1,0 = ( − 1)( + a) − z ( + a) ( + a + 1),
T,U(0)
1,0 = ( − 1)( − N − 1) − z ( + a) ( − N).
R ECURRENCE RELATIONS FOR THE MOMENTS 147
U(N+1),T
1,0 = ( − N − 1)( − N − 2) − z ( + a)( − N) ( − N − 1),
U(N),T
1,0 = ( − N − 1)( − N) − z ( + a)( − N) ( − N + 1).
Generalized Charlier
Differential operator:
0,1 = ( + b) − z.
Transformations (of type 1, 2):
(i) Christoffel
C0,1 = ( + b)( − − 1) − z ( − + 1),
(ii) Geronimus
G
0,1 = ( + b)( − ) − z ( − ) ,
(iii) Reduced-Uvarov
U(0)
0,1 = ( + b)( − 1) − z ,
U(−b)
0,1 = ( + b)( + b − 1)− z ( + b),
(iv) Truncation
T0,1 = ( + b)( − N − 1) − z ( − N) .
148 D. D OMINICI
Generalized Meixner
Differential operator:
1,1 = ( + b) − z ( + a).
Transformations (of type 2, 2):
(i) Christoffel
(ii) Geronimus
G
1,1 = ( + b)( − ) − z ( + a)( − ),
(iii) Reduced-Uvarov
U(0)
1,1 = ( + b)( − 1) − z ( + a) ,
U(−b)
1,1 = ( + b)( + b − 1) − z ( + a)( + b),
U(−a)
1,1 = ( + b)( + a) − z ( + a)( + a + 1),
(iv) Truncation
Generalized Krawtchouk
Differential operator:
2,0;N = − z ( − N) ( + a)
Transformations (of type 3, 1; N):
(i) Christoffel
(ii) Geronimus
G
2,0;N = ( − ) − z ( − N) ( + a)( − ),
(iii) Reduced-Uvarov
U(0)
2,0;N = ( − 1) − z ( − N) ( + a) ,
U(N)
2,0;N = ( − N) − z ( − N) ( + a)( − N + 1),
U(−a)
2,0;N = ( + a) − z ( − N) ( + a)( + a + 1).
R ECURRENCE RELATIONS FOR THE MOMENTS 149
Differential operator:
2,1 = ( + b) − z ( + a1 ) ( + a2) .
Transformations (of type 3, 2):
(i) Christoffel
(ii) Geronimus
G
2,1 = ( + b)( − ) − z ( + a1 ) ( + a2 ) ( − ),
(iii) Reduced-Uvarov
U(0)
2,1 = ( + b)( − 1) − z ( + a1 ) ( + a2) ,
U(−b)
2,1 = ( + b)( + b − 1)− z ( + a1 ) ( + a2) ( + b),
U(−a )
2,1 1 = ( + b)( + a1) − z ( + a1 ) ( + a2) ( + a1 + 1),
U(−a )
2,1 2 = ( + b)( + a2) − z ( + a1 ) ( + a2) ( + a2 + 1),
(iv) Truncation
5. Conclusion
We obtained recurrence relations for the moments on the monomial and falling factorial
polynomial bases, and gave examples for all functionals of class s 2, where s =
max {p − 1, q}.
We are currently working on further applications of our results to study some prop-
erties of the orthogonal polynomials themselves (representations, recurrence-relation
coefficients, generating functions, etc.).
REFERENCES
[1] M. A LFARO AND R. Á LVAREZ N ODARSE, A characterization of the classical orthogonal discrete
and q -polynomials, J. Comput. Appl. Math. 201 (1), 48–54 (2007).
[2] R. Á LVAREZ N ODARSE , J. A RVES Ú , AND F. M ARCELL ÁN, Modifications of quasi-definite linear
functionals via addition of delta and derivatives of delta Dirac functions, Indag. Math. (N.S.) 15 (1),
1–20 (2004).
[3] R. Á LVAREZ N ODARSE , A. G. G ARC ÍA , AND F. M ARCELL ÁN, On the properties for modifications
of classical orthogonal polynomials of discrete variables, In “Proceedings of the International Confer-
ence on Orthogonality, Moment Problems and Continued Fractions (Delft, 1994)”, vol. 65, pp. 3–18
(1995).
[4] R. Á LVAREZ N ODARSE AND F. M ARCELL ÁN, Difference equation for modifications of Meixner poly-
nomials, J. Math. Anal. Appl. 194 (1), 250–258 (1995).
[5] R. Á LVAREZ N ODARSE AND F. M ARCELL ÁN, The modification of classical Hahn polynomials of a
discrete variable, Integral Transform. Spec. Funct. 3 (4), 243–262 (1995).
[6] R. Á LVAREZ N ODARSE AND J. P ETRONILHO, On the Krall-type discrete polynomials, J. Math. Anal.
Appl. 295 (1), 55–69 (2004).
[7] H. BAVINCK AND R. K OEKOEK, On a difference equation for generalizations of Charlier polynomi-
als, J. Approx. Theory 81 (2), 195–206 (1995).
[8] H. BAVINCK AND H. VAN H AERINGEN, Difference equations for generalized Meixner polynomials,
J. Math. Anal. Appl. 184 (3), 453–463 (1994).
[9] S. B ELMEHDI , On semi-classical linear functionals of class s = 1 . Classification and integral repre-
sentations, Indag. Math. (N.S.) 3 (3), 253–275 (1992).
[10] L. B OELEN , G. F ILIPUK , AND W. VAN A SSCHE, Recurrence coefficients of generalized Meixner
polynomials and Painlevé equations, J. Phys. A 44 (3), 035202, 19 (2011).
[11] M. I. B UENO AND F. M ARCELL ÁN, Darboux transformation and perturbation of linear functionals,
Linear Algebra Appl. 384, 215–242 (2004).
[12] P. L. B UTZER AND T. H. K OORNWINDER, Josef Meixner: his life and his orthogonal polynomials,
Indag. Math. (N.S.) 30 (1), 250–264 (2019).
[13] K. C ASTILLO , F. R. R AFAELI , AND A. S UZUKI , Stieltjes’ theorem for classical discrete orthogonal
polynomials, J. Math. Phys. 61 (10), 103505, 16 (2020).
[14] C. V. L. C HARLIER, Über die Darstellung willkürlicher Funktionen, Ark. Mat., Astr. Fys. 2 (20), 35
(1905–1906).
[15] T. S. C HIHARA, Orthogonal polynomials and measures with end point masses, Rocky Mountain J.
Math. 15 (3), 705–719 (1985).
[16] P. A. C LARKSON, Recurrence coefficients for discrete orthonormal polynomials and the Painlevé
equations, J. Phys. A 46 (18), 185205, 18 (2013).
[17] D. D ICKINSON, On Lommel and Bessel polynomials, Proc. Amer. Math. Soc. 5, 946–956 (1954).
[18] D. D OMINICI , Asymptotic analysis of the Askey-scheme. I. From Krawtchouk to Charlier, Cent. Eur.
J. Math. 5 (2), 280–304 (2007).
[19] D. D OMINICI , Asymptotic analysis of the Krawtchouk polynomials by the WKB method, Ramanujan
J. 15 (3), 303–338 (2008).
[20] D. D OMINICI , Laguerre-Freud equations for generalized Hahn polynomials of type I, J. Difference
Equ. Appl. 24 (6), 916–940 (2018).
[21] D. D OMINICI , Recurrence relations for the moments of discrete semiclassical orthogonal polynomials,
J. Class. Anal. 20 (2), 143–180 (2022).
[22] D. D OMINICI AND F. M ARCELL ÁN, Discrete semiclassical orthogonal polynomials of class one,
Pacific J. Math. 268 (2), 389–411 (2014).
[23] D. D OMINICI AND F. M ARCELL ÁN, Discrete semiclassical orthogonal polynomials of class 2, In
“Orthogonal polynomials: current trends and applications”, vol. 22 of “SEMA SIMAI Springer Ser.”,
pp. 103–169, Springer, Cham ([2021] 2021).
c
[24] A. J. D UR ÁN, Orthogonal polynomials satisfying higher-order difference equations, Constr. Approx.
36 (3), 459–486 (2012).
[25] A. J. D UR ÁN, From Krall discrete orthogonal polynomials to Krall polynomials, J. Math. Anal. Appl.
450 (2), 888–900 (2017).
R ECURRENCE RELATIONS FOR THE MOMENTS 151
[26] A. J. D UR ÁN, The algebras of difference operators associated to Krall-Charlier orthogonal polyno-
mials, J. Approx. Theory 234, 64–81 (2018).
[27] A. J. D UR ÁN, Christoffel transform of classical discrete measures and invariance of determinants
of classical and classical discrete polynomials, J. Math. Anal. Appl. 503 (2), Paper No. 125306, 29
(2021).
[28] P. F EINSILVER , J. M C S ORLEY, AND R. S CHOTT, Combinatorial interpretation and operator calcu-
lus of Lommel polynomials, J. Combin. Theory Ser. A 75 (1), 163–171 (1996).
[29] G. F ILIPUK AND W. VAN A SSCHE, Recurrence coefficients of a new generalization of the Meixner
polynomials, SIGMA Symmetry Integrability Geom. Methods Appl. 7, Paper 068, 11 (2011).
[30] G. F ILIPUK AND W. VAN A SSCHE, Discrete orthogonal polynomials with hypergeometric weights
and Painlevé VI, SIGMA Symmetry Integrability Geom. Methods Appl. 14, Paper No. 088, 19 (2018).
[31] M. F OUPOUAGNIGNI , M. N. H OUNKONNOU , AND A. RONVEAUX, Laguerre-Freud equations for
the recurrence coefficients of D semi-classical orthogonal polynomials of class one, In “Proceedings
of the VIIIth Symposium on Orthogonal Polynomials and Their Applications (Seville, 1997)”, vol. 99,
pp. 143–154 (1998).
[32] A. G. G ARC ÍA , F. M ARCELL ÁN , AND L. S ALTO, A distributional study of discrete classical orthog-
onal polynomials, In “Proceedings of the Fourth International Symposium on Orthogonal Polynomials
and their Applications (Evian-Les-Bains, 1992)”, vol. 57, pp. 147–162 (1995).
[33] W. G AUTSCHI , On generating orthogonal polynomials, SIAM J. Sci. Statist. Comput. 3 (3), 289–317
(1982).
[34] E. G ODOY, F. M ARCELL ÁN , L. S ALTO , AND A. Z ARZO, Perturbations of discrete semiclassical
functionals by Dirac masses, Integral Transform. Spec. Funct. 5 (1–2), 19–46 (1997).
[35] M. N. H OUNKONNOU , C. H OUNGA , AND A. RONVEAUX, Discrete semi-classical orthogonal poly-
nomials: generalized Charlier, J. Comput. Appl. Math. 114 (2), 361–366 (2000).
[36] S. K ARLIN AND J. L. M C G REGOR, The Hahn polynomials, formulas and an application, Scripta
Math. 26, 33–46 (1961).
[37] M. P. K RAVCHUK, Sur une généralisation des polynômes d’Hermite, C. R. Acad. Sci. Paris 189,
620–622 (1929).
[38] K. F. L EE AND R. W ONG, Asymptotic expansion of the modified Lommel polynomials hn, (x) and
their zeros, Proc. Amer. Math. Soc. 142 (11), 3953–3964 (2014).
[39] M. M A ÑAS , Pearson equations for discrete orthogonal polynomials: III – Christoffel and Geronimus
transformations, Rev. R. Acad. Cienc. Exactas Fı́s. Nat. Ser. A Mat. RACSAM 116 (4), Paper No.
168, 23 (2022).
[40] M. M A ÑAS , I. F ERN ÁNDEZ -I RISARRI , AND O. F. G ONZ ÁLEZ -H ERN ÁNDEZ, Pearson equations for
discrete orthogonal polynomials: I. Generalized hypergeometric functions and Toda equations, Stud.
Appl. Math. 148 (3), 1141–1179 (2022).
[41] I. G. M ACDONALD, Symmetric functions and Hall polynomials, Oxford Classic Texts in the Physical
Sciences. The Clarendon Press, Oxford University Press, New York, second ed. (2015).
[42] D. M AKI , On constructing distribution functions: With applications to Lommel polynomials and
Bessel functions, Trans. Amer. Math. Soc. 130, 281–297 (1968).
[43] F. M ARCELL ÁN AND L. S ALTO, Discrete semi-classical orthogonal polynomials, J. Differ. Equations
Appl. 4 (5), 463–496 (1998).
[44] P. M ARONI , Sur la suite de polynômes orthogonaux associée à la forme u = c d + (x − c)−1 L ,
Period. Math. Hungar. 21 (3), 223–248 (1990).
[45] P. M ARONI AND M. M EJRI , The symmetric D -semi-classical orthogonal polynomials of class one,
Numer. Algorithms 49 (1–4), 251–282 (2008).
[46] J. M EIXNER, Orthogonale Polynomsysteme Mit Einer Besonderen Gestalt Der Erzeugenden Funktion,
J. London Math. Soc. 9 (1), 6–13 (1934).
[47] A. F. N IKIFOROV, S. K. S USLOV, AND V. B. U VAROV, Classical orthogonal polynomials of a dis-
crete variable, Springer Series in Computational Physics. Springer-Verlag, Berlin (1991).
[48] P. N JIONOU S ADJANG , W. K OEPF, AND M. F OUPOUAGNIGNI , On moments of classical orthogonal
polynomials, J. Math. Anal. Appl. 424 (1), 122–151 (2015).
[49] F. W. J. O LVER , D. W. L OZIER , R. F. B OISVERT, AND C. W. C LARK (eds.), NIST handbook of
mathematical functions, U.S. Department of Commerce, National Institute of Standards and Technol-
ogy, Washington, DC; Cambridge University Press, Cambridge (2010).
152 D. D OMINICI
[50] A. RONVEAUX , S. B ELMEHDI , E. G ODOY, AND A. Z ARZO, Recurrence relation approach for con-
nection coefficients. Applications to classical discrete orthogonal polynomials, In “Symmetries and
integrability of difference equations (Estérel, PQ, 1994)”, vol. 9 of “CRM Proc. Lecture Notes”, pp.
319–335. Amer. Math. Soc., Providence, RI (1996).
[51] A. RONVEAUX AND L. S ALTO, Discrete orthogonal polynomials – polynomial modification of a
classical functional, J. Differ. Equations Appl. 7 (3), 323–344 (2001).
[52] M. S GHAIER AND M. Z AATRA, A large family of Dw -semiclassical polynomials of class one, Integral
Transforms Spec. Funct. 28 (5), 386–402 (2017).
[53] M. S GHAIER AND M. Z AATRA, A family of discrete semi-classical orthogonal polynomials of class
one, Period. Math. Hungar. 76 (1), 68–87 (2018).
[54] C. S MET AND W. VAN A SSCHE, Orthogonal polynomials on a bi-lattice, Constr. Approx. 36 (2),
215–242 (2012).
[55] L. T RUKSA, Hypergeometric orthogonal systems of polynomials. i, Aktuárské Vědy 2 (2), 65–84
(1931).
[56] L. T RUKSA, Hypergeometric orthogonal systems of polynomials. ii, Aktuárské Vědy 2 (3), 113–144
(1931).
[57] L. T RUKSA, Hypergeometric orthogonal systems of polynomials. iii, Aktuárské Vědy 2 (4), 177–203
(1931).
[58] W. VAN A SSCHE AND M. F OUPOUAGNIGNI, Analysis of non-linear recurrence relations for the
recurrence coefficients of generalized Charlier polynomials, J. Nonlinear Math. Phys. 10 (suppl. 2),
231–237 (2003).
[59] L. V ERDE -S TAR, Recurrence coefficients and difference equations of classical discrete orthogonal
and q -orthogonal polynomial sequences, Linear Algebra Appl. 440, 293–306 (2014).
[60] L. V ERDE -S TAR, A unified construction of all the hypergeometric and basic hypergeometric families
of orthogonal polynomial sequences, Linear Algebra Appl. 627, 242–274 (2021).
[61] L. V INET AND A. Z HEDANOV, Hypergeometric orthogonal polynomials with respect to Newtonian
bases, SIGMA Symmetry Integrability Geom. Methods Appl. 12, Paper No. 048, 14 (2016).