Page:Laplace - Œuvres complètes, Gauthier-Villars, 1878, tome 11.djvu/385

371

THÉORIE DES SATELLITES DE JUPITER.

puissance de ces forces doivent disparaître d’eux-mêmes, on aura

${\begin{aligned}0=&{\frac {d^{2}\left(r\delta 'r+{\frac {1}{2}}\delta r^{2}\right)}{dt^{2}}}+{\frac {(1+m)r\delta 'r}{r^{3}}}-{\frac {(1+m)\delta r^{2}}{r^{3}}}\\&+2\int \operatorname {d} \mathrm {R} +x{\frac {\partial \mathrm {R} }{\partial x}}+y{\frac {\partial \mathrm {R} }{\partial y}}+z{\frac {\partial \mathrm {R} }{\partial z}},\\\\0=&r{\frac {d^{2}\delta 'r}{dt^{2}}}+\delta 'r{\frac {d^{2}r}{dt^{2}}}+\delta r{\frac {d^{2}\delta r}{dt^{2}}}-2r^{2}{\frac {dv}{dt}}{\frac {d\delta 'v}{dt}}-2r\delta 'r{\frac {d^{2}v}{dt^{2}}}\\&-r^{2}\left({\frac {d\delta v}{dt}}\right)^{2}-4r\delta r{\frac {dv}{dt}}{\frac {d\delta v}{dt}}-\delta r^{2}\left({\frac {dv}{dt}}\right)^{2}\\&-{\frac {(1+m)r\delta 'r}{r^{3}}}+{\frac {(1+m)\delta r^{2}}{r^{3}}}+x{\frac {\partial \mathrm {R} }{\partial x}}+y{\frac {\partial \mathrm {R} }{\partial y}}+z{\frac {\partial \mathrm {R} }{\partial z}},\end{aligned}}$

en observant de ne conserver dans les fonctions

$\int \operatorname {d} \mathrm {R} \quad {\text{et}}\quad x{\frac {\partial \mathrm {R} }{\partial x}}+y{\frac {\partial \mathrm {R} }{\partial y}}+z{\frac {\partial \mathrm {R} }{\partial z}}$

que les termes dépendants des carrés et des produits des forces perturbatrices.

Les propriétés du mouvement elliptique donnent

${\begin{aligned}{\frac {r^{2}dv}{dt}}=&{\sqrt {(1+m)a\left(1-e^{2}\right)}}\\{\frac {rd^{2}v}{dt^{2}}}=&{\frac {d^{2}r}{dt^{2}}}+{\frac {1+m}{r^{2}}}\,;\end{aligned}}$

en substituant ces valeurs dans la seconde des deux équations différentielles précédentes, et en y mettant, au lieu de ${\frac {(1+m)r\delta 'r}{r^{3}}},$ sa valeur tirée de la première, on aura

${\begin{aligned}{\frac {2d\delta 'v}{dt}}&{\sqrt {(1+m)a\left(1-e^{2}\right)}}\\=&{\frac {d(rd\delta 'r-\delta 'rdr)}{dt^{2}}}+3{\frac {d^{2}\left(r\delta 'r+{\frac {1}{2}}\delta r^{2}\right)}{dt^{2}}}\\&-{\frac {2(1+m)\delta r^{2}}{r^{3}}}+6\int \operatorname {d} \mathrm {R} +4\left(x{\frac {\partial \mathrm {R} }{\partial x}}+y{\frac {\partial \mathrm {R} }{\partial y}}+z{\frac {\partial \mathrm {R} }{\partial z}}\right)\\&+{\frac {\delta rd^{2}\delta r-r^{2}(d\delta v)^{2}-4rdv\delta rd\delta v-dv^{2}\delta r^{2}}{dt^{2}}}.\end{aligned}}$