Algorithm 911: Multiple-Precision Exponential Integral and Related Functions
This article describes a collection of Fortran-95 routines for evaluating the exponential integral function, error function, sine and cosine integrals, Fresnel integrals, Bessel functions, and related mathematical special functions using the FM multiple-precision arithmetic package.
algorithms exponential integral function accuracy function evaluation floating point fortran mathematical library portable software multiple precision