# On classical orthogonal polynomials and the Cholesky factorization of a
class of Hankel matrices

Misael E. Marriaga*, Guillermo Vera de Salas, Marta Latorre, Rubén Muñoz
Alcázar  Departamento de Matemática Aplicada, Ciencia e Ingeniería de
Materiales y Tecnología Electrónica, Universidad Rey Juan Carlos (Spain)
<EMAIL_ADDRESS>Departamento de Matemática Aplicada, Ciencia e
Ingeniería de Materiales y Tecnología Electrónica, Universidad Rey Juan Carlos
(Spain<EMAIL_ADDRESS>Departamento de Matemática Aplicada, Ciencia e
Ingeniería de Materiales y Tecnología Electrónica, Universidad Rey Juan Carlos
(Spain<EMAIL_ADDRESS>Departamento de Álgebra, Geometría y Topología,
Universidad de Málaga (Spain<EMAIL_ADDRESS>
###### Abstract.

Classical moment functionals (Hermite, Laguerre, Jacobi, Bessel) can be
characterized as those linear functionals whose moments satisfy a second order
linear recurrence relation. In this work, we use this characterization to link
the theory of classical orthogonal polynomials and the study of Hankel
matrices whose entries satisfy a second order linear recurrence relation.
Using the recurrent character of the entries of such Hankel matrices, we give
several characterizations of the triangular and diagonal matrices involved in
their Cholesky factorization and connect them with a corresponding
characterization of classical orthogonal polynomials.

###### Key words and phrases:

Orthogonal polynomials, Hankel matrices, Cholesky factorization.

###### 2010 Mathematics Subject Classification:

15A23, 15A23, 15A63, 33C45, 33D45

*Corresponding Author: Misael E. Marriaga

##  1\. Introduction

Classical orthogonal polynomials (Hermite, Laguerre, Jacobi, and Bessel) have
been characterized using different approaches. For instance, they can be
characterized in terms of differential equations [3], their derivatives ([10,
14]), structure relations ([2, 6, 16]), and a Rodrigues formula ([21]), among
others (see [7, 17] and the references therein). More recently, an approach
that uses linear functionals and duality was introduced by P. Maroni ([18]).
In all of these approaches, the starting point is to use the basis of
monomials to represent polynomials and state the results.

However, an interesting (and most recent) approach is to start from the theory
of semi-infinite matrices. The bibliography on this subject has grown greatly
in the last years, and it has become increasingly difficult to do a
comprehensive review of all the references. Hence, we refer the reader to [22,
23] (and the refereces therein) where the algebra of infinite triangular
matrices and the algebra of infinite Hessenberg matrices are used to study
some aspects of orthogonal polynomials, and to [5, 15] (and the references
therein) where the main tool is the Cholesky factorization of Gram matrices of
bilinear forms. We remark that the Cholesky factorization proves to be quite
fruitful in the study of non standard orthogonality such as multiple, matrix,
Sobolev, and multivariate orthogonality as well as orthogonality on the unit
circle of the complex plane, and have successfully found its way into
applications in random matrices, Toda lattices, integrable systems, Riemann-
Hilbert problems, Painlevé equations, and Darboux transformations, among
others topics.

Our goal is to contribute to the link between matrix factorization and
orthogonal polynomials. In particular, we deal with several characterizations
of classical orthogonal polynomials. However, we shift our paradigm from
infinite matrices to the finite Gram matrix $G_{n}$ associated with a bilinear
form defined on the linear space of polynomials of degree at most $n\geqslant
0$. For standard orthogonality, $G_{n}$ is a Hankel matrix (all of its
antidiagonals are constant). Taking into account that the moments of a linear
functional associated with a family of classical orthogonal polynomials
satisfy a second order linear recurrence relation ([16]), we can say that this
paper deals with Hankel matrices with an additional structure: the entries of
$G_{n}$ satisfy such recurrence relation. In this way, we can extend the
bilinear form to the linear space of polynomials of degree at most $n+1$ by
constructing a new Gram matrix $G_{n+1}$ by means of bordering $G_{n}$ with a
new row and column whose entries are obtained using the recurrence relation
and the entries of $G_{n}$. The resulting matrix $G_{n+1}$ will also be a
Hankel matrix with the additional structure mentioned above. Consequently, it
will be possible to prove by induction that the properties satisfied by
$G_{n}$ are also satisfied by $G_{n+1}$.

The change from infinite matrices to subsequently bordering finite matrices is
motivated by an alternative proof of a classical result about the interlacing
of zeros of orthogonal polynomials of consecutive degrees found in [7]. This
classical result states that if the zeros of any two polynomials of
consecutive degrees interlace, then these polynomials are elements of a
sequence of orthogonal polynomials associated with a positive definite moment
functional. This result can be proved using the Euclidean division algorithm
for polynomials. However, this result can also be deduced using the following
theorem about interlacing eigenvalues of Hermitian matrices found in [12, p.
185]:

######  Theorem 1.1.

Let $n$ be a given positive integer, and let $\\{x_{n,k}\\}_{k=1}^{n}$ and
$\\{x_{n-1,k}\\}_{k=1}^{n-1},$ be two given sequences of numbers such that

$x_{n,1}<x_{n-1,1}<\cdots<x_{n,k}<x_{n-1,k}<x_{n,k+1}<\cdots<x_{n-1,n-1}<x_{n,n}.$

Let $\Lambda=\operatorname{diag}(x_{n-1,1},x_{n-1,2},\cdots,x_{n-1,n-1})$.
Then there exists real number $b$ and a real vector
$\bar{y}=(y_{1},\ldots,y_{n-1})^{\top}\in\mathbb{R}^{n-1}$ such that
$(x_{n,k})_{k=1}^{n}$ is the set of eigenvalues of the real symmetric matrix

$B=\left(\begin{array}[]{c|c}\Lambda&\bar{y}\\\ \hline\cr\\\\[-6.0pt]
\bar{y}^{\top}&b\end{array}\right).$

From this, it seems reasonable to think that the procedure of bordering
matrices encoding information about polynomial sequences is well suited for
presenting and deducing results about orthogonal polynomials most likely due
to the fact that for all $n\geqslant 0$, the linear space of polynomials of
degree at most $n$ is a subspace of the space of polynomials of degree at most
$n+1$. In this way, the Gram matrices of bilinear forms associated with
classical orthogonal polynomials posses the adequate structure to start
exploring our proposed paradigm.

The paper is organized as follows. Section 2 presents basic background on
classical orthogonal polynomials and their associated linear functionals, and
we introduce classical sequences of real numbers in Section 3. In Section 4,
we discuss the Cholesky factorization of Hankel matrices obtained from given
sequences of real number and its relation to orthogonal polynomials. We
present several characterizations of classical sequences of real numbers in
Section 5.

##  2\. Orthogonal polynomials and linear functionals

For $n\geqslant 0$, let $\Pi_{n}$ be the linear space of polynomials of degree
at most $n$ of a real variable and real coefficients, and let
$\Pi=\bigcup_{n\geqslant 0}\Pi_{n}$.

Let $\Pi^{*}$ denote the algebraic dual space of $\Pi$. That is, $\Pi^{*}$ is
the linear space of linear functionals defined on $\Pi$,

$\Pi^{*}=\left\\{\mathbf{u}:\Pi\rightarrow\mathbb{R}\ \ :\ \ \mathbf{u}\text{
is linear}\right\\}.$

We denote by $\langle\mathbf{u},p\rangle$ the image of the polynomials $p$
under the linear functional $\mathbf{u}$.

Any linear functional $\mathbf{u}$ is completely defined by the values

$\mu_{n}\,:=\,\langle\mathbf{u},x^{n}\rangle,\quad n\geqslant 0,$

and extended by linearity to all polynomials, where $\mu_{n}$ is called the
$n$-th moment of $\mathbf{u}$. Therefore, we refer to $\mathbf{u}$ as a moment
functional.

A moment functional $\mathbf{u}$ is called positive definite if
$\langle\mathbf{u},p^{2}\rangle>0$ for every non zero polynomial $p\in\Pi$.

Let $\mathbf{u}$ be a moment functional. A sequence of polynomials
$\\{P_{n}(x)\\}_{n\geqslant 0}$ is called an orthogonal polynomial sequence
(OPS) with respect to $\mathbf{u}$ if

  1. (1)

$\deg\,P_{n}=n$,

  2. (2)

$\langle\mathbf{u},P_{n}\,P_{m}\rangle=h_{n}\,\delta_{n,m}$, with $h_{n}\neq
0$.

Here $\delta_{n,m}$ denotes the Kronecker delta defined as

$\delta_{n,m}=\left\\{\begin{array}[]{ll}1,&n=m,\\\ 0,&n\neq
m.\end{array}\right.$

If there exists an OPS associated with $\mathbf{u}$, then $\mathbf{u}$ is
called quasi-definite. Positive definite moment functionals are quasi-
definite.

Observe that an OPS $\\{P_{n}(x)\\}_{n\geqslant 0}$ constitutes a basis for
$\Pi$. If for all $n\geqslant 0$, the leading coefficient of $P_{n}(x)$ is 1,
then $\\{P_{n}(x)\\}_{n\geqslant 0}$ is called a monic orthogonal polynomial
sequence (MOPS).

Given a moment functional $\mathbf{u}$ and a polynomial $q(x)$, we define the
left multiplication of $\mathbf{u}$ by $q(x)$ as the moment functional
$q\,\mathbf{u}$ such that

$\langle q\,\mathbf{u},p\rangle\,=\,\langle\mathbf{u},q\,p\rangle,\quad\forall
p\in\Pi,$

and we define the distributional derivative $D\mathbf{u}$ by

$\langle
D\mathbf{u},p\rangle\,=\,-\langle\mathbf{u},p^{\prime}\rangle,\quad\forall
p\in\Pi.$

Moreover, the product rule is satisfied, that is,

$D(q\,\mathbf{u})\,=\,q^{\prime}\,\mathbf{u}+q\,D\mathbf{u}.$

######  Definition 2.1.

Let $\mathbf{u}$ be a quasi-definite moment functional, and let
$\\{P_{n}(x)\\}_{n\geqslant 0}$ be an OPS with respect to $\mathbf{u}$. Then
$\mathbf{u}$ is classical if there are nonzero polynomials $\phi(x)$ and
$\psi(x)$ with $\deg\phi\leqslant 2$ and $\deg\psi=1$, such that $\mathbf{u}$
satisfies the distributional Pearson equation

$D(\phi\,\mathbf{u})\,=\,\psi\,\mathbf{u}.$ (2.1)

The sequence $\\{P_{n}(x)\\}_{n\geqslant 0}$ is called a classical OPS.

The following characterizations of classical moment functionals and OPS will
be of central importance in the sequel.

######  Theorem 2.2.

Let $\mathbf{u}$ be a quasi-definite moment functional, and
$\\{P_{n}(x)\\}_{n\geqslant 0}$ its associated MOPS. The following statements
are equivalent:

  1. 1.

$\mathbf{u}$ is a classical moment functional.

  2. 2.

(Bochner, [3]) There are nonzero polynomials $\phi(x)$ and $\psi(x)$ with
$\deg\phi\leqslant 2$ and $\deg\psi=1$ such that, for $n\geqslant 0$,
$P_{n}(x)$ satisfies

$\phi(x)\,P_{n}^{\prime\prime}(x)+\psi(x)\,P_{n}^{\prime}(x)\,=\,\lambda_{n}\,P_{n}(x),$
(2.2)

where $\lambda_{n}=n\,(\frac{n-1}{2}\phi^{\prime\prime}+\psi^{\prime})$.

  3. 3.

(Hahn,[10]) There is a nonzero polynomial $\phi(x)$ with $\deg\phi\leqslant
2$, such that $\left\\{\frac{P_{n+1}^{\prime}(x)}{n+1}\right\\}_{n\geqslant
0}$ is the MOPS associated with the moment functional
$\mathbf{v}=\phi(x)\,\mathbf{u}$.

  4. 4.

(First structure relation, [2]) There is a nonzero polynomial with
$\deg\phi\leqslant 2$, and real numbers $a_{n}$, $b_{n}$, $c_{n}$, $n\geqslant
1$, with $c_{n}\neq 0$, such that

$\phi(x)\,P^{\prime}_{n}(x)\,=\,a_{n}\,P_{n+1}(x)+b_{n}\,P_{n}(x)+c_{n}\,P_{n-1}(x),\quad
n\geqslant 1.$

  5. 5.

(Second structure relation, [6, 16]) There are real numbers $\alpha_{n}$ and
$\beta_{n}$, $n\geqslant 2$, such that

$P_{n}(x)\,=\,\frac{P^{\prime}_{n+1}(x)}{n+1}+\alpha_{n}\,\frac{P_{n}^{\prime}(x)}{n}+\beta_{n}\,\frac{P^{\prime}_{n-1}(x)}{n-1},\quad
n\geqslant 2.$ (2.3)

  6. 6.

(Rodrigues formula, [21]) There is a non zero polynomial $\phi(x)$ with
$\deg\phi\leqslant 2$, and a non zero real number $k_{n}\neq 0$ such that

$D^{n}(\phi^{n}(x)\,\mathbf{u})\,=\,k_{n}\,P_{n}(x)\,\mathbf{u},\quad
n\geqslant 0.$

It is well known (see [3] as well as [13]) that, up to affine transformations
of the independent variable, the only families of positive definite classical
orthogonal polynomials are the Hermite, Laguerre, and Jacobi polynomials. The
corresponding moment functionals admit an integral representation of the form

$\langle\mathbf{u},p\rangle\,=\,\int_{I}p(x)\,w(x)\,dx,\quad p\in\Pi,$

where $I=\mathbb{R}$ and $w(x)=e^{-x^{2}}$ in the Hermite case,
$I=(0,+\infty)$ and $w(x)=x^{\alpha}e^{-x}$ with $\alpha>-1$ in the Laguerre
case, and $I=(-1,1)$ and $w(x)=(1-x)^{\alpha}(1+x)^{\beta}$ with
$\alpha,\beta>-1$ in the Jacobi case. We note that in each case, $w(x)>0$ on
$I$ and, thus, we say that $w(x)$ is a weight function.

The definition of classical moment functionals in terms of the distributional
Pearson equation not only encompasses positive definite moment functionals
associated with weight functions, but includes the non positive case as well.
Considering the non positive definite case gives rise to the Bessel classical
moment functional satisfying the distributional Pearson equation (2.1) with
$\phi(x)=x^{2}$ and $\psi(x)=ax+2$. The Bessel functional is quasi-definite
when $a\neq-1,-2,\ldots$ Moreover, it has the following integral
representation

$\langle\mathbf{u},p\rangle=\int_{c}p(z)\,w(z)\,dz,\quad p\in\Pi,$

where $w(z)=(2\,\pi\,i)^{-1}z^{a-2}e^{-2/z}$, and $c$ is the unit circle
oriented in the counter-clockwise direction.

Observe that from Theorem 2.2, if $\mathbf{u}$ is a classical moment
functional satisfying (2.1), then $\mathbf{v}=\phi(x)\,\mathbf{u}$ is a
classical moment functional satisfying the Pearson equation

$D(\phi\,\mathbf{v})\,=\,(\psi+\phi^{\prime})\,\mathbf{v}.$

Iterating this idea, we get that the high-order derivatives of classical
orthogonal polynomials are again classical orthogonal polynomials of the same
type.

######  Theorem 2.3 ([10, 13, 14]).

Let $\mathbf{u}$ be a classical moment functional satisfying (2.1), and
$\\{P_{n}(x)\\}_{n\geqslant 0}$ its corresponding MOPS. For $k\geqslant 0$,
let $\mathbf{v}_{k}=\phi^{k}(x)\,\mathbf{u}$ and $\\{Q_{n,k}(x)\\}_{n\geqslant
0}$ be the sequence of polynomials given by

$Q_{n,k}(x):=\dfrac{1}{(n+1)_{k}}P_{n+k}^{(k)}(x),\quad n\geqslant 0,$ (2.4)

where $p^{(k)}$ is the $k$-th derivative of $p$, and
$(\nu)_{k}=\nu\,(\nu+1)\cdots(\nu+k-1)$, $(\nu)_{0}=1$, denotes the Pochhammer
symbol. Then, for each $k\geqslant 0$, $\\{Q_{n,k}(x)\\}_{n\geqslant 0}$ is a
MOPS associated with the moment functional $\mathbf{v}_{k}$, satisfying

$D(\phi\,\mathbf{v}_{k})\,=\,\psi_{k}\,\mathbf{v}_{k},$

where $\psi_{k}(x)=\psi(x)+k\,\phi^{\prime}(x)$. Hence, $\mathbf{v}_{k}$ is a
classical moment functional.

##  3\. Classical sequences of numbers

This section is devoted to presenting the definition of classical moment
functionals from a different approach. We start by introducing sequences of
real numbers that satisfy a second order recurrence relation and use them to
construct linear functionals defined on $\Pi$.

######  Definition 3.1.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of real numbers with
$\mu_{0}\neq 0$. Then $\\{\mu_{n}\\}_{n\geqslant 0}$ is a pre-classical
sequence if there are real numbers $a,b,c,d,e$ satisfying

$|a|+|b|+|c|>0,\quad n\,a+d\neq 0\quad n\geqslant 0,$

such that the following holds

$(n\,a+d)\,\mu_{n+1}+(n\,b+e)\,\mu_{n}+n\,c\,\mu_{n-1}=0,\quad n\geqslant 0.$
(3.1)

By convention, $\mu_{n}=0$ whenever $n<0$.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a pre-classical sequence of real
numbers. Then it is possible to define a functional $\mathbf{u}$ as

$\mu_{n}\,:=\,\langle\mathbf{u},x^{n}\rangle,\quad n\geqslant 0,$

and extend it by linearity to all polynomials, where $\mu_{n}$ is called the
$n$-th moment of $\mathbf{u}$. Therefore, we refer to $\mathbf{u}$ as a pre-
classical moment functional. Observe that the condition $n\,a+d\neq 0$,
$n\geqslant 0$, guarantees that $\mathbf{u}$ is completely defined since each
moment

$\mu_{n+1}=-\frac{1}{n\,a+d}[(n\,b+e)\,\mu_{n}+n\,c\,\mu_{n-1}],\quad
n\geqslant 0,$

is well-defined.

The recurrence relation (3.1) can be passed down to the pre-classical moment
functional associated with $\\{\mu_{n}\\}_{n\geqslant 0}$.

######  Theorem 3.2.

A sequence $\\{\mu_{n}\\}_{n\geqslant 0}$ is pre-classical if and only if
there are non zero polynomials $\phi(x)$ and $\psi(x)$ with $\deg\phi\leqslant
2$, $\deg\psi=1$, and $\frac{n}{2}\phi^{\prime\prime}+\psi^{\prime}\neq 0$ for
$n\geqslant 0$, such that the moment functional $\mathbf{u}$ defined by
$\mu_{n}=\langle\mathbf{u},x^{n}\rangle$ satisfies the distributional Pearson
equation

$D(\phi\,\mathbf{u})\,=\,\psi\,\mathbf{u}.$

###### Proof.

Suppose that $\mathbf{u}$ satisfies (2.1) with non zero
$\phi(x)=a\,x^{2}+b\,x+c$ and $\psi(x)=d\,x+e$ with $\deg\psi=1$ such that
$0\neq\frac{n}{2}\phi^{\prime\prime}+\psi^{\prime}=n\,a+d$ for $n\geqslant 0$.
Then

$\langle\mathbf{u},\phi\,p^{\prime}+\psi\,p\rangle=0,\quad\forall p\in\Pi.$

In particular,

$0=\langle\mathbf{u},n\,\phi\,x^{n-1}+\psi\,x^{n}\rangle=(n\,a+d)\,\mu_{n+1}+(n\,b+e)\,\mu_{n}+n\,c\,\mu_{n-1},\quad
n\geqslant 0.$

Therefore $\\{\mu_{n}\\}_{n\geqslant 0}$ is a pre-classical sequence of real
numbers and, thus, $\mathbf{u}$ is a pre-classical moment functional. It is
easy to verify that the implications in the opposite direction holds by
inverting each of the previous steps. ∎

For any sequence $\\{\mu_{n}\\}_{n\geqslant 0}$, we can define the sequence of
matrices $\\{G_{n}\\}_{n\geqslant 0}$ where $G_{n}$ is an $(n+1)\times(n+1)$
matrix given by $G_{0}=\mu_{0}$ and

$G_{n}=\left[\begin{array}[]{ccc|c}&&&\mu_{n}\\\ &G_{n-1}&&\vdots\\\
&&&\mu_{2n-1}\\\
\hline\cr\mu_{n}&\ldots&\mu_{2n-1}&\mu_{2n}\end{array}\right]=\begin{bmatrix}\mu_{0}&\mu_{1}&\ldots&\mu_{n}\\\
\mu_{1}&\mu_{2}&\ldots&\mu_{n+1}\\\ \vdots&\vdots&\ddots&\vdots\\\
\mu_{n}&\mu_{n+1}&\dots&\mu_{2n}\end{bmatrix},\quad n\geqslant 0.$ (3.2)

In particular, consider the sequence $\\{\mu_{n}\\}_{n\geqslant 0}$ with
$\mu_{0}=1$ and $\mu_{n}=0$ for $n\geqslant 1$. Observe that this sequence
corresponds to the linear functional $\boldsymbol{\delta}\in\Pi^{*}$, known as
the Dirac delta, defined as

$\langle\boldsymbol{\delta},p\rangle=p(0),\quad\forall p\in\Pi.$

In this case, we have $G_{0}=1$,

$G_{n}=\begin{bmatrix}1&0&\ldots&0\\\ 0&0&\ldots&0\\\
\vdots&\vdots&\ddots&\vdots\\\ 0&0&\dots&0\end{bmatrix},$

and $\det G_{n}=0$ for $n\geqslant 1$. In the sequel, we will need to exclude
this and other similar cases and, therefore, we impose that $\det G_{n}\neq 0$
for $n\geqslant 0$. Hence, we have the following definition.

######  Definition 3.3.

A pre-classical sequence $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical if the
sequence of matrices $\\{G_{n}\\}_{n\geqslant 0}$ defined as in (3.2) satisfy
$\det G_{n}\neq 0$ for $n\geqslant 0$. The moment functional defined by
$\mu_{n}=\langle\mathbf{u},x^{n}\rangle$ is called a classical moment
functional.

In the sequel, some orthogonal bases for $\mathbb{R}^{n+1}$ associated with a
classical sequence of numbers will play an important role in our study of such
sequences. Therefore, in the following section we discuss the orthogonal
structure of $\mathbb{R}^{n+1}$ induced by general sequences of numbers.

##  4\. Sequences of numbers and orthogonality

Given an integer $n\geqslant 0$, a sequence of numbers induces a bilinear form
in $\mathbb{R}^{n+1}$ whose Gram matrix is $G_{n}$ defined in (3.2). In this
section, we explore orthogonal bases of $\mathbb{R}^{n+1}$ associated with a
sequence of numbers and use it to construct bases of polynomials in $\Pi_{n}$
orthogonal with respect to its corresponding moment functional.

For $n\geqslant 0$, we will denote by $\bar{e}_{0},\ldots,\bar{e}_{n}$ the
columns of the identity matrix $I_{n+1}$, that is, $I_{n+1}=[\bar{e}_{0}\
\bar{e}_{1}\ \ldots\ \bar{e}_{n}]$. The set of column vectors
$\mathcal{E}_{n}=\\{\bar{e}_{0},\ldots,\bar{e}_{n}\\}$ is called the canonical
basis for $\mathbb{R}^{n+1}$.

######  Definition 4.1.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of numbers. For $n\geqslant
0$, $\mathfrak{B}_{n}(\cdot,\cdot)$ denotes the bilinear form defined by

$\mathfrak{B}_{n}(\bar{u},\bar{v})\,=\,\bar{u}^{\top}\,G_{n}\,\bar{v},\quad\forall\bar{u},\bar{v}\in\mathbb{R}^{n+1},$

and is called the bilinear form associated with $\\{\mu_{n}\\}_{n\geqslant 0}$
(relative to the canonical basis of $\mathbb{R}^{n+1}$).

For $n\geqslant 0$, if $\det G_{n}>0$, then $\mathfrak{B}_{n}(\cdot,\cdot)$ is
an inner product on $\mathbb{R}^{n+1}$. In this case, we can define the norm

$\|\bar{u}\|_{n}=\sqrt{\mathfrak{B}_{n}(\bar{u},\bar{u})},\quad\forall\bar{u}\in\mathbb{R}^{n+1}.$

Of course, there are many orthogonal bases for $\mathbb{R}^{n+1}$ associated
with $\mathfrak{B}_{n}$. However, we are interested in orthogonal bases
obtained from the Cholesky factorization of $G_{n}$. Recall that if $\det
G_{n}\neq 0$, there is an $(n+1)\times(n+1)$ unit lower triangular matrix
(that is, with $1$’s in its main diagonal) $S_{n}^{-1}$ and an
$(n+1)\times(n+1)$ non singular diagonal matrix $H_{n}$ such that

$G_{n}=S_{n}^{-1}\,H_{n}\,S_{n}^{-\top},$

where $S_{n}^{-\top}=(S_{n}^{-1})^{\top}$ (see Theorem 4.1.3 in [9]).
Moreover, this matrix factorization is unique. We can immediately observe that
if we write the above identity as

$S_{n}\,G_{n}\,S_{n}^{\top}\,=\,H_{n},$

then we have an orthogonality relation for the columns of $S_{n}^{\top}$ as we
show in the following theorem.

######  Theorem 4.2.

Let $G_{n}=S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the Cholesky decomposition of
$G_{n}$. Then the columns of $S_{n}^{\top}$ form an orthogonal basis for
$\mathbb{R}^{n+1}$ with respect to the bilinear form $\mathfrak{B}_{n}$.

###### Proof.

If $\bar{s}_{0},\bar{s}_{1},\ldots,\bar{s}_{n}$ are the columns of
$S_{n}^{\top}$, then

$H_{n}=S_{n}\,G_{n}\,S_{n}^{\top}=\begin{bmatrix}\bar{s}_{0}^{\top}\\\\[3.0pt]
\bar{s}_{1}^{\top}\\\ \vdots\\\
\bar{s}_{n}^{\top}\end{bmatrix}\,G_{n}\,\begin{bmatrix}\bar{s}_{0}&\bar{s}_{1}&\ldots&\bar{s}_{n}\end{bmatrix}.$

Since $H_{n}$ is a diagonal matrix and $G_{n}$ is the representation of the
bilinear form $\mathfrak{B}_{n}$, it follows that

$\mathfrak{B}_{n}(\bar{s}_{i},\bar{s}_{j})=0\quad\text{for}\quad i\neq j,$

and

$\mathfrak{B}_{n}(\bar{s}_{i},\bar{s}_{i})=h_{i},\quad i=0,1,\ldots,n,$

where $h_{i}$ is the $i$-th non zero entry of $H_{n}$, that is,
$H_{n}=\text{diag}(h_{0},h_{1},\ldots,h_{n})$ ∎

The following theorem shows that the Cholesky factorization of $G_{n}$ and
$G_{n+1}$ are related. In fact, the Cholesky factorization of $G_{n+1}$ is
obtained by bordering $S_{n}$ and $H_{n}$ with a new row and a new column. The
proof relies heavily on the expressions for $G_{n}^{-1}$ and $\det G_{n}$ in
terms of cofactors. If $(G_{n})_{i,j}$ is the $(i,j)$-cofactor of $G_{n}$ with
$i,j=0,1,\ldots,n$, then

$\det G_{n}=\sum_{k=0}^{n}\,\mu_{i+k}\,(G_{n})_{i,k},\quad
G_{n}^{-1}=\frac{1}{\det
G_{n}}\begin{bmatrix}(G_{n})_{0,0}&(G_{n})_{1,0}&\ldots&(G_{n})_{n,0}\\\
(G_{n})_{0,1}&(G_{n})_{1,1}&\ldots&(G_{n})_{n,1}\\\
\vdots&\vdots&\ddots&\vdots\\\
(G_{n})_{0,n}&(G_{n})_{1,n}&\ldots&(G_{n})_{n,n}\end{bmatrix}.$

Hereon, $\mathtt{0}$ will denote the zero matrix of appropriate size.

######  Theorem 4.3.

Let $n\geqslant 0$, and let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the Cholesky
factorization of $G_{n}$. Then,
$G_{n+1}=S_{n+1}^{-1}\,H_{n+1}\,S_{n+1}^{-\top}$ with

$S_{n+1}^{\top}=\left[\begin{array}[]{c|c}S_{n}^{\top}&\hat{s}_{n+1}\\\
\hline\cr\mathtt{0}&1\end{array}\right],\quad
H_{n+1}=\left[\begin{array}[]{c|c}H_{n}&\mathtt{0}\\\
\hline\cr\mathtt{0}&h_{n+1}\end{array}\right],$

where $[\hat{s}_{n+1}^{\top}\ 1]^{\top}$ is a vector given by the formal
determinant

$\begin{bmatrix}\hat{s}_{n+1}\\\ 1\end{bmatrix}=\frac{1}{\det
G_{n}}\det\left[\begin{array}[]{ccc|c}&&&\mu_{n+1}\\\ &G_{n}&&\vdots\\\
&&&\mu_{2n+1}\\\
\hline\cr\bar{e}_{0}&\ldots&\bar{e}_{n}&\bar{e}_{n+1}\end{array}\right],$
(4.1)

and

$h_{n+1}=\frac{\det G_{n+1}}{\det G_{n}}.$ (4.2)

###### Proof.

Observe that the columns of $[S_{n}^{\top}\ \mathtt{0}]^{\top}$ are linear
combinations of the vectors
$\bar{e}_{0},\ldots,\bar{e}_{n}\in\mathbb{R}^{n+2}$ (that is, the first $n+1$
columns of $I_{n+2}$). Then $H_{n+1}=S_{n+1}\,G_{n+1}\,S_{n+1}^{\top}$ with
$S_{n+1}^{\top}$ and $H_{n+1}$ as in (4.1) and (4.2) if and only if

$\mathfrak{B}_{n+1}([\hat{s}_{n+1}^{\top}\ 1]^{\top},\bar{e}_{i})=0,\quad
i=0,1,\ldots,n,$

and

$\mathfrak{B}_{n+1}([\hat{s}_{n+1}^{\top}\ 1]^{\top},[\hat{s}_{n+1}^{\top}\
1]^{\top})=\frac{\det G_{n+1}}{\det G_{n}}.$

The above conditions can be written as a system of linear equations:

$G_{n}\,\hat{s}_{n+1}=-\begin{bmatrix}\mu_{n+1}\\\ \vdots\\\
\mu_{2n+1}\end{bmatrix}.$

By Cramer’s rule,

$\hat{s}_{n+1}=\frac{1}{\det
G_{n}}\bigg{(}(G_{n+1})_{n+1,0}\,\bar{e}_{0}+\cdots+(G_{n+1})_{n+1,n}\,\bar{e}_{n}\bigg{)},$
(4.3)

which is the expansion of the determinant (4.1) across the last row.

Moreover, we also have

$\hat{s}_{n+1}=-G_{n}^{-1}\begin{bmatrix}\mu_{n+1}\\\ \vdots\\\
\mu_{2n+1}\end{bmatrix},$ (4.4)

and

$\displaystyle S_{n+1}\,G_{n+1}\,S_{n+1}^{\top}$
$\displaystyle=\left[\begin{array}[]{c|c}S_{n}&\mathtt{0}\\\
\hline\cr\\\\[-8.0pt]
\hat{s}_{n+1}^{\top}&1\end{array}\right]\,\left[\begin{array}[]{ccc|c}&&&\mu_{n+1}\\\
&G_{n}&&\vdots\\\ &&&\mu_{2n+1}\\\
\hline\cr\mu_{n+1}&\ldots&\mu_{2n+1}&\mu_{2n+2}\end{array}\right]\,\left[\begin{array}[]{c|c}S_{n}^{\top}&\hat{s}_{n+1}\\\
\hline\cr\mathtt{0}&1\end{array}\right]$
$\displaystyle=\left[\begin{array}[]{c|c}H_{n}&S_{n}\,B_{n}\\\
\hline\cr\\\\[-8.0pt] C_{n}\,S_{n}^{\top}&h_{n+1}\end{array}\right],$

where

$B_{n}=G_{n}\,\hat{s}_{n+1}+\begin{bmatrix}\mu_{n+1}\\\ \vdots\\\
\mu_{2n+1}\end{bmatrix},\quad
C_{n}=\hat{s}_{n+1}^{\top}\,G_{n}+\begin{bmatrix}\mu_{n+1}\\\ \vdots\\\
\mu_{2n+1}\end{bmatrix}^{\top},\quad
h_{n+1}=[\hat{s}_{n+1}^{\top}\,|\,1]\begin{bmatrix}\mu_{n+1}\\\ \vdots\\\
\mu_{2n+1}\\\ \hline\cr\mu_{2n+2}\end{bmatrix}.$

It follows from (4.4) that $B_{n}=\mathtt{0}$ and $C_{n}=\mathtt{0}$. Finally,
using (4.3) and the fact that $(G_{n+1})_{n+1,n+1}=\det G_{n}$, we obtain

$h_{n+1}=\frac{1}{\det
G_{n}}\sum_{j=0}^{n+1}\mu_{n+1+j}(G_{n+1})_{n+1,j}=\frac{\det G_{n+1}}{\det
G_{n}},$

which proves (4.2). ∎

We can reformulate the above discussion in the context of $\Pi_{n}$ as
follows. Given a sequence of numbers $\\{\mu_{n}\\}_{n\geqslant 0}$, denote by
$\mathbf{u}$ the moment functional defined by
$\mu_{n}=\langle\mathbf{u},x^{n}\rangle$. If $\det G_{n}\neq 0$, then

$\langle\mathbf{u},p\,q\rangle,\quad p,q\in\Pi_{n},$

is a (non degenerate) bilinear form defined on $\Pi_{n}$. It is easy to see
that its Gram matrix relative to the basis of monomials is $G_{n}$. Let
$G_{n}=S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the Cholesky decomposition of
$G_{n}$ ($\det G_{n}\neq 0$) and let

$S_{n}=\begin{bmatrix}1&0&0&\ldots&0\\\ s_{1,0}&1&0&\ldots&0\\\
s_{2,0}&s_{2,1}&1&\ldots&0\\\ \vdots&\vdots&\vdots&\ddots&\vdots\\\
s_{n,0}&s_{n,1}&s_{n,2}&\ldots&1\end{bmatrix}$

be the explicit expression of the matrix $S_{n}$. It follows from Theorem 4.2
that the set of polynomials $\\{P_{0}(x),P_{1}(x),\ldots,P_{n}(x)\\}$ where

$P_{k}(x)=s_{k,0}+s_{k,1}x+\cdots+s_{k,k-1}x^{k-1}+x^{k},\quad
k=0,1,2,\ldots,n,$

form an orthogonal basis for $\Pi_{n}$ with respect to $\mathbf{u}$; that is,
$\deg P_{k}=k$ and

$\langle\mathbf{u},P_{j}\,P_{i}\rangle=h_{j}\,\delta_{i,j},\quad 0\leqslant
i,j\leqslant n,$

with $h_{j}\neq 0$. Moreover, by Theorem 4.3 the polynomial

$P_{n+1}(x)=\frac{1}{\det
G_{n}}\det\left[\begin{array}[]{ccc|c}&&&\mu_{n+1}\\\ &G_{n}&&\vdots\\\
&&&\mu_{2n+1}\\\ \hline\cr 1&\ldots&x^{n}&x^{n+1}\end{array}\right],$

has degree exactly $n+1$, is orthogonal to every polynomial in $\Pi_{n}$ and

$h_{n+1}=\langle\mathbf{u},P_{n+1}^{2}\rangle=\frac{\det G_{n+1}}{\det
G_{n}}.$

##  5\. Characterizations of classical sequences

Our goal for this section is to recast Theorem 2.2 in terms of classical
sequences by shifting our point of view from the moment functional
$\mathbf{u}$ to the Gram matrix $G_{n}$ associated with a bilinear form
defined on $\Pi_{n}$. Recall that $G_{n}$ is a Hankel matrix (all of its
antidiagonals are constant). Hence, we can say that this section deals with
Hankel matrices with an additional structure: the entries of $G_{n}$ satisfy
the recurrence relation (3.1). In this way, we can extend the bilinear form to
$\Pi_{n+1}$ by constructing a new Gram matrix $G_{n+1}$ by means of bordering
$G_{n}$ with a new row and column whose entries are obtained with (3.1) from
the entries of $G_{n}$. The resulting matrix $G_{n+1}$ will also be a Hankel
matrix with the additional structure mentioned above. Consequently, it will be
possible to prove by induction that the properties satisfied by $G_{n}$ are
also satisfied by $G_{n+1}$.

The following matrices will play an important role in the sequel.

######  Definition 5.1.

Let $a,b,c,d$, and $e$ be real numbers such that $|a|+|b|+|c|>0$ and
$n\,a+d\neq 0$ for $n\geqslant 0$. For $n\geqslant 1$, we define the
$n\times(n+1)$ matrix $R_{n}$ recursively as follows:

$R_{n}=\left[\begin{array}[]{ccccc|c}&&&&&0\\\ &&&R_{n-1}&&\vdots\\\ &&&&&0\\\
\hline\cr 0&\ldots&0&(n-1)\,c&(n-1)\,b+e&(n-1)\,a+d\end{array}\right],\quad
R_{1}=\begin{bmatrix}e&d\end{bmatrix}.$

Consider the differential operator $\mathcal{R}:\Pi\rightarrow\Pi$ defined as

$\mathcal{R}[p]=\phi(x)\,p^{\prime}+\psi(x)\,p,\quad\forall p\in\Pi,$

where $\phi(x)=a\,x^{2}+b\,x+c$ and $\psi(x)=d\,x+e$. Observe that

$\mathcal{R}[x^{n}]\,=\,(n\,a+d)\,x^{n+1}+(n\,b+e)\,x^{n}+n\,c\,x^{n-1},\quad
n\geqslant 0.$

In this way, the matrix $R_{n+1}^{\top}$ is the matrix representation relative
to the basis of monomials of $\mathcal{R}$ restricted to $\Pi_{n}$.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of real numbers. Define the
vector of moments

$\mathtt{M}_{n}=\begin{bmatrix}\mu_{0}&\mu_{1}&\ldots&\mu_{n}\end{bmatrix}^{\top},\quad
n\geqslant 0.$

If $\\{\mu_{n}\\}_{n\geqslant 0}$ is pre-classical, then equation (3.1) for
$\\{\mu_{0},\ldots,\mu_{n}\\}$ can be written as

$R_{n}\,\mathtt{M}_{n}=\mathtt{0}.$

This implies that if $\mathbf{u}$ is the moment functional defined as
$\mu_{n}:=\langle\mathbf{u},x^{n}\rangle$, then by (3.1), the following holds

$\left\langle\mathbf{u},\mathcal{R}[x^{n}]\right\rangle=0,\quad n\geqslant 0.$

Now, consider the vector whose entries are the monomials in $\Pi_{n}$:

$\mathtt{X}_{n}=\begin{bmatrix}1&x&\ldots&x^{n}\end{bmatrix}^{\top}.$

If we define the $n\times(n+1)$ matrix

$N_{n}=\left[\begin{array}[]{c|c}N_{n-1}&\mathtt{0}\\\
\hline\cr\mathtt{0}&n\end{array}\right],\quad
N_{1}=\begin{bmatrix}0&1\end{bmatrix},\quad N_{0}=0,$

then

$\frac{d}{dx}\mathtt{X}_{n}\,=\,N_{n}^{\top}\,\mathtt{X}_{n-1}.$

###  5.1. Bochner-type characterization

For $n\geqslant 0$, let us introduce the $(n+1)\times(n+1)$ matrix
$\mathtt{D}_{n}$ defined as

$\mathtt{D}_{0}=0,\qquad\mathtt{D}_{n}:=R_{n}^{\top}N_{n},\quad n\geqslant 1.$

It is possible to express $\mathtt{D}_{n}$ recursively as follows:

$\displaystyle\mathtt{D}_{n}=\left[\begin{array}[]{ccc|c}&&&0\\\ &&&\vdots\\\
&\mathtt{D}_{n-1}&&0\\\ &&&n\,(n-1)\,c\\\\[6.0pt]
&&&n\,((n-1)\,b+e)\\\\[6.0pt] \hline\cr
0&\ldots&0&n\,((n-1)\,a+d)\end{array}\right],\quad n\geqslant 1.$

Observe that $\mathtt{D}_{n}$ is the matrix representation relative the basis
of monomials of the operator

$\mathcal{D}[p]:=\mathcal{R}\left[\frac{d}{dx}p\right]=\phi(x)\,p^{\prime\prime}+\psi(x)\,p^{\prime},\quad\forall
p\in\Pi,$ (5.1)

restricted to $\Pi_{n}$.

In the following theorem, we show that $\mathtt{D}_{n}$ is a self-adjoint
matrix with respect to the bilinear form $\mathfrak{B}_{n}$ given in
Definition 4.1.

######  Theorem 5.2.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence satisfying (3.1),
and let $\mathfrak{B}_{n}$ be the operator defined in Definition 4.1. Then,
for $n\geqslant 0$, the matrix $\mathtt{D}_{n}$ satisfies

$\mathfrak{B}_{n}(\mathtt{D}_{n}\bar{u},\bar{v})=\mathfrak{B}_{n}(\bar{u},\mathtt{D}_{n}\bar{v}),\quad\bar{u},\bar{v}\in\mathbb{R}^{n+1}.$
(5.2)

###### Proof.

Observe that proving (5.2) is equivalent to proving

$\mathtt{D}_{n}^{\top}\,G_{n}=G_{n}\,\mathtt{D}_{n}.$

We prove this for $n\geqslant 0$ by induction.

It is obvious that $\mathtt{D}_{0}^{\top}\,G_{0}=G_{0}\,\mathtt{D}_{0}$. We
also prove the case $n=1$ for the sake of clarity since it is the first non
trivial case. We compute

$\mathtt{D}_{1}^{\top}\,G_{1}=\left[\begin{array}[]{c|c}0&0\\\ \hline\cr
e&d\end{array}\right]\left[\begin{array}[]{c|c}\mu_{0}&\mu_{1}\\\
\hline\cr\mu_{1}&\mu_{2}\end{array}\right]=\left[\begin{array}[]{c|c}\mathtt{D}_{0}^{\top}&0\\\
\hline\cr e&d\end{array}\right]\left[\begin{array}[]{c|c}G_{0}&\mu_{1}\\\
\hline\cr\mu_{1}&\mu_{2}\end{array}\right].$

Again, multiplying by blocks, we get

$\mathtt{D}_{1}^{\top}\,G_{1}=\left[\begin{array}[]{c|c}\mathtt{D}_{0}^{\top}\,G_{0}&\mathtt{D}_{0}^{\top}\,\mu_{1}\\\
\hline\cr
d\,\mu_{1}+e\,\mu_{0}&e\,\mu_{1}+d\,\mu_{2}\end{array}\right]=\left[\begin{array}[]{c|c}G_{0}\,\mathtt{D}_{0}^{\top}&0\\\
\hline\cr d\,\mu_{1}+e\,\mu_{0}&e\,\mu_{1}+d\,\mu_{2}\end{array}\right].$

Since $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence, by condition
(3.1) with $n=0$, we have

$\mu_{1}\,\mathtt{D}_{0}=0=d\,\mu_{1}+e\,\mu_{0}.$

Therefore, we can write

$\mathtt{D}_{1}^{\top}\,G_{1}=\left[\begin{array}[]{c|c}G_{0}\,\mathtt{D}_{0}^{\top}&e\,\mu_{0}+d\,\mu_{1}\\\
\hline\cr\mu_{1}\,\mathtt{D}_{0}&e\,\mu_{1}+d\,\mu_{2}\end{array}\right]=\left[\begin{array}[]{c|c}G_{0}&\mu_{1}\\\
\hline\cr\mu_{1}&\mu_{2}\end{array}\right]\left[\begin{array}[]{c|c}\mathtt{D}_{0}&e\\\
\hline\cr 0&d\end{array}\right]=G_{1}\,\mathtt{D}_{1}.$

This proves $\mathtt{D}_{1}^{\top}\,G_{1}=G_{1}\,\mathtt{D}_{1}$.

Now, suppose that $\mathtt{D}_{k}^{\top}\,G_{k}=G_{k}\,\mathtt{D}_{k}$ holds
for some $k\geqslant 0$. We compute

$\displaystyle\mathtt{D}_{k+1}^{\top}\,G_{k+1}$
$\displaystyle=\left[\begin{array}[]{ccccc|c}&&&&&0\\\
&&&\mathtt{D}_{k}^{\top}&&\vdots\\\ &&&&&0\\\ \hline\cr
0&\ldots&0&(k+1)\,k\,c&(k+1)\,(k\,b+e)&(k+1)\,(k\,a+d)\end{array}\right]$
$\displaystyle\qquad\times\left[\begin{array}[]{ccc|c}&&&\mu_{k+1}\\\
&G_{k}&&\vdots\\\ &&&\mu_{2k+1}\\\
\hline\cr\mu_{k+1}&\ldots&\mu_{2k+1}&\mu_{2k+2}\end{array}\right].$

Multiplying by blocks, we get

$\mathtt{D}_{k+1}^{\top}\,G_{k+1}=\left[\begin{array}[]{c|c}\mathtt{D}_{k}^{\top}\,G_{k}&\bar{x}_{k}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt] \bar{y}_{k}^{\top}&\bar{z}_{k}\end{array}\right],$

where

$\bar{x}_{k}=\mathtt{D}_{k}^{\top}\begin{bmatrix}\mu_{k+1}\\\ \vdots\\\
\mu_{2k+1}\end{bmatrix},$ $\bar{y}_{k}^{\top}=\left[(k+1)\,k\,c\ \ \
(k+1)\,(k\,b+e)\ \ \
(k+1)\,(k\,a+d)\right]\left[\begin{array}[]{ccc}\mu_{k-1}&\ldots&\mu_{2k-1}\\\
\mu_{k}&\ldots&\mu_{2k}\\\
\hline\cr\mu_{k+1}&\ldots&\mu_{2k+1}\end{array}\right],$

and

$\bar{z}_{k}=(k+1)\,k\,c\,\mu_{2k}+(k+1)\,(k\,b+e)\,\mu_{2k+1}+(k+1)\,(k\,a+d)\,\mu_{2k+2}=\bar{z}_{k}^{\top}.$

Since by our induction hypothesis we have

$\mathtt{D}_{k+1}^{\top}\,G_{k+1}=\left[\begin{array}[]{c|c}G_{k}\,\mathtt{D}_{k}&\bar{x}_{k}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt]
\bar{y}_{k}^{\top}&\bar{z}_{k}^{\top}\end{array}\right],$

it is clear that our main efforts should focus on showing that
$\bar{x}_{k}=\bar{y}_{k}$. Observe that

$\bar{x}_{k}=\left[\begin{array}[]{ccccc|c}&&&&&0\\\
&&&\mathtt{D}_{k-1}^{\top}&&\vdots\\\ &&&&&0\\\ \hline\cr
0&\ldots&0&k\,(k-1)\,c&k\,((k-1)\,b+e)&k\,((k-1)\,a+d)\end{array}\right]\begin{bmatrix}\mu_{k+1}\\\
\vdots\\\ \mu_{2k+1}\end{bmatrix}.$

Then, the $k$-th entry of $\bar{x}_{k}$ is

$k\,(k-1)\,c\,\mu_{2k-1}+k\,((k-1)\,b+e)\,\mu_{2k}+k\,((k-1)\,a+d)\,\mu_{2k+1},$

while the $k$-th entry of $\bar{y}_{k}$ is

$(k+1)\,k\,c\,\mu_{2k-1}+(k+1)\,(k\,b+e)\,\mu_{2k}+(k+1)\,(k\,a+d)\,\mu_{2k+1},$

which means that the $k$-th entry of $\bar{y}_{k}-\bar{x}_{k}$ is

$(2k\,a+d)\,\mu_{2k+1}+(2k\,b+e)\,\mu_{2k}+2k\,c\,\mu_{2k-1}=0,$

where we have used condition (3.1).

Now, noticing that the last row of $\mathtt{D}_{k-1}^{\top}$ is

$\begin{bmatrix}0&\ldots&0&(k-1)\,(k-2)\,c&(k-1)\,((k-2)\,b+2)&(k-1)\,((k-2)\,a+d)\end{bmatrix},$

we have that the $(k-1)$-th entry of $\bar{x}_{k}$ is

$(k-1)\,(k-2)\,c\,\mu_{2k-2}+(k-1)\,((k-2)\,b+2)\,\mu_{2k-1}+(k-1)\,((k-2)\,a+d)\,\mu_{2k},$

while the $(k-1)$-th entry of $\bar{y}_{k}$ is

$(k+1)\,k\,c\,\mu_{2k-2}+(k+1)\,(k\,b+e)\,\mu_{2k-1}+(k+1)\,(k\,a+d)\,\mu_{2k}.$

Then the $(k-1)$-th entry of $\bar{y}_{k}-\bar{x}_{k}$ is

$2\bigg{[}((2k-1)\,a+d)\,\mu_{2k}+((2k-1)\,b+e)\,\mu_{2k-1}+(2k-1)\,c\,\mu_{2k-2}\bigg{]}=0,$

where, again, we have used (3.1) with $n=k-1$.

If we continue in this way, then for $i=0,1,\ldots,k$, the $(i+1)$-th entry of
$\bar{y}_{k}-\bar{x}_{k}$ is

$(k+1-i)\bigg{[}((k+i)\,a+d)\mu_{k+1+i}+((k+i)\,b+e)\,\mu_{k+i}+(k+i)\,c\,\mu_{k-1+i}\bigg{]}=0.$

This means that $\bar{x}_{k}=\bar{y}_{k}$ and, consequently,

$\mathtt{D}_{k+1}^{\top}\,G_{k+1}=\left[\begin{array}[]{c|c}G_{k}\,\mathtt{D}_{k}&\bar{y}_{k}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt]
\bar{x}_{k}^{\top}&\bar{z}_{k}^{\top}\end{array}\right]=G_{k+1}\,\mathtt{D}_{k+1},$

which proves that $\mathtt{D}_{n}^{\top}\,G_{n}=G_{n}\,\mathtt{D}_{n}$ holds
for $n\geqslant 0$. ∎

Let us now consider the eigenvectors of $\mathtt{D}_{n}$. Suppose that
$\bar{u},\bar{v}\in\mathbb{R}^{n+1}$ are eigenvectors corresponding to
distinct eigenvalues $\lambda$ and $\overline{\lambda}$, respectively. Then,

$\lambda\,\mathfrak{B}_{n}(\bar{u},\bar{v})=\mathfrak{B}_{n}(\lambda\,\bar{u},\bar{v})=\mathfrak{B}_{n}(\mathtt{D}_{n}\bar{u},\bar{v})=\mathfrak{B}_{n}(\bar{u},\mathtt{D}_{n}\bar{v})=\mathfrak{B}_{n}(\bar{u},\overline{\lambda}\,\bar{v})=\overline{\lambda}\,\mathfrak{B}_{n}(\bar{u},\bar{v}),$

which implies

$(\lambda-\overline{\lambda})\,\mathfrak{B}_{n}(\bar{u},\bar{v})=0.$

Since $\lambda\neq\overline{\lambda}$, we must have that
$\mathfrak{B}_{n}(\bar{u},\bar{v})=0$ and, therefore, $\bar{u}$ and $\bar{v}$
are orthogonal with respect to $\mathfrak{B}_{n}$. We have already encountered
the eigenvectors of $\mathtt{D}_{n}$, as the following theorem shows. This
theorem is, in fact, a characterization of classical sequences in terms of
$\mathtt{D}_{n}$.

######  Theorem 5.3.

For $n\geqslant 0$, let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the Cholesky
factorization of $G_{n}$ and let
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns of
$S_{n}^{\top}$. Then $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence if
and only if

$\mathtt{D}_{n}\,\bar{s}_{n,j}\,=\,\lambda_{j}\,\bar{s}_{n,j},\quad n\geqslant
0,\quad 0\leqslant j\leqslant n,$

where $\lambda_{j}=j\,[(j-1)\,a+d]$.

###### Proof.

Suppose that $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence. For
$n=0$, then $\mathtt{D}_{0}=0$ and $\bar{s}_{0,0}=1$. Thus it is obvious that
$\mathtt{D}_{0}\,\bar{s}_{0,0}=\lambda_{0}\,\bar{s}_{0,0}$.

For $n=1$,

$\mathtt{D}_{1}=\left[\begin{array}[]{c|c}\mathtt{D}_{0}&e\\\ \hline\cr
0&d\end{array}\right],\quad\bar{s}_{1,0}=\left[\begin{array}[]{c}\bar{s}_{0,0}\\\
\hline\cr 0\end{array}\right].$

Then, multiplying by blocks, we have

$\mathtt{D}_{1}\,\bar{s}_{1,0}=\left[\begin{array}[]{c}\mathtt{D}_{0}\,\bar{s}_{0,0}\\\
\hline\cr
0\end{array}\right]=\lambda_{0}\,\left[\begin{array}[]{c}\bar{s}_{0,0}\\\
\hline\cr 0\end{array}\right]=\lambda_{0}\,\bar{s}_{1,0}.$

Since $S_{1}^{\top}$ is upper triangular with $1$’s on its diagonal, then
$\\{\bar{s}_{1,0},\bar{s}_{1,1}\\}$ constitutes a basis for $\mathbb{R}^{2}$.
Then we can write

$\mathtt{D}_{1}\,\bar{s}_{1,1}=a_{1,1}\,\bar{s}_{1,1}+a_{1,0}\,\bar{s}_{1,0},$

for some constants $a_{1,1}$ and $a_{1,0}$. Using the orthogonality of the
columns of $S_{1}^{\top}$ with respect to $\mathfrak{B}_{1}$ and Theorem 5.2,
we obtain

$a_{1,j}=\frac{\mathfrak{B}_{1}(\mathtt{D}_{1}\,\bar{s}_{1,1},\bar{s}_{1,j})}{h_{j}}=\frac{\mathfrak{B}_{1}(\bar{s}_{1,1},\mathtt{D}_{1}\,\bar{s}_{1,j})}{h_{j}},\quad
j=0,1.$

It follows that $a_{1,0}=0$ and, thus,

$\mathtt{D}_{1}\,\bar{s}_{1,1}=a_{1,1}\,\bar{s}_{1,1}=\begin{bmatrix}*\\\
a_{1,1}\end{bmatrix},$

where we have taken into account that $\bar{s}_{1,1}=[*\ 1]^{\top}$ (the value
of the entry denoted by $*$ has no relevance). If we multiply by blocks, we
get

$\mathtt{D}_{1}\,\bar{s}_{1,1}=\begin{bmatrix}*\\\ d\end{bmatrix},$

which implies that $a_{1,1}=d=\lambda_{1}$.

Now, suppose that

$\mathtt{D}_{k}\,\bar{s}_{k,j}\,=\,\lambda_{j}\,\bar{s}_{k,j},\quad 0\leqslant
j\leqslant k,$

holds for some $k\geqslant 0$. Recall that

$\displaystyle\mathtt{D}_{k+1}=\left[\begin{array}[]{ccc|c}&&&0\\\
&&&\vdots\\\ &\mathtt{D}_{k}&&0\\\ &&&(k+1)\,k\,c\\\\[6.0pt]
&&&(k+1)\,(k\,b+e)\\\\[6.0pt] \hline\cr
0&\ldots&0&(k+1)\,(k\,a+d)\end{array}\right],$

and that, by Theorem (4.3),

$\bar{s}_{k+1,j}=\left[\begin{array}[]{c}\bar{s}_{k,j}\\\ \hline\cr
0\end{array}\right],\quad 0\leqslant j\leqslant
k,\quad\text{and}\quad\bar{s}_{k+1,k+1}=\left[\begin{array}[]{c}*\\\ \hline\cr
1\end{array}\right],$

where the values of the entries denoted by $*$ are not relevant here.
Multiplying by blocks and using the induction hypothesis, we get

$\mathtt{D}_{k+1}\,\bar{s}_{k+1,j}=\left[\begin{array}[]{c}\mathtt{D}_{k}\,\bar{s}_{k,j}\\\
\hline\cr 0\end{array}\right]=\lambda_{j}\,\bar{s}_{k+1,j},\quad 0\leqslant
j\leqslant k.$

Since $S_{k+1}^{\top}$ is an upper triangular matrix with $1$’s on its
diagonal, its columns constitute a basis for $\mathbb{R}^{k+1}$. Then, we can
write

$\mathtt{D}_{k+1}\,\bar{s}_{k+1,k+1}=\sum_{j=0}^{k+1}\,a_{k+1,j}\,\bar{s}_{k+1,j},$

where, by the orthogonality of the columns of $S_{k+1}^{\top}$ with respect to
$\mathfrak{B}_{k+1}$ and Theorem 5.2, we have

$\displaystyle a_{k+1,j}$
$\displaystyle=\frac{\mathfrak{B}_{k+1}(\mathtt{D}_{k+1}\,\bar{s}_{k+1,k+1},\bar{s}_{k+1,j})}{h_{j}}$
$\displaystyle=\frac{\mathfrak{B}_{k+1}(\bar{s}_{k+1,k+1},\mathtt{D}_{k+1}\,\bar{s}_{k+1,j})}{h_{j}},\quad
0\leqslant j\leqslant k+1.$

It follows that

$\displaystyle
a_{k+1,j}=\,\lambda_{j}\,\frac{\mathfrak{B}_{k+1}(\bar{s}_{k+1,k+1},\bar{s}_{k+1,j})}{h_{j}}\,=\,0,\quad
0\leqslant j\leqslant k,$

and, consequently,

$\mathtt{D}_{k+1}\,\bar{s}_{k+1,k+1}=a_{k+1,k+1}\,\bar{s}_{k+1,k+1}=\left[\begin{array}[]{c}*\\\
\hline\cr a_{k+1,k+1}\end{array}\right].$

Moreover, if we multiply by blocks, we get

$\mathtt{D}_{k+1}\,\bar{s}_{k+1,k+1}=\left[\begin{array}[]{c}*\\\
\hline\cr\lambda_{k+1}\end{array}\right],$

which implies that $a_{k+1,k+1}=\lambda_{k+1}$. Then the sufficient condition
is proved by the principle of induction.

Conversely, if $\mathfrak{B}_{n}$ is the bilinear form defined in Definition
4.1, then

$\mathfrak{B}_{n}(\mathtt{D}_{n}\,\bar{s}_{n,j},\bar{s}_{n,k})\,=\,\mathfrak{B}_{n}(\bar{s}_{n,j},\mathtt{D}_{n}\,\bar{s}_{n,k}),\quad
n\geqslant 0,\quad 0\leqslant j,k\leqslant n.$

Indeed,

$\displaystyle\mathfrak{B}_{n}(\mathtt{D}_{n}\,\bar{s}_{n,j},\bar{s}_{n,k})-\mathfrak{B}_{n}(\bar{s}_{n,j},\mathtt{D}_{n}\,\bar{s}_{n,k})=\,(\lambda_{j}-\lambda_{k})\,\mathfrak{B}_{n}(\bar{s}_{n,j},\bar{s}_{n,k})=\,(\lambda_{j}-\lambda_{k})\,h_{j}\,\delta_{j,k},$

which vanishes for all $0\leqslant j,k\leqslant n$. This implies that

$S_{n}\,(\mathtt{D}_{n}^{\top}\,G_{n}-G_{n}\,\mathtt{D}_{n})\,S_{n}^{\top}\,=\,\mathtt{0},\quad
n\geqslant 0,$

or, equivalently,

$\mathtt{D}_{n}^{\top}\,G_{n}-G_{n}\,\mathtt{D}_{n}\,=\,\mathtt{0}.$

Since $\mathtt{D}_{n}=R_{n}^{\top}N_{n}$, the first column of the above matrix
identity reads $N_{n}^{\top}\,R_{n}\,\mathtt{M}_{n}=\mathtt{0}$, or,
equivalently,

$k\,\left[\big{(}(k-1)\,a+d\big{)}\,\mu_{k}+\big{(}(k-1)\,b+e\big{)}\,\mu_{k-1}+(k-1)\,\mu_{k-2}\right]=0,$

for $n\geqslant 0$ and $0\leqslant k\leqslant n$. It follows that
$\\{\mu_{n}\\}_{n\geqslant 0}$ satisfies (3.1); hence, it is pre-classical.
Furthermore, $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical since, otherwise,
$G_{n}$ would not have a Cholesky factorization for some $n\geqslant 0$. ∎

The above results can be passed down to the operator $\mathcal{D}$ defined in
(5.1). Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence of real
numbers, and let $\mathbf{u}$ be the moment functional defined as
$\mu_{n}=\langle\mathbf{u},x^{n}\rangle$, $n\geqslant 0$. Then (5.2) implies
that

$\left\langle\mathbf{u},\mathcal{D}[p]\,q\right\rangle\,=\,\left\langle\mathbf{u},p\,\mathcal{D}[q]\right\rangle,\quad\forall
p,q\in\Pi.$

That is, $\mathcal{D}$ is a self-adjoint operator on polynomials. Moreover,
for $n\geqslant 0$, let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the Cholesky
factorization of $G_{n}$ and let
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns of
$S_{n}^{\top}$. From Theorem 5.3 we deduce that the sequence of polynomials
$\\{P_{n}\\}_{n\geqslant 0}$ with

$P_{n}(x)\,=\,\bar{s}_{n,n}^{\top}\,\mathtt{X}_{n},\quad n\geqslant 0,$

are eigenfunctions of the operator $\mathcal{D}$. That is,

$\mathcal{D}[P_{n}]\,=\,\lambda_{n}\,P_{n},\quad n\geqslant 0,$

with $\lambda_{n}=n\,[(n-1)\,a+d]$. Note that $\\{P_{n}\\}_{n\geqslant 0}$ is
a sequence of polynomials orthogonal with respect to $\mathbf{u}$.

###  5.2. Hahn-type characterization

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of real numbers and let
$a,b,c\in\mathbb{R}$ such that $|a|+|b|+|c|>0$. We can define a new sequence
$\\{\sigma_{n}\\}_{n\geqslant 0}$ as follows:

$\sigma_{n}=a\,\mu_{n+2}+b\,\mu_{n+1}+c\,\mu_{n},\quad n\geqslant 0.$

Notice that if $\mathbf{u}$ is the moment functional defined as
$\mu_{n}=\langle\mathbf{u},x^{n}\rangle$, then $\\{\sigma_{n}\\}_{n\geqslant
0}$ is the sequence of moments of the functional given by
$\mathbf{v}=\phi(x)\,\mathbf{u}$ where $\phi(x)=a\,x^{2}+b\,x+c$. Indeed, for
$n\geqslant 0$,

$\langle\mathbf{v},x^{n}\rangle=\langle\mathbf{u},\phi\,x^{n}\rangle=\langle\mathbf{u},a\,x^{n+2}+b\,x^{n+1}+c\,x^{n}\rangle=a\mu_{n+2}+b\mu_{n+1}+c\mu_{n}=\sigma_{n}.$

We denote by $\\{G_{n}^{(1)}\\}_{n\geqslant 0}$ the sequence of
$(n+1)\times(n+1)$ matrices with

$G_{0}^{(1)}=\sigma_{0},\quad\text{and}\quad
G_{n}^{(1)}=\left[\begin{array}[]{ccc|c}&&&\sigma_{n}\\\
&G_{n-1}^{(1)}&&\vdots\\\ &&&\sigma_{2n-1}\\\
\hline\cr\sigma_{n}&\ldots&\sigma_{2n-1}&\sigma_{2n}\end{array}\right],\quad
n\geqslant 1.$ (5.3)

The following theorem shows that the pre-classical character is inherited by
$\\{\sigma_{n}\\}_{n\geqslant 0}$.

######  Theorem 5.4.

If $\\{\mu_{n}\\}_{n\geqslant 0}$ is pre-classical satisfying (3.1), then
$\\{\sigma_{n}\\}_{n\geqslant 0}$ is pre-classical satisfying

$(n\,a+d_{1})\,\sigma_{n+1}+(n\,b+e_{1})\,\sigma_{n}+n\,c\,\sigma_{n-1}\,=\,0,\quad
n\geqslant 0,$

where $d_{1}=d+2\,a$ and $e_{1}=e+b$. Moreover,

$\sigma_{n}\,=\,-(n\,a+d)\,\mu_{n+2}-(n\,b+e)\,\mu_{n+1}-n\,c\,\mu_{n},\quad
n\geqslant 0.$ (5.4)

###### Proof.

For $n\geqslant 0$, we compute

$\displaystyle(n\,a+d_{1})\,\sigma_{n+1}+(n\,b+e_{1})\,\sigma_{n}+n\,c\,\sigma_{n-1}$
$\displaystyle=a\,\big{(}[(n+2)a+d]\mu_{n+3}+[(n+2)b+e]\mu_{n+2}+(n+2)c\mu_{n+1}-b\,\mu_{n+2}-2\,c\,\mu_{n+1}\big{)}$
$\displaystyle\ \
+b\,\big{(}[(n+1)a+d]\mu_{n+2}+[(n+1)b+e]\mu_{n+1}+(n+1)c\,\mu_{n}+a\,\mu_{n+2}-c\,\mu_{n}\big{)}$
$\displaystyle\ \ \
+c\,\big{(}(n\,a+d)\,\mu_{n+1}+(n\,b+e)\,\mu_{n}+n\,c\,\mu_{n-1}+2\,a\,\mu_{n+1}+b\,\mu_{n}\big{)},$

where we have used $\sigma_{n}=a\,\mu_{n+2}+b\,\mu_{n+1}+c\,\mu_{n}$. By
(3.1), we have

$\displaystyle(n\,a+d_{1})$
$\displaystyle\,\sigma_{n+1}+(n\,b+e_{1})\,\sigma_{n}+n\,c\,\sigma_{n-1}$
$\displaystyle=-a\,b\,\mu_{n+2}-2\,a\,c\,\mu_{n+1}+a\,b\,\mu_{n+2}-b\,c\,\mu_{n}+2\,a\,c\,\mu_{n+1}+b\,c\,\mu_{n}=0.$

Finally, (5.4) follows from the fact that (3.1) can be written as

$a\,\mu_{n+2}+b\,\mu_{n+1}+c\,\mu_{n}\,=\,-(n\,a+d)\,\mu_{n+2}-(n\,b+e)\,\mu_{n+1}-n\,c\,\mu_{n}.$

∎

When $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence of real numbers,
the matrix $G_{n}^{(1)}$ satisfies an interesting and useful relation
involving the matrices $G_{n}$ and $\mathtt{D}_{n}$.

######  Proposition 5.5.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence of real numbers
satisfying (3.1). Then, for $n\geqslant 0$,

$N_{n+1}^{\top}\,G_{n}^{(1)}\,N_{n+1}\,=\,-\mathtt{D}_{n+1}^{\top}\,G_{n+1}.$

###### Proof.

We prove this theorem by induction. For $n=0$, on one hand we have

$N_{1}^{\top}\,G_{0}^{(1)}\,N_{1}\,=\,\begin{bmatrix}0\\\
1\end{bmatrix}\,\sigma_{0}\,\begin{bmatrix}0&1\end{bmatrix}\,=\,\begin{bmatrix}0&0\\\
0&\sigma_{0}\end{bmatrix}.$

On the other hand, we have

$-\mathtt{D}_{1}^{\top}\,G_{1}\,=\,\begin{bmatrix}0&0\\\
-e&-d\end{bmatrix}\,\begin{bmatrix}\mu_{0}&\mu_{1}\\\
\mu_{1}&\mu_{2}\end{bmatrix}\,=\,\begin{bmatrix}0&0\\\
-e\,\mu_{0}-d\,\mu_{1}&-e\,\mu_{1}-d\,\mu_{2}\end{bmatrix}.$

Since $\\{\mu_{n}\\}_{n\geqslant 0}$ satisfies (3.1), we have that
$-e\,\mu_{0}-d\,\mu_{1}=0$ and, by (5.4), $-e\,\mu_{1}-d\,\mu_{2}=\sigma_{0}$.
Therefore,

$-\mathtt{D}_{1}^{\top}\,G_{1}\,=\,\begin{bmatrix}0&0\\\
0&\sigma_{0}\end{bmatrix},$

which proves that
$N_{1}^{\top}\,G_{0}^{(1)}\,N_{1}=-\mathtt{D}_{1}^{\top}\,G_{1}$.

Now, suppose that
$N_{k+1}^{\top}\,G_{k}^{(1)}\,N_{k+1}\,=\,-\mathtt{D}_{k+1}^{\top}\,G_{k+1}$
holds for $k\geqslant 0$. On one hand, we compute

$\displaystyle N_{k+2}^{\top}\,$ $\displaystyle G_{k+1}^{(1)}\,N_{k+2}$
$\displaystyle=\,\left[\begin{array}[]{c|c}N_{k+1}^{\top}&\mathtt{0}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt]
\mathtt{0}&k+2\end{array}\right]\,\left[\begin{array}[]{ccc|c}&&&\sigma_{k+1}\\\
&G_{k}^{(1)}&&\vdots\\\ &&&\sigma_{2k+1}\\\
\hline\cr\sigma_{k+1}&\ldots&\sigma_{2k+1}&\sigma_{2k+2}\end{array}\right]\,\left[\begin{array}[]{c|c}N_{k+1}&\mathtt{0}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt] \mathtt{0}&k+2\end{array}\right].$

