Consider the integral
$$\displaystyle \int_{0}^{2\pi} \int_{0}^{2\pi} (x^2 + y^2)^{-3/4}((b-x)^2 + (c-y)^2)^{-3/4}\mathrm{d}x\mathrm{d}y,$$
where $b, c \in \mathbb{Z}$ are non-zero constants. It seems like there should be an obvious substitution to make here, but I can't spot one that would turn the integral into something nicer. In particular, I'm having trouble finding a substitution that works with the constants. If this were an integral in one dimension, a substitution like $x = bt$ would work rather easily, given that we could just pull the constant out. In this case, however, the fact that we have two constants here makes it tricky. Converting to polar co-ordinates would certainly simplify the $(x^2 + y^2)^{-3/4}$ term but then dealing with the second term is still an issue, as are the constants. Can anyone recommend a substitution to use here?
For anyone wondering why I need this result, see here. The NIntegrate function in Mathematica quickly gives results if we plug in values for $b,c$, but I need a general result for $b,c$. I am also willing to accept a bound for the integral provided the resulting expression in $b,c$ is square-summable (see the link).