Airyai Airyai_sub Airyaid Airyaid_sub Airybi Airybi_sub Airybid Airybid_sub Cairy Cbiry

Cbiry() Sub Cbiry ( Z As  Complex, Id As  Long, Y As  Complex, Info As  Long, Optional Kode As  Long = 1  ) Computes the Airy function Bi(z) or its derivative Bi'(z) (complex argument) Purpose This routine computes the Airy function Bi(z) or its derivative Bi'(z). Bi(z) = (1/π) * ∫ (sin((1/3)t^3 + zt) + exp(-(1/3)t^3 + zt)) dt [0, ∞] A scaling option is available to help avoid overflow. Parameters [in] Z Argument z. [in] Id Order of derivative. = 0: Bi(z) is computed. = 1: Bi'(z) is computed. [out] Y Airy function Bi(z) or its derivative Bi'(z). If Kode = 2, Y is scaled by exp(-|re((2/3)*z^(3/2))|). [out] Info = 0: Successful exit. = -2: The argument Id had an illegal value. (Id <> 0 and Id <> 1) = -5: The argument Kode had an illegal value. (Kode <> 1 and Kode <> 2) = 2: Precision warning. (Half precision or less) = 3: Precision error. (No precision) = 4: Overflow re(Z) too large on Kode = 1. [in] Kode (Optional) A parameter to indicate the scaling option. (default = 1) = 1: No scaling. = 2: Returns exponentially scaled values. Reference SLATEC