Multiplying by blocks and using the induction hypothesis, we get

$N_{k+2}^{\top}\,G_{k+1}^{(1)}\,N_{k+2}\,=\,\left[\begin{array}[]{c|c}-\mathtt{D}_{k+1}^{\top}\,G_{k+1}&\bar{x}_{k}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt]
\bar{x}_{k}^{\top}&(k+2)^{2}\,\sigma_{2k+2}\end{array}\right],$

where

$\bar{x}_{k}\,=\,(k+2)\,N_{k+1}^{\top}\,\begin{bmatrix}\sigma_{k+1}\\\
\vdots\\\ \sigma_{2k+1}\end{bmatrix}.$

On the other hand,

$\displaystyle-\mathtt{D}_{k+2}^{\top}\,G_{k+2}$
$\displaystyle=\left[\begin{array}[]{ccccc|c}&&&&&0\\\
&&&-\mathtt{D}_{k+1}^{\top}&&\vdots\\\ &&&&&0\\\ \hline\cr
0&\ldots&0&-(k+2)\,(k+1)\,c&-(k+2)\,[(k+1)\,b+e]&-(k+2)\,[(k+1)\,a+d]\end{array}\right]$
$\displaystyle\qquad\times\left[\begin{array}[]{ccc|c}&&&\mu_{k+2}\\\
&G_{k+1}&&\vdots\\\ &&&\mu_{2k+3}\\\
\hline\cr\mu_{k+2}&\ldots&\mu_{2k+3}&\mu_{2k+4}\end{array}\right].$

