| Pré. | Proc. |
iv
Modified Bessel function of noninteger order.SYNOPSIS:
double v, x, y, iv();
y = iv(v, x);
DESCRIPTION:
Returns the modified Bessel function of order v of the argument. If x is negative, v must be integer-valued.
The function is defined as Iv(x) = Jv(ix). Here, it is computed in terms of the confluent hypergeometric function, according to the formula:
v -x
Iv(x) = (x/2) e hyperg(v+0.5, 2v+1, 2x) / gamma(v+1)
If v is a negative integer, then v is replaced by -v.
ACCURACY:
Tested at random points (v, x), with v between 0 and 30, x between 0 and 28.
Relative error:
arithmetic domain # trials peak rms
DEC 0,30 2000 3.1e-15 5.4e-16
IEEE 0,30 10000 1.7e-14 2.7e-15
Accuracy is diminished if v is near a negative integer.