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