integral - sqrt(r0 / (2 G M)) sqrt(r / (r0-r)) dr