(ij0,iy0)=itj0y0(x) returns simple integrals from 0 to x of the zeroth order Bessel functions j0 and y0.
scipy.special.k1e
scipy.special.it2j0y0