 Solution

The relativistic orbit-shape equation can be derived without dividing by $\dot{r}$ at a turning point. Substitute the norm constraint into the radial geodesic equation to eliminate $\dot{t}$ and ${\dot{r}}^2$. Using $f^{\prime }=2\mu /r^2$ givesFor a massive particle with proper time as parameter, $\kappa =c^2$. Put $u=1/r$, and use primes in this paragraph for derivatives with respect to $\varphi$. The constant specific angular momentum gives $\dot{\varphi }=h{u}^2$, so the chain rule yieldsSubstituting these expressions into the radial equation and dividing by $h^2{u}^2$ for a nonradial orbit givesThe extra $3\mu u^2$ term is responsible for Schwarzschild perihelion precession. It vanishes in the Newtonian weak-field limit $\mu /r\to 0$. The derivation uses a massive particle; for a null geodesic, $\kappa =0$ and the constant $GM/h^2$ on the right is absent.