Two Fortran 77 routines for the evaluation of Airy functions of complex arguments $Ai(z)$, $Bi(z)$ and their first derivatives are presented. The routines are based on the use of Gaussian quadrature, Maclaurin series and asymptotic expansions. Comparison with a previous code by D.E. Amos (ACM Trans. Math. Soft. 12 (1986)) is provided.

Bessel and Airy functions, cylinder functions, ${}_0F_1$ (msc 33C10), Computation of special functions, construction of tables (msc 65D20), Quadrature and cubature formulas (msc 65D32), Approximation in the complex domain (msc 30E10)
CWI
Modelling, Analysis and Simulation [MAS]
Computational Dynamics

Gil, A, Segura, J, & Temme, N.M. (2002). AIZ, BIZ: two Fortran 77 routines for the computation of complex Airy functions. Modelling, Analysis and Simulation [MAS]. CWI.