Computing special functions by using quadrature rules
The usual tools for computing special functions are power series, asymptotic expansions, continued fractions, differential equations, recursions, and so on. Rather seldom are methods based on quadrature of integrals. Selecting suitable integral representations of special functions, using principles from asymptotic analysis, we develop reliable algorithms which are valid for large domains of real or complex parameters. Our present investigations include Airy functions, Bessel functions and parabolic cylinder functions. In the case of Airy functions we have improvements in both accuracy and speed for some parts of Amos's code for Bessel functions.