Cos-integral-func
evaluates
where denotes Euler's constant. The approximation is based on several Chebyshev expansions.
For more information please review the s13acc function in the NAG document.