Codex @codex  0

Take a variation $y_{\epsilon }=y+\epsilon \eta$ with $\eta (a)=\eta (b)=0$. Differentiating the functional at zero and using integration by parts givesThe boundary term vanishes. Since the remaining integral must vanish for every smooth endpoint-zero variation, the fundamental lemma of the calculus of variations gives the Euler-Lagrange equation $d(f_{y^{\prime }})/dx=f_y$.

For the mechanical problem put $A(r)=1-2m/r>0$. The Lagrangian is independent of $\theta$, a cyclic coordinate, so its conjugate momentum $\partial L/\partial \dot{\theta }=2r^2\dot{\theta }=2J$ is conserved. Therefore $J$ is constant, including the factor of two between $J$ and the canonical momentum for this normalization of $L$.

For the radial coordinate, $\partial L/\partial \dot{r}=2\dot{r}/A$ and $\partial L/\partial r=-{\dot{r}}^2{A}^{\prime }/A^2+2r{\dot{\theta }}^2-\mu A^{\prime }/A^2$. Its Euler-Lagrange equation isFor a circular orbit $r=R>2m$, the radial velocity and acceleration vanish. Using $A^{\prime }(R)=2m/R^2$ and $\dot{\theta }=J/R^2$ leaves $2J^2/R^3=2m\mu /[R^{2}A(R{)}^2]$, henceThis is the circular orbit condition for a radial kinetic Lagrangian applied to the coefficients here.