Airy functions are solutions of the differential equation
Two linearly independent solutions that are real for real values of are denoted by and . They have the integral representation
where we assume that is real. In this paper the authors are concerned with the numerical evaluation of and for complex values of by numerical quadrature. 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.