toms726

toms726, a Fortran77 code which computes recursion coefficients for orthogonal polynomials, and the abscissa and weights for related quadrature rules, by Walter Gautschi. This is commonly called orthpol(). This is a version of ACM TOMS algorithm 726.

Licensing:

The computer code and data files made available on this web page are distributed under the MIT license

Languages:

toms726 is available in a FORTRAN77 version and a FORTRAN90 version.

Related Data and Programs:

toms726_test

toms655, a FORTRAN77 code which computes the weights for interpolatory quadrature rules.

toms793, a FORTRAN77 library which carries out Gauss quadrature for rational functions, by Walter Gautschi; this is ACM toms algorithm 793.

Reference:

1. William Cody, Kenneth Hillstrom,
   Chebyshev Approximations for the Natural Logarithm of the Gamma Function,
   Mathematics of Computation,
   Volume 21, Number 98, April 1967, pages 198-203.
2. Walter Gautschi,
   On Generating Orthogonal Polynomials,
   SIAM Journal on Scientific and Statistical Computing,
   Volume 3, Number 3, 1982, pages 289-317.
3. Walter Gautschi,
   Algorithm 726: ORTHPOL - A Package of Routines for Generating Orthogonal Polynomials and Gauss-Type Quadrature Rules,
   ACM Transactions on Mathematical Software,
   Volume 20, Number 1, March 1994, pages 21-62.

Source Code:

• toms726.f, the source code.
• toms726.sh, commands to compile the source code.