scipy.special.itj0y0

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

(ij0,iy0)=itj0y0(x) returns simple integrals from 0 to x of the zeroth order bessel functions j0 and y0.