I'm having trouble understanding how Gelfand & Fomin go from
$$ \begin{align} \Delta J[h] &= J[y+h] - J[y] \\ &= \int_a^b \left( F_y h + F_y' h' \right) \, dx + \frac{1}{2}\int_a^b \left( \bar{F}_{yy}h^2 + 2\bar{F}_{yy'}hh' + \bar{F}_{y'y'}h'^2 \right) \, dx \tag{1}\label{eq1} \end{align} $$
to
$$ \Delta J[h] = \int_a^b \left( F_y h + F_y' h' \right) \, dx + \frac{1}{2}\int_a^b \left( F_{yy}h^2 + 2F_{yy'}hh' + F_{y'y'}h'^2 \right) \, dx + \epsilon \tag{2}\label{eq2} $$
where
$$ \epsilon = \int_a^b \left( \epsilon_1h^2 + \epsilon_2hh' + \epsilon_3h'^2 \right) \, dx $$
in Section 25. Step \ref{eq1} is just Taylor's theorem with remainder from David Widder's Advanced Calculus (Section 9.2, Eq. 3). Going from \ref{eq1} to \ref{eq2} however, Gelfand and Fomin state:
"If we replace $\bar{F}_{yy}$, $\bar{F}_{yy'}$, and $\bar{F}_{y'y'}$ by the derivatives $F_{yy}$, $F_{yy'}$, and $F_{y'y'}$ evaluated at the point $(x, y(x), y'(x))$, then [$\ref{eq1}$] becomes [$\ref{eq2}$]."
How do they do this? (Continuity or some application of mean value theorem perhaps?)