scipy.special.ive

scipy.special.ive(x1, x2[, out])

y=ive(v,z) returns the exponentially scaled modified Bessel function of real order v and complex z: ive(v,z) = iv(v,z) * exp(-abs(z.real))