integrate tan(arccos(1-(GM/(rc^2)))) dr