Заряд Нётер для лагранжиана с производными высшего порядка

Я пытаюсь найти заряд Нётер для симметрии

$x\rightarrow x+f\left(x\right)$

Это преобразование должно оставить действие инвариантным, поэтому

$$\begin{align*} dS&=S\left(x+f\left(x\right),\dots\right)-S\left(x\right)=0\\ &=\int dt\ \mathcal{L}\left(x+f\left(x\right),\dot{x}+\dot{f},\dots,t\right)-\mathcal{L}\left(x,\dot{x},\dots,t\right) \end{align*}$$

С использованием$f\left(x+\epsilon\right)-f\left(x\right)\approx \epsilon \frac{df}{dx}$

$$\begin{align*} dS=\int dt\ \frac{\partial \mathcal{L}}{\partial x}f+\frac{\partial \mathcal{L}}{\partial \dot{x}}\dot{f}+\frac{\partial \mathcal{L}}{\partial \ddot{x}}\ddot{f}+\dots \end{align*}$$ Запись второго члена в виде полной производной $$\begin{align*} dS&=\int dt\ \frac{\partial \mathcal{L}}{\partial x}f+\frac{d}{dt}\left[\frac{\partial L}{\partial \dot{x}} f\right]-f\frac{d}{dt}\left(\frac{\partial\mathcal{L}}{\partial\dot{x}}\right)+\frac{\partial\mathcal{L}}{\partial \ddot{x}}\ddot{f}+\dots\\ &=\frac{\partial\mathcal{L}}{\partial \dot{x}}f\bigg|_0^T+\int dt\ \frac{\partial \mathcal{L}}{\partial x}f-f\frac{d}{dt}\left(\frac{\partial\mathcal{L}}{\partial\dot{x}}\right)+\frac{\partial\mathcal{L}}{\partial \ddot{x}}\ddot{f}+\dots \end{align*}$$ Для членов более высокого порядка мы можем сделать то же самое $$\begin{align*} \frac{\partial \mathcal{L}}{\partial \ddot{x}}\ddot{f}&=\frac{d}{dt}\left(\frac{\partial \mathcal{L}}{\partial \ddot{x}}\dot{f}\right)-\dot{f}\frac{d}{dt}\left(\frac{\partial \mathcal{L}}{\partial \ddot{x}}\right)\\ &=\frac{d}{dt}\left(\frac{\partial \mathcal{L}}{\partial \ddot{x}}\dot{f}\right)-\boxed{\frac{d}{dt}\left(f\frac{d}{dt}\left(\frac{\partial\mathcal{L}}{\partial\ddot{x}}\right)\right)}+f\frac{d^2}{dt^2}\left(\frac{\partial\mathcal{L}}{d\ddot{x}}\right) \end{align*}$$ Итак, теперь интеграл становится $$\begin{align*} dS=\sum_{n=0}^N\frac{\partial\mathcal{L}}{\partial x^{\left(n+1\right)}}f^{(n)}\bigg|_0^T+\int dt\ f\left[\sum_{n=0}^N\left(-1\right)^n\frac{d^n}{dt^n}\left(\frac{\partial \mathcal{L}}{\partial x^{(n)}}\right)\right]-\boxed{\frac{d}{dt}\left(f\frac{d}{dt}\left(\frac{\partial\mathcal{L}}{\partial\ddot{x}}\right)\right)+\dots}=0\end{align*}$$

Где сумма под интегралом представляет уравнения Эйлера-Лагранжа для невозмущенного действия. Я как бы ожидаю, что термины в штучной упаковке также исчезнут, оставив только первый термин. Правильно ли я сделал, какие шаги я пропустил?