scipy.special.expm1

scipy.special.expm1(x[, out])

y=expm1(x) calculates exp(x) - 1 for use when x is near zero.