SciPy

scipy.special.it2i0k0

scipy.special.it2i0k0(x[, out1, out2]) = <ufunc 'it2i0k0'>

(ii0,ik0)=it2i0k0(x) returns the integrals int((i0(t)-1)/t,t=0..x) and int(k0(t)/t,t=x..infinitity).