On the Polynomial Solution of Divided-Difference Equations of the Hypergeometric Type on Nonuniform Lattices

: In this paper, we provide a formal proof of the existence of a polynomial solution of ﬁxed degree for a second-order divided-difference equation of the hypergeometric type on non-uniform lattices, generalizing therefore previous work proving existence of the polynomial solution for second-order differential, difference or q -difference equation of hypergeometric type. This is achieved by studying the properties of the mean operator and the divided-difference operator as well as by deﬁning explicitly, the right and the “left” inverse for the second operator. The method constructed to provide this formal proof is likely to play an important role in the characterization of orthogonal polynomials on non-uniform lattices and might also be used to provide hypergeometric representation (when it does exist) of the second solution—non polynomial solution—of a second-order divided-difference equation of hypergeometric type.


Introduction
Classical orthogonal polynomials of a continuous variable (P n ) 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 λ = λ n . This is achieved mainly by showing that : the nth derivative y (n) of any solution y of (1) satisfies an equation of the same type (hypergeometric aspect), that is, an equation of the form σ(x) y (x) + τ n (x) y (x) + λ n y(x) = 0, where τ n is a first degree polynomial and λ n is a constant given by τ n (x) = τ + n σ (x), λ n = λ + n τ + n(n − 1) 2 σ .
The fact that the constant solution of (2) when λ n = 0 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 B n 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 D q is the Hahn operator [8] defined by provided that f (0) 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]) ∇y(x(s)) ∇x(s) ∆y(x(s)) ∆x(s) + ∇y(x(s)) ∇x(s) + λ n y(x(s)) = 0, where ψ and φ are polynomials of degree 1 and at most 2, respectively; λ n is a constant depending on n and on the leading coefficients of φ and ψ. The lattice x(s) is defined by [9,10] x(s) = is known as non-uniform lattice and fulfills various important properties. Equation (6) can be transformed into equation [12] φ(x(s))D 2 x y(x(s)) + ψ(x(s))S x D x y(x(s)) + λ n y(x(s)) = 0, called divided-difference equation of the hypergeometric type by means of the two companion operators D x (called divided-difference operator) and S x (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 D x and S x 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 D x y(x(s)) = 2wq 1 4 1 − q y(x(s)), 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: 1. redto define the right and the "left" inverses of the operator D x ; 2. to provide a formal proof of the existence of a polynomial solution of a preassigned degree of the divided-difference equation of hypergeometric type (8), extending and generalising therefore-by means of specialisation and limiting situations on the lattice x(s)-similar results obtained for second-order differential, difference or q-difference equation of hypergeometric type.

Preliminary Results: Known and New Properties
Since the main result of this paper is based on the operators D x and S x which are defined by using the lattice x(s), we will provide in this section some known and basic properties of x(s), D x and S x . We will also derive new properties such as the right and the "left" inverses of the operator D x , required in the next section.

Known Properties of the Operators D x and S x
The operators S x and D x fulfil the so-called Product rules I [13,14]: where U 2 is a polynomial of degree 2 and η x is a constant given by [14] The operators D x and S x also satisfy the so-called Product Rules II [13,14]: where I is the identity operator I f (x) = f (x), and

An Appropriate Basis for the Operators D x and S x
Searching for a polynomial basis on which the action of the companion operators will give a linear combination of at most two elements of the basis, Foupouagnigni et al. proved in [14] that the polynomial F n defined by where ε is the unique solution (provided that the lattice x(s) is quadratic or q-quadratic: i.e., the constants c j in (7) satisfy is the right basis for the operators S x and D x because it satisfies the following properties where the constants α n and γ n are given in (14) and (15). After reviewing some properties of the operators D x and S x , we now state and prove the following proposition providing the left and right inverse of the operator D x , to be used in the next section to complete the proof of the main theorem of this paper.

Proposition 1.
Let F x be a linear operator defined on the basis (F n ) n by Then, F x satisfies the following relations: where I is the identity operator and δ x(ε) is the Dirac delta distribution defined by δ x(ε) , P = P(x(ε)), ∀P, with ε is defined in (21).
Proof. For all positive integer n, F n defined by (21) is a polynomial of degree exactly n.
We have Hence, the first part of Relation (24) holds. Using (26), we have The second part of (25) is therefore satisfied.

Existence of the Polynomial Solution of the Divided-Difference Equation of the Hypergeometric Type
Having stated and proved required properties of the operators D x and S x , 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 that Then, the divided-difference equation φ(x(s)) D 2 x y(x(s)) + ψ(x(s)) S x D x y(x(s)) + λ n,0 y(x(s)) = 0, with where ψ 1 and φ 2 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 y 0 is a solution of (28), then the function y 1 = D x y 0 satisfies φ [1] (x(s)) D 2 x y(x(s)) + ψ [1] (x(s)) S x D x y(x(s)) + λ n,1 y(x(s)) = 0, where Proof. Assume that y 0 satisfies (28). Applying the operator D x to (28) in which y is replaced by y 0 and using the product rule I in (16) and (17), we obtain Using the product rules II in (20) to replace S 2 x and D x S x in the previous equation, we have where φ [1] (x(s)), ψ [1] (x(s)) and λ n,1 are defined by (31). Therefore, y 1 = D x y 0 is a solution of Equation (30).

Lemma 2.
If the function y 0 is a solution of (28), then the function y k = D k x y 0 is a solution of the equation where the polynomials φ [k] , ψ [k] and the constant λ n,k satisfy with the following initial values: Proof. Lemma 1 assures the validity of the result for k = 1. Let k be a positive integer. Assume that y k is solution of Equation (32). Applying the operator D x to (32) in which y is replaced by y k and using the Product Rules I, we obtain Using the products rule II to replace S 2 x and D x S x in the previous equation, we have where φ [k+1] (x(s)), ψ [k+1] (x(s)) and λ n,k+1 are defined by (33). Thus, y k+1 = D x y k satisfies φ [k+1] (x(s)) D 2 x y k+1 (x(s)) + ψ [k+1] (x(s)) S x D x y k+1 (x(s)) + λ n,k+1 y k+1 (x(s)) = 0.
We prove that v 0 is solution of (28) by replacing y 1 in the Equation (35) by D x v 0 .
0 , then from (33), we have the following system of recurrence equation 1 . (12) and (13), the previous system of equations becomes
Using the definition of λ n,0 which of course coincides with the one of λ n,k for k = 0, we derive (37) from (39).
Solving the following the equation in terms of the unknown k keeping in mind (27), gives a unique solution k = n. 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 0 ≤ k ≤ n. Then, if y k is a solution of Equation (32), then there exists y 0 solution of Equation (28) such that Proof. Let k be a nonnegative integer with k ≤ n. Assume that y k satisfies (32). Then, we obtain that there exists a function y k−1 solution of the equation obtained by replacing k in (32) by k − 1, namely, such that y k = D x y k−1 .
This is achieved using the fact that λ n,k−1 = 0 thanks to Lemma 4, and also using Lemma 3 but with the functions y 0 and y 1 replaced, respectively, by the functions y k−1 and y k 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 y k−1 , y k−2 , . . . , y 1 and using Lemmas 3 and 4.
Next, we apply the operator F x on both members of the previous equation and deduce by applying the second relation of Equation (25) of Proposition 1 that Hence, where C n−1 = D n−1 x v 0 (x(s)) | s=ε . By applying again the operator F x on both members of the previous equation and using the second relation of Equation (25), we get where C n−2 = D n−2 x v 0 (x(s)) | s=ε . Repeating the same process, we express v 0 as v 0 (x(s)) = F n x F 0 (x(s)) + where C j = D j x v 0 (x(s)) | s=ε . Therefore, v 0 (x(s)) is a polynomial of degree exactly n in x(s).

Conclusions and Perspectives
In this work, we have derived the right and the "left" inverse of the operator D x and used the properties of the inverse operators, as well as those of the operators D x and S x , 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 b = a q  (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 D x 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).
Author Contributions: These authors contributed equally to this work, in terms of investigation, conceptualization and methodology.
Funding: The Research of the first author was partially supported by the AIMS-Cameroon Research Allowance 2018-2019.