Abstract
Integral representations are considered of solutions of the Airy differential equation w ′′−zw=0 for computing Airy functions for complex values of z. In a first method contour integral representations of the Airy functions are written as non-oscillating integrals for obtaining stable representations, which are evaluated by the trapezoidal rule. In a second method an integral representation is evaluated by using generalized Gauss–Laguerre quadrature; this approach provides a fast method for computing Airy functions to a predetermined accuracy. Comparisons are made with well-known algorithms of Amos, designed for computing Bessel functions of complex argument. Several discrepancies with Amos' code are detected, and it is pointed out for which regions of the complex plane Amos' code is less accurate than the quadrature algorithms. Hints are given in order to build reliable software for complex Airy functions.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
M. Abramowitz and I.A. Stegun (eds.), Handbook of Mathematical Functions, National Bureau of Standards Applied Mathematics Series No. 55 (US Government Printing Office, Washington, DC, 1964).
D.E. Amos, Algorithm 644: A portable package for Bessel functions of a complex argument and nonnegative order, ACM Trans. Math. Softw. 12 (1986) 265-273.
R.M. Corless, D.J. Jeffrey and H. Rasmussen, Numerical evaluation of Airy functions with complex arguments, J. Comput. Phys. 99 (1992) 106-114.
W. Gautschi, Computation of Bessel and Airy functions and of related Gaussian quadrature, BIT 42 (2002) 110-118.
A. Gil, J. Segura and N.M. Temme, On non-oscillating integrals for computing inhomogeneous Airy functions, Math. Comput. 70 (2001) 1183-1194.
A. Gil, J. Segura and N.M. Temme, AiZ: a Fortran program for the computation of complex Airy functions (submitted).
E.T. Goodwin, The evaluation of integrals of the form \(\smallint _{ - \infty }^\infty f(x){\text{e}}^{ - x^2 } {\text{d}}x\), Cambridge Phil. Soc. Proc. 45 (1949) 241-245.
D.W. Lozier and F.W.J. Olver, Numerical evaluation of special functions, in: AMS Proceedings of Symposia in Applied Mathematics, Vol. 48, ed. W. Gautschi (1994) pp. 79-125.
Y.L. Luke, The Special Functions and Their Approximations, Vol. I-II (Academic Press, New York, 1969).
F.W.J. Olver, Asymptotics and Special Functions (Academic Press, New York, 1974, reprinted by A.K. Peters, 1997).
Z. Schulten, D.G.M. Anderson and R.G. Gordon, An algorithm for the evaluation of complex Airy functions, J. Comput. Phys. 31 (1979) 60-75.
N.M. Temme, Special Functions: An Introduction to the Classical Functions of Mathematical Physics (Wiley, New York, 1996).
R. Wong, Asymptotic Approximations of Integrals (Academic Press, New York, 1989).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Gil, A., Segura, J. & Temme, N.M. Computing Complex Airy Functions by Numerical Quadrature. Numerical Algorithms 30, 11–23 (2002). https://doi.org/10.1023/A:1015636825525
Issue Date:
DOI: https://doi.org/10.1023/A:1015636825525