There is a step in what seems to be the standard derivation of the Euler-Lagrange equation that I do not see how can be justified (it is the first derivation in the Wikipedia article on the equation).
We have that $g(\epsilon) = f(t) + \epsilon \lambda(t)$
We are computing the total derivative $\frac{\partial F}{\partial\epsilon}(t,g(\epsilon),\dot{g}(\epsilon))$ and get that this is equal to $\frac{\partial F}{\partial g}\lambda(t) + \frac{\partial F}{\partial \dot{g}}\dot{\lambda}(t)$. The thing that does not make sense to me is that this is supposed to equal $\frac{\partial F}{\partial f}\lambda(t) + \frac{\partial F}{\partial \dot{f}}\dot{\lambda}(t)$ at $\epsilon = 0$.
To make some sense of it, I tried using the definition of the derivative on $\frac{\partial F}{\partial g}\lambda(t)$ to see if I can make it equal $\frac{\partial F}{\partial f}\lambda(t)$. However, what I do get is: $$\lim_{\Delta \epsilon \to 0}\frac{F(t,f(t)+(\epsilon+\Delta \epsilon) \lambda(t),\dot{g}) - F(t,f(t)+ \epsilon \lambda(t),\dot{g})}{\Delta \epsilon }$$ $$=\lim_{\Delta \epsilon \to 0}\frac{F(t,f(t)+(\epsilon+\Delta \epsilon) \lambda(t),\dot{g}) - F(t,f(t)+ \epsilon \lambda(t),\dot{g})}{\Delta g } \frac{\Delta g}{\Delta\epsilon}$$ Here we have that $\Delta g = \Delta \epsilon \lambda(t)$. This gives me: $$=\lim_{\Delta \epsilon \to 0}\frac{F(t,f(t)+(\epsilon+\Delta \epsilon) \lambda(t),\dot{g}) - F(t,f(t)+ \epsilon \lambda(t),\dot{g})}{\Delta g } \lambda(t)$$
Now I am just stuck. I cannot convince myself that $\Delta g = \Delta f$. I mean, the derivation is with respect to $\epsilon$, not t. So, why does this lead to the Euler-Lagrange equation?
Addendum
I think I may have to clarify my problem: $\frac{\partial F}{\partial f} \lambda$ looks to me to be equal to: $$\lim_{\Delta f \to 0}\frac{F(t,f(t)+\Delta f(t),\dot{g}) - F(t,f(t),\dot{g})}{\Delta f(t) } \lambda(t)$$
However, my computations above imply that the solution is $$\lim_{\Delta \epsilon \to 0}\frac{F(t,f(t)+(\epsilon+\Delta \epsilon) \lambda(t),\dot{g}) - F(t,f(t)+ \epsilon \lambda(t),\dot{g})}{\Delta g } \lambda(t) = \lim_{\Delta \epsilon \to 0}\frac{F(t,f(t)+\Delta \epsilon \lambda(t),\dot{g}) - F(t,f(t),\dot{g})}{\Delta \epsilon \lambda(t) } \lambda(t) $$
Clearly, my reasoning and computations are wrong. However, I don't see how I can get the solution above by computing $\frac{dF}{d\epsilon}$.