/* e2, second exponential integral, E2(x),
 * first variable is argument x,
 * second is exp(-x)
 */
double e2(double,double);