int 1/(b+sqrt(r^2-x^2)) dx