/// <summary>
/// Airy function Ai(z).
/// <p/>
/// If expScaled is true, returns Exp(zta) * Ai(z), where zta = (2/3) * z * Sqrt(z).
/// </summary>
/// <param name="z">The value to compute the Airy function of.</param>
/// <param name="expScaled">If true, returns exponentially-scaled Airy function</param>
/// <returns></returns>
public static double AiryAi(double z, bool expScaled = false)
{
if (expScaled)
{
return(Amos.ScaledCairy(z));
}
else
{
return(AiryAi(new Complex(z, 0), expScaled).Real);
}
}
/// <summary>
/// Returns the exponentially scaled Airy function Ai.
/// <para>ScaledAiryAi(z) is given by Exp(zta) * AiryAi(z), where zta = (2/3) * z * Sqrt(z).</para>
/// </summary>
/// <param name="z">The value to compute the Airy function of.</param>
/// <returns>The exponentially scaled Airy function Ai.</returns>
public static double AiryAiScaled(double z)
{
return(Amos.ScaledCairy(z));
}
/// <summary>
/// Returns the exponentially scaled Airy function Ai.
/// <para>ScaledAiryAi(z) is given by Exp(zta) * AiryAi(z), where zta = (2/3) * z * Sqrt(z).</para>
/// </summary>
/// <param name="z">The value to compute the Airy function of.</param>
/// <returns>The exponentially scaled Airy function Ai.</returns>
public static Complex AiryAiScaled(Complex z)
{
return(Amos.ScaledCairy(z));
}
/// <summary>
/// Airy function Ai(z).
/// <p/>
/// If expScaled is true, returns Exp(zta) * Ai(z), where zta = (2/3) * z * Sqrt(z).
/// </summary>
/// <param name="z">The value to compute the Airy function of.</param>
/// <param name="expScaled">If true, returns exponentially-scaled Airy function</param>
/// <returns></returns>
public static Complex AiryAi(Complex z, bool expScaled = false)
{
return((expScaled) ? Amos.ScaledCairy(z) : Amos.Cairy(z));
}
