/* expnf, returns exponential integral,
 * n is order, 1 for first integral integral
 * xx is argument, must be positive
 */
double expnf( int n, double xx );