— 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 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 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 S_A(x) Ai(x). For x>0 the scaling factor S_A(x) is \exp(+(2/3) x^(3/2)), and is 1 for 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 S_B(x) Bi(x). For x>0 the scaling factor S_B(x) is exp(-(2/3) x^(3/2)), and is 1 for x<0.