On the exponential integral computation
Matematičeskoe modelirovanie, Tome 20 (2008) no. 1, pp. 87-91
Citer cet article
Voir la notice de l'article provenant de la source Math-Net.Ru
New high-precision algorithm for exponential integral calculation was developed. It is based on the representation of exponential integral in form of convergent series when $x$-argument is not large and in form of asymptotically convergent continued fraction when $x$ is large. It was shown that the optimal bound between these representations is $x=1$. At the same time using of 18 series members and 220 continued fraction members guarantees relative pre-cision lower than $2\cdot 10^{-15}$, that exceeds practical needs.
[1] E. S. Kuznetsov, Izbrannye nauchnye trudy, FIZMATLIT, M., 2003
[2] M. Abramowitz, I. A. Stegun (ed.), Handbook of mathematical functions with formulas, graphs, and mathematical tables, John Wiley Sons, New York, 1972 | MR | Zbl
[3] S. Chepmen, T. Kauling, Matematicheskaya teoriya neodnorodnykh gazov, Izdatelstvo inostrannoi literatury, M., 1960 | MR