integrate 1/(sqrt(v+2*G*M*(1/r-1/c)))dr