s15asc calculates approximate values for the error function
for an array of arguments
, for
.
Let
be the root of the equation
(then
). For
the value of
is based on the following rational Chebyshev expansion for
:
where
denotes a rational function of degree
in the numerator and
in the denominator.
For
the value of
is based on a rational Chebyshev expansion for
: for
the value is based on the expansion
and for
it is based on the expansion
For each expansion, the specific values of
and
are selected to be minimal such that the maximum relative error in the expansion is of the order
, where
is the maximum number of decimal digits that can be accurately represented for the particular implementation (see
X02BEC).
For
there is a danger of setting underflow in
(the value of
is given in the
Users' Note for your implementation). For
,
s15asc returns
; for
it returns
.
None.
This example reads values of
x from a file, evaluates the function at each value of
and prints the results.