double expint(double x); float expintf(float x); long double expintl(long double x);
Ei(x)=−∫∞−xe−ttdt
where x is x.