After multiplying by blocks, we have

$-\mathtt{D}_{k+2}^{\top}\,G_{k+2}\,=\,\left[\begin{array}[]{c|c}-\mathtt{D}_{k+1}^{\top}\,G_{k+1}&\bar{y}_{k}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt] \bar{z}_{k}&w_{k}\end{array}\right],$

where

$\displaystyle\bar{y}_{k}=-\mathtt{D}_{k+1}^{\top}\,\begin{bmatrix}\mu_{k+2}\\\
\vdots\\\ \mu_{2k+3}\end{bmatrix},$
$\displaystyle\bar{z}_{k}=\begin{bmatrix}0&\ldots&0&-(k+2)\,(k+1)\,c&-(k+2)\,[(k+1)\,b+e]\end{bmatrix}\,G_{k+1}$
$\displaystyle\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad-(k+2)\,[(k+1)\,a+d]\,\begin{bmatrix}\mu_{k+2}&\ldots&\mu_{2k+3}\end{bmatrix},$
$\displaystyle
w_{k}=-(k+2)\bigg{[}(k+1)\,c\,\mu_{2k+2}+[(k+1)\,b+e]\,\mu_{2k+3}+[(k+1)\,a+d]\,\mu_{2k+4}\bigg{]}.$

By Theorem 5.2, we have that $\bar{z}_{k}=\bar{y}_{k}^{\top}$. Hence, our
efforts should focus on showing that $\bar{x}_{k}=\bar{y}_{k}$ and
$w_{k}=(k+2)^{2}\,\sigma_{2k+2}$. Let us start by proving the second identity.

