◆ Cbiry()

エアリー関数 Bi(z) またはその導関数 Bi'(z) (複素引数)

目的: 本ルーチンはエアリー関数 Bi(z) またはその導関数 Bi'(z) を求める.
Bi(z) = (1/π) * ∫ (sin((1/3)t^3 + zt) + exp(-(1/3)t^3 + zt)) dt [0, ∞]

オーバーフローを避けるためにスケーリング・オプションを使用することができる.

引数

[in]	Z	引数 z.
[in]	Id	導関数の次数. = 0: Bi(z) を求める. = 1: Bi'(z) を求める.
[out]	Y	Airy関数 Bi(z) またはその導関数 Bi'(z). Kode = 2 の場合, exp(-\|re((2/3)*z^(3/2))\|) でスケーリングする.
[out]	Info	= 0: 正常終了. = -2: パラメータ Id の誤り. (Id <> 0 かつ Id <> 1) = -5: パラメータ Kode の誤り. (Kode <> 1 かつ Kode <> 2) = 2: 精度警告. (精度が半分以下になった) = 3: 精度エラー. (正しく計算できない) = 4: Kode = 1においてre(Z)が大きすぎるためオーバーフローが起きた.
[in]	Kode	(省略可) スケーリングの指定. (省略時 = 1) = 1: スケーリングなし. = 2: 指数スケーリングあり.