| Pré. | Proc. |
k1
Modified Bessel function, third kind, order one.SYNOPSIS:
double x, y, k1();
y = k1(x);
DESCRIPTION:
Computes the modified Bessel function of the third kind, of order one of the argument.
The range is partitioned into the two intervals [0,2] and (2, infinity). Chebyshev polynomial expansions are employed in each interval.
ACCURACY:
Relative error:
arithmetic domain # trials peak rms
DEC 0, 30 3300 8.9e-17 2.2e-17
IEEE 0, 30 30000 1.2e-15 1.6e-16
ERROR MESSAGES:
message condition value returned
domain x <= 0 MAXNUM