/*ee1 first exponential integral */
double ee1(double x);