scipy.special.it2i0k0

scipy.special.it2i0k0(x[, out1, out2])

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