Function File: y = expint_E1 (x)

Compute the exponential integral,

                        infinity
                       /
           expint(x) = | exp(t)/t dt
                       /
                      x

See also: expint expint_Ei