Observe that

$\displaystyle w_{k}=$
$\displaystyle-(k+2)\,\bigg{[}[(2k+2)\,a+d]\,\mu_{2k+4}+[(2k+2)\,b+e]\,\mu_{2k+3}+(2k+2)\,c\,\mu_{2k+2}\bigg{]}$
$\displaystyle+(k+2)\,(k+1)\,(a\,\mu_{2k+4}+b\,\mu_{2k+3}+c\,\mu_{2k+2}).$

From (5.4) and the fact that
$\sigma_{2k+2}=a\,\mu_{2k+4}+b\,\mu_{2k+3}+c\,\mu_{2k+2}$, we deduce that

$w_{k}=(k+2)\,\sigma_{2k+2}+(k+2)\,(k+1)\,\sigma_{2k+2}=(k+2)^{2}\,\sigma_{2k+2}.$

In order to prove that $\bar{x}_{k}=\bar{y}_{k}$, we note that

$\bar{x}_{k}\,=\,(k+2)\,\left[\begin{array}[]{c|c}N_{k}^{\top}&\mathtt{0}\\\\[3.0pt]
\hline\cr\\\\[-6.0pt]
\mathtt{0}&k+1\end{array}\right]\,\begin{bmatrix}\sigma_{k+1}\\\ \vdots\\\
\sigma_{2k+1}\end{bmatrix},$

