THÉORIE
DES VARIATIONS PÉRIODIQUES
DES MOUVEMENTS DES PLANÈTES.
(Nouveaux Mémoires de l’Académie royale des Sciences et Belles-Lettres
de Berlin, année 1783.)
Après avoir donné dans la Théorie des variations séculaires la partie la plus importante et la plus difficile du Calcul des perturbations des Planètes, j’ai cru devoir compléter ce Calcul, en donnant aussi l’autre partie qui concerne les variations périodiques.
Plusieurs Géomètres se sont déjà occupés de ce dernier objet ; mais leurs travaux se trouvent épars dans divers Ouvrages, et les résultats de ces travaux, dépendant de méthodes et de données différentes, présentent une variété embarrassante pour les Astronomes qui voudraient en faire usage. D’ailleurs ils n’ont point tenu compte de l’effet des variations séculaires, et, quoique cet effet ne puisse être que très-petit, la rigueur ne permet pas de le négliger, sans avoir auparavant démontré qu’il n’en saurait résulter que des altérations insensibles dans le mouvement des Planètes.
L’Analyse par laquelle nous avons déterminé les variations séculaires donne aussi directement les variations périodiques ; car nous avons commencé par réduire tout l’effet des perturbations à la variation des éléments des orbites ; nous avons ensuite séparé ce qui, dans ces variations, n’était que périodique de ce qui était indépendant des lieux des Planètes et cette dernière Partie nous a donné les variations séculaires qui affectent immédiatement et continuellement les dimensions et la position des orbites elliptiques. Il ne s’agit donc que d’avoir égard à la partie périodique des variations de ces mêmes éléments, et de tenir compte des termes négligés dans le premier calcul. On aura par là les variations totales des éléments des Planètes causées par leur attraction mutuelle ; et, en employant les éléments corrigés par ces variations, on pourra calculer pour chaque instant les lieux des Planètes par les règles ordinaires.
À la vérité, comme les variations périodiques demeurent toujours très-petites et n’ont, pour ainsi dire, qu’un effet passager et alternatif, il est peu important pour l’Astronomie de connaître en particulier les altérations qui en résultent dans chacun des éléments des orbites ; et il suffit d’avoir l’effet total de ces variations sur les lieux des Planètes. Aussi les Astronomes sont-ils dans l’usage de regarder les orbites comme constantes relativement aux variations périodiques, et de ne traiter ces variations que comme des corrections à faire aux lieux elliptiques calculés par les règles ordinaires. Mais après avoir déterminé par nos formules générales la partie périodique des variations des éléments, rien ne sera plus facile que d’en conclure les inégalités du rayon vecteur, de la longitude et de la latitude. Il est vrai qu’on peut calculer ces inégalités d’une manière plus simple, en les déduisant immédiatement des équations différentielles de l’orbite ; mais cette méthode a, d’un autre côté, l’inconvénient de ne donner les variations séculaires qu’indirectement et par des réductions particulières ; et ne doit-il pas être plus satisfaisant de trouver toutes les inégalités du mouvement des Planètes par une même analyse et par des procédés directs et uniformes ?
1. Nous nous proposons donc ici d’appliquer à la détermination des inégalités périodiques des Planètes les formules générales trouvées dans la première Partie de la Théorie des variations séculaires ; mais auparavant il est bon de faire à ces formules un changement qui servira à les rendre plus simples et plus exactes à quelques égards, sans influer d’ailleurs en rien sur les résultats trouvés pour les variations séculaires. Ce changement est relatif à la manière dont nous avons déterminé l’anomalie vraie par l’anomalie moyenne dans une ellipse dont les éléments sont variables ; voici en quoi il consiste.
Après être parvenus dans le no 30 de la Partie citée à l’équation différentielle
![{\displaystyle dq+\beta \sin qdq+\gamma \cos qdq+\delta \sin 2qdq+\varepsilon \cos 2qdq+\ldots =dp,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e5912b92142d8ed71fa55f6eab85e17c5358b42)
dans laquelle
est l’angle du mouvement moyen,
celui du mouvement vrai, et
des coefficients dépendant des éléments de l’ellipse, nous avons intégré le premier membre par approximation, et nous avons obtenu la formule
![{\displaystyle q-(\beta )\cos q+(\gamma )\sin q-{\frac {1}{2}}(\delta )\cos 2q+{\frac {1}{2}}(\varepsilon )\sin 2q-\ldots =p,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eb917700960e1ff3eb3e55b8e898219d3c925b8b)
d’où nous avons ensuite déduit la valeur de
en ![{\displaystyle p.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/88532f4eab1d4cef71ef96c0f8c98cac36fd9257)
J’ai reconnu depuis qu’au lieu de cette approximation il vaut mieux laisser le premier membre de l’intégrale sous la forme qu’il aurait si les éléments étaient constants, et appliquer à la valeur de
qui constitue le second membre, les corrections résultantes de la variabilité des mêmes éléments.
Ainsi, en faisant
![{\displaystyle d\Sigma =-\cos qd\beta +\sin qd\gamma -{\frac {1}{2}}\cos 2q\,d\delta +{\frac {1}{2}}\sin 2q\,d\varepsilon -\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/82bc69ab173673556957ccc9081b7f5aceffec26)
et ajoutant cette équation à l’équation différentielle ci-dessus, on aura, après l’intégration, cette formule
![{\displaystyle q-\beta \cos q+\gamma \sin q-{\frac {1}{2}}\delta \cos 2q+{\frac {1}{2}}\varepsilon \sin 2q-\ldots =p+\Sigma ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f0fccb732d15cb34bfd9521264a4d751ab965956)
qu’on voit être de la même forme que la précédente, les coefficients ![{\displaystyle (\beta ),(\gamma ),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a11a47b8d09126d6f5ae792ad10ae7914a9d73eb)
étant changés en
et la variable
en ![{\displaystyle p+\Sigma .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/200e92e17371918f9f8216c7a7b4b3d1a329ebcf)
On fera donc ces mêmes changements dans les autres formules dépendantes de celle-ci ; et par conséquent les quantités
et
des nos 38 et suivants se réduiront d’abord à
et
réduction que nous avons faite aussi dans le no 49, mais d’après la considération de la petitesse des forces perturbatrices dont nous avons négligé les carrés. Par cette considération, on pourra donc aussi négliger dans les formules des variations des éléments la quantité, ou du moins la partie périodique de cette quantité, l’autre partie proportionnelle à
se confondant avec le mouvement moyen ; de sorte que les équations différentielles des nos 40 et 50 demeureront les mêmes ; et toute la Théorie des variations séculaires des éléments des orbites planétaires subsistera en son entier.
À l’égard de la quantité
il est clair qu’elle représentera la variation du mouvement moyen due aux forces perturbatrices, variation indépendante du grand axe de l’ellipse, et qu’on pourra rapporter à l’époque du moyen mouvement, laquelle, dans les orbites invariables, contient la sixième constante arbitraire des intégrales, et par conséquent le sixième élément du mouvement elliptique.
2. Pour avoir donc la formule de cette variation, il faut d’abord connaître les valeurs des coefficients,
en
or, par le no 30 de l’Ouvrage cité, on voit que la série
![{\displaystyle \alpha \left(1+\beta \sin q+\gamma \cos q+\delta \sin 2q+\varepsilon \cos 2q+\zeta \sin 3q+\eta \cos 3q+\ldots \right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/63a5af781d3695db32fe677bb6fa0559d1d3dd64)
résulte du développement de la fraction
![{\displaystyle \left[g{\sqrt {\mathrm {R} ^{2}+(\mathrm {Q} \sin q-\mathrm {P} \cos q)^{2}}}+\mathrm {C} \sin q-\mathrm {B} \cos q\right]^{-2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b9cbf8a3c87470d72b31f6cbdff33c84f015a71)
dans laquelle (nos 8, 29)
![{\displaystyle \mathrm {B=RN-PL,\quad C=-RM-QL} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/fdbf45140c0d1cabaca223ba7530d82a4e7267bf)
et
![{\displaystyle \mathrm {L={\frac {MQ-NP}{R}}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/40b4bc713aee3275b4e4eb892efca72fea7b1b2c)
l’on trouve par les méthodes ordinaires, en poussant l’approximation
jusqu’aux troisièmes dimensions de
![{\displaystyle \mathrm {M,N,P,Q} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/cfe8638669527d4b6d44439109526dd1da1e6692)
et faisant
![{\displaystyle g=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ba0027652d8e5f694d4aa1c71ce16c9380ce1186)
conformément au
no 37 du même Ouvrage,
![{\displaystyle {\begin{aligned}\beta =&2\mathrm {M-{\frac {NPQ}{2R^{2}}}+{\frac {M\left(3Q^{2}+P^{2}\right)}{4R^{2}}}} ,\\\gamma =&2\mathrm {N-{\frac {MQP}{2R^{2}}}+{\frac {N\left(3P^{2}+Q^{2}\right)}{4R^{2}}}} ,\\\delta =&3\mathrm {MN+{\frac {PQ}{2R^{2}}}} ,\\\varepsilon =&\mathrm {{\frac {3\left(N^{2}-M^{2}\right)}{2}}+{\frac {Q^{2}-P^{2}}{2R^{2}}}} ,\\\zeta =&\mathrm {3MN^{2}-M^{3}+{\frac {3M\left(Q^{2}-P^{2}\right)}{4R^{2}}}+{\frac {3NPQ}{2R^{2}}}} ,\\\eta =&\mathrm {N^{3}-3NM^{2}+{\frac {3N\left(Q^{2}-P^{2}\right)}{4R^{2}}}-{\frac {3MPQ}{2R^{2}}}} .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/049af88962ac3a3dd36daa7df0eb0a8cf0f5bdbf)
Après avoir différentié ces valeurs, on y substituera pour ![{\displaystyle d\mathrm {M} ,d\mathrm {N} ,d\mathrm {P} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/63a5783bbce0c462acd76686a563b207059c96b4)
celles que nous avons données dans les nos 17 et 41 de l’Ouvrage cité ; on aura ainsi les expressions de
qu’il faudra substituer dans celle de
du numéro précédent.
3. Si dans cette expression de
on veut se contenter d’avoir égard aux premières dimensions des quantités
ainsi que nous l’avons fait dans les formules des variations des éléments de l’orbite, il suffira alors d’avoir égard aux secondes dimensions de ces quantités dans les valeurs de
ce qui les réduira à celles-ci
![{\displaystyle \beta =2\mathrm {M} ,\quad \gamma =2\mathrm {N} ,\quad \delta =\mathrm {3MN+{\frac {PQ}{R^{2}}}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c83c07d357d55af143266680521b0109b3faa8cd)
![{\displaystyle \varepsilon =\mathrm {{\frac {3\left(N^{2}-M^{2}\right)}{2}}+{\frac {Q^{2}-P^{2}}{2R^{2}}}} ,\quad \zeta =0,\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/45cf411ea52f38657648c857037d0f1db3b1743b)
et, comme les valeurs de
et
contiennent déjà les premières dimensions de
et
(ainsi qu’on le voit par les formules du no 39), il est clair qu’on pourra négliger les différences des termes
et
On aura
ainsi, après les substitutions, et en n’employant pour
![{\displaystyle d\mathrm {N} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/4a63a2fa5d72f6383a41c7ba41c5423f10f67ed8)
et
![{\displaystyle d\mathrm {M} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/7d9aa26d4c9da82341023553a2ce27fc9f9d483a)
que les valeurs du
no 42 du même Ouvrage,
![{\displaystyle {\begin{aligned}d\Sigma =&2r^{2}{\frac {d\Omega }{dr}}dq+2{\frac {d\Omega }{dq}}dr+3r{\frac {d\Omega }{dq}}(\mathrm {N} \sin q-\mathrm {M} \cos q)dq\\&+{\frac {3r^{2}}{2}}{\frac {d\Omega }{dr}}(\mathrm {M} \sin q+\mathrm {N} \cos q)+{\frac {3}{2}}{\frac {d\Omega }{dq}}(\mathrm {M} \sin q+\mathrm {N} \cos q)dr\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3fe43b08d0ac94eb917969b3152037ef01e747e6)
et il n’y aura plus qu’à faire les autres substitutions de ce dernier numéro, en se souvenant d’y changer
en
et
en ![{\displaystyle p+\Sigma .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/200e92e17371918f9f8216c7a7b4b3d1a329ebcf)
4. Dénotons par
la valeur moyenne de l’angle
c’est-à-dire la partie indépendante des sinus et cosinus, et proportionnelle au temps ; et par
le demi-grand axe de l’ellipse invariable dans laquelle le mouvement moyen serait
en sorte que l’on ait
il est clair que, dans les termes dus aux forces perturbatrices, on pourra mettre partout
au lieu de
et
au lieu de
demi-grand axe de l’ellipse variable (37), du moins en tant qu’on néglige les carrés de ces forces, comme on l’a toujours pratiqué dans la Théorie des perturbations des corps célestes.
Ainsi les substitutions dont il s’agit se réduiront à faire
![{\displaystyle {\begin{alignedat}{2}q=&{\overline {p}}+2\mathrm {M} \cos {\overline {p}}-2\mathrm {N} \sin {\overline {p}},\qquad &dq=&\left(1-2\mathrm {M} \sin {\overline {p}}-2\mathrm {N} \cos {\overline {p}}\right)d{\overline {p}},\\r=&{\overline {r}}\left(1+\mathrm {M} \sin {\overline {p}}+\mathrm {N} \cos {\overline {p}}\right),&dr=&{\overline {r}}\left(\mathrm {M} \cos {\overline {p}}-\mathrm {N} \sin {\overline {p}}\right)d{\overline {p}},\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3134d36b0d9a1f7293a72beb7b508a7f857d59b2)
![{\displaystyle {\begin{aligned}{\frac {d\Omega }{dr}}=&{\frac {d\Omega }{d{\overline {r}}}}+{\overline {r}}{\frac {d^{2}\Omega }{d{\overline {r}}^{2}}}\left(\mathrm {M} \sin {\overline {p}}+\mathrm {N} \cos {\overline {p}}\right)\\&+{\overline {r}}'{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {r}}'}}\left(\mathrm {M} '\sin {\overline {p}}'+\mathrm {N} '\cos {\overline {p}}'\right)+\ldots \\&+2{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}}}\left(\mathrm {M} \cos {\overline {p}}-\mathrm {N} \sin {\overline {p}}\right)\\&+2{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}'}}\left(\mathrm {M} '\cos {\overline {p}}'-\mathrm {N} '\sin {\overline {p}}'\right)+\ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fef05c510f5ff77316b31ef5a6f987db71e32d7e)
![{\displaystyle {\begin{aligned}{\frac {d\Omega }{dq}}=&{\frac {d\Omega }{d{\overline {p}}}}+{\overline {r}}{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}}}\left(\mathrm {M} \sin {\overline {p}}+\mathrm {N} \cos {\overline {p}}\right)\\&+{\overline {r}}'{\frac {d^{2}\Omega }{d{\overline {r}}'d{\overline {p}}}}\left(\mathrm {M} '\sin {\overline {p}}'+\mathrm {N} '\cos {\overline {p}}'\right)+\ldots \\&+2{\frac {d^{2}\Omega }{d{\overline {p}}^{2}}}\left(\mathrm {M} \cos {\overline {p}}-\mathrm {N} \sin {\overline {p}}\right)\\&+2{\frac {d^{2}\Omega }{d{\overline {p}}d{\overline {p}}'}}\left(\mathrm {M} '\cos {\overline {p}}'-\mathrm {N} '\sin {\overline {p}}'\right)+\ldots ,\\\Omega =&\mathrm {T} '\ \left[{\frac {{\overline {r}}\cos \left({\overline {p}}-{\overline {p}}'\ \right)}{{\overline {r}}'^{2}\ }}-{\frac {1}{\sqrt {{\overline {r}}^{2}-2{\overline {r}}{\overline {r}}'\ \cos \left({\overline {p}}-{\overline {p}}'\ \right)+{\overline {r}}'^{2}}}}\right]\\+&\mathrm {T} ''\left[{\frac {{\overline {r}}\cos \left({\overline {p}}-{\overline {p}}''\right)}{{\overline {r}}''^{2}}}-{\frac {1}{\sqrt {{\overline {r}}^{2}-2{\overline {r}}{\overline {r}}''\cos \left({\overline {p}}-{\overline {p}}''\right)+{\overline {r}}''^{2}}}}\right]\\+&\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/c3d0d0ff6e07aa250755d3ac993b331accec34d6)
Négligeant donc dans ces substitutions les secondes dimensions des quantités
et changeant les lettres
en
ainsi que nous l’avons fait dans le no 50 de l’Ouvrage cité, on aura pour la variation de
cette formule
![{\displaystyle {\begin{aligned}d\Sigma =&2{\overline {r}}^{2}{\frac {d\Omega }{d{\overline {r}}}}d{\overline {p}}+\left(2{\overline {r}}^{3}{\frac {d^{2}\Omega }{d{\overline {r}}^{2}}}+{\frac {3{\overline {r}}^{2}}{2}}{\frac {d\Omega }{d{\overline {r}}}}\right)\left(x\sin {\overline {p}}+y\cos {\overline {p}}\right)d{\overline {p}}\\&+2{\overline {r}}^{2}{\overline {r}}'{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {r}}'}}\left(x'\sin {\overline {p}}'+y'\cos {\overline {p}}'\right)d{\overline {p}}\\&+2{\overline {r}}^{2}{\overline {r}}''{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {r}}''}}\left(x''\sin {\overline {p}}''+y''\cos {\overline {p}}''\right)d{\overline {p}}=\ldots \\&+\left(4{\overline {r}}^{2}{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}}}-{\overline {r}}{\frac {d\Omega }{d{\overline {p}}}}\right)\left(x\cos {\overline {p}}+y\cos {\overline {p}}\right)d{\overline {p}}\\&+4{\overline {r}}^{2}{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}'}}\left(x'\cos {\overline {p}}'+y'\cos {\overline {p}}'\right)d{\overline {p}}\\&+4{\overline {r}}^{2}{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}''}}\left(x''\cos {\overline {p}}''+y''\cos {\overline {p}}''\right)d{\overline {p}}+\ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d092cb3090c3b34257070dfa8c5ce3a1f488763)
dans laquelle il ne s’agira plus que de substituer la valeur précédente de ![{\displaystyle \Omega .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5dc55a0f828ae7d92ad5426ce7348d84ee62a70)
5. Mais, pour avoir une formule intégrable, il est nécessaire de réduire en série les radicaux qui entrent dans l’expression de
on supposera donc
![{\displaystyle \left[{\overline {r}}^{2}-2{\overline {r}}{\overline {r}}'\cos({\overline {p}}-{\overline {p}}')+{\overline {r}}'^{2}\right]^{-{\frac {1}{2}}}=\left[{\overline {r}},{\overline {r}}'\right]+\left[{\overline {r}},{\overline {r}}'\right]_{1}\cos \left({\overline {p}}-{\overline {p}}'\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4ab9d05c2bcd3242efd927f43e09be53d4b60782)
![{\displaystyle +\left[{\overline {r}},{\overline {r}}'\right]_{2}\cos 2\left({\overline {p}}-{\overline {p}}'\right)+\left[{\overline {r}},{\overline {r}}'\right]_{3}\cos 3\left({\overline {p}}-{\overline {p}}'\right)+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a2b85a6d81d32da6e4b554471c61f72ccc4e5c09)
et ainsi des autres radicaux semblables ; alors l’expression de
se trouvera composée de termes tous intégrables, dont les uns, simplement proportionnels à
seront
![{\displaystyle -2\mathrm {T} '{\overline {r}}^{2}{\frac {d[{\overline {r}},{\overline {r}}']}{d{\overline {r}}'}}d{\overline {p}}-2\mathrm {T} ''{\overline {r}}^{2}{\frac {d[{\overline {r}},{\overline {r}}'']}{d{\overline {r}}''}}d{\overline {p}}-\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/750db2cb8e5e6dc041fe50b6ef07f5932e43538b)
et dont les autres contiendront tous des sinus ou cosinus des angles ![{\displaystyle {\overline {p}},{\overline {p}}',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/da5d8c611d095129fa6f045ce7e7377dc217cac9)
et donneront par conséquent la partié périodique de la valeur de ![{\displaystyle \Sigma .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/203943be81c35e85bc0f08a1f8822bf725dc43a4)
Quant aux premiers termes, nous verrons ci-après comment, étant joints à un terme semblable provenant de la valeur de
il en résulte le terme unique
que nous avons supposé représenter la partie uniforme du mouvement moyen
dans l’ellipse variable.
Ainsi, tant qu’on n’a égard qu’aux premières dimensions des excentricités et des inclinaisons, auxquelles les quantités
ou ![{\displaystyle x,y,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4ad6024a5a4f15ab696de3d6e4866da50c8a315d)
sont proportionnelles, on ne trouve point de variation séculaire dans le mouvement moyen des Planètes ; mais, si l’on poussait l’approximation jusqu’aux, secondes dimensions de ces quantités dans l’expression de
on trouverait alors des termes proportionnels à
et indépendants des sinus et cosinus de
lesquels auraient pour coefficients les quantités
![{\displaystyle x^{2}+y^{2},\ \ x'^{2}+y'^{2},\ \ xx'+yy',\ldots ,\ \ s^{2}+u^{2},\ \ s'^{2}+u'^{2},\ \ ss'+uu',\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a705345178757ec3093db5d34dd0dc8e176d3917)
ces termes donneraient donc, par la substitution des valeurs de ![{\displaystyle x,y,x',y',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5d74fabc023d020f5e287a1d6c216c5d2a9fe2fb)
![{\displaystyle \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a3dbfc5796975effdfc4a5e30c7b0ce9e80e0d5f)
relatives à chaque Planète, des variations séculaires dans leur mouvement moyen, variations qui n’affecteraient que le
mouvement circulatoire dans l’ellipse variable, mais nullement le demigrand axe ou la distance moyenne. Comme cet objet est de la plus grande importance dans la Théorie des Planètes, je me propose de le discuter à fond dans une autre occasion
[1].
6. Si dans l’expression de
du no 4 on néglige les termes où se trouvent les quantités
dues aux excentricités des orbites, elle se réduit à celle-ci
![{\displaystyle {\begin{aligned}d\Sigma =&2{\overline {r}}^{2}{\frac {d\Omega }{d{\overline {r}}}}d{\overline {p}}\\=&-2{\overline {r}}^{2}\mathrm {T} '\ \left[{\frac {d[{\overline {r}},{\overline {r}}'\ ]}{d{\overline {r}}}}+\left({\frac {d[{\overline {r}},{\overline {r}}'\ ]_{1}}{d{\overline {r}}}}-{\frac {1}{{\overline {r}}'^{2}\ }}\right)\cos({\overline {p}}-{\overline {p}}'\ )\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad \left.+{\frac {d[{\overline {r}},{\overline {r}}'\ ]_{2}}{d{\overline {r}}}}\cos 2({\overline {p}}-{\overline {p}}'\ )+\ldots \right]d{\overline {p}}\\=&-2{\overline {r}}^{2}\mathrm {T} ''\ \left[{\frac {d[{\overline {r}},{\overline {r}}'']}{d{\overline {r}}}}+\left({\frac {d[{\overline {r}},{\overline {r}}'']_{1}}{d{\overline {r}}}}-{\frac {1}{{\overline {r}}''^{2}}}\right)\cos({\overline {p}}-{\overline {p}}'')\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad \left.+{\frac {d[{\overline {r}},{\overline {r}}'']_{2}}{d{\overline {r}}}}\cos 2({\overline {p}}-{\overline {p}}'')+\ldots \right]d{\overline {p}}\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b964311bac5e5dae3c93a48a273ead8cba1e8248)
Soient
![{\displaystyle d{\overline {p}}-d{\overline {p}}'=n'd{\overline {p}},\quad d{\overline {p}}-d{\overline {p}}''=n''d{\overline {p}},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a041a3e25b232656a2352da14de4b9f22b11bac)
les coefficients
étant constants et donnés par les rapports connus des mouvements moyens des Planètes
puisque ![{\displaystyle d{\overline {p}}={\frac {dt}{{\overline {r}}^{\frac {3}{2}}}},d{\overline {p}}'={\frac {dt}{{\overline {r}}'^{\frac {3}{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/68808e70e98d46b192eee4e8f398e9af66b593c5)
![{\displaystyle \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a3dbfc5796975effdfc4a5e30c7b0ce9e80e0d5f)
sera
on aura, en substituant pour
ou
ou
une équation intégrable, et dont l’intégrale sera
![{\displaystyle {\begin{aligned}\Sigma =&-2{\overline {r}}^{2}\left(\mathrm {T} '{\frac {d[{\overline {r}},{\overline {r}}']_{1}}{d{\overline {r}}}}+\mathrm {T} ''{\frac {d[{\overline {r}},{\overline {r}}'']_{1}}{d{\overline {r}}}}+\ldots \right){\overline {p}}\\&-{\frac {2{\overline {r}}^{2}\mathrm {T} '\ }{n'}}\left[\left({\frac {d[{\overline {r}},{\overline {r}}'\ ]_{1}}{d{\overline {r}}}}-{\frac {1}{{\overline {r}}'^{2}}}\right)\sin({\overline {p}}-{\overline {p}}'\ )\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad \left.+{\frac {1}{2}}{\frac {d[{\overline {r}},{\overline {r}}'\ ]_{2}}{d{\overline {r}}}}\sin 2({\overline {p}}-{\overline {p}}'\ )+\ldots \right]\\&-{\frac {2{\overline {r}}^{2}\mathrm {T} ''\ }{n''}}\left[\left({\frac {d[{\overline {r}},{\overline {r}}''\ ]_{1}}{d{\overline {r}}}}-{\frac {1}{{\overline {r}}''^{2}}}\right)\sin({\overline {p}}-{\overline {p}}''\ )\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad \left.+{\frac {1}{2}}{\frac {d[{\overline {r}},{\overline {r}}''\ ]_{2}}{d{\overline {r}}}}\sin 2({\overline {p}}-{\overline {p}}''\ )+\ldots \right]\\&-\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/682dece2e794c8e9417056a32debdd3a01347afe)
7. Si l’on voulait aussi tenir compte des excentricités dans la valeur de
il n’y aurait qu’à substituer dans l’expression de
les valeurs de
déterminées dans la Théorie des variations séculaires ; ce qui ne donnerait que des termes en sinus et cosinus d’angles proportionnels à
ou
et par conséquent tous intégrables.
Mais pour simplifier ce calcul on remarquera :
1o Que l’on peut regarder les quantités
comme constantes, puisque leurs différentielles étant de l’ordre des forces perturbatrices doivent être négligées, par la raison que ces quantités se trouvent ici déjà multipliées par la quantité
qui est aussi du même ordre.
2o Qu’à cause de la petitesse des termes dont il s’agit, il suffira d’lavoir égard à ceux qui pourront augmenter beaucoup par l’intégration ; et il est clair que ce sont les termes qui contiendront des sinus ou cosinus d’angles dont la variation sera très-petite vis-à-vis de celle de l’angle
car si, par exemple,
désigne un de ces angles, et que
étant un coefficients fort petit, la première intégration introduira au dénominateur le coefficient
et la seconde y introduira le carré
Comme l’angle
ne peut être composé que de multiples des angles
joints ensemble par les signes
ou
ce n’est que par les rapports connus des mouvements moyens des Planètes qu’on pourra juger dans chaque cas de la valeur du coefficient
il faudra donc un examen particulier pour chaque Planète dont on voudra calculer les perturbations.
3o Que, s’il arrivait que le coefficient
fût très-petit du même ordre que les coefficients de
dans les sinus ou cosinus des valeurs de
il faudrait alors, dans les termes qui contiendraient les sinus ou cosinus de
substituer immédiatement ces valeurs, et intégrer ensuite, après avoir réduit les produits des sinus et cosinus en sinus et cosinus simples. Mais ce cas ne paraît pas pouvoir exister dans notre Système planétaire.
8. Considérons maintenant les variations du demi-grand axe de l’ellipse. Nous avons trouvé, dans le no 35 de la Théorie citée, que ce grand axe est constant, relativement aux variations séculaires ; mais il ne l’est plus lorsqu’on a égard aux variations périodiques. Car la quantité
qui est en raison inverse du demi-grand axe, est telle que
![{\displaystyle d\Delta =2(d\Omega ),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6f00032d876b393928d367d16e501e56fc532d85)
en désignant par
la différentielle partielle de la fonction
relativement aux seules-variables de la Planète troublée ; par conséquent cette quantité sera essentiellement variable, mais ses variations ne seront que périodiques, puisque l’expression de
ne peut contenir que des termes affectés de sinus et cosinus.
Prenons pour un moment
pour dénoter le demi-grand axe de l’ellipse variable dans laquelle
est l’angle du mouvement moyen ; on aura par les nos 34 et 37 de la Théorie citée
![{\displaystyle r={\frac {1}{\Delta }},\quad {\frac {dp}{dt}}=\Delta ^{\frac {3}{2}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ab6a87b5904ee58d3be2dea81816078e969e117e)
Donc : 1o
![{\displaystyle d{\frac {1}{r}}=2(d\Omega ),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a025911b1560ae8bb83ec06b3cc2e0e10a05946)
et, intégrant,
![{\displaystyle {\frac {1}{r}}={\frac {1}{\overline {r}}}+2\int (d\Omega ).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9836146c0ede9e450851528eac6b869ef3580062)
La quantité
est une constante arbitraire qui serait égale à la valeur moyenne de
si l’intégrale
ne contenait aucun terme constant. Mais nous supposerons que cette intégrale contienne encore une petite constante arbitraire, que nous déterminerons en sorte que le mouvement moyen dans l’ellipse qui aurait
pour demi-grand axe soit le même que la partie uniforme
de l’angle
conformément à ce que nous avons supposé dans le no 4.
Si donc on fait
la quantité
exprimant ainsi les variations périodiques du demi-grand axe, et qu’on néglige les secondes puissances de cette quantité, on aura
![{\displaystyle \rho =-2{\overline {r}}^{2}\int (d\Omega ).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2dc2a78a21d454adda078e2941a34f16d4e160ee)
2o On aura par conséquent
![{\displaystyle dp={\frac {dt}{r^{\frac {3}{2}}}}={\frac {dt}{{\overline {r}}^{\frac {3}{2}}}}-{\frac {3\rho dt}{2r^{\frac {5}{2}}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/29f3fad6fb9cb028172bd662d72b86b89dbe4c61)
Mais
![{\displaystyle {\frac {dt}{{\overline {r}}^{\frac {3}{2}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/406c1f8d8ef2361bc0eb3286feec47dd84b63c9c)
est la valeur de
![{\displaystyle dp}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fbfc7a78f0e04743a4bcdc18bf6dd926e4f215a7)
qui aurait lieu dans l’ellipse dont le demi-grand axe serait
![{\displaystyle {\overline {r}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cf414b822992e16b6735a3c9b5173868e6db9226)
et que nous avons déjà désignée par
![{\displaystyle d\,{\overline {p}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e25df5f1563a913f54e0b82a71f37fdf74b56e5c)
; donc, intégrant,
![{\displaystyle p={\overline {p}}-{\frac {3}{2{\overline {r}}^{\frac {5}{2}}}}\int \rho dt={\overline {p}}+3{\overline {r}}\int d\,{\overline {p}}\int (d\Omega ).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78083b08b003b602506fc149eef3e9d98d80428f)
Donc, si l’on fait
la quantité
exprimant la partie périodique de l’angle
on aura
![{\displaystyle \varpi =\Sigma +3{\overline {r}}\int d\,{\overline {p}}\int (d\Omega )\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a61c65598d846d2d19f89503034ab91b3697c00)
et il faudra que dans cette expression il n’entre aucun terme proportionnel à
.
9. Il ne s’agit maintenant que d’évaluer la quantité
Or on a vu dans le no 41 de la Théorie citée que cette quantité est réductible à la forme
![{\displaystyle {\frac {d\Omega }{dr}}dr+{\frac {d\Omega }{dq}}dq+{\frac {d\Omega }{dz}}dz,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3f774033adc622a4ce968e814750b42fb2980e36)
étant ici le rayon vecteur de l’orbite ; donc, négligeant le terme
parce qu’il contiendrait les secondes puissances des inclinaisons des orbites, et faisant dans les autres les mêmes substitutions et réductions que dans le no 4 ci-dessus, on aura
![{\displaystyle {\begin{aligned}(d\Omega )=&{\frac {d\Omega }{d\,{\overline {p}}}}d\,{\overline {p}}+\left({\overline {r}}{\frac {d^{2}\Omega }{d\,{\overline {r}}d\,{\overline {p}}}}-2{\frac {d\Omega }{d\,{\overline {p}}}}\right)\left(x\sin {\overline {p}}+y\cos {\overline {p}}\right)d\,{\overline {p}}\\&+{\overline {r}}'{\frac {d^{2}\Omega }{d\,{\overline {r}}'d\,{\overline {p}}}}\left(x'\sin {\overline {p}}'+y'\cos {\overline {p}}'\right)d\,{\overline {p}}\\&+{\overline {r}}''{\frac {d^{2}\Omega }{d\,{\overline {r}}''d\,{\overline {p}}}}\left(x''\sin {\overline {p}}''+y''\cos {\overline {p}}''\right)d\,{\overline {p}}+\ldots \\&+\left(2{\frac {d^{2}\Omega }{d\,{\overline {p}}^{2}}}+{\overline {r}}'{\frac {d\Omega }{d{\overline {r}}}}\right)\left(x\cos {\overline {p}}-y\sin {\overline {p}}\right)d\,{\overline {p}}\\&+2{\frac {d^{2}\Omega }{d\,{\overline {p}}d\,{\overline {p}}'}}\left(x'\cos {\overline {p}}'-y'\sin {\overline {p}}'\right)d\,{\overline {p}}\\&+2{\frac {d^{2}\Omega }{d\,{\overline {p}}d\,{\overline {p}}''}}\left(x''\cos {\overline {p}}''-y''\sin {\overline {p}}''\right)d\,{\overline {p}}+\ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e81d3e48d27ed2119047a287f69baef77f2d1c31)
Et l’on intégrera cette formule d’une manière semblable à celle que nous avons indiquée plus haut pour l’intégration de la différentielle
10. Si l’on néglige les excentricités et par conséquent les quantités
qui en dépendent, on aura simplement
![{\displaystyle (d\Omega )={\frac {d\Omega }{d\,{\overline {p}}}}d\,{\overline {p}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c7cb2489f39d6074fbdf6cc6df9afa37f2780fab)
Cette différentielle est intégrable rigoureusement ; car en substituant la valeur de
du no 4 ci-dessus, et faisant comme dans le no 6
![{\displaystyle d\,{\overline {p}}={\frac {d\,{\overline {p}}-d\,{\overline {p}}'}{n'}}={\frac {d\,{\overline {p}}-d\,{\overline {p}}''}{n''}}=\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc0bb1fe82249ede7cc6c514d4fb337283ed497d)
on aura
![{\displaystyle {\begin{aligned}\int (d\Omega )&={\frac {\mathrm {T} '}{n'}}\left[{\frac {{\overline {r}}\cos \left({\overline {p}}-{\overline {p}}'\ \right)}{{\overline {r}}'^{2}}}-{\frac {1}{\sqrt {{\overline {r}}^{2}-2{\overline {r}}{\overline {r}}'\ \cos \left({\overline {p}}-{\overline {p}}'\ \right)+{\overline {r}}'^{2}}}}\right]\\&+{\frac {\mathrm {T} ''}{n''}}\left[{\frac {{\overline {r}}\cos \left({\overline {p}}-{\overline {p}}''\right)}{{\overline {r}}''^{2}}}-{\frac {1}{\sqrt {{\overline {r}}^{2}-2{\overline {r}}{\overline {r}}''\cos \left({\overline {p}}-{\overline {p}}''\right)+{\overline {r}}''^{2}}}}\right]\\&+\ldots +\chi .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d95d23b227c19d22efe83226d7e4645c9a552c07)
étant une constante arbitraire.
Cette expression multipliée par
et prise négativement donnera la valeur de
c’est-à-dire les variations du demi-grand axe en tant qu’on néglige les excentricités des orbites. Mais pour en déduire ensuite celles du mouvement moyen qui en dépendent, il faudra une nouvelle intégration qui, n’étant pas possible en général, oblige d’avoir recours aux séries. On pourrait, à la vérité, construire l’intégrale par les méthodes connues pour la quadrature arithmétique des courbes, et l’on pourrait en faire de même pour la valeur de
dans ce même cas ; mais on n’aurait pas de cette manière des formules générales qui fassent connaître la marche des variations, et l’on n’aurait pas même une plus grande exactitude que par les séries. Employant donc les séries du no 5 ci-dessus, on aura
![{\displaystyle {\begin{aligned}\int (d\Omega )=&-{\frac {\mathrm {T} '}{n'}}\left[\left[{\overline {r}},{\overline {r}}'\ \right]+\left(\left[{\overline {r}},{\overline {r}}'\ \right]_{1}-{\frac {\overline {r}}{{\overline {r}}'^{2}}}\right)\cos \left({\overline {p}}-{\overline {p}}'\ \right)\right.\\&\quad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+\left[{\overline {r}},{\overline {r}}'\ \right]_{2}\cos 2\left({\overline {p}}-{\overline {p}}'\ \right)+\ldots {\Bigr ]}\\&-{\frac {\mathrm {T} ''}{n''}}\left[\left[{\overline {r}},{\overline {r}}''\right]+\left(\left[{\overline {r}},{\overline {r}}''\right]_{1}-{\frac {\overline {r}}{{\overline {r}}''^{2}}}\right)\cos \left({\overline {p}}-{\overline {p}}''\ \right)\right.\\&\quad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+\left[{\overline {r}},{\overline {r}}''\right]_{2}\cos 2\left({\overline {p}}-{\overline {p}}''\ \right)+\ldots {\Bigr ]}\\&-\ldots +\chi .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7bb0be49c7e1fafc7144dbd43e15bb4497cd2920)
et de là
![{\displaystyle {\begin{aligned}\int d\,{\overline {p}}\int (d\Omega )=&-{\frac {\mathrm {T} '}{n'^{2}}}\left[\left(\left[{\overline {r}},{\overline {r}}'\ \right]_{1}-{\frac {\overline {r}}{{\overline {r}}'^{2}\ }}\right)\sin \left({\overline {p}}-{\overline {p}}'\ \right)\right.\\&\qquad \qquad \qquad \qquad \quad \left.+{\frac {1}{2}}\left[{\overline {r}},{\overline {r}}'\ \right]_{2}\sin 2\left({\overline {p}}-{\overline {p}}'\ \right)+\ldots \right]\\&-{\frac {\mathrm {T} ''}{n''^{2}}}\left[\left(\left[{\overline {r}},{\overline {r}}''\right]_{1}-{\frac {\overline {r}}{{\overline {r}}''^{2}}}\right)\sin \left({\overline {p}}-{\overline {p}}''\ \right)\right.\\&\qquad \qquad \qquad \qquad \quad \left.+{\frac {1}{2}}\left[{\overline {r}},{\overline {r}}''\right]_{2}\sin 2\left({\overline {p}}-{\overline {p}}''\ \right)+\ldots \right]\\&-\ldots +\left(\chi -{\frac {\mathrm {T} '}{n'}}\left[{\overline {r}},{\overline {r}}'\right]-{\frac {\mathrm {T} ''}{n''}}\left[{\overline {r}},{\overline {r}}''\right]-\ldots \right){\overline {p}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d244d1152c6d6eec4a2f0a4235af993e8b1174cf)
Or il faut (8) que la valeur de
c’est-à-dire
![{\displaystyle \Sigma +3{\overline {r}}\int d\,{\overline {p}}\int (d\Omega ),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/705c77a84933691b59f8e1a55b6403afa17a3040)
ne contienne aucun terme proportionnel à
égalant donc à zéro la somme des coefficients de ces sortes de termes dans
et dans
![{\displaystyle 3{\overline {r}}\int d\,{\overline {p}}\int (d\Omega ),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de34fbdef4de98730f8aa4565d210585af530722)
on aura l’équation
![{\displaystyle {\begin{aligned}&-2{\overline {r}}^{2}\left(\mathrm {T} '{\frac {d\left[{\overline {r}},{\overline {r}}'\right]}{d\,{\overline {r}}}}+\mathrm {T} ''{\frac {d\left[{\overline {r}},{\overline {r}}''\right]}{d\,{\overline {r}}}}+\ldots \right)\\&+3{\overline {r}}\left(\chi -{\frac {\mathrm {T} '}{n'}}\left[{\overline {r}},{\overline {r}}'\right]-{\frac {\mathrm {T} ''}{n''}}\left[{\overline {r}},{\overline {r}}''\right]-\ldots \right)=0,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4678891cf12afcf673466b1f60eff4770f71f00c)
d’où l’on tire
![{\displaystyle \chi =\mathrm {T} '\left({\frac {\left[{\overline {r}},{\overline {r}}'\right]}{n'}}+{\frac {2{\overline {r}}}{3}}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]}{d\,{\overline {r}}}}\right)+\mathrm {T} ''\left({\frac {\left[{\overline {r}},{\overline {r}}''\right]}{n''}}+{\frac {2{\overline {r}}}{3}}{\frac {d\left[{\overline {r}},{\overline {r}}''\right]}{d\,{\overline {r}}}}\right)+\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/684407a4c68237ca69fac66c54f5d1c1c4697476)
Par le moyen de ces substitutions, les expressions de
et de
du no 8 deviendront
![{\displaystyle {\begin{aligned}\rho =&-\mathrm {T} '{\frac {4{\overline {r}}^{2}}{3}}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]}{d\,{\overline {r}}}}-\mathrm {T} ''{\frac {4{\overline {r}}^{2}}{3}}{\frac {d\left[{\overline {r}},{\overline {r}}''\right]}{d\,{\overline {r}}}}-\ldots \\&+2\mathrm {T} '{\frac {{\overline {r}}^{2}\ }{n'}}\left[\left(\left[{\overline {r}},{\overline {r}}'\ \right]_{1}-{\frac {\overline {r}}{{\overline {r}}'^{2}}}\right)\cos({\overline {p}}-{\overline {p}}'\ )\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+\left[{\overline {r}},{\overline {r}}'\ \right]_{2}\cos 2({\overline {p}}-{\overline {p}}')+\ldots {\Bigr ]}\\&+2\mathrm {T} ''{\frac {{\overline {r}}^{2}}{n''}}\left[\left(\left[{\overline {r}},{\overline {r}}''\right]_{1}-{\frac {\overline {r}}{{\overline {r}}''^{2}}}\right)\cos({\overline {p}}-{\overline {p}}'')\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+\left[{\overline {r}},{\overline {r}}''\right]_{2}\cos 2({\overline {p}}-{\overline {p}}'')+\ldots {\Bigr ]}\\&-\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .\\\varpi =&\mathrm {T} '\ {\Bigl [}\left[\left[{\overline {r}},{\overline {r}}'\,\right]_{1}\right]\sin \left({\overline {p}}-{\overline {p}}'\ \right)+\left[\left[{\overline {r}},{\overline {r}}'\ \right]_{2}\right]\sin 2\left({\overline {p}}-{\overline {p}}'\right)+\ldots {\Bigr ]}\\+&\mathrm {T} ''{\Bigl [}\left[\left[{\overline {r}},{\overline {r}}''\right]_{1}\right]\sin \left({\overline {p}}-{\overline {p}}''\right)+\left[\left[{\overline {r}},{\overline {r}}''\right]_{2}\right]\sin 2\left({\overline {p}}-{\overline {p}}''\right)+\ldots {\Bigr ]}\\+&\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/6480fec13e374cb0b084f5247be8302703e026a9)
en supposant, pour abréger,
![{\displaystyle {\begin{aligned}&\left[\left[{\overline {r}},{\overline {r}}'\right]_{1}\right]={\frac {\overline {r}}{{\overline {r}}'^{2}}}\left({\frac {3}{n'^{2}}}+{\frac {2}{n'}}\right)-{\frac {3{\overline {r}}}{n'^{2}}}\left[{\overline {r}},{\overline {r}}'\right]_{1}-{\frac {2{\overline {r}}^{2}}{n'}}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{1}}{d\,{\overline {r}}}},\\&\left[\left[{\overline {r}},{\overline {r}}'\right]_{2}\right]=-{\frac {3{\overline {r}}}{2n'^{2}}}\left[{\overline {r}},{\overline {r}}'\right]_{2}-{\frac {2{\overline {r}}^{2}}{2n'}}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{2}}{d\,{\overline {r}}}},\\&\left[\left[{\overline {r}},{\overline {r}}'\right]_{3}\right]=-{\frac {3{\overline {r}}}{3n'^{2}}}\left[{\overline {r}},{\overline {r}}'\right]_{3}-{\frac {2{\overline {r}}^{2}}{3n'}}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{3}}{d\,{\overline {r}}}},\\&\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/61dcf4dc957b442256491f23e0a0ed6d75b6ceef)
et ainsi des autres expressions semblables.
11. Après avoir vu comment on doit déterminer les variations périodiques de la distance moyenne et du mouvement moyen, passons à la recherche de celles des excentricités et des aphélies.
Celles-ci sont contenues dans les formules du no 41 de la Théorie citée, lesquelles donnent les variations des quantités
et
ou
et
car en nommant
l’excentricité,
la longitude de l’aphélie sur le plan de projection et
sa latitude, nous avons supposé
![{\displaystyle x=\lambda \cos \eta \sin \varphi ,\quad \operatorname {y} =\lambda \cos \eta \cos \varphi .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d458f7d81241ef4797d3a3b3530685cc20a0721)
Or on a vu dans le no 42 qu’en négligeant les carrés des excentricités et des inclinaisons, les formules dont il s’agit se réduisent à celles-ci
![{\displaystyle {\begin{aligned}d\mathrm {M} =&-r^{2}{\frac {d\Omega }{dr}}\cos qdq+{\frac {d\Omega }{dq}}(2r\sin qdq-\cos qdr),\\d\mathrm {N} =&r^{2}{\frac {d\Omega }{dr}}\sin qdq+{\frac {d\Omega }{dq}}(2r\cos qdq+\sin qdr).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f85647a27dedea16138fe0d604eefd016e11777a)
Faisant donc dans ces formules les substitutions du no 4 ci-dessus, et changeant
en
on aura des équations de cette forme
![{\displaystyle dx=\left(\mathrm {X} \sin {\overline {p}}-\mathrm {Y} \cos {\overline {p}}\right)d{\overline {p}},\quad dy=\left(\mathrm {X} \cos {\overline {p}}+\mathrm {Y} \sin {\overline {p}}\right)d{\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/13a00e7b7695528b0afa7a15bd6a7b50fd9434aa)
dans lesquelles
![{\displaystyle {\begin{aligned}\mathrm {X} =&2{\overline {r}}{\frac {d\Omega }{d{\overline {p}}}}+\left(-2{\overline {r}}{\frac {d\Omega }{d{\overline {p}}}}+2{\overline {r}}^{2}{\frac {d^{2}\Omega }{d{\overline {p}}d{\overline {r}}}}\right)\left(x\sin {\overline {p}}+y\cos {\overline {p}}\right)\\&+2{\overline {r}}{\overline {r}}'\ {\frac {d^{2}\Omega }{d{\overline {p}}d{\overline {r}}'}}\left(x'\ \sin {\overline {p}}'\ +y'\ \cos {\overline {p}}'\right)\\&+2{\overline {r}}{\overline {r}}''{\frac {d^{2}\Omega }{d{\overline {p}}d{\overline {r}}''}}\left(x''\sin {\overline {p}}''+y''\cos {\overline {p}}''\right)+\ldots \\&+\left(2{\overline {r}}^{2}{\frac {d\Omega }{d{\overline {r}}}}+4{\overline {r}}{\frac {d^{2}\Omega }{d{\overline {p}}^{2}}}\right)\left(x\cos {\overline {p}}-y\sin {\overline {p}}\right)\\&+4{\overline {r}}{\frac {d^{2}\Omega }{d{\overline {p}}d{\overline {p}}'\,}}\left(x'\,\cos {\overline {p}}'\,-y'\ \sin {\overline {p}}'\ \right)\\&+4{\overline {r}}{\frac {d^{2}\Omega }{d{\overline {p}}d{\overline {p}}''}}\left(x''\cos {\overline {p}}''-y''\sin {\overline {p}}''\right)\\&+\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\\\mathrm {Y} =&{\overline {r}}^{2}{\frac {d\Omega }{d{\overline {r}}}}+{\overline {r}}^{3}{\frac {d^{2}\Omega }{d{\overline {r^{2}}}}}\left(x\sin {\overline {p}}+y\cos {\overline {p}}\right)\\&+{\overline {r}}^{2}{\overline {r}}'\ {\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {r}}'}}\left(x'\ \sin {\overline {p}}'\ +y'\ \cos {\overline {p}}'\right)\\&+{\overline {r}}^{2}{\overline {r}}''{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {r}}''}}\left(x''\sin {\overline {p}}''+y''\cos {\overline {p}}''\right)+\ldots \\&+\left(-3{\overline {r}}{\frac {d\Omega }{d{\overline {p}}}}+2{\overline {r}}^{2}{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}}}\right)\left(x\cos {\overline {p}}-y\sin {\overline {p}}\right)\\&+2{\overline {r}}^{2}{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}'}}\left(x'\,\cos {\overline {p}}'\ -y'\ \sin {\overline {p}}'\ \right)\\&+2{\overline {r}}^{2}{\frac {d^{2}\Omega }{d{\overline {r}}d{\overline {p}}''}}\left(x''\cos {\overline {p}}''-y''\sin {\overline {p}}''\right)\\&+\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e05c0b1d4b590571113add12dc0cf501cf983732)
12. Si maintenant on substitue dans ces expressions de
et
la valeur de
en série (4 et 5), qu’ensuite on développe les différents produits des sinus et cosinus en sinus et cosinus simples, on aura des termes proportionnels à
et
lesquels étant ensuite multipliés par
et
donneront, dans les valeurs de
et
des termes sans sinus ni cosinus ; ce sont ceux que nous avons déterminés à part dans les nos 43 et 49 de la Théorie citée, et auxquels nous avons eu uniquement égard dans les équations différentielles en
et
du no 50, parce que nous faisions alors abstraction des inégalités périodiques. Les autres termes de
et
ne pourront donner dans
et
que des termes affectés de sinus ou cosinus d’angles composés de multiples de
et il faudra maintenant tenir compte aussi de ces termes, pour pouvoir déterminer la partie périodique des valeurs de
et
Désignons, pour abréger, par
la totalité des termes affectés de sinus et cosinus dans la valeur de
![{\displaystyle \mathrm {X} \sin {\overline {p}}-\mathrm {Y} \cos {\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10ca904d4361353dc56f81ddf22dc456ad6a9969)
et par
la totalité des termes pareils dans la valeur de
![{\displaystyle \mathrm {X} \cos {\overline {p}}+\mathrm {Y} \sin {\overline {p}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9d1542d78e89881c73697367706069bcb3de764a)
il faudra donc ajouter respectivement aux valeurs de
des équations différentielles citées les quantités
de sorte que ces quantités formeront maintenant les seconds membres des deux premières équations différentielles dont il s’agit.
Ainsi, puisque (4)
![{\displaystyle d\,{\overline {p}}={\frac {dt}{{\overline {r}}^{\frac {3}{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b6054a8a138c79d36094218407ca1232a04ffd8d)
les équations complètes seront
![{\displaystyle {\begin{aligned}&{\frac {dx}{dt}}-\left((0,1)+(0,2)+\ldots \right)y+[0,1]y'+[0,2]y''+\ldots ={\frac {(\mathrm {X} )}{{\overline {r}}^{\frac {3}{2}}}},\\&{\frac {dy}{dt}}+\left((0,1)+(0,2)+\ldots \right)x-[0,1]x'-[0,2]x''+\ldots ={\frac {(\mathrm {Y} )}{{\overline {r}}^{\frac {3}{2}}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e6b7a494a1ca419e3202351ea3cc64e9d20c5d76)
De même, en dénotant par
ce que deviennent les expressions de
et
lorsqu’on y change
en
et vice versâ, et par
la totalité des termes affectés de sinus et cosinus dans les valeurs de
![{\displaystyle \mathrm {X} '\sin {\overline {p}}'-\mathrm {Y} '\cos {\overline {p}}',\quad \mathrm {X} '\cos {\overline {p}}'+\mathrm {Y} '\sin {\overline {p}}',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/813bc438a904e12a34639f1d31cce2f05b237d2c)
on aura les équations
![{\displaystyle {\begin{aligned}&{\frac {dx'}{dt}}-{\Bigl (}(1,0)+(1,2)+\ldots {\Bigr )}y'+[1,0]y+[1,2]y''+\ldots ={\frac {(\mathrm {X} ')}{{\overline {r}}^{\frac {3}{2}}}},\\&{\frac {dy'}{dt}}+{\Bigl (}(1,0)+(1,2)+\ldots {\Bigr )}x'-[1,0]x+[1,2]x''-\ldots ={\frac {(\mathrm {Y} ')}{{\overline {r}}^{\frac {3}{2}}}},\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0068e83a2a852b9ecad4287ebbe62724156bef1f)
et ainsi de suite.
13. Comme dans les premiers membres de ces équations les variations sont linéaires, et que les seconds membres peuvent être regardés comme des fonctions connues de la variable
l’intégration est toujours possible par les méthodes connues. Dans la Théorie précédente nous avons déjà donné les intégrales complètes pour le cas où les seconds membres seraient nuls ; et ces mêmes intégrales, en y faisant varier les constantes arbitraires, donneront celles des équations dont il s’agit par la méthode indiquée dans le no 27 de la même Théorie.
Désignons, en général, par
les expressions de ![{\displaystyle x,y,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4ad6024a5a4f15ab696de3d6e4866da50c8a315d)
trouvées dans cette Théorie (51), et soient ![{\displaystyle {\overline {x}}+\xi ,{\overline {y}}+\psi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/458e8d0b14943ea109c6f6d7b050670f5da0441c)
les valeurs de
dans les équations ci-dessus, en ayant égard à leurs seconds membres ; il est clair qu’en y substituant ces valeurs, les termes en
s’en iront d’eux-mêmes, et que les transformées ne seront autre chose que les mêmes équations, en y changeant
en
et comme les valeurs de
contiennent déjà toutes les constantes arbitraires nécessaires pour l’intégration complète, il suffira que celles de
satisfassent aux équations d’une manière quelconque.
Or les quantités qui forment les seconds membres des équations dont il s’agit étant dues aux forces perturbatrices, on peut supposer que les valeurs des variables
soient très-petites du même ordre, et, dans cette supposition, les termes qui renferment ces mêmes variables sous une forme finie deviendront très-petits de l’ordre des carrés de ces forces, puisque les coefficients
sont eux-mêmes très-petits de l’ordre des mêmes forces. Ainsi, comme nous négligeons dans les recherches présentes les carrés des forces perturbetrices, et toutes les quantités du même ordre, les équations en
se réduiront simplement à
![{\displaystyle {\frac {d\xi }{dt}}={\frac {(\mathrm {X} )}{{\overline {r}}^{\frac {3}{2}}}},\quad {\frac {d\psi }{dt}}={\frac {(\mathrm {Y} )}{{\overline {r}}^{\frac {3}{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b4ae489838a13ce7209713b8b0b9df6974585321)
d’où l’on tire, en remettant
pour ![{\displaystyle {\frac {dt}{{\overline {r}}^{\frac {3}{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e06de38be3c5ab2a891de215374f3d3e3aebda1)
![{\displaystyle \xi =\int (\mathrm {X} )d\,{\overline {p}},\quad \psi =\int (\mathrm {Y} )d\,{\overline {p}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d42087c3b90ca38accaef92697d31866b606f21a)
Et il en sera de même pour ![{\displaystyle \xi ',\psi ',\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e7c0d941290bf76c36c835216cf1eb803b65f1f6)
14. Si dans les expressions de
du no 11 on fait d’abord abstraction des excentricités, et que par conséquent on y néglige les termes multipliés par
elles se réduisent aux quantités
et comme par la substitutions de la valeur de
en série il ne vient aucun terme proportionnel à
et
on aura simplement, dans ce cas,
![{\displaystyle (\mathrm {X} )=\mathrm {X} \sin {\overline {p}}-\mathrm {Y} \cos {\overline {p}},\quad (\mathrm {Y} )=\mathrm {X} \cos {\overline {p}}+\mathrm {Y} \sin {\overline {p}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47fc7dd6bdbfe1c6dec012829497f654592b3933)
L’intégration de
et
ne présente aucune difficulté, puisqu’on a (6)
![{\displaystyle d\,{\overline {p}}-d\,{\overline {p}}'=n'd\,{\overline {p}},\quad d\,{\overline {p}}-d\,{\overline {p}}''=n''d\,{\overline {p}},\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/82eba5b9dc4aa538a8c619a2474593daceda6d4a)
mais on peut la simplilier beaucoup en remarquant, en général, que si
est un terme quelconque de
et
le terme correspondant de
étant un angle tel que
![{\displaystyle d\Pi =\mu d\,{\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b476b81d484beb486e888b46448215ad3c3e95bd)
ces termes donneront dans l’intégrale de
![{\displaystyle (\mathrm {X} )d\,{\overline {p}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0004a83e92dbb2bc528ca5d48e22a7fd0ab44fe5)
deux termes de la forme
![{\displaystyle \mathrm {P} \cos \Pi \sin {\overline {p}}-\mathrm {Q} \sin \Pi \cos {\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/467840f1de686d68e99f863a0f5b8266fb3b3dae)
et dans l’intégrale de
les termes correspondants
![{\displaystyle \mathrm {P} \cos \Pi \cos {\overline {p}}+\mathrm {Q} \sin \Pi \sin {\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/abbbab66f6ab9380ff0c7f06af58374df4cbd0f7)
En effet, en différentiant et comparant les termes analogues, on parvient aux équations
![{\displaystyle \mathrm {P-\mu Q=-N,\quad -\mu P+Q=M} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff39351c64b7b774da093de66192f1f850f3f6fc)
Or dans le cas présent on a
![{\displaystyle {\begin{aligned}\mathrm {X} =&2{\overline {r}}{\frac {d\Omega }{d\,{\overline {p}}}}\\=&2\mathrm {T} '{\overline {r}}\left[\left(\left[{\overline {r}},{\overline {r}}'\ \right]_{1}-{\frac {\overline {r}}{{\overline {r}}'^{2}}}\right)\sin({\overline {p}}-{\overline {p}}'\ )\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+2\left[{\overline {r}},{\overline {r}}'\ \right]_{2}\sin 2({\overline {p}}-{\overline {p}}')+\ldots {\Bigr ]}\\=&2\mathrm {T} ''{\overline {r}}\left[\left(\left[{\overline {r}},{\overline {r}}''\right]_{1}-{\frac {\overline {r}}{{\overline {r}}''^{2}}}\right)\sin({\overline {p}}-{\overline {p}}'')\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+2\left[{\overline {r}},{\overline {r}}''\right]_{2}\sin 2({\overline {p}}-{\overline {p}}'')+\ldots {\Bigr ]}\\&+\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\\\mathrm {Y} =&{\overline {r}}^{2}{\frac {d\Omega }{d\,{\overline {r}}}}\\=&-\mathrm {T} '{\overline {r}}^{2}\left[{\frac {d\left[{\overline {r}},{\overline {r}}'\ \right]}{d\,{\overline {r}}}}+\left({\frac {d\left[{\overline {r}},{\overline {r}}'\ \right]_{1}}{d\,{\overline {r}}}}-{\frac {1}{{\overline {r}}'^{2}}}\right)\cos({\overline {p}}-{\overline {p}}'\ )\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+{\frac {d\left[{\overline {r}},{\overline {r}}'\ \right]_{2}}{d\,{\overline {r}}}}\cos 2({\overline {p}}-{\overline {p}}')+\ldots {\Bigr ]}\\=&-\mathrm {T} ''{\overline {r}}^{2}\left[{\frac {d\left[{\overline {r}},{\overline {r}}''\right]}{d\,{\overline {r}}}}+\left({\frac {d\left[{\overline {r}},{\overline {r}}''\right]_{1}}{d\,{\overline {r}}}}-{\frac {1}{{\overline {r}}''^{2}}}\right)\cos({\overline {p}}-{\overline {p}}'')\right.\\&\qquad \qquad \qquad \qquad \qquad \qquad {\Bigl .}+{\frac {d\left[{\overline {r}},{\overline {r}}''\right]_{2}}{d\,{\overline {r}}}}\cos 2({\overline {p}}-{\overline {p}}'')+\ldots {\Bigr ]}\\&-\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/7bd2a07f4f0b0c6111f5d99f39024d33ecd11f6c)
Si donc on fait successivement
![{\displaystyle \mu =0,\ \ n',\ \ 2n',\ldots \ \ n'',\ \ 2n'',\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec27c1621d3b04273f0e3e7d8ee09b5ae09dedbf)
et qu’on prenne pour
et
les coefficients respectifs des sinus et cosinus dans les séries précédentes, on aura pour
![{\displaystyle \int (X)d\,{\overline {p}}\quad {\text{et}}\quad \int (Y)d\,{\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/623c09ee0bb44bacb1183e5ca61877df73708b63)
c’est-à-dire pour
![{\displaystyle \xi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e0b461aaf61091abd5d2c808931c48b8ff9647db)
et
![{\displaystyle \psi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bfc4043b55bade492740e58cba74198873db1464)
des expressions de cette forme
![{\displaystyle \xi =\Xi \sin {\overline {p}}-\Psi \cos {\overline {p}},\quad \psi =\Xi \cos {\overline {p}}+\Psi \sin {\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e88d89cd7c9a8326695103789cf0a4388dfd2c65)
dans lesquelles
![{\displaystyle {\begin{aligned}\Xi =&\mathrm {T} '\ {\Bigl [}\left[\left[{\overline {r}},{\overline {r}}'\ \right]\right.\left[\left[{\overline {r}},{\overline {r}}'\,\ \right]\right._{1}\cos \left({\overline {p}}-{\overline {p}}'\ \right)+\left[\left[{\overline {r}},{\overline {r}}'\ \right]_{2}\right.\cos 2\left({\overline {p}}-{\overline {p}}'\right)+\ldots {\Bigr ]}\\+&\mathrm {T} ''{\Bigl [}\left[\left[{\overline {r}},{\overline {r}}''\,\right]\right.\left[\left[{\overline {r}},{\overline {r}}''\,\right]\right._{1}\cos \left({\overline {p}}-{\overline {p}}''\right)+\left[\left[{\overline {r}},{\overline {r}}''\right]_{2}\right.\cos 2\left({\overline {p}}-{\overline {p}}''\right)+\ldots {\Bigr ]}\\+&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\\psi =&\mathrm {T} '\,\ {\Bigl [}\left.\left[{\overline {r}},{\overline {r}}'\,\ \right]\right]_{1}\sin \left({\overline {p}}-{\overline {p}}'\ \ \right)+\left.\left[{\overline {r}},{\overline {r}}'\ \right]_{2}\right]\sin 2\left({\overline {p}}-{\overline {p}}'\ \right)+\ldots {\Bigr ]}\\+&\mathrm {T} ''\ {\Bigl [}\left.\left[{\overline {r}},{\overline {r}}''\,\right]\right]_{1}\sin \left({\overline {p}}-{\overline {p}}''\ \right)+\left.\left[{\overline {r}},{\overline {r}}''\ \right]_{2}\right]\sin 2\left({\overline {p}}-{\overline {p}}''\right)+\ldots {\Bigr ]}\\+&\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/0e9ca84d0a27fdb11637bc7ca68e62a880148ff2)
en supposant
![{\displaystyle {\begin{aligned}&\left[\left[{\overline {r}},{\overline {r}}'\right]_{1}\right.={\overline {r}}^{2}{\frac {\left[{\overline {r}},{\overline {r}}'\right]}{dr}},\\&\left[\left[{\overline {r}},{\overline {r}}'\right]_{1}\right.={\frac {{\overline {r}}^{2}{\cfrac {d\left[{\overline {r}},{\overline {r}}'\right]_{1}}{d\,{\overline {r}}}}+2n'{\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{1}-{\frac {{\overline {r}}^{2}}{{\overline {r}}'^{2}}}(1+2n')}{1-n'^{2}}},\\&\left[\left[{\overline {r}},{\overline {r}}'\right]_{2}\right.={\frac {{\overline {r}}^{2}{\cfrac {d\left[{\overline {r}},{\overline {r}}'\right]_{2}}{d\,{\overline {r}}}}+2.4n'{\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{2}}{1-4n'^{2}}},\\&\left[\left[{\overline {r}},{\overline {r}}'\right]_{3}\right.={\frac {{\overline {r}}^{2}{\cfrac {d\left[{\overline {r}},{\overline {r}}'\right]_{3}}{d\,{\overline {r}}}}+2.9n'{\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{3}}{1-9n'^{2}}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\\&\left.\left[{\overline {r}},{\overline {r}}'\right]_{1}\right]={\frac {n'{\overline {r}}^{2}{\cfrac {d\left[{\overline {r}},{\overline {r}}'\right]_{1}}{d\,{\overline {r}}}}+2{\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{1}-{\frac {{\overline {r}}^{2}}{{\overline {r}}'^{2}}}(n'+2)}{1-n'^{2}}},\\&\left.\left[{\overline {r}},{\overline {r}}'\right]_{2}\right]=2{\frac {n'{\overline {r}}^{2}{\cfrac {d\left[{\overline {r}},{\overline {r}}'\right]_{2}}{d\,{\overline {r}}}}+2{\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{2}}{1-4n'^{2}}},\\&\left.\left[{\overline {r}},{\overline {r}}'\right]_{3}\right]=3{\frac {n'{\overline {r}}^{2}{\cfrac {d\left[{\overline {r}},{\overline {r}}'\right]_{3}}{d\,{\overline {r}}}}+2{\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{3}}{1-9n'^{2}}},\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad8bf26bff11d1ca090dec03370c0a816cbe96bb)
et ainsi des autres fonctions semblables.
On pourra avec la même facilité tenir compte des termes multipliés par les variables
dans les expressions de
et
pour cela on y changera d’abord ces variables en
et l’on y substituera leurs valeurs trouvées dans la Théorie des variations séculaires on n’aura ainsi à intégrer que des termes en sinus et cosinus, et l’on pourra même simplifier beaucoup le calcul par des remarques semblables à celles que nous avons faites (7).
15. Il ne reste plus qu’à déterminer les inégalités des nœuds et des inclinaisons, lesquelles dépendent des équations différentielles du no 39 de la Théorie citée ; mais nous réduirons auparavant ces mêmes équations à une forme plus simple et plus générale à quelques égards.
En supposant, comme nous l’avons fait,
![{\displaystyle s=\mathrm {\frac {P}{R}} =\theta \sin \omega ,\quad u=\mathrm {\frac {Q}{R}} =\theta \cos \omega ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bea06ae1f23a7303439caf3f890c7050a5d5984f)
où
représente la tangente de l’inclinaison et
la longitude du nœud ascendant comptée depuis un point fixe, on a
![{\displaystyle ds={\frac {d\mathrm {P} }{\mathrm {R} }}-{\frac {\mathrm {P} d\mathrm {R} }{\mathrm {R} ^{2}}},\quad du={\frac {dQ}{R}}-{\frac {\mathrm {Q} d\mathrm {R} }{\mathrm {R} ^{2}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f8ba2e86bc0cd5e786812f8903c1a2d24bbb4d0)
Or
![{\displaystyle {\begin{aligned}d\mathrm {P} =&\mathrm {T} '\left({\frac {1}{\rho '^{3}}}-{\frac {1}{\sigma '^{3}}}\right)(y'z-yz')dt+\mathrm {T} ''\left({\frac {1}{\rho ''^{3}}}-{\frac {1}{\sigma ''^{3}}}\right)(y''z-yz'')dt+\ldots ,\\d\mathrm {Q} =&\mathrm {T} '\left({\frac {1}{\rho '^{3}}}-{\frac {1}{\sigma '^{3}}}\right)(x'z-xz')dt+\mathrm {T} ''\left({\frac {1}{\rho ''^{3}}}-{\frac {1}{\sigma ''^{3}}}\right)(x''z-xz'')dt+\ldots ,\\d\mathrm {R} =&\mathrm {T} '\left({\frac {1}{\rho '^{3}}}-{\frac {1}{\sigma '^{3}}}\right)(x'y-xy')dt+\mathrm {T} ''\left({\frac {1}{\rho ''^{3}}}-{\frac {1}{\sigma ''^{3}}}\right)(x''y-xy'')dt+\ldots \,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/62bf5b08bc82cd1cee6c7dd59227ebf444aee506)
mais
![{\displaystyle z={\frac {\mathrm {Q} y-\mathrm {P} x}{\mathrm {R} }},\quad z'={\frac {\mathrm {Q} 'y'-\mathrm {P} 'x'}{\mathrm {R} '}},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2646102b8519dca77a14f92d973c7685a3a7e75a)
donc
![{\displaystyle {\begin{aligned}y'z-yz'=&\mathrm {\left({\frac {Q}{R}}-\mathrm {\frac {Q'}{R'}} \right)} yy'-\mathrm {\frac {P}{R}} xy'+\mathrm {\frac {P'}{R'}} x'y,\\x'z-xz'=&-\mathrm {\left({\frac {P}{R}}-\mathrm {\frac {P'}{R'}} \right)} xx'+\mathrm {\frac {Q}{R}} x'y-\mathrm {\frac {Q'}{R'}} xy'\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec392df766ab9975420ba4aa724a9f46a8bb6376)
faisant ces substitutions dans les valeurs de
![{\displaystyle ds}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5f0fb36e4308227d3e4a1f809c2571ec02527100)
et de
![{\displaystyle du,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6d009e975bd6131b0512861d72dbcea41a712f0)
et supposant, pour abréger,
![{\displaystyle {\begin{aligned}\mathrm {Z} =&\mathrm {T} '\ \left({\frac {1}{\rho '^{3}\,}}-{\frac {1}{\sigma '^{3}\,}}\right)\left[(u-u'\ )y'\,-(s-s'\ )x'\ \right]\\+&\mathrm {T} ''\left({\frac {1}{\rho ''^{3}}}-{\frac {1}{\sigma ''^{3}}}\right)\left[(u-u'')y''-(s-s'')x''\right]\\+&\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/3c2d9587d70149b90720e2b4d844f2fb55c25fe4)
on aura
![{\displaystyle ds=\mathrm {Z} y{\frac {dt}{\mathrm {R} }},\quad du=\mathrm {Z} x{\frac {dt}{\mathrm {R} }}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e2b9b30ef3c4c853489cbbc5de66ad1ac784271f)
et il n’y aura plus qu’à mettre pour
leurs valeurs
![{\displaystyle r\cos q,\quad r\sin q,\quad r'\cos q',\quad r'\sin q',\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/05a7d5aee5fe342f90e849b17974064c395d3826)
et pour
![{\displaystyle {\frac {1}{\rho '^{3}}}-{\frac {1}{\sigma '^{3}}},\quad {\frac {1}{\rho ''^{3}}}-{\frac {1}{\sigma ''^{3}}},\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb81ba29444b28efa140631d0ae53c44b00f5fc8)
les expressions en séries données dans le numéro cité, ensuite substituer partout les valeurs de
en ![{\displaystyle {\overline {p}},{\overline {p}}',\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4661d2aceb921acb85eb892f11353aa952697c42)
En négligeant les quantités du second ordre (
étant regardées comme du premier, ainsi que
) il suffira de changer
en
en
et
en
ainsi, à cause de
![{\displaystyle d{\overline {p}}={\frac {d\varepsilon }{{\overline {r}}^{\frac {3}{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ae7a39cd8cc6420f1b02020888872a2ff561c75a)
on aura alors simplement
![{\displaystyle ds={\overline {r}}^{2}\mathrm {Z} \sin {\overline {p}}d{\overline {p}},\quad du={\overline {r}}^{2}\mathrm {Z} \cos {\overline {p}}d{\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/65807f812b91de8e80cd52982ae3c13f94de5eb8)
et
![{\displaystyle {\begin{aligned}\mathrm {Z} &=\mathrm {T} '\ {\overline {r}}'\,\left[(u-u'\,)\sin {\overline {p}}'-(s-s'\ )\cos {\overline {p}}'\ \right]\\&\quad \times \left[{\frac {1}{{\overline {r}}'^{2}}}-({\overline {r}},{\overline {r}}'\ )-({\overline {r}},{\overline {r}}'\ )_{1}\cos({\overline {p}}-{\overline {p}}'\,)-({\overline {r}},{\overline {r}}'\ )_{2}\cos 2({\overline {p}}-{\overline {p}}'\ )-\ldots \right]\\&+\mathrm {T} ''{\overline {r}}''\left[(u-u'')\sin {\overline {p}}''-(s-s'')\cos {\overline {p}}''\right]\\&\quad \times \left[{\frac {1}{{\overline {r}}''^{2}}}-({\overline {r}},{\overline {r}}'')-({\overline {r}},{\overline {r}}'')_{1}\cos({\overline {p}}-{\overline {p}}'')-({\overline {r}},{\overline {r}}'')_{2}\cos 2({\overline {p}}-{\overline {p}}'')-\ldots \right]\\&+\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/dda2be7377e0ea43aa3f9042b04eb86f93b8c8b2)
les quantités
étant les coefficients de la série
qui représente le radical
![{\displaystyle \left[{\overline {r}}^{2}-2{\overline {r}}{\overline {r}}'\cos \left({\overline {p}}-{\overline {p}}'\right)+{\overline {r}}'^{2}\right]^{-{\frac {3}{2}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f41257ade3f0eed51dc55dee8ba68061618a4d1b)
et ainsi des autres fonctions semblables.
16. Si dans cette expression de
on développe les produits des sinus et cosinus, il viendra des termes proportionnels à
et
lesquels donneront par conséquent, dans les valeurs de
et
des termes exempts de sinus et cosinus. Ce seront les mêmes que nous avons trouvés dans l’endroit cité de la Théorie des variations séculaires, et qui nous ont donné les équations pour les variations séculaires de
et
et comme l’analyse précédente est indépendante de la condition
que nous avions employée dans le même endroit, il s’ensuit que les équations dont il s’agit auraient également lieu quand même le grand axe de l’orbite serait aussi sujet à des variations séculaires ; qu’ainsi notre Théorie des nœuds et des inclinaisons des orbites planétaires subsiste indépendamment de l’inaltérabilité des distances moyennes.
À l’égard des autres termes, ils ne pourront donner que des sinus et \inftysinus dans les valeurs de
et
nous les avons négligés dans les équations du no 40 de la Théorie citée, parce qu’il n’était question alors que des variations séculaires ; mais, pour avoir maintenant les valeurs complètes de
et
il faudra ajouter ces termes aux équations dont il s’agit.
Désignons par
la totalité des termes en sinus et cosinus de la valeur de
et par
la totalité des termes pareils dans la valeur de
il faudra donc ajouter respectivement aux valeurs de ![{\displaystyle {\frac {ds}{dt}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3befa0f22887a83e01c0c4302c4d03a97bd94907)
les quantités
de sorte que ces quantités formeront les seconds inembres des équations du numéro cité, lesquelles deviendront de cette manière
![{\displaystyle {\begin{aligned}&{\frac {ds}{dt}}+(0,1)(u-u')+(0,2)(u-u'')+\ldots ={\frac {(\mathrm {S} )}{{\overline {r}}^{\frac {3}{2}}}},\\&{\frac {du}{dt}}-(0,1)(s-s'\,)-(0,2)(s\,-s''\,)-\ldots ={\frac {(\mathrm {U} )}{{\overline {r}}^{\frac {3}{2}}}},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b46e0951397244c66bba1dc579209d7a78f07ac7)
et ainsi des autres équations semblables.
17. On appliquera à ces équations ce que nous avons dit (13) relativement aux équations en
et
et l’on en conclura que si
désignent les valeurs de
déterminées dans la Théorie des variations séculaires, et qu’on suppose
![{\displaystyle s={\overline {s}}+\sigma ,\quad u={\overline {u}}+\nu ,\quad s'={\overline {s}}'+\sigma ',\quad u'={\overline {u}}'+\nu ',\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fbe0992dc33de35883fa7c088d9d66e0607fad36)
on aura avec assez d’exactitude
![{\displaystyle \sigma =\int (\mathrm {S} )d{\overline {p}},\quad \nu =\int (\mathrm {U} )d{\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4a0818f39b0020b272f67d55b1e7f549028043d5)
et ainsi de ![{\displaystyle \sigma ',\nu ',\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5eb98731d3e95abdacc7f9947708344cba22e7f3)
Au reste, comme la quantité
du no 15 ne contient que des termes déjà multipliés par
il s’ensuit qu’en faisant abstraction, dans l’effet des forces perturbatrices, des inclinaisons des orbites, les quantités dont il s’agit seront nulles ; par conséquent on aura aussi dans cette hypothèse
et
nuls, et de là
![{\displaystyle \sigma =0,\quad \nu =0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8334aff00bd5769a659a2e8fc8ffb27959bbd54)
et de même
![{\displaystyle \sigma '=0,\quad \nu '=0,\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/36aea93495e37403cc8b2ae355e7c1b4115eb65e)
Lorsqu’on voudra pousser l’exactitude plus loin, ce qui cependant ne paraît guère nécessaire dans la Théorie du système planétaire, il n’y aura qu’à substituer dans
pour
les valeurs déjà connues ![{\displaystyle {\overline {s}},{\overline {u}},{\overline {s}}',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/14b3a45efd35c5abbc533e3d1eba4e2c828880be)
et intégrer les différents termes qui en proviendront, suivant les méthodes ordinaires. Sur quoi voyez ce que nous avons dit (7).
18. La pratique ordinaire des Astronomes dans le calcul des Planètes est de chercher d’abord le lieu dans l’orbite, et de le réduire ensuite à l’écliptique. Ce lieu dépend de la longitude moyenne
et des éléments de l’orbite, c’est-à-dire du demi-grand axe
de l’excentricité
et de la longitude de l’aphélie ; or, quoique dans les orbites peu inclinées, comme celles des Planètes, la longitude de l’aphélie dans l’orbite soit presque la même que sa longitude dans l’écliptique, que nous avons désignée par
si l’on voulait néanmoins la déterminer rigoureusement par nos formules, il faudrait faire attention que,
étant la longitude du nœud dans le plan de projection, sa longitude doit être exprimée par
étant l’angle d’inclinaison ; car il est évident que, tandis que le nœud avance de l’angle
dans le plan de projection, il n’avance dans le plan de l’orbite que de l’angle
qui est le côté adjacent à l’angle
dans le triangle rectangle dont
est l’hypoténuse.
Ainsi, nommant
la longitude
du nœud dans l’orbite,
la longitude de l’aphélie dans la même orbite, et
la latitude de cet aphélie comme dans le no 9 de la Théorie des variations séculaires, et considérant le triangle sphérique rectangle dont
est l’hypoténuse,
les deux côtés, et
l’angle opposé au côté
on aura par les formules connues
![{\displaystyle \sin(\Phi -\Omega )={\frac {\sin \eta }{\sin i}},\quad \cos(\Phi -\Omega )=\cos(\varphi -\omega )\cos \eta \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3c2980e2f716d9f371f3c9c5b959c893062e7ace)
faisant les substitutions du numéro cité, on aura
![{\displaystyle \lambda \sin(\Phi -\Omega )={\frac {\mathrm {L} }{\sin i}},\quad \lambda \cos(\Phi -\Omega )=\mathrm {M} \sin \omega +\mathrm {N} \cos \omega .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/43a81388bd22363046a8e2dc7229c3bbb511e23d)
Or
![{\displaystyle \mathrm {L=(M\cos \omega -N\sin \omega )} \operatorname {tang} i\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cc636278437c60212ff1e9fb05461656861d5af4)
donc, puisque
![{\displaystyle {\frac {1}{\sin i}}={\frac {1}{\operatorname {tang} i}}+\operatorname {tang} {\frac {i}{2}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7670f78517faaa5c974df5f57bd8d289c53441e5)
la première équation deviendra
![{\displaystyle \lambda \sin(\Phi -\Omega )=\mathrm {M} \cos \omega -\mathrm {N} \sin \omega +\mathrm {L} \operatorname {tang} {\frac {i}{2}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1879cd71f37066f6795e343f9ce97adb8254516a)
laquelle, étant combinée avec la seconde, donnera
![{\displaystyle {\begin{aligned}\lambda \sin \Phi =&\mathrm {M} \cos(\omega -\Omega )-\mathrm {N} \sin(\omega -\Omega )+\mathrm {L} \operatorname {tang} {\frac {i}{2}}\cos \Omega ,\\\lambda \cos \Phi =&\mathrm {M} \sin(\omega -\Omega )+\mathrm {N} \cos(\omega -\Omega )-\mathrm {L} \operatorname {tang} {\frac {i}{2}}\sin \Omega .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/790f1cdd298189f08378648e70e8fba07054150c)
En ajoutant ensemble les carrés de ces équations, on aura
![{\displaystyle \lambda ^{2}=\mathrm {M^{2}+N^{2}+L^{2}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6b10fa6a0075d6cac82a3aad001a2d57913e9c1e)
comme cela doit être d’après nos formules primitives, et en les divisant l’une par l’autre on aura la valeur de
en quantités toutes connues, puisque les angles
et
sont donnés par les formules
![{\displaystyle s=\operatorname {tang} i\sin \omega ,\quad u=\operatorname {tang} i\cos \omega .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b4dd803941158a579eb9c3f6bf8633c6786ea079)
Si l’on voulait au contraire déterminer les valeurs de
en ![{\displaystyle \lambda ,\Phi ,\Omega ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f8a51636afd2830d4844de7a9619b7ec7de94f04)
ce qui pourrait être utile dans quelques occasions, on trouverait par les formules précédentes
![{\displaystyle {\begin{aligned}\mathrm {N} \ =&\lambda \left[\cos(\Phi -\Omega )\cos \omega -\sin(\Phi -\Omega )\sin \omega \cos i\right],\\\mathrm {M} =&\lambda \left[\cos(\Phi -\Omega )\sin \omega +\sin(\Phi -\Omega )\cos \omega \cos i\right],\\\mathrm {L} \ =&\lambda \sin(\Phi -\Omega )\sin i.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8afadbbdd90b8311cfe06e521a7331749c4c26e7)
19. Comme l’objet final des observations et des calculs astronomiques relativement aux Planètes se réduit toujours à la détermination de leurs lieux rapportés à l’écliptique, et que cette détermination ne dépend que de trois éléments, de la distance accourcie au rayon vecteur de l’orbite projetée, de la longitude vraie ou angle de ce rayon avec une ligne fixe, et de la latitude ou angle du rayon vecteur de l’orbite vraie avec celui de l’orbite projetée, il suffira de connaître les altérations que ces éléments doivent subir par l’effet des variations périodiques que nous avons déterminées.
Cet effet consiste à changer les éléments de l’orbite
![{\displaystyle {\overline {r}},\quad {\overline {x}},\quad {\overline {y}},\quad {\overline {s}},\quad {\overline {u}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/996aebec775d47cc736f6a279aaa6a428a9b2ed1)
en
![{\displaystyle {\overline {r}}+\rho ,\quad {\overline {x}}+\xi ,\quad {\overline {y}}+\psi ,\quad {\overline {s}}+\sigma ,\quad {\overline {u}}+\nu }](https://wikimedia.org/api/rest_v1/media/math/render/svg/a846e6b74e9a0b8b5fdf5cab9516d624ee4432c0)
et la longitude moyenne
en
ainsi il n’y aura qu’à faire ces changements dans les expressions connues de la longitude vraie, du rayon vecteur et de la latitude, calculées pour une ellipse invariable. Par conséquent, si
est la longitude vraie,
le rayon vecteur,
Ia tangente de
la latitude, ces quantités étant exprimées par l’angle
![{\displaystyle {\overline {p}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3884eccbf972d435831da42dd4154e192d550e16)
du mouvement moyen, et par les éléments
![{\displaystyle {\overline {r}},{\overline {x}},{\overline {y}},{\overline {s}},{\overline {u}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c8eaff378dae5e72381b613d91978601cae55590)
supposés constants, on aura, en négligeant les carrés des corrections
![{\displaystyle \rho ,\xi ,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/b320f19c02f5314a43f2cbc6a7f98d5ef201d910)
dues aux forces perturbatrices,
Correction de la longitude
![{\displaystyle {\frac {dq}{d{\overline {p}}}}\varpi +{\frac {dq}{d{\overline {x}}}}\xi +{\frac {dq}{d{\overline {y}}}}\psi +{\frac {dq}{d{\overline {s}}}}\sigma +{\frac {dq}{d{\overline {u}}}}\nu ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5fd53af67ab8edbeb8a2a14e3a9811822e776a8)
Correction du rayon vecteur
![{\displaystyle {\frac {dr}{d{\overline {p}}}}\varpi +{\frac {dr}{d{\overline {x}}}}\xi +{\frac {dr}{d{\overline {y}}}}\psi +{\frac {dr}{d{\overline {s}}}}\sigma +{\frac {dr}{d{\overline {u}}}}\nu +{\frac {dr}{d{\overline {r}}}}\rho ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e53fb025b705fff449b9c015e57be24f481533d0)
Correction de la tangente de latitude
![{\displaystyle {\frac {dl}{d{\overline {p}}}}\varpi +{\frac {dl}{d{\overline {x}}}}\xi +{\frac {dl}{d{\overline {y}}}}\psi +{\frac {dl}{d{\overline {s}}}}\sigma +{\frac {dl}{d{\overline {u}}}}\nu .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/57935ecff125e3ed1b66db9078ccc8b902663716)
20. Les expressions de
en
peuvent se déduire immédiatement des formules de la Théorie citée (29 et suivants). Car, en n’ayant égard qu’aux quantités très-petites du premier et du second ordre, ce qui suffit dans la recherche présente, on aura
![{\displaystyle q=p+(\mathrm {B} )\cos {\overline {p}}-(\mathrm {C} )\sin {\overline {p}}+(\mathrm {D} )\cos 2{\overline {p}}-(\mathrm {E} )\sin 2{\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e3175f5254b60d60bddef494c1d1a5b008e06ae)
où
![{\displaystyle {\begin{aligned}(\mathrm {B} )=&\beta -{\frac {1}{4}}\beta \varepsilon +{\frac {1}{4}}\gamma \delta ,\\(\mathrm {C} )=&\gamma +{\frac {1}{4}}\beta \delta +{\frac {1}{4}}\gamma \varepsilon ,\\(\mathrm {D} )=&{\frac {1}{2}}\delta -\beta \gamma ,\\(\mathrm {E} )=&{\frac {1}{2}}\varepsilon +{\frac {2}{1}}\beta ^{2}-{\frac {1}{2}}\gamma ^{2}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9fa4ac6a58a38cfef42396ab61711d1a5310ccba)
de sorte qu’en faisant les substitutions du no 2, et changeant
en
on aura
![{\displaystyle q={\overline {p}}+2{\overline {x}}\cos {\overline {p}}-2{\overline {y}}\sin {\overline {p}}-{\frac {5{\overline {x}}{\overline {y}}-{\overline {s}}{\overline {u}}}{2}}\cos 2{\overline {p}}-{\frac {5\left({\overline {x}}^{2}-{\overline {y}}^{2}\right)+{\overline {u}}^{2}-{\overline {s}}^{2}}{4}}\sin 2{\overline {p}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7cfe64c90474e589650511c827edb91727382e92)
Pour le rayon vecteur
de l’orbite projetée, on aura d’abord cette série en
![{\displaystyle {\overline {r}}\left[1+{\overline {x}}\sin q+{\overline {y}}\cos q-\left({\frac {{\overline {x}}^{2}-{\overline {y}}^{2}}{2}}+{\frac {{\overline {s}}^{2}-{\overline {u}}^{2}}{4}}\right)\cos 2q\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/527ed5a1f9b96d3fd3e14b0566073a016f5647db)
![{\displaystyle \left.+\left({\overline {x}}{\overline {y}}+{\frac {{\overline {s}}{\overline {u}}}{2}}\right)\sin 2q-{\frac {{\overline {x}}^{2}+{\overline {y}}^{2}}{2}}-{\frac {{\overline {s}}^{2}+{\overline {u}}^{2}}{4}}\right],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eb29b7ac4dc596746c8a0d35f3eeb45b79f6499a)
laquelle par la substitution de la valeur précédente de
en
donnera
![{\displaystyle r={\overline {r}}\left[1+{\overline {x}}\sin {\overline {p}}+{\overline {y}}\cos {\overline {p}}+\left({\frac {{\overline {x}}^{2}-{\overline {y}}^{2}}{2}}-{\frac {{\overline {s}}^{2}-{\overline {u}}^{2}}{4}}\right)\cos 2{\overline {p}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/074034858fe80eb30cdc3b3ae1ed37a9eb3f8624)
![{\displaystyle \left.-\left({\overline {x}}{\overline {y}}+{\frac {{\overline {s}}{\overline {u}}}{2}}\right)\sin 2{\overline {p}}+{\frac {{\overline {x}}^{2}+{\overline {y}}^{2}}{2}}-{\frac {{\overline {s}}^{2}+{\overline {u}}^{2}}{4}}\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e9d4645c97d58a076f2ec78d6d0ded81312b0ec9)
Enfin, pour la tangente de la latitude, on a l’expression
![{\displaystyle {\frac {z}{r}}={\overline {u}}\sin q-{\overline {s}}\cos q\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/393d45b676c60894a12e7a09446f3c3bc803f39f)
de sorte qu’en substituant aussi la valeur de
en
on aura
![{\displaystyle l={\overline {u}}\sin {\overline {p}}-{\overline {s}}\cos {\overline {p}}+({\overline {x}}{\overline {u}}-{\overline {y}}{\overline {s}})\cos 2{\overline {p}}+({\overline {x}}{\overline {s}}-{\overline {y}}{\overline {u}})\sin 2{\overline {p}}+{\overline {x}}{\overline {u}}-{\overline {y}}{\overline {s}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ecf3997ece85ec5c877af1b42d141cc653e23077)
Ainsi l’on trouvera par des différentiations partielles les valeurs des quantités d’où dépendent les corrections cherchées.
21. Si dans ces corrections on néglige d’abord les excentricités et les inclinaisons des orbites, on aura simplement par les formules précédentes
![{\displaystyle {\begin{alignedat}{2}{\frac {dq}{d\,{\overline {p}}}}=&1,\qquad &{\frac {dq}{d\,{\overline {x}}}}=&2\cos {\overline {p}},&{\frac {dq}{d\,{\overline {y}}}}=&-2\sin {\overline {p}},\\{\frac {dr}{d\,{\overline {r}}}}=&1,&{\frac {dr}{d\,{\overline {x}}}}=&{\overline {r}}\sin {\overline {p}},&{\frac {dr}{d\,{\overline {y}}}}=&{\overline {r}}\cos {\overline {p}},\\{\frac {dl}{d\,{\overline {s}}}}=&-\cos {\overline {p}},&&{\frac {dl}{d\,{\overline {u}}}}=\sin {\overline {p}},\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4d4b83bdf3d6e07410073230e2e836cf1217ef2f)
et toutes les autres différences partielles seront nulles.
Donc
1o La correction de la longitude deviendra
![{\displaystyle \varpi +2\xi \cos {\overline {p}}-2\psi \sin {\overline {p}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a7681d28ec64cb771f40d0120632813b403d905)
et, mettant pour
et
les valeurs trouvées dans le no 14, elle se réduira à ![{\displaystyle \varpi -2\Psi \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/79bd24edbfb5ec1a2072a76e63a4c97bc9949f79)
2o La correction de la latitude sera nulle, puisque les quantités
et
deviennent nulles dans le cas dont il s’agit (17) ;
3o La correction du rayon vecteur sera
![{\displaystyle \rho +{\overline {r}}\left(\xi \sin {\overline {p}}+\psi \cos {\overline {p}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bb3eee20c3af7cc48cee6c373be29e476f67409d)
laquelle par la substitution des valeurs de
et
se réduit à ![{\displaystyle \rho +{\overline {r}}\Xi .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2dc0a8552833d00ba6f69d474226a020e4526b9c)
Les valeurs de
ont été données dans les nos 10, 14 ; en les réunissant, on aura ces formules très-simples
Correction de la longitude
![{\displaystyle {\begin{aligned}&-\mathrm {T} '\left[{\frac {2}{n'}}{\overline {r}}^{2}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{1}}{d{\overline {r}}}}+\left(1+{\frac {3}{n'^{2}}}\right){\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{1}-\left(1+{\frac {2}{n'}}+{\frac {3}{n'^{2}}}\right){\frac {{\overline {r}}^{2}}{{\overline {r}}'^{2}}}\right]\\&\ \ \quad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \times {\frac {\sin \left({\overline {p}}-{\overline {p}}'\right)}{1-n'^{2}}}\\&-\mathrm {T} '\left[{\frac {2}{2n'}}{\overline {r}}^{2}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{2}}{d{\overline {r}}}}+\left(2+{\frac {3}{2n'^{2}}}\right){\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{2}\right]{\frac {\sin 2\left({\overline {p}}-{\overline {p}}'\right)}{1-4n'^{2}}}\\&-\mathrm {T} '\left[{\frac {2}{3n'}}{\overline {r}}^{2}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{3}}{d{\overline {r}}}}+\left(3+{\frac {3}{3n'^{2}}}\right){\overline {r}}\left[{\overline {r}},{\overline {r}}'\right]_{3}\right]{\frac {\sin 3\left({\overline {p}}-{\overline {p}}'\right)}{1-9n'^{2}}}\\&+\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/d98422d33f8acd71ff954dbe4b71fce0c9de524f)
Correction du rayon vecteur
![{\displaystyle {\begin{aligned}-\mathrm {T} '{\overline {r}}^{3}{\frac {1}{3}}{\frac {\left[{\overline {r}},{\overline {r}}'\right]}{\overline {r}}}&+\mathrm {T} '\left[{\overline {r}}^{3}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{1}}{d{\overline {r}}}}+{\frac {2}{n'}}{\overline {r}}^{2}\left[{\overline {r}},{\overline {r}}'\right]_{1}-\left(1+{\frac {2}{n'}}\right){\frac {{\overline {r}}^{3}}{{\overline {r}}'^{2}}}\right]\\&\ \qquad \qquad \qquad \qquad \qquad \qquad \qquad \times {\frac {\cos \left({\overline {p}}-{\overline {p}}'\right)}{1-n'^{2}}}\\&+\mathrm {T} '\left({\overline {r}}^{3}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{2}}{d{\overline {r}}}}+{\frac {2}{n'}}{\overline {r}}^{2}\left[{\overline {r}},{\overline {r}}'\right]_{2}\right){\frac {\cos 2\left({\overline {p}}-{\overline {p}}'\right)}{1-4n'^{2}}}\\&+\mathrm {T} '\left({\overline {r}}^{3}{\frac {d\left[{\overline {r}},{\overline {r}}'\right]_{3}}{d{\overline {r}}}}+{\frac {2}{n'}}{\overline {r}}^{2}\left[{\overline {r}},{\overline {r}}'\right]_{3}\right){\frac {\cos 3\left({\overline {p}}-{\overline {p}}'\right)}{1-9n'^{2}}}\\&+\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/19c215370cd50ce1943c286e078f1e8957c2e032)
Ces formules représentent les inégalités périodiques d’une Planète quelconque
troublée par une autre Planète
la masse du Soleil étant prise pour l’unité ; chaque Planète perturbatrice en donnera de pareilles, et leur somme exprimera l’effet total des perturbations, en tant qu’on fait abstraction des excentricités et des inclinaisons des orbites ; ce qui suffit dans la plupart des cas. Le calcul ne sera pas plus difficile, mais seulement un peu plus long, lorsqu’on voudra tenir cornpte aussi des termes dus aux excentricités et aux inclinaisons.
Quant à l’application de ces formules aux Planètes, elle n’aura d’autre difficulté que celle qui consiste dans l’évaluation des fonctions ![{\displaystyle \left[{\overline {r}},{\overline {r}}'\right],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c619ac470080600fd7bab415cc749873c37ded43)
et de leurs différences ; et nous avons aussi déjà donné tous les secours nécessaires pour ce calcul dans la Théorie des variations séculaires, à laquelle nous renvoyons.