I would like to numerically integrate the following function:
$I = \int^{a}_{-a}\int^{b}_{-b}Y_0(r)dx dy$
Where: $r = \sqrt{x^2+y^2}$, and $Y_0$ is a 0 order bessel function of the second kind.
But there is a singularity at r=0 so the calculation dosn't converge. Is there a way to deal with this sort of problem?