scipy.special.iv

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

y=iv(v,z) returns the modified Bessel function of real order v of z. If z is of real type and negative, v must be integer valued.