Solution
= Solution
Take endpoint-preserving variations $y+\varepsilon\eta$, with $\eta(a)=\eta(b)=0$. The first variation of the <functional> is $\delta I=\int_a^b(f_y\eta+f_{y'}\eta')dx$; <integration by parts> supplies the differential equation and its conserved quantity below.