Quarterly of Applied Mathematics

Quarterly of Applied Mathematics

Online ISSN 1552-4485; Print ISSN 0033-569X



Accurate calculation of the modified Mathieu functions of integer order

Authors: Arnie L. Van Buren and Jeffrey E. Boisvert
Journal: Quart. Appl. Math. 65 (2007), 1-23
MSC (2000): Primary 33E10
DOI: https://doi.org/10.1090/S0033-569X-07-01039-5
Published electronically: January 9, 2007
MathSciNet review: 2313146
Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: Several different expressions involving infinite series are available for calculating the radial (i.e., modified) Mathieu functions of integer order. Mathieu functions depend on three parameters: order, radial coordinate, and a size parameter that is chosen here to be either real or imaginary. Expressions traditionally used to calculate the radial functions result in inaccurate function values over some parameter ranges due to unavoidable subtraction errors that occur in the series evaluation. For many of the expressions the errors for small orders increase without bound as the size parameter increases. In the present paper, the subtraction error obtained using traditional expressions is explored with regard to parameter values. Included is a discussion of the Bessel function product series, which has an integer offset for the order of the Bessel functions that is traditionally chosen to be zero (or one). It is shown here that the use of larger offset values that tend to increase with increasing radial function order usually eliminates the subtraction errors. This paper identifies the expressions and evaluation procedures that provide accurate radial Mathieu function values. A brief discussion of the calculation of the angular functions of the first kind that appear in many of these expressions is included. The paper also gives a description of a Fortran computer program that provides accurate values of radial Mathieu functions together with the associated angular functions over extremely wide parameter ranges. This effort was guided by recent advancements in the calculation of prolate spheroidal functions. The Mathieu functions are a special case of spheroidal functions, resulting in a similarity of behavior in their evaluation.

References [Enhancements On Off] (What's this?)

Similar Articles

Retrieve articles in Quarterly of Applied Mathematics with MSC (2000): 33E10

Retrieve articles in all journals with MSC (2000): 33E10

Additional Information

Arnie L. Van Buren
Affiliation: Cary, North Carolina

Jeffrey E. Boisvert
Affiliation: NAVSEA Newport, Newport, Rhode Island

DOI: https://doi.org/10.1090/S0033-569X-07-01039-5
Received by editor(s): October 4, 2005
Published electronically: January 9, 2007
Article copyright: © Copyright 2007 Brown University

American Mathematical Society