rexpm1 function

TOMS 708 Approximation REXP(x) of expm1(x) = exp(x) - 1