Summary: Alternative expressions for calculating the prolate spheroidal radial functions of the second kind and their first derivatives with respect to are shown to provide accurate values over wide parameter ranges where the traditional expressions fail to do so.
The first alternative expression is obtained from the expansion of the product of and the prolate spheroidal angular function of the first kind in a series of products of the corresponding spherical functions. A similar expression for the radial functions of the first kind was shown previously to provide accurate values for the prolate spheroidal radial functions of the first kind and their first derivatives over all parameter ranges.
The second alternative expression for involves an integral of the product of and a spherical Neumann function kernel. It provides accurate values when is near unity and is not too large, even when becomes large and traditional expressions fail. The improvement in accuracy using the alternative expressions is quantified and discussed.