integral (1 - (r)/x)^(-1/2) dx