One obtains for a functional for single function of single variable with higher derivatives
$ I[f] = \int_{x_0}^{x_1} \mathcal{L}(x, f, f', f'')~\mathrm{d}x ~;~~ f' := \cfrac{\mathrm{d}f}{\mathrm{d}x}, ~f'' := \cfrac{\mathrm{d}^2f}{\mathrm{d}x^2} $
the Euler-Lagrange equation as $\cfrac{\partial \mathcal{L}}{\partial f} - \cfrac{\mathrm{d}}{\mathrm{d} x}\left(\cfrac{\partial \mathcal{L}}{\partial f'}\right) + \cfrac{\mathrm{d}^2}{\mathrm{d} x^2}\left(\cfrac{\partial \mathcal{L}}{\partial f''}\right) = 0 $
I want to prove the above equation extending Euler-Lagrange for several functions of one variable i.e.
For a given functional of the form :
$ I[f_1,f_2] = \int_{x_0}^{x_1} \mathcal{L}(x, f_1, f_2,f_1', f_2')~\mathrm{d}x ~;~~ f_i' := \cfrac{\mathrm{d}f_i}{\mathrm{d}x}$
the Euler-Lagrange Equation is
$\frac{\partial \mathcal{L}}{\partial f_i} - \frac{\mathrm{d}}{\mathrm{d}x}\left(\frac{\partial \mathcal{L}}{\partial f_i'}\right) = 0 , i=1,2$
This can be achieved by defining $f_1 = f; f_2 =f'$.
But I'm not able to achieve the same equation as for single function of single variable with higher derivatives.