 |
This paper not only reviews the various methodologies for evaluating the angular and radial prolate and oblate spheroidal functions and their eigenvalues, but also presents an efficient algorithm which is developed with the software package Mathematica. Two algorithms are developed for computation of the eigenvalues lmn and coefficients d^mn/r. Important steps in programming are provided for estimating eigenvalues of the spheroidal harmonics with a complex argument c. Furthermore, the starting and ending points for searching for the eigenvalues by Newton's method are successfully obtained. As compared with the published data by Caldwell or Press etal (for real argument) and Oguchi (for complex argument), the spheroidal harmonics and their eigenvalues estimated using this algorithm are of a much higher accuracy. In particular, a lot of tabulated data for the intermediate coefficients d^mn/p|r, the prolate and oblate radial spheroidal functions of the second king, and their first-order derivatives, as obtained by Flammer are found to be inaccurate, although these tabulated data have been considered as exact referenced results for about half a century. The algorithm developed here for evaluating the spheroidal harmonics with the Mathematica program is also found to be simple, fast, and numerically efficient, and of a much better accuracy than the other results tabulated by Flammer and others, being able to produce results of 100 significant digits or more.
|
|
For the newest resources, visit Wolfram Repositories and Archives »
