Partial wave expansion is encountered in a wide range of problems in theoretical physics. Typically, spherical Bessel functions are encountered in scattering of nuclei, atoms, and molecules. A routine is described for computing spherical Bessel and Neumann functions jl( x) and nl( x), respectively, by backward and forward recurrence relations. The method is fast, very accurate, not code intensive, and easily adopted for a personal computer. These features make the code well suited for the physics student research environment.

This content is only available via PDF.
You do not currently have access to this content.