On the computation of Mathieu functions - Delft Numerical Analysis Group (J. Engrg. Math., 1973)
Variable-precision Algol program, maximum precision 15D