Ei_3(x)

Function: double gsl_sf_expint_3 (double x)
Function: int gsl_sf_expint_3_e (double x, gsl_sf_result * result)
These routines compute the exponential integral Ei_3(x) = \int_0^x dt \exp(-t^3) for x >= 0.