integrate r^2 / (sqrt(1 - a/r))