scipy.special.it2j0y0

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

(ij0,iy0)=it2j0y0(x) returns the integrals int((1-j0(t))/t,t=0..x) and int(y0(t)/t,t=x..infinitity).