RECHERCHES
SUR LES
INÉGALITÉS DES SATELLITES DE JUPITER
CAUSÉES PAR LEUR ATTRACTION MUTUELLE.
Multùm adhuc restat operis.
Sen., Epist. 64.
(Prix de l’Académie Royale des Sciences de Paris, tome IX, 1766.)
I.
Soient nommés :

le rayon vecteur de l’orbite d’un satellite quelconque projetée sur le plan de l’orbite de Jupiter ;

la tangente de la latitude du satellite par rapport à ce même plan ;

la force que Jupiter exerce sur le satellite à la distance 1.
On aura la distance du satellite au plan de l’orbite de Jupiter égale à

Donc la distance du satellite au centre de Jupiter sera
Par conséquent la force par laquelle le satellite est poussé vers Jupiter sera

Cette force peut être regardée comme composée de deux autres : l’une parallèle au rayon vecteur et égale à
l’autre perpendiculaire au plan de l’orbite de Jupiter et égale à
Or on peut, en général, réduire les forces perturbatrices du satellite à trois forces uniques, dont :
La première, que j’appelle
soit parallèle au rayon
La seconde, que j’appelle
soit perpendiculaire au rayon vecteur, et parallèle au plan de l’orbite de Jupiter ;
La troisième, que j’appelle
soit perpendiculaire à ce même plan.
Donc le satellite sera sollicité, dans les directions dont nous parlons, par les forces

dont les deux premières déterminent le mouvement que le satellite doit avoir dans le plan de l’orbite de Jupiter, ou pour mieux dire, parallèlement à ce plan.
II.
Cela posé, soient
le temps écoulé depuis le commencement du mouvement ;
l’angle décrit par le rayon
durant ce temps ; l’élément du temps
constant, c’est-à-dire, 
On aura pour la vitesse circulatoire du satellite, parallèlement au plan de l’orbite de Jupiter,
d’où résulte la force centrifuge
laquelle étant retranchée de la force
on aura la véritable force qui tend à diminuer le rayon
.
Donc, par le principe des forces accélératrices, on aura
(A)
|
|
|
Maintenant on sait que, si la force perpendiculaire

était nulle, le rayon

décrirait des aires proportionnelles aux temps, de sorte que l’on aurait, à cause de

constant,

mais la force
fait parcourir perpendiculairement à
l’espace
pendant le temps
donc le secteur
croîtra pendant ce temps de la quantité
par conséquent on aura l’équation

dont l’intégrale, en ajoutant
est

d’où l’on tire
(B)
|
|
|
Enfin on aura, en vertu de la force perpendiculaire au plan de l’orbite de Jupiter,

ou bien

d’où, en mettant pour
sa valeur

tirée de l’équation (A), et effaçant ce qui se détruit, on aura l’équation suivante
(C)
|
|
|
III.
Les équations (A), (B), (C) donneront
et
en
ce qui suffira pour laire connaître le lieu du satellite à chaque instant. Que si l’on voulait connaître la figure même de l’orbite qu’il décrit, il faudrait éliminer des équations (A), (C) l’élément
Or de l’équation (B) on tire, après quelques réductions fort simples,

donc, si l’on substitue cette valeur dans (A), (C), et qu’on fasse pour plus de simplicité
on aura, en prenant
constant,

Supposons pour un moment que les forces perturbatrices
soient nulles ; on aura par l’équation

dont l’intégrale est, comme on sait,

ou bien

et
étant deux constantes arbitraires. Cette dernière expression de
fait voir que l’orbite est toute dans un plan fixe, dont la position dépend des quantités
qui expriment, la première, la tangente de l’inclinaison, et la seconde, la longitude du nœud. Retenons maintenant cette même expression de
et supposons, à cause des forces perturbatrices
et
variables ; on aura

Or, afin que le corps puisse être regardé comme se mouvant réellement dans le plan déterminé par
et
il faut que la valeur de
soit la même que si ces quantités demeuraient constantes, c’est-à-dire, que

donc

par conséquent, à cause de
constant,

et

On réduira ainsi l’équation
ci-dessus à deux équations du premier degré, qui donneront
et
en
d’où l’on connaîtra la variation de l’inclinaison de l’orbite et le mouvement de la ligne des nœuds. C’est ainsi que la plupart des Géomètres en ont usé jusqu’ici dans la recherche des orbites des Planètes ; mais il nous paraît plus court de chercher directement la latitude
par une seule équation, d’autant plus que les quantités
et
s’en déduiront plus aisément ; car, puisque

on aura

On pourrait faire une pareille transformation sur l’équation
ce qui réduirait l’orbite à une ellipse dont l’excentricité et la position de la ligne des apsides seraient variables, ainsi que M. Newton l’a pratiqué par rapport à la Lune. En effet, si l’on suppose d’abord

l’équation
devient

dont l’intégrale, étant mise sous cette forme

donne une ellipse dans laquelle
est le demi-paramètre,
l’excentricité, et
la longitude de l’apside inférieure. Qu’on regarde maintenant
et
comme variables, et qu’on suppose, par une raison analogue à celle que nous avons expliquée ci-dessus,

on trouvera


Ainsi l’équation
se réduira à deux équations du premier degré, d’où l’on tirera aisément
et 
IV.
Les observations nous apprennent que les inégalités des mouvements des satellites de Jupiter sont très-petites, aussi bien que les inclinaisons de leurs orbites, par rapport à l’orbite de cette Planète ; d’où il suit que, si l’on nomme
la valeur moyenne de
la valeur moyenne de
c’est-à-dire la vitesse angulaire moyenne, et qu’on dénote par
un coefficient très-petit, et par
des quantités variables, on aura les expressions suivantes

où l’on remarquera que les valeurs de

et de

ne doivent contenir aucun terme constant ; autrement

et

ne seraient plus les valeurs moyennes de

et de

ce qui est contre l’hypothèse.
V.
Substituons maintenant ces expressions de
dans les équations de l’Article II et négligeons les termes qui se trouveraient multipliés par des puissances de
plus hautes que
parce qu’une plus grande exactitude serait superflue dans le sujet que nous traitons ; nous changerons d’abord l’équation (A) en celle-ci

ou bien

Si
était
on aurait

donc,
étant très-petite, la quantité

devra l’être aussi ; de sorte qu’on pourra supposer

Cette substitution faite, on divisera toute l’équation par

et l’on aura, en mettant pour plus de simplicité

au lieu de

