It would be useful, given the general importance of orthogonal polynomials in applications, to have some shortcuts to evaluate the most important ones in "by name" in the subdirectory poly, instead of defining them by a list of values that are interpolated. A proposal/implementation for the Hermite, Laguerre and Chebyshev polynomials is attached. I can implement more of these if this finds support. -- --- Richard J Mathar Tel (+31) (0) 71 527 8459 Sterrewacht Universiteit Leiden Fax (+31) (0) 71 527 5819 Postbus 9513 The Netherlands URL http://www.strw.leidenuniv.nl/~mathar office: Niels Bohrweg 2, 2333 CA Leiden, #456