A combined symbolic and numerical algorithm for the computation of zeros of orthogonal polynomials and special functions

A Maple algorithm for the computation of the zeros of orthogonal polynomials (OPs) and special functions (SFs) in a given interval [x(1), x(2)] is presented. The program combines symbolic and numerical calculations and it is based on fixed point iterations. The program uses as inputs the analytic expressions for the coefficients of the three-term recurrence relation and a difference-differential relation satisfied by the set of ON or SFs. The performance of the method is illustrated with several examples: Hermite, Chebyshev, Legendre, Jacobi and Gegenbauer polynomials, Bessel, Coulomb and Conical functions. (C) 2003 Elsevier Science Ltd. All rights reserved.