This routine uses Steedâs method to compute the values of the regular spherical Bessel functions j_l(x) for l from 0 to lmax inclusive for lmax >= 0 and x >= 0. The Steed/Barnett algorithm is described in Comp. Phys. Comm. 21, 297 (1981). Steedâs method is more stable than the recurrence used in the other functions but is also slower.
This function is from the GNU Scientific Library, see http://www.gnu.org/software/gsl/ for documentation.
Package: gsl