(1/(2 b) (Exp[a b] ExpIntegralE[1, a b] - Exp[-a b] ExpIntegralE[1, -a b]))/.{a->1,b->2.0}