y=i1(x) returns the modified Bessel function of the second kind (sometimes called the third kind) of order 1 at x.
scipy.special.k0e
scipy.special.k1e