scipy.special.kn

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

y=kn(n,x) returns the modified Bessel function of the second kind (sometimes called the third kind) for integer order n at x.