### Airy Functions

Function: double gsl_sf_airy_Ai (double x, gsl_mode_t mode)
Function: int gsl_sf_airy_Ai_e (double x, gsl_mode_t mode, gsl_sf_result * result)
These routines compute the Airy function @math{Ai(x)} with an accuracy specified by mode.

Function: double gsl_sf_airy_Bi (double x, gsl_mode_t mode)
Function: int gsl_sf_airy_Bi_e (double x, gsl_mode_t mode, gsl_sf_result * result)
These routines compute the Airy function @math{Bi(x)} with an accuracy specified by mode.

Function: double gsl_sf_airy_Ai_scaled (double x, gsl_mode_t mode)
Function: int gsl_sf_airy_Ai_scaled_e (double x, gsl_mode_t mode, gsl_sf_result * result)
These routines compute a scaled version of the Airy function @math{S_A(x) Ai(x)}. For @math{x>0} the scaling factor @math{S_A(x)} is @c{\$\exp(+(2/3) x^{3/2})\$} @math{\exp(+(2/3) x^(3/2))}, and is 1 for @math{x<0}.

Function: double gsl_sf_airy_Bi_scaled (double x, gsl_mode_t mode)
Function: int gsl_sf_airy_Bi_scaled_e (double x, gsl_mode_t mode, gsl_sf_result * result)
These routines compute a scaled version of the Airy function @math{S_B(x) Bi(x)}. For @math{x>0} the scaling factor @math{S_B(x)} is @c{\$\exp(-(2/3) x^{3/2})\$} @math{exp(-(2/3) x^(3/2))}, and is 1 for @math{x<0}.