| Pré. | Proc. |
erf
Error function.SYNOPSIS:
double x, y, erf();
y = erf(x);
DESCRIPTION:
The integral is
x
-
2 | | 2
erf(x) = -------- | exp(- t) dt.
sqrt(pi) | |
-
0
The magnitude of x is limited to 9.231948545 for DEC arithmetic; 1 or -1 is returned outside this range.
For 0 <= |x| < 1, erf(x) = x * P4(x**2)/Q5(x**2); otherwise erf(x) = 1 - erfc(x).
ACCURACY:
Relative error:
arithmetic domain # trials peak rms
DEC 0,1 14000 4.7e-17 1.5e-17
IEEE 0,1 30000 3.7e-16 1.0e-16