1. Introduction
Classical orthogonal polynomials of a continuous variable
are known to satisfy a second-order differential equation of hypergeometric type
where
is a polynomial of degree at most 2,
is a first degree polynomial and
is a constant with respect to
x.
In [
1,
2], it is shown that Equation (
1) has a polynomial solution of exactly
n degree for a specific given constant
. This is achieved mainly by showing that:
- -
the
nth derivative
of any solution
y of (
1) satisfies an equation of the same type (hypergeometric aspect), that is, an equation of the form
where
is a first degree polynomial and
is a constant given by
- -
Any solution of (
2) can be written as the
nth derivative of a solution of (
1), provided that
.
The fact that the constant solution of (
2) when
is the
nth derivative of a solution of (
1) leads to the existence of a polynomial solution of (
1), of exactly
n degree, when
This result proves not only the existence of a polynomial solution for Equation (
1) but also allows for establishing the Rodrigues formula expressing the polynomial solution in the term of the
nth derivative:
where
is a constant and
is the weight function satisfying the Pearson equation
It is worth mentioning that Hermite, Laguerre, Jacobi and Bessel polynomials are the polynomial eigenfunctions of the second-order linear differential operation given in (
1).
Using the same approach, similar results have been established in [
3] (See also [
4]) for the classical orthogonal polynomials of a discrete variable satisfying instead a second-order difference equation of hypergeometric type
where
and ∇ are the forward and the backward operators defined by
Furthermore, it should be noticed that Charlier, Krawtchuk, Meixner and Hahn polynomials are the polynomial eigenfunctions of the second-order linear difference operation given in (
4).
The same result can be established in the same way to for the classical orthogonal polynomials of a
q-discrete variable satisfying a second-order
q-difference equation of hypergeometric type [
5] (See also [
6,
7])
where
is the Hahn operator [
8] defined by
provided that
exists. Orthogonal polynomials which are eigenfunctions of the second-order
q-difference operator given defined by (
5) are [
5]: Big q-Jacobi, Big q-Laguerre, Little q-Jacobi, Little q-Laguerre (Wall), q-Laguerre, Alternative q-Charlier, Al-Salam–Carlitz I, Al-Salam–Carlitz II, Stieltjes–Wigert, Discrete q–Hermite I, Discrete q–Hermite II, q-Hahn, q-Meixner, Quantum q-Krawtchouk, q-Krawtchouk, Affine q-Krawtchouk, the q-Charlier and the q-Charlier II polynomials.
Classical orthogonal polynomials on non-uniform lattices (including but not limited to Askey–Wilson polynomials, Racah and
q-Racah polynomials), are known to satisfy a second-order divided-difference equation of the form [
9,
10] (see also [
11])
where
and
are polynomials of degree 1 and at most 2, respectively;
is a constant depending on
n and on the leading coefficients of
and
. The lattice
is defined by [
9,
10]
is known as non-uniform lattice and fulfills various important properties.
Equation (
6) can be transformed into equation [
12]
called divided-difference equation of the hypergeometric type by means of the two companion operators
(called divided-difference operator) and
(called mean operator) defined as [
9,
10,
12,
13]
Using appropriate bases, computer algebra software has been used to solve divided-difference Equation (
8) for specific families of classical orthogonal polynomials on a non-uniform lattice. For some special values of the parameter for the specific case of Askey–Wilson polynomials, non-polynomial solution has been recovered together with the polynomial one [
14] (see page 15, Equations (62) and (63)). In addition, the operators
and
have played a decisive role not only for establishing the functional approach of the characterization theorem of classical orthogonal polynomials on non-uniform lattices, but also for providing algorithmic solution to linear homogeneous divided-difference equations with polynomial coefficients, allowing to solve explicitly [
13] the first-order divided-difference equations satisfied by the basic exponential function
and the second-order divided-difference equation satisfied by the basic trigonometric functions
where
w is a given constant.
The aim of this work is:
3. Existence of the Polynomial Solution of the Divided-Difference Equation of the Hypergeometric Type
Having stated and proved required properties of the operators and , we will now state and prove the main theorem of this paper.
Theorem 1. Let n be a nonnegative integer, ψ and ϕ be two polynomials of degree 1 and at most 2, respectively, such thatThen, the divided-difference equationwithwhere and are leading coefficients of polynomials ψ and ϕ respectively, has a polynomial solution of exactly n degree. The proof of Theorem 1 will be organized as follows: we split the proof in five lemmas which we first state, prove, and then put these lemmas together in combination with Proposition 1 to deduce the proof of this theorem.
Lemma 1.
If the function is a solution of (28), then the function satisfieswhere Proof. Assume that
satisfies (
28). Applying the operator
to (
28) in which
y is replaced by
and using the product rule I in (
16) and (17), we obtain
Using the product rules II in (
20) to replace
and
in the previous equation, we have
where
,
and
are defined by (
31). Therefore,
is a solution of Equation (
30). □
Lemma 2.
If the function is a solution of (28), then the function is a solution of the equationwhere the polynomials and the constant satisfywith the following initial values: Proof. Lemma 1 assures the validity of the result for .
Let
k be a positive integer. Assume that
is solution of Equation (
32). Applying the operator
to (
32) in which
y is replaced by
and using the Product Rules I, we obtain
Using the products rule II to replace
and
in the previous equation, we have
where
,
and
are defined by (
33). Thus,
satisfies
□
Lemma 3.
If a given function satisfies (30) with , then there exists a function satisfying (28) such that Proof. Let
be a solution of (
30) with
. If there would exist a solution
of (
28) such that
, then from (
28) we can express
as:
Now, it remains to verify that the function
defined in terms of
by (
35) satisfies Equation (
28) with
By applying
to (
35) and using product rules I, II and the fact that
is solution of (
30), we get
Therefore,
since
.
We prove that
is solution of (
28) by replacing
in the Equation (
35) by
. □
Lemma 4.
For any positive integer n, the coefficients defined by relation (33) satisfywhere Proof. If we denote by
and
, then from (
33), we have the following system of recurrence equation
Using relations
derived from Equations (
12) and (
13), the previous system of equations becomes
Solving this system of recurrence equations with the initial values
,
, we obtain for the
q-quadratic lattice
Using the definition of
which of course coincides with the one of
for
, we derive (
37) from (
39).
Solving the following the equation
in terms of the unknown
k keeping in mind (
27), gives a unique solution
. Thus, relation (38) is satisfied. It can easily be proved in the same way that relation (38) is satisfied for the quadratic lattice. □
Lemma 5.
Let n be a fixed positive integer and let k be an integer such that . Then, if is a solution of Equation (32), then there exists solution of Equation (28) such that Proof. Let
k be a nonnegative integer with
. Assume that
satisfies (
32). Then, we obtain that there exists a function
solution of the equation obtained by replacing
k in (
32) by
, namely,
such that
This is achieved using the fact that
thanks to Lemma 4, and also using Lemma 3 but with the functions
and
replaced, respectively, by the functions
and
while Equation (
28) is replaced by Equation (
40). In addition, Equation (
30) is replaced by the Equation (
32). □
The proof is completed by repeating the same process for and using Lemmas 3 and 4.
Proof of Theorem 1. Since, for
thanks to (
37), Equation (
32) admits a constant solution, namely
. We therefore deduce from Lemma 5 that there exists a function
solution of (
28) such that
Next, we apply the operator
on both members of the previous equation and deduce by applying the second relation of Equation (
25) of Proposition 1 that
Hence,
where
.
By applying again the operator
on both members of the previous equation and using the second relation of Equation (
25), we get
where
. Repeating the same process, we express
as
where
. Therefore,
is a polynomial of degree exactly
n in
. □
4. Conclusions and Perspectives
In this work, we have derived the right and the “left” inverse of the operator
and used the properties of the inverse operators, as well as those of the operators
and
, to provide a formal proof that the divided-difference equation of hypergeometric type (
28) has a polynomial solution of degree exactly
n.
The novelty of our work is the formal proof of the existence of this polynomial solution, confirming therefore the fact that, in [
14], by solving divided-difference (
8) on a case by case basis and using most appropriate polynomial basis for each case, we have obtained for each family of classical orthogonal polynomials on non-uniform lattice, a hypergeometric or
q-hypergeometric solution which happens to be a polynomial because of the form of one of the upper parameters obtained in the hypergeometric (or
q-hypergeometric) representation of the obtained solution.
Finding hypergeometric representation of the non polynomial solution of (
8) is not obvious and this was obtained unexpectedly for the Askey–Wilson polynomials when the parameters fulfill
[
14] (see page 15, Equations (62) and (63)). The method developed here might help to understand when and why such a hypergeometric representation exists for non-polynomial solutions.
As an additional potential application of our paper, the right and the “left” inverse of the operator
are likely to play important role in the study of the properties of orthogonal polynomials on the non-uniform latices, and on the search of the solutions of divided-difference equations on non-uniform lattices, as well as on the hypergeometric representation (when they exist) of the second-solution—non polynomial solution—of Equation (
28).