integrate 2/pi^2 * 1/sqrt((r^2-x^2)(R^2-(R-r)^2)) for r from 0 to R