It has been known for a long time that there is a close connection
between stochastic processes and orthogonal polynomials. For example, N.
Wiener [112] and K. Ito [56] knew that Hermite polynomials play an
important role in the integration theory with respect to Brownian
motion. In the 1950s D. G. Kendall [66], W. Ledermann and G. E. H.
Reuter [67] [74], and S. Kar- lin and J. L. McGregor [59]
established another important connection. They expressed the transition
probabilities of a birth and death process by means of a spectral
representation, the so-called Karlin-McGregor representation, in terms
of orthogonal polynomials. In the following years these relation- ships
were developed further. Many birth and death models were related to
specific orthogonal polynomials. H. Ogura [87], in 1972, and D. D. En-
gel [45], in 1982, found an integral relation between the Poisson
process and the Charlier polynomials. Some people clearly felt the
potential im- portance of orthogonal polynomials in probability theory.
For example, P. Diaconis and S. Zabell [29] related Stein equations
for some well-known distributions, including Pearson's class, with the
corresponding orthogonal polynomials. The most important orthogonal
polynomials are brought together in the so-called Askey scheme of
orthogonal polynomials. This scheme classifies the hypergeometric
orthogonal polynomials that satisfy some type of differ- ential or
difference equation and stresses the limit relations between them.