scipy.special.kn

scipy.special.kn(x1, x2[, out]) = <ufunc 'kn'>

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