(ii0,ik0)=iti0k0(x) returns simple integrals from 0 to x of the zeroth order modified Bessel functions i0 and k0.
scipy.special.it2j0y0
scipy.special.it2i0k0