Předpokládejme, že máme funkcionál
nazývaný akce, kde
je konfigurační prostor závisející obecně na všech veličinách
, kterými popisujeme daný systém. (Multiindex A tyto veličiny čísluje.)
Předpokládejme dále, že akci můžeme vyjádřit v obvyklém tvaru, jako integrál hustoty lagrangiánu

přes celý prostor
![{\displaystyle {\mathcal {S}}[\phi _{A}]\,=\,\int _{M}\mathrm {d} ^{n}x\,{\mathcal {L}}(x,\phi _{A}(x),\partial _{\mu }\phi _{A}(x)).}](//wikimedia.org/api/rest_v1/media/math/render/svg/dea02386ec788494460acb4490ce0ecb6b2709ea)
(Přitom předpokládáme, že lagrangián závisí jen na prvních derivacích zkoumaných proměnných. Pro vyšší derivace je zobecnění přímočaré.) Podle principu stacionární akce platí
![{\displaystyle {\frac {\delta {\mathcal {S}}[\phi _{A}]}{\delta \phi _{A}}}\,=\,0,}](//wikimedia.org/api/rest_v1/media/math/render/svg/563967126d554098b71a6ab60a9c9986508bd871)
a to tak, že na okraji M jsou veličiny
nulové. (Jde o úlohu s pevnými okraji.) Provedením variace
![{\displaystyle \delta {\mathcal {S}}[\phi ]\,=\,\int _{M}\mathrm {d} ^{n}x\,{\frac {\partial {\mathcal {L}}}{\partial \phi _{A}}}\delta \phi _{A}+{\frac {\partial {\mathcal {L}}}{\partial (\partial _{\mu }\phi _{A})}}\delta (\partial _{\mu }\phi _{A})=\,\int _{M}\mathrm {d} ^{n}x\,{\frac {\partial {\mathcal {L}}}{\partial \phi _{A}}}\delta \phi _{A}-\partial _{\mu }\left({\frac {\partial {\mathcal {L}}}{\partial (\partial _{\mu }\phi _{A})}}\right)\delta \phi _{A}=\,\int _{M}\mathrm {d} ^{n}x\,{\mathcal {L}}^{*}\delta \phi _{A},}](//wikimedia.org/api/rest_v1/media/math/render/svg/08e1c3c851d6de37fc164d64a8f281659c3a0b3d)
kde

Protože
jsou na M libovolné, musí platit
což jsou rovnice popisující daný systém.
Nyní mějme spojitou k-parametrickou transformaci souřadnic

vůči které jsou fyzikální zákony invariantní.
jsou obecně libovolné diferencovatelné funkce,
jsou infinitezimální parametry, které generují příslušnou Lieovu grupu symetrií a i probíhá hodnoty 1..k. Tato transformace indukuje transformaci zkoumaných veličin



Protože se při ní (z předpokladů věty) nemění tvar pohybových rovnic, platí

takže pro akci systému platí
![{\displaystyle {\mathcal {S'}}[\phi ]\,=\,{\mathcal {S}}[\phi ]\implies \,\int _{M'}\mathrm {d} ^{n}x'\,{\mathcal {L'}}(x',\phi '_{A},\partial _{\mu }\phi '_{A})\,=\,\int _{M}\mathrm {d} ^{n}x\,{\mathcal {L}}(x,\phi _{A},\partial _{\mu }\phi _{A}).\;\;\;\;\;\;\;\;\;\;\;\;\;(\spadesuit )}](//wikimedia.org/api/rest_v1/media/math/render/svg/53e3832ea7c590838f8cbda886c4120bd2006775)
Protože transformace symetrie je infinitezimální, můžeme nahradit integrování přes M' v souřadnicích x' integrováním přes M v souřadnicích x, pokud při tom zároveň přičteme povrchový člen, o který se liší na hranici M. Ten vypočteme jako plošný integrál hustoty lagrangiánu krát skalární součin normály hranice s
.

což dále upravíme podle Gaussovy věty a pro hraniční členy dosadíme
, (variace na hranici je nulová,) čímž obdržíme tvar

Dosadíme-li tento tvar do rovnice (♠), obdržíme tvar

což přepíšeme jako

kde jsme zavedli

je tzv. Lieova derivace
podle pole
. Veličina
se v užším kontextu, kde uvažujeme jako souřadnici jenom čas, označuje jako izochronní variace.
Diferencováním získáme

což lze pomocí integrace per partes a Gaussovy věty jako

Pokud se dá
vyjádřit jako
a má platit pro všechna
a zároveň platí
, získáme k rovnic

což jsou hledané zákony zachování.
Zákon zachování energie
Nyní uvažujme pouze soustavu hmotných bodů, které se nacházejí v potenciálu, který závisí jen na jejich vzájemné poloze. Lagrangián je tedy dán jako
![{\displaystyle {\mathcal {S}}[{\vec {x}}]\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/2f6adb06761bd005f020e7ce1a8f979ba86c49c3) |
 |
|
![{\displaystyle =\int \mathrm {d} t\left[\sum _{R=1}^{N}{\frac {m_{R}}{2}}\left({\dot {\vec {X}}}_{R}\right)^{2}-\sum _{R<S}V_{RS}({\vec {X}}_{S}-{\vec {X}}_{R})\right]}](//wikimedia.org/api/rest_v1/media/math/render/svg/f6a611fa87691c4011b7cb1ff81a0a59404b6d71) |
Roli souřadnice zde hraje jen čas t, zatímco polohy
hrají roli zkoumaných veličin výše označených jako
. Protože nás zajímá infinitezimální transformace spojená s posunem v čase, zvolíme



Z toho dopočteme

Povšimněme si, že
je nenulové, zatímco
jsou nulové - o symetrie spojené s posunem v prostoru se nyní nezajímáme. Zákon zachování tedy dostáváme ve tvaru

což po dalších úpravách
![{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left(\sum _{R=1}^{N}m_{R}{\dot {\vec {X}}}_{R}({\vec {\dot {X}}}_{R})-\left[\sum _{R=1}^{N}{\frac {m_{R}}{2}}\left({\dot {\vec {X}}}_{R}\right)^{2}-\sum _{R<S}V_{RS}({\vec {X}}_{S}-{\vec {X}}_{R})\right]\right)=0,}](//wikimedia.org/api/rest_v1/media/math/render/svg/33866421ade8f69149793aea87fac3ecac6d9975)

přejde na hledaný zákon zachování energie.