Sin-integral-func
evaluates
The approximation is based on several Chebyshev expansions.
For more information please review the s13adc function in the NAG document.