returns the value of the sine integral
via the function name.
it is based on the Chebyshev expansion
is an implementation-dependent number,
There are no failure exits from s13adf
. The argument ifail
has been included for consistency with other routines in this chapter.
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
is always less than
, the accuracy will be limited by machine precision