s13adf returns the value of the sine integral
via the function name.
For
it is based on the Chebyshev expansion
For
, where
is an implementation-dependent number,
where
and
,
.
There are no failure exits from
s13adf. The argument
ifail has been included for consistency with other routines in this chapter.
If
and
are the relative errors in the argument and result, respectively, then in principle
The equality may hold if
is greater than the
machine precision (
due to data errors etc.) but if
is simply due to round-off in the machine representation, then since the factor relating
to
is always less than
, the accuracy will be limited by
machine precision.
None.