Computes the exponential integral E_1(x),
E_1(x) := Re \int_1^\infty dt \exp(-xt)/t.
err contains an estimate of the absolute error in the value z.
This function is from the GNU Scientific Library, see http://www.gnu.org/software/gsl/ for documentation.
Package: gsl