Computation of Special Functions

Abstract : Documentation is given for some subroutines which compute potentials and other functions. A set of subroutines uses rational approximations to compute Bessel functions of integral order. One subroutine uses the Debye approximation for the efficient computation of Bessel functions of complex argument and complex order. Empirical formulae have been developed to express the limiting boundaries of the modes of computation. (Author)