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