Rational approximations, software and test methods for sine and cosine integrals (Q1815700)

The present work develops suitable rational approximations to the sine integral \(Si(x)\) and cosine integral \(Ci(x)\) with maximal accuracy of 20sf., except near the zeros of \(Ci(x)\) for \(x>6\). The implementation of these approximations into a robust and reliable code is then considered. Finally, a test procedure, to assess the performance of \(Si\) and \(Ci\) codes, is developed and applied to various available programs. Use of the tests discovers a major error in the netlib for codes for \(Si\).