legen/jacobl.f:   compute the Gauss-Lobatto points of Jacobi
                  polynomials

legen/leinit.f:   initialization for Legendre transforms

cheby/chtrfm1.f:  Subroutine for Chebyshev transform (fftpack is needed)

dfft:             contains the double-precision version of fftpac