Pré. | Proc. |
j0
Bessel function of order zero.SYNOPSIS:
double x, y, j0();
y = j0(x);
DESCRIPTION:
Returns a Bessel function of order zero of the argument. The domain is divided into the intervals [0, 5] and (5, infinity). In the first interval this rational approximation is used:
2 2
(w - r ) (w - r ) P (w) / Q (w)
1 2 3 8
2
where w = x and each r is a zero of the function.
In the second interval, the Hankel asymptotic expansion is employed with two rational functions of degree 6/6 and 7/7.
ACCURACY:
Absolute error:
arithmetic domain # trials peak rms
DEC 0, 30 10000 4.4e-17 6.3e-18
IEEE 0, 30 60000 4.2e-16 1.1e-16