j’ai fondé, dans le no 12 du Livre cinquième de la Mécanique céleste, la démonstration de ce théorème.
Considérons la Terre et la Lune comme formant un système. Si l’on prend son centre de gravité pour origine des coordonnées ; si l’on nomme
les coordonnées du centre de gravité de la Terre,
celles d’une quelconque de ses molécules
les coordonnées de la Lune et
la masse de la Terre, on aura
![{\displaystyle {\begin{aligned}c=&\int dm{\frac {x_{1}dy_{1}-y_{1}dx_{1}}{dt}}+\mathrm {L} {\frac {xdy-ydx}{dt}}\\=&\int dm{\frac {\left(x_{1}-\mathrm {X} \right)\left(dy_{1}-d\mathrm {Y} \right)-\left(y_{1}-\mathrm {Y} \right)\left(dx_{1}-d\mathrm {X} \right)}{dt}}\\&+\mathrm {L} {\frac {\left(x-\mathrm {X} \right)\left(dy-d\mathrm {Y} \right)-\left(y-\mathrm {Y} \right)\left(dx-d\mathrm {X} \right)}{dt}}\\&-(\mathrm {T+L} ){\frac {\mathrm {X} d\mathrm {Y} -\mathrm {Y} d\mathrm {X} }{dt}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4427596d9a6037b247c58dd2017ab43fbf4ed8da)
car on a, par la propriété du centre de gravité du système,
![{\displaystyle {\begin{alignedat}{2}\int x_{1}dm+\mathrm {L} x=&0,\qquad &\int dm{\frac {dx_{1}}{dt}}+\mathrm {L} {\frac {dx}{dt}}=&0,\\\int y_{1}dm+\mathrm {L} y=&0,&\int dm{\frac {dy_{1}}{dt}}+\mathrm {L} {\frac {dy}{dt}}=&0,\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f2430f071c445cbde4d7e70b8112b3185e02883)
étant la vitesse de rotation de la Terre et
désignant la somme des produits des molécules terrestres par le carré de leurs distances à l’axe de rotation, on a
![{\displaystyle \int dm{\frac {\left(x_{1}-\mathrm {X} \right)\left(dy_{1}-d\mathrm {Y} \right)-\left(y_{1}-\mathrm {Y} \right)\left(dx_{1}-d\mathrm {X} \right)}{dt}}=\mathrm {C} n\cos \theta ',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83b9039896d12bac56da97a839d25fc4eb91b05f)
étant l’inclinaison de l’équateur terrestre au plan fixe. On a pareillement
![{\displaystyle \mathrm {L} {\frac {\left(x-\mathrm {X} \right)\left(dy-d\mathrm {Y} \right)-\left(y-\mathrm {Y} \right)\left(dx-d\mathrm {X} \right)}{dt}}=\mathrm {L} n'a^{2}{\sqrt {1-e^{2}}}\cos \gamma ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3550603521166ebf2d9ef9a333cdd7a27b6626c0)
étant la vitesse moyenne angulaire du mouvement de la Lune autour de la Terre,
étant sa moyenne distance à la Terre,
étant le rapport de l’excentricité au demi-grand axe de son orbite, et
étant