and

$\bar{y}_{k}=\left[\begin{array}[]{ccccc|c}&&&&&0\\\
&&&-\mathtt{D}_{k}^{\top}&&\vdots\\\ &&&&&0\\\ \hline\cr
0&\ldots&0&-(k+1)\,k\,c&-(k+1)\,(k\,b+e)&-(k+1)\,(k\,a+d)\end{array}\right]\,\begin{bmatrix}\mu_{k+2}\\\
\vdots\\\ \mu_{2k+3}\end{bmatrix}.$

Then, the last entry of $\bar{x}_{k}-\bar{y}_{k}$ is

$\displaystyle(k+2)\,(k+1)\,\sigma_{2k+1}+(k+1)\,[k\,c\,\mu_{2k+1}+(k\,b+e)\,\mu_{2k+2}+(k\,a+d)\,\mu_{2k+3}]$
$\displaystyle=(k+2)\,(k+1)\,\sigma_{2k+1}-(k+1)\,\sigma_{2k+1}-(k+1)^{2}\,\sigma_{2k+1}=0,$

where we have used (5.4) and the fact that
$\sigma_{2k+1}=a\,\mu_{2k+3}+b\,\mu_{2k+2}+c\,\mu_{2k+1}$. In general, the
$i$-th entry of $\bar{x}_{k}-\bar{y}_{k}$, with $1\leqslant i\leqslant k+2$,
is

$\displaystyle(i-1)\,\bigg{[}(k+2)\,\sigma_{k+i-1}+[(i-2)\,c\,\mu_{k+i-1}+[(i-2)\,b+e]\,\mu_{k+i}+[(i-2)\,a+d]\,\mu_{k+1+i}]\bigg{]}$
$\displaystyle=(i-1)\,\bigg{[}(k+2)\,\sigma_{k+i-1}-\sigma_{k+i-1}-(k+1)\,\sigma_{k+i-1}\bigg{]}=0,$

where we have used (5.4) and the fact that
$\sigma_{k+i-1}=a\,\mu_{k+1+i}+b\,\mu_{k+i}+c\,\mu_{k+i-1}$. This proves that
$\bar{x}_{k}-\bar{y}_{k}=\mathtt{0}$ and, in turn, that
$N_{k+2}^{\top}\,G_{k+1}^{(1)}\,N_{k+2}\,=\,-\mathtt{D}_{k+2}^{\top}\,G_{k+2}$.

It follows from the Principle of Induction that
$N_{n+1}^{\top}\,G_{n}^{(1)}\,N_{n+1}\,=\,-\mathtt{D}_{n+1}^{\top}\,G_{n+1}$
holds for $n\geqslant 0$. ∎

Now, let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence of real
numbers. For $n\geqslant 1$, let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the
Cholesky factorization of $G_{n}$ and let
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns of
$S_{n}^{\top}$. Consider the set
$\\{\bar{s}_{n,0}^{(1)},\ldots,\bar{s}_{n,n}^{(1)}\\}$ of vectors in
$\mathbb{R}^{n+1}$ with

$\bar{s}_{n,j}^{(1)}=\frac{1}{j+1}\,N_{n+1}\,\bar{s}_{n+1,j+1},\quad
0\leqslant j\leqslant n.$

Observe that $\\{\bar{s}_{n,0}^{(1)},\ldots,\bar{s}_{n,n}^{(1)}\\}$
constitutes a basis for $\mathbb{R}^{n+1}$. Furthermore, since $S_{n}^{\top}$
is a unit upper triangular matrix, the $(n+1)\times(n+1)$ matrix
$S_{n,1}^{\top}$ defined as

$S_{n,1}^{\top}=\left[\bar{s}_{n,0}^{(1)}\
\bar{s}_{n,1}^{(1)}\,\ldots\,\bar{s}_{n,n}^{(1)}\right],$ (5.5)

is also a unit upper triangular matrix. We show that $G_{n}^{(1)}$ defined in
(5.3) admits a Cholesky factorization with $S_{n,1}^{-1}$ its triangular
matrix factor.

######  Theorem 5.6.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence of real numbers
satisfying (3.1). For $n\geqslant 0$, let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$
be the Cholesky factorization of $G_{n}$. Then $G_{n}^{(1)}$ admits the
Cholesky factorization given by

$G_{n}^{(1)}=S_{n,1}^{-1}\,H_{n,1}\,S_{n,1}^{-\top},$

where $S_{n,1}^{\top}$ is the matrix defined in (5.5) and
$H_{n,1}=\textnormal{diag}[h_{0}^{(1)},\ldots,h_{n}^{(1)}]$ with

$h^{(1)}_{j}=-\frac{\lambda_{j+1}}{(j+1)^{2}}\,h_{j+1},\quad j\geqslant 0,$

and $\lambda_{j}=j\,[(j-1)\,a+d]$.

###### Proof.

For $0\leqslant j,k\leqslant n$, we compute

$\left(\bar{s}_{n,j}^{(1)}\right)^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,k}^{(1)}\,=\,\frac{1}{(j+1)\,(k+1)}\,\bar{s}_{n+1,j+1}^{\top}\,N_{n+1}^{\top}\,G_{n}^{(1)}\,N_{n+1}\,\bar{s}_{n+1,k+1}.$

Using Theorem 5.3 and Proposition 5.5, we obtain

$\displaystyle\left(\bar{s}_{n,j}^{(1)}\right)^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,k}^{(1)}\,$
$\displaystyle=\,-\frac{1}{(j+1)\,(k+1)}\,\bar{s}_{n+1,j+1}^{\top}\,\mathtt{D}_{n+1}^{\top}\,G_{n+1}\,\bar{s}_{n+1,k+1}$
$\displaystyle=-\frac{\lambda_{j+1}}{(j+1)\,(k+1)}\,\bar{s}_{n+1,j+1}^{\top}\,G_{n+1}\,\bar{s}_{n+1,k+1}$
$\displaystyle=-\frac{\lambda_{j+1}}{(j+1)\,(k+1)}\,\mathcal{B}_{n+1}(\bar{s}_{n+1,j+1},\,\bar{s}_{n+1,k+1})$
$\displaystyle=-\frac{\lambda_{j+1}}{(j+1)^{2}}\,h_{j+1}\,\delta_{j,k}.$

This implies that $S_{n,1}\,G_{n}^{(1)}\,S_{n,1}^{\top}=H_{n,1}$ and, hence,
$G_{n}^{(1)}=S_{n,1}^{-1}\,H_{n,1}\,S_{n,1}^{-\top}$. ∎

######  Corollary 5.7.

If $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence, then so is
$\\{\sigma_{n}\\}_{n\geqslant 0}$.

###### Proof.

By Theorem (5.4), $\\{\sigma_{n}\\}_{n\geqslant 0}$ is a pre-classical
sequence.

Now, we must show that $\det G_{n}^{(1)}\neq 0$ for $n\geqslant 0$ (Definition
3.3). From Theorem 5.6 we deduce that

$\det G_{n}^{(1)}=\det H_{n,1}=h^{(1)}_{0}\cdots h^{(1)}_{n},\quad n\geqslant
0,$

with

$h^{(1)}_{j}=-\frac{\lambda_{j+1}}{(j+1)^{2}}\,h_{j+1},\quad j\geqslant 0,$

and $\lambda_{j}=j\,[(j-1)\,a+d]$. Since $\\{\mu_{n}\\}_{n\geqslant 0}$ is
classical, then we have

$h_{n}\neq 0\quad\text{and}\quad n\,a+d\neq 0,\quad n\geqslant 0,$

(see equality (4.2) and Definition 3.1). This implies that $h^{(1)}_{j}\neq 0$
for $j\geqslant 0$. Therefore $\det G_{n}^{(1)}\neq 0$ for $n\geqslant 0$ and,
thus, $\\{\sigma_{n}\\}_{n\geqslant 0}$ is classical. ∎

Theorem 5.6 implies that for $n\geqslant 0$, the columns of $S_{n,1}^{\top}$
constitute an orthogonal basis for $\mathbb{R}^{n+1}$ with respect to the
bilinear form associated with $\\{\sigma_{n}\\}_{n\geqslant 0}$ (see
Definition 4.1), which we denote by

$\mathfrak{B}^{(1)}_{n}(\bar{u},\bar{v})\,=\,\bar{u}^{\top}\,G_{n}^{(1)}\,\bar{v},\quad\forall\bar{u},\bar{v}\in\mathbb{R}^{n+1}.$

We are ready for the following characterizations of classical sequences.

######  Theorem 5.8.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of real numbers such that
$\det G_{n}\neq 0$ for $n\geqslant 0$. Let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$
be the Cholesky factorization of $G_{n}$ and let
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns of
$S_{n}^{\top}$. Then $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical if and only
if there are real numbers $a,b,c$ satisfying

$|a|+|b|+|c|>0,$

such that the set $\\{\bar{s}_{n,0}^{(1)},\ldots,\bar{s}_{n,n}^{(1)}\\}$ of
vectors in $\mathbb{R}^{n+1}$ with

$\bar{s}_{n,j}^{(1)}=\frac{1}{j+1}\,N_{n+1}\,\bar{s}_{n+1,j+1},\quad
0\leqslant j\leqslant n,$

constitutes an orthogonal basis for $\mathbb{R}^{n+1}$ with respect to the
bilinear form associated with $\\{\sigma_{n}\\}_{n\geqslant 0}$, where

$\sigma_{n}\,=\,a\,\mu_{n+2}+b\,\mu_{n+1}+c\,\mu_{n},\qquad n\geqslant 0.$

###### Proof.

If $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical, then it follows from Theorem
5.6 that

$S_{n,1}\,G_{n}^{(1)}\,S_{n,1}^{\top}\,=\,H_{n,1},\quad n\geqslant 0,$

where

$S_{n,1}^{\top}=\left[\bar{s}_{n,0}^{(1)}\
\bar{s}_{n,1}^{(1)}\,\ldots\,\bar{s}_{n,n}^{(1)}\right],$

and $H_{n,1}=\textnormal{diag}[h_{0}^{(1)},\ldots,h_{n}^{(1)}]$ with

$h^{(1)}_{j}=-\frac{\lambda_{j+1}}{(j+1)^{2}}\,h_{j+1},\quad j\geqslant 0,$

and $\lambda_{j}=j\,[(j-1)\,a+d]$. This implies that

$\mathfrak{B}_{n}^{(1)}(\bar{s}_{n,i}^{(1)},\bar{s}_{n,j}^{(1)})=h^{(1)}_{j}\,\delta_{i,j},\quad
0\leqslant i,j\leqslant n.$

This proves the necessary condition.

Conversely, for $n\geqslant 0$, on one hand we have

$\mathfrak{B}_{n}^{(1)}(\bar{s}_{n,k}^{(1)},\bar{s}_{n,0}^{(1)})=h_{0}^{(1)}\,\delta_{k,0}$

or, equivalently,

$\left(\bar{s}_{n,k}^{(1)}\right)^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,0}^{(1)}\,=\,\frac{1}{k+1}\,\bar{s}_{n+1,k+1}^{\top}\,N_{n+1}^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,0}^{(1)}=-\lambda_{1}\,h_{1}\,\delta_{k,0}.$

Using the fact that $s_{n,0}^{(1)}=\bar{e}_{0}$ where, recall, $\bar{e}_{0}$
is the first column of the identity matrix of order $n+1$, we write

$\bar{s}_{n+1,k+1}^{\top}\,N_{n+1}^{\top}\,\begin{bmatrix}\sigma_{0}\\\
\sigma_{1}\\\ \vdots\\\
\sigma_{n}\end{bmatrix}=-(k+1)\,\lambda_{1}\,h_{1}\,\delta_{k,0}.$

On the other hand

$\mathfrak{B}_{n+1}(\bar{s}_{n+1,k+1},\bar{s}_{n+1,1})=\bar{s}_{n+1,k+1}^{\top}\,G_{n+1}\,\bar{s}_{n+1,1}=h_{1}\,\delta_{k,0}=h_{1}\,(k+1)\,\delta_{k,0}.$

Therefore,

$\bar{s}_{n+1,k+1}^{\top}\,N_{n+1}^{\top}\,\begin{bmatrix}\sigma_{0}\\\
\sigma_{1}\\\ \vdots\\\
\sigma_{n}\end{bmatrix}=-\lambda_{1}\,\bar{s}_{n+1,k+1}^{\top}\,G_{n+1}\,\bar{s}_{n+1,1},$

and, thus,

$\bar{s}_{n+1,k+1}^{\top}\left(N_{n+1}^{\top}\,\begin{bmatrix}\sigma_{0}\\\
\sigma_{1}\\\ \vdots\\\
\sigma_{n}\end{bmatrix}+\lambda_{1}\,G_{n+1}\,\bar{s}_{n+1,1}\right)=0,\quad
0\leqslant k\leqslant n.$

Since the first row of $N_{n+1}^{\top}$ is $[0\ 0\ldots 0]$ and
$\bar{s}_{n+1,0}^{\top}\,G_{n+1}\,\bar{s}_{n+1,1}=0$, we have

$S_{n+1}^{\top}\left(N_{n+1}^{\top}\,\begin{bmatrix}\sigma_{0}\\\
\sigma_{1}\\\ \vdots\\\
\sigma_{n}\end{bmatrix}+\lambda_{1}\,G_{n+1}\,\bar{s}_{n+1,1}\right)=\mathtt{0},$

and, consequently,

$N_{n+1}^{\top}\,\begin{bmatrix}\sigma_{0}\\\ \sigma_{1}\\\ \vdots\\\
\sigma_{n}\end{bmatrix}+\lambda_{1}\,G_{n+1}\,\bar{s}_{n+1,1}=\mathtt{0},$
(5.6)

where we have used the fact that $S_{n}^{\top}$ is an invertible matrix. If

$\bar{s}_{n+1,1}=[s_{1,0}\ 1\ 0\ldots 0]^{\top},$

then let $d=\lambda_{1}$ and $e=\lambda_{1}\,s_{1,0}$. Observe that by Theorem
4.3, $s_{1,0}$ is independent of $n$. In this way, the entries of (5.6) read

$\displaystyle d\,\mu_{1}+e\,\mu_{0}=0,$ $\displaystyle
k\,\sigma_{k-1}+d\,\mu_{k+1}+e\,\mu_{k}=(k\,a+d)\,\mu_{k+1}+(k\,b+e)\,\mu_{k}+k\,c\,\mu_{k-1}=0,\quad
1\leqslant k\leqslant n,$

which proves that $\\{\mu_{n}\\}_{n\geqslant 0}$ satisfies (3.1) and, thus,
$\\{\mu_{n}\\}_{n\geqslant 0}$ is classical. ∎

Corollary 5.7 shows that the classical character of $\\{\mu_{n}\\}_{n\geqslant
0}$ is inherited by $\\{\sigma_{n}\\}_{n\geqslant 0}$ allowing us to apply all
of our previous results about classical sequences to
$\\{\sigma_{n}\\}_{n\geqslant 0}$ and $\\{G_{n}^{(1)}\\}_{n\geqslant 0}$.
Therefore, we can define a new sequence $\\{\sigma_{n}^{(2)}\\}_{n\geqslant
0}$ as follows:

$\sigma_{n}^{(2)}=a\,\sigma_{n+2}+b\,\sigma_{n+1}+c\,\sigma_{n},\quad
n\geqslant 0.$

We denote by $\\{G_{n}^{(2)}\\}_{n\geqslant 0}$ the sequence of
$(n+1)\times(n+1)$ matrices with

$G_{0}^{(2)}=\sigma_{0}^{(2)},\quad\text{and}\quad
G_{n}^{(2)}=\left[\begin{array}[]{ccc|c}&&&\sigma_{n}^{(2)}\\\
&G_{n-1}^{(2)}&&\vdots\\\ &&&\sigma_{2n-1}^{(2)}\\\
\hline\cr\sigma_{n}^{(2)}&\ldots&\sigma_{2n-1}^{(2)}&\sigma_{2n}^{(2)}\end{array}\right],\quad
n\geqslant 1.$

By Theorem 5.4, $\\{\sigma_{n}^{(2)}\\}_{n\geqslant 0}$ is pre-classical
satisfying

$(n\,a+d_{2})\,\sigma_{n+1}^{(2)}+(n\,b+e_{2})\,\sigma_{n}^{(2)}+n\,c\,\sigma_{n-1}^{(2)}\,=\,0,\quad
n\geqslant 0,$

where $d_{2}=d+4\,a$ and $e_{2}=e+2\,b$. Moreover, if $\det G_{n}^{(2)}\neq 0$
for $n\geqslant 0$, then $\\{\sigma_{n}^{(2)}\\}_{n\geqslant 0}$ is classical
and, by Theorem 5.8, the set
$\\{\bar{s}_{n,0}^{(2)},\ldots,\bar{s}_{n,n}^{(2)}\\}$ of vectors in
$\mathbb{R}^{n+1}$ with

$\bar{s}_{n,j}^{(2)}=\frac{1}{j+1}\,N_{n+1}\,\bar{s}_{n+1,j+1}^{(1)},\quad
0\leqslant j\leqslant n,$

constitutes an orthogonal basis for $\mathbb{R}^{n+1}$ with respect to the
bilinear form $\mathfrak{B}_{n}^{(2)}$ associated with
$\\{\sigma_{n}^{(2)}\\}_{n\geqslant 0}$, that is,

$\mathfrak{B}^{(2)}_{n}(\bar{s}_{n,j}^{(2)},\bar{s}_{n,k}^{(2)})=h_{j}^{(2)}\,\delta_{j,k},\quad
0\leqslant j,k\leqslant n,$

where

$h^{(2)}_{j}=-\frac{\lambda_{j+1}^{(1)}}{(j+1)^{2}}\,h_{j+1}^{(2)},\quad
j\geqslant 0,$

and $\lambda_{j}^{(1)}=j\,[(j-1)\,a+d_{1}]$.

Iterating this idea, we obtain the following result.

######  Corollary 5.9.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence of real numbers.
Let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the Cholesky factorization of
$G_{n}$ and let $\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the
columns of $S_{n}^{\top}$. For each $k\geqslant 1$, define the sequence of
real numbers $\\{\sigma_{n}^{(k)}\\}_{n\geqslant 0}$ by

$\sigma_{n}^{(k)}\,=\,a\,\sigma_{n+2}^{(k-1)}+b\,\sigma_{n+1}^{(k-1)}+c\,\sigma_{n}^{(k-1)},\quad
n\geqslant 0,$

where $\sigma_{n}^{(0)}=\mu_{n}$ for $n\geqslant 0$. Then
$\\{\sigma_{n}^{(k)}\\}_{n\geqslant 0}$ is classical satisfying

$(n\,a+d_{k})\,\sigma_{n+1}^{(k)}+(n\,b+e_{k})\,\sigma_{n}^{(k)}+n\,c\,\sigma_{n-1}^{(k)}\,=\,0,\quad
n\geqslant 0,$

where $d_{k}=d+2\,k\,a$ and $e_{k}=e+k\,b$. Moreover, the set
$\\{\bar{s}_{n,0}^{(k)},\ldots,\bar{s}_{n,n}^{(k)}\\}$ of vectors in
$\mathbb{R}^{n+1}$ with

$\bar{s}_{n,j}^{(k)}=\frac{1}{j+1}\,N_{n+1}\,\bar{s}_{n+1,j+1}^{(k-1)},\quad
0\leqslant j\leqslant n,$

where $s_{n,j}^{(0)}=\bar{s}_{n,j}$, constitutes an orthogonal basis for
$\mathbb{R}^{n+1}$ with respect to the bilinear form associated with
$\\{\sigma_{n}^{(k)}\\}_{n\geqslant 0}$.

Observe that the vectors
$\\{\bar{s}_{n,0}^{(k)},\ldots,\bar{s}_{n,n}^{(k)}\\}$ in Corollary 5.9 can be
written in terms of the vectors
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ as follows: for each
$k\geqslant 1$,

$\bar{s}_{n,j}^{(k)}=\frac{1}{(j+1)_{k}}N_{n+1}\,N_{n+2}\cdots
N_{n+k}\,\bar{s}_{n+k,j+k},\quad 0\leqslant j\leqslant n,$ (5.7)

where $(\nu)_{k}=\nu\,(\nu+1)\cdots(\nu+k-1)$, $(\nu)_{0}=1$, denotes the
Pochhammer symbol. If $\\{G_{n}^{(k)}\\}_{n\geqslant 0}$ dentoes the sequence
of $(n+1)\times(n+1)$ matrices with

$G_{0}^{(k)}=\sigma_{0}^{(k)},\quad\text{and}\quad
G_{n}^{(k)}=\left[\begin{array}[]{ccc|c}&&&\sigma_{n}^{(k)}\\\
&G_{n-1}^{(k)}&&\vdots\\\ &&&\sigma_{2n-1}^{(k)}\\\
\hline\cr\sigma_{n}^{(k)}&\ldots&\sigma_{2n-1}^{(k)}&\sigma_{2n}^{(k)}\end{array}\right],\quad
n\geqslant 1,$ (5.8)

then the orthogonality of
$\\{\bar{s}_{n,0}^{(k)},\ldots,\bar{s}_{n,n}^{(k)}\\}$ with respect to the
bilinear form $\mathfrak{B}_{n}^{(k)}$ associated with
$\\{\sigma_{n}^{(k)}\\}_{n\geqslant 0}$ is given by

$\mathfrak{B}^{(k)}_{n}(\bar{s}_{n,j}^{(k)},\bar{s}_{n,i}^{(k)})\,=\,(\bar{s}_{n,j}^{(k)})^{\top}\,G_{n}^{(k)}\,\bar{s}_{n,i}^{(k)}\,=\,h_{j}^{(k)}\,\delta_{j,i},\quad
0\leqslant i,j\leqslant n,$ (5.9)

where

$h^{(k)}_{j}=-\frac{\lambda_{j+1}^{(k-1)}}{(j+1)^{2}}\,h_{j+1}^{(k-1)},\quad
j\geqslant 0,$

and $\lambda_{j}^{(k)}=j\,[(j-1)\,a+d_{k}]$. Note that we can write

$h^{(k)}_{j}=(-1)^{k}\frac{\prod_{i=0}^{k-1}\lambda_{j+k-i}^{(i)}}{[(j+1)_{k}]^{2}}\,h_{j+k},\quad
j\geqslant 0.$

Moreover, (5.9) implies the Cholesky factorization of $G_{n}^{(k)}$:

$S_{n,k}\,G_{n}^{(k)}\,S_{n,k}^{\top}\,=\,H_{n,k},\quad n\geqslant 0,$

where

$S_{n,k}^{\top}=\left[\bar{s}_{n,0}^{(k)}\
\bar{s}_{n,1}^{(k)}\,\ldots\,\bar{s}_{n,n}^{(k)}\right],$

and $H_{n,k}=\textnormal{diag}[h_{0}^{(k)},\ldots,h_{n}^{(k)}]$.

Let us reformulate the above results in terms of polynomials and moment
functionals. Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence of
real numbers satisfying (3.1), and let $\mathbf{u}$ be the moment functional
defined as $\mu_{n}=\langle\mathbf{u},x^{n}\rangle$, $n\geqslant 0$. For each
$k\geqslant 1$, the sequence of real numbers
$\\{\sigma_{n}^{(k)}\\}_{n\geqslant 0}$ defined by

$\sigma_{n}^{(k)}\,=\,a\,\sigma_{n+2}^{(k-1)}+b\,\sigma_{n+1}^{(k-1)}+c\,\sigma_{n}^{(k-1)},\quad
n\geqslant 0,$

where $\sigma_{n}^{(0)}=\mu_{n}$ for $n\geqslant 0$, is the sequence of
moments of the functional given by $\mathbf{v}_{k}\,=\,\phi^{k}\,\mathbf{u}$
where $\phi(x)=a\,x^{2}+b\,x+c$. Moreover, for $n\geqslant 0$, let
$S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be the Cholesky factorization of $G_{n}$
and let $\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns
of $S_{n}^{\top}$. For $n\geqslant 0$, Corollary 5.9 implies that the
polynomials $\\{Q_{0,k}(x),Q_{1,k}(x),\ldots,Q_{n,k}(x)\\}$ given by

$Q_{j,k}(x)\,=\,(\bar{s}_{n,j}^{(k)})^{\top}\,\mathtt{X}_{n},\quad 0\leqslant
j\leqslant n,$

satisfy
$\langle\mathbf{v}_{k},Q_{i,k}\,Q_{j,k}\rangle=h_{j}^{(k)}\,\delta_{i,j}$.
Therefore, $\\{Q_{0,k}(x),Q_{1,k}(x),\ldots,Q_{n,k}(x)\\}$ constitutes an
orthogonal basis for $\Pi_{n}$ with respect to $\mathbf{v}_{k}$. Furthermore,
Theorem 4.3 allows us to write

$Q_{n,k}(x)\,=\,(\bar{s}_{n,n}^{(k)})^{\top}\,\mathtt{X}_{n},\quad n\geqslant
0,$

and, in this way, $\\{Q_{n,k}(x)\\}_{n\geqslant 0}$ is an MOPS associated with
$\mathbf{v}_{k}$. We note that if

$P_{n}(x)\,=\,\bar{s}_{n,n}^{\top}\,\mathtt{X}_{n},\quad n\geqslant 0,$

then, by (5.7), we have

$Q_{n,k}(x)\,=\,\frac{P_{n+k}^{(k)}(x)}{(n+1)_{k}},\quad n\geqslant 0,$

where $P_{n}^{(k)}(x)$ denotes the $k$-th order derivative of $P_{n}(x)$.

###  5.3. First structure relation

For $n\geqslant 1$ and given real numbers $a,b$, and $c$, we define the
$(n+3)\times(n+1)$ matrices

$\Phi_{n}=\left[\begin{array}[]{c|c}&\mathtt{0}\\\ \Phi_{n-1}&c\\\ &b\\\
\hline\cr\mathtt{0}&a\end{array}\right],\quad\Phi_{0}=\begin{bmatrix}c\\\ b\\\
a\end{bmatrix}.$ (5.10)

######  Lemma 5.10.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ and $\\{\sigma_{n}\\}_{n\geqslant 0}$ be
sequences of real numbers satisfying

$\sigma_{n}=a\,\mu_{n+2}+b\,\mu_{n+1}+c\,\mu_{n},\quad n\geqslant 0,$

where $a,b,c\in\mathbb{R}$. Then, for $n\geqslant 2$,

$G_{n}^{(1)}\,\begin{bmatrix}\bar{v}\\\ 0\\\
0\end{bmatrix}\,=\,G_{n}\,\Phi_{n-2}\,\bar{v},\quad\forall\,\bar{v}\in\mathbb{R}^{n-1}.$
(5.11)

###### Proof.

We use induction to prove this result. For $n=2$, observe that, for all
$v\in\mathbb{R}$,

$G_{2}^{(1)}\,\begin{bmatrix}v\\\ 0\\\
0\end{bmatrix}=\begin{bmatrix}\sigma_{0}\\\ \sigma_{1}\\\
\sigma_{2}\end{bmatrix}v=a\,v\,\begin{bmatrix}\mu_{2}\\\ \mu_{3}\\\
\mu_{4}\end{bmatrix}+b\,v\,\begin{bmatrix}\mu_{1}\\\ \mu_{2}\\\
\mu_{3}\end{bmatrix}+c\,v\,\begin{bmatrix}\mu_{0}\\\ \mu_{1}\\\
\mu_{2}\end{bmatrix}=G_{2}\,\Phi_{0}\,v.$

This proves the base case.

Suppose that (5.11) holds for some $k\geqslant 2$. Let

$\bar{v}=\begin{bmatrix}\nu_{1}\\\ \vdots\\\
\nu_{k}\end{bmatrix}\in\mathbb{R}^{k}.$

We compute

$\displaystyle G_{k+1}^{(1)}\begin{bmatrix}\bar{v}\\\ 0\\\
0\end{bmatrix}=G_{k+1}^{(1)}\left(\begin{bmatrix}\nu_{1}\\\ \vdots\\\
\nu_{k-1}\\\ 0\\\ 0\\\ 0\end{bmatrix}+\begin{bmatrix}0\\\ \vdots\\\ 0\\\
\nu_{k}\\\ 0\\\ 0\end{bmatrix}\right).$

Multiplying by blocks and using the induction hypothesis, we get

$\displaystyle G_{k+1}^{(1)}\begin{bmatrix}\bar{v}\\\ 0\\\ 0\end{bmatrix}$
$\displaystyle=\left[\begin{array}[]{ccc|c}&&&\sigma_{k+1}\\\
&G_{k}^{(1)}&&\vdots\\\ &&&\sigma_{2k+1}\\\
\hline\cr\sigma_{k+1}&\ldots&\sigma_{2k+1}&\sigma_{2k+2}\end{array}\right]\begin{bmatrix}\nu_{1}\\\
\vdots\\\ \nu_{k-1}\\\ 0\\\ 0\\\
0\end{bmatrix}+\nu_{k}\begin{bmatrix}\sigma_{k-1}\\\ \vdots\\\
\sigma_{2k-1}\\\ \sigma_{2k}\end{bmatrix}$
$\displaystyle=\left[\begin{array}[]{ccc}&G_{k}\,\Phi_{k-2}&\\\\[-6.0pt] \\\
\hline\cr\sigma_{k+1}&\cdots&\sigma_{2k-1}\end{array}\right]\begin{bmatrix}\nu_{1}\\\
\vdots\\\ \nu_{k-1}\end{bmatrix}+\nu_{k}\begin{bmatrix}\sigma_{k-1}\\\
\vdots\\\ \sigma_{2k-1}\\\ \sigma_{2k}\end{bmatrix}$
$\displaystyle=\left[\begin{array}[]{ccc|c}&&&\sigma_{k-1}\\\
&G_{k}\,\Phi_{k-2}&&\vdots\\\ &&&\sigma_{2k-1}\\\
\hline\cr\sigma_{k-1}&\ldots&\sigma_{2k-1}&\sigma_{2k}\end{array}\right]\,\bar{v}.$

Then, it is straightforward to verify that

$\left[\begin{array}[]{ccc|c}&&&\sigma_{k-1}\\\ &G_{k}\,\Phi_{k-2}&&\vdots\\\
&&&\sigma_{2k-1}\\\
\hline\cr\sigma_{k-1}&\ldots&\sigma_{2k-1}&\sigma_{2k}\end{array}\right]=\left[\begin{array}[]{ccc|c}&&&\mu_{k+1}\\\
&G_{k}&&\vdots\\\ &&&\mu_{2k+1}\\\
\hline\cr\mu_{k+1}&\ldots&\mu_{2k+1}&\mu_{2k+2}\end{array}\right]\left[\begin{array}[]{c|c}&\mathtt{0}\\\
\Phi_{k-2}&c\\\ &b\\\ \hline\cr\mathtt{0}&a\end{array}\right].$

This proves that

$G_{k+1}^{(1)}\,\begin{bmatrix}\bar{v}\\\ 0\\\
0\end{bmatrix}\,=\,G_{k+1}\,\Phi_{k-1}\,\bar{v},$

and, thus, (5.11) holds for $n\geqslant 2$. ∎

We have the following characterization of classical sequences of real numbers.

######  Theorem 5.11.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of real numbers such that
$\det G_{n}\neq 0$ for $n\geqslant 0$. Let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$
be the Cholesky factorization of $G_{n}$ and let
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns of
$S_{n}^{\top}$. Then $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical if and only
if there are real numbers $a,b,c$ satisfying

$|a|+|b|+|c|>0,$

and real numbers $a_{j},b_{j},c_{j}$, $j\geqslant 0$, with $c_{j}\neq 0$, such
that the set $\\{\bar{s}_{n,0}^{(1)},\ldots,\bar{s}_{n,n}^{(1)}\\}$ of vectors
in $\mathbb{R}^{n+1}$ with

$\bar{s}_{n,j}^{(1)}=\frac{1}{j+1}\,N_{n+1}\,\bar{s}_{n+1,j+1},\quad
0\leqslant j\leqslant n,$

satisfy

$\Phi_{n-2}\,\bar{s}_{n-2,j}^{(1)}=a_{j}\,\bar{s}_{n,j+2}+b_{j}\,\bar{s}_{n,j+1}+c_{j}\,\bar{s}_{n,j},\quad
0\leqslant j\leqslant n-2,$ (5.12)

with $\Phi_{n}$ as defined in (5.10)

###### Proof.

We will use the fact that

$S_{n,1}^{\top}=\left[\bar{s}_{n,0}^{(1)}\
\bar{s}_{n,1}^{(1)}\,\ldots\,\bar{s}_{n,n}^{(1)}\right],$

is a unit upper triangular matrix.

Suppose that $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence satisfying
(3.1). Since $\Phi_{n-2}\,\bar{s}_{n-2,j}^{(1)}\in\mathbb{R}^{n+1}$ and
$\\{\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}\\}$ constitutes an
orthogonal basis for $\mathbb{R}^{n+1}$, we can write

$\Phi_{n-2}\,\bar{s}_{n-2,j}^{(1)}=\sum_{k=0}^{n}a_{j,k}\,\bar{s}_{n,k},$

with

$a_{j,k}\,h_{k}\,=\,\bar{s}_{n,k}^{\top}\,G_{n}\,\Phi_{n-2}\,\bar{s}_{n-2,j}^{(1)},\quad
0\leqslant k\leqslant n.$

The absence of a subindex indicating some dependence on $n$ is justified by
Theorem 4.3 and Theorem 5.6, which imply that the expression for
$\bar{s}_{n,k}$ and $\bar{s}_{n-2,j}^{(1)}$ are independent of $n$.

By (5.11) and since

$S_{n,1}^{\top}=\left[\begin{array}[]{c|c}S_{n-2,1}^{\top}&\ast\\\
\hline\cr\mathtt{0}&\begin{matrix}1&\ast\\\
0&1\end{matrix}\end{array}\right],$

we have

$a_{j,k}\,h_{k}\,=\,\bar{s}_{n,k}^{\top}\,G_{n}^{(1)}\,\begin{bmatrix}\bar{s}_{n-2,j}^{(1)}\\\
0\\\
0\end{bmatrix}=\bar{s}_{n,k}^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,j}^{(1)},\quad
0\leqslant k\leqslant n,\quad 0\leqslant j\leqslant n-2.$

On one hand, by Corollary 5.9, the set
$\\{\bar{s}_{n,0}^{(1)},\ldots,\bar{s}_{n,n}^{(1)}\\}$ is an orthogonal basis
for $\mathbb{R}^{n+1}$, and since
$\bar{s}_{n,k}=\bar{s}_{n,k}^{(1)}+\alpha_{k,k-1}\,\bar{s}^{(1)}_{n,k-1}+\cdots+\alpha_{k,0}\,\bar{s}_{n,0}^{(1)}$,
we get that $a_{j,k}=0$ for $0\leqslant k\leqslant j-1$, and

$a_{j,j}\,h_{j}=\bar{s}_{n,j}^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,j}^{(1)}=(\bar{s}_{n,j}^{(1)})^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,j}^{(1)}=h_{j}^{(1)}.$

Therefore,

$a_{j,j}=\frac{h_{j}^{(1)}}{h_{j}}=-\frac{\lambda_{j+1}}{(j+1)^{2}}\neq 0.$

On the other hand, since $S_{n-2,1}^{\top}$ is upper triangular, we have

$\Phi_{n-2}\,\bar{s}_{n-2,j}^{(1)}=\begin{bmatrix}\ast\\\ a\\\
\mathtt{0}\end{bmatrix}\in\mathbb{R}^{n+1},\quad 0\leqslant j\leqslant n-2,$

where the last $n-j-2$ entries are zero. This implies that $a_{j,k}=0$ for
$k\geqslant j+3$. Therefore, (5.12) holds with $a_{j}=a_{j,j+2}$,
$b_{j}=a_{j,j+1}$ and $c_{j}=a_{j,j}\neq 0$.

Conversely, suppose that (5.12) holds with the entries of $\Phi_{n-2}$ are
real numbers $a,b,c$ satisfying $|a|+|b|+|c|>0$. On one hand, multiplying both
sides of (5.12) by $\bar{s}_{n,0}^{\top}\,G_{n}$, we obtain for $0\leqslant
j\leqslant n-2$,

$\bar{s}_{n,0}^{\top}\,G_{n}\left(a_{j}\,\bar{s}_{n,j+2}+b_{j}\,\bar{s}_{n,j+1}+c_{j}\,\bar{s}_{n,j}\right)=\bar{s}_{n,0}^{\top}\,G_{n}\Phi_{n-2}\,\bar{s}_{n-2,j}^{(1)}.$

By the orthogonality of
$\\{\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}\\}$, we have

$\bar{s}_{n,0}^{\top}\,G_{n}\Phi_{n-2}\,N_{n-1}\,\bar{s}_{n-1,j+1}\,=\,(j+1)\,c_{0}\,h_{0}\,\delta_{j,0}.$

On the other hand, we have

$\bar{s}_{n-1,1}^{\top}\,G_{n-1}\,\bar{s}_{n-1,j+1}\,=\,h_{1}\,\delta_{j,0}=(j+1)\,h_{1}\,\delta_{j,0},\quad
0\leqslant j\leqslant n-2.$

Therefore,

$\bar{s}_{n,0}^{\top}\,G_{n}\Phi_{n-2}\,N_{n-1}\,\bar{s}_{n-1,j+1}\,=\,c_{0}\,\frac{h_{0}}{h_{1}}\,\bar{s}_{n-1,1}^{\top}\,G_{n-1}\,\bar{s}_{n-1,j+1},$

or, equivalently,

$\left(\bar{s}_{n,0}^{\top}\,G_{n}\,\Phi_{n-2}\,N_{n-1}-c_{0}\,\frac{h_{0}}{h_{1}}\,\bar{s}_{n-1,1}^{\top}\,G_{n-1}\right)\,\bar{s}_{n-1,j+1}=0,\quad
0\leqslant j\leqslant n-2.$

Since $N_{n-1}\bar{s}_{n-1,0}=\mathtt{0}$ and
$\bar{s}_{n-1,1}^{\top}\,G_{n-1}\,\bar{s}_{n-1,0}=0$, we have

$\left(\bar{s}_{n,0}^{\top}\,G_{n}\,\Phi_{n-2}\,N_{n-1}-c_{0}\,\frac{h_{0}}{h_{1}}\,\bar{s}_{n-1,1}^{\top}\,G_{n-1}\right)\,S_{n-1}^{\top}=\mathtt{0},$

and, therefore,

$\bar{s}_{n,0}^{\top}\,G_{n}\,\Phi_{n-2}\,N_{n-1}-c_{0}\,\frac{h_{0}}{h_{1}}\,\bar{s}_{n-1,1}^{\top}\,G_{n-1}=\mathtt{0},$
(5.13)

where we have used the fact that $S_{n-1}^{\top}$ is an invertible matrix. If

$\bar{s}_{n-1,1}^{\top}=[s_{1,0}\ 1\ 0\ldots 0],$

then let $d=-c_{0}\frac{h_{0}}{h_{1}}$ and $e=d\,s_{1,0}$. In this way, the
entries of (5.13) read

$\displaystyle d\,\mu_{1}+e\,\mu_{0}=0,$
$\displaystyle(k\,a+d)\,\mu_{k+1}+(k\,b+e)\,\mu_{k}+k\,c\,\mu_{k-1}=0,\quad
1\leqslant k\leqslant n-1,$

which proves that $\\{\mu_{n}\\}_{n\geqslant 0}$ satisfies (3.1) and, thus,
$\\{\mu_{n}\\}_{n\geqslant 0}$ is classical. ∎

We briefly recast Theorem 5.11 in terms of polynomials. Let
$\\{\mu_{n}\\}_{n\geqslant 0}$ be a classical sequence of real numbers, and
let $\mathbf{u}$ be the moment functional defined as
$\mu_{n}=\langle\mathbf{u},x^{n}\rangle$, $n\geqslant 0$. Let
$\\{P_{n}(x)\\}_{n\geqslant 0}$ and $\\{Q_{n}(x)\\}_{n\geqslant 0}$ let be
sequences of polynomials with

$P_{n}(x)\,=\,\bar{s}_{n,n}^{\top}\,\mathtt{X}_{n}\quad\text{and}\quad
Q_{n}(x)\,=\,(\bar{s}_{n,n}^{(1)})^{\top}\,\mathtt{X}_{n}\quad n\geqslant 0.$

Then $\\{P_{n}(x)\\}_{n\geqslant 0}$ is an OPS associated with $\mathbf{u}$,
and

$Q_{j}(x)\,=\,\frac{P_{j+1}^{\prime}(x)}{j+1},\quad 0\leqslant j\leqslant n.$

For $a,b,c\in\mathbb{R}$ such that $|a|+|b|+|c|>0$, and $n\geqslant 0$, the
matrix $\Phi_{n}$ is the matrix representation of the linear mapping from
$\Pi_{n}$ to $\Pi_{n+2}$ defined by $p(x)\mapsto\phi(x)\,p(x)$ with
$\phi(x)=a\,x^{2}+b\,x+c$. Therefore,

$(\Phi_{n}\,\bar{s}_{n,n}^{(1)})^{\top}\,\mathtt{X}_{n+2}=\phi(x)\,Q_{n}(x).$

In this way, by Theorem 5.11, $\mathbf{u}$ is a classical moment functional if
and only if there is a non zero polynomial $\phi(x)$ with
$\deg\phi(x)\leqslant 2$, and real numbers $a_{n},b_{n},c_{n}$, $n\geqslant
0$, with $c_{n}\neq 0$,

$\phi(x)\,Q_{n}(x)=a_{n}\,p_{n+2}(x)+b_{n}\,p_{n+2}(x)+c_{n}\,p_{n}(x),\quad
n\geqslant 0.$

###  5.4. Second structure relation

The following characterization is similar to (5.12) but it has a dual flavor
in the sense that the roles of $\\{\bar{s}_{n,j}\\}_{j=0}^{n}$ and
$\\{\bar{s}_{n,j}^{(1)}\\}_{j=0}^{n}$ are interchanged.

######  Theorem 5.12.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of real numbers such that
$\det G_{n}\neq 0$ for $n\geqslant 0$. Let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$
be the Cholesky factorization of $G_{n}$ and let
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns of
$S_{n}^{\top}$. Then $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical if and only
if there are real numbers $\kappa_{j},\xi_{j}$, $j\geqslant 0$, such that

$\bar{s}_{n,j}\,=\,\bar{s}_{n,j}^{(1)}+\kappa_{j}\,\bar{s}_{n,j-1}^{(1)}+\xi_{j}\,\bar{s}_{n,j-2}^{(1)},\quad
0\leqslant j\leqslant n,$ (5.14)

where, by convention, we set
$\bar{s}_{n,-2}^{(1)}\,=\,\bar{s}_{n,-1}^{(1)}\,=\,\mathtt{0}$.

###### Proof.

We will use the fact that

$S_{n,1}^{\top}=\left[\bar{s}_{n,0}^{(1)}\
\bar{s}_{n,1}^{(1)}\,\ldots\,\bar{s}_{n,n}^{(1)}\right],$

is a unit upper triangular matrix.

Suppose that $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence. By
Corollary 5.9, $\\{\bar{s}_{n,0}^{(1)},\ldots,\bar{s}_{n,n}^{(1)}\\}$ is an
orthogonal basis for $\mathbb{R}^{n+1}$. Since $S_{n,1}^{\top}$ is a unit
upper triangular matrix, we can write

$\bar{s}_{n,j}\,=\,\bar{s}_{n,j}^{(1)}+\sum_{k=0}^{j-1}a_{j,k}\,\bar{s}_{n,k}^{(1)},\quad
0\leqslant j\leqslant n,$

where

$a_{j,k}\,h_{k}^{(1)}\,=\,(\bar{s}_{n,k}^{(1)})^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,j},\quad
0\leqslant k\leqslant j-1.$

As before, the absence of a subindex indicating some dependence on $n$ is
justified by Theorem 4.3 and Theorem 5.6, which imply that the expression for
$\bar{s}_{n,k}$ and $\bar{s}_{n-2,j}^{(1)}$ are independent of $n$. Using
(5.11) and Theorem 5.11, we get

$a_{j,k}\,h_{k}^{(1)}\,=\,(a_{k}\,\bar{s}_{n,k+2}+b_{k}\,\bar{s}_{n,k+1}+c_{k}\,\bar{s}_{n,k})^{\top}\,G_{n}\,\bar{s}_{n,j},\quad
0\leqslant k\leqslant n-2.$

From the orthogonality of $\\{\bar{s}_{n,0},\ldots,\bar{s}_{n,n}\\}$, we
obtain $a_{j,k}\,=\,0$ for $k\leqslant j-3$. Therefore, (5.14) holds with
$\kappa_{j}=a_{j,j-1}$ and $\xi_{j}=a_{j,j-2}$.

Conversely, suppose that there are real numbers $\kappa_{j},\xi_{j}$,
$j\geqslant 0$, such that (5.14) holds. For $n\geqslant 0$, define the
$(n+1)\times(n+1)$ unit upper triangular matrices

$U_{0}=1,\quad U_{1}=\begin{bmatrix}1&\kappa_{1}\\\ 0&1\end{bmatrix},\quad
U_{n}=\left[\begin{array}[]{c|c}&\mathtt{0}\\\ U_{n-1}&\xi_{n}\\\
&\kappa_{n}\\\ \hline\cr\mathtt{0}&1\end{array}\right],\quad n\geqslant 2.$

Then (5.14) can be written as

$S_{n}^{\top}\,=\,S_{n,1}^{\top}\,U_{n}.$

Using $H_{n}=S_{n}\,G_{n}\,S_{n}^{\top}$, we get

$U_{n}\,=\,U_{n}\,H_{n}^{-1}\,S_{n}\,G_{n}\,S_{n}^{\top}\,=\,U_{n}\,H_{n}^{-1}\,S_{n}\,G_{n}\,S_{n,1}^{\top}\,U_{n},$

which implies

$I_{n+1}\,=\,U_{n}\,H_{n}^{-1}\,S_{n}\,G_{n}\,S_{n,1}^{\top}.$

From this equality we get

$S_{n,1}^{-\top}\,=\,U_{n}\,H_{n}^{-1}\,S_{n}\,G_{n}.$ (5.15)

On one hand, observe that

$\displaystyle N_{n+1}\,S_{n+1}^{\top}$
$\displaystyle=\,\begin{bmatrix}\mathtt{0}&N_{n+1}\,\bar{s}_{n+1,1}&N_{n+1}\,\bar{s}_{n+1,2}&\ldots&N_{n+1}\,\bar{s}_{n+1,n+1}\end{bmatrix}$
(5.16)
$\displaystyle=\,\begin{bmatrix}\mathtt{0}&\bar{s}_{n,0}^{(1)}&2\,\bar{s}_{n,1}^{(1)}&\ldots&(n+1)\,\bar{s}_{n,n}^{(1)}\end{bmatrix}$
$\displaystyle=\,S_{n,1}^{\top}\,N_{n+1},$

Combining (5.15) and (5.16), we get

$N_{n+1}\,=\,U_{n}\,H_{n}^{-1}\,S_{n}\,G_{n}\,N_{n+1}\,S_{n+1}^{\top}.$

On the other hand, we have

$\frac{1}{h_{1}}\,\bar{s}_{n+1,1}^{\top}\,G_{n+1}\,S_{n+1}^{\top}\,=\,\bar{s}_{n,0}^{\top}\,N_{n+1}.$

Therefore,

$\frac{1}{h_{1}}\,\bar{s}_{n+1,1}^{\top}\,G_{n+1}\,=\,\bar{s}_{n,0}^{\top}\,U_{n}\,H_{n}^{-1}\,S_{n}\,G_{n}\,N_{n+1}.$
(5.17)

If we let

$\begin{bmatrix}c\\\ b\\\ a\\\
\mathtt{0}\end{bmatrix}=\frac{\xi_{2}}{h_{2}}\,\bar{s}_{n,2}+\frac{\kappa_{1}}{h_{1}}\,\bar{s}_{n,1}+\frac{1}{h_{0}}\,\bar{s}_{n,0},\quad\text{and}\quad\begin{bmatrix}e\\\
d\\\ \mathtt{0}\end{bmatrix}=-\frac{1}{h_{1}}\,\bar{s}_{n+1,1},$ (5.18)

then the entries of (5.17) read

$\displaystyle-e\,\mu_{0}-d\,\mu_{1}$ $\displaystyle=0,$
$\displaystyle-e\,\mu_{k}-d\,\mu_{k+1}$
$\displaystyle=k\,c\,\mu_{k-1}+k\,b\,\mu_{k}+k\,a\,\mu_{k+1},\quad 1\leqslant
k\leqslant n+1,$

which proves that $\\{\mu_{n}\\}_{n\geqslant 0}$ satisfies (3.1) and, thus,
$\\{\mu_{n}\\}_{n\geqslant 0}$ is classical. ∎

For a classical sequence $\\{\mu_{n}\\}_{n\geqslant 0}$, let $\mathbf{u}$ be
the moment functional defined as $\mu_{n}=\langle\mathbf{u},x^{n}\rangle$,
$n\geqslant 0$, and let $\\{P_{n}(x)\\}_{n\geqslant 0}$ and
$\\{Q_{n}(x)\\}_{n\geqslant 0}$ let be sequences of polynomials with

$P_{n}(x)\,=\,\bar{s}_{n,n}^{\top}\,\mathtt{X}_{n}\quad\text{and}\quad
Q_{n}(x)\,=\,(\bar{s}_{n,n}^{(1)})^{\top}\,\mathtt{X}_{n}\quad n\geqslant 0.$

Theorem 5.12 implies that there are real numbers $\kappa_{n},\xi_{n}$,
$n\geqslant 0$, such that

$P_{n}(x)\,=\,Q_{n}(x)+\kappa_{n}\,Q_{n-1}(x)+\xi_{n}\,Q_{n-2}(x),\quad
n\geqslant 0,$

where, by convention, $Q_{-2}(x)=Q_{-1}(x)=0$. Moreover, we deduce from (5.18)
and Theorem 3.2 that $\mathbf{u}$ satisfies
$D(\phi\,\mathbf{u})=\psi\,\mathbf{u}$ with

$\phi(x)=\frac{\xi_{2}}{h_{2}}P_{2}(x)+\frac{\kappa_{1}}{h_{1}}P_{1}(x)+\frac{1}{h_{0}}P_{0}(x)\quad\text{and}\quad\psi(x)=-\frac{1}{h_{1}}P_{1}(x).$

###  5.5. Rodrigues-type formula

Recall that classical sequences of real numbers $\\{\mu_{n}\\}_{n\geqslant 0}$
satsisfy the three-term recurrence relation (3.1) which can be written in
matrix form as

$N_{n+1}^{\top}\,G_{n}^{(1)}\bar{s}_{n,0}+G_{n}\,\begin{bmatrix}e\\\ d\\\
\mathtt{0}\end{bmatrix}\,=\,\mathtt{0},\quad n\geqslant 1,$

with $G_{n}^{(1)}$ as defined in (5.3). The following characterization shows
that classical sequences satisfy higher order recurrence relations which can
be written in matrix form as well.

######  Theorem 5.13.

Let $\\{\mu_{n}\\}_{n\geqslant 0}$ be a sequence of real numbers such that
$\det G_{n}\neq 0$, $n\geqslant 0$. Let $S_{n}^{-1}\,H_{n}\,S_{n}^{-\top}$ be
the Cholesky factorization of $G_{n}$ and let
$\bar{s}_{n,0},\bar{s}_{n,1},\ldots,\bar{s}_{n,n}$ denote the columns of
$S_{n}^{\top}$. Then $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical if and only
if there are $a,b,c\in\mathbb{R}$ such that $|a|+|b|+|c|>0$, and non zero real
numbers $\varpi_{k}$, $k\geqslant 1$, such that

$N_{n+k}^{\top}\cdots
N_{n+1}^{\top}\,G_{n}^{(k)}\,\bar{s}_{n,0}\,=\,\varpi_{k}\,\,G_{n+k}\,\bar{s}_{n+k,k},\quad
n\geqslant 1,\quad 1\leqslant k\leqslant n,$ (5.19)

with $G_{n}^{(k)}$ as defined in (5.8).

###### Proof.

Suppose that $\\{\mu_{n}\\}_{n\geqslant 0}$ is classical satisfying (3.1). For
$k\geqslant 0$, by Corollary 5.9, the set
$\\{\bar{s}_{n,0}^{(k)},\ldots,\bar{s}_{n,n}^{(k)}\\}$ of vectors in
$\mathbb{R}^{n+1}$ with

$\bar{s}_{n,j}^{(k)}=\frac{1}{j+1}\,N_{n+1}\,\bar{s}_{n+1,j+1}^{(k-1)},\quad
0\leqslant j\leqslant n,$

where $s_{n,j}^{(0)}=\bar{s}_{n,j}$, constitutes an orthogonal basis for
$\mathbb{R}^{n+1}$ with respect to the bilinear form associated with
$\\{\sigma_{n}^{(k)}\\}_{n\geqslant 0}$. Let

$S_{n,k}^{\top}=\left[\bar{s}_{n,0}^{(k)}\
\bar{s}_{n,1}^{(k)}\,\ldots\,\bar{s}_{n,n}^{(k)}\right].$

On one hand, observe that for $k\geqslant 1$,

$\displaystyle N_{n+1}\,S_{n+1,k-1}^{\top}$
$\displaystyle=\,\begin{bmatrix}\mathtt{0}&N_{n+1}\,\bar{s}_{n+1,1}^{(k-1)}&N_{n+1}\,\bar{s}_{n+1,2}^{(k-1)}&\ldots&N_{n+1}\,\bar{s}_{n+1,n+1}^{(k-1)}\end{bmatrix}$
$\displaystyle=\,\begin{bmatrix}\mathtt{0}&\bar{s}_{n,0}^{(k)}&2\,\bar{s}_{n,1}^{(k)}&\ldots&(n+1)\,\bar{s}_{n,n}^{(k)}\end{bmatrix}$
$\displaystyle=\,S_{n,k}^{\top}\,N_{n+1},$

Then, since $\bar{s}^{(k)}_{n,0}=\bar{s}_{n,0}$ for $k\geqslant 1$,

$\displaystyle S_{n+k}\,N_{n+k}^{\top}\cdots
N_{n+1}^{\top}\,G_{n}^{(k)}\,\bar{s}_{n,0}$
$\displaystyle=\,N_{n+k}^{\top}\cdots
N_{n+1}^{\top}\,S_{n,k}\,G_{n}^{(k)}\,\bar{s}_{n,0}^{(k)}$
$\displaystyle=\,k!\,h_{0}^{(k)}\,\bar{e}_{k}\in\mathbb{R}^{n+k+1},$

where $\bar{e}_{k}$ denotes the $k+1$-th column of the identity matrix
$I_{n+k+1}$. On the other hand,

$\frac{1}{h_{k}}\,S_{n+k}\,G_{n+k}\,\bar{s}_{n+k,k}\,=\,\bar{e}_{k}.$

Since $S_{n+k}$ is invertible, it follows that

$N_{n+k}^{\top}\cdots
N_{n+1}^{\top}\,G_{n}^{(k)}\,\bar{s}_{n,0}\,=\,k!\,\frac{h_{0}^{(k)}}{h_{k}}\,\,G_{n+k}\,\bar{s}_{n+k,k},\quad
k\geqslant 1.$

Hence, (5.19) holds with $\varpi_{k}=k!\,\frac{h_{0}^{(k)}}{h_{k}}$.

Conversely, suppose that there are $a,b,c\in\mathbb{R}$ such that
$|a|+|b|+|c|>0$, and non zero real numbers $\varpi_{k}$, $k\geqslant 1$, such
that (5.19) holds. In particular, for $k=1$, we have

$N_{n+1}^{\top}\,G_{n}^{(1)}\,\bar{s}_{n,0}\,=\,\varpi_{1}\,\,G_{n+1}\,\bar{s}_{n+1,1},\quad
n\geqslant 0.$ (5.20)

If we let

$\begin{bmatrix}e\\\ d\end{bmatrix}=-\varpi_{1}\,\bar{s}_{1,1},$

and since $\bar{s}_{n+1,1}=\begin{bmatrix}\bar{s}_{1,1}\\\
\mathtt{0}\end{bmatrix}$, then the entries of (5.20) read

$\displaystyle-e\,\mu_{0}-d\,\mu_{1}$ $\displaystyle=0,$
$\displaystyle-e\,\mu_{k}-d\,\mu_{k+1}$
$\displaystyle=k\,c\,\mu_{k-1}+k\,b\,\mu_{k}+k\,a\,\mu_{k+1},\quad 1\leqslant
k\leqslant n+2,$

which proves that $\\{\mu_{n}\\}_{n\geqslant 0}$ satisfies (3.1) and, thus,
$\\{\mu_{n}\\}_{n\geqslant 0}$ is classical. ∎

We remark that if $\\{\mu_{n}\\}_{n\geqslant 0}$ is a classical sequence of
real numbers, and $\mathbf{u}$ be the moment functional defined as
$\mu_{n}=\langle\mathbf{u},x^{n}\rangle$, $n\geqslant 0$, then the entries of
(5.20) can be written as

$\langle
D(\phi\,\mathbf{u}),x^{k}\rangle\,=\,\langle\psi\,\mathbf{u},x^{k}\rangle,\quad
0\leqslant k\leqslant n+2,$

where $\phi(x)=a\,x^{2}+b\,x+c$ and $\psi(x)=d\,x+e$, which holds for all
$n\geqslant 0$. Hence, $\mathbf{u}$ satisfies (2.1). Moreover, it is
straightforward, but tedious, to verify that the entries of (5.19) can be
written as

$\langle
D^{k}(\phi^{k}\,\mathbf{u}),x^{j}\rangle\,=\,\langle\varpi_{k}\,P_{k}\,\mathbf{u},x^{j}\rangle,\quad
0\leqslant j\leqslant n+k+1,$

where $P_{k}(x)=\bar{s}_{n+k,k}^{\top}\,\mathtt{X}_{n+k}$, $k\geqslant 1$,
which holds for $n\geqslant 1$. Hence, $\mathbf{u}$ satisfies
$D^{k}(\phi^{k}\,\mathbf{u})=\varpi_{k}\,P_{k}\,\mathbf{u}$ for $k\geqslant 1$
(and holds for $k=0$ with $\varpi_{0}=1$).

## References

  * [1] R. Álvarez-Nodarse, On characterizations of classical polynomials, J. Comput. Appl. Math. 196 (2006), 320–337.
  * [2] W. Al Salam, T. S. Chihara, Another characterization of the classical orthogonal polynomials, SIAM J. Math. Anal. 3(1), (1972), 65-70.
  * [3] S. Bochner, Über Sturm-Liouvillesche Polynomsysteme, Math. Z. 29 (1929), no. 1, 730–736.
  * [4] F. A. Costabile, M. I. Gualtieri, A. Napoli, Matrix Calculus-Based Approach to Orthogonal Polynomial Sequences, Mediterr. J. Math. 17 (118) (2020).
  * [5] D. Dominici, Matrix factorizations and orthogonal polynomials, Random Matrices Theor. 9 (2020), 2040003, 33pp.
  * [6] J. Geronimus, On polynomials orthogonal with regard to a given sequence of numbers, Comm. Inst. Math. Mec. Univ. Kharkoff [Zapiski Inst. Math. Mech.] 17 (4) (1940), 3–18 (Russian).
  * [7] J. C. Garcia-Ardila, F. Marcellán, M. E. Marriaga, Orthogonal Polynomials and Linear Functionals. An Algebraic Approach and Applications, European Mathematical Society (EMS), Zürich, 2021.
  * [8] J. C. Garcia-Ardila, M. E. Marriaga, On Sobolev bilinear forms and polynomial solutions of second-order differential equations, Rev. Real Acad. Cienc. Exactas Fis. Nat. Ser. A-Mat. (2021).
  * [9] G. H. Golub, C. F. Van Loan, Matrix computations, Johns Hopkins Studies in Mathematical Sciences, Johns Hopkins University Press, Baltimore, MD, 4th ed., 2013.
  * [10] W. Hahn, Über die Jacobischen Polynome und Zwei Verwandte Polynomklassen, Math. Z. 39 (1935), 634–638.
  * [11] W. Hahn. Über Differentialgleichunqen für Orthoqonalpolynome. Monat. Math. 95, (1983), p. 269–274.
  * [12] R. A. Horn, C. R. Johnson, Matrix Analysis, Second ed., Cambridge University Presss, Cambridge, 2013.
  * [13] H. L. Krall, On derivatives of orthogonal polynomials II, Bull. Amer. Math. Soc. 47 (1941), 261–264.
  * [14] H. L. Krall, On higher derivatives orthogonal polynomials, Bull. Amer. Math. Soc. 42 (1936), 423–428.
  * [15] M. Mañas, Revisiting Biorthogonal Polynomials: An LU Factorization Discussion, In: Marcellán, F., Huertas, E.J. (eds) Orthogonal Polynomials: Current Trends and Applications. SEMA SIMAI Springer Series, vol 22. Springer, Cham.
  * [16] F. Marcellán, A. Branquinho, J. Petronilho, Classical orthogonal polynomials: a functional approach, Acta Appl. Math. 34 (1994), 283-303.
  * [17] F. Marcellán, J. Petronilho, On the solution of some distributional differential equations: existence and characterizations of the classical moment functionals, Integral Transform. Spec. Funct. 2 (3) (1994), 185–218.
  * [18] P. Maroni, Prolégomènes à l’étude des polynômes semiclassiques, Ann. Mat. Pura Appl. (4) 149 (1987), 165–184.
  * [19] M. Piñar, V. Ramírez, Matrix interpretation of formal orthogonal polynomials for non-definite functionals, J. Comput. Appl. Math. 18 (1987) 265–277.
  * [20] G. Szegő, Orthogonal polynomials, 4th edition, vol. 23. Amer. Math. Soc. Colloq. Publ., Amer. Math. Soc., Providence RI, 1975.
  * [21] F. Tricomi, Vorlesungen über orthogonalreihen, 2nd edn, Springer-Verlag, Berlin, 1970.
  * [22] L. Verde-Star, Characterization and construction of classical orthogonal polynomials using a matrix approach, Linear Algebra Appl. 438 (2013), 3635–3648.
  * [23] L. Verde-Star, Infinite Matrices in the Theory of Orthogonal Polynomials, In: Marcellán, F., Huertas, E.J. (eds) Orthogonal Polynomials: Current Trends and Applications. SEMA SIMAI Springer Series, vol 22. Springer, Cham.