Algorithm 861: Fortran 90 Subroutines for Computing the Expansion Coefficients of Mathieu Functions Using Blanch's Algorithm

摘要

A translation to Fortran 90 of Gertrude Blanch's algorithm for computing the expansion coefficients of the series that represent Mathieu functions is presented. Its advantages are portability, higher precision, practicality of use, and extended documentation. In addition, numerical validations and comparisons with other existing methods are presented.