VI.
L’équation (B) deviendra, par les mêmes substitutions,
![{\displaystyle \mu +n{\frac {dy}{dt}}=\left[{\frac {c}{a^{2}}}+{\frac {\int \mathrm {Q} (1+nx)dt}{a}}\right]\left(1-2nx+3n^{2}x^{2}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/03dce5386d2bdd098cffa87b145810b397d06b89)
Si
était
on aurait

supposons donc

on aura, après les réductions,

VII.
Enfin l’équation (C) se changera en celle-ci

et l’on prouvera ici, comme on a fait ci-dessus, qu’il faut que la quantité
soit très-petite de l’ordre
c’est pourquoi nous supposerons

d’où nous aurons l’équation

VIII.
Voilà les formules par lesquelles on pourra déterminer les inégalités des satellites de Jupiter, dès qu’on aura trouvé les valeurs des quantités
qui résultent de leur action mutuelle.
Pour rendre ces formules encore plus commodes pour le calcul, nous substituerons dans celles des Articles V et VII la valeur de
tirée de l’Article VI.
De cette manière, on aura, en négligeant toujours les termes affectés de
(D)
|
|
|
(E)
|
|
|
(F)
|
|
|
Nous avons dit que les valeurs de
et
ne doivent renfermer aucun terme constant ; on remplira ces deux conditions par le moyen des constantes
et
CHAPITRE II.
DÉTERMINATION DES FORCES PERTURBATRICES DES SATELLITES DE JUPITER.
IX.
Soient
la masse de Jupiter,
la masse du premier satellite,
la masse du second satellite,
la masse du troisième satellite,
la masse du quatrième satellite.
Supposons de plus que toutes les quantités que nous avons nommées

dans le Chapitre précédent, soient désignées ici, relativement au premier satellite, par

relativement au second satellite, par

relativement au troisième, par

et relativement au quatrième, par

En général, nous conserverons toujours dans la suite les noms donnés dans les Articles précédents, avec cette seule différence que nous marquerons les lettres d’un trait pour le premier satellite, de deux traits pour le second satellite, etc.[1].
Enfin nous dénoterons, pour plus de simplicité, la distance entre deux satellites quelconques, c’est-à-dire, la ligne droite qui joint leurs centres, par
étant les rayons vecteurs des deux satellites ; ainsi la distance entre le premier et le second satellite sera désignée par
la distance entre le premier et le troisième par
et ainsi des autres.
X.
Cela posé, il est visible
1o Que le satellite
est attiré vers Jupiter avec une force

et qu’en même temps Jupiter est attiré lui-même vers le satellite avec une force

d’où il suit que la force totale qui tend à rapprocher le satellite de Jupiter est

Cette expression doit être comparée avec l’expression de la force centrale
(Article I), c’est-à-dire, en la rapportant au premier satellite, avec
ce qui donne d’abord

2o Que le satellite
est attiré vers le satellite
avec une force

laquelle peut se décomposer en deux autres l’une dans la direction du rayon mené du satellite
à Jupiter, qui sera

l’autre parallèle au rayon mené du satellite
à Jupiter, et qui sera

De plus le même satellite
doit être regardé comme attiré par une force égale, et en sens contraire, à celle avec laquelle Jupiter est attiré par le satellite
c’est-à-dire par une force

et dirigée parallèlement au rayon mené de ce dernier satellite à Jupiter.
Donc l’action du satellite

produit dans le satellite

deux forces l’une

dirigée vers Jupiter, l’autre
![{\displaystyle {\mathfrak {S}}_{2}\left[{\frac {1}{r_{2}^{2}\left(1+p_{2}^{2}\right)}}-{\frac {r_{2}{\sqrt {1+p_{1}^{2}}}}{\Delta (r_{1},r_{2})^{3}}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ca7367870a82b984d8064fcc256d21fad203eea)
dans une direction parallèle à celle qui va du satellite
à Jupiter.
3o Or la force

se décompose en deux autres l’une perpendiculaire au plan de l’orbite de Jupiter

l’autre parallèle au même plan dans la direction du rayon
qui sera

Pareillement la force
![{\displaystyle {\mathfrak {S}}_{2}\left[{\frac {1}{r_{2}^{2}\left(1+p_{2}^{2}\right)}}-{\frac {r_{2}{\sqrt {1+p_{2}^{2}}}}{\Delta (r_{1},r_{2})^{3}}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c0f7a90a302e4be989943e878799cc1feb60f58f)
se change en deux autres forces : l’une perpendiculaire au plan de l’orbite de Jupiter
![{\displaystyle {\mathfrak {S}}_{2}\left[{\frac {p_{2}}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}-{\frac {r_{2}p_{2}}{\Delta (r_{1},r_{2})^{3}}}\right],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1441358255a666b2e12444f81bc3e649bad56287)
et l’autre parallèle à ce plan, dans la direction du rayon 
![{\displaystyle {\mathfrak {S}}_{2}\left[{\frac {1}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}-{\frac {r_{2}}{\Delta (r_{1},r_{2})^{3}}}\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10aeab453cc0f3764995d9f07e09d0e873e8dc62)
Enfin cette dernière force se décompose encore en deux autres : l’une dans la direction du rayon
avec le rayon
fait l’angle
l’autre
perpendiculaire à cette direction la première sera exprimée par
![{\displaystyle {\mathfrak {S}}_{2}\left[{\frac {1}{r_{2}^{2}\left(1+p_{2}^{2}\right)}}-{\frac {r_{2}}{\Delta (r_{1},r_{2})^{3}}}\right]\cos(\varphi _{2}-\varphi _{1}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/715d83448c1100a750377d23a44aaa4a6b44fa91)
la seconde par
![{\displaystyle {\mathfrak {S}}_{2}\left[{\frac {1}{r_{2}^{2}\left(1+p_{2}^{2}\right)}}-{\frac {r_{2}}{\Delta (r_{1},r_{2})^{3}}}\right]\sin(\varphi _{2}-\varphi _{1}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fb828cb30487a313014931e3a7818751309906f5)
et tendra à diminuer l’angle
au lieu que nous avons supposé (Article II) que la force perpendiculaire
tendait à augmenter l’angle
c’est pourquoi il faudra la prendre négativement.
4o Comparant donc toutes ces forces avec les forces
(Article I), ou bien
(Article IX), on aura en conséquence de l’action du satellite
les expressions suivantes
![{\displaystyle {\begin{aligned}\mathrm {R} _{1}=&{\mathfrak {S}}_{2}\left[{\frac {r_{1}-r_{2}\cos(\varphi _{2}-\varphi _{1})}{\Delta (r_{1},r_{2})^{3}}}+{\frac {\cos(\varphi _{2}-\varphi _{1})}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right],\\\mathrm {Q} _{1}=&{\mathfrak {S}}_{2}\left[{\frac {r_{2}\sin(\varphi _{2}-\varphi _{1})}{\Delta (r_{1},r_{2})^{3}}}-{\frac {\sin(\varphi _{2}-\varphi _{1})}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right],\\\mathrm {P} _{1}=&{\mathfrak {S}}_{2}\left[{\frac {r_{1}p_{1}-r_{2}p_{2}}{\Delta (r_{1},r_{2})^{3}}}+{\frac {p_{2}}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6b1a00b2e5ea268c4f36522c7a1383f7faf09bcd)
On trouvera de la même manière les expressions des forces
en tant qu’elles résultent de l’action des satellites
et
et il est clair que l’on aura les mêmes formules que ci-dessus, en marquant seulement de trois traits ou de quatre traits les lettres qui sont marquées de deux traits[2].
XI.
Si l’on veut avoir égard aussi à l’action du Soleil sur le satellite
on nommera :
la masse du Soleil,
la distance du satellite
au Soleil,
le rayon vecteur de l’orbite du Soleil autour de Jupiter,
la longitude du Soleil vu du centre de
:
et il n’y aura qu’à mettre, dans les expressions de

de l’Article X,

au lieu de

au lieu de

au lieu de

au lieu de

et supposer

De cette manière on aura, en vertu de l’action du Soleil,
![{\displaystyle {\begin{aligned}\mathrm {R} _{1}=&\circledast \left[{\frac {r_{1}-\rho _{1}\cos(\psi -\varphi _{1})}{\delta _{1}^{3}}}+{\frac {\cos(\psi -\varphi _{1})}{\rho _{1}^{2}}}\right],\\\mathrm {Q} _{1}=&\circledast \left[{\frac {\rho _{1}\sin(\psi -\varphi _{1})}{\delta _{1}^{3}}}-{\frac {\sin(\psi -\varphi _{1})}{\rho _{1}^{2}}}\right],\\\mathrm {P} _{1}=&\circledast {\frac {r_{1}p_{1}}{\delta _{1}^{3}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/725ec90e11cf1f31204eadc855e4f66c58854aef)
XII.
Donc, en joignant ensemble les forces qui proviennent de l’action des trois satellites
et du Soleil sur le satellite
on aura les valeurs complètes de
exprimées de la manière suivante
![{\displaystyle {\begin{aligned}\mathrm {R} _{1}&={\mathfrak {S}}_{2}\left[{\frac {r_{1}-r_{2}\cos(\varphi _{2}-\varphi _{1})}{\Delta (r_{1},r_{2})^{3}}}+{\frac {\cos(\varphi _{2}-\varphi _{1})}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{3}\left[{\frac {r_{1}-r_{3}\cos(\varphi _{3}-\varphi _{1})}{\Delta (r_{1},r_{3})^{3}}}+{\frac {\cos(\varphi _{3}-\varphi _{1})}{r_{3}^{2}\left(1+p_{3}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{4}\left[{\frac {r_{1}-r_{4}\cos(\varphi _{4}-\varphi _{1})}{\Delta (r_{1},r_{4})^{3}}}+{\frac {\cos(\varphi _{4}-\varphi _{1})}{r_{4}^{2}\left(1+p_{4}^{2}\right)^{\frac {3}{2}}}}\right]\\&+\circledast \ \ \left[{\frac {r_{1}-\rho _{1}\cos(\psi \ -\varphi _{1})}{\delta _{1}^{3}}}+{\frac {\cos(\ \psi \ -\varphi _{1})}{\rho _{1}^{2}}}\right],\\\\\mathrm {Q} _{1}&={\mathfrak {S}}_{2}\left[{\frac {r_{2}\sin(\varphi _{2}-\varphi _{1})}{\Delta (r_{1},r_{2})^{3}}}-{\frac {\sin(\varphi _{2}-\varphi _{1})}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{3}\left[{\frac {r_{3}\sin(\varphi _{3}-\varphi _{1})}{\Delta (r_{1},r_{3})^{3}}}-{\frac {\sin(\varphi _{3}-\varphi _{1})}{r_{3}^{2}\left(1+p_{3}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{4}\left[{\frac {r_{4}\sin(\varphi _{4}-\varphi _{1})}{\Delta (r_{1},r_{4})^{3}}}-{\frac {\sin(\varphi _{4}-\varphi _{1})}{r_{4}^{2}\left(1+p_{4}^{2}\right)^{\frac {3}{2}}}}\right]\\&+\circledast \ \ \left[{\frac {r_{4}\sin(\psi \ -\varphi _{1})}{\delta _{1}^{3}}}-{\frac {\sin(\ \psi \ -\varphi _{1})}{\rho _{1}^{2}}}\right],\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ceeea161559c59cd68af543763b33a020c5d58cc)
![{\displaystyle {\begin{aligned}\mathrm {P} _{1}&={\mathfrak {S}}_{2}\left[{\frac {r_{1}p_{1}-r_{2}p_{2}}{\Delta (r_{1},r_{2})^{3}}}+{\frac {p_{2}}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{3}\left[{\frac {r_{1}p_{1}-r_{3}p_{3}}{\Delta (r_{1},r_{3})^{3}}}-{\frac {p_{3}}{r_{3}^{2}\left(1+p_{3}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{4}\left[{\frac {r_{1}p_{1}-r_{4}p_{4}}{\Delta (r_{1},r_{4})^{3}}}-{\frac {p_{4}}{r_{4}^{2}\left(1+p_{4}^{2}\right)^{\frac {3}{2}}}}\right]\\&+\circledast {\frac {r_{1}p_{1}}{\delta _{1}^{3}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d024aa0e55a68b9353584393921bcfcb9565a5d)
XIII.
Telles sont les expressions des forces perturbatrices du satellite
d’où il est facile de déduire celles des trois autres satellites
En effet, un peu de réflexion suffit pour faire voir que les quantités
deviendront
en marquant seulement de deux traits les lettres qui sont marquées d’un trait, et réciproquement[3] ; ainsi l’on aura pour les forces perturbatrices du second satellite les expressions suivantes
![{\displaystyle {\begin{aligned}\mathrm {R} _{2}&={\mathfrak {S}}_{1}\left[{\frac {r_{2}-r_{1}\cos(\varphi _{1}-\varphi _{2})}{\Delta (r_{2},r_{1})^{3}}}+{\frac {\cos(\varphi _{1}-\varphi _{2})}{r_{1}^{2}\left(1+p_{1}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{3}\left[{\frac {r_{2}-r_{3}\cos(\varphi _{3}-\varphi _{2})}{\Delta (r_{2},r_{3})^{3}}}+{\frac {\cos(\varphi _{3}-\varphi _{2})}{r_{3}^{2}\left(1+p_{3}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{4}\left[{\frac {r_{2}-r_{4}\cos(\varphi _{4}-\varphi _{2})}{\Delta (r_{2},r_{4})^{3}}}+{\frac {\cos(\varphi _{4}-\varphi _{2})}{r_{4}^{2}\left(1+p_{4}^{2}\right)^{\frac {3}{2}}}}\right]\\&+\circledast \ \ \left[{\frac {r_{2}-\rho _{1}\cos(\psi \ -\varphi _{2})}{\delta _{2}^{3}}}+{\frac {\cos(\ \psi \ -\varphi _{2})}{\rho _{1}^{2}}}\right],\\\\\mathrm {Q} _{2}&={\mathfrak {S}}_{1}\left[{\frac {r_{1}\sin(\varphi _{1}-\varphi _{2})}{\Delta (r_{2},r_{1})^{3}}}-{\frac {\sin(\varphi _{1}-\varphi _{2})}{r_{1}^{2}\left(1+p_{1}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{3}\left[{\frac {r_{3}\sin(\varphi _{3}-\varphi _{2})}{\Delta (r_{2},r_{3})^{3}}}-{\frac {\sin(\varphi _{3}-\varphi _{2})}{r_{3}^{2}\left(1+p_{3}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{4}\left[{\frac {r_{4}\sin(\varphi _{4}-\varphi _{2})}{\Delta (r_{2},r_{4})^{3}}}-{\frac {\sin(\varphi _{4}-\varphi _{2})}{r_{4}^{2}\left(1+p_{4}^{2}\right)^{\frac {3}{2}}}}\right]\\&+\circledast \ \ \left[{\frac {r_{4}\sin(\psi \ -\varphi _{2})}{\delta _{2}^{3}}}-{\frac {\sin(\ \psi \ -\varphi _{2})}{\rho _{1}^{2}}}\right],\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2b0a9a6d3678b482189730bb55294d134161c3a1)
![{\displaystyle {\begin{aligned}\mathrm {P} _{2}&={\mathfrak {S}}_{1}\left[{\frac {r_{2}p_{2}-r_{1}p_{1}}{\Delta (r_{2},r_{1})^{3}}}+{\frac {p_{1}}{r_{1}^{2}\left(1+p_{1}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{3}\left[{\frac {r_{2}p_{2}-r_{3}p_{3}}{\Delta (r_{2},r_{3})^{3}}}-{\frac {p_{3}}{r_{3}^{2}\left(1+p_{3}^{2}\right)^{\frac {3}{2}}}}\right]\\&+{\mathfrak {S}}_{4}\left[{\frac {r_{2}p_{2}-r_{4}p_{4}}{\Delta (r_{2},r_{4})^{3}}}-{\frac {p_{4}}{r_{4}^{2}\left(1+p_{4}^{2}\right)^{\frac {3}{2}}}}\right]\\&+\circledast {\frac {r_{2}p_{2}}{\delta _{2}^{3}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eaed0cdffa032e91a0920295c7de45d0a7c6b63a)
On aura pareillement les expressions de
et de
en marquant successivement de trois et de quatre traits les lettres qui ne sont marquées que d’un seul trait dans les expressions de
et réciproquement[4].
XIV.
Il reste à chercher les valeurs des quantités
qui expriment les distances entre le premier satellite et le second, entre le premier et le troisième, etc. (Article IX). Or il est facile de trouver qu’on aura
![{\displaystyle {\begin{aligned}\Delta (r_{1},r_{2})^{2}=&(r_{2}p_{2}-r_{1}p_{1})^{2}+\left[r_{1}\sin(\varphi _{2}-\varphi _{1})\right]^{2}+\left[r_{2}-r_{1}\cos(\varphi _{2}-\varphi _{1})\right]^{2}\\=&r_{1}^{2}\left(1+p_{1}^{2}\right)-2r_{1}r_{2}\left[\cos(\varphi _{2}-\varphi _{1})+p_{1}p_{2}\right]+r_{2}^{2}\left(1+p_{2}^{2}\right)\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/286f86c9eb415e5e29b4a602d094ac65dd40dbd9)
donc, tirant la racine carrée,
![{\displaystyle \Delta (r_{1},r_{2})={\sqrt {r_{1}^{2}\left(1+p_{1}^{2}\right)-2r_{1}r_{2}\left[\cos(\varphi _{2}-\varphi _{1})+p_{1}p_{2}\right]+r_{2}^{2}\left(1+p_{2}^{2}\right)}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e1fed9c7428f13f26da580c4af5cff39a6c7f83)
On trouvera pareillement
![{\displaystyle \Delta (r_{1},r_{3})={\sqrt {r_{1}^{2}\left(1+p_{1}^{2}\right)-2r_{1}r_{3}\left[\cos(\varphi _{3}-\varphi _{1})+p_{1}p_{3}\right]+r_{3}^{2}\left(1+p_{3}^{2}\right)}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d0724dead4d7f7e4009fbd4ee02b45f9850efde7)
et ainsi des autres. On voit par là que

car l’expression de cette dernière quantité demeure la même, en changeant
en
et réciproquement ; ce qui est d’ailleurs évident.
XV.
Pour avoir maintenant la valeur de
il n’y aura qu’à changer, dans celle de
en
en
et effacer la quantité
(Article XI) ; on aura donc ainsi

on trouvera pareillement

et ainsi des autres.
XVI.
Nous avons supposé (Article X) que l’attraction de Jupiter sur les satellites était exactement en raison inverse du carré des distances ; c’est ce qui n’est rigoureusement vrai qu’en regardant Jupiter comme un globe de densité uniforme.
Or on sait par les observations et par la Théorie que cette Planète est considérablement aplatie ; de plus il peut se faire qu’elle ne soit pas partout de la même densité deux circonstances qui peuvent aussi influer sur le mouvement des satellites, et auxquelles il est bon par conséquent d’avoir égard ici. Pour cela nous supposerons 1o que la figure de Jupiter soit celle d’un sphéroïde elliptiqùe peu différent d’une sphère ; 2o que ce sphéroïde soit formé d’une infinité de couches toutes sphéroïdiques, et de densités différentes ; 3o que l’équateur de Jupiter soit dans le plan de l’orbite de cette Planète.
Cette dernière supposition n’est pas tout à fait exacte ; car on sait que l’équateur de Jupiter est incliné d’environ
degrés sur le plan de son orbite ; mais l’erreur qui en résulte est si petite qu’il serait superflu d’en tenir compte.
Cela posé soient
le demi-axe d’une couche quelconque,
son ellipticité et
sa densité ; on trouvera par les Théorèmes de la figure de la Terre de M. Clairaut (§§ XXVI et XLVI, seconde Partie) que l’attraction de Jupiter sur un satellitequelconque produit deux forces l’une, dirigée au centre de Jupiter, égale à
![{\displaystyle {\frac {2\varpi }{r^{2}\left(1+p^{2}\right)}}\left[\int \mathrm {DA^{2}} d\mathrm {A} +{\frac {2}{3}}\int \mathrm {D} d\mathrm {\left(A^{3}E\right)} \right]+{\frac {2\varpi \left(1-2p^{2}\right)}{5r^{4}\left(1+p^{2}\right)^{3}}}\int \mathrm {D} d\mathrm {\left(A^{3}E\right)} \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eb25d328b5a4dc5e1675df8f17feac00c048f946)
l’autre ; perpendiculaire à cette direction dans le plan d’un méridien, égale à

(
dénote ici la périphérie d’un cercle dont le rayon est égal à
). La partie
![{\displaystyle {\frac {2\varpi p}{r^{2}\left(1+p^{2}\right)}}\left[\int \mathrm {DA^{2}} d\mathrm {A} +{\frac {2}{3}}\int \mathrm {D} d\mathrm {\left(A^{3}E\right)} \right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/04122e0d8b4fb5ab01136151916b134a7db73427)
de la première de ces deux forces, étant réciproquement proportionnelle au carré de la distance, doit être comparée avec la force
(Article X) ; d’où l’on aura
![{\displaystyle \mathbb {Z} \!^{\upsilon }=2\varpi \left[\int \mathrm {DA^{2}} d\mathrm {A} +{\frac {2}{3}}\int \mathrm {D} d\mathrm {\left(A^{3}E\right)} \right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/84a03fd366bbc17e5c2514e460bb37aa689eb7a8)
L’autre partie de la même force

aussi bien que la force perpendiculaire

devront être regardées comme des forces perturbatrices, et par conséquent décomposées suivant les directions de
cette décomposi-
tion étant faite, on aura les deux forces suivantes

dans la direction de la force
et

dans la direction de la force
donc, si l’on suppose
![{\displaystyle \nu ={\frac {\int \mathrm {D} d\mathrm {\left(A^{5}E\right)} }{\mathrm {A} ^{2}\left[\int \mathrm {DA^{2}} d\mathrm {A} +{\cfrac {2}{3}}\int \mathrm {D} d\mathrm {\left(A^{3}E\right)} \right]}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8e8153777e7733157e9fcad98871abce497170cb)
les forces perturbatrices
et
qui résultent de l’aplatissement de Jupiter et de l’hétérogénéité de ses couches, seront, en général,

d’où l’on tire : par rapport au premier satellite,

par rapport au second satellite,

et ainsi des autres.
Il n’y aura donc qu’à ajouter ces valeurs à celles des Articles XII et XIII. Au reste, comme l’aplatissement de Jupiter n’est que d’environ
suivant les dernières observations, la quantité
sera fort petite, aussi bien que la quantité
de plus le rapport de
à
sera toujours exprimé par une fraction fort petite ; de sorte que les forces perturbatrices dont nous venons de parler seront nécessairement très-petites.
Si l’on suppose
constante, on aura

En général, quelle que soit
on aura, par les conditions de l’équilibre,

(ϐ étant le rapport de la force centrifuge à la pesanteur, sous l’équateur) ; donc

ϐ
à très-peu près.
XVII.
Il faut maintenant développer les expressions des forces perturbatrices
en employant les suppositions de l’article V. Pour cela nous remarquerons d’abord que nous pouvons négliger dans ce calcul tous les termes de l’ordre
parce que les quantités
sont déjà elles-mêmes de l’ordre
comme nous le verrons plus bas. Donc, mettant premièrement dans la valeur de
[Article XIV], au lieu de
au lieu de
et de même, au lieu de
et au lieu de
suivant ce que nous avons dit à l’Article IX on aura

d’où l’on tire, par les séries,
![{\displaystyle {\begin{aligned}{\frac {1}{\Delta (r_{1},r_{2})^{3}}}=&\left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {3}{2}}}\\&-3n\left[a_{1}^{2}x_{1}+a_{2}^{2}x_{2}-a_{1}a_{2}(x_{1}+x_{2})\cos(\varphi _{2}-\varphi _{1})\right]\\&\qquad \qquad \times \left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {5}{2}}}+\ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe39b5a214cc2437b2f8e695ec4049e379add88c)
On trouvera de même
![{\displaystyle {\begin{aligned}{\frac {1}{\Delta (r_{1},r_{3})^{3}}}=&\left[a_{1}^{2}-2a_{1}a_{3}\cos(\varphi _{3}-\varphi _{1})+a_{3}^{2}\right]^{-{\frac {3}{2}}}\\&-3n\left[a_{1}^{2}x_{1}+a_{3}^{2}x_{3}-a_{1}a_{3}(x_{1}+x_{3})\cos(\varphi _{3}-\varphi _{1})\right]\\&\qquad \qquad \times \left[a_{1}^{2}-2a_{1}a_{3}\cos(\varphi _{3}-\varphi _{1})+a_{3}^{2}\right]^{-{\frac {5}{2}}}+\ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f79cd4d1339abb9932817ad073010d0618f2b3c)
Et pareillement
![{\displaystyle {\begin{aligned}{\frac {1}{\Delta (r_{1},r_{4})^{3}}}=&\left[a_{1}^{2}-2a_{1}a_{4}\cos(\varphi _{4}-\varphi _{1})+a_{4}^{2}\right]^{-{\frac {3}{2}}}\\&-3n\left[a_{1}^{2}x_{1}+a_{4}^{2}x_{4}-a_{1}a_{4}(x_{1}+x_{4})\cos(\varphi _{4}-\varphi _{1})\right]\\&\qquad \qquad \times \left[a_{1}^{2}-2a_{1}a_{4}\cos(\varphi _{4}-\varphi _{1})+a_{4}^{2}\right]^{-{\frac {5}{2}}}+\ldots \,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0d645574bf48dbfd93e3de9b9a0249fb0beeddc7)
et ainsi des autres.
XVIII.
Mais il se présente ici une difficulté, par rapport aux quantités
![{\displaystyle \left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {3}{2}}},\quad \left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {5}{2}}},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/433b09868e573c3dadcd0921b09102cc72d47d13)
c’est de pouvoir les réduire à une forme rationnelle, condition absolument nécessaire pour l’intégration des équations des satellites.
Pour résoudre cette difliculté, on écrira d’abord les radicaux proposés ainsi
![{\displaystyle a_{2}^{-3}\left[1-{\frac {2a_{1}}{a_{2}}}\cos(\varphi _{2}-\varphi _{1})+{\frac {a_{1}^{2}}{a_{2}^{2}}}\right]^{-{\frac {3}{2}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/58be9fa790243f9322510b250575b8e27f106564)
![{\displaystyle a_{2}^{-5}\left[1-{\frac {2a_{1}}{a_{2}}}\cos(\varphi _{2}-\varphi _{1})+{\frac {a_{1}^{2}}{a_{2}^{2}}}\right]^{-{\frac {5}{2}}},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a46b72261b07178400be75f2df805a41bbdaf4a)
et la question se réduira à changer en une fonction rationnelle une quantité de cette forme

dans laquelle
est un nombre moindre que l’unité.
Pour y parvenir, je remarque que la quantité
est égale au produit de ces deux quantités

je les élève donc l’une et l’autre à la puissance
en écrivant au lieu du carré de
et ainsi de suite ;
j’ai
![{\displaystyle {\begin{aligned}&\left[1-q\left(\cos \theta \pm \sin \theta {\sqrt {-1}}\right)\right]^{-\lambda }=\\&\qquad \qquad \qquad \qquad 1+\lambda q\left(\cos \theta \pm \sin \theta {\sqrt {-1}}\right)\\&\qquad \qquad \qquad \qquad +{\frac {\lambda (\lambda +1)}{2}}q^{2}\left(\cos 2\theta \pm \sin 2\theta {\sqrt {-1}}\right)\\&\qquad \qquad \qquad \qquad +{\frac {\lambda (\lambda +1)(\lambda +2)}{2.3}}q^{3}\left(\cos 3\theta \pm \sin 3\theta {\sqrt {-1}}\right)+\ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19fdbc14712a8eb25b8a8b13f4e010ffa5b69e3e)
Soit, pour abréger,

on aura
![{\displaystyle {\begin{aligned}&\left[1-q\left(\cos \theta +\sin \theta {\sqrt {-1}}\right)\right]^{-\lambda }=\mathrm {M+N} {\sqrt {-1}},\\&\left[1-q\left(\cos \theta -\sin \theta {\sqrt {-1}}\right)\right]^{-\lambda }=\mathrm {M-N} {\sqrt {-1}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be56ef8a6d1781a80980d8a2142fb1b793913587)
donc

Or, si l’on fait les carrés des deux séries
et
qu’on ajoute ensemble les termes qui ont le même coefficient, et qu’on remarque que

et
étant des nombres quelconques, on trouvera

Et les coefficients
seront exprimés de la manière suivante

et ainsi de suite.
Au reste il ne sera nécessaire que de connaître les deux premiers coefficients
pour avoir tous les autres
car on trouvera par les formules de l’Article XXVI de la Pièce sur le mouvement de Saturne [Prix 1748][5]

et ainsi de suite.
XIX.
Tout consiste donc à déterminer les valeurs de
et
Or, dans la Théorie des satellites de Jupiter, la plus grande valeur de
est d’environ
comme on le verra plus bas ; donc
sera toujours moindre que
donc, si l’on fait
les suites
et
seront assez convergentes pour qu’on puisse se contenter d’un petit nombre de termes. Ces suites seront représentées, en général, par celles-ci

dont les coefficients numériques sont très-aisés à calculer.
Voici les logarithmes de ces coefficients pour les différentes puissances de
qui entrent dans les deux séries dont il s’agit ; les logarithmes qui répondent aux puissances paires de
sont ceux des coefficients des termes de la série
et les logarithmes qui répondent aux puissances impaires de
sont ceux des coefficients des termes de la série

Il ne s’agira donc plus que d’ajouter à chacun de ces logarithmes celui de la puissance correspondante de
et de chercher ensuite le nombre qui répond à chaque somme ; on aura ainsi les valeurs d’autant de termes des deux séries
et
qu’on voudra ; d’où l’on pourra tirer pour
et
des valeurs aussi approchées qu’on le croira nécessaire. Pour juger de la quantité de l’approximation, on remarquera que les différences des logarithmes de la Table précédente forment une progression décroissante ; d’où il suit que, si après avoir pris la somme d’un nombre quelconque de termes de la série
ou
on regarde le reste de la série comme une propression géométrique, l’erreur sera toujours moindre que la somme de cette progression. Au reste, dans le cas même où
sera la plus grande (ce cas est celui où
comme on le verra dans la suites), il suffira de prendre les dix premiers termes des séries
et
pour avoir les valeurs de ces coefficients en millièmes, c’est-à-dire aux dix-millièmes près, et en prenant encore trois ou quatre termes, on poussera l’exactitude jusqu’aux dix-millièmes et au delà.
XX.
Ayant ainsi les valeurs des coefficients
de la suite qui représente

on trouvera aisément ceux de la suite qui exprime

car, dénotant ces derniers par
il faudra que la série

étant multipliée par

devienne égale à la série

La multiplication faite, on trouvera, en comparant les deux premiers termes,

Or
est donné en
et
de la même manière que
est donné en
et
il suffira pour cela de mettre dans l’expression de
Article XVIII,
au lieu de
au lieu de
au lieu de
et
au lieu de
ce qui donnera

donc, si l’on substitue cette valeur de
on aura deux équations en
d’où l’on tirera

Connaissant
et
on connaîtra tous les suivants (Article XVIII).
XXI.
De ce qu’on vient de démontrer, il suit qu’on peut supposer
![{\displaystyle {\begin{aligned}&\left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {3}{2}}}\\&\quad =\Gamma (a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})\\&\qquad +\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\Gamma _{3}(a_{1},a_{2})\cos 3(\varphi _{2}-\varphi _{1})+\ldots ,\\\\&\left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {5}{2}}}\\&\quad =\Lambda (a_{1},a_{2})+\Lambda _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})\\&\qquad +\Lambda _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\Lambda _{3}(a_{1},a_{2})\cos 3(\varphi _{2}-\varphi _{1})+\ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf8f843a882258a8a40e46ca43f2d9791fc9f2e7)
J’entends par

des fonctions données de

dont on trouvera la valeur par les méthodes des Articles précédents.
Donc, si l’on fait ces substitutions dans la quantité
(Article XVII), et que l’on développe les produits des sinus et des cosinus, on trouvera
![{\displaystyle {\begin{aligned}&{\frac {1}{\Delta (r_{1},r_{2})^{3}}}\\&=\Gamma (a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})\\&\quad +\Gamma _{3}(a_{1},a_{2})\cos 3(\varphi _{2}-\varphi _{1})+\ldots \\&\quad -3nx_{1}[a_{1}^{2}\Lambda (a_{1},a_{2})-a_{1}a_{2}{\frac {\Lambda _{1}(a_{1},a_{2})}{2}}]\\&\quad -3nx_{1}[a_{1}^{2}\Lambda _{1}(a_{1},a_{2})-a_{1}a_{2}{\frac {2\Lambda (a_{1},a_{2})+\Lambda _{2}(a_{1},a_{2})}{2}}]\cos(\varphi _{2}-\varphi _{1})\\&\quad -3nx_{1}[a_{1}^{2}\Lambda _{2}(a_{1},a_{2})-a_{1}a_{2}{\frac {2\Lambda _{1}(a_{1},a_{2})+\Lambda _{3}(a_{1},a_{2})}{2}}]\cos 2(\varphi _{2}-\varphi _{1})\\&\quad -\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&\quad -3nx_{2}[a_{2}^{2}\Lambda (a_{1},a_{2})-a_{1}a_{2}{\frac {\Lambda _{1}(a_{1},a_{2})}{2}}]\\&\quad -3nx_{2}[a_{2}^{2}\Lambda _{1}(a_{1},a_{2})-a_{1}a_{2}{\frac {2\Lambda (a_{1},a_{2})+\Lambda _{2}(a_{1},a_{2})}{2}}]\cos(\varphi _{2}-\varphi _{1})\\&\quad -3nx_{2}[a_{2}^{2}\Lambda _{2}(a_{1},a_{2})-a_{1}a_{2}{\frac {2\Lambda _{1}(a_{1},a_{2})+\Lambda _{3}(a_{1},a_{2})}{2}}]\cos 2(\varphi _{2}-\varphi _{1})\\&\quad -\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8996526d7696b6e30f2ca85431da861ebc3d3629)
XXII.
Soit fait, pour plus de simplicité,

Soit aussi

On aura la quantité
exprimée de la manière suivante
![{\displaystyle {\begin{aligned}&{\frac {1}{\Delta (r_{1},r_{2})^{3}}}\\&=\Gamma (a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \\&\quad +3nx_{1}\left[\Pi (a_{1},a_{2})+\Pi _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Pi _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&\quad +3nx_{2}\left[\Psi (a_{1},a_{2})+\Psi _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Psi _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/84e805334af002788ccd1230439eff1656ead58b)
On trouvera de même, en changeant simplement
en
en
en
et
en
![{\displaystyle {\begin{aligned}&{\frac {1}{\Delta (r_{1},r_{3})^{3}}}\\&=\Gamma (a_{1},a_{3})+\Gamma _{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\Gamma _{2}(a_{1},a_{3})\cos 2(\varphi _{3}-\varphi _{1})+\ldots \\&\quad +3nx_{1}\left[\Pi (a_{1},a_{3})+\Pi _{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\Pi _{2}(a_{1},a_{3})\cos 2(\varphi _{3}-\varphi _{1})+\ldots \right]\\&\quad +3nx_{3}\left[\Psi (a_{1},a_{3})+\Psi _{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\Psi _{2}(a_{1},a_{3})\cos 2(\varphi _{3}-\varphi _{1})+\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11327b6e941a9bd028286d1fbcfbaf3399dcad06)
et pareillement
![{\displaystyle {\begin{aligned}&{\frac {1}{\Delta (r_{1},r_{4})^{3}}}\\&=\Gamma (a_{1},a_{4})+\Gamma _{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\Gamma _{2}(a_{1},a_{4})\cos 2(\varphi _{4}-\varphi _{1})+\ldots \\&\quad +3nx_{1}\left[\Pi (a_{1},a_{4})+\Pi _{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\Pi _{2}(a_{1},a_{4})\cos 2(\varphi _{4}-\varphi _{1})+\ldots \right]\\&\quad +3nx_{4}\left[\Psi (a_{1},a_{4})+\Psi _{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\Psi _{2}(a_{1},a_{4})\cos 2(\varphi _{4}-\varphi _{1})+\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/71f63544b9ff572d4820f4ab52df31b722a41dc6)
XXIII.
Cela posé, on aura
![{\displaystyle {\begin{aligned}&{\frac {r_{1}}{\Delta (r_{1},r_{2})^{3}}}\\&=a_{1}\left[\Gamma (a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&\quad +nx_{1}a_{1}\left[3\Pi (a_{1},a_{2})+\Gamma (a_{1},a_{2})+\left[3\Pi _{1}(a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})\right]\cos(\varphi _{2}-\varphi _{1})\right.\\&\qquad \qquad \qquad \qquad \left.+\left[3\Pi _{2}(a_{1},a_{2})+\Gamma _{2}(a_{1},a_{2})\right]\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&\quad +3nx_{2}a_{1}\left[\Psi (a_{1},a_{2})+\Psi _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Psi _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f0542b39214d68e2b4d6c95ff3675a8c781286ee)
On trouvera de la même manière
![{\displaystyle {\begin{aligned}&{\frac {r_{2}}{\Delta (r_{1},r_{2})^{3}}}\\&=a_{2}\left[\Gamma (a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&\quad +3nx_{1}a_{2}\left[\Pi (a_{1},a_{2})+\Pi _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Pi _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&\quad +nx_{2}a_{2}\left[3\Psi (a_{1},a_{2})+\Gamma (a_{1},a_{2})+\left[3\Psi _{1}(a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})\right]\cos(\varphi _{2}-\varphi _{1})\right.\\&\qquad \qquad \qquad \qquad \left.+\left[3\Psi _{2}(a_{1},a_{2})+\Gamma _{2}(a_{1},a_{2})\right]\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7bf75e4d8cd918215c8557cfdce037e834bbe862)
On aura ensuite

Donc

![{\displaystyle {\begin{aligned}&={\frac {1}{a_{2}^{2}}}-a_{2}\Gamma (a_{1},a_{2})-a_{2}\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})-a_{2}\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})-\ldots \\&\quad -3nx_{1}a_{2}\left[\Pi (a_{1},a_{2})+\Pi _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\Pi _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&\quad -nx_{2}\left[{\frac {2}{a_{2}^{2}}}+3a_{2}\Psi (a_{1},a_{2})+a_{2}\Gamma (a_{1},a_{2})\right.\\&\qquad \qquad \qquad \qquad +\left[3a_{2}\Psi _{1}(a_{1},a_{2})+a_{2}\Gamma _{1}(a_{1},a_{2})\right]\cos(\varphi _{2}-\varphi _{1})\\&\qquad \qquad \qquad \qquad +{\biggl .}\left[3a_{2}\Psi _{2}(a_{1},a_{2})+a_{2}\Gamma _{2}(a_{1},a_{2})\right]\cos 2(\varphi _{2}-\varphi _{1})+\ldots {\biggr ]}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1444790b4950b6aff23e7bbb5fec897f57809f39)
Cette quantité étant multipliée par

on aura

![{\displaystyle {\begin{aligned}&=-{\frac {a_{2}\Gamma _{1}(a_{1},a_{2})}{2}}-\left[{\frac {a_{2}\Gamma _{2}(a_{1},a_{2})+2a_{2}\Gamma (a_{1},a_{2})}{2}}-{\frac {1}{a_{2}^{2}}}\right]\cos(\varphi _{2}-\varphi _{1})\\&\quad -{\frac {a_{2}\Gamma _{3}(a_{1},a_{2})+a_{2}\Gamma _{1}(a_{1},a_{2})}{2}}\cos(\varphi _{2}-\varphi _{1})-\ldots \\&\quad -nx_{1}\left[3a_{2}{\frac {\Pi _{1}(a_{1},a_{2})}{2}}+3a_{2}{\frac {\Pi _{2}(a_{1},a_{2})+2\Pi (a_{1},a_{2})}{2}}\cos(\varphi _{2}-\varphi _{1})\right.\\&\qquad \qquad \qquad \qquad \qquad \quad +\left.3a_{2}{\frac {\Pi _{3}(a_{1},a_{2})+\Pi _{1}(a_{1},a_{2})}{2}}\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&\quad -nx_{2}\left[{\frac {3a_{2}\Psi _{1}(a_{1},a_{2})+a_{2}\Gamma _{1}(a_{1},a_{2})}{2}}\right.\\&\qquad \qquad \quad +\left(3a_{2}{\frac {\Psi _{2}(a_{1},a_{2})+2\Psi (a_{1},a_{2})}{2}}+a_{2}{\frac {\Gamma _{2}(a_{1},a_{2})+2\Gamma (a_{1},a_{2})}{2}}+{\frac {2}{a_{2}^{2}}}\right)\\&\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \times \cos(\varphi _{2}-\varphi _{1})\\\\&\qquad \qquad \quad +\left(3a_{2}{\frac {\Psi _{3}(a_{1},a_{2})+\Psi _{1}(a_{1},a_{2})}{2}}+a_{2}{\frac {\Gamma _{3}(a_{1},a_{2})+\Gamma _{1}(a_{1},a_{2})}{2}}\right)\\&\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad {\biggl .}\times \cos 2(\varphi _{2}-\varphi _{1})+\ldots {\biggr ]}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3afbba2db36695c31a30e71c6b07842af6fede20)
Enfin, multipliant la même quantité pa
on aura

![{\displaystyle {\begin{aligned}&=-\left[a_{2}{\frac {\Gamma _{2}(a_{1},a_{2})-2\Gamma (a_{1},a_{2})}{2}}+{\frac {1}{a_{2}^{2}}}\right]\sin(\varphi _{2}-\varphi _{1})\\&\quad -a_{2}{\frac {\Gamma _{3}(a_{1},a_{2})-\Gamma _{1}(a_{1},a_{2})}{2}}\sin 2(\varphi _{2}-\varphi _{1})-\ldots \\&\quad -nx_{1}\left[3a_{2}{\frac {\Pi _{2}(a_{1},a_{2})-2\Pi (a_{1},a_{2})}{2}}\sin(\varphi _{2}-\varphi _{1})\right.\\&\qquad \qquad \qquad \qquad +\left.3a_{2}{\frac {\Pi _{3}(a_{1},a_{2})-\Pi _{1}(a_{1},a_{2})}{2}}\sin 2(\varphi _{2}-\varphi _{1})+\ldots \right]\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e47a06c502c57dba81a7e3c91b76bbfb3ac22de5)
![{\displaystyle {\begin{aligned}&\quad -nx_{2}\left[\left(3a_{2}{\frac {\Psi _{2}(a_{1},a_{2})-2\Psi (a_{1},a_{2})}{2}}+a_{2}{\frac {\Gamma _{2}(a_{1},a_{2})-2\Gamma (a_{1},a_{2})}{2}}-{\frac {2}{a_{2}^{2}}}\right)\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \times \sin(\varphi _{2}-\varphi _{1})\\\\&\qquad \qquad \quad +\left(3a_{2}{\frac {\Psi _{3}(a_{2},a_{1})-\Psi _{1}(a_{1},a_{2})}{2}}+a_{2}{\frac {\Gamma _{3}(a_{2},a_{1})-\Gamma _{1}(a_{1},a_{2})}{2}}\right)\\&\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad {\Biggl .}\times \sin(\varphi _{2}-\varphi _{1}){\Biggr ]}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/da43af4c47614802ffa2f7273b7a8624ead6f97b)
XXIV.
Soit maintenant


On aura
![{\displaystyle {\begin{aligned}&{\mathfrak {S}}_{2}\left[{\frac {r_{1}-r_{2}\cos(\varphi _{2}-\varphi _{1})}{\Delta (r_{1},r_{2})^{3}}}+{\frac {\cos(\varphi _{2}-\varphi _{1})}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right]=\\&-{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}\left[{\breve {\Gamma }}(a_{1},a_{2})+{\breve {\Gamma }}_{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+{\breve {\Gamma }}_{2}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}x_{1}\left[{\breve {\Pi }}(a_{1},a_{2})+{\breve {\Pi }}_{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+{\breve {\Pi }}_{2}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}x_{2}\left[{\breve {\Psi }}(a_{1},a_{2})+{\breve {\Psi }}_{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+{\breve {\Psi }}_{2}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb310f33b027c75730d73c3015388458594ab872)
C’est la partie de la force
qui résulte de l’action du satellite
(Article X).
XXV.
Soit de plus

On aura
![{\displaystyle {\begin{aligned}&{\mathfrak {S}}_{2}\left[{\frac {r_{2}\sin(\varphi _{2}-\varphi _{1})}{\Delta (r_{1},r_{2})^{3}}}-{\frac {\sin(\varphi _{2}-\varphi _{1})}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right]=\\&-{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}\left[{\widehat {\Gamma }}_{1}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+{\widehat {\Gamma }}_{2}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}x_{1}\left[{\widehat {\Pi }}_{1}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+{\widehat {\Pi }}_{2}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}x_{2}\left[{\widehat {\Psi }}_{1}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+{\widehat {\Psi }}_{2}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/422d0e212cf53ef343661767ae65a687b9fc94f4)
C’est la valeur de la force
en tant qu’elle vient de l’action du satellite
XXVI.
Enfin on trouvera
![{\displaystyle {\begin{aligned}{\mathfrak {S}}_{2}&\left[{\frac {r_{1}p_{1}-r_{2}p_{2}}{\Delta (r_{1},r_{2})^{3}}}+{\frac {p_{2}}{r_{2}^{2}\left(1+p_{2}^{2}\right)^{\frac {3}{2}}}}\right]\\=&n{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}z_{1}\left[a_{1}^{3}\Gamma (a_{1},a_{2})+a_{1}^{3}\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})\right.\\&\quad \qquad \qquad \qquad \qquad \qquad \qquad \left.+a_{1}^{3}\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{2}}{a_{1}^{2}}}z_{2}\left[a_{1}^{2}a_{2}\Gamma (a_{1},a_{2})-{\frac {a_{1}^{2}}{a_{2}^{2}}}+a_{1}^{2}a_{2}\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad \qquad {\biggl .}+a_{1}^{2}a_{2}\Gamma _{2}(a_{1},a_{2})\cos 2(\varphi _{2}-\varphi _{1}){\biggr ]}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d9b14b131fb93e451b9b36a1242fb884558e015a)
C’est la partie de force
qui vient de l’action du même satellite
On changera maintenant dans les expressions précédentes les quantités
en
et en
successivement, et l’on aura les valeurs de
dues à l’action des satellites
Il ne restera donc plus qu’à chercher les valeurs de ces mêmes forces, en tant qu’elles viennent de l’action du Soleil.
Pour cela nous remarquerons d’abord que le rayon
de l’orbite de Jupiter est considérablement plus grand que le rayon
de l’orbite d’un satellite quelconque ; d’où il suit que la valeur de
qui est exprimée généralement (Article XV) par

se réduira en une suite très-convergente, dont il suffira de prendre les premiers termes ; on aura donc

donc
![{\displaystyle {\begin{aligned}-{\frac {\rho _{1}\cos(\psi -\varphi )}{\delta ^{3}}}+{\frac {\cos(\psi -\varphi )}{\rho _{1}^{2}}}=&-{\frac {3r}{2\rho _{1}^{3}}}\left[1+\cos 2(\psi -\varphi )\right],\\{\frac {\rho _{1}\sin(\psi -\varphi )}{\delta ^{3}}}-{\frac {\sin(\psi -\varphi )}{\rho _{1}^{2}}}=&{\frac {3r}{2\rho _{1}^{3}}}\sin 2(\psi -\varphi ),\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2738d65cd60dd0af1f34538b1f895634c0b1ce53)
et

XXVII.
Soient à présent
la valeur moyenne de
et
la valeur moyenne de
c’est-à-dire la vitesse angulaire moyenne de Jupiter autour du Soleil.
On supposera, à l’imitation de ce que nous avons fait (Article IV),

Dans ces formules,
représente l’équation de la distance de Jupiter au Soleil, et
l’équation du centre de Jupiter ; lesquelles sont connues par la théorie de cette Planète. En effet, en n’ayant égard qu’aux équations elliptiques, et supposant que
soit l’excentricité et
l’anomalie moyenne, on a à très-peu près

On aura donc

donc enfin
![{\displaystyle {\begin{aligned}\circledast &\left[r-{\frac {\rho _{1}\cos(\psi -\varphi _{1})}{\delta _{1}^{3}}}+{\frac {\cos(\psi -\varphi _{1})}{\rho _{1}^{2}}}\right]\\&=-{\frac {a_{1}\circledast }{2\alpha ^{3}}}\left[1+3\cos 2(\psi -\varphi _{1})\right]-n{\frac {a_{1}\circledast }{2\alpha ^{3}}}(x_{1}-3\xi )\left[1+3\cos 2(\psi -\varphi _{1})\right],\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e8caaab49b02a128d1c47afc24dc3105890cb941)
![{\displaystyle {\begin{aligned}\circledast &\left[{\frac {\rho _{1}\sin(\psi -\varphi _{1})}{\delta _{1}^{3}}}+{\frac {\sin(\psi -\varphi _{1})}{\rho _{1}^{2}}}\right]\\&={\frac {3a_{1}\circledast }{2\alpha ^{3}}}\sin 2(\psi -\varphi _{1})+n{\frac {3a_{1}\circledast }{2\alpha ^{3}}}(x_{1}-3\xi )\sin 2(\psi -\varphi _{1}),\\\circledast &{\frac {r_{1}p_{1}}{\delta _{1}^{3}}}=n{\frac {a_{1}\circledast }{\alpha _{1}^{3}}}z_{1}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd6afb1e7bff7ec82ba76b2cc5f19084bece6c47)
Ce sont les valeurs des forces
qui viennent de l’action du Soleil (Article XI).
XXVIII.
En joignant ensemble toutes ces différentes valeurs, on aura les valeurs complètes des forces
exprimées de la manière suivante
![{\displaystyle {\begin{aligned}\mathrm {R} _{1}=&-{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\breve {\Gamma }}(a_{1},a_{2})+{\breve {\Gamma }}_{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\breve {\Gamma }}(a_{1},a_{3})+{\breve {\Gamma }}_{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\ldots \right]\\&-{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\breve {\Gamma }}(a_{1},a_{4})+{\breve {\Gamma }}_{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\ldots \right]\\&-{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\frac {1}{2}}+{\frac {3}{2}}\cos 2(\psi -\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\left[{\breve {\Pi }}(a_{1},a_{2})+{\breve {\Pi }}_{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\left[{\breve {\Pi }}(a_{1},a_{3})+{\breve {\Pi }}_{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\left[{\breve {\Pi }}(a_{1},a_{4})+{\breve {\Pi }}_{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\ldots \right]\\&-n{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\left[{\frac {1}{2}}+{\frac {3}{2}}\cos 2(\psi -\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{2}\left[{\breve {\Psi }}(a_{1},a_{2})+{\breve {\Psi }}_{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{3}\left[{\breve {\Psi }}(a_{1},a_{3})+{\breve {\Psi }}_{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\ldots \right]\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b68c1e2fdccf7cbf23ffb7bce02ea454c4b3b249)
![{\displaystyle {\begin{aligned}&-n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\widehat {\Psi }}(a_{1},a_{4})+{\widehat {\Psi }}_{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\ldots \right]\\&-n{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\mathbb {Z} \!^{\upsilon }}{\circledast }}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\xi \left[{\frac {3}{2}}+{\frac {9}{2}}\cos 2(\psi -\varphi _{1})+\ldots \right],\\\mathrm {Q} _{1}=&-{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\widehat {\Gamma }}_{1}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\widehat {\Gamma }}_{1}(a_{1},a_{3})\sin(\varphi _{3}-\varphi _{1})+\ldots \right]\\&-{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\left[{\widehat {\Gamma }}_{1}(a_{1},a_{4})\sin(\varphi _{4}-\varphi _{1})+\ldots \right]\\&-{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{\alpha _{1}^{2}}}\times {\frac {3}{2}}\sin 2(\psi -\varphi _{1})\\\\&-n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\left[{\widehat {\Pi }}_{1}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\left[{\widehat {\Pi }}_{1}(a_{1},a_{3})\sin(\varphi _{3}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\left[{\widehat {\Pi }}_{1}(a_{1},a_{4})\sin(\varphi _{4}-\varphi _{1})+\ldots \right]\\&-n{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{1}\times {\frac {3}{2}}\sin 2(\psi -\varphi _{1})+\ldots \\&-n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{2}\left[{\widehat {\Psi }}_{1}(a_{1},a_{2})\sin(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{3}\left[{\widehat {\Psi }}_{1}(a_{1},a_{3})\sin(\varphi _{3}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}x_{4}\left[{\widehat {\Psi }}_{1}(a_{1},a_{4})\sin(\varphi _{4}-\varphi _{1})+\ldots \right]\\&-n{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}\xi \times {\frac {9}{2}}\sin 2(\psi -\varphi _{1}),\\\\\mathrm {P} _{1}=&n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}z_{1}\left[a_{1}^{3}\Gamma (a_{1},a_{2})+a_{1}^{3}\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\ldots \right]\\&+n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}z_{1}\left[a_{1}^{3}\Gamma (a_{1},a_{2})+a_{1}^{3}\Gamma _{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\ldots \right]\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6bda30b9f790f62663c3898722f635e038f55bb4)
![{\displaystyle {\begin{aligned}&+n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}z_{1}\left[a_{1}^{3}\Gamma (a_{1},a_{4})+a_{1}^{3}\Gamma _{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\ldots \right]\\&+n{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}z_{1}\\&-n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}z_{2}\left[a_{1}^{2}a_{2}\Gamma (a_{1},a_{2})-{\frac {a_{1}^{2}}{a_{2}^{2}}}+a_{1}^{2}a_{2}\Gamma _{1}(a_{1},a_{2})\cos(\varphi _{2}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}z_{3}\left[a_{1}^{2}a_{3}\Gamma (a_{1},a_{3})-{\frac {a_{1}^{2}}{a_{3}^{2}}}+a_{1}^{2}a_{3}\Gamma _{1}(a_{1},a_{3})\cos(\varphi _{3}-\varphi _{1})+\ldots \right]\\&-n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}z_{4}\left[a_{1}^{2}a_{4}\Gamma (a_{1},a_{4})-{\frac {a_{1}^{2}}{a_{4}^{2}}}+a_{1}^{2}a_{4}\Gamma _{1}(a_{1},a_{4})\cos(\varphi _{4}-\varphi _{1})+\ldots \right]\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c5cd5bda41f97b832d41686f3a7406bff63f2813)
Il ne reste plus qu’à substituer pour
et
leurs valeurs
et
ce qui est très-facile, car il n’y aura qu’à changera,
en
et ajouter ensuite aux expressions de
et
les quantités suivantes
![{\displaystyle {\begin{aligned}&n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(y_{2}-y_{1})\left[{\breve {\Gamma }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{2})+\ldots \right]\\+&n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(y_{3}-y_{1})\left[{\breve {\Gamma }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{3})+\ldots \right]\\+&n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(y_{4}-y_{1})\left[{\breve {\Gamma }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{4})+\ldots \right]\\+&n{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(\mathrm {J} -y_{1})\times 3\sin 2(m-\mu _{1})t,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be0fd9f402344dfb0d623e741b7517e90a9f5141)
et
![{\displaystyle {\begin{aligned}-&n{\frac {{\mathfrak {S}}_{2}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(y_{2}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{2})+\ldots \right]\\-&n{\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(y_{3}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{3})+\ldots \right]\\-&n{\frac {{\mathfrak {S}}_{4}}{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(y_{4}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{4})+\ldots \right]\\+&n{\frac {a_{1}^{3}}{\alpha ^{3}}}{\frac {\circledast }{\mathbb {Z} \!^{\upsilon }}}{\frac {\mathbb {Z} \!^{\upsilon }}{a_{1}^{2}}}(\mathrm {J} -y_{1})\times 3\cos 2(m-\mu _{1})t.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/74dcaacd05a578f272ecd605f11ab52e2adb76c3)
Pour trouver les valeurs de
il ne faudra qu’ajouter un trait aux lettres qui n’en ont qu’un, et en ôter un à celles qui en ont deux ; et ainsi des autres quantités
[Article XIII][6].
À l’égard des forces perturbatrices qui viennent de la non-sphéricité de Jupiter, on trouvera, en négligeant dans les formules de l’Article XVI ce qu’on y doit négliger, qu’il faut ajouter aux valeurs de
et
les quantités suivantes

et de même aux valeurs de
les quantités

et ainsi de suite.
CHAPITRE III.
CALCUL DES PERTURBATIONS DES SATELLITES DE JUPITER.
XXIX.
Nous nous contenterons ici de chercher les formules qui déterminent le mouvement du premier satellite, parce que les autres s’en déduiront aisément par les remarques des Articles IX et XIII.
Pour appliquer au mouvement du premier satellite les équations générales de l’Article VIII il est visible qu’il ne faut que marquer les lettres d’un trait[6], et substituer ensuite pour
leurs valeurs tirées de celles de
(Article précédent). Mais, avant que de faire cette substitution, nous remarquerons que les équations (Articles V et VI)

ne peuvent subsister dans l’hypothèse de

très-petit, à moins que les quantités constantes

et les quantités variables

et

ne soient chacune très-petites de l’ordre

.
Or, en examinant les valeurs de
et
(Article précédent), il est facile de voir qu’elles ne sauraient être supposées très-petites qu’en regardant comme telles les quantités constantes
et
Soit donc, en général,

Soit, de plus,

Nous aurons, à cause de
(Article V),

c’est-à-dire, à cause de
(Article X),

À l’égard de la quantité
elle sera déterminée par l’équation (Article VII)

laquelle se réduit à

où l’on remarquera que
est déjà toute multipliée par
(Article XXVIII).
XXX.
Appliquons maintenant ces formules au premier satellite. Nous aurons d’abord en substituant la valeur de
(Article XXIX)
![{\displaystyle {\begin{aligned}&\mathrm {X} _{1}=g_{1}\\&-{\frac {\chi _{2}}{1+n\chi _{1}}}\left[{\breve {\Gamma }}(a_{1},a_{2})+{\breve {\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+{\breve {\Gamma }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&-{\frac {\chi _{3}}{1+n\chi _{1}}}\left[{\breve {\Gamma }}(a_{1},a_{3})+{\breve {\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+{\breve {\Gamma }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&-{\frac {\chi _{4}}{1+n\chi _{1}}}\left[{\breve {\Gamma }}(a_{1},a_{4})+{\breve {\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+{\breve {\Gamma }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&-{\frac {\mathrm {K} _{1}}{1+n\chi _{1}}}\left[{\frac {1}{2}}+{\frac {3}{2}}\cos 2(m-\mu _{1})t\right]\\&-{\frac {n\chi _{2}}{1+n\chi _{1}}}x_{1}\left[{\breve {\Pi }}(a_{1},a_{2})+{\breve {\Pi }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+{\breve {\Pi }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&-{\frac {n\chi _{3}}{1+n\chi _{1}}}x_{1}\left[{\breve {\Pi }}(a_{1},a_{3})+{\breve {\Pi }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+{\breve {\Pi }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&-{\frac {n\chi _{4}}{1+n\chi _{1}}}x_{1}\left[{\breve {\Pi }}(a_{1},a_{4})+{\breve {\Pi }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+{\breve {\Pi }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&-{\frac {n\mathrm {K} _{1}}{1+n\chi _{1}}}x_{1}\left[{\frac {1}{2}}+{\frac {3}{2}}\cos 2(m-\mu _{1})\right]\\&-{\frac {n\chi _{2}}{1+n\chi _{1}}}x_{2}\left[{\breve {\Psi }}(a_{1},a_{2})+{\breve {\Psi }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+{\breve {\Psi }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&-{\frac {n\chi _{3}}{1+n\chi _{1}}}x_{3}\left[{\breve {\Psi }}(a_{1},a_{3})+{\breve {\Psi }}_{1}(a_{1},a_{3})(\cos \mu _{3}-\mu _{1})t+{\breve {\Psi }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&-{\frac {n\chi _{4}}{1+n\chi _{1}}}x_{4}\left[{\breve {\Psi }}(a_{1},a_{4})+{\breve {\Psi }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+{\breve {\Psi }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&+{\frac {n\mathrm {K} _{1}}{1+n\chi _{1}}}\xi \left[{\frac {1}{2}}+{\frac {3}{2}}\cos 2(m-\mu _{1})t\right]\\&+{\frac {n\chi _{2}}{1+n\chi _{1}}}(y_{2}-y_{1})\left[{\breve {\Gamma }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&+{\frac {n\chi _{3}}{1+n\chi _{1}}}(y_{3}-y_{1})\left[{\breve {\Gamma }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&+{\frac {n\chi _{4}}{1+n\chi _{1}}}(y_{4}-y_{1})\left[{\breve {\Gamma }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&+{\frac {n\mathrm {K} _{1}}{1+n\chi _{1}}}(\mathrm {J} -y)3\sin 2(m-\mu _{1})t\\&+{\frac {\varkappa _{1}}{1+n\chi _{1}}}\left({\frac {1}{5}}-{\frac {4}{5}}nx_{1}\right).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e7d1c775d458ebcf727a68d09c89155ce98434ae)
XXXI.
Multipliant la valeur de
par
et faisant, pour abréger,

on aura, après l’intégration et la substitution,
![{\displaystyle {\begin{aligned}&\mathrm {Y} _{1}=\mathrm {H} _{1}\\&+{\frac {\chi _{2}}{1+n\chi _{1}}}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{2})}{\mu _{2}-\mu _{1}}}\cos(\mu _{2}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{2})}{2(\mu _{2}-\mu _{1})}}\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&+{\frac {\chi _{3}}{1+n\chi _{1}}}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{3})}{\mu _{3}-\mu _{1}}}\cos(\mu _{3}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{3})}{2(\mu _{3}-\mu _{1})}}\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&+{\frac {\chi _{4}}{1+n\chi _{1}}}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{4})}{\mu _{4}-\mu _{1}}}\cos(\mu _{4}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{4})}{2(\mu _{4}-\mu _{1})}}\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&-{\frac {\mathrm {K} _{1}}{1+n\chi _{1}}}{\frac {3}{4(m-\mu _{1})}}\cos 2(m-\mu _{1})t\\&-{\frac {n\chi _{2}}{1+n\chi _{1}}}\int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\&-{\frac {n\chi _{3}}{1+n\chi _{1}}}\int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\&-{\frac {n\chi _{4}}{1+n\chi _{1}}}\int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\&+{\frac {n\mathrm {K} _{1}}{1+n\chi _{1}}}\int x_{1}\times 3\sin 2(m-\mu _{1})tdt\\&-{\frac {n\chi _{2}}{1+n\chi _{1}}}\int x_{2}\left[{\widehat {\Psi }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\&-{\frac {n\chi _{3}}{1+n\chi _{1}}}\int x_{3}\left[{\widehat {\Psi }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\&-{\frac {n\chi _{4}}{1+n\chi _{1}}}\int x_{4}\left[{\widehat {\Psi }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\&-{\frac {n\mathrm {K} _{1}}{1+n\chi _{1}}}\int \xi \times {\frac {9}{2}}\sin 2(m-\mu _{1})tdt\\&-{\frac {n\chi _{2}}{1+n\chi _{1}}}\int (y_{2}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f8ad2f70c41d9df0b1d295878df4c9b57d82ba02)
![{\displaystyle {\begin{aligned}&-{\frac {n\chi _{3}}{1+n\chi _{1}}}\int (y_{3}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\&-{\frac {n\chi _{4}}{1+n\chi _{1}}}\int (y_{4}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\&+{\frac {n\mathrm {K} _{1}}{1+n\chi _{1}}}\int (\mathrm {J} -y_{1})\times 3\cos 2(m-\mu _{1})tdt.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/13876c51e88d555b8014c462dbf94a124a643976)
XXXII.
Enfin, si l’on fait

et ainsi de suite, on trouvera
![{\displaystyle {\begin{aligned}&\mathrm {Z} _{1}=\\&\quad {\frac {n\chi _{2}}{1+n\chi _{1}}}z_{1}\left[{\overset {\backsim }{\Gamma }}(a_{1},a_{2})+{\overset {\backsim }{\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+{\overset {\backsim }{\Gamma }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&+{\frac {n\chi _{3}}{1+n\chi _{1}}}z_{1}\left[{\overset {\backsim }{\Gamma }}(a_{1},a_{3})+{\overset {\backsim }{\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+{\overset {\backsim }{\Gamma }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&+{\frac {n\chi _{4}}{1+n\chi _{1}}}z_{1}\left[{\overset {\backsim }{\Gamma }}(a_{1},a_{4})+{\overset {\backsim }{\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+{\overset {\backsim }{\Gamma }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&+{\frac {n\mathrm {K} _{1}}{1+n\chi _{1}}}{\frac {2}{5}}z_{1}\\&-{\frac {n\chi _{2}}{1+n\chi _{1}}}z_{2}\left[{\overset {\circ }{\Gamma }}(a_{1},a_{2})+{\overset {\circ }{\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+{\overset {\circ }{\Gamma }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&-{\frac {n\chi _{3}}{1+n\chi _{1}}}z_{3}\left[{\overset {\circ }{\Gamma }}(a_{1},a_{3})+{\overset {\circ }{\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+{\overset {\circ }{\Gamma }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&-{\frac {n\chi _{4}}{1+n\chi _{1}}}z_{4}\left[{\overset {\circ }{\Gamma }}(a_{1},a_{4})+{\overset {\circ }{\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+{\overset {\circ }{\Gamma }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right].\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc6e5086417c592cd8ebb4d4aa5ce15e0cbf2c67)
XXXIII.
Et le mouvement du satellite
sera déterminé par les équations suivantes (Article VIII)
1
o

2
o
3
o
On se souviendra que les quantités
et
ne doivent renfermer aucun terme constant, suivant la remarque de l’Article IV.
XXXIV.
Il ne s’agit donc plus que d’intégrer les équations que nous venons de donner ; pour cela, on commencera par rejeter tous les termes affectés de
et l’on cherchera par l’intégration les valeurs de
on substituera ensuite ces valeurs dans les termes qu’on avait négligés, et l’on en tirera de nouvelles valeurs de
plus approchées que les premières. On opérerait ainsi de suite, si nous avions eu égard aux termes affectés de
Par ce moyen, l’intégration de la première et de la troisième équation de l’Article précédent se réduira à celle d’une équation de cette forme

étant une fonction composée de sinus et de cosinus d’angles multiples de
or l’intégrale de cette équation est, comme l’on sait,

de sorte qu’en supposant

on aura

et
sont les valeurs de
et de
lorsque 
On voit de là que, pour avoir la valeur de
il n’y a qu’à diviser chacun des sinus et des cosinus qui entrent dans
par
étant le coefficient de
et y ajouter encore deux autres termes, qui renferment
et
avec des coefficients arbitraires.
Il ne peut y avoir de difficulté que dans le cas où
serait égal à
car alors le diviseur
sera nul, et les termes

aussi bien que les termes

deviendraient
ce qui ne fait rien connaitre.
Pour résoudre cette difficulté, on supposera que
ne soit pas tout à fait égale à
mais qu’elle en diffère d’une quantité infiniment petite ; et l’on trouvera que les deux premiers termes se réduisent à

et les deux autres à

d’où l’on voit que la valeur de
contiendra des termes multipliés par l’angle 
Au reste, si dans la quantité
il se trouve des termes de cette forme
ou
étant égal à
il est visible que ces termes augmenteront beaucoup par l’intégration, puisqu’ils se trouveront divisés par la quantité très-petite
Donc, si ces sortes de termes ont des coefficients finis dans l’équation différentielle, ils deviendront comme infinis dans l’intégrale ; et, s’ils n’ont que des coefficients très-petits de l’ordre
dans la différentielle, ils deviendront finis dans l’intégrale, et appartiendront à la première valeur de
.
§ I. — Premières formules du mouvement des satellites
de Jupiter autour de cette Planète.
XXXV.
Si l’on substitue dans les trois équations de l’Article XXXIII les valeurs
et
qu’on rejette d’abord tous les termes affectés de
et que l’on fasse, pour abréger,


on trouvera les trois équations suivantes
![{\displaystyle \left.{\begin{aligned}1^{\circ }&\ \ {\frac {d^{2}x_{1}}{dt^{2}}}+\mathrm {M} _{1}^{2}x_{1}+f_{1}\mathrm {L} _{1}\\&-\chi _{2}f_{1}\left[{\breve {\Gamma }}_{1}(a_{1},a_{2})+{\frac {2\mu _{1}}{\mu _{2}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\right]\cos(\mu _{2}-\mu _{1})t\\&-\chi _{2}f_{1}\left[{\breve {\Gamma }}_{2}(a_{1},a_{2})+{\frac {2\mu _{1}}{2(\mu _{2}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{2})\right]\cos 2(\mu _{2}-\mu _{1})t+\ldots \\&-\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&-\chi _{3}f_{1}\left[{\breve {\Gamma }}_{1}(a_{1},a_{3})+{\frac {2\mu _{1}}{\mu _{3}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\right]\cos(\mu _{3}-\mu _{1})t\\&-\chi _{3}f_{1}\left[{\breve {\Gamma }}_{2}(a_{1},a_{3})+{\frac {2\mu _{1}}{2(\mu _{3}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{3})\right]\cos 2(\mu _{3}-\mu _{1})t+\ldots \\&-\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&-\chi _{4}f_{1}\left[{\breve {\Gamma }}_{1}(a_{1},a_{4})+{\frac {2\mu _{1}}{\mu _{4}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{4})\right]\cos(\mu _{4}-\mu _{1})t\\&-\chi _{4}f_{1}\left[{\breve {\Gamma }}_{2}(a_{1},a_{4})+{\frac {2\mu _{1}}{2(\mu _{4}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{4})\right]\cos 2(\mu _{4}-\mu _{1})t+\ldots \\&-\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\\end{aligned}}\right\}=0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/af2dc2fd1e5d9f9324908ce6aff7c764070892aa)
![{\displaystyle \left.{\begin{aligned}2^{\circ }&\ \ {\frac {d^{2}y_{1}}{dt^{2}}}+2\mu _{1}x_{1}-f_{1}h_{1}\\&-\chi _{2}f_{1}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{2})}{\mu _{2}-\mu _{1}}}\cos(\mu _{2}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{2})}{2(\mu _{2}-\mu _{1})}}\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&-\chi _{3}f_{1}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{3})}{\mu _{3}-\mu _{1}}}\cos(\mu _{3}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{3})}{2(\mu _{3}-\mu _{1})}}\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&-\chi _{4}f_{1}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{4})}{\mu _{4}-\mu _{1}}}\cos(\mu _{4}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{4})}{2(\mu _{4}-\mu _{1})}}\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&+\mathrm {K} _{1}{\frac {3f_{1}}{4(m-\mu _{1})}}\cos 2(m-\mu _{1})t\\\end{aligned}}\right\}=0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/214f0faf410d88e94c2f4220be211377b38ea530)

XXXVI.
La première équation étant intégrée par la méthode de l’Article XXXIV, on trouvera que la valeur de
renferme premièrement le terme constant
lequel devant être nul (Article XXXIII), on aura l’équation
Ensuite la valeur de
renfermera deux termes, tels que
et
avec des coefficients arbitraires, lesquels pourront se réduire à un seul terme représenté par
étant pareillement des constantes arbitraires.
De cette manière, si l’on fait
![{\displaystyle {\begin{aligned}&\Xi _{1}(a_{1},a_{2})=\left[{\breve {\Gamma }}_{1}(a_{1},a_{2})+{\frac {2\mu _{1}}{\mu _{2}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\right]{\frac {f_{1}}{(\mu _{2}-\mu _{1})^{2}-\mathrm {M} _{1}^{2}}},\\&\Xi _{2}(a_{1},a_{2})=\left[{\breve {\Gamma }}_{2}(a_{1},a_{2})+{\frac {2\mu _{1}}{2(\mu _{2}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{2})\right]{\frac {f_{1}}{4(\mu _{2}-\mu _{1})^{2}-\mathrm {M} _{1}^{2}}},\\&\Xi _{3}(a_{1},a_{2})=\left[{\breve {\Gamma }}_{3}(a_{1},a_{2})+{\frac {2\mu _{1}}{3(\mu _{2}-\mu _{1})}}{\widehat {\Gamma }}_{3}(a_{1},a_{2})\right]{\frac {f_{1}}{9(\mu _{2}-\mu _{1})^{2}-\mathrm {M} _{1}^{2}}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/edf6f045beee245aa588686bb7abbbc85dcf2b4f)
et de même
![{\displaystyle {\begin{aligned}&\Xi _{1}(a_{1},a_{3})=\left[{\breve {\Gamma }}_{1}(a_{1},a_{3})+{\frac {2\mu _{1}}{\mu _{3}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\right]{\frac {f_{1}}{(\mu _{3}-\mu _{1})^{2}-\mathrm {M} _{1}^{2}}},\\&\Xi _{2}(a_{1},a_{3})=\left[{\breve {\Gamma }}_{2}(a_{1},a_{3})+{\frac {2\mu _{1}}{2(\mu _{3}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{3})\right]{\frac {f_{1}}{4(\mu _{3}-\mu _{1})^{2}-\mathrm {M} _{1}^{2}}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a16eb8d236027cc0addc8eb2952dd6273ff04708)
et ainsi des autres, et qu’on suppose de plus

on aura
![{\displaystyle {\begin{aligned}x_{1}=&\varepsilon _{1}\cos(M_{1}t+\omega _{1})\\&-\chi _{2}\left[\Xi _{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+\Xi _{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&-\chi _{3}\left[\Xi _{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+\Xi _{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&-\chi _{4}\left[\Xi _{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+\Xi _{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&-\mathrm {K} _{1}\beta _{1}\cos 2(m-\mu _{1})t.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7d9b1cf55aaab120df82a71e8424424fbfca64c)
XXXVII.
Ayant trouvé la valeur
on aura l’expression du rayon vecteur
de l’orbite du premier satellite rapportée au plan de l’orbite de Jupiter, par la formule
[Article IV].
Or, en examinant cette expression de
on reconnaîtra aisément que le terme
représente l’équation elliptique qui vient de l’excentricité de l’orbite, de sorte que
exprimera la valeur de l’excentricité, et
sera l’anomalie moyenne ; d’où l’on voit que le mouvement de cette anomalie sera au mouvement moyen du satellite comme
à
par conséquent le mouvement moyen de la ligne des apsides sera au mouvement moyen du satellite comme
à
Nous verrons plus bas (Article XLV), qu’en négligeant les quantités de l’ordre
on a
de sorte que la ligne des apsides sera fixe, au moins par cette première approximation.
À l’égard de
on le déterminera par le moyen d’une époque quelconque donnée de l’anomalie moyenne ; ainsi les quantités
et
dépendent entièrement des observations.
Les autres termes de la valeur de
expriment les inégalités qui viennent de l’action des trois satellites
et du Soleil sur le satellite
XXXVIII.
On substituera maintenant la valeur de
dans la seconde équation de l’Article XXXV, et l’on tirera par l’intégration la valeur de
mais on aura attention de faire évanouir auparavant (Article XXXIII) le terme constant
ce qui donnera
Soit, pour abréger,

et pareillement

et ainsi de suite.
Soit de plus

on trouvera
![{\displaystyle {\begin{aligned}y_{1}=&-{\frac {2\mu _{1}\varepsilon _{1}}{\mathrm {M} _{1}}}\sin(\mathrm {M} _{1}t+\omega _{1})\\&+\chi _{2}\left[\Phi _{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+\Phi _{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&+\chi _{3}\left[\Phi _{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+\Phi _{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&+\chi _{4}\left[\Phi _{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+\Phi _{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&+\mathrm {K} _{1}\gamma _{1}\sin 2(m-\mu _{1})t.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f782638f08da810d7a9b9e20a6f94fa67702316)
XXXIX.
Puisque
(Article IV), on aura, en connaissant
l’expression du mouvement vrai du premier satellite par son mouvement moyen.
Le terme
représentera l’équation du centre qui vient de la figure elliptique de l’orbite, et les termes suivants exprimeront les inégalités causées par l’action des trois autres satellites et du Soleil.
XL.
Enfin l’équation

donnera

et
étant des quantités arbitraires ; car il est visible que cette expression
laquelle représente généralement la valeur de
(Article XXXIV), peut se réduire à celle-ci : 
XLI.
On aura donc, à cause de
(Article IV),


tangente de la latitude du satellite par rapport au plan de l’orbite de Jupiter ;
d’où l’on voit que l’orbite réelle du satellite sera toute dans un plan passant par le centre de Jupiter, et dont on connaîtra la position en remarquant : 1o que
étant la plus grande valeur de
exprimera la tangente de l’inclinaison ; 2o que
sera la distance du satellite au nœud ascendant, comptée sur l’orbite de Jupiter, laquelle étant retranchée de la longitude moyenne
on aura
pour la longitude moyenne du nœud.
Au reste, puisque l’on a ici
(Article XXXV), le mouvement du nœud sera nul, et sa longitude moyenne sera
ou plutôt
quantité qui dépend des observations ; mais il faut se souvenir que ce résultat n’est exact qu’aux quantités de l’ordre
près.
XLII.
On trouvera de même, pour le second satellites, les formules suivantes
![{\displaystyle {\begin{aligned}x_{2}=&\varepsilon _{2}\cos(M_{2}t+\omega _{2})\\&-\chi _{1}\left[\Xi _{1}(a_{2},a_{1})\cos(\mu _{1}-\mu _{2})t+\Xi _{2}(a_{2},a_{1})\cos 2(\mu _{1}-\mu _{2})t+\ldots \right]\\&-\chi _{3}\left[\Xi _{1}(a_{2},a_{3})\cos(\mu _{3}-\mu _{2})t+\Xi _{2}(a_{2},a_{3})\cos 2(\mu _{3}-\mu _{2})t+\ldots \right]\\&-\chi _{4}\left[\Xi _{1}(a_{2},a_{4})\cos(\mu _{4}-\mu _{2})t+\Xi _{2}(a_{2},a_{4})\cos 2(\mu _{4}-\mu _{2})t+\ldots \right]\\&-\mathrm {K} _{2}\beta _{2}\cos 2(m-\mu _{2})t\,;\\\\y_{2}=&-{\frac {2\mu _{2}\varepsilon _{2}}{\mathrm {M} _{2}}}\sin(\mathrm {M} _{2}t+\omega _{2})\\&+\chi _{1}\left[\Phi _{1}(a_{2},a_{1})\sin(\mu _{1}-\mu _{2})t+\Phi _{2}(a_{2},a_{1})\sin 2(\mu _{1}-\mu _{2})t+\ldots \right]\\&+\chi _{3}\left[\Phi _{1}(a_{2},a_{3})\sin(\mu _{3}-\mu _{2})t+\Phi _{2}(a_{2},a_{3})\sin 2(\mu _{3}-\mu _{2})t+\ldots \right]\\&+\chi _{4}\left[\Phi _{1}(a_{2},a_{4})\sin(\mu _{4}-\mu _{2})t+\Phi _{2}(a_{2},a_{4})\sin 2(\mu _{4}-\mu _{2})t+\ldots \right]\\&+\mathrm {K} _{2}\gamma _{2}\sin 2(m-\mu _{2})t\,;\\\\z_{2}=&\lambda _{2}\sin(\mu _{2}t+\eta _{2}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/429d3d990896771c519034c55b15c5f796c72a80)
Et l’on aura ensuite

Quant aux quantités marquées par
et
on aura
![{\displaystyle {\begin{aligned}&\Xi _{1}(a_{2},a_{1})=\left[{\breve {\Gamma }}_{1}(a_{2},a_{1})+{\frac {2\mu _{2}}{\mu _{1}-\mu _{2}}}{\widehat {\Gamma }}_{1}(a_{2},a_{1})\right]{\frac {f_{2}}{(\mu _{1}-\mu _{2})^{2}-\mathrm {M} _{2}^{2}}},\\&\Xi _{2}(a_{2},a_{1})=\left[{\breve {\Gamma }}_{2}(a_{2},a_{1})+{\frac {2\mu _{2}}{2(\mu _{1}-\mu _{2})}}{\widehat {\Gamma }}_{2}(a_{2},a_{1})\right]{\frac {f_{2}}{4(\mu _{1}-\mu _{2})^{2}-\mathrm {M} _{2}^{2}}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \,;\\&\Phi _{1}(a_{2},a_{1})={\frac {2\mu _{2}}{\mu _{1}-\mu _{2}}}\Xi _{1}(a_{2},a_{1})+{\frac {f_{2}}{(\mu _{1}-\mu _{2})^{2}}}{\widehat {\Gamma }}_{1}(a_{2},a_{1}),\\&\Phi _{2}(a_{2},a_{1})={\frac {2\mu _{2}}{2(\mu _{1}-\mu _{2})}}\Xi _{2}(a_{2},a_{1})+{\frac {f_{2}}{4(\mu _{1}-\mu _{2})^{2}}}{\widehat {\Gamma }}_{2}(a_{2},a_{1}),\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4d337afe57eae78310d9697525c7d90ea48a446e)
Outre cela, on aura

Enfin on trouvera les deux conditions
et
qui servirout à déterminer les deux constantes
et
(Article XIX).
On aura des formules analogues pour le troisième et le quatrième satellite, que nous nous dispenserons de rappeler ici, parce qu’elles se déduisent à l’œil de celles que nous venons de donner.
§ II. — Valeurs numériques des coefficients des formules précédentes.
XLIII.
Soit, en général, suivant l’Article XVIII,

étant une fraction moindre que l’unité ; on aura, en faisant
et 
![{\displaystyle \left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {3}{2}}}=}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e0cba3c89d6688add60aa89f771de73a9e19b8d1)
![{\displaystyle a_{2}^{3}\left[\mathrm {A+B\cos(\varphi _{2}-\varphi _{1})+C\cos 2(\varphi _{2}-\varphi _{1})+\ldots } \right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c5d4e9149043cf3fd00fc01c7286dc96919828fc)
donc (Article XXI)

On trouvera de même, en faisant successivement

on trouvera, dis-je,

et ainsi de suite.
À l’égard des quantités
il est évident qu’elles doivent être égales à leurs réciproques
car les fonctions

demeurent les mêmes, en changeant
en
et
en
ou bien
en
et
en 
XLIV.
De là on trouvera (Article XXIV),
étant égal à

Et, par l’Article XXV on aura

et ces mêmes formules serviront aussi pour les quantités

en faisant successivement

Mais, pour les quantités réciproques
on aura les formules suivantes

lesquelles auront lieu pareillement pour les quantités
en faisant comme ci-devant 
XLV.
On formera ensuite les quantités marquées par
et par
(Articles XXXIV et XXXVIII), ce qui n’aura aucune difficulté. On remarquera seulement que, à cause de
(Article XXIX), on aura, en négligeant la quantité
qui est de l’ordre

et de là

par conséquent (Article XXXV)

et pareillement

Donc, supposant
on aura
![{\displaystyle {\begin{aligned}&\Xi _{1}(a_{1},a_{2})=\left[{\breve {\Gamma }}_{1}(a_{1},a_{2})+{\frac {2}{s-1}}{\breve {\Gamma }}_{1}(a_{1},a_{2})\right]{\frac {1}{(s-1)^{2}-1}},\\&\Xi _{2}(a_{1},a_{2})=\left[{\breve {\Gamma }}_{2}(a_{1},a_{2})+{\frac {2}{2(s-1)}}{\breve {\Gamma }}_{2}(a_{1},a_{2})\right]{\frac {1}{4(s-1)^{2}-1}},\\&\Xi _{3}(a_{1},a_{2})=\left[{\breve {\Gamma }}_{3}(a_{1},a_{2})+{\frac {2}{3(s-1)}}{\breve {\Gamma }}_{3}(a_{1},a_{2})\right]{\frac {1}{9(s-1)^{2}-1}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \,;\\&\Phi _{1}(a_{1},a_{2})={\frac {2}{s-1}}\Xi _{1}(a_{1},a_{2})+{\frac {1}{(s-1)^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2}),\\&\Phi _{2}(a_{1},a_{2})={\frac {2}{2(s-1)}}\Xi _{2}(a_{1},a_{2})+{\frac {1}{4(s-1)^{2}}}{\widehat {\Gamma }}_{2}(a_{1},a_{2}),\\&\Phi _{3}(a_{1},a_{2})={\frac {2}{3(s-1)}}\Xi _{3}(a_{1},a_{2})+{\frac {1}{9(s-1)^{2}}}{\widehat {\Gamma }}_{3}(a_{1},a_{2}),\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cff7635ae82c1ebc65cbd5bba883a2c5fc35db65)
De méme, en faisant
on aura
![{\displaystyle {\begin{aligned}&\Xi _{1}(a_{1},a_{3})=\left[{\breve {\Gamma }}_{1}(a_{1},a_{3})+{\frac {2}{s-1}}{\widehat {\Gamma }}(a_{1},a_{3})\right]{\frac {1}{(s-1)^{2}-1}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \,;\\&\Phi _{1}(a_{1},a_{3})={\frac {2}{s-1}}\Xi _{1}(a_{1},a_{3})+{\frac {1}{(s-1)^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3}),\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/afacce161b78f2189c066ff96e79d27af59a583b)
Pareillement on trouvera,

étant égal à
![{\displaystyle {\begin{aligned}&\Xi _{1}(a_{2},a_{1})=\left[{\breve {\Gamma }}_{1}(a_{2},a_{1})+{\frac {2}{s-1}}{\widehat {\Gamma }}(a_{2},a_{1})\right]{\frac {1}{(s-1)^{2}-1}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \,;\\&\Phi _{1}(a_{2},a_{1})={\frac {2}{s-1}}\Xi _{1}(a_{2},a_{1})+{\frac {1}{(s-1)^{2}}}{\widehat {\Gamma }}_{1}(a_{2},a_{1}),\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f928475e3ac34b7be99e66bffb26c47dcbd80bd6)
et ainsi des autres.
XLVI.
Cela posé, on remarquera :
1o Que les quantités
expriment les vitesses angulaires moyennes des satellites
autour de Jupiter (Articles IV et IX) ; d’où il suit que ces quantités sont réciproquement proportionnelles aux temps périodiques des mêmes satellites.
Or on a, par les observations, en négligeant les secondes,
Révolutions périodiques.

Donc, réduisant ces quantités en minutes, on aura

2o Que l’on a généralement (Article XLV)

c’est-à-dire, à cause de
(Article V) et de
(Article X) 

et par conséquent

d’où l’on voit que les quantités
sont entre elles comme les quantités
ainsi l’on trouvera les valeurs de ces quantités, ou plutôt de leurs rapports, qui sont les seules dont nous ayons besoin ici.
Au reste, comine ces valeurs ne sont exactes qu’aux quantités de l’ordre
près, nous emploierons, pour les distances moyennes des satellites, les déterminations que M. Cassini a tirées des observations, lesquelles ne s’écartent d’ailleurs que très-peu de la loi de Kepler ; on aura donc, en prenant le demi-diamètre de Jupiter pour l’unité,

XLVII.
Par le moyen de ces valeurs numériques et des formules des Articles XVIII et XIX, j’ai trouvé les déterminations suivantes

Et de là (Article
XLIV)


D’où enfin (Article
XLV)


En consultant cette dernière Table, on voit qu’il y a des quantités dont les valeurs numériques sont fort grandes ; telles sont les quantités

et leurs correspondantes

La raison pour laquelle ces quantités ont des valeurs si considérables, c’est que le diviseur
se trouve fort petit dans le cas où
et
et que pareillement le diviseur
est fort petit lorsque
et
comme il est facile de s’en assurer, au moyen des valeurs de
données ci-dessus.
Cette remarque est d’autant plus essentielle qu’elle sert à expliquer pourquoi les équations empiriques des satellites de Jupiter sont en effet les seules qui puissent être bien sensibles (voir plus bas Articles LVIII et suivants).
XLVIII.
Il ne reste plus maintenant qu’à chercher les valeurs des quantités et (Articles XXXVI et XXXVIII).
Pour cela, on remarquera que la quantité
qui représente la vitesse angulaire moyenne du Soleil autour de Jupiter (Article XXVII), est extrêmement petite par rapport aux quantités
vitesses moyennes des satellites ; d’où il suit qu’elle pourra être négligée vis-à-vis de ces dernières quantités ; or on a généralement (Articles cités)

c’est-à-dire, à cause de
et
(Article XLV),

donc, en négligeant la quantité

on aura

À l’égard des quantités
et
qui doivent être déterminées par les équations
et
(Articles XXIX, XXXV et suivants), il est inutile d’en chercher la valeur, puisqu’elles ne se trouvent point dans l’expression des coefficients de nos formules.
§ III. — Formules des rayons vecteurs et des longitudes vraies des satellites de Jupiter par rapport au plan de l’orbite de cette Planète.
XLIX.
Dans les formules suivantes, j’ai remis, au lieu de
les quantités
(Article XXIX) ; et j’ai substitué pour
leurs valeurs en nombres ; car, puisque
(Article cité) et que
(Article XLVI), et par conséquent aussi
on aura
donc

Or,
étant la vitesse moyenne angulaire du Soleil autour de Jupiter et
les vitesses moyennes angulaires des satellites, on aura


d’où l’on tire à très-peu près


Outre cela, au lieu de
et
qui représentent les valeurs moyennes des angles
et
c’est-à-dire des longitudes moyennes des quatre satellites et du Soleil vus de Jupiter, et rapportés au plan de son orbite, j’ai mis les lettres
et
De même, au lieu de
anomalies moyennes des satellites, j’ai substitué, pour plus de simplicité,
Enfin j’ai exprimé les rayons vecteurs en demi-diamètres de Jupiter, et les longitudes en minutes. De cette manière j’ai trouvé
L.
Rayon vecteur du premier satellite.
![{\displaystyle {\begin{aligned}r_{1}=&5{,}67(1+n\varepsilon _{1}\cos \mathrm {V} _{1})\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[4{,}19\cos(u_{2}-u_{1})-1014{,}93\cos 2(u_{2}-u_{1})\right.\\&\qquad \qquad \qquad \left.-3{,}87\cos 3(u_{2}-u_{1})-1{,}02\cos 4(u_{2}-u_{1})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[0{,}75\cos(u_{3}-u_{1})-1{,}06\cos 2(u_{3}-u_{1})\right.\\&\qquad \qquad \qquad \left.-0{,}13\cos 3(u_{3}-u_{1})-0{,}02\cos 4(u_{3}-u_{1})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[0{,}14\cos(u_{4}-u_{1})-0{,}09\cos 2(u_{4}-u_{1})\right.\\&\qquad \qquad \qquad \left.-0{,}00\cos 3(u_{4}-u_{1})-0{,}00\cos 4(u_{4}-u_{1})\ldots \right]\\&-0{,}00000092\cos 2(v-u_{1}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a7a47eb3f2b76e3d6b08ae7f6e4b3ea4bd0be22)
Rayon vecteur du deuxième satellite.
![{\displaystyle {\begin{aligned}r_{2}=&9{,}00(1+n\varepsilon _{2}\cos \mathrm {V} _{2})\\&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[518{,}78\cos(u_{1}-u_{2})+5{,}73\cos 2(u_{1}-u_{2})\right.\\&\qquad \qquad \qquad \left.+1{,}36\cos 3(u_{1}-u_{2})+0{,}49\cos 4(u_{1}-u_{2})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[7{,}16\cos(u_{3}-u_{2})-824{,}07\cos 2(u_{3}-u_{2})\right.\\&\qquad \qquad \qquad \left.-6{,}29\cos 3(u_{3}-u_{2})-1{,}66\cos 4(u_{3}-u_{2})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[0{,}52\cos(u_{4}-u_{1})-4{,}85\cos 2(u_{4}-u_{1})\right.\\&\qquad \qquad \qquad \left.-0{,}11\cos 3(u_{4}-u_{1})-0{,}04\cos 4(u_{4}-u_{1})\ldots \right]\\&-0{,}0000060\cos 2(v-u_{2}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/266efa7b6fac44a941b28328465f68b864b8e9a3)
Rayon vecteur du troisième satellite.
![{\displaystyle {\begin{aligned}r_{3}=&14{,}38(1+n\varepsilon _{3}\cos \mathrm {V} _{3})\\&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[5{,}88\cos(u_{1}-u_{3})+0{,}19\cos 2(u_{1}-u_{3})\right.\\&\qquad \qquad \qquad \left.+0{,}03\cos 3(u_{1}-u_{3})+0{,}00\cos 4(u_{2}-u_{1})\ldots \right]\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[452{,}98\cos(u_{2}-u_{3})+9{,}13\cos 2(u_{2}-u_{3})\right.\\&\qquad \qquad \qquad \left.+2{,}16\cos 3(u_{2}-u_{3})+0{,}59\cos 4(u_{2}-u_{3})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[7{,}46\cos(u_{4}-u_{3})-33{,}62\cos 2(u_{4}-u_{3})\right.\\&\qquad \qquad \qquad \left.-3{,}95\cos 3(u_{4}-u_{3})+0{,}95\cos 4(u_{4}-u_{3})\ldots \right]\\&-0{,}0000392\cos 2(v-u_{3}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69f2bde76a70cc9c6c7b617b38f51851f59dd352)
Rayon vecteur du quatrième satellite.
![{\displaystyle {\begin{aligned}r_{4}=&25{,}30(1+n\varepsilon _{4}\cos \mathrm {V} _{4})\\&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[5{,}66\cos(u_{1}-u_{4})+0{,}01\cos 2(u_{1}-u_{4})\right.\\&\qquad \qquad \qquad \left.+0{,}00\cos 3(u_{1}-u_{4})+0{,}00\cos 4(u_{1}-u_{4})\ldots \right]\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[8{,}18\cos(u_{2}-u_{4})+0{,}17\cos 2(u_{2}-u_{4})\right.\\&\qquad \qquad \qquad \left.+0{,}07\cos 3(u_{2}-u_{4})+0{,}00\cos 4(u_{3}-u_{2})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[31{,}55\cos(u_{3}-u_{4})+1{,}62\cos 2(u_{3}-u_{4})\right.\\&\qquad \qquad \qquad \left.+1{,}35\cos 3(u_{3}-u_{4})+0{,}43\cos 4(u_{3}-u_{4})\ldots \right]\\&-0{,}0003754\cos 2(v-u_{4}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f32576a636dcc1b94e3c90fe801d32456d28fae)
LI.
Longitude vraie du premier satellite.
![{\displaystyle {\begin{aligned}\varphi _{1}=&u_{1}-\left(114^{\circ }35'\right)n\varepsilon _{1}\sin \mathrm {V} _{1})\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[9500'\sin(u_{2}-u_{1})-1227214'\sin 2(u_{2}-u_{1})\right.\\&\qquad \qquad \qquad \left.-3526'\sin 3(u_{2}-u_{1})-705'\sin 4(u_{2}-u_{1})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[1158'\sin(u_{3}-u_{1})-1007'\sin 2(u_{3}-u_{1})\right.\\&\qquad \qquad \qquad \left.-101'\sin 3(u_{3}-u_{1})-18'\sin 4(u_{3}-u_{1})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[188'\sin(u_{4}-u_{1})-81'\sin 2(u_{4}-u_{1})\right.\\&\qquad \qquad \qquad \left.-3'\sin 3(u_{4}-u_{1})-0'\sin 4(u_{4}-u_{1})\ldots \right]\\&-0{,}047\sin 2(v-u_{1}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11a2dcbd5c86cd766bf1a8a69d817f7252c2655f)
Longitude vraie du deuxième satellite.
![{\displaystyle {\begin{aligned}\varphi _{2}=&u_{2}-\left(114^{\circ }35'\right)n\varepsilon _{2}\sin \mathrm {V} _{2})\\&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[-387482'\sin(u_{1}-u_{2})-2727'\sin 2(u_{1}-u_{2})\right.\\&\qquad \qquad \qquad \left.-509'\sin 3(u_{1}-u_{2})-12'\sin 4(u_{1}-u_{2})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[10067'\sin(u_{3}-u_{2})-626246'\sin 2(u_{3}-u_{2})\right.\\&\qquad \qquad \qquad \left.-3717'\sin 3(u_{3}-u_{2})-825'\sin 4(u_{3}-u_{2})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[506'\sin(u_{4}-u_{2})-3123'\sin 2(u_{4}-u_{2})\right.\\&\qquad \qquad \qquad \left.-52'\sin 3(u_{4}-u_{2})-17'\sin 4(u_{4}-u_{2})\ldots \right]\\&-0{,}190\sin 2(v-u_{2}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/33e3c383b79d93e7465c4110e3019ecb96dd8290)
Longitude vraie du troisième satellite.
![{\displaystyle {\begin{aligned}\varphi _{3}=&u_{3}-\left(114^{\circ }35'\right)n\varepsilon _{3}\sin \mathrm {V} _{3})\\&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[1306'\sin(u_{1}-u_{3})-38'\sin 2(u_{1}-u_{3})\right.\\&\qquad \qquad \qquad \left.-8'\sin 3(u_{1}-u_{3})-1'\sin 4(u_{1}-u_{3})\ldots \right]\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[-207375'\sin(u_{2}-u_{3})-2760'\sin 2(u_{2}-u_{3})\right.\\&\qquad \qquad \qquad \left.-559'\sin 3(u_{2}-u_{3})-142'\sin 4(u_{2}-u_{3})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[5703'\sin(u_{4}-u_{3})-14911'\sin 2(u_{4}-u_{3})\right.\\&\qquad \qquad \qquad \left.-1376'\sin 3(u_{4}-u_{3})-306'\sin 4(u_{4}-u_{3})\ldots \right]\\&-0{,}977\sin 2(v-u_{3}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d491acb54f0272989df35b24b9fa687dd848a06c)
Longitude vraie du quatrième satellite.
![{\displaystyle {\begin{aligned}\varphi _{4}=&u_{4}-\left(114^{\circ }35'\right)n\varepsilon _{4}\sin \mathrm {V} _{4})\\&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[768'\sin(u_{1}-u_{4})-0'\sin 2(u_{1}-u_{4})\right.\\&\qquad \qquad \qquad \left.-0'\sin 3(u_{1}-u_{4})-0'\sin 4(u_{1}-u_{4})\ldots \right]\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[1219'\sin(u_{2}-u_{4})-15'\sin 2(u_{2}-u_{4})\right.\\&\qquad \qquad \qquad \left.-3'\sin 3(u_{2}-u_{4})-0'\sin 4(u_{2}-u_{4})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[-1691'\sin(u_{3}-u_{4})-444'\sin 2(u_{3}-u_{4})\right.\\&\qquad \qquad \qquad \left.-180'\sin 3(u_{3}-u_{4})-59'\sin 4(u_{3}-u_{4})\ldots \right]\\&-4{,}208\sin 2(v-u_{4}).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/64ced92eddffd734816a8c11212faf7680d6e4f7)
§ IV. — Où l’on donne les inégalités des satellites qui dépendent de leurs configurations, et, qui ont lieu au temps des éclipses.
LII.
Il est visible que les éclipses des satellites, c’est-à-dire leurs conjonctions avec Jupiter, arrivent lorsque leurs longitudes
diffèrent de
degrés de la longitude
du Soleil vu de Jupiter ; de sorte qu’on aura généralement l’équation

ou bien, en mettant pour
et
leurs valeurs
et 

où
exprime l’équation de Jupiter,
l’équation, ou plutôt la somme des équations du satellite, et
la distance, ou bien l’élongation moyenne du satellite ; donc, pour avoir la conjonction vraie, il n’y aura qu’à ajouter au temps de la conjonction moyenne la quantité
convertie en temps, à raison du mouvement moyen du satellite au Soleil, conversion qu’on fera aisément en multipliant la quantité proposée par
pour le premier satellite, par
pour le deuxième, par
pour le troisième et par
pour le quatrième ; et changeant ensuite les degrés en heures, les minutes de degré en minutes de temps, etc. ces nombres se trouvent, en divisant les durées des révolutions synodiques des satellites, lesquelles sont de

par
degrés, après avoit réduit le tout en secondes.
LIII.
Nous allons donc donner ici les équations des conjonctions des satellites ; mais nous ferons abstraction de celles qui viennent, de l’excentricité de Jupiter, et des excentricités particulières des satellites, parce que les unes sont assez connues des Astronomes, et que les autres ne sont pas assez exactes pour qu’on puisse s’y fier.
LIV.
Équation du premier satellite.
![{\displaystyle {\begin{aligned}&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[-1097^{\text{m}}\sin(u_{2}-u_{1})+144800^{\text{m}}\sin 2(u_{2}-u_{1})\right.\\&\qquad \qquad \qquad \qquad \left.+416^{\text{m}}\sin 3(u_{2}-u_{1})+83^{\text{m}}\sin 4(u_{2}-u_{1})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[-137^{\text{m}}\sin(u_{3}-u_{1})+119^{\text{m}}\sin 2(u_{3}-u_{1})\right.\\&\qquad \qquad \qquad \qquad \left.+12^{\text{m}}\sin 3(u_{3}-u_{1})-1^{\text{m}}\sin 4(u_{3}-u_{1})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[-22^{\text{m}}\sin(u_{4}-u_{1})+9^{\text{m}}\sin 2(u_{4}-u_{1})\right.\\&\qquad \qquad \qquad \qquad \left.+0^{\text{m}}\sin 3(u_{4}-u_{1})+0^{\text{m}}\sin 4(u_{4}-u_{1})\ldots \right]\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7553d38ac0b805f7cc0d321f64b24a64af906b43)
LV.
Équation du deuxième satellite.
![{\displaystyle {\begin{aligned}&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[91810^{\text{m}}\sin(u_{1}-u_{2})+646^{\text{m}}\sin 2(u_{1}-u_{2})\right.\\&\qquad \qquad \qquad \qquad \left.+121^{\text{m}}\sin 3(u_{1}-u_{2})+3^{\text{m}}\sin 4(u_{1}-u_{2})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[-2385^{\text{m}}\sin(u_{3}-u_{2})+148383^{\text{m}}\sin 2(u_{3}-u_{2})\right.\\&\qquad \qquad \qquad \qquad \left.+881^{\text{m}}\sin 3(u_{3}-u_{2})+195^{\text{m}}\sin 4(u_{3}-u_{2})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[-119^{\text{m}}\sin(u_{4}-u_{2})+740^{\text{m}}\sin 2(u_{4}-u_{2})\right.\\&\qquad \qquad \qquad \qquad \left.+12^{\text{m}}\sin 3(u_{4}-u_{2})+4^{\text{m}}\sin 4(u_{4}-u_{2})\ldots \right]\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/70ae7cec6c2af45bbc3210cbe5d54e7f685ef507)
LV.
Équation du troisième satellite.
![{\displaystyle {\begin{aligned}&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[-624^{\text{m}}\sin(u_{1}-u_{3})+18^{\text{m}}\sin 2(u_{1}-u_{3})\right.\\&\qquad \qquad \qquad \qquad \left.+4^{\text{m}}\sin 3(u_{1}-u_{3})+0^{\text{m}}\sin 4(u_{1}-u_{3})\ldots \right]\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[99075^{\text{m}}\sin(u_{2}-u_{3})+1319s^{\text{m}}\sin 2(u_{2}-u_{3})\right.\\&\qquad \qquad \qquad \qquad \left.+267^{\text{m}}\sin 3(u_{2}-u_{3})+68^{\text{m}}\sin 4(u_{2}-u_{3})\ldots \right]\\&+{\frac {{\text{☾}}_{4}}{\mathbb {Z} \!^{\upsilon }}}\left[-2725^{\text{m}}\sin(u_{4}-u_{3})+7124^{\text{m}}\sin 2(u_{4}-u_{3})\right.\\&\qquad \qquad \qquad \qquad \left.+627^{\text{m}}\sin 3(u_{4}-u_{3})+146^{\text{m}}\sin 4(u_{4}-u_{3})\ldots \right]\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e6a7a35b8cd4ecd39eeedffdae54d4b8f9437629)
LV.
Équation du quatrième satellite.
![{\displaystyle {\begin{aligned}&+{\frac {{\text{☾}}_{1}}{\mathbb {Z} \!^{\upsilon }}}\left[-858^{\text{m}}\sin(u_{1}-u_{4})+0^{\text{m}}\sin 2(u_{1}-u_{4})\ldots \right]\\&+{\frac {{\text{☾}}_{2}}{\mathbb {Z} \!^{\upsilon }}}\left[-1362^{\text{m}}\sin(u_{2}-u_{4})+18^{\text{m}}\sin 2(u_{2}-u_{4})\right.\\&\qquad \qquad \qquad \qquad \left.+4^{\text{m}}\sin 3(u_{2}-u_{4})-0^{\text{m}}\sin 4(u_{2}-u_{4})\ldots \right]\\&+{\frac {{\text{☾}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[1888^{\text{m}}\sin(u_{3}-u_{4})+496^{\text{m}}\sin 2(u_{3}-u_{4})\right.\\&\qquad \qquad \qquad \qquad \left.+201^{\text{m}}\sin 3(u_{3}-u_{4})+66^{\text{m}}\sin 4(u_{3}-u_{4})\ldots \right]\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dcbe7536e9678ed8fab99be4b7ea8915614ecbae)
J’ai négligé dans ces formules les termes dus à l’action du Soleil, et qui sont de la forme de
parce que ces termes deviennent nuls au temps des conjonctions où l’on a
§ V. — Comparaison des formules précédentes avec les observations, et conséquences qui en résultent par rapport aux masses des satellites.
LVIII.
Nous nous contenterons ici de comparer nos formoles avec les Tables de M. Wargentin, qui sont, comme l’on sait, le résultat d’un grand nombres d’observations ; mais, avant d’entreprendre ce parallèle, il est bon d’avertir que les Tables de ce grand Astronome sont dressées de manièce qu’il n’y a aucune équation soustractive, quoique les équations qu’il a employées soient de nature à être tantôt additives et tantôt soustractives ; cela vient de ce que l’Auteur a retranché, par avance des époques, chacune des plus grandes équations soustractives ; de sorte que les équations des Tables se trouvent nulles dans le cas où elles auraient été les plus grandes à soustraire, et que leur plus grande valeur est double de ce qu’elle aurait dû être.
LIX.
En examinant les différents termes de la formule de l’Article LIV on voit qu’il y en a un dont le coefficient numérique est très-grand, et vis-à-vis duquel tous les autres termes ne sont presque d’aucune considération c’est le terme

d’où résulte une équation qui a pour argument
savoir le double de la distance moyenne du second satellite au premier, au temps des conjonctions de celui-ci, et dont la plus grande valeur est de
,
exprimant le rapport de la masse du second satellite à celle de Jupiter.
Pour mieux connaître la nature de cette inégalité qui doit avoir lieu dans les conjonctions du premier satellite, il faut cherchec sa période, laquelle dépend du rapport des révolutions synodiques des deux premiers satellites. Or, suivant M. Wargentin, on a pour la durée de la révolution synodique du premier

et pour celle du second

d’où l’on trouve, en additionnant successivement ces nombres, que
révolutions du premier font

et que
révolutions du second font

ainsi, pendant que le premier fait une révolution par rapport au Soleil, le second ne fait que
d’une pareille révolution ; d’où il suit que la distance
u, du second satellite au premier augmente, dans l’intervalle d’une conjonction à l’autre, de
pour avoir une exactitude plus grande, on additionnera de nouveau les périodes du premier et du second satellite que nous venons de trouver, jusqu’à ce qu’ils fassent des sommes à peu près égales, et l’on trouvera que
révolutions du premier font

et que
révolutions du second font

c’est pourquoi on aura, au lieu de la fraction
celle-ci beaucoup plus exacte 
Soit maintenant
la distance du second satellite au premier au temps d’une conjonction de celui-ci, cette distance deviendra, après
n
révolutions[7],

donc on aura

et

donc, pour que cette quantité redevienne
il faut que

c’est le nombre des révolutions du premier satellite qui exprime la période de l’équation 
Or
révolutions font à très-peu près

et
de révolution font
donc la période cherchée sera de

LX.
Voyons à présent quelle doit être la marche de cette équation ; pour cela, nous supposerons
c’est-à-dire que les deux satellites se trouvent à la fois en conjonction, et nous aurons, après un nombre quelconque
de révolutions du premier satellite,

ou bien en faisant, pour abréger,

égal au nombre des révolutions qui forment la période de l’équation

De là on voit que l’équation
sera nulle au commencement de la période, qu’ensuite elle deviendra soustractive, et qu’elle sera la plus grande à soustraire lorsque
c’est-à-dire, au quart de la période ; après quoi elle redeviendra nulle à la moitié de la période, ensuite se changera en additive croissante jusqu’aux trois quarts de la période, où elle sera la plus grande, et enfin décroîtra pendant le dernier quart, pour se retrouver nulle au commencement de la période suivante.
LXI.
Je dis maintenant que l’équation que nous venons d’examiner est la même que celle qui se trouve dans les Tables du premier satellite, désignée par la lettre
et qui est la seule que les observations aient fait connaître jusqu’ici. En effet : 1o la période de cette équation est, selon M. Wargentin, de
environ, ce qui s’accorde admirablement bien avec ce que nous avons trouvé dans l’Article LIX ; car la différence de
qui s’y trouve, n’est d’aucune considération par rapport à un intervalle de
jours ; 2o si l’on examine l’équation
on verra qu’en ôtant toujours
(moitié de la plus grande valeur de cette équation, selon la remarque de l’Article LVIII), et établissant le commencement de la période (qui est divisée en
parties) au nombre
on verra, dis-je, que la marche de cette équation est la même que celle de l’équation
de l’Article précédent. De plus on trouvera, par les Tables du premier et du second satellite, que, dans les conjonctions du premier satellite qui répondent exactement au nombre
l’élongation du second satellite est nulle. Donc, etc.
LXII.
De là il suit que les nombres
des Tables du premier satellite ne sont autre chose que les distances, c’est-à-dire, les élongations du second satellite au premier, au temps des conjonctions de celui-ci, le cercle étant supposé divisé en
parties, de sorte que le nombre
réponde aux conjonctions des deux satellites et
à leurs oppositions. Cette remarque fournit un moyen de rectifier les époques de ces nombres, si elles en avaient besoin, et de les prolonger autant qu’on voudra, sans craindre de s’égarer.
LXIII.
La plus grande valeur de l’équation
du premier satellite est de
dont il ne faut prendre que la moitié (Article LVIII) ; donc, comparant cette valeur avec le coefficient de l’équation
lequel est
on aura

d’où l’on tire

à peu près ;
c’est le rapport de la masse, du second satellite à celle de Jupiter. Si l’on prend la masse de la Terre pour l’unité, on a

ce qui donne

à peu près.
Supposons que la densité de ce satellite soit la même que celle de Jupiter, ou au moins qu’elle n’en diffère que très-peu, ce qui est trèsnaturel, on trouvera, en prenant le demi-diamètre de Jupiter pour l’unité, que celui du satellite est
c’est-à-dire, environ
ce qui donnerait pour le temps que le satellite doit employer à entrer dans l’ombre de Jupiter
ce qui est à peu près le milieu entre les résultats des observations de M. Maraldi et de M. Whiston (Mémoires d’e l’Académie, 1734).
LXIV.
Il serait tout à fait inutile d’examiner les autres termes de la formule de l’Article LIV ; car il est clair qu’il n’en pourrait résulter que des équations extrêmement petites, et par conséquent insensibles, à moins qu’on ne voulût supposer les masses du troisième et du quatrième satellite énormément grandes par rapport à celle du second, ce qui ne paraît guère naturel ; d’ailleurs l’équation que nous avons examinée est la seule qu’on ait jusqu’ici déduite des observations.
LXV.
Passons donc à la formule de l’Article LV, qui renferme les équations des conjonctions du second satellite. Parmi tous les termes dont cette formule est composée, j’en distingue d’abord deux qui sont beaucoup plus considérables que les autres par les coefficients numériques dont ils sont affectés, savoir

dont l’un vient de l’action du premier satellite, et l’autre de l’action du troisième. Ces deux termes produisent, comme l’on voit, deux équations dont les arguments sont
distance moyenne du premier satellite au second, et
double de la distance moyenne du troisième satellite au second au temps des conjonctions de celui-ci.
Je remarque maintenant que la durée de la révolution synodique du troisième satellite est de

selon M. Wargentin ; ce qui donne, pour
révolutions,

et, pour
111021
révolutions,

or nous avons déjà trouvé que
révolutions du premier sont

et que
révolutions du second sont (Article LIX)

donc les mouvements des trois premiers satellites au Soleil sont entre eux comme les nombres
et les différences entre les mouvements des deux premiers et les mouvements du second et du troisième sont au mouvement du second comme les nombres
au nombre
donc, pendant que le second achève une révolution au Soleil, les angles
et
croissent de
et
donc l’angle
diminue à chaque révolution du second de la même quantité dont l’angle
augmente, savoir de
donc la quantité

est toujours la même dans les conjonctions du second satellite.
Examinons donc une conjonction quelconque de ce satellite, et voyons quelles sont les élongations du premier et du troisième, c’est-à-dire, les valeurs de
et de
je prends pour exemple la première conjonction de l’année 1760, laquelle est marquée dans les Tables à
à quoi ajoutant la moitié des plus grandes équations, savoir,
(Article LVIII), on a

pour le temps moyen de la conjonction moyenne du second satellite ; je trouve de la même manière que les premières conjonctions moyennes
du premier et du troisième satellite ont dû arriver à

de temps moyen ; d’où je conclus qu’au temps de la conjonction du second satellite, le premier était plus avancé de
ce qui fait
et que le troisième était en arrière de
ce qui vaut
donc

par conséquent

à très-peu près.
On aura donc, en général,

ainsi les deux termes

peuvent se réduire à un terme unique, tel que

lequel ne donne qu’une équation dépendante de l’élongation du premier satellite au second.
LXVI.
Soit, dans une conjonction du second satellite,
on aura, par ce que nous avons démontré dans l’Article précédent, après
révolutions de ce même satellite,

et par conséquent

d’où l’on voit que cette quantité ne peut redevenir

à moins que l’on n’ait

ce qui donne

c’est le nombre des révolutions du second satellite qui forment la période de l’équation
et l’on trouvera que cette période est la même que celle de l’équation du premier satellite, savoir (Article LIX)

Mettons
au lieu du nombre
nous aurons

donc, supposant au commencement de la période
c’est-à-dire, les deux premiers satellites en conjonction à la fois, et faisant successivement
et
on trouvera que l’équation dont il s’agit est nulle au commencement de la période, qu’ensuite elle augmente jusqu’au quart de la période, où elle est la plus grande ; que de là elle diminue et redevient nulle à la moitié de la période, après quoi elle se change en négative, etc.
LXVII.
Si l’on compare maintenant la marche de cette équation avec celle de l’équation
des Tables du second satellite, on verra qu’elles s’accordent parfaitement, pourvu que l’on ait attention d’ôter constamment de cette dernière équation
moitié de sa plus grande valeur, et qu’on fixe le commencement de la période au nombre
Ainsi les nombres
des Tables du second satellite indiquent les élongations du premier au temps des conjonctions du second, de sorte que le nombre
répond aux conjonctions des deux satellites, et le nombre
à leurs oppositions. (Voyez là-dessus la dissertation de M. Wargentin qui est à la tête des observations du second satellite, dans les Mémoires de la Société d’Upsal pour l’année 1743.)
LXVIII.
Il ne reste donc plus qu’à égaler le coefficient de l’équation
à la plus grande valeur de l’équation
des Tables, ce qui donne

de sorte qu’en supposant
on aura

Soit par exemple
c’esl-à-dire, les masses du premier et du troisième satellite égales entre elles, on aura

environ ;
d’où, en supposant les densités des satellites égales à celles de Jupiter, on tire leurs demi-diamètres
environ
de celui de Jupiter ; ce qui donne pour le temps que le premier devrait employer à entrer dans l’ombre
et pour le temps que devrait employer le troisième 
Au reste, quel que soit le nombre
comme il ne saurait être ni infini ni nul, il est clair que les quantités
sont toujours nécessairement moindres que la fraction
c’est-à-dire, en prenant la masse de la Terre pour l’unité,

et

LXIX.
À l’égard des autres termes de la formule de l’Article LV il est facile de voir qu’ils ne donnent que des équations extrêmement petites, et qui peuvent par conséquent être négligées ; en effet, le terme qui a le plus grand coefficient numérique, après ceux que nous venons d’examiner, est
![{\displaystyle {\frac {{\mathfrak {S}}_{3}}{\mathbb {Z} \!^{\upsilon }}}\left[-2385^{\text{m}}\sin(u_{3}-u_{2})\right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/415188ca542b8e67f9daa41a394b8c17983107db)
or nous avons trouvé que
donc la plus grande équation sera 
LXX
L’équation

trouvée dans l’Article LXV et d’où nous avons tiré

à très-peu près,
est une suite du rapport que nous avons établi entre les révolutions synodiques des trois premiers satellites ; ce rapport n’est pas exact à la rigueur, mais il ne s’écarte de la vérité que d’une quantité infiniment petite, de sorte qu’au bout de
ans l’erreur qui en pourra résulter sera encore presque insensible.
En effet, on trouve qu’il faudrait environ
ans pour que l’équation dont nous parlons devînt

pourvu que les moyens mouvements des satellites fussent assez exacts pour pouvoir être employés dans une si longue suite de siècles. (Voyez l’Ouvrage de M. Wargentin cité ci-dessus.)
LXXI.
La formule de l’Article LVI, qui renferme les équations du troisième satellite, ne nous présente qu’un terme qui puisse être de quelque considération c’est le terme

dont l’argument est l’élongation du second satellite au troisième au temps des conjonctions de celui-ci.
Avant d’entrer dans le détail de l’inégalité qui en résulte, voyons si elle est assez considérable pour qu’on doive en tenir compte. Pour cela on substituera, au lieu de
sa valeur trouvée ci-dessus (Article LXIII), savoir
et l’on trouvera

de sorte que l’inégalité dont il s’agit montera à
à cause que l’équation
est tantôt additive, tantôt soustractive.
Maintenant on sait que les mouvements moyens du second et du troisième satellite sont entre eux comme les nombres
et
(Article LXV), d’où il suit que, pendant que le troisième achève une révolution au Soleil, la distance
augmente de
de sorte que, si l’on appelle
l’élongation du second satellite au troisième au temps d’une conjonction quelconque de ce dernier, on aura, après
révolutions,

et de là

d’où l’on voit : 1o que la période de cette équation sera de
révolutions du troisième satellite, ce qui revient au même que celles du premier et du second satellite (Articles LIX et LXVI) ; 2o que si l’on prend pour le commencement de la période une conjonction du troisième sa-
tellite dans laquelle

c’est-à-dire, que le second satellite soit aussi en conjonction, on trouvera que la marche de l’équation dont il s’agit sera entièrement analogue à celle de l’équation du second satellite (Article
LXVI).
LXXII.
L’équation que nous venons d’examiner ne se trouve point dans les Tables du troisième satellite ; M. Wargentin s’est contenté de l’indiquer dans la Préface de ses Tables (Mémoires de la Société d’Upsal, pour l’année 1741), où il dit Multæ etiam observationes satis manifeste indicant tertium æquatione alia indigere cujus fere eadem est quantitas et natura cum æquatione nova primi ; sed quoniam observationes non paucas habeam que eam vel minorem, vel nullam arguant, hujus æquationis in Tabulis nullam habere rationem satius judicavi ; et ailleurs (dans la Dissertation qui est à la tête des observations du second satellite) : In motibus tertii satellitis deprehenditur inæqualitas quædam quæ indicat eum esse retardatum in conjunctionibus, sed acceleratum in oppositionibus secundi ; et plus bas : Est etiam hæc inæqualitas tertii similis inaequalitati supra descriptæ secundi ; ce qui s’accorde parfaitement avec ce que nous avons trouvé dans l’Article précédent. Il est vrai que cette équation a paru à M. Wargentin de la même quantité que celle du premier satellite, au lieu qu’elle n’en est qu’environ les deux tiers, suivant notre Théorie ; mais ce savant Astronome avoue lui-même qu’il ne regarde pas son résultat comme fort exact, l’ayant trouvé quelquefois moindre, et même nul, et que c’est pour cette raison qu’il a cru devoir s’abstenir d’en faire usage dans ses Tables.
LXXIII.
Avant de quitter la formule de l’Article LVI, nous dirons deux mots des termes qui dépendent de
élongation du quatrième satellite au troisième, et dont le plus considérable est celui-ci

Supposons d’abord que l’équation qui en provient soit, lorsqu’elle est la plus grande, de
minutes ; on aura

donc

d’où l’on voit que, pour que
soit au moins
=1,
il faut que la masse du quatrième satellite surpasse de beaucoup celles des trois premiers.
Si l’on veut que la densité de ce satellite soit la même que celle de Jupiter, on trouvera son diamètre
de celui de Jupiter ; et par conséquent le temps qu’il doit employer à entrer dans l’ombre
ce qui, en faisant
est assez conforme au résultat des observations de M. Maraldi.
Cette équation, au reste, supposé qu’elle montât à quelques minutes, ce qui ne serait nullement impossible, mériterait d’autant plus l’attention des Astronomes qu’elle varie beaucoup d’une conjonction à l’autre ; en effet, les révolutions synodiques du troisième et du quatrième satellite étant de
et
on trouve que l’angle
doit augmenter pendant une révolution du troisième de
c’est-à-dire, diminuer de
donc, nommant
l’angle
dans le temps d’une conjonction quelconque de ce satellite, on aura, après
révolutions,

d’où

par conséquent la période de cette équation ne sera que de
révolutions, c’est-à-dire, de
révolutions, ce qui fait
à peu
près. Ne serait-ce point là la source de ces inégalités qu’on observe dans les conjonctions du troisième satellite, et qui font des sauts considérables d’une conjonction à l’autre ? C’est une vue que nous proposons aux Astronomes qui s’occupent de la Théorie des satellites.
LXXIV.
Il ne resterait plus qu’à examiner les équations des conjonctions du quatrième satellite, contenues dans la formule de l’Article LVIII ; mais ayant déjà trouvé (Articles LXIII et LXVIII)

et

on verra aisément que les coefficients de-ces équations ne s’étendent point au delà d’un petit nombre de secondes ; ce qui est trop peu de chose pour qu’on doive en tenir compte dans le mouvement de ce satellite, surtout vu l’imperfection qui règne encore dans les Tables de Jupiter.
CHAPITRE IV.
SUITE DU CALCUL DES PERTURBATIONS DES SATELLITES DE JUPITER.
LXXV.
Ayant trouvé les premières valeurs de
(Articles XXXVI et suivants), on reprendra les équations de l’Article XXXII ; et après y avoir mis au lieu de
leurs expressions (Articles XXX et suivants), sans négliger les termes de l’ordre
on substituera dans tous les termes de cet ordre les valeurs trouvées de
et l’on aura de nouvelles équations en
plus exactes que celles de l’Article XXXV, et qui s’intégreront encore par la même méthode.
LXXVI.
Soit pour le premier satellite (ces formules s’appliquent également aux trois autres, suivant les remarques des Articles IX et XIII)
![{\displaystyle {\begin{aligned}\mathrm {L} _{1}\ =&g_{1}-2\mu _{1}\mathrm {H} _{1}-{\frac {\chi _{2}}{1+n\chi _{1}}}{\breve {\Gamma }}(a_{1},a_{2})-{\frac {\chi _{3}}{1+n\chi _{1}}}{\breve {\Gamma }}(a_{1},a_{3})-{\frac {\chi _{4}}{1+n\chi _{1}}}{\breve {\Gamma }}(a_{1},a_{4})\\&-{\frac {1}{2}}{\frac {\mathrm {K} _{1}}{1+n\chi _{1}}}+{\frac {1}{5}}{\frac {\varkappa _{1}}{1+n\chi _{1}}},\\\mathrm {M} _{1}^{2}=&3\mu _{1}^{2}-2f_{1}-nf_{1}\left[\chi _{2}{\breve {\Pi }}(a_{1},a_{2})+\chi _{3}{\breve {\Pi }}(a_{1},a_{3})+\chi _{4}{\breve {\Pi }}(a_{1},a_{4})+{\frac {1}{2}}\mathrm {K} _{1}+{\frac {1}{5}}\varkappa _{1}\right],\\\mathrm {N} _{1}^{2}\,=&\mu _{1}^{2}\\+&nf_{1}\left[\chi _{2}{\overset {\backsim }{\Gamma }}(a_{1},a_{2})+\chi _{3}{\overset {\backsim }{\Gamma }}(a_{1},a_{3})+\chi _{4}{\overset {\backsim }{\Gamma }}(a_{1},a_{4})+{\frac {3}{2}}\mathrm {K} _{1}+{\frac {2}{5}}\varkappa _{1}\right]+2n\mu _{1}f_{1}\mathrm {H} _{1}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19467210fa56856575b9089aae3887543eba4d53)
Supposons de plus (Articles XXXVI et XXXVIII)
![{\displaystyle {\begin{aligned}\theta _{1}=&-{\frac {\chi _{2}}{1+n\chi _{1}}}\left[\Xi _{1}(a_{1},a_{2})\cos(u_{2}-u_{1})t+\Xi _{2}(a_{1},a_{2})\cos 2(u_{2}-u_{1})t+\ldots \right]\\&-{\frac {\chi _{3}}{1+n\chi _{1}}}\left[\Xi _{1}(a_{1},a_{3})\cos(u_{3}-u_{1})t+\Xi _{2}(a_{1},a_{3})\cos 2(u_{3}-u_{1})t+\ldots \right]\\&-{\frac {\chi _{4}}{1+n\chi _{1}}}\left[\Xi _{1}(a_{1},a_{4})\cos(u_{4}-u_{1})t+\Xi _{2}(a_{1},a_{4})\cos 2(u_{4}-u_{1})t+\ldots \right]\\&-{\frac {\mathrm {K} _{1}}{1+n\chi _{1}}}\gamma _{1}\sin 2(m-\mu _{1})t,\\\\\vartheta _{1}=&{\frac {\chi _{2}}{1+n\chi _{1}}}\left[\Phi _{1}(a_{1},a_{2})\sin(u_{2}-u_{1})t+\Phi _{2}(a_{1},a_{2})\sin 2(u_{2}-u_{1})t+\ldots \right]\\+&{\frac {\chi _{3}}{1+n\chi _{1}}}\left[\Phi _{1}(a_{1},a_{3})\sin(u_{3}-u_{1})t+\Phi _{2}(a_{1},a_{3})\sin 2(u_{3}-u_{1})t+\ldots \right]\\+&{\frac {\chi _{4}}{1+n\chi _{1}}}\left[\Phi _{1}(a_{1},a_{4})\sin(u_{4}-u_{1})t+\Phi _{2}(a_{1},a_{4})\sin 2(u_{4}-u_{1})t+\ldots \right]\\+&{\frac {\mathrm {K} _{1}}{1+n\chi _{1}}}\gamma _{1}\sin 2(m-\mu _{1})t,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78349be95eebb4ac7582061c3d05b532638e4cc1)
et faisons
(nous verrons bientôt la raison de ces substitutions). Les équations de l’Article XXXIII se changeront en
celles-ci, dans lesquelles nous avons négligé les termes affectés de
![{\displaystyle {\begin{aligned}(\mathrm {G} )&\ \ {\frac {d^{2}x_{1}}{dt^{2}}}+\mathrm {M} _{1}^{2}x_{1}+f_{1}\mathrm {L} _{1}-n\left(6\mu _{1}^{2}-3f_{1}\right)x_{1}^{2}-{\frac {1}{2}}nf_{1}z_{1}^{2}-nf_{1}^{2}\mathrm {Y} _{1}^{2}\\&-n\chi _{2}f_{1}x_{1}\left[{\breve {\Pi }}_{1}(a_{1},a_{2})+{\frac {2\mu _{1}}{\mu _{2}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\right]\cos(\mu _{2}-\mu _{1})t\\&-n\chi _{2}f_{1}x_{1}\left[{\breve {\Pi }}_{2}(a_{1},a_{2})+{\frac {2\mu _{1}}{2(\mu _{2}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{2})\right]\cos 2(\mu _{2}-\mu _{1})t\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&-n\chi _{3}f_{1}x_{1}\left[{\breve {\Pi }}_{1}(a_{1},a_{3})+{\frac {2\mu _{1}}{\mu _{3}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\right]\cos(\mu _{3}-\mu _{1})t\\&-n\chi _{3}f_{1}x_{1}\left[{\breve {\Pi }}_{2}(a_{1},a_{3})+{\frac {2\mu _{1}}{2(\mu _{3}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{3})\right]\cos 2(\mu _{3}-\mu _{1})t\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&-n\chi _{4}f_{1}x_{1}\left[{\breve {\Pi }}_{1}(a_{1},a_{4})+{\frac {2\mu _{1}}{\mu _{4}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{4})\right]\cos(\mu _{4}-\mu _{1})t\\&-n\chi _{4}f_{1}x_{1}\left[{\breve {\Pi }}_{2}(a_{1},a_{4})+{\frac {2\mu _{1}}{2(\mu _{4}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{4})\right]\cos 2(\mu _{4}-\mu _{1})t\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\\\&-n\mathrm {K} _{1}f_{1}x_{1}\left[{\frac {3}{2}}-{\frac {3\mu _{1}}{2(m-\mu _{1})}}\cos 2(m-\mu _{1})t\right]\\&-n\chi _{2}f_{1}x_{2}\\&\quad \times \left[{\breve {\Psi }}(a_{1},a_{2})+{\breve {\Psi }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+{\breve {\Psi }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&-n\chi _{3}f_{1}x_{3}\\&\quad \times \left[{\breve {\Psi }}(a_{1},a_{3})+{\breve {\Psi }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+{\breve {\Psi }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&-n\chi _{4}f_{1}x_{4}\\&\quad \times \left[{\breve {\Psi }}(a_{1},a_{4})+{\breve {\Psi }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+{\breve {\Psi }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\\&+n\mathrm {K} _{1}f_{1}\xi \left[{\frac {3}{2}}+{\frac {9}{2}}\cos 2(m-\mu _{1})t\right]\\&+n\chi _{2}f_{1}(y_{2}-y_{1})\\&\quad \times \left[{\breve {\Gamma }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t-\ldots \right]\\&+n\chi _{3}f_{1}(y_{3}-y_{1})\\&\quad \times \left[{\breve {\Gamma }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t-\ldots \right]\\&+n\chi _{4}f_{1}(y_{4}-y_{1})\\&\quad \times \left[{\breve {\Gamma }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+2{\breve {\Gamma }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t-\ldots \right]\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/428f8a5098d6a7326458a08a0a12e2bb9cdc4d1d)
![{\displaystyle {\begin{aligned}&+n\mathrm {K} _{1}f_{1}(\mathrm {J} -y_{1})\times 3\sin 2(m-\mu _{1})t\\&+2n\chi _{2}f_{1}\mu _{1}\\&\quad \times \int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\&+2n\chi _{3}f_{1}\mu _{1}\\&\quad \times \int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\&+2n\chi _{4}f_{1}\mu _{1}\\&\quad \times \int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\\\&+2n\mathrm {K} _{1}f_{1}\mu _{1}\int x_{1}\times 3\sin 2(m-\mu _{1})tdt\\&+2n\chi _{2}f_{1}\mu _{1}\\&\quad \times \int x_{2}\left[{\widehat {\Psi }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\&+2n\chi _{3}f_{1}\mu _{1}\\&\quad \times \int x_{3}\left[{\widehat {\Psi }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\&+2n\chi _{4}f_{1}\mu _{1}\\&\quad \times \int x_{4}\left[{\widehat {\Psi }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\\\&+2n\mathrm {K} _{1}f_{1}\mu _{1}\int \xi \times {\frac {9}{2}}\sin 2(m-\mu _{1})tdt\\+&2n\chi _{2}f_{1}\mu _{1}\\&\ \ \times \int (y_{2}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\+&2n\chi _{3}f_{1}\mu _{1}\\&\ \ \times \int (y_{3}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\+&2n\chi _{4}f_{1}\mu _{1}\\&\ \ \times \int (y_{4}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\\\&-2n\mathrm {K} _{1}f_{1}x_{1}\int (\mathrm {J} -y_{1})\times 3\cos 2(m-\mu _{1})tdt\\(\mathrm {H} )&\ \ {\frac {dy_{1}}{dt}}+2\mu _{1}x_{1}-f_{1}\mathrm {H} _{1}-3n\mu x_{1}^{2}\\&+2n\chi _{2}f_{1}x_{1}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{2})}{\mu _{2}-\mu _{1}}}\cos(\mu _{2}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{2})}{2(\mu _{2}-\mu _{1})}}\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]\\&+2n\chi _{3}f_{1}x_{1}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{3})}{\mu _{3}-\mu _{1}}}\cos(\mu _{3}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{3})}{2(\mu _{3}-\mu _{1})}}\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]\\&+2n\chi _{4}f_{1}x_{1}\left[{\frac {{\widehat {\Gamma }}_{1}(a_{1},a_{4})}{\mu _{4}-\mu _{1}}}\cos(\mu _{4}-\mu _{1})t+{\frac {{\widehat {\Gamma }}_{2}(a_{1},a_{4})}{2(\mu _{4}-\mu _{1})}}\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4c7b674932f7df9e91ae10e02ed886ee088f6581)
![{\displaystyle {\begin{aligned}&-2n\mathrm {K} _{1}f_{1}x_{1}\times {\frac {3}{4(m-\mu _{1})}}\cos 2(m-\mu _{1})t\\&+n\chi _{2}f_{1}\\&\quad \times \int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\&+n\chi _{3}f_{1}\\&\quad \times \int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\&+n\chi _{4}f_{1}\\&\quad \times \int x_{1}\left[{\overset {\backsim }{\Pi }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+{\overset {\backsim }{\Pi }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\\\&+n\mathrm {K} _{1}f_{1}\int x_{1}\times 3\sin 2(m-\mu _{1})tdt\\&+n\chi _{2}f_{1}\\&\quad \times \int x_{2}\left[{\widehat {\Psi }}_{1}(a_{1},a_{2})\sin(\mu _{2}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{2})\sin 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\&+n\chi _{3}f_{1}\\&\quad \times \int x_{3}\left[{\widehat {\Psi }}_{1}(a_{1},a_{3})\sin(\mu _{3}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{3})\sin 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\&+n\chi _{4}f_{1}\\&\quad \times \int x_{4}\left[{\widehat {\Psi }}_{1}(a_{1},a_{4})\sin(\mu _{4}-\mu _{1})t+{\widehat {\Psi }}_{2}(a_{1},a_{4})\sin 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\\\&+n\mathrm {K} _{1}f_{1}\int \xi \times {\frac {9}{2}}\sin 2(m-\mu _{1})tdt\\+&n\chi _{2}f_{1}\\&\ \ \times \int (y_{2}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t+\ldots \right]dt\\+&n\chi _{3}f_{1}\\&\ \ \times \int (y_{3}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t+\ldots \right]dt\\+&n\chi _{4}f_{1}\\&\ \ \times \int (y_{4}-y_{1})\left[{\widehat {\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+2{\widehat {\Gamma }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t+\ldots \right]dt\\\\&-n\mathrm {K} _{1}f_{1}\int (\mathrm {J} -y_{1})\times 3\cos 2(m-\mu _{1})tdt=0,\\\\(\mathrm {K} )&\ \ {\frac {d^{2}z_{1}}{dt^{2}}}+\mathrm {N} _{1}^{2}-4n\mu _{1}^{2}z_{1}x_{1}+2n{\frac {dx_{1}dz_{1}}{dt^{2}}}\\&+n\chi _{2}f_{1}z_{1}\left[{\overset {\backsim }{\Gamma }}_{1}(a_{1},a_{2})+{\frac {2\mu _{1}}{\mu _{2}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\right]\cos(\mu _{2}-\mu _{1})t\\&+n\chi _{2}f_{1}z_{1}\left[{\overset {\backsim }{\Gamma }}_{2}(a_{1},a_{2})+{\frac {2\mu _{1}}{2(\mu _{2}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{2})\right]\cos 2(\mu _{2}-\mu _{1})t\ldots \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de9b6ba2b57fb9f1f76f0db970cf322628040386)
![{\displaystyle {\begin{aligned}&+n\chi _{3}f_{1}z_{1}\left[{\overset {\backsim }{\Gamma }}_{1}(a_{1},a_{3})+{\frac {2\mu _{1}}{\mu _{3}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\right]\cos(\mu _{3}-\mu _{1})t\\&+n\chi _{3}f_{1}z_{1}\left[{\overset {\backsim }{\Gamma }}_{2}(a_{1},a_{3})+{\frac {2\mu _{1}}{2(\mu _{3}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{3})\right]\cos 2(\mu _{3}-\mu _{1})t\ldots \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&+n\chi _{4}f_{1}z_{1}\left[{\overset {\backsim }{\Gamma }}_{1}(a_{1},a_{4})+{\frac {2\mu _{1}}{\mu _{4}-\mu _{1}}}{\widehat {\Gamma }}_{1}(a_{1},a_{4})\right]\cos(\mu _{4}-\mu _{1})t\\&+n\chi _{4}f_{1}z_{1}\left[{\overset {\backsim }{\Gamma }}_{2}(a_{1},a_{4})+{\frac {2\mu _{1}}{2(\mu _{4}-\mu _{1})}}{\widehat {\Gamma }}_{2}(a_{1},a_{4})\right]\cos 2(\mu _{4}-\mu _{1})t\ldots \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&+n\mathrm {K} _{1}f_{1}z_{1}\left[{\frac {3}{2}}-{\frac {3\mu _{1}}{2(m-\mu _{1})}}\right]\cos 2(m-\mu _{1})t\\-&n\chi _{2}f_{1}z_{2}\\&\ \ \times \left[{\overset {\circ }{\Gamma }}(a_{1},a_{2})+{\overset {\circ }{\Gamma }}_{1}(a_{1},a_{2})\cos(\mu _{2}-\mu _{1})t+{\overset {\circ }{\Gamma }}_{2}(a_{1},a_{2})\cos 2(\mu _{2}-\mu _{1})t\ldots \right]\\-&n\chi _{3}f_{1}z_{3}\\&\ \ \times \left[{\overset {\circ }{\Gamma }}(a_{1},a_{3})+{\overset {\circ }{\Gamma }}_{1}(a_{1},a_{3})\cos(\mu _{3}-\mu _{1})t+{\overset {\circ }{\Gamma }}_{2}(a_{1},a_{3})\cos 2(\mu _{3}-\mu _{1})t\ldots \right]\\-&n\chi _{4}f_{1}z_{4}\\&\ \ \times \left[{\overset {\circ }{\Gamma }}(a_{1},a_{4})+{\overset {\circ }{\Gamma }}_{1}(a_{1},a_{4})\cos(\mu _{4}-\mu _{1})t+{\overset {\circ }{\Gamma }}_{2}(a_{1},a_{4})\cos 2(\mu _{4}-\mu _{1})t\ldots \right]=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a1d88b94371888b98f686748553320a5b77ab30e)
LXXVII.
Si l’on rejette dans les équations
et
tous les termes affectés de
comme aussi tous les termes constants qui doivent être nuls par les conditions de l’Article XXXII on a

D’où l’on tire

ce qui donne pour
et
les mêmes valeurs que nous avons déjà trouvées (Article LXXV).
La quantité
n’est que le premier terme de l’équation du centre calculée dans une ellipse mobile (Article XXXVIII) ; si l’on voulait avoir le terme suivant, c’est-à-dire celui qui contient le carré de l’excentricité, il n’y aurait qu’à mettre au lieu de
dans les termes
et
des équations
la valeur de
qu’on vient de trouver.
On aurait donc, en négligeant toujours les termes constants,

La première de ces équations donne (Article XXXIV)

c’est-à-dire, en mettant au lieu de
et de
leurs valeurs approchées
(Article XLV),

Donc, substituant cette valeur de
dans la seconde, et l’intégrant ensuite, on aura

Ce qui s’accorde avec ce que l’on sait d’ailleurs ; mais nous verrons plus bas qu’il y a dans l’équation
d’autres termes qui influent considérablement sur l’équation du centre, et qui empêchent qu’on ne puisse regarder l’expression précédente comme assez exacte, même dans le cas où l’on néglige les quantités de l’ordre
.
Il en faut dire autant de l’expression de la latitude que nous avons déjà trouvée (Article XL) ; mais avant que d’entrer dans cette discussion, il est bon de voir ce que donnent les nouvelles valeurs de
et de
(Article précédent) pour le mouvement des apsides et des nœuds.
§ I. — Premières valeurs du mouvement des apsides et des nœuds des satellites.
LXXVIII.
Nous avons trouvé (Article LXXVI)

![{\displaystyle -nf_{1}\left[\chi _{2}{\breve {\Pi }}(a_{1},a_{2})+\chi _{3}{\breve {\Pi }}(a_{1},a_{3})+\chi _{4}{\breve {\Pi }}(a_{1},a_{4})+{\frac {1}{2}}\mathrm {K} _{1}+{\frac {1}{5}}\varkappa _{1}\right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1bfa6bed0211342415554e018e69fba68bc1f232)
or, on a généralement (Article XXIX)

d’où

et par conséquent

donc, négligeant les termes affectés de
on aura
![{\displaystyle \mathrm {M} _{1}^{2}=\mu _{1}^{2}\left[1-ng_{1}-n\chi _{2}{\breve {\Pi }}(a_{1},a_{2})-n\chi _{3}{\breve {\Pi }}(a_{1},a_{3})-n\chi _{4}{\breve {\Pi }}(a_{1},a_{4})-{\frac {1}{2}}\mathrm {K} _{1}-{\frac {4}{5}}\varkappa _{1}\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aea2c4b1487ed1a45329824b1f88554ff44e28d7)
Maintenant on a par l’équation
de l’Article LXXVI

![{\displaystyle +{\frac {1}{1+n\chi _{1}}}\left[\chi _{2}{\breve {\Gamma }}(a_{1},a_{2})+\chi _{3}{\breve {\Gamma }}(a_{1},a_{3})+\chi _{4}{\breve {\Gamma }}(a_{1},a_{4})+{\frac {1}{2}}\mathrm {K} _{1}-{\frac {1}{5}}\varkappa _{1}\right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a28c7b06be6e6b7f5bc95f0ecff9bd00110eb3bd)
donc
![{\displaystyle {\begin{aligned}\mathrm {M} _{1}^{2}=\mu _{1}^{2}&{\biggl [}1-n\chi _{2}\left[{\breve {\Pi }}(a_{1},a_{2})+2{\breve {\Gamma }}(a_{1},a_{2})\right]{\biggr .}\\&-n\chi _{3}\left[{\breve {\Pi }}(a_{1},a_{3})+2{\breve {\Gamma }}(a_{1},a_{3})\right]-n\chi _{4}\left[{\breve {\Pi }}(a_{1},a_{4})+2{\breve {\Gamma }}(a_{1},a_{4})\right]\\&{\biggl .}-{\frac {3}{2}}n\mathrm {K} _{1}-{\frac {2}{5}}n\varkappa _{1}-2n\mathrm {L} _{1}-4\mu _{1}\mathrm {H} _{1}{\biggr ]},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/705867c9e505453af778ea09238dd3dbc5d0de1c)
les quantités
devant être déterminées par la condition que les équations
ne renferment aucun terme constant.
Pour remplir ces deux conditions, on supposera que
soit la quantité constante qui entre dans la valeur de
car la valeur de
étant composée de sinus et de cosinus, il est évident que le carré
contiendra nécessairement des termes constants, quoique
n’en contienne point ; de même, soient
les quantités constantes qui entreront dans les valeurs de
on aura donc


d’où l’on tirera
et
qui seront de l’ordre de
c’est pourquoi on peut négliger dans la valeur de
les quantités
et
qui seraient de l’ordre 
LXXIX.
Donc, si l’on fait, pour abréger,
![{\displaystyle {\begin{aligned}{\text{ϐ}}_{1}&=\chi _{2}\left[{\breve {\Pi }}(a_{1},a_{2})+2{\breve {\Gamma }}(a_{1},a_{2})\right]+\chi _{3}\left[{\breve {\Pi }}(a_{1},a_{3})+2{\breve {\Gamma }}(a_{1},a_{3})\right]\\&+\chi _{4}\left[{\breve {\Pi }}(a_{1},a_{4})+2{\breve {\Gamma }}(a_{1},a_{4})\right]+{\frac {3}{2}}\mathrm {K} _{1}+{\frac {2}{5}}\varkappa _{1},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f6712af6b0fd544cc748be74a8937a1876dda35e)
et de même (Article IX)
![{\displaystyle {\begin{aligned}{\text{ϐ}}_{2}&=\chi _{1}\left[{\breve {\Pi }}(a_{2},a_{1})+2{\breve {\Gamma }}(a_{2},a_{1})\right]+\chi _{3}\left[{\breve {\Pi }}(a_{2},a_{3})+2{\breve {\Gamma }}(a_{2},a_{3})\right]\\&+\chi _{4}\left[{\breve {\Pi }}(a_{2},a_{4})+2{\breve {\Gamma }}(a_{2},a_{4})\right]+{\frac {3}{2}}\mathrm {K} _{2}+{\frac {2}{5}}\varkappa _{2},\\\\{\text{ϐ}}_{3}&=\chi _{1}\left[{\breve {\Pi }}(a_{3},a_{1})+2{\breve {\Gamma }}(a_{3},a_{1})\right]+\chi _{2}\left[{\breve {\Pi }}(a_{3},a_{2})+2{\breve {\Gamma }}(a_{3},a_{2})\right]\\&+\chi _{4}\left[{\breve {\Pi }}(a_{3},a_{4})+2{\breve {\Gamma }}(a_{3},a_{4})\right]+{\frac {3}{2}}\mathrm {K} _{3}+{\frac {2}{5}}\varkappa _{3},\\\\{\text{ϐ}}_{4}&=\chi _{1}\left[{\breve {\Pi }}(a_{4},a_{1})+2{\breve {\Gamma }}(a_{4},a_{1})\right]+\chi _{2}\left[{\breve {\Pi }}(a_{4},a_{2})+2{\breve {\Gamma }}(a_{4},a_{2})\right]\\&+\chi _{3}\left[{\breve {\Pi }}(a_{4},a_{3})+2{\breve {\Gamma }}(a_{4},a_{3})\right]+{\frac {3}{2}}\mathrm {K} _{4}+{\frac {2}{5}}\varkappa _{4},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1edff8cbd47eff062631ffe1a2ee9a7aac70b019)
on aura

et ainsi des autres.
Or le mouvement de la ligne des apsides étant au mouvement moyen comme
à
(Article XXXVII), cette ligne avancera pendant une révolution de
d’où l’on connaîtra le mouvement des apsides de tous les satellites.
LXXX.
Pour évaluer en nombres les quantités
il faut commencer par chercher les valeurs des quantités
lesquelles dépendent des quantités
c’est-à-dire des coefficients de la série qui représente la quantité radicale
(Articles XX et suivants).
Soit donc, comme dans cet Article,

on aura, en faisant
et 
![{\displaystyle \left[a_{1}^{2}-2a_{1}a_{2}\cos(\varphi _{2}-\varphi _{1})+a_{2}^{2}\right]^{-{\frac {5}{2}}}=}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ff98c2879fb83b4bdf1c3e6ecd8b8f5d627dac2)
![{\displaystyle a_{2}^{-5}\left[(\mathrm {A} )+(\mathrm {B} )\cos(\varphi _{2}-\varphi _{1})+(\mathrm {C} )\cos 2(\varphi _{2}-\varphi _{1})+\ldots \right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3e914f42db0b0082704e715c2485432e94e75112)
donc(Article XXI)

De même, en faisant
on aura

et ainsi de suite ; où l’on remarquera que les quantités réciproques
sont les mêmes que les quantités
(Article XLIII).
Cela posé, on trouvera (Article XXII),
étant égal à

et de là, en faisant, pour abréger,

on aura par l’Article XXIV

On trouvera des expressions semblables pour les fonctions
de 
en faisant
De la même manière on trouvera que l’on a,
étant encore égal à

d’où, en faisant

on tire

expressions qui serviront aussi pour les quantités
en faisant successivement 
LXXXI.
Ayant donc fait le calcul de ces différentes quantités, j’ai trouvé les valeurs suivantes :



LXXXII.
À l’égard des valeurs de
nous les avons données ci-dessus (Article XLVII) aussi bien que celles de
(Article XLIX) ; et pour ce qui est des quantités
(Article XXIX) on aura, en faisant
demi-diamètre de Jupiter, égal à
et mettant pour
leurs valeurs (Article XLVI), on aura, dis-je,

la quantité
dépendant de la figure et de la constitution intérieure de Jupiter, comme on l’a vu (Article XVI).
LXXXIII.
Toutes ces substitutions faites, on aura, après avoir remis au lieu de
les quantités

LXXXIV.
Passons maintenant aux formules qui donnent le mouvement des nœuds, et nous trouvons d’abord pour le premier satellite (Article LXXV)
![{\displaystyle \mathrm {N} _{1}^{2}=\mu _{1}^{2}+nf_{1}\left[\chi _{2}{\overset {\backsim }{\Gamma }}(a_{1},a_{2})+\chi _{3}{\overset {\backsim }{\Gamma }}(a_{1},a_{3})+\chi _{4}{\overset {\backsim }{\Gamma }}(a_{1},a_{4})+{\frac {3}{2}}\mathrm {K} _{1}+{\frac {3}{5}}\varkappa _{1}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/405b0406ea870f985fcc1c3ecb30c75a901fda28)

c’est-à-dire, en mettant
au lieu de
et négligeant le terme 
qui est du second ordre, à cause que

est déjà de l’ordre de
(Article LXXVIII),
![{\displaystyle \mathrm {N} _{1}^{2}=\mu _{1}^{2}\left[1+n\chi _{2}{\overset {\backsim }{\Gamma }}(a_{1},a_{2})+n\chi _{3}{\overset {\backsim }{\Gamma }}(a_{1},a_{3})+n\chi _{4}{\overset {\backsim }{\Gamma }}(a_{1},a_{4})+{\frac {3}{2}}n\mathrm {K} _{1}+{\frac {2}{5}}n\varkappa _{1}\right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/121a57f92c0219a1d8f358dfacfc1d3a2b7207d8)
donc si l’on fait, pour abréger,

de même

On aura pour tous les quatre satellites

Pour tirer de là le mouvement des nœuds, on remarquera que
exprime, en général, la longitude moyenne du nœud ascendant (Article XLI) ; d’où il suit que le mouvement de la ligne des nœuds sera au mouvement moyen comme
à
c’est-à-dire, comme
à
par conséquent les nœuds reculeront à chaque révolution de
LXXXV.
Or, on trouvé par l’Article XXXII en faisant successivement 

ensuite,

ce qui donne (Article XLVII)


LXXXVI.
Donc, faisant ces substitutions, et remettant
au lieu de
on aura

LXXXVII.
Nous avons trouvé (Article LXXVII)

on trouvera de même

et ainsi des autres. Cela posé, si l’on reprend l’équation
de l’Article LXXVI, et qu’on substitue dans le terme

au lieu de
sa valeur
on verra que la quantité
renfermera un terme de cette forme
[8] ;
lequel étant intégré deviendra
![{\displaystyle {\frac {\cos \left[\left(\mu _{1}-{\cfrac {n{\text{ϐ}}_{2}}{2}}\mu _{2}\right)t+\omega _{2}\right]}{\left(\mu _{1}-{\cfrac {n{\text{ϐ}}_{2}}{2}}\mu _{2}\right)^{2}-\left(\mu _{1}-{\cfrac {n{\text{ϐ}}_{1}}{2}}\mu _{1}\right)^{2}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0c43d80870e31bd485a75a6da4c20f9daad13f1c)
ainsi le terme

de l’équation différentielle donnera dans la valeur de

le terme suivant
![{\displaystyle {\frac {\varepsilon _{2}}{2}}{\frac {n\chi _{2}f_{1}{\breve {\Psi }}_{1}(a_{1},a_{2})\cos \left[\left(\mu _{1}-{\frac {n{\text{ϐ}}_{2}}{2}}\mu _{2}\right)t+\omega _{2}\right]}{\left(\mu _{1}-{\cfrac {n{\text{ϐ}}_{2}}{2}}\mu _{2}\right)^{2}-\left(\mu _{1}-{\cfrac {n{\text{ϐ}}_{1}}{2}}\mu _{1}\right)^{2}}}=}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95f3522326b164543ac99c7cd94064c1ba6c9a7a)
![{\displaystyle {\frac {\varepsilon }{2}}{\frac {\chi _{2}f_{1}{\breve {\Psi }}_{1}(a_{1},a_{2})}{\mu _{1}({\text{ϐ}}_{1}\mu _{1}-{\text{ϐ}}_{2}\mu _{2})}}\cos \left[\left(\mu _{1}-{\frac {n{\text{ϐ}}_{2}}{2}}\mu _{2}\right)t+\omega _{2}\right],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d2876df31f2a9539bd859e2a8ebbefbb8d5ac5c6)
lequel appartient, comme on voit, à la première valeur de
Pareillement le terme

donnera dans la valeur de
le terme
![{\displaystyle {\frac {\varepsilon _{3}}{2}}{\frac {\chi _{3}f_{1}{\breve {\Psi }}_{1}(a_{1},a_{3})}{\mu _{1}({\text{ϐ}}_{1}\mu _{1}-{\text{ϐ}}_{3}\mu _{3})}}\cos \left[\left(\mu _{1}-{\frac {n{\text{ϐ}}_{3}}{2}}\mu _{3}\right)t+\omega _{3}\right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f8b7178ef63c82c877bda95a215aaf07940a8639)
et il en sera de même de quelques autres termes de l’équation
dont nous parlerons plus bas.
On trouvera de la même manière dans la valeur de
les termes
![{\displaystyle {\begin{aligned}{\frac {\varepsilon _{1}}{2}}{\frac {\chi _{1}f_{2}{\breve {\Psi }}_{1}(a_{2},a_{1})}{\mu _{2}({\text{ϐ}}_{2}\mu _{2}-{\text{ϐ}}_{1}\mu _{1})}}\cos \left[\left(\mu _{2}-{\frac {n{\text{ϐ}}_{1}}{2}}\mu _{1}\right)t+\omega _{1}\right],\\{\frac {\varepsilon _{3}}{2}}{\frac {\chi _{3}f_{2}{\breve {\Psi }}_{1}(a_{2},a_{3})}{\mu _{2}({\text{ϐ}}_{2}\mu _{2}-{\text{ϐ}}_{3}\mu _{3})}}\cos \left[\left(\mu _{2}-{\frac {n{\text{ϐ}}_{3}}{2}}\mu _{3}\right)t+\omega _{3}\right],\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3bf8352128573194e8e7d2d61617712aab407496)
lesquels étant de nouveau substitués dans le terme

de l’équation
en donneront deux autres de cette forme
![{\displaystyle {\begin{aligned}&{\frac {1}{2}}n\chi _{2}f_{1}{\breve {\Psi }}_{1}(a_{1},a_{2}){\frac {\varepsilon _{1}}{2}}{\frac {\chi _{1}f_{2}{\breve {\Psi }}_{1}(a_{2},a_{1})}{\mu _{2}({\text{ϐ}}_{2}\mu _{2}-{\text{ϐ}}_{1}\mu _{1})}}\cos \left[\left(\mu _{1}-{\frac {n{\text{ϐ}}_{1}}{2}}\mu _{1}\right)t+\omega _{1}\right],\\&{\frac {1}{2}}n\chi _{2}f_{1}{\breve {\Psi }}_{1}(a_{1},a_{2}){\frac {\varepsilon _{3}}{2}}{\frac {\chi _{3}f_{2}{\breve {\Psi }}_{1}(a_{2},a_{3})}{\mu _{2}({\text{ϐ}}_{2}\mu _{2}-{\text{ϐ}}_{3}\mu _{3})}}\cos \left[\left(\mu _{1}-{\frac {n{\text{ϐ}}_{3}}{2}}\mu _{3}\right)t+\omega _{3}\right]\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/74a55211fca89990c6e6a0da00c778bc39579bd5)
le premier de ces deux termes produira

à cause des
![{\displaystyle {\Biggl .}\mu \left(1-{\frac {n{\text{ϐ}}}{2}}\right)=\mathrm {M} {\Biggr ]}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/29917fbb475afd5248f6e5129590b6f737caa758)
dans la valeur de

un terme qui sera multiplié par l’angle

(Article
XXXIV) ; ce qui donnera des arcs de cercle dans le rayon vecteur de l’orbite ; le second y produira le terme
![{\displaystyle {\frac {1}{2}}{\frac {\chi _{2}f_{1}{\breve {\Psi }}_{1}(a_{1},a_{2})}{\mu _{2}({\text{ϐ}}_{3}\mu _{3}-{\text{ϐ}}_{1}\mu _{1})}}{\frac {\varepsilon _{3}}{2}}{\frac {\chi _{3}f_{2}{\breve {\Psi }}_{1}(a_{2},a_{3})}{\mu _{2}({\text{ϐ}}_{2}\mu _{2}-{\text{ϐ}}_{3}\mu _{3})}}\cos \left[\left(\mu _{1}-{\frac {n{\text{ϐ}}_{3}}{2}}\mu _{3}\right)t+\omega _{3}\right],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4a3c29e4eed90d394aa835760998cd4dcb173a9d)
qui est de la même forme que celui que nous avons déjà trouvé.
Ces termes en reproduiront d’autres dans la valeur de
de la même forme que ceux que nous venons d’examiner, d’où il renaîtra encore dans la valeur de
d’autres termes de même espèce que les précédents, et ainsi de suite à l’infini.
LXXXVIII.
De là je tire ces deux conséquences fort importantes : 1o que les termes dont il s’agit, quoique de l’ordre
dans l’équation différentielle, appartiennent cependant à la première approximation et ne doivent point être négligés dans les premières valeurs de
2o que la méthode ordinaire d’approximation, suivant laquelle on emploie à chaque nouvelle correction les valeurs trouvées dans la correction précédente, est absolument insuffisante pour calculer ces sortes de termes.
On appliquera le même raisonnement à l’équation
et l’on en tirera des conclusions analogues par rapport à la valeur de
.
LXXXIX.
Il est donc nécessaire d’avoir une méthode particulière pour intégrer les équations
on verra dans le paragraphe suivant comment je m’y suis pris pour arriver à ce but ; mais il faut commencer ici par voir quels sont les termes de ces équations, auxquels on doit avoir égard.
Pour peu qu’on examine l’équation (
), on reconnaîtra aisément que les termes dont il s’agit viennent uniquement des termes qui renferment


en tant qu’on y substitue
à la place de
de sorte qu’on pourra réduire cette équation à celle-ci

À l’égard de l’équation (
), on trouvera qu’elle se réduit de même à celle-ci

XC.
Comme notre dessein n’est pas d’avoir égard dans les valeurs de
et
aux termes de l’ordre
mais seulement à ceux qui ont des coefficients finis, nous pourrons négliger dans les équations
et
tous les termes qui se trouveront affectés de
parce que ces termes seront encore de l’ordre
après l’intégration.
Or les équations
et
donnent, en rejetant les termes affectés de

d’où l’on tire.

et, intégrant,

il ne faut point ici de constante, ce qui est évident par la nature de nos formules ; on trouvera de même

donc, substituant ces valeurs de

dans l’équation

de l’Article précédent, on changera les termes

en ceux-ci

et les termes

en ceux-ci

que l’on peut encore changer en ceux-ci

Par ce moyen, l’équation
ne contiendra plus que des termes de la forme de

Je reprends maintenant l’équation

laquelle, étant rapportée au second satellite, devient

je multiplie cette dernière par
je l’intègre, j’ai

je change l’expression

en son équivalente


et il me vient l’équation

![{\displaystyle +\left[\mathrm {M} _{2}^{2}-(\mu _{2}-\mu _{1})^{2}\right]\int \mathrm {x} _{2}\sin(\mu _{2}-\mu _{1})tdt=0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ea0afc0ba54e7aa02be6a29764e6c61f709b549)
d’où je tire

Je trouve, de la même manière,

On fera toutes ces substitutions dans l’équation
moyennant quoi elle n’aura plus que des termes de la forme de

XCI.
Donc si l’on fait, pour abréger,
![{\displaystyle {\begin{aligned}{\overset {\circ }{\Psi }}_{1}=&{\widehat {\Psi }}_{1}(a_{1},a_{2})-4{\frac {\mu _{1}\mu _{2}}{\mathrm {M} _{2}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\\&\quad -\left[2{\widehat {\Psi }}_{1}(a_{1},a_{2})+4{\frac {\mu _{2}(\mu _{2}-\mu _{1})}{\mathrm {M} _{2}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\right]{\frac {\mu _{1}(\mu _{2}-\mu _{1})}{\mathrm {M} _{2}^{2}-(\mu _{2}-\mu _{1})^{2}}},\\{\overset {\circ }{\Psi }}_{1}=&{\widehat {\Psi }}_{1}(a_{1},a_{3})-4{\frac {\mu _{1}\mu _{3}}{\mathrm {M} _{3}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\\&\quad -\left[2{\widehat {\Psi }}_{1}(a_{1},a_{3})+4{\frac {\mu _{3}(\mu _{3}-\mu _{1})}{\mathrm {M} _{3}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\right]{\frac {\mu _{1}(\mu _{3}-\mu _{1})}{\mathrm {M} _{3}^{2}-(\mu _{3}-\mu _{1})^{2}}},\\{\overset {\circ }{\Psi }}_{1}=&{\widehat {\Psi }}_{1}(a_{1},a_{4})-4{\frac {\mu _{1}\mu _{4}}{\mathrm {M} _{4}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{4})\\&\quad -\left[2{\widehat {\Psi }}_{1}(a_{1},a_{4})+4{\frac {\mu _{4}(\mu _{4}-\mu _{1})}{\mathrm {M} _{4}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{4})\right]{\frac {\mu _{1}(\mu _{4}-\mu _{1})}{\mathrm {M} _{4}^{2}-(\mu _{4}-\mu _{1})^{2}}},\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a3ae1e32fbfe542d2b96a27cc69e335920baf0c)
et de plus
![{\displaystyle {\begin{aligned}{\overset {\circ }{\Pi }}_{1}=&{\frac {-2\mu _{2}}{\mathrm {M} _{2}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\\&\quad +\left[2{\widehat {\Psi }}_{1}(a_{1},a_{2})+4{\frac {\mu _{2}(\mu _{2}-\mu _{1})}{\mathrm {M} _{2}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{2})\right]{\frac {\mu _{1}}{\mathrm {M} _{2}^{2}-(\mu _{2}-\mu _{1})^{2}}},\\{\overset {\circ }{\Pi }}_{1}=&{\frac {-2\mu _{3}}{\mathrm {M} _{3}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\\&\quad +\left[2{\widehat {\Psi }}_{1}(a_{1},a_{3})+4{\frac {\mu _{3}(\mu _{3}-\mu _{1})}{\mathrm {M} _{3}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{3})\right]{\frac {\mu _{1}}{\mathrm {M} _{3}^{2}-(\mu _{3}-\mu _{1})^{2}}},\\{\overset {\circ }{\Pi }}_{1}=&{\frac {-2\mu _{4}}{\mathrm {M} _{4}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{4})\\&\quad +\left[2{\widehat {\Psi }}_{1}(a_{1},a_{4})+4{\frac {\mu _{4}(\mu _{4}-\mu _{1})}{\mathrm {M} _{4}^{2}}}{\widehat {\Gamma }}_{1}(a_{1},a_{4})\right]{\frac {\mu _{1}}{\mathrm {M} _{4}^{2}-(\mu _{4}-\mu _{1})^{2}}},\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/860a8238e4632f6f3d2d69fc8bf5cb04a1ac97fc)
on aura, pour le premier satellite, l’équation

et de même pour les trois autres satellites


Pareillement on aura, par rapport aux variables
ces quatre équations (Article LXXXIX)


§ III. — Où l’on donne une nouvelle méthode pour intégrer les équations précédentes.
XCII.
Je fais

d’où je tire

Je substitue ces valeurs dans l’équation
ce qui la change en celle-ci
![{\displaystyle {\begin{aligned}&{\frac {d^{2}\mathrm {x} _{1}}{dt^{2}}}+\mathrm {M} _{1}^{2}\mathrm {x} _{1}\\&-nf_{1}\chi _{2}\left[{\overset {\circ }{\Psi }}_{1}(a_{1},a_{2})-(\mu _{2}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{1},a_{2})\right]\mathrm {P} -nf_{1}\chi _{2}{\overset {\circ }{\Pi }}_{1}(a_{1},a_{2}){\frac {dp}{dt}}\\&-nf_{1}\chi _{3}\left[{\overset {\circ }{\Psi }}_{1}(a_{1},a_{3})-(\mu _{3}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{1},a_{3})\right]\mathrm {Q} -nf_{1}\chi _{3}{\overset {\circ }{\Pi }}_{1}(a_{1},a_{3}){\frac {dq}{dt}}\\&-nf_{1}\chi _{4}\left[{\overset {\circ }{\Psi }}_{1}(a_{1},a_{4})-(\mu _{4}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{1},a_{4})\right]\mathrm {R} -nf_{1}\chi _{4}{\overset {\circ }{\Pi }}_{1}(a_{1},a_{4}){\frac {dr}{dt}}=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cc9f3b17732143ddf9c0308241f243e9eb29d74e)
C’est l’équation qu’il s’agit maintenant d’intégrer, en regardant les quantités
chacune comme une variable particulière. Pour y parvenir, voici comment je m’y prends.
XCIII.
Je reprends les formules

et je trouve de même

De là je tire, par la différentiation, les formules suivantes

Cela posé, je multiplie d’abord l’équation
par
j’ai
![{\displaystyle {\begin{aligned}{\frac {d^{2}\mathrm {x} _{2}}{dt^{2}}}&\sin(\mu _{2}-\mu _{1})t+\mathrm {M} _{2}^{2}\mathrm {x} _{2}\sin(\mu _{2}-\mu _{1})t\\&-{\frac {n}{2}}f_{2}\chi _{1}{\overset {\circ }{\Psi }}_{1}(a_{2},a_{1})\mathrm {x} _{1}\sin 2(\mu _{2}-\mu _{1})t\\&-{\frac {n}{2}}f_{2}\chi _{1}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{1}){\frac {d\mathrm {x} _{1}}{dt}}\left[-1+\cos 2(\mu _{2}-\mu _{1})t\right]\\&-{\frac {n}{2}}f_{2}\chi _{3}{\overset {\circ }{\Psi }}_{1}(a_{2},a_{3})\mathrm {x} _{3}\left[\sin(\mu _{3}-\mu _{1})t-\sin(\mu _{3}-2\mu _{2}+\mu _{1})t\right]\\&-{\frac {n}{2}}f_{2}\chi _{3}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{3}){\frac {d\mathrm {x} _{3}}{dt}}\left[-\cos(\mu _{3}-\mu _{1})t+\cos(\mu _{3}-2\mu _{2}+\mu _{1})t\right]\\&-{\frac {n}{2}}f_{2}\chi _{4}{\overset {\circ }{\Psi }}_{1}(a_{2},a_{4})\mathrm {x} _{4}\left[\sin(\mu _{4}-\mu _{1})t-\sin(\mu _{4}-2\mu _{2}+\mu _{1})t\right]\\&-{\frac {n}{2}}f_{2}\chi _{4}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{4}){\frac {d\mathrm {x} _{4}}{dt}}\left[-\cos(\mu _{4}-\mu _{1})t+\cos(\mu _{4}-2\mu _{2}+\mu _{1})t\right]=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/855ec6edf86aa2163df25e4ed3bc75ac69f42ce6)
Je ne conserve dans cette équation que les termes analogues à ceux de l’équation
c’est-à-dire, les termes qui, en faisant pour
les substitutions de l’Article LXXXVII, en donneraient d’autres où le coefficient de
serait presque égal à
et qui sont les seuls auxquels nous devions avoir égard dans l’intégration de l’équation de l’Article XCII.
J’aurai donc simplement

Je substitue au lieu de

leurs valeurs en
j’ai
(1o)
![{\displaystyle {\begin{aligned}{\frac {d^{2}p}{dt^{2}}}&-2(\mu _{2}-\mu _{1}){\frac {d\mathrm {P} }{dt}}\left[\mathrm {M} _{2}^{2}-(\mu _{2}-\mu _{1})^{2}\right]p+{\frac {n}{2}}f_{2}\chi _{1}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{1}){\frac {d\mathrm {x} _{1}}{dt}}\\&-{\frac {n}{2}}f_{2}\chi _{3}\left[{\overset {\circ }{\Psi }}_{1}(a_{2},a_{3})-(\mu _{3}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{2},a_{3})\right]q+{\frac {n}{2}}f_{2}\chi _{3}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{3}){\frac {d\mathrm {Q} }{dt}}\\&-{\frac {n}{2}}f_{2}\chi _{4}\left[{\overset {\circ }{\Psi }}_{1}(a_{2},a_{4})-(\mu _{4}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{2},a_{4})\right]r+{\frac {n}{2}}f_{2}\chi _{4}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{4}){\frac {d\mathrm {R} }{dt}}=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b36759645b19b266c6b21e5d95217f7ac7dabd04)
Je multiplie en second lieu la même équation
par
j’ai
![{\displaystyle {\begin{aligned}{\frac {d^{2}\mathrm {x} _{2}}{dt^{2}}}&\cos(\mu _{2}-\mu _{1})t+\mathrm {M} _{2}^{2}\mathrm {x} _{2}\cos(\mu _{2}-\mu _{1})t\\&-{\frac {n}{2}}f_{2}\chi _{1}{\overset {\circ }{\Psi }}_{1}(a_{2},a_{1})\mathrm {x} _{1}\left[1+\cos 2(\mu _{2}-\mu _{1})t\right]\\&+{\frac {n}{2}}f_{2}\chi _{1}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{1}){\frac {d\mathrm {x} _{1}}{dt}}\sin 2(\mu _{2}-\mu _{1})t\\&-{\frac {n}{2}}f_{2}\chi _{3}{\overset {\circ }{\Psi }}_{1}(a_{2},a_{3})\ \mathrm {x} _{3}\ \left[\cos(\mu _{3}-\mu _{1})t+\cos(\mu _{3}-2\mu _{2}+\mu _{1})t\right]\\&-{\frac {n}{2}}f_{2}\chi _{3}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{3}){\frac {d\mathrm {x} _{3}}{dt}}\left[\sin(\mu _{3}-\mu _{1})t+\sin(\mu _{3}-2\mu _{2}+\mu _{1})t\right]\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e4e8e196ae0232a495b7df71f0f5aa21e8e05be2)
![{\displaystyle {\begin{aligned}&-{\frac {n}{2}}f_{2}\chi _{4}{\overset {\circ }{\Psi }}_{1}(a_{2},a_{4})\ \mathrm {x} _{4}\ \left[\cos(\mu _{4}-\mu _{1})t+\cos(\mu _{4}-2\mu _{2}+\mu _{1})t\right]\\&-{\frac {n}{2}}f_{2}\chi _{4}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{4}){\frac {d\mathrm {x} _{4}}{dt}}\left[\sin(\mu _{4}-\mu _{1})t+\sin(\mu _{4}-2\mu _{2}+\mu _{1})t\right]=0,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/32daa1f22e529c24009ea5426fad641e5ead5045)
équation que je réduis, par la raison que j’ai dite tantôt, à celle-ci

laquelle me donne, après les substitutions,
(2o).
![{\displaystyle {\begin{aligned}{\frac {d^{2}\mathrm {P} }{dt^{2}}}&+2(\mu _{2}-\mu _{1}){\frac {dp}{dt}}+\left[\mathrm {M} _{2}^{2}-(\mu _{2}-\mu _{1})^{2}\right]\mathrm {P} -{\frac {n}{2}}f_{2}\chi _{1}{\overset {\circ }{\Psi }}_{1}(a_{2},a_{1})\mathrm {x} _{1}\\&-{\frac {n}{2}}f_{2}\chi _{3}\left[{\overset {\circ }{\Psi }}_{1}(a_{2},a_{3})-(\mu _{3}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{2},a_{3})\right]\mathrm {Q} -{\frac {n}{2}}f_{2}\chi _{3}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{3}){\frac {dq}{dt}}\\&-{\frac {n}{2}}f_{2}\chi _{4}\left[{\overset {\circ }{\Psi }}_{1}(a_{2},a_{4})-(\mu _{4}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{2},a_{4})\right]\mathrm {R} -{\frac {n}{2}}f_{2}\chi _{4}{\overset {\circ }{\Pi }}_{1}(a_{2},a_{4}){\frac {dr}{dt}}=0..\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/46b4dc99734f1f793f0c445bd325fa53144721e8)
En troisième lieu, je multiplie l’équation
par
j’aurai, après les réductions et les substitutions,
(3o).
![{\displaystyle {\begin{aligned}{\frac {d^{2}q}{dt^{2}}}&-2(\mu _{3}-\mu _{1}){\frac {d\mathrm {Q} }{dt}}+\left[\mathrm {M} _{3}^{2}-(\mu _{3}-\mu _{1})^{2}\right]q+{\frac {n}{2}}f_{3}\chi _{1}{\overset {\circ }{\Psi }}_{1}(a_{3},a_{1}){\frac {d\mathrm {x} _{1}}{dt}}\\&-{\frac {n}{2}}f_{3}\chi _{2}\left[{\overset {\circ }{\Psi }}_{1}(a_{3},a_{2})-(\mu _{2}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{3},a_{2})\right]p+{\frac {n}{2}}f_{3}\chi _{2}{\overset {\circ }{\Pi }}_{1}(a_{3},a_{2}){\frac {d\mathrm {P} }{dt}}\\&-{\frac {n}{2}}f_{3}\chi _{4}\left[{\overset {\circ }{\Psi }}_{1}(a_{3},a_{4})-(\mu _{4}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{3},a_{4})\right]r+{\frac {n}{2}}f_{3}\chi _{4}{\overset {\circ }{\Pi }}_{1}(a_{3},a_{4}){\frac {d\mathrm {R} }{dt}}=0..\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8be6925a98bd7cb826d6ffcc6cacde222b49fbf7)
En quatrième lieu, je multiplie la même équation par
et je trouve
(4o)
![{\displaystyle {\begin{aligned}{\frac {d^{2}\mathrm {Q} }{dt^{2}}}&+2(\mu _{3}-\mu _{1}){\frac {dq}{dt}}+\left[\mathrm {M} _{3}^{2}-(\mu _{3}-\mu _{1})^{2}\right]\mathrm {Q} -{\frac {n}{2}}f_{3}\chi _{1}{\overset {\circ }{\Psi }}_{1}(a_{3},a_{1})\mathrm {x} _{1}\\&-{\frac {n}{2}}f_{3}\chi _{2}\left[{\overset {\circ }{\Psi }}_{1}(a_{3},a_{2})-(\mu _{2}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{3},a_{2})\right]\mathrm {P} -{\frac {n}{2}}f_{3}\chi _{2}{\overset {\circ }{\Pi }}_{1}(a_{3},a_{2}){\frac {dp}{dt}}\\&-{\frac {n}{2}}f_{3}\chi _{4}\left[{\overset {\circ }{\Psi }}_{1}(a_{3},a_{4})-(\mu _{4}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{3},a_{4})\right]\mathrm {R} -{\frac {n}{2}}f_{3}\chi _{4}{\overset {\circ }{\Pi }}_{1}(a_{3},a_{4}){\frac {dr}{dt}}=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/729667d9b06492d41dc1e6d48aa4faa9c4ff73f4)
En cinquième lieu, je multiplie l’équation
par
j’ai
(5o)
![{\displaystyle {\begin{aligned}{\frac {d^{2}r}{dt^{2}}}&-2(\mu _{4}-\mu _{1}){\frac {d\mathrm {R} }{dt}}+\left[\mathrm {M} _{4}^{2}-(\mu _{4}-\mu _{1})^{2}\right]r+{\frac {n}{2}}f_{4}\chi _{1}{\overset {\circ }{\Psi }}_{1}(a_{4},a_{1}){\frac {d\mathrm {x} _{1}}{dt}}\\&-{\frac {n}{2}}f_{4}\chi _{2}\left[{\overset {\circ }{\Psi }}_{1}(a_{4},a_{2})-(\mu _{2}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{4},a_{2})\right]p+{\frac {n}{2}}f_{4}\chi _{2}{\overset {\circ }{\Pi }}_{1}(a_{4},a_{2}){\frac {d\mathrm {P} }{dt}}\\&-{\frac {n}{2}}f_{4}\chi _{3}\left[{\overset {\circ }{\Psi }}_{1}(a_{4},a_{3})-(\mu _{3}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{4},a_{3})\right]q+{\frac {n}{2}}f_{4}\chi _{3}{\overset {\circ }{\Pi }}_{1}(a_{4},a_{3}){\frac {d\mathrm {Q} }{dt}}=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a3304dbd22fd4f4de031d1d4349004f9b7397f2)
En sixième et dernier lieu, je multiplie la même équation par
et je trouve
(6o)
![{\displaystyle {\begin{aligned}{\frac {d^{2}\mathrm {R} }{dt^{2}}}&+2(\mu _{4}-\mu _{1}){\frac {dr}{dt}}+\left[\mathrm {M} _{4}^{2}+(\mu _{4}-\mu _{1})^{2}\right]\mathrm {R} -{\frac {n}{2}}f_{4}\chi _{1}{\overset {\circ }{\Psi }}_{1}(a_{4},a_{1})\mathrm {x} _{1}\\&-{\frac {n}{2}}f_{4}\chi _{2}\left[{\overset {\circ }{\Psi }}_{1}(a_{4},a_{2})-(\mu _{2}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{4},a_{2})\right]\mathrm {P} -{\frac {n}{2}}f_{4}\chi _{2}{\overset {\circ }{\Pi }}_{1}(a_{4},a_{2}){\frac {dp}{dt}}\\&-{\frac {n}{2}}f_{4}\chi _{3}\left[{\overset {\circ }{\Psi }}_{1}(a_{4},a_{3})-(\mu _{3}-\mu _{1}){\overset {\circ }{\Pi }}_{1}(a_{4},a_{3})\right]\mathrm {Q} -{\frac {n}{2}}f_{4}\chi _{3}{\overset {\circ }{\Pi }}_{1}(a_{4},a_{3}){\frac {dq}{dt}}=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ed0da5f5044229133f0f2bbab51552637870213b)
Voilà, comme on voit, six équations différentielles de la même nature que l’équation de l’Article XCII, et qui, étant combinées avec cette dernière équation, suffiront pour déterminer les sept variables
XCIV.
Pour cet effet, je multiplie l’équation de l’Article XCII par
l’équation (1o) de l’Article précédent par
l’équation (2o) par
l’équation (3o) par
l’équation (4o) par
l’équation (5o) par
l’équation (6o) par
(
et
sont des constantes indéterminées), et après en avoir fait une somme, j’en prends l’intégrale ; j’ai