RECHERCHES
SUR LA
THÉORIE DES PERTURBATIONS
QUE LES COMÈTES PEUVENT ÉPROUVER PAR L’ACTION DES PLANÈTES.
(Mémoires des Savants étrangers, t. X, 1785. Prix de l’Académie pour l’année 1778.)
Ces Recherches sont divisées en quatre Sections, que je vais parcourir sommairement.
Dans la première Section, je donne d’abord les équations générales du mouvement d’une comète autour du Soleil, en ayant égard aux perturbations qu’elle peut éprouver par l’action d’une ou de plusieurs planètes, et en rapportant les lieux tant de la comète que des planètes à des coordonnées rectangles. Je simplifie ensuite ces équations en partageant chacune d’elles en deux, dont l’une appartienne à l’orbite non altérée et dont l’autre renferme l’effet des perturbations, et je fais voir qu’en négligeant, à l’exemple des grands Géomètres qui ont déjà traité la Théorie des comètes, les carrés et les produits des forces perturbatrices, on peut considérer à part l’action de chaque planète, et prendre la somme des effets de leurs différentes actions pour l’effet total de leurs actions réunies. Enfin je montre comment on peut satisfaire aux équations différentielles des perturbations, dans le cas où la comète serait à une distance du Soleil infiniment grande par rapport à la distance de la planète au Soleil d’où résulte naturellement une transformation de ces mêmes équations, laquelle en facilite beaucoup l’intégration relativement à la partie supérieure de l’orbite de la comète. Cette transformation tient lieu des méthodes synthétiques proposées jusqu’ici pour simplifier le calcul des perturbations dans les régions supérieures de l’orbite, et elle a en même temps l’avantage de conserver l’uniformité dans la marche du calcul.
La deuxième Section est destinée uniquement à l’intégration des équations différentielles de l’orbite non altérée, et contient une solution complète du fameux Problème que Newton a résolu le premier, et une foule d’Auteurs après lui. Je me flatte que mon analyse pourra paraître encore digne de l’attention des Géomètres par sa simplicité et par sa généralité elle est d’ailleurs nécessaire pour les calculs de la Section suivante, et fournit différentes formules qui sont d’un grand usage dans tout le cours de cet Ouvrage.
Dans la troisième Section, je m’occupe de l’intégration des équations différentielles des perturbations. Je fais voir comment leurs intégrales se déduisent naturellement de celles des équations de l’orbite non altérée, en y faisant varier les constantes arbitraires qui représentent les éléments de l’orbite ce qui conduit directement à exprimer l’effet des perturbations par la variation des éléments de l’orbite considérée comme elliptique et ces variations se trouvent déterminées par des formules différentielles assez simples, dont chacune ne demande qu’une seule intégration. Je fais ensuite usage des transformations proposées dans la première Section pour les parties supérieures de l’orbite ; les formules différentielles dont il s’agit deviennent par là composées d’une partie absolument intégrable et d’une partie non intégrable, mais qui est toujours d’autant plus petite que la comète est plus éloignée du Soleil, en sorte qu’elle devient insensible lorsque la comète est à une très-grande distance du Soleil. Je termine cette Section par les formules générales qui expriment l’altération de la durée des révolutions anomalistiques et périodiques de la comète.
La quatrième Section contient l’application des méthodes et des formules données dans les Sections précédentes aux perturbations des comètes, et en particulier à celles de la comète de 1532 et de 1661. Toute la difficulté de cette application consiste dans l’intégration des formules différentielles qui déterminent les variations des éléments de l’orbite.
Après avoir mis ces formules sous une forme plus simple et plus commode pour le calcul, je montre les obstacles qui s’opposent à leur intégration générale, et qui obligent d’avoir recours aux quadratures des courbes mécaniques. Comme la méthode de ces quadratures est assez connue par les Ouvrages de Cotes et de Stirling, je n’entre là-dessus dans aucun détail ; mais je remarque qu’il y a des cas où l’usage de cette méthode cesse d’être légitime c’est lorsque la distance entre la comète et la planète perturbatrice est fort petite et approche de son minimum. Je donne pour ces sortes de cas une méthode particulière, qui réduit l’intégration aux logarithmes ou aux arcs de cercles, et ne peut jamais être sujette à aucun inconvénient. Tout ce que nous venons de dire ne regarde que la partie inférieure de l’orbite de la comète ; car, pour la partie supérieure de cette orbite, dans laquelle la distance de la comète au Soleil sera beaucoup plus grande que la distance de la planète au Soleil, je fais voir que la partie des formules différentielles qu’il reste à enregistrer se partage de nouveau en deux parties l’une indépendante du lieu de la planète, et qui est absolument intégrable ; l’autre qui contient les sinus ou cosinus de l’angle du moyen mouvement de la planète, et qui n’est intégrable par aucune méthode connue, mais dont je démontre que l’intégrale est nécessairement beaucoup plus petite que celle de la première partie, en sorte qu’on peut la négliger entièrement ; et, au cas qu’on voulût pousser l’exactitude plus loin, je donne un moyen d’approcher de plus en plus de la vraie valeur de cette intégrale. D’où il s’ensuit que, dans les régions supérieures de l’orbite des comètes, on peut déterminer leurs perturbations par des formules analytiques, qui ne demandent que des substitutions numériques pour donner les résultats cherchés, comme dans le cas des planètes. Je considère enfin la comète des années 1532 et 1661, que les Astronomes attendent vers 1789 ou 1790, et je déduis des éléments de cette comète toutes les données nécessaires pour le calcul de ses perturbations. Comme, dans le Programme de 1778, on n’exige pas que les concurrents donnent les résultats numériques de ce calcul, je m’abstiens d’entrer dans aucun détail à cet égard ; mais je me flatte qu’il n’y aura point de calculateur tant soit peu intelligent qui ne soit en état d’appliquer à la comète dont il s’agit la Théorie exposée dans cet Ouvrage.
Tels sont les principaux objets du travail que je soumets au jugement de l’Académie ; j’en serai suffisamment récompensé si cette illustre Compagnie daigne l’honorer de quelque attention.
section première.
équations différentielles du mouvement d’une comète autour du soleil, en ayant égard aux perturbations qu’elle peut éprouver par l’action des planètes.
1. Je prends la masse du Soleil pour l’unité, et je nomme
la masse de la comète,
les masses des planètes perturbatrices. Il est clair que ces quantités
doivent être des fractions très-petites, puisqu’elles expriment les rapports des masses de la comète et des planètes à la masse du Soleil ; en effet on sait que Jupiter, la plus grosse de toutes les planètes, a environ mille fois moins de masse que le Soleil ; et quant aux masses des comètes, quoiqu’elles soient inconnues, on ne peut guère les supposer plus grandes que celle de Jupiter, autrement il pourrait résulter de leur attraction des dérangements sensibles dans les orbites des planètes ; ce que les observations n’ont pas encore fait connaître, et ce qu’on ne suppose pas d’ailleurs qui arrive dans le Problème des comètes, tel qu’on l’a envisagé jusqu’à présent.
Nous regarderons donc dans la suite et nous traiterons les quantités
comme des quantités très-petites, dont il sera permis de négliger les puissances et les produits de deux ou de plusieurs dimensions. Cette supposition est conforme à ce que les Géomètres ont pratiqué jusqu’ici dans la Théorie des planètes principales et dans celle des comètes, et une plus grande exactitude ne serait peut-être d’aucune utilité.
2. Je rapporte les orbites que la comète et les planètes décrivent autour du Soleil à des coordonnées rectangles, prises du centre de cet astre, et parallèles à trois droites fixes et perpendiculaires entre elles.
Et je nomme ces coordonnées
pour l’orbite de la comète ;
pour l’orbite de la planète
pour l’orbite de la planète
etc.
Je nomme de plus
la distance de la comète au Soleil, ou le rayon vecteur de son orbite ;
le rayon vecteur de l’orbite de la planète
le rayon vecteur de l’orbite de la planète
Enfin je désigne par
la distance de la planète
à la comète ; par
la distance de la planète
à la comète, etc.
Il est clair qu’on aura
![{\displaystyle r={\sqrt {x^{2}+y^{2}+z^{2}}},\quad \rho ={\sqrt {\xi ^{2}+\eta ^{2}+\zeta ^{2}}},\quad \rho '={\sqrt {\xi '^{2}+\eta '^{2}+\zeta '^{2}}},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc4db582dffd68f7e454a2faac64bd46aea49683)
![{\displaystyle \mathrm {R} ={\sqrt {(x-\xi )^{2}+(y-\eta )^{2}+(z-\zeta )^{2}}},\quad \mathrm {R} '={\sqrt {(x-\xi ')^{2}+(y-\eta ')^{2}+(z-\zeta ')^{2}}},\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fbbaa531898810a234c71bc0688d9b4937b6d18c)
3. Cela posé, si l’on décompose, suivant les directions des trois coordonnées rectangles
toutes les forces qui agissent sur la comète pour lui faire décrire son orbite autour du Soleil, savoir les attractions
exercées par le Soleil et par les planètes sur la comète, et les attractions
exercées par la comète et par les planètes sur le Soleil, et qui doivent être transportées à la comète en sens contraire et qu’on égale la somme de toutes les forces qui agissent suivant la ligne
et qui tendent à diminuer cette ligne à
la somme de toutes les forces qui agissent suivant
à
et la somme de toutes les forces qui agissent suivant
à
étant les éléments du temps supposés constants, on aura ces trois équations
![{\displaystyle {\begin{aligned}&{\frac {d^{2}x}{dt^{2}}}+{\frac {(1+m)x}{r^{3}}}+\mu \left({\frac {x-\xi }{\mathrm {R} ^{3}}}+{\frac {\xi }{\rho ^{3}}}\right)+\mu '\left({\frac {x-\xi '}{\mathrm {R} '^{3}}}+{\frac {\xi '}{\rho '^{3}}}\right)+\ldots =0,\\&{\frac {d^{2}y}{dt^{2}}}\,+{\frac {(1+m)y}{r^{3}}}+\mu \left({\frac {y-\eta }{\mathrm {R} ^{3}}}+{\frac {\eta }{\rho ^{3}}}\right)+\mu '\left({\frac {y-\eta '}{\mathrm {R} '^{3}}}+{\frac {\eta '}{\rho '^{3}}}\right)+\ldots =0,\\&{\frac {d^{2}z}{dt^{2}}}\,+{\frac {(1+m)z}{r^{3}}}+\mu \left({\frac {z-\zeta }{\mathrm {R} ^{3}}}+{\frac {\zeta }{\rho ^{3}}}\right)\,+\mu '\left({\frac {z-\zeta '}{\mathrm {R} '^{3}}}\,+{\frac {\zeta '}{\rho '^{3}}}\right)+\ldots =0,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe834da5fa0a575199758a7ccb49ca8bb8d970a9)
lesquelles, d’après les principes connus de la Dynamique, serviront à déterminer le mouvement de la comète par rapport au Soleil regardé comme immobile.
4. On aura des équations semblables pour le mouvement de la planète
autour du Soleil, en tant qu’elle est dérangée par l’action de la comète et des autres planètes ; pour cela, il n’y aura qu’à changer, dans les équations précédentes, les quantités
appartenant à l’orbite de la comète, dans les quantités analogues
appartenant à l’orbite de la planète, et vice versâ, celles-ci en celles-là.
Mais il faut considérer que pour notre objet il n’est pas nécessaire de tenir compte des termes affectés des quantités très-petites
dans les équations de la planète, parce que les quantités
dépendant de ces équations ne se trouvent dans les équations de la comète que dans des termes déjà affectés de la quantité très-petite
On peut donc réduire les équations de la planète
aux deux premiers termes, savoir
![{\displaystyle {\begin{aligned}&{\frac {d^{2}\xi }{dt^{2}}}+{\frac {(1+\mu )\xi }{\rho ^{3}}}=0,\\&{\frac {d^{2}\eta }{dt^{2}}}+{\frac {(1+\mu )\eta }{\rho ^{3}}}=0,\\&{\frac {d^{2}\zeta }{dt^{2}}}+{\frac {(1+\mu )\zeta }{\rho ^{3}}}=0\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0713199937f42d4a55014e367b0573815ebc2f59)
et l’on réduira, par des raisons semblables, les équations du mouvement
de la planète
![{\displaystyle \mu '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53b9ddd6bbc6f1ef9ca40b1c2e6e2c1c8d141aed)
à celles-ci
![{\displaystyle {\begin{aligned}&{\frac {d^{2}\xi '}{dt^{2}}}+{\frac {(1+\mu ')\xi '}{\rho '^{3}}}=0,\\&{\frac {d^{2}\eta '}{dt^{2}}}+{\frac {(1+\mu ')\eta '}{\rho '^{3}}}=0,\\&{\frac {d^{2}\zeta '}{dt^{2}}}+{\frac {(1+\mu ')\zeta '}{\rho '^{3}}}=0\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5397fcbf449d99dd0da3ac930223f32eff950ea5)
et ainsi pour les autres planètes perturbatrices.
Ces réductions sont fondées, comme l’on voit, sur la supposition que, dans le calcul des perturbations des comètes, on néglige les perturbations des planètes perturbatrices. Si cette supposition n’est pas rigoureusement exacte, elle est du moins permise dans la première approximation, à laquelle nous nous contenterons ici de borner nos Recherches, à l’exemple des grands Géomètres qui ont traité avant nous le Problème des comètes.
5. En considérant les expressions des quantités
il est aisé de voir qu’on peut mettre les équations précédentes sous une forme plus simple que voici :
Pour la comète,
![{\displaystyle {\begin{aligned}&{\frac {d^{2}x}{dt^{2}}}=(1+m){\frac {d{\cfrac {1}{r}}}{dx}}+\mu {\frac {d\left({\cfrac {1}{\rho }}-{\cfrac {1}{\mathrm {R} }}\right)}{d\xi }}+\mu '{\frac {d\left({\cfrac {1}{\rho '}}-{\cfrac {1}{\mathrm {R} '}}\right)}{d\xi '}}+\ldots ,\\&{\frac {d^{2}y}{dt^{2}}}=(1+m){\frac {d{\cfrac {1}{r}}}{dy}}+\mu {\frac {d\left({\cfrac {1}{\rho }}-{\cfrac {1}{\mathrm {R} }}\right)}{d\eta }}+\mu '{\frac {d\left({\cfrac {1}{\rho '}}-{\cfrac {1}{\mathrm {R} '}}\right)}{d\eta '}}+\ldots ,\\&{\frac {d^{2}z}{dt^{2}}}=(1+m){\frac {d{\cfrac {1}{r}}}{dz}}+\mu {\frac {d\left({\cfrac {1}{\rho }}-{\cfrac {1}{\mathrm {R} }}\right)}{d\zeta }}+\mu '{\frac {d\left({\cfrac {1}{\rho '}}-{\cfrac {1}{\mathrm {R} '}}\right)}{d\zeta '}}+\ldots \,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/72c9f2d2b72f90d420ac0b705e0b9101c45c97bc)
Pour la planète ![{\displaystyle \mu ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1e7e1ef161a49a22b500d63307460ad92eeb6a16)
![{\displaystyle {\frac {d^{2}\xi }{dt^{2}}}=(1+\mu ){\frac {d{\cfrac {1}{\rho }}}{d\xi }},\quad {\frac {d^{2}\eta }{dt^{2}}}=(1+\mu ){\frac {d{\cfrac {1}{\rho }}}{d\eta }},\quad {\frac {d^{2}\zeta }{dt^{2}}}=(1+\mu ){\frac {d{\cfrac {1}{\rho }}}{d\zeta }}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f14184b9c42190617b5300526207267d1ed9e4d5)
Pour la planète ![{\displaystyle \mu ',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/edb713e4931434716e2aae1c44e2b80607293b34)
![{\displaystyle {\frac {d^{2}\xi '}{dt^{2}}}=(1+\mu '){\frac {d{\cfrac {1}{\rho '}}}{d\xi '}},\quad {\frac {d^{2}\eta '}{dt^{2}}}=(1+\mu '){\frac {d{\cfrac {1}{\rho '}}}{d\eta '}},\quad {\frac {d^{2}\zeta '}{dt^{2}}}=(1+\mu '){\frac {d{\cfrac {1}{\rho '}}}{d\zeta '}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/149ba333933c2a126afbafa252f8d0ab56213679)
et ainsi des autres.
Dans ces formules, les expressions
dénotent, suivant la notation reçue parmi les Géomètres, les coefficients de
dans la différentielle de
et ainsi des autres expressions semblables.
6. Si l’on suppose que dans le mouvement de la comète on fasse abstraction des forces perturbatrices, il faudra rejeter, dans les équations de la comète, les termes affectés de
on aura ainsi
Pour l’orbite non altérée de la comète,
![{\displaystyle {\frac {d^{2}x}{dt^{2}}}=(1+m){\frac {d{\cfrac {1}{r}}}{dx}},\quad {\frac {d^{2}y}{dt^{2}}}=(1+m){\frac {d{\cfrac {1}{r}}}{dy}},\quad {\frac {d^{2}z}{dt^{2}}}=(1+m){\frac {d{\cfrac {1}{r}}}{dz}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11fff5133a3cde4e2398dd8d3b09c0b87a5c6217)
Nous pouvons supposer que les quantités
se rapportent à l’orbite non altérée, et sont par conséquent déterminées par les équations précédentes ; dans cette supposition, il est clair que les vraies valeurs des quantités
dans l’orbite troublée, ne peuvent différer des précédentes que par des quantités très-petites de l’ordre de
qu’on peut désigner, pour plus de simplicité, par la caractéristique
à la manière des différences ordinaires ; et la recherche des perturbations de la comète se réduira à déterminer les valeurs des différences
7. Dorénavant donc les quantités
appartiendront toujours à l’orbite non altérée de la comète, et devront par conséquent se déterminer par les équations du numéro précédent. Dans l’orbite troublée, ces quantités deviendront
et devront être déterminées par les équations du no 5, en mettant dans ces équations ces nouvelles quantités à la place des premières
Or, comme les différences
sont très-petites de l’ordre
il suffira de conserver, dans cette substitution, les premières dimensions de ces différences (par l’hypothèse du no 1), dans les termes non affectés de
et, dans les termes affectés de ces quantités, on pourra négliger tout à fait les différences dont il s’agit.
Ainsi donc le terme
de la première équation du no 5 deviendra
![{\displaystyle {\frac {d^{2}x}{dt^{2}}}+{\frac {d^{2}\delta x}{dt^{2}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d0fea87c84f5b16a5a564889e7bd42d6005585e1)
le terme
de la même équation deviendra
![{\displaystyle (1+m){\frac {d{\cfrac {1}{r}}}{dx}}+(1+m)\left({\frac {d^{2}{\cfrac {1}{r}}}{dx^{2}}}\delta x+{\frac {d^{2}{\cfrac {1}{r}}}{dxdy}}\delta y+{\frac {d^{2}{\cfrac {1}{r}}}{dxdz}}\delta z\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eb8ddbcce9b0e089c8f70873d841f7e058227540)
et les autres termes demeureront les mêmes.
On transformera de même la deuxième et la troisième équation du mouvement de la comète, et, effaçant ensuite les termes qui se détruisent en vertu des équations du no 6, on aura ces trois-ci, qui serviront à déterminer les valeurs des quantités
, dues aux perturbations de la comète :
Pour les perturbations de la comète,
![{\displaystyle {\begin{aligned}{\frac {d^{2}\delta x}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\cfrac {1}{r}}}{dx^{2}}}\delta x+{\frac {d^{2}{\cfrac {1}{r}}}{dxdy}}\delta y+{\frac {d^{2}{\cfrac {1}{r}}}{dxdz}}\delta z\right)\\&\qquad \qquad \qquad +\mu {\frac {d\left({\cfrac {1}{\rho }}-{\cfrac {1}{\mathrm {R} }}\right)}{d\xi }}+\mu '{\frac {d\left({\cfrac {1}{\rho '}}-{\cfrac {1}{\mathrm {R} '}}\right)}{d\xi '}}+\ldots ,\\{\frac {d^{2}\delta y}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\cfrac {1}{r}}}{dydx}}\delta x+{\frac {d^{2}{\cfrac {1}{r}}}{dy^{2}}}\delta y+{\frac {d^{2}{\cfrac {1}{r}}}{dydz}}\delta z\right)\\&\qquad \qquad \qquad +\mu {\frac {d\left({\cfrac {1}{\rho }}-{\cfrac {1}{\mathrm {R} }}\right)}{d\eta }}+\mu '{\frac {d\left({\cfrac {1}{\rho '}}-{\cfrac {1}{\mathrm {R} '}}\right)}{d\eta '}}+\ldots ,\\{\frac {d^{2}\delta z}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\cfrac {1}{r}}}{dzdx}}\delta x+{\frac {d^{2}{\cfrac {1}{r}}}{dzdy}}\delta y+{\frac {d^{2}{\cfrac {1}{r}}}{dz^{2}}}\delta z\right)\\&\qquad \qquad \qquad +\mu {\frac {d\left({\cfrac {1}{\rho }}-{\cfrac {1}{\mathrm {R} }}\right)}{d\zeta }}+\mu '{\frac {d\left({\cfrac {1}{\rho '}}-{\cfrac {1}{\mathrm {R} '}}\right)}{d\zeta '}}+\ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bee1648f597a6ff220421f7e90476cac970fe3eb)
C’est donc de l’intégration de ces équations que dépend la solution du Problème des perturbations des comètes. Nous allons nous en occuper, après avoir fait quelques remarques générales sur la nature de ces équations.
8. J’observe d’abord que, comme ces équations ne renferment que les premières dimensions des variables
on peut chercher à part les valeurs de ces variables pour les différents termes affectés des quantités
et qui viennent de l’action des différentes planètes dont ces quantités sont les masses ; car il est visible que, si l’on réunit ensuite ces différentes valeurs, on aura les valeurs complètes des variables
qui satisfont aux équations proposées.
En général, il est facile de concevoir que, lorsqu’on néglige, ainsi que nous l’avons fait, les carrés et les produits des forces perturbatrices, l’effet total de ces forces doit être égal à la somme des effets que chacune en particulier produirait si elle était seule.
9. Je remarque ensuite que les termes multipliés par les masses
des planètes perturbatrices deviennent d’autant plus petits que les quantités
sont plus petites, c’est-à-dire, que la comète est plus près du Soleil. En effet, en supposant
des quantités très-petites vis-à-vis de
on a à très-peu près (no 2)
![{\displaystyle {\frac {1}{\mathrm {R} }}={\frac {1}{\rho }}-{\frac {d{\cfrac {1}{\rho }}}{d\xi }}x-{\frac {d{\cfrac {1}{\rho }}}{d\eta }}y-{\frac {d{\cfrac {1}{\rho }}}{d\zeta }}z\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc28d0469f53f8230fe27cb4034f8e7956ad8a27)
d’où l’on voit que la quantité
ainsi que ses différences divisées par
seront du même ordre de petitesse que les quantités
Par conséquent, si l’on suppose que ces quantités soient devenues de l’ordre des quantités
il est clair que les termes dont il s’agit seront pour lors de l’ordre de
de sorte qu’on pourra les négliger, d’autant plus que, dans ce cas, la quantité
devient d’autant
plus grande. Ces termes disparaissant, il est visible qu’on pourra satisfaire aux équations proposées par la supposition de
![{\displaystyle \delta z=0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/130ce912071a8a376514e2ba2b89bc67f4922ed2)
Ainsi l’on peut regarder ces valeurs comme les limites des variables
![{\displaystyle \delta x,\delta y,\delta z,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/91994651e02e93f207cf9de0e9cdad63a668c60d)
du côté du Soleil.
10. Voyons maintenant quelles sont les limites des mêmes variables du côté opposé.
Supposons donc les quantités
infiniment grandes vis-à-vis de
on aura ici (no 2)
![{\displaystyle {\frac {1}{\mathrm {R} }}={\frac {1}{r}}-{\frac {d{\frac {1}{r}}}{dx}}\xi -{\frac {d{\frac {1}{r}}}{dy}}\eta -{\frac {d{\frac {1}{r}}}{dz}}\zeta +{\frac {1}{2}}{\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\xi ^{2}+{\frac {1}{2}}{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\eta ^{2}+{\frac {1}{2}}{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\zeta ^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5560813349362a041e614e109d23424fbf23f457)
![{\displaystyle +{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\xi \eta +{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\xi \zeta +{\frac {d^{2}{\frac {1}{r}}}{dydz}}\eta \zeta +\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb7d50063808bc3cb5000fd3d2717fa70a8715aa)
J’ai poussé ici l’approximation jusqu’à la seconde dimension des quantités
parce que la différentiation par
fait disparaître une dimension de ces quantités.
On aura donc
![{\displaystyle {\begin{aligned}{\frac {d{\frac {1}{\mathrm {R} }}}{d\xi }}=-{\frac {d{\frac {1}{r}}}{dx}}+{\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\xi +{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\eta +{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\zeta ,\\{\frac {d{\frac {1}{\mathrm {R} }}}{d\eta }}=-{\frac {d{\frac {1}{r}}}{dy}}+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\xi +{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\eta +{\frac {d^{2}{\frac {1}{r}}}{dydz}}\zeta ,\\{\frac {d{\frac {1}{\mathrm {R} }}}{d\zeta }}=-{\frac {d{\frac {1}{r}}}{dz}}+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\xi +{\frac {d^{2}{\frac {1}{r}}}{dydz}}\eta +{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\zeta .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/961646c94aa7dc0435ea1312147791894152dacd)
Qu’on substitue ces valeurs dans les équations du no 7, en n’ayant égard qu’aux termes affectés de
par la remarque ci-dessus (no 8), on aura les équations suivantes
![{\displaystyle {\frac {d^{2}\delta x}{dt^{2}}}=(1+m)\left[{\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\left(\delta x-{\frac {\mu }{1+m}}\xi \right)+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\left(\delta y-{\frac {\mu }{1+m}}\eta \right)\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eeef94cdcf737d2f3db9d422d7b2bb92bef1b20a)
![{\displaystyle \left.+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\left(\delta z-{\frac {\mu }{1+m}}\zeta \right)\right]+\mu \left({\frac {d{\frac {1}{\rho }}}{d\xi }}+{\frac {d{\frac {1}{r}}}{dx}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/33bbf0d7427dde0495050f0f30b8eccd890bc512)
![{\displaystyle {\frac {d^{2}\delta y}{dt^{2}}}=(1+m)\left[{\frac {d^{2}{\frac {1}{r}}}{dydx}}\left(\delta x-{\frac {\mu }{1+m}}\xi \right)+{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\left(\delta y-{\frac {\mu }{1+m}}\eta \right)\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f37b38b851fb848c7679b5c4b7ab58846d2b29b)
![{\displaystyle \left.+{\frac {d^{2}{\frac {1}{r}}}{dydz}}\left(\delta z-{\frac {\mu }{1+m}}\zeta \right)\right]+\mu \left({\frac {d{\frac {1}{\rho }}}{d\eta }}+{\frac {d{\frac {1}{r}}}{dy}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/59de7d12cc7e187b434a307ba9e8946d9903e2cd)
![{\displaystyle {\frac {d^{2}\delta z}{dt^{2}}}=(1+m)\left[{\frac {d^{2}{\frac {1}{r}}}{dzdx}}\left(\delta x-{\frac {\mu }{1+m}}\xi \right)+{\frac {d^{2}{\frac {1}{r}}}{dzdy}}\left(\delta y-{\frac {\mu }{1+m}}\eta \right)\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b13df67f68bd2e3ec5befb0288e85e6698ff783)
![{\displaystyle \left.+{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\left(\delta z-{\frac {\mu }{1+m}}\zeta \right)\right]+\mu \left({\frac {d{\frac {1}{\rho }}}{d\zeta }}+{\frac {d{\frac {1}{r}}}{dz}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b08cb29cde293bdcbaa0dfd792daad40ba1e861f)
Or on a, par les équations de la planète
(no 5),
![{\displaystyle {\frac {d{\frac {1}{\rho }}}{d\xi }}={\frac {1}{1+\mu }}{\frac {d^{2}\xi }{dt^{2}}},\quad {\frac {d{\frac {1}{\rho }}}{d\eta }}={\frac {1}{1+\mu }}{\frac {d^{2}\eta }{dt^{2}}},\quad {\frac {d{\frac {1}{\rho }}}{d\zeta }}={\frac {1}{1+\mu }}{\frac {d^{2}\zeta }{dt^{2}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/725e1b8e31f78df2bba4728803c5247b41357357)
faisant ces substitutions dans les pénultièmes termes des équations précédentes, on verra incontinent que, si les derniers termes
![{\displaystyle \mu {\frac {d{\frac {1}{r}}}{dx}},\quad \mu {\frac {d{\frac {1}{r}}}{dy}},\quad \mu {\frac {d{\frac {1}{r}}}{dz}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/331d34c5f10b570fb4f250065d59b018d1017e43)
n’existaient pas, et que l’on eût
on satisferait exactement à ces équations, en faisant
![{\displaystyle \delta x={\frac {\mu }{1+\mu }}\xi ,\quad \delta y={\frac {\mu }{1+\mu }}\eta ,\quad \delta z={\frac {\mu }{1+\mu }}\zeta .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/84a34f15b2e946ff59dd075c4596acdc9249e8ab)
Supposons donc
![{\displaystyle {\begin{aligned}\delta x=&{\frac {\mu }{1+\mu }}\xi +\alpha ,\\\delta y=&{\frac {\mu }{1+\mu }}\eta +\beta ,\\\delta z=&{\frac {\mu }{1+\mu }}\zeta +\gamma \,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/013bda6ebda7fc6aa25546ba03fc8d7ca33a15a6)
si l’on substitue ces valeurs, et qu’on rejette les termes qui auraient pour coefficient
d’après l’hypothèse éta-
blie dans le
no 1, on aura ces nouvelles équations
![{\displaystyle {\begin{aligned}{\frac {d^{2}\alpha }{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\alpha +{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\beta +{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\gamma \right)+\mu {\frac {d{\frac {1}{r}}}{dx}},\\{\frac {d^{2}\beta }{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dydx}}\alpha +{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\beta +{\frac {d^{2}{\frac {1}{r}}}{dydz}}\gamma \right)+\mu {\frac {d{\frac {1}{r}}}{dy}},\\{\frac {d^{2}\gamma }{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dzdx}}\alpha +{\frac {d^{2}{\frac {1}{r}}}{dzdy}}\beta +{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\gamma \right)+\mu {\frac {d{\frac {1}{r}}}{dz}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/17e86ebfdd547b0d83292cc608f91ebc72beb4ca)
J’observe qu’on peut satisfaire à ces équations en faisant
étant un coefficient constant ; car elles deviennent par là
![{\displaystyle {\begin{aligned}\mathrm {K} {\frac {d^{2}x}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}x+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}y+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}z\right)\mathrm {K} +\mu {\frac {d{\frac {1}{r}}}{dx}},\\\mathrm {K} {\frac {d^{2}y}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dydx}}x+{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}y+{\frac {d^{2}{\frac {1}{r}}}{dydz}}z\right)\mathrm {K} +\mu {\frac {d{\frac {1}{r}}}{dy}},\\\mathrm {K} {\frac {d^{2}z}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dzdx}}x+{\frac {d^{2}{\frac {1}{r}}}{dzdy}}y+{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}z\right)\mathrm {K} +\mu {\frac {d{\frac {1}{r}}}{dz}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8c1c5554738e39d3df54326d6ee2dc239d978cbf)
Mais les équations de l’orbite non altérée (no 6) donnent
![{\displaystyle {\frac {d^{2}x}{dt^{2}}}=(1+m){\frac {d{\frac {1}{r}}}{dx}},\quad {\frac {d^{2}y}{dt^{2}}}=(1+m){\frac {d{\frac {1}{r}}}{dy}},\quad {\frac {d^{2}z}{dt^{2}}}=(1+m){\frac {d{\frac {1}{r}}}{dz}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/531f853c5c7b6a03f50cf31ce479870c6bb538e4)
Ensuite je remarque que la quantité
est une fonction homogène de
de la dimension
qu’ainsi les quantités
sont aussi
des fonctions homogènes de
![{\displaystyle x,y,z,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d08d690d7e19ea7aee8574fc6abd6a15d97fa026)
mais de la dimension
![{\displaystyle -2.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/535176f9d062cd724aa815b4e014dd014c6537a7)
Donc, par le Théorème connu concernant ces sortes de fonctions, on aura
![{\displaystyle {\begin{aligned}{\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}x+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}y+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}z=-2{\frac {d{\frac {1}{r}}}{dx}},\\{\frac {d^{2}{\frac {1}{r}}}{dydx}}x+{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}y+{\frac {d^{2}{\frac {1}{r}}}{dydz}}z=-2{\frac {d{\frac {1}{r}}}{dy}},\\{\frac {d^{2}{\frac {1}{r}}}{dzdx}}x+{\frac {d^{2}{\frac {1}{r}}}{dzdy}}y+{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}z=-2{\frac {d{\frac {1}{r}}}{dz}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7dcaf33ce300d6e7a2119d42b0034e8c338672f9)
C’est de quoi on peut d’ailleurs s’assurer par la différentiation actuelle. Ces substitutions faites, on verra d’abord que, pour satisfaire aux trois équations dont il s’agit, il suffit de satisfaire à celle-ci
![{\displaystyle \mathrm {K} (1+m)=-2(1+m)\mathrm {K} +\mu ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9e25a5da864e89b7a6b5383fd42759932bf890d3)
laquelle donne
![{\displaystyle \mathrm {K} ={\frac {\mu }{3(1+m)}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/293ab7b24a7981feab3e862132bb810a3ee1a9d1)
Donc enfin on aura
![{\displaystyle {\begin{aligned}\delta x=&{\frac {\mu }{1+\mu }}\xi +{\frac {\mu }{3(1+m)}}x,\\\delta y=&{\frac {\mu }{1+\mu }}\eta +{\frac {\mu }{3(1+m)}}y,\\\delta z=&{\frac {\mu }{1+\mu }}\zeta +{\frac {\mu }{3(1+m)}}z.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/28af69151a230779e3de61c5c7dee45077269c47)
Ce sont les limites cherchées, dont les quantités
s’approchent d’autant plus que la comète s’éloigne davantage du Soleil.
11. De là il s’ensuit que si l’on suppose, en général,
![{\displaystyle {\begin{aligned}\delta x=&\mu \left({\frac {x}{3(1+m)}}+{\frac {\xi }{1+\mu }}\right)+\delta x',\\\delta y=&\mu \left({\frac {y}{3(1+m)}}+{\frac {\eta }{1+\mu }}\right)+\delta y',\\\delta z=&\mu \left({\frac {z}{3(1+m)}}+{\frac {\zeta }{1+\mu }}\right)+\delta z',\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19cb77857bfffc1a1d2d17d8e48529a67f558376)
qu’on substitue ces valeurs dans les équations du no 7 et qu’on y fasse les réductions enseignées ci-dessus, on aura, en n’ayant égard qu’aux termes affectés de
et faisant, pour abréger,
![{\displaystyle {\frac {1}{\mathrm {S} }}={\frac {1}{r}}-{\frac {1+m}{1+\mu }}\left({\frac {d{\cfrac {1}{r}}}{dx}}\xi +{\frac {d{\cfrac {1}{r}}}{dy}}\eta +{\frac {d{\cfrac {1}{r}}}{dz}}\zeta \right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5d044c6d6a83be49d935eaff41c407c5d714f0be)
on aura, dis-je, ces transformées
![{\displaystyle {\begin{aligned}{\frac {d^{2}\delta x'}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\delta x'+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\delta y'+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\delta z'\right)-\mu {\frac {d\mathrm {\left({\frac {1}{S}}-{\frac {1}{R}}\right)} }{dx}},\\{\frac {d^{2}\delta y'}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dydx}}\delta x'+{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\delta y'+{\frac {d^{2}{\frac {1}{r}}}{dydz}}\delta z'\right)-\mu {\frac {d\mathrm {\left({\frac {1}{S}}-{\frac {1}{R}}\right)} }{dy}},\\{\frac {d^{2}\delta z'}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dzdx}}\delta x'+{\frac {d^{2}{\frac {1}{r}}}{dzdy}}\delta y'+{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\delta z'\right)-\mu {\frac {d\mathrm {\left({\frac {1}{S}}-{\frac {1}{R}}\right)} }{dz}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e2a762e179534138a02257ebf20d642fb35fde58)
Dans ces équations, j’ai mis à la place des quantités
leurs équivalentes
pour mettre plus d’uniformité dans les formules.
12. On peut aussi donner une forme semblable aux équations primitives du no 7. En effet, si l’on fait
![{\displaystyle {\frac {1}{\sigma }}={\frac {1}{\rho }}-\left({\frac {d{\cfrac {1}{\rho }}}{d\xi }}x+{\frac {d{\cfrac {1}{\rho }}}{d\eta }}y+{\frac {d{\cfrac {1}{\rho }}}{d\zeta }}z\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b6c93a19500657573426fe4589872c853289b2e8)
et que l’on fasse abstraction des termes affectés de
on aura
![{\displaystyle {\begin{aligned}{\frac {d^{2}\delta x}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\delta z\right)-\mu {\frac {d\left({\frac {1}{\sigma }}-{\frac {1}{\mathrm {R} }}\right)}{dx}},\\{\frac {d^{2}\delta y}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dydx}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dydz}}\delta z\right)-\mu {\frac {d\left({\frac {1}{\sigma }}-{\frac {1}{\mathrm {R} }}\right)}{dy}},\\{\frac {d^{2}\delta z}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dzdx}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dzdy}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\delta z\right)-\mu {\frac {d\left({\frac {1}{\sigma }}-{\frac {1}{\mathrm {R} }}\right)}{dz}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/539feaa8dfc41e9ce9388b2d298f9ac053ed2a67)
13. On voit que la quantité
contient les deux premiers termes de la quantité
développée en suite ascendante par rapport aux puissances de
comme la quantité
contient les deux premiers termes de la même quantité
développée en suite ascendante par rapport aux puissances de
en négligeant (ce qui est permis ici) la différence infiniment petite entre
et l’unité : d’où résultent naturellement les conclusions que nous avons trouvées plus haut (nos 10, 11).
Il s’ensuit aussi de là que, tant que
il est plus simple de se servir des formules du no 12, et qu’au contraire, lorsque
il est plus avantageux d’employer celles du no 11 ; d’autant plus que dans celles-ci les termes affectés de
et qui sont l’effet des forces perturbatrices, deviennent presque nuls lorsque la comète est à une grande distance du Soleil.
section deuxième.
intégration des équations différentielles de l’orbite non altérée.
14. Ayant décomposé les équations générales du mouvement de la comète en équations de l’orbite non troublée (no 6) et en équations des perturbations (no 7), nous allons nous occuper, dans cette Section, de l’intégration des premières. Nous pourrions, à la vérité, nous en dispenser, puisqu’on sait d’avance, par les Théorèmes de Newton, que, sans les forces perturbatrices, la comète doit décrire autour du Soleil une section conique dont cet astre occupe le foyer, et que le temps doit être proportionnel à l’aire parcourue, divisée par la racine carrée du paramètre. Mais, comme nous avons besoin de connaître les intégrales mêmes des équations dont il s’agit, il est beaucoup plus court et en même temps plus direct de chercher ces intégrales par l’intégration effective que de les déduire des propriétés des sections coniques.
Les équations qu’il s’agit d’intégrer sont celles-ci, en mettant pour
leurs valeurs
![{\displaystyle {\begin{aligned}&{\frac {d^{2}x}{dt^{2}}}+{\frac {(1+m)x}{r^{3}}}=0,\\&{\frac {d^{2}y}{dt^{2}}}+{\frac {(1+m)y}{r^{3}}}=0,\\&{\frac {d^{2}z}{dt^{2}}}+{\frac {(1+m)z}{r^{3}}}=0.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2356c01cc177ea4ebdcdfe0ed8c091e92fd6fbc1)
On peut intégrer ces équations par différentes méthodes ; celle dont je vais faire usage m’a paru une des plus simples.
Je remarque d’abord que, en supposant les deux premières équations, on peut satisfaire à la troisième en faisant
(A)
|
|
|
![{\displaystyle b}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f11423fbb2e967f986e36804a8ae4271734917c3)
et
![{\displaystyle c}](https://wikimedia.org/api/rest_v1/media/math/render/svg/86a67b81c2de995bd608d5b2df50cd8cd7d92455)
étant deux constantes arbitraires ; et il est visible que cette valeur de
![{\displaystyle z}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf368e72c009decd9b6686ee84a375632e11de98)
est en même temps l’intégrale complète de la troisième équation, puisqu’elle renferme deux constantes arbitraires.
Je multiplie maintenant la première des trois équations différentielles proposées par
la seconde par
la troisième par
ensuite je les ajoute ensemble, et j’intègre j’ai
(B)
|
|
|
étant une constante arbitraire.
Je multiplie ensuite les mêmes équations par
et j’ajoute à leur somme l’équation précédente j’ai, à cause de
(C)
|
|
|
Cette équation étant multipliée par
et ensuite intégrée, donne celle-ci
(D)
|
|
|
étant une nouvelle constante arbitraire. Or
![{\displaystyle {\frac {1}{2}}d^{2}\left(r^{2}\right)=rd^{2}r+dr^{2}\quad {\text{et}}\quad \left[{\frac {1}{2}}d\left(r^{2}\right)\right]^{2}=r^{2}dr^{2}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/880fd4307abe6001a0d49cba31e3823d693212ce)
donc, si l’on divise l’équation (D) par
et qu’on la retranche de l’équation (C), on aura, après avoir divisé par ![{\displaystyle r,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/250644a0f511e9078be6f89ba78a606a0e08c0a0)
![{\displaystyle {\frac {d^{2}r}{dt^{2}}}+(1+m)\left({\frac {1}{r^{2}}}-{\frac {2h}{r^{3}}}\right)=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f3aa4c81bb2c306cda7752a1ee0d5548c4dfcf15)
équation qui, en faisant
prend cette forme
![{\displaystyle {\frac {d^{2}p}{dt^{2}}}+(1+m){\frac {p}{r^{3}}}=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7dffd078a13cb4d39295e09e76657831d0d8e0ed)
qui est, comme l’on voit, semblable aux équations primitives.
C’est pourquoi on aura sur-le-champ cette intégrale
ou bien
(E)
|
|
|
et
étant deux nouvelles constantes arbitraires, en sorte que l’intégrale est complète.
Les équations (A) et (E) offrent déjà, comme l’on voit, deux intégrales finies. On trouvera la troisième au moyen de l’équation (D), laquelle se réduit à
![{\displaystyle {\frac {rdr}{\sqrt {r-{\cfrac {r^{2}}{a}}-h}}}=dt{\sqrt {2(1+m)}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1386c1b0e1b015cb782e0fda36aff53f8039c884)
dont l’intégrale est
(F)
|
|
|
étant encore une constante arbitraire.
Cette équation détermine
en
et les équations (A) et (E), combinées avec celle-ci
servent à déterminer
en
ainsi l’on aura
en
Mais ces valeurs, pour être complètes, doivent renfermer six constantes, parce que les équations différentielles proposées sont chacune du second ordre. Or l’équation (A) renferme deux constantes arbitraires
et
l’équation (E) en renferme trois
et
et l’équation (F) en renferme encore deux autres
et
Il y en a donc en tout sept, et par conséquent une de plus qu’il ne faut.
En examinant la chose de plus près, il est aisé de s’apercevoir que cela vient de ce que la constante
a été introduite par l’intégration qui a donné l’équation (B), équation dont nous n’avons point tenu compte dans la suite du calcul comme d’une équation intégrale. Il est donc nécessaire d’avoir égard à cette équation, et il en doit résulter une équation de condition entre les constantes ; en sorte qu’il n’en restera plus que six d’arbitraires, comme le Problème le demande
15. Commençons par déterminer
![{\displaystyle x,y,z}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bbeca34b28f569a407ef74a955d041df9f360268)
en
![{\displaystyle r.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10110093812676dd04a92ce4c8b75940c366330a)
Les équations (A) et (E) donnent, en retenant
![{\displaystyle p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/81eac1e205430d1f40810df36a0edffdc367af36)
à la place de
![{\displaystyle x={\frac {gz-cp}{bg-cf}},\quad y={\frac {fz-bp}{cf-bg}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37a2b892e5eb06fe8f6253aff7f8c2a10c8c4380)
substituant ces valeurs dans
et ordonnant par rapport à
on a
![{\displaystyle \left[(cf-bg)^{2}+f^{2}+g^{2}\right]z^{2}-2(bf+cg)pz+\left(b^{2}+c^{2}\right)p^{2}-(cf-bg)^{2}r^{2}=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/666c9b94ff534f3593b69495832bb6e2b03ec504)
d’où l’on tirera
et ensuite
et ![{\displaystyle y.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83f72471aff7c6fbb27df0f971283a068efe091f)
Si l’on fait, pour plus de simplicité, on trouvera
![{\displaystyle q={\sqrt {\left[(cf-bg)^{2}+f^{2}+g^{2}\right]r^{2}-\left(1+b^{2}+c^{2}\right)p^{2}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f42e7f1bd28182865f6564ce91b7d093d901291)
on trouvera
(G)
|
|
|
16. Maintenant l’équation (B) donne, en chassant
par le moyen de l’équation (D),
![{\displaystyle dx^{2}+dy^{2}+dz^{2}={\frac {r-{\cfrac {r^{2}}{a}}}{r-{\cfrac {r^{2}}{a}}-h}}dr^{2}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/469757b3e78d26e74e98d29db7b1168bf8b19d33)
mais les équations précédentes donnent
![{\displaystyle dx^{2}+dy^{2}+dz^{2}={\frac {\left(1+b^{2}+c^{2}\right)dp^{2}+dq^{2}}{(cf-bg)^{2}+f^{2}+g^{2}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f75c13a690608dd6d31f15ccd1720c314d8d164e)
il faut donc que ces deux expressions de
deviennent identiques après qu’on aura substitué dans la dernière les valeurs de
et
en ![{\displaystyle r.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10110093812676dd04a92ce4c8b75940c366330a)
Ces substitutions faites, on verra que l’identité aura lieu en effet, en faisant
(H)
|
|
|
C’est l’équation de condition cherchée.
Si, dans l’expression de
du numéro précédent, on substitue la valeur de
donnée par l’équation (G), que nous venons de trouver, et qu’on y mette de plus pour
sa valeur
elle deviendra
![{\displaystyle q=2{\sqrt {h\left(1+b^{2}+c^{2}\right)}}{\sqrt {r-{\frac {r^{2}}{a}}-h}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3f7119c3d605458e5f07f73a1042b21fca500445)
17. Pour pouvoir appliquer les formules précédentes au mouvement des comètes, il faut connaître les valeurs des constantes que ces formules renferment.
Pour cet effet, je remarque d’abord que l’équation (A) est celle d’un plan dont la position, à l’égard du plan des coordonnées
et
dépend des constantes
et
Ce plan sera donc celui de l’orbite de la comète, et qui est déterminé par les observations.
Soient
l’angle que l’intersection des deux plans, c’est-à-dire, la ligne des nœuds de l’orbite sur le plan des
et
fait avec l’axe des
et
l’inclinaison de l’orbite sur ce dernier plan ; il est facile de prouver qu’on aura
![{\displaystyle c=\operatorname {tang} \psi \cos \omega ,\quad b=-\operatorname {tang} \psi \sin \omega .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0cfc5c0b142fac52bf4f3bdddedfcb332186f843)
L’équation (E) servira ensuite à déterminer la figure de l’orbite ; et il est aisé de conclure de la forme même de cette équation que l’orbite ne peut être qu’une section conique, ayant le foyer dans l’origine des coordonnées, en sorte que
sera le rayon vecteur de l’orbite.
Les deux apsides seront donc aux points où
or, dans ce cas, l’équation (D) donne
![{\displaystyle r-{\frac {r^{2}}{a}}-h=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42e3e03b44b19604a00af5edfd3007980f5bee91)
équation dont les deux racines sont
![{\displaystyle {\frac {a\pm {\sqrt {a^{2}-4ah}}}{2}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e3cf1f7a0aa294361d744dc5fa606ce5cb9d1e3)
La somme de ces deux racines sera le grand axe, et leur différence, divisée par la somme, sera l’excentricité. Donc le grand axe de l’orbite sera
et l’excentricité sera
que je désignerai dans la suite par
.
Puis donc que
![{\displaystyle e={\sqrt {1-{\frac {4h}{a}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be3de7189e418c24de5c99d4f077cb54954490f4)
on aura
![{\displaystyle a{\sqrt {1-e^{2}}}={\sqrt {4ha}}=}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f648810034d1e91ecb7ae49e2a9e530e56302bce)
au petit axe de l’orbite ;
par conséquent,
sera le paramètre du grand axe.
Or on sait que le rayon vecteur qui répond à
degrés d’anomalie, c’est-à-dire, qui est perpendiculaire au grand axe, est égal au demi-paramètre. Donc on aura à
degrés d’anomalie
et l’équation (E) donnera
d’où l’on tire
égal par conséquent à la tangente de l’angle que fait avec l’axe des
dans le plan des
et
la projection du rayon vecteur qui répond à
degrés d’anomalie dans l’orbite.
Soit cet angle
on aura donc
donc, faisant
on aura
ces valeurs étant substituées dans l’équation (H) du no 16, ainsi que celles de
et
trouvées ci-dessus, et mettant
à la place de
elle deviendra
![{\displaystyle e^{2}=l^{2}{\frac {1+\operatorname {tang} ^{2}\psi \cos ^{2}(\omega -\varepsilon )}{1+\operatorname {tang} ^{2}\psi }}=l^{2}\left[1-\sin ^{2}\psi \sin ^{2}(\omega -\varepsilon )\right]\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1e617f921cd5565559beeaa23f9852651cb198de)
d’où l’on tire
![{\displaystyle l={\frac {e}{\sqrt {1-\sin ^{2}\psi \sin ^{2}(\omega -\varepsilon )}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/89ea5e2d6c076b4a40d36db2d983722d3d06ba97)
de sorte qu’on aura
![{\displaystyle f={\frac {e\cos \varepsilon }{\sqrt {1-\sin ^{2}\psi \sin ^{2}(\omega -\varepsilon )}}},\quad g={\frac {e\sin \varepsilon }{\sqrt {1-\sin ^{2}\psi \sin ^{2}(\omega -\varepsilon )}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/317e3aa2ddaf9cc15458aff65c1cde38ce64354b)
18. Si l’on fait coïncider le plan de l’orbite avec celui de
et
on aura
et l’angle
sera évidemment celui que le grand axe de l’orbite fait avec l’axe des
Donc, si l’on suppose de plus que ces deux axes coïncident, on aura aussi
de sorte que, dans cette hypothèse,
et les formules (G) du no 15 donneront
![{\displaystyle x={\frac {p}{e}},\quad u={\frac {q}{e}},\quad z=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2236803862465e7a3256fb822f6ba6e758e20a80)
savoir
![{\displaystyle x={\frac {2h-r}{e}},\quad y={\frac {2{\sqrt {h}}}{e}}{\sqrt {r-{\frac {r^{2}}{a}}-h}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5653403661b83d90770535d21cd61b35021682af)
Or il est visible que, dans ce cas,
et
deviennent les coordonnées de l’orbite dans le plan même de cette orbite ; et comme ces coordonnées doivent être indépendantes de la position du plan de l’orbite, il s’ensuit que les valeurs précédentes de
et
exprimeront toujours, l’une l’abscisse prise dans le grand axe depuis le foyer, et l’autre l’ordonnée rectangle dans le plan même de l’orbite, quelle que soit d’ailleurs la position de ce plan.
Donc, si l’on nomme
l’angle du rayon vecteur
avec le grand axe, on aura, dans la supposition précédente,
![{\displaystyle x=r\cos \varphi ,\quad y=r\sin \varphi \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37e6af01a14aa50d65d5314c8aa7206bf83621f1)
savoir
![{\displaystyle r\cos \varphi ={\frac {2h-r}{e}},\quad r\sin \varphi ={\frac {2{\sqrt {h}}}{e}}{\sqrt {r-{\frac {r^{2}}{a}}-h}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83bb10f071f26ba0c890df96fec48fb789f914fc)
d’où l’on tire
![{\displaystyle r={\frac {2h}{1+e\cos \varphi }},\quad {\sqrt {r-{\frac {r^{2}}{a}}-h}}={\frac {e{\sqrt {h}}\sin \varphi }{1+e\cos \varphi }}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69f0f985a0238bdfa4315391424bc830ce3d9ac3)
cette expression de
fait voir que
est l’anomalie vraie de l’orbite, comptée de son périhélie.
On aura donc, en général,
![{\displaystyle p=er\cos \varphi ,\quad q={\frac {er\sin \varphi }{\cos \psi }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1f9f8dc04fbf4ee25301e8b9c01d7c01c2302574)
et l’on pourra, par ces substitutions, dans les formules (F) et (G), avoir les valeurs de
exprimées par la seule anomalie ![{\displaystyle \varphi .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b0b6c90c1e9984232aed2d530ac2fb2660ea000a)
19. Dans le nœud, on a
donc (équation G)
![{\displaystyle (bf+cg)p+(cf-bg)q=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/386cc573ddac1cf51f48a5bc3fd3590db7260433)
savoir, en substituant les valeurs précédentes de
et ![{\displaystyle q,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1aa74b799849683cad6a0b79ebd9bf58bdf9890a)
![{\displaystyle (bf+cg)\cos \varphi +(cf-bg){\frac {\sin \varphi }{\cos \psi }}=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bdb8a6f5f2daaf23a56ad378d6df4866d664d79e)
où
est l’anomalie qui répond au nœud.
Dénotons cette anomalie par
on aura donc
![{\displaystyle (bf+cg)\cos \alpha \cos \psi +(cf-bg)\sin \alpha =0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bbace99e3b5187b2b9fd2c9b6ebd71298f55de50)
d’où l’on tire
![{\displaystyle {\frac {g}{f}}={\frac {c\sin \alpha +b\cos \psi \cos \alpha }{b\sin \alpha -c\cos \psi \cos \alpha }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3577df892976279f87756b0bb10fa5f0f1f2892f)
et, en mettant pour
et
leurs valeurs (no 17),
![{\displaystyle {\frac {g}{f}}={\frac {-\cos \omega \sin \alpha +\cos \psi \sin \omega \cos \alpha }{\sin \omega \sin \alpha +\cos \psi \cos \omega \cos \alpha }}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d21d775bbb28c7495d0bad038b6763000423257)
mais on a trouvé plus haut (no 17)
ainsi l’on aura l’équation
![{\displaystyle \operatorname {tang} \varepsilon ={\frac {-\cos \omega \sin \alpha +\cos \psi \sin \omega \cos \alpha }{\sin \omega \sin \alpha +\cos \psi \cos \omega \cos \alpha }}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7eb11b40661b7e0620d8b1ff4435a3708ca267ed)
d’où il est aisé de tirer
![{\displaystyle \sin \varepsilon ={\frac {-\cos \omega \sin \alpha +\cos \psi \sin \omega \cos \alpha }{\sqrt {1-\sin ^{2}\psi \cos ^{2}\alpha }}},\quad \cos \varepsilon ={\frac {\sin \omega \sin \alpha +\cos \psi \cos \omega \cos \alpha }{\sqrt {1-\sin ^{2}\psi \cos ^{2}\alpha }}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ecaaa7641233b5645583e1aad2b0ed70b3272ee0)
et, en substituant ces valeurs dans les expressions de
![{\displaystyle f}](https://wikimedia.org/api/rest_v1/media/math/render/svg/132e57acb643253e7810ee9702d9581f159a1c61)
et
![{\displaystyle g}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d3556280e66fe2c0d0140df20935a6f057381d77)
du
no 17, on trouvera
![{\displaystyle f=e\left(\cos \omega \cos \alpha +{\frac {\sin \omega \sin \alpha }{\cos \psi }}\right),\quad g=e\left(\sin \omega \cos \alpha -{\frac {\cos \omega \sin \alpha }{\cos \psi }}\right)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8be0b280de495baea5c83e0049aa8c5e4d829104)
par là, et par les valeurs de
et
on aura
![{\displaystyle cf-bg=e\operatorname {tang} \psi \cos \alpha ,\quad bf+cg=-{\frac {e\operatorname {tang} \psi \sin \alpha }{\cos \psi }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/460edc153c110e7da630d70c06d71f8de27934e2)
![{\displaystyle f+c(cf-bg)={\frac {e(\cos \omega \cos \alpha +\sin \omega \sin \alpha \cos \psi )}{\cos ^{2}\psi }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7563930d306d21e496660ab1882fa0991ed6e11a)
![{\displaystyle g-b(cf-bg)={\frac {e(\sin \omega \cos \alpha -\cos \omega \sin \alpha \cos \psi )}{\cos ^{2}\psi }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ecf336e2ee1db61c2002211f193d8bf1bc1c5cc9)
![{\displaystyle (cf-bg)^{2}+f^{2}+g^{2}={\frac {e^{2}}{\cos ^{2}\psi }}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/772c07e5679ef523307a4d9d4bd28b0cf486985f)
de sorte que, à cause de
les formules (G) du no 15 deviendront
![{\displaystyle {\begin{aligned}x=&r\left[\cos \omega \cos(\varphi -\alpha )-\sin \omega \sin(\varphi -\alpha )\cos \psi \right],\\y=&r\left[\sin \omega \cos(\varphi -\alpha )+\cos \omega \sin(\varphi -\alpha )\cos \psi \right],\\z=&r\,\ \sin \psi \sin(\varphi -\alpha ),\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/36fbfa7e77cce2dc57bd1ad73cad0547c860cff1)
dans lesquelles
est ce qu’on nomme l’argument de latitude.
20. Si l’on fait
![{\displaystyle {\frac {2{\sqrt {r-{\cfrac {r^{2}}{a}}-h}}}{\sqrt {a-4h}}}=\sin u,\quad 2t{\sqrt {\frac {2(1+m)}{a^{3}}}}+i=\theta ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/779d6b0276bc964814ba119afc13430851ca048d)
et qu’on se souvienne que
(no 17), il est clair que l’équation (F) du no 15 prendra cette forme très-simple
![{\displaystyle u-e\sin u=\theta ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/22eaf1259bad89b6df7df04fde64faf61885d61a)
dans laquelle
sera évidemment ce que l’on nomme, d’après Kepler, l’anomalie excentrique, mais comptée du périhélie, et où
sera, par conséquent, l’anomalie moyenne.
Donc, comme
lorsque
on voit que la constante
n’est autre chose que l’époque de l’anomalie moyenne.
En appliquant les formules au mouvement de la Terre autour du Soleil, et prenant la distance moyenne
pour l’unité, on aura, en négligeant
vis-à-vis de
![{\displaystyle t+i=\theta =}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d0b6dc3658df15ecbd032d9d1806fe5b09116c2b)
l’anomalie moyenne du Soleil ;
d’où il s’ensuit que
exprimé en angles, représentera proprement le mouvement moyen du Soleil pendant le temps écoulé depuis l’époque d’où l’on part ; et qu’ainsi, en divisant la valeur de
par
degrés, ou bien, si
est exprimé en nombres (la distance moyenne du Soleil étant l’unité), en divisant la valeur de
par le rapport de la circonférence au rayon, on aura le temps exprimé en années sidérales, puisque nous pouvons faire abstraction ici du mouvement de l’apogée du Soleil.
Or, puisque
![{\displaystyle 1-{\frac {4h}{a}}-4\left({\frac {r}{a}}-{\frac {r^{2}}{a^{2}}}-{\frac {h}{a}}\right)=\left(1-{\frac {2r}{a}}\right)^{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e069a4f1e3a5306245bc1f03009bf08a9a3b7c66)
il est clair qu’on aura
![{\displaystyle 1-{\frac {2r}{a}}=e\cos u\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/232f6f02b09e8725973b8bd4480a7b8d93c797ef)
donc
![{\displaystyle r={\frac {a}{2}}(1-e\cos u),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7253e17e0445cb88187c1873f01d171d1a648181)
et comme (nos 15 et 16)
![{\displaystyle p=2h-r,\quad q=2{\sqrt {h\left(1+b^{2}+c^{2}\right)}}{\sqrt {r-{\frac {r^{2}}{a}}-h}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67c8f5e094411d8ada5f56858595f1db6e1e2b2c)
on aura, à cause de ![{\displaystyle h={\frac {a\left(1-e^{2}\right)}{4}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cfb1b491cd7e6355d496347beeca6970aee159b1)
![{\displaystyle {\begin{aligned}p=&{\frac {ae}{2}}(\cos u-e),\\q=&{\frac {ae}{2}}{\sqrt {\left(1-e^{2}\right)\left(1+b^{2}+c^{2}\right)}}\sin u.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a60522c12da792d4066e1ecde1ec159925490b68)
De sorte que, en substituant ces valeurs dans les formules (G), on aura aussi
exprimées en ![{\displaystyle u.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/edd5636410da69bac33da075162221527401793c)
Dans la parabole, le grand axe
devient infini, par conséquent l’angle
est infiniment petit. Dans les ellipses très-allongées, telles que sont les orbites des comètes, la quantité
est seulement très-grande ; donc l’angle
sera très-petit, du moins tant que
n’est pas fort grand.
Dans ce cas donc on aura
![{\displaystyle u=\sin u+{\frac {1}{6}}\sin ^{3}u+{\frac {3}{40}}\sin ^{5}+{\frac {5}{112}}\sin ^{7}u+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6f0e7089dfa61feefe2bc23c5ebbb6b7c264e093)
et l’équation entre
et
deviendra
![{\displaystyle \theta =(1-e)\sin u+{\frac {1}{6}}\sin ^{3}u+{\frac {3}{40}}\sin ^{5}+{\frac {5}{112}}\sin ^{7}u+\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/20873aa78c9eaaae6a70c2f4a6e824486403dec6)
mais
![{\displaystyle 1-e={\frac {1-e^{2}}{1+e}}={\frac {4h}{a(1+e)}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b42977f300684277b95b36b13e332676a2e2c3f0)
donc, si l’on met pour
sa valeur
![{\displaystyle 2t{\sqrt {\frac {2(1+m)}{a^{3}}}}+i,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b28409b2e7fe8aefe8262d5905f18dd6ee3240af)
et qu’on fasse
![{\displaystyle i{\sqrt {\frac {a^{3}}{8}}}=j,\quad \sin u={\frac {w}{\sqrt {\cfrac {a}{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1d826f97ef563b5ca2260b23ad5ef700e4b4c6e2)
on aura, après avoir tout multiplié par ![{\displaystyle {\sqrt {\frac {a^{3}}{8}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b980f4f29ffb7c56ccfa843e6a558d84523af8f2)
![{\displaystyle t{\sqrt {1+m}}+j={\frac {2h}{1+e}}w+{\frac {1}{6}}w^{3}+{\frac {3}{20a}}w^{5}+{\frac {5}{28a^{2}}}w^{7}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc66208be06e4a823af851737a4ef316730634cd)
où
sera une quantité finie, puisqu’on aura
![{\displaystyle w={\frac {{\sqrt {2}}{\sqrt {r-{\cfrac {r^{2}}{a}}-h}}}{\sqrt {1-{\cfrac {4h}{a}}}}}={\frac {{\sqrt {2}}{\sqrt {r-{\cfrac {r^{2}}{a}}-h}}}{e}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11aa2b1cc2b115cf985f80c3ef3bd687947ce8dc)
et, substituant pour
sa valeur en
trouvée ci-dessus, il
viendra
![{\displaystyle w={\frac {{\sqrt {2h}}\sin \varphi }{1+e\cos \varphi }}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94316238c00e3704e8934a9f52ccd6660710653f)
Pour la parabole, on fera
et l’on aura
![{\displaystyle t{\sqrt {1+m}}+j=hw+{\frac {w^{3}}{6}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e24123cca3a5b9a0a7742b6e733fbdbcdf7ce96)
où (à cause de
)
![{\displaystyle w={\sqrt {2(r-h)}}={\sqrt {2h}}\operatorname {tang} {\frac {\varphi }{2}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0fad3916350f97cfa7fd489c69729b92f9906020)
et
sera pour lors égal à la distance périhélie.
21. Nous remarquerons encore que, si, dans l’équation différentielle entre
et
du no 15, on substitue pour
et pour
leurs valeurs en
(no 18), on trouve
![{\displaystyle dt={\sqrt {\frac {1}{2h(1+m)}}}r^{2}d\varphi \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92465ae7003080b3617a852ab6444db3049dee4b)
et, si l’on différentie l’équation qui donne la relation entre
et
(no 20), et qu’on y mette
pour
il vient
![{\displaystyle dt={\sqrt {\frac {a}{2(1+m)}}}rdu\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8b9dd7a1a4ac04337d4222b574bcc8e859962531)
dans la première formule,
est l’anomalie vraie, et dans la seconde
est l’anomalie excentrique.
section troisième.
intégration des équations différentielles des perturbations.
22. Nous avons vu, dans la première Section, que
étant les coordonnées de l’orbite non altérée, et
celles de l’orbite troublée par l’action d’une planète
on a pour la détermination des quantités
les équations suivantes
![{\displaystyle {\begin{aligned}{\frac {d^{2}\delta x}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\delta z\right)-\mu \mathrm {X} ,\\{\frac {d^{2}\delta y}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dydx}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dydz}}\delta z\right)-\mu \mathrm {Y} ,\\{\frac {d^{2}\delta z}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dzdx}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dzdy}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\delta z\right)-\mu \mathrm {Z} ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4feba46d95948a0013026b068b3698d6431aa30e)
en faisant, pour abréger (no 12),
![{\displaystyle \mathrm {X} ={\frac {d\left({\cfrac {1}{\sigma }}-{\cfrac {1}{\mathrm {R} }}\right)}{dx}},\quad \mathrm {Y} ={\frac {d\left({\cfrac {1}{\sigma }}-{\cfrac {1}{\mathrm {R} }}\right)}{dy}},\quad \mathrm {Z} ={\frac {d\left({\cfrac {1}{\sigma }}-{\cfrac {1}{\mathrm {R} }}\right)}{dz}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7feff74bd86646893ddbf0e4b0081ce0f8033b38)
23. C’est donc de l’intégration de ces équations que dépend la recherche des perturbations causées par l’action d’une planète quelconque sur la comète. Or cette intégration dépend, comme l’on sait, de celle du cas où il n’y aurait aucun terme tout connu, à cause que les variables inconnues
ne paraissent que sous la forme linéaire. Ainsi toute la difficulté se réduit à intégrer des équations de la forme suivante
![{\displaystyle {\begin{aligned}{\frac {d^{2}\delta x}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dx^{2}}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dxdy}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dxdz}}\delta z\right),\\{\frac {d^{2}\delta y}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dydx}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dy^{2}}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dydz}}\delta z\right),\\{\frac {d^{2}\delta z}{dt^{2}}}=&(1+m)\left({\frac {d^{2}{\frac {1}{r}}}{dzdx}}\delta x+{\frac {d^{2}{\frac {1}{r}}}{dzdy}}\delta y+{\frac {d^{2}{\frac {1}{r}}}{dz^{2}}}\delta z\right).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/311f38cb4a63525b67e15ad0add85350c76f3750)
24. Si l’on se rappelle les calculs du no 7, on doit voir que les équations précédentes résultent des équations de l’orbite non altérée, en y faisant varier les quantités
des différences
regardées comme infiniment petites. Donc les intégrales des équations dont il s’agit doivent résulter aussi des intégrales des mêmes équations de l’orbite non altérée, en y faisant varier non-seulement ces mêmes quantités, mais encore les constantes arbitraires introduites par les différentes intégrations, et qui, n’existant point dans les équations différentielles, peuvent, à leur égard, être aussi regardées comme variables.
Ainsi donc, pour avoir les intégrales des trois équations différentielles du numéro précédent, il n’y aura qu’à différentier à l’ordinaire les intégrales de l’orbite non altérée, trouvées dans la seconde Section, en y regardant les trois indéterminées
et les six arbitraires
comme variables à la fois, et marquant leurs différences par la caractéristique
[à l’égard de
elle doit aussi être traitée comme variable, parce que c’est une fonction de
donnée par l’équation (H) du no 16] ; les différences de ces arbitraires seront elles-mêmes les nouvelles constantes arbitraires que les intégrales cherchées doivent contenir pour être complètes.
25. Comme les formules (G) du no 15 donnent
en
et que la formule (F) du no 14 donner en
on pourra tirer directement de la différentiation des premières les valeurs de
en
ensuite on aura
par la différentiation de la dernière ; mais, à la place de
il sera plus simple d’introduire l’angle
au moyen des formules du no 20 ; et, pour donner à notre calcul toute la simplicité dont il est susceptible, nous remarquerons de plus que, la position du plan des
et
auquel nous avons jusqu’ici rapporté l’orbite de la comète, étant arbitraire, on peut, sans nuire à la généralité du Problème, supposer que ce plan coïncide avec celui de l’orbite non altérée de la comète ; et l’on peut, par la même raison, supposer aussi que l’axe des
coïncide avec le grand axe de la même orbite, en sorte que les abscisses
soient prises depuis le foyer et soient positives en allant vers l’apside inférieure.
Ces deux suppositions donneront (nos 17 et 19)
![{\displaystyle \psi =0,\quad \alpha =\omega \,;\quad {\text{donc}}\quad b=0,\quad c=0,\quad f=e,\quad g=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/261100ad0ab8e61ecf2e15febe9f0ba30ef64917)
ce qui simplifiera beaucoup les expressions finies de
![{\displaystyle x,y,z\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/41656228f52d6c7b999c2944e813136d460b0102)
mais, comme les différences
![{\displaystyle \delta b,\delta c,\delta g}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e8ea7caed61ff0a690300c2ecc172fe2a6de886)
ne sont pas nulles, il ne faudra pas faire disparaître entièrement les quantités
![{\displaystyle b,c,g\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ce366091e1e8f1a726ceafb04f88e62b4853e30)
mais il en faudra conserver lespremières dimensions dans les expressions de
![{\displaystyle x,y,z,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d08d690d7e19ea7aee8574fc6abd6a15d97fa026)
afin de pouvoir en tirer par la différentiation les valeurs complètes de
26. De cette manière, on aura donc (no 15, formule G)
![{\displaystyle x={\frac {fp-gq}{f^{2}}},\quad y={\frac {gp+fq}{f^{2}}},\quad z={\frac {bp+cq}{f}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/daafdf850bd9efabf3d47a05bcf07dd93d3ab081)
et par les formules du no 20 on aura, à cause de ![{\displaystyle e=f,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e2276f0378d4922fe51732593d7a1c3aa4af48c)
![{\displaystyle p={\frac {af}{2}}(\cos u-f),\quad q={\frac {af}{2}}{\sqrt {1-f^{2}}}\sin u,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/320150889d365e68a05518e5cf64aaf62a65c477)
de sorte qu’en substituant il viendra
![{\displaystyle {\begin{aligned}x=&{\frac {a}{2}}(\cos u-f)-{\frac {ag}{2f}}{\sqrt {1-f^{2}}}\sin u,\\y=&{\frac {a}{2}}{\sqrt {1-f^{2}}}\sin u+{\frac {ag}{2f}}(\cos u-f),\\z=&{\frac {ab}{2}}(\cos u-f)+{\frac {ac}{2}}{\sqrt {1-f^{2}}}\sin u.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a04a5867184cd50b7dc9abc8a120a30609ddf23c)
Différentiant donc suivant la caractéristique
en faisant tout varier, et supposant ensuite les constantes
nulles, on aura
![{\displaystyle {\begin{aligned}\delta x=&{\frac {\cos u-f}{2}}\delta a-{\frac {a}{2}}\delta f-{\frac {a{\sqrt {1-f^{2}}}}{2f}}\sin u\,\delta g-{\frac {a}{2}}\sin u\,\delta u,\\\delta y=&{\frac {\sqrt {1-f^{2}}}{2}}\sin u\,\delta a-{\frac {af}{2{\sqrt {1-f}}}}\sin u\,\delta f+{\frac {a(\cos u-f)}{2f}}\delta g\\&\ \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad +{\frac {a{\sqrt {1-f^{2}}}}{2}}\cos u\,\delta u,\\\delta z=&{\frac {a(\cos u-f)}{2}}\delta b+{\frac {a{\sqrt {1-f^{2}}}}{2}}\sin u\,\delta c.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/59155b22f977772136de7e5b30dcc4c029f75628)
Mais, par le no 20, on a, en mettant
à la place de ![{\displaystyle e,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1cfc3e7aeae6ef256412eaa9f68d097804943b5a)
![{\displaystyle u-f\sin u=\theta =2t{\sqrt {\frac {2(1+m)}{a^{3}}}}+i\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9c72c1177fc5f57135a1618bf12964caa862bcad)
donc, faisant varier
![{\displaystyle f,a,i}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7425c23485a62fbf159d527eb36ac376b15f770d)
et
![{\displaystyle u,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/30dcc93e14b40416ed2d1391bc6c08ee99fa5ff6)
on en tirera la valeur de
![{\displaystyle \delta u,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c8e9b0f8b7ee6e8dfe0dc8d413f192c348d4c90f)
laquelle sera
![{\displaystyle \delta u={\frac {-3t{\sqrt {\cfrac {2(1+m)}{a^{3}}}}{\cfrac {\delta a}{a}}+\sin u\,\delta f+\delta i}{1-f\cos u}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47def6c144cc494c522380dc75a3d876e08ed3bb)
Substituant donc cette valeur de
dans les expressions précédentes de
on aura les valeurs cherchées, lesquelles seront évidemment de cette forme
![{\displaystyle {\begin{aligned}\delta x=&\mathrm {A} \delta a+\mathrm {B} \delta f+\mathrm {C} \delta g+\mathrm {D} \delta i,\\\delta y=&\mathrm {E} \delta a\,+\mathrm {F} \delta f+\mathrm {G} \delta g+\mathrm {H} \delta i,\\\delta z=&\mathrm {K} \delta b\,+\mathrm {L} \delta c,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/970504ea26c4f5fdb64ffe4b3f71741cdb750961)
en supposant, pour abréger,
![{\displaystyle {\begin{aligned}\mathrm {A} =&{\frac {\cos u-f}{2}}+{\frac {3t}{2}}{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\sin u}{1-f\cos u}},\\\mathrm {B} =&-{\frac {a}{2}}\left(1+{\frac {\sin ^{2}u}{1-f\cos u}}\right),\\\mathrm {C} =&-{\frac {a{\sqrt {1-f^{2}}}}{2f}}\sin u,\\\mathrm {D} =&-{\frac {a}{2}}{\frac {\sin u}{1-f\cos u}},\\\mathrm {E} =&{\frac {\sqrt {1-f^{2}}}{2}}\sin u-{\frac {3t}{2}}{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {{\sqrt {1-f^{2}}}\cos u}{1-f\cos u}},\\\mathrm {F} =&-{\frac {af}{2{\sqrt {1-f^{2}}}}}\sin u+{\frac {a{\sqrt {1-f^{2}}}}{2}}{\frac {\sin u\cos u}{1-f\cos u}},\\\mathrm {G} =&{\frac {a}{2f}}(\cos u-f),\\\mathrm {H} =&{\frac {a{\sqrt {1-f^{2}}}}{2}}{\frac {\cos u}{1-f\cos u}},\\\mathrm {K} =&{\frac {a}{2}}(\cos u-f),\\\mathrm {L} =&{\frac {a{\sqrt {1-f^{2}}}}{2}}\sin u.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/031ed59808325dae9d62970e8536e82826bb41a6)
Telles sont les valeurs complètes des quantités
en tant qu’elles résultent des trois équations différentielles du no 23 ; et les quantités
sont les six constantes arbitraires que ces valeurs doivent contenir à raison des six intégrations qu’elles supposent.
27. Voyons maintenant comment on doit déterminer ces nouvelles arbitraires il est clair qu’elles dépendent des valeurs des quantités
et de leurs différences premières
pour un instant quelconque donné. Il faudra donc différentier les expressions de
trouvées ci-dessus, en y regardant les arbitraires
comme constantes, c’est-à-dire, en y faisant varier seulement les quantités qui sont des fonctions du temps
pour avoir les valeurs de
lesquelles seront représentées, en général, par les formules suivantes
![{\displaystyle {\begin{aligned}{\frac {d\delta x}{dt}}=&{\frac {d\mathrm {A} }{dt}}\delta a+{\frac {d\mathrm {B} }{dt}}\delta f+{\frac {d\mathrm {C} }{dt}}\delta g+{\frac {d\mathrm {D} }{dt}}\delta i,\\{\frac {d\delta y}{dt}}=&{\frac {d\mathrm {E} }{dt}}\delta a\,+{\frac {d\mathrm {F} }{dt}}\delta f+{\frac {d\mathrm {G} }{dt}}\delta g+{\frac {d\mathrm {H} }{dt}}\delta i,\\{\frac {d\delta z}{dt}}=&{\frac {d\mathrm {K} }{dt}}\delta b\,+{\frac {d\mathrm {L} }{dt}}\delta c,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b6338fff66ea4a3c781cc54b5c0a4e1e01e21ab9)
Ces trois équations étant combinées avec les trois du numéro précédent, on en tirera, par la méthode ordinaire d’élimination, les valeurs des six inconnues
et il est aisé de voir que ces valeurs seront de la forme suivante
![{\displaystyle {\begin{aligned}\delta a=&\mathrm {A} '\delta x+\mathrm {B} '\delta y\,+\mathrm {C} '{\frac {d\delta x}{dt}}\,+\mathrm {D} '{\frac {d\delta y}{dt}},\\\delta f=&\mathrm {E} '\delta x\,+\mathrm {F} '\delta y\,+\mathrm {G} '{\frac {d\delta x}{dt}}\,+\mathrm {H} '{\frac {d\delta y}{dt}},\\\delta g=&\mathrm {A} ''\delta x+\mathrm {B} ''\delta y+\mathrm {C} ''{\frac {d\delta x}{dt}}+\mathrm {D} ''{\frac {d\delta y}{dt}},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce1dec29a4a5d171bbb91e7949a888946b552b1a)
![{\displaystyle {\begin{aligned}\delta i=&\mathrm {E} ''\delta x+\mathrm {F} ''\delta y\,+\mathrm {G} ''{\frac {d\delta x}{dt}}\,+\mathrm {H} ''{\frac {d\delta y}{dt}},\\\delta b=&\mathrm {K} '\delta z\,+\mathrm {L} '{\frac {d\delta z}{dt}},\\\delta c=&\mathrm {K} ''\delta z+\mathrm {L} ''{\frac {d\delta z}{dt}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1bcf78dbb7fc5a78a1b97f52e5c64e4a84efd77e)
les quantités
étant des fonctions rationnelles de
et de ![{\displaystyle {\frac {d\mathrm {A} }{dt}},{\frac {d\mathrm {B} }{dt}},\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/97c138802a59524c92fe7bf2f1b19d054fbb6ade)
28. Quoique la détermination de ces quantités
ne soit pas difficile, elle pourrait néanmoins entraîner dans des calculs très-longs, si on l’entreprenait par la méthode ordinaire ; voici un moyen de la simplifier beaucoup.
Ce moyen consiste à chercher d’abord les valeurs des constantes
en
et en
à quoi on parviendra facilement par le moyen des formules du no 14 ; ensuite à différentier ces valeurs relativement à la caractéristique
c’est-à-dire, en faisant varier seulement les constantes dont il s’agit et les indéterminées
et marquant les variations par
et comme les différentiations relatives aux deux caractéristiques différentes
et
sont totalement indépendantes entre elles, on voit aisément que
sera la même chose que
de sorte qu’on aura ainsi directement les valeurs de
en
Nous allons donner ici les résultats de ce calcul, parce qu’ils nous seront utiles dans la suite.
29. On voit d’abord (no 14) que l’équation (B) donnera la valeur de
et que l’équation (D) donnera celle de
ensuite l’équation finie (E), combinée avec sa différentielle, donnera les valeurs de
et
et de même l’équation (A), combinée avec sa différentielle, donnera celle de
et
enfin l’équation (F) donnera la valeur de
on aura donc d’abord
![{\displaystyle {\frac {1}{a}}={\frac {1}{r}}-{\frac {dx^{2}+dy^{2}+dz^{2}}{2(1+m)dt^{2}}},\quad h=r-{\frac {r^{2}}{a}}-{\frac {\left(dr^{2}\right)^{2}}{8(1+m)dt^{2}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f5e35193b2ed7bc4f6ac2276199a690658fbb2f)
ensuite on trouvera
![{\displaystyle {\begin{alignedat}{2}f=&{\frac {(2h-r)dy+ydr}{xdy-ydx}},\qquad &g=&{\frac {(2h-r)dx+xdr}{ydx-xdy}},\\b=&{\frac {zdy-ydz}{xdy-ydx}},&c=&{\frac {zdx-xdz}{ydx-xdy}}\,;\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ef159d0347e1973ed85114e44806c7623809ab3)
or, si dans la valeur de
on substitue celle de
et qu’on y mette
à la place de
on a
![{\displaystyle h={\frac {r^{2}\left(dx^{2}+dy^{2}+dz^{2}\right)-(xdx+ydy+zdz)^{2}}{2(1+m)dt^{2}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f36a86b18e0014dcb507734f69e58e516111e2ce)
ce qui, à cause de
peut se réduire à cette forme
![{\displaystyle h={\frac {(xdy-ydx)^{2}+(zdy-ydz)^{2}+(xdz-zdx)^{2}}{2(1+m)dt^{2}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cf13e62cc316ccf8eba1ce04785c31adb18ec32b)
mais on vient de trouver
![{\displaystyle zdy-ydz=b(xdy-ydx),\quad xdz-zdx=c(xdy-ydx)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/06d4a63489f4ab29add5b713f8cd10f6189312fa)
donc, faisant ces substitutions, extrayant la racine carrée et supposant, pour abréger,
![{\displaystyle {\frac {h}{1+b^{2}+c^{2}}}=\mathrm {K} ^{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7dd128c7987e64262158d50d5dd210d65ef0b7a7)
on aura
![{\displaystyle \mathrm {K} ={\frac {xdy-ydx}{dt{\sqrt {2(1+m)}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/df8d5da24247d83d8301fbca374e352cff64b8d5)
et les autres formules deviendront, étant multipliées par ![{\displaystyle \mathrm {K} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ea6c3420085e49c56993db435a6d9fbd04a5300d)
![{\displaystyle {\begin{alignedat}{2}\mathrm {K} f=&{\frac {(2h-r)dy+ydr}{dt{\sqrt {2(1+m)}}}},\qquad &\mathrm {K} g=&-{\frac {(2h-r)dx+xdr}{dt{\sqrt {2(1+m)}}}},\\\mathrm {K} b=&{\frac {zdy-ydz}{dt{\sqrt {2(1+m)}}}},&\mathrm {K} c=&-{\frac {zdx-xdz}{dt{\sqrt {2(1+m)}}}}.\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d7254e3468af0a3d68f102a8501a09a7e0b4ab2e)
Enfin on aura (formule F)
![{\displaystyle i=-2t{\sqrt {\frac {2(1+m)}{a^{3}}}}+\operatorname {arc} \,\sin {\frac {2{\sqrt {r-{\cfrac {r^{2}}{a}}-h}}}{\sqrt {a-4h}}}-{\frac {2{\sqrt {r-{\cfrac {r^{2}}{a}}-h}}}{\sqrt {a}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4ee24fad3bbe6f535d92c0834395386aa91c40a9)
30. En différentiant ces équations par rapport à la caractéristique
et changeant partout
en
on trouvera les formules suivantes
![{\displaystyle {\begin{aligned}\delta a=&a^{2}\left({\frac {\delta r}{r^{2}}}+{\frac {dx\,d\delta x+dy\,d\delta y+dz\,d\delta z}{(1+m)dt^{2}}}\right),\\\delta \mathrm {K} =&\quad {\frac {dy\,\delta x-dx\,\delta y+x\,d\,\delta y-y\,d\,\delta x}{dt{\sqrt {2(1+m)}}}},\\\delta f=&\quad {\frac {dy(2\delta h-\delta r)+dr\,\delta y+(2h-r)d\,\delta y+y\,d\,\delta r}{\mathrm {K} dt{\sqrt {2(1+m)}}}}-{\frac {f\delta \mathrm {K} }{\mathrm {K} }},\\\delta g=&-{\frac {dx(2\delta h-\delta r)+dr\,\delta x+(2h-r)d\,\delta x+x\,d\,\delta r}{\mathrm {K} dt{\sqrt {2(1+m)}}}}-{\frac {g\delta \mathrm {K} }{\mathrm {K} }},\\\delta b=&\quad {\frac {dy\delta z-dz\delta y+z\,d\,\delta y-y\,d\,\delta z}{\mathrm {K} dt{\sqrt {2(1+m)}}}}-{\frac {b\delta \mathrm {K} }{\mathrm {K} }},\\\delta c=&-{\frac {dx\delta z-dz\delta x+z\,d\,\delta x-x\,d\,\delta z}{\mathrm {K} dt{\sqrt {2(1+m)}}}}-{\frac {c\delta \mathrm {K} }{\mathrm {K} }},\\\delta h=&2\left(1+b^{2}+c^{2}\right)\mathrm {K\delta K+2K^{2}} (b\delta b+c\delta c),\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5dc5f7ebbffd502bfc76c0246633d7ce7b0f0ed0)
dans lesquelles il faudra substituer pour
sa valeur
![{\displaystyle {\frac {x\delta x+y\delta y+z\delta z}{r}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b3cff08dba50aed1e909a3276e92b5706f069f4a)
et pour ![{\displaystyle ds\,\delta r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c9d52dd26a22f5c878465e733071ac339f0decd4)
![{\displaystyle {\frac {x\,d\,\delta x+y\,d\,\delta y+z\,d\,\delta z}{r}}+d{\frac {x}{r}}\delta x+d{\frac {y}{r}}\delta y+d{\frac {z}{r}}\delta z.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/807e44438cefa58d0e8bfe5e3ca73f63ef165692)
Quant à la valeur de
pour la trouver plus aisément, on supposera
![{\displaystyle {\frac {r}{a}}=v,\quad {\frac {h}{a}}=n,\quad {\frac {2{\sqrt {v-v^{2}-n}}}{\sqrt {1-4n}}}=\mathrm {V} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ca549e4f98fefd2a0f91060fe0abbc7e7936bc45)
ee qui réduira la valeur de
![{\displaystyle i}](https://wikimedia.org/api/rest_v1/media/math/render/svg/add78d8608ad86e54951b8c8bd6c8d8416533d20)
à cette forme
![{\displaystyle i=-2t{\sqrt {\frac {2(1+m)}{a^{3}}}}+\operatorname {arc} \,\sin \mathrm {V-V} {\sqrt {1-4n}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f5aa5dcc2c82a7738747f5cb0f0bd6f683a8e0d)
de sorte qu’on aura, en différentiant suivant
et faisant tout varier, excepté ![{\displaystyle t,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4ea3ad87830a1055c7b85c04cf940cfd3b847ae6)
![{\displaystyle \delta i=3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}+\mathrm {\frac {\delta V}{\sqrt {1-V^{2}}}} -{\sqrt {1-4n}}\,\delta \mathrm {V-V} \delta {\sqrt {1-4n}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e7b10744d862d37104857d9e0aae33ef3a030695)
or
![{\displaystyle {\sqrt {1-\mathrm {V} ^{2}}}={\frac {1-2v}{\sqrt {1-4n}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3eff8f2d0e28232c63cba0f57a99d2c5d5d4d1c0)
donc, substituant cette valeur, ainsi que celles de
et de
et réduisant, il viendra
![{\displaystyle \delta i=3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}+{\frac {2v\delta v+{\cfrac {2\delta n}{1-4n}}(v-2n)}{\sqrt {v-v^{2}-n}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67e348daa0b9471856e491c0c8c46e36122fed89)
où il n’y aura plus qu’à remettre pour
et
leurs valeurs
et par conséquent pour
et
les quantités ![{\displaystyle {\frac {a\delta r-r\delta a}{a^{2}}},\ {\frac {a\delta h-h\delta a}{a^{2}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7ede8773c1e3c303350357c8caeefb4d5e9efe5a)
Après avoir trouvé cette expression de
j’ai remarqué qu’elle avait l’inconvénient de contenir au dénominateur le radical
savoir
![{\displaystyle {\frac {\sqrt {r-{\cfrac {r^{2}}{a}}-h}}{\sqrt {a}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f7ebb8d0252882eef773faf9378b577db2ad6f65)
lequel, comme on l’a vu dans le, no 17, devient nul dans les deux apsides de l’orbite ; de sorte que, comme
ne saurait devenir infini, il faut nécessairement que le numérateur soit alors pareillement nul ; d’où il s’ensuit que la formule sera en défaut dans ces deux points.
Pour éviter cet inconvénient, il faut tâcher de donner une autre forme à la valeur de
et qui soit telle, qu’il n’y ait aucune fonction des variables au dénominateur ; voici comment je suis parvenu à ce but.
Je considère que la quantité
![{\displaystyle \mathrm {\frac {\delta V}{\sqrt {1-V^{2}}}} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/81e43e779f68221edea1eb555dd419b2859e304d)
est la même chose que celle-ci
![{\displaystyle \mathrm {{\sqrt {1-V^{2}}}\,\delta V-V\delta {\sqrt {1-V^{2}}}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/99e33391d668aac19149f414a187bfd12f5b7d88)
et qu’ainsi on peut réduire la première expression de
à celle-ci
![{\displaystyle \delta i=3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}+\left({\sqrt {1-\mathrm {V} ^{2}}}-{\sqrt {1-4n}}\right)\delta \mathrm {V} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/19d13f2489c5f9b474bb0214a8da520c6e2b84ab)
![{\displaystyle -\mathrm {V} \delta \left({\sqrt {1-\mathrm {V} ^{2}}}+{\sqrt {1-4n}}\right)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/599021e941cc16f4691edd7b34e8d89eb7c6bfd5)
mais
![{\displaystyle \mathrm {V} ={\frac {2{\sqrt {v-v^{2}-n}}}{\sqrt {1-4n}}}={\frac {2{\sqrt {r-{\cfrac {r^{2}}{a}}-h}}}{{\sqrt {a}}{\sqrt {1-4n}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53c2764adec02908fbc68eefd85fb5990866433b)
et, par le no 16 et à cause de ![{\displaystyle h=na,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7bc91bb895ed8b0233cdd28d25a6fe97dbe8d68d)
![{\displaystyle \mathrm {V} ={\frac {q}{a{\sqrt {1-4n}}{\sqrt {n\left(1+b^{2}+c^{2}\right)}}}},\quad {\sqrt {1-\mathrm {V} ^{2}}}={\frac {1-2v}{\sqrt {1-4n}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/00796042125965d39d78a712f08f58e5e6f9c455)
par conséquent (no 14)
![{\displaystyle {\sqrt {1-\mathrm {V} ^{2}}}-{\sqrt {1-4n}}={\frac {4n-2v}{\sqrt {1-4n}}}={\frac {4h-2r}{a{\sqrt {1-4n}}}}={\frac {2p}{a{\sqrt {1-4n}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/85d6a43d9a9edbc44dffab63d433b52b25b5cfe5)
![{\displaystyle {\sqrt {1-\mathrm {V} ^{2}}}+{\sqrt {1-4n}}={\frac {2p}{a{\sqrt {1-4n}}}}+2{\sqrt {1-4n}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c3dcef0341475562f59c4a23c96eee43687da681)
donc, faisant ces substitutions dans l’équation précédente, et effaçant les termes qui viendront de la différentiation de la quantité
![{\displaystyle a{\sqrt {1-4n}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ea243bf44e06122e3dba769abe1860cf794159a)
laquelle divise
et
parce que ces termes se détruisent mutuellement, on aura
![{\displaystyle \delta i=3t{\frac {\sqrt {2(1+m)}}{a^{3}}}{\frac {\delta a}{a}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7b9ea44644ef836c01b2835b2d782b0563ec9592)
![{\displaystyle +{\frac {2}{a^{2}(1-4n)}}\left[p\delta {\frac {q}{\sqrt {n\left(1+b^{2}+c^{2}\right)}}}-{\frac {q}{\sqrt {n\left(1+b^{2}+c^{2}\right)}}}(\delta p-2a\,\delta n)\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/21ef126c7540e75d19469ccab719ff7c55af3f44)
Or les formules (G) du no 15 donnent
![{\displaystyle {\begin{aligned}p=&fx+gy,\\q=&\left[f+c(cf-bg)\right]y-\left[g-b(cf-bg)\right]x\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42daefcf29f97069a5000cc88e1807bb8d71aa46)
donc, substituant ces valeurs, on aura pour
une expression toute rationnelle et entière, et qui ne sera par conséquent sujette à aucun inconvénient.
On remarquera encore, à l’égard de
qu’on peut aussi les exprimer d’une manière plus simple et plus commode à quelques égards, en les déduisant directement des équations (A) et (E), différentiées d’abord relativement à la caractéristique
et ensuite par rapport à la caractéristique ordinaire
ce qui, dans le fond, revient au même que si on les différentie d’abord par rapport à cette dernière caractéristique, et ensuite par rapport à la première, ainsi que nous en avons usé plus haut.
De cette manière, l’équation (E) donnera ces deux-ci
![{\displaystyle 2\delta h-\delta r=x\delta f+y\delta g+f\delta x+g\delta y,\quad -d\,\delta r=dx\delta f+dy\delta g+fd\,\delta x+gd\,\delta y\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc2bbb8764f347bafadcd44d7cefcf7ea5452555)
d’où l’on tire, en mettant pour
sa valeur ![{\displaystyle \mathrm {K} dt{\sqrt {2(1+m)}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/25618776762c07b96e9743f09cb3403fc164f918)
![{\displaystyle {\begin{aligned}\delta f=&\quad {\frac {dy(2\delta h-\delta r)+yd\,\delta r-f(dy\delta x-yd\,\delta x)-g(dy\delta y-yd\,\delta y)}{\mathrm {K} dt{\sqrt {2(1+m)}}}}\\\delta g=&-{\frac {dx(2\delta h-\delta r)+xd\,\delta r-f(dx\delta x-xd\,\delta x)-g(dx\delta y-xd\,\delta y)}{\mathrm {K} dt{\sqrt {2(1+m)}}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/22f2712225cbd3a1b48d476b80104b29ca3d2a37)
De même l’équation (A) donnera ces deux-ci
![{\displaystyle \delta z=x\delta b+y\delta c+b\delta x+c\delta y,\quad d\,\delta z=dx\delta b+dy\delta c+bd\,\delta x+cd\,\delta y\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6b45f7556d869c1e99d76df76b9d738074161c1c)
d’où l’on tire
![{\displaystyle {\begin{aligned}\delta b=&{\frac {dy\delta z-yd\,\delta z-b(dy\delta x-yd\,\delta x)-c(dy\delta y-yd\,\delta y)}{\mathrm {K} dt{\sqrt {2(1+m)}}}}\\\delta c=&-{\frac {dx\delta z-xd\,\delta z-b(dx\delta x-xd\,\delta x)-c(dx\delta y-xd\,\delta y)}{\mathrm {K} dt{\sqrt {2(1+m)}}}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/721263583a451d951af21009cc03daac679234c7)
ce sont les formules que nous emploierons dans la suite, par préférence aux autres trouvées ci-dessus.
Enfin on observera que, comme
on aura, par la formule (H) du no 16,
![{\displaystyle 4n=1-{\frac {(cf-bg)^{2}+f^{2}+g^{2}}{1+b^{2}+c^{2}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b7e66da182ebc3d1d041db245f58458cfc4a87d5)
de sorte que, en différentiant suivant
on aura la valeur de
exprimée à volonté par
et
ou par ![{\displaystyle \delta f,\delta g,\delta b,\delta c.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/721e25bc7dd26f9cd8f82f1c66cffbde19227d51)
31. Les formules précédentes ont toute la généralité possiblé ; mais, pour les appliquer à notre cas, il y faut supposer (no 25)
![{\displaystyle b=0,\quad c=0,\quad g=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/254f615d6292689b120ba78ae4fda8ecee8dba90)
ce qui donne [équation (G)]
![{\displaystyle z=0,\quad {\frac {dz}{dt}}=0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/643f001eab8edd5f4128553b68400cc9bca78bce)
par conséquent
![{\displaystyle d{\frac {x}{r}}=-{\frac {(x\,dy-y\,dx)y}{r^{3}}},\quad d{\frac {y}{r}}={\frac {(x\,dy-y\,dx)x}{r^{3}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9be1beb50ddcc576002404b15efb710f9d0a56d3)
donc
![{\displaystyle \delta r={\frac {x\delta x+y\delta y}{r}},\quad d\,\delta r={\frac {x\,d\,\delta x+y\,d\,\delta y}{r}}+{\frac {(x\delta y-y\delta x)(x\,dy-y\,dx)}{r^{3}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a9af0470527a224a3fc7ef41f76ca2ff155f893)
de plus on aura
![{\displaystyle \mathrm {K} ={\sqrt {h}},\quad \delta \mathrm {K} ={\frac {\delta h}{2{\sqrt {h}}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bafa283f9dd66dba76af2ab6c848b85697c5244e)
donc, faisant ces différentes réductions dans les formules ci-dessus, elles deviendront
![{\displaystyle {\begin{aligned}\delta h=&2h{\frac {dy\delta x-dx\delta y+x\,d\,\delta y-y\,d\,\delta x}{dt{\sqrt {2h(1+m)}}}}\\\delta a=&a^{2}\left({\frac {x\delta x+y\delta y}{r^{3}}}+{\frac {dx\,d\,\delta x+dy\,d\,\delta y}{(1+m)dt^{2}}}\right),\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/35281e15a6a3de97c8a86f79871b6d24a1c3e84f)
![{\displaystyle {\begin{aligned}\delta f=&{\frac {y}{r^{3}}}(x\delta y-y\delta x)-\left(f+{\frac {x}{r}}\right){\frac {dy\delta x-y\,d\,\delta x}{dt{\sqrt {2h(1+m)}}}}\\&-{\frac {y}{r}}{\frac {dy\delta y-y\,d\,\delta y}{dt{\sqrt {2h(1+m)}}}}+{\frac {2dy\delta h}{dt{\sqrt {2h(1+m)}}}},\\\\\delta g=&-{\frac {x}{r^{3}}}(x\delta y-y\delta x)+\left(f+{\frac {x}{r}}\right){\frac {dx\delta x-x\,d\,\delta x}{dt{\sqrt {2h(1+m)}}}}\\&+{\frac {y}{r}}{\frac {dx\delta y-x\,d\,\delta y}{dt{\sqrt {2h(1+m)}}}}-{\frac {2dx\delta h}{dt{\sqrt {2h(1+m)}}}},\\\\\delta b=&{\frac {dy\delta z-y\,d\,\delta z}{dt{\sqrt {2h(1+m)}}}}\\\\\delta c=&{\frac {-dx\delta z+x\,d\,\delta z}{dt{\sqrt {2h(1+m)}}}}\\\\\delta i=&3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}+{\frac {2(p\delta q-q\delta p)}{a^{2}(1-4n){\sqrt {n}}}}-{\frac {q(p-4an)\delta n}{a^{2}(1-4n)n^{\frac {3}{2}}}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b70bbb52f417c0d76b685f97b42fb99c34e6ec6f)
mais, à cause de
on aura
![{\displaystyle p=fx,\quad q=fy,\quad \delta p=f\delta x+x\delta f+y\delta g,\quad \delta q=f\delta y+y\delta f-x\delta g\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f50c12c689ccf01b51bd251358d767ffe52756a0)
donc
![{\displaystyle p\delta q-q\delta p=f^{2}(x\delta y-y\delta x)-\left(x^{2}+y^{2}\right)f\delta g\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0b8550b076fd517b19e31e3d31588f7cd6baf626)
de plus
![{\displaystyle n={\frac {h}{a}}={\frac {1-f^{2}}{4}},\quad \delta n=-{\frac {f\delta f}{2}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d11f96037983af11a63deaeb1d14bc0d9194de99)
donc la dernière formule deviendra
![{\displaystyle \delta i=3{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}+{\frac {2\left(x\delta y-y\delta x-{\cfrac {r^{2}\delta g}{f}}\right)}{\sqrt {a^{3}h}}}+{\frac {y(fx-4h)\delta f}{2{\sqrt {ah^{3}}}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4745996587fcfbc075b857384a393cec3f6b8659)
Et l’on remarquera que, à cause de
et de
on aura encore cette équation entre
et
savoir
![{\displaystyle a\delta h-h\delta a+{\frac {a^{2}f\delta f}{2}}=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aa0b8abb3b3b625ecd285b4d55081c74a245bd05)
laquelle pourra tenir lieu d’une quelconque des trois premières formules.
Telles sont donc les valeurs des quantités
en
par conséquent, si l’on met dans ces formules les valeurs de
et
en
et
savoir
![{\displaystyle x={\frac {a}{2}}(\cos u-f),\quad y={\frac {a}{2}}{\sqrt {1-f^{2}}}\cos u,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d789f731c6dcc9439ec02633850e164194282f09)
et (no 21)
![{\displaystyle dt={\sqrt {\frac {a}{2(1+m)}}}rdu={\frac {a}{2}}{\sqrt {\frac {a}{2(1+m)}}}(1-f\cos u)du,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c2dc1c07df8d680ab54751f96657727fef9a0611)
à cause de
elles doivent devenir identiques avec celles du no 27 ; de sorte qu’on pourra trouver, par la comparaison des coefficients de
dans les expressions de
les valeurs des quantités
des formules de ce numéro ; lesquelles valeurs seront nécessairement les mêmes que si on les avait déduites des formules du no 26, au moyen de l’élimination.
32. Revenons maintenant aux expressions de
trouvées dans le numéro que nous venons de citer ; comme ces expressions satisfont aux équations différentielies du no 23, les quantités
demeurant constantes et indéterminées, il s’ensuit que, par la substitution de ces valeurs de
dans les mêmes équations, tous les termes doivent se détruire d’eux-mêmes, indépendamment des quantités
donc, en général, les termes qui renfermeront les quantités finies
se détruiront toujours dans les équations dont il s’agit, soit que ces quantités soient constantes ou non.
Donc, si l’on suppose, ce qui est permis, que les mêmes expressions de
satisfassent aux équations du no 22 (lesquelles ne diffèrent, comme l’on voit, de celles du no 23 que par les termes
ajoutés à leurs seconds membres), mais en y regardant les six quantités
comme des variables indéterminées, et qu’on en fasse la substitution, il est visible que les termes qui renfermeront ces variables finies s’en iront aussi, et que les équations résultantes seront (no 26)
![{\displaystyle {\begin{aligned}&{\frac {\mathrm {A} d^{2}\delta a+\mathrm {B} d^{2}\delta f+\mathrm {C} d^{2}\delta g+\mathrm {D} d^{2}\delta i}{dt^{2}}}\\&\qquad \qquad \qquad \qquad +2{\frac {d\mathrm {A} \,d\,\delta a+d\mathrm {B} \,d\,\delta f+d\mathrm {C} \,d\,\delta g+d\mathrm {D} \,d\,\delta i}{dt^{2}}}=-\mu \mathrm {X} ,\\\\&{\frac {\mathrm {E} d^{2}\delta a+\mathrm {F} d^{2}\delta f+\mathrm {G} d^{2}\delta g+\mathrm {H} d^{2}\delta i}{dt^{2}}}\\&\qquad \qquad \qquad \qquad +2{\frac {d\mathrm {E} \,d\,\delta a+d\mathrm {F} \,d\,\delta f+d\mathrm {G} \,d\,\delta g+d\mathrm {H} \,d\,\delta i}{dt^{2}}}=-\mu \mathrm {Y} ,\\\\&\qquad \qquad \qquad \qquad \quad \ \ \ {\frac {\mathrm {K} d^{2}\delta b+\mathrm {L} d^{2}\delta c+2\left(d\mathrm {K} d\,\delta b+d\mathrm {L} d\,\delta c\right)}{dt^{2}}}=-\mu \mathrm {Z} .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e6d6d2a6aecc6fd5ed368370a54967d61912cd2)
33. Comme il y a ici six variables indéterminées, et qu’il n’y a que trois équations pour la détermination de ces variables, il est clair qu’on peut supposer a volonté trois autres équations entre ces mêmes variables, et il sera à propos de prendre ces équations en sorte que les différences secondes des variables
disparaissent d’elles-mèmes ; c’est de quoi on viendra à bout en supposant
![{\displaystyle {\begin{aligned}{\frac {\mathrm {A} \,d\,\delta a+\mathrm {B} \,d\,\delta f+\mathrm {C} \,d\,\delta g+\mathrm {D} \,d\,\delta i}{dt}}=&0,\\{\frac {\mathrm {E} \,d\,\delta a+\mathrm {F} \,d\,\delta f+\mathrm {G} \,d\,\delta g+\mathrm {H} \,d\,\delta i}{dt}}=&0,\\{\frac {\mathrm {K} d\,\delta b+\mathrm {L} d\,\delta c}{dt}}=0\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e704571ca189dcd17cd6f91b4434cfa50482de31)
car, en retranchant respectivement des équations précédentes les différences de celles-ci, on aura
![{\displaystyle {\begin{aligned}{\frac {d\mathrm {A} \,d\,\delta a+d\mathrm {B} \,d\,\delta f+d\mathrm {C} \,d\,\delta g+d\mathrm {D} \,d\,\delta i}{dt^{2}}}=&-\mu \mathrm {X} ,\\{\frac {d\mathrm {E} \,d\,\delta a+d\mathrm {F} \,d\,\delta f+d\mathrm {G} \,d\,\delta g+d\mathrm {H} \,d\,\delta i}{dt^{2}}}=&-\mu \mathrm {Y} ,\\{\frac {\left(d\mathrm {K} d\,\delta b+d\mathrm {L} d\,\delta c\right)}{dt^{2}}}=&-\mu \mathrm {Z} .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5b839eeb875896026bb565aa9da88e490041c953)
Ayant ainsi six équations entre les six quantités
![{\displaystyle {\frac {d\,\delta a}{dt}},\quad {\frac {d\,\delta f}{dt}},\quad {\frac {d\,\delta g}{dt}},\quad {\frac {d\,\delta i}{dt}},\quad {\frac {d\,\delta b}{dt}},\quad {\frac {d\,\delta c}{dt}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3e2499a755194900191fc42e95eb9afd3d71078b)
on déterminera, par l’élimination, la valeur de chacune de ces quantités ensuite il n’y aura plus qu’à multiplier ces différentes valeurs par
![{\displaystyle dt,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/993271528b71660c0c4282c0d0b0f4cbf88e5e74)
et les intégrer ; on aura de cette manière les valeurs des variables
![{\displaystyle \delta a,\delta b,\delta c,\delta f,\delta g,\delta i}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c377ef4e26fc92b96d366cdd2e51ca8f017e246f)
qu’il faudra substituer dans les expressions de
![{\displaystyle \delta x,\delta y,\delta z\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a006a33687288da44ecc7c1eec2e5cbdec682233)
et les équations du
no 22, qui expriment les perturbations de la comète, seront résolues.
34. Il est important de remarquer que, si l’on différentie les expressions de
on aura, en vertu des équations supposées ci-dessus,
![{\displaystyle {\begin{aligned}d\,\delta x=&d\mathrm {A} \delta a+d\mathrm {B} \delta f+d\mathrm {C} \delta g+d\mathrm {D} \delta i,\\d\,\delta y=&d\mathrm {E} \delta a+d\mathrm {F} \delta f\,+d\mathrm {G} \delta g+d\mathrm {H} \delta i,\\d\,\delta z=&d\mathrm {K} \delta b+d\mathrm {L} \delta c,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eeadd0f6383985a0c588bdc884b40d154ea28e86)
précisément comme si les quantités
étaient constantes, parce que les termes dépendant des variations de ces quantités sont précisément ceux qui forment les équations supposées. D’où il est facile de conclure que, si les équations différentielles du no 22 contenaient aussi les différences premières de
elles s’intégreraient également par la méthode du numéro précédent, et l’on parviendrait aux mêmes résultats.
Il y a plus, et c’est ici le point essentiel, dans l’orbite non altérée on a pour coordonnées
fonctions du temps
et des six constantes arbitraires
lesquelles déterminent les six éléments de l’orbite, savoir, le grand axe, l’excentricité, la position du périhélie, l’époque du passage par le périhélie, le lieu du nœud et l’inclinaison (nos 17, 19, 20). Dans l’orbite troublée, les coordonnées sont
les quantités
n’étant autre chose que les variations de
provenant des variations
des six constantes
comme on l’a vu ci-dessus. Ainsi, dans l’orbite troublée, les coordonnées sont exprimées de la même manière que dans l’orbite non troublée, c’est-à-dire, qu’elles sont les mêmes fonctions de
et de
qu’elles le sont de
dans l’orbite non troublée. Par conséquent on peut à chaque instant regarder l’orbite troublée comme étant de la même forme que l’orbite non troublée, mais dont les éléments dépendent des quantités
lesquelles étant variables, il s’ensuit que les éléments de l’orbite troublée seront variables aussi, et que les quantités
serviront à déterminer leurs variations. Or, comme nous venons de voir que les valeurs de ces quantités sont telles, que les différences premières de
sont les mêmes que si ces quantités étaient constantes, il est aisé d’en conclure que les éléments de l’orbite troublée, quoique essentiellement variables, peuvent néanmoins être regardés et traités comme constants pendant un instant, et qu’ainsi non-seulement le lieu de la comète dans l’orbite troublée, mais encore son mouvement instantané, c’est-à-dire, sa vitesse et sa direction, seront dans chaque instant les mêmes que l’on trouverait en les déterminant à l’ordinaire dans une orbite fixe dont les éléments seraient ceux qui répondent à ce même instant.
35. La difficulté est donc réduite à déterminer les valeurs des quantités
au moyen des six équations du no 33.
Or, en examinant ces équations et en les comparant avec les formules qui donnent les valeurs de
et de leurs différences
(nos 26, 27), il est aisé de voir qu’elles sont semblables, et qu’elles peuvent se déduire de ces mêmes formules en y changeant seulement les quantités
en leurs différences
et en y supposant en même temps
![{\displaystyle \delta x=0,\ \ \delta y=0,\ \ \delta z=0,\ \ {\frac {d\,\delta x}{dt}}=-\mu \mathrm {X} ,\ \ {\frac {d\,\delta y}{dt}}=-\mu \mathrm {Y} ,\ \ {\frac {d\,\delta z}{dt}}=-\mu \mathrm {Z} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/172d0577c3c96f960a8f9e47025480cb97bcc9d4)
Donc, en faisant ces mêmes suppositions dans les expressions de
, en
on aura les valeurs des différences ![{\displaystyle {\frac {d\,\delta a}{dt}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4b899157b5dc38b3334f0f02583c4775ca4b1f98)
![{\displaystyle {\frac {d\,\delta f}{dt}},\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/774f928c098484d24981e0bc2b3425e0cc5cf799)
ainsi l’on aura (
no 23)
![{\displaystyle {\begin{alignedat}{2}{\frac {d\,\delta a}{dt}}=&-\mu \mathrm {\left(C'X+D'Y\right)} ,\qquad &{\frac {d\,\delta f}{dt}}=&-\mu \mathrm {\left(G'X+H'Y\right)} ,\\{\frac {d\,\delta g}{dt}}=&-\mu \mathrm {\left(C''X+D''Y\right)} ,&{\frac {d\,\delta i}{dt}}=&-\mu \mathrm {\left(G''X+H''Y\right)} ,\\{\frac {d\,\delta b}{dt}}=&-\mu \mathrm {L'Z} ,&{\frac {d\,\delta c}{dt}}=&-\mu \mathrm {L''Z} .\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/390f2c6de04618b69cc8205ae3e78e641638435b)
36. En général, il est visible que les équations du no 33 ne sont autre chose que les différentielles de celles qui donnent les valeurs de
en
en y faisant varier seulement ces dernières quantités, ainsi que les différences premières
et mettant à la place des différences secondes
les quantités
de sorte qu’en faisant les mêmes opérations sur les équations qui donnent directement les valeurs de
en
on aura sur-le-champ les valeurs cherchées de
C’est ce qu’on peut aussi démontrer à priori par le raisonnement suivant.
Soit, en général,
![{\displaystyle \Delta =\Phi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/764ce568787106d1096e96cde3a9e32bcb252c85)
une quelconque des équations dont il s’agit, étant une des six constantes arbitraires
et
la fonction de
et de
qui lui est égale, il est clair que cette équation considérée en elle-même n’est autre chose qu’une intégrale première, ou du premier ordre, des équations du no 23, dans laquelle à est la constante arbitraire introduite par l’intégration ; donc, en différentiant, on aura cette équation du second ordre
![{\displaystyle d\Phi =0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8d569615904dc4905b58fb10f401e46f82b9a5d)
laquelle, ne contenant plus de constantes arbitraires, devra être identique, c’est-à-dire, avoir lieu en même temps que les équations du numéro cité ; de sorte que la différentielle
![{\displaystyle d\Phi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/eaeade320fd3cfaadbdf57739c3edf0e59a5ddc9)
devra être telle que, si l’on y substitue à la place des différences secondes
![{\displaystyle {\frac {d^{2}\delta x}{dt^{2}}},\ {\frac {d^{2}\delta y}{dt^{2}}},\ {\frac {d^{2}\delta z}{dt^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6047d70785ac05f65708eab304640656787ab6d5)
leurs valeurs données par ces mêmes équations, tous ses termes se détruisent d’eux-mêmes ; c’est aussi de quoi on pourra se convaincre
à posteriori par le calcul.
Or, comme les équations du no 22 ne diffèrent de celles du no 23 que parce que les valeurs de
ont les termes
de plus, il s’ensuit que si, au lieu de substituer dans l’expression de
les valeurs de
déduites des équations du no 23, on y substituait les valeurs de ces mêmes quantités, déduites des équations du no 22, on aurait nécessairement le même résultat que si l’on y substituait simplement
à la place de
et qu’on y effaçât en même temps tous les autres termes. Soit
ce que devient alors la valeur de
(
étant ici regardée comme variable) on aura donc, pour les équations du no 22,
![{\displaystyle d\Phi =d\Delta ,\quad \mathrm {et\ de\ l{\grave {a}}} \quad \Phi =\Delta ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d3189c81828615132c8cc1fa2f39d4cfee87ff6c)
comme pour celles du no 23, mais avec cette différence, que
ne sera plus ici constante, mais une fonction donnée de
et cette équation
sera par conséquent aussi une intégrale première des équations du no 22.
D’où il est aisé de conclure, en général, que, pour trouver les intégrales de ces dernières équations, qui sont proprement celles qui déterminent les perturbations de la comète, il n’y aura qu’à différentier chacune des formules
trouvées plus haut (nos 30, 31), en n’y faisant varier que la constante
et les différences premières
et y substituer ensuite à la place de
les quantités
on aura par ce moyen la valeur de
dont l’intégrale sera celle de
Ayant déterminé ainsi les valeurs des différentes quantités
qui étaient auparavant constantes, et qui sont devenues maintenant des fonctions de
on aura des intégrales premières de la même forme qu’auparavant ; par conséquent les intégrales secondes ou finies qui résulteront de celles-là par l’élimination des différences premières
seront encore de la même forme ; d’où il s’ensuit que tant ces différences que les variables finies
seront aussi de la même forme, c’est-à-dire, les mêmes fonctions de
et des différentes quantités
que dans le cas où ces quantités seraient constantes.
Et il est facile de se convaincre qu’il n’est pas nécessaire, pour l’exactitude de cette méthode, que les différentes constantes
soient dégagées tout à fait des variables dans les intégrales premières des équations du no 23, ainsi que nous l’avons supposé il suffit de les imaginer dégagées, ce qui est toujours possible, et de les traiter comme toutes variables à la fois dans la différentiation des mêmes équations intégrales ; on éliminera ensuite successivement les différentielles de ces différentes quantités
pour avoir la valeur de chacune de ces différentielles.
Voilà, comme l’on voit, un moyen aussi simple que direct pour déduire les intégrales des équations du no 22 de celles des équations plus simples du no 23, et, en général, pour intégrer toutes sortes d’équations linéaires, en supposant qu’on sache déjà intégrer ces mêmes équations dans le cas où elles ne contiendraient aucun terme tout connu.
37. Qu’on différentie donc, d’après la méthode précédente, les formules du no 31, en y faisant varier seulement les quantités
ainsi que les trois differences premières
et qu’on y mette ensuite, à la place des différences secondes
les quantités
c’est-à-dire,
à la place de
on aura les équations suivantes
![{\displaystyle {\begin{aligned}d\,\delta h=&-\mu {\frac {\sqrt {2h}}{1+m}}(x\mathrm {Y} -y\mathrm {X} )dt,\\d\,\delta a=&-{\frac {\mu }{1+m}}a^{2}(\mathrm {X} dx+\mathrm {Y} dy),\\d\,\delta f=&-{\frac {\mu ydt}{\sqrt {2h(1+m)}}}\left[\left(f+{\frac {x}{r}}\right)\mathrm {X} +{\frac {y}{r}}\mathrm {Y} \right]+{\frac {2dyd\,\delta h}{dt{\sqrt {2h(1+m)}}}},\\d\,\delta g=&\quad {\frac {\mu xdt}{\sqrt {2h(1+m)}}}\left[\left(f+{\frac {x}{r}}\right)\mathrm {X} +{\frac {y}{r}}\mathrm {Y} \right]-{\frac {2dxd\,\delta h}{dt{\sqrt {2h(1+m)}}}},\\d\,\delta b=&\quad {\frac {\mu y}{\sqrt {2h(1+m)}}}\mathrm {Z} dt,\\d\,\delta c=&-{\frac {\mu x}{\sqrt {2h(1+m)}}}\mathrm {Z} dt,\\d\,\delta i=&3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {d\,\delta a}{a}}-{\frac {2r^{2}d\,\delta g}{f{\sqrt {a^{3}h}}}}+{\frac {y(fx-4h)d\,\delta f}{2{\sqrt {ah^{3}}}}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a1fb4944a734aa200ef273be65e11bc61b074dbc)
et l’équation entre
étant différentiée aussi, donnera
![{\displaystyle ad\,\delta h-hd\,\delta a+{\frac {a^{2}fd\,\delta f}{2}}=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc9be984c7b53e2b5f144b39616febbf05e8d8a9)
qui servira à déterminer, si l’on veut,
en connaissant
et ![{\displaystyle d\,\delta a.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94259102e81ef6c536fc4809f82c877ee2e57042)
Or je remarque qu’on a cette combinaison
![{\displaystyle xd\,\delta f+yd\,\delta g=2{\frac {xdy-ydx}{dt{\sqrt {2h(1+m)}}}}d\,\delta h=2d\,\delta h\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d4784f7f3f696b3e53519c73b8ebf6f0408fb8a6)
de sorte qu’on aura
![{\displaystyle d\,\delta g={\frac {2d\,\delta h-xd\,\delta f}{y}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c48732c3d1df78ced7ae469861bd69cce502bcb7)
ainsi, comme
![{\displaystyle d\,\delta f={\frac {2(hd\,\delta a-ad\,\delta h)}{a^{2}f}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9cd7f98e0d0c18e947bbb14952c42da0a8a87d42)
on aura
![{\displaystyle d\,\delta g={\frac {2}{afy}}\left[(af+x)d\,\delta h-xh\right]{\frac {d\,\delta a}{a}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2b7013f5631099939a7116c12f950e62e773f2d3)
valeurs que l’on pourra employer à la place des précédentes.
Telles sont les formules par l’intégration desquelles il faudra déterminer les valeurs des quantités
et il est visible que ces intégrations ne demandent que de simples quadratures, puisque les quantités
et
sont censées données en
d’après les mouvements supposés connus de la comète dans l’orbite non altérée, et de la planète perturbatrice dans son orbite.
38. Connaissant ces différentes quantités, on aura les éléments de l’orbite troublée, au moyen desquels on pourra calculer, par les méthodes ordinaires, tant le lieu que la vitesse et la direction de la comète dans un instant quelconque, ainsi que nous l’avons démontré plus haut (no 34).
Pour cet effet, on se ressouviendra que
est le grand axe de l’orbite non altérée,
le paramètre du grand axe, et
l’excentricité (no 17).
Ainsi
sera le grand axe de l’orbite troublée,
le paramètre de cette orbite, et
son excentricité.
Ensuite, en différentiant suivant
les valeurs de
et de
de ce même no 17, et faisant, suivant l’hypothèse du no 25,
on aura
![{\displaystyle \delta b=-\sin \omega \delta \psi ,\quad \delta c=\cos \omega \delta \psi \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/683907d23f6e5059582e850a8053902213ab42af)
ainsi
sera l’inclinaison du plan de l’orbite troublée sur le plan de l’orbite non troublée, et
sera l’angle que la ligne des nœuds de ces deux plans fait avec l’axe des
lequel est en même temps le grand axe de l’orbite non altérée (no 25) ; de sorte que sera proprement la longitude du nœud ascendant de l’orbite troublée, comptée sur le plan de l’orbite non troublée depuis le périhélie de cette dernière orbite.
En différentiant de même les valeurs de
et de
du no 17, et faisant, d’après le no 25,
et
on aura
![{\displaystyle \delta f=\delta e,\quad \delta g=e\delta \varepsilon \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d720f6b6738211d456ee059e967dc26ce69317f7)
et il est clair, par les dénominations du no 13, que
sera la lon-
gitude du point de l’orbite troublée qui est à
![{\displaystyle 90}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8bda6ec8d2b197921d42d52bd57b265cf2c7b392)
degrés du périhélie, comptée sur le plan de l’orbite non troublée, depuis le périhélie de celle-ci ; mais, à cause que ces deux orbites ne font entre elles qu’un très-petit angle
![{\displaystyle \delta \psi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0858c07501362600bd7d707b8d212de1cf955b1a)
et que nous négligeons ici les
![{\displaystyle \delta \psi ^{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fcdc9d53e2d87b0d4f63b43b397cddfc1e5b1aa5)
il est très-facile de prouver que sera la longitude même du périhélie de l’orbite troublée, la projection d’un arc de
![{\displaystyle 90}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8bda6ec8d2b197921d42d52bd57b265cf2c7b392)
degrés ne pouvant différer de
![{\displaystyle 90}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8bda6ec8d2b197921d42d52bd57b265cf2c7b392)
degrés que par des quantités de l’ordre de
![{\displaystyle \delta \psi ^{2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0736089a8a1c66013d7d24f777ddac271bc38817)
Ainsi le petit angle
![{\displaystyle \delta \varepsilon }](https://wikimedia.org/api/rest_v1/media/math/render/svg/844848a786c7840e72bfe5ea141639e5606d37a1)
exprimera proprement le mouvement du périhélie en longitude, en vertu des perturbations.
Enfin on se rappellera que
est l’époque de l’anomalie moyenne dans l’orbite non troublée, c’est-à-dire, la valeur de cette anomalie lorsque
(no 20) ; donc
sera aussi l’époque de la même anomalie dans l’orbite troublée ; en sorte qu’ajoutant à cette époque le mouvement moyen pendant le temps
dans une orbite dont le grand axe serait
on aura l’anomalie moyenne qui servira à déterminer le lieu de la comète dans l’orbite troublée.
Ainsi
![{\displaystyle \theta =2t{\sqrt {\frac {2(1+m)}{a^{3}}}}+i}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c77b7781f530d85edab5a5381d7bcb8d7814a929)
étant (numéro cité) l’anomalie moyenne dans l’orbite non troublée, on aura
pour l’anomalie moyenne dans l’orbite troublée, et l’on trouvera la valeur de
par la différentiation de l’équation précédente, en y faisant varier
et
seulement, en sorte qu’on aura
![{\displaystyle \delta \theta =-3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}+\delta i.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7eb977b3f2150e7975a0e26864f0bdb1ac00399c)
Comme
et
sont ici des quantités variables, si l’on différentie à l’ordinaire cette valeur de
on aura
![{\displaystyle d\,\delta \theta =-3dt{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}-3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {d\,\delta a}{a}}+d\,\delta i\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c8f2f3d34932323b52242e5ad269fb04ed02e3c0)
et substituant pour
sa valeur trouvée dans le numéro précédent,
on aura
![{\displaystyle d\,\delta \theta =-3dt{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {\delta a}{a}}-{\frac {2r^{2}d\,\delta g}{f{\sqrt {a^{3}h}}}}+{\frac {y(xf-4h)d\,\delta f}{2{\sqrt {ah^{3}}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3deb9453c1d2ff9d2245c234c729a06c5a3f87fa)
dont l’intégrale donnera directement la valeur de
, qui est l’altération de l’anomalie moyenne causée par les perturbations.
39. Nous avons donné, dans la première Section (nos 10, 11), une manière de transformer les équations générales des perturbations, en sorte que les forces perturbatrices deviennent très-petites lorsque la comète est à une grande distance du Soleil ; comme cette transformation est d’une grande utilité pour le calcul des perturbations dans la partie supérieure de l’orbite, il faut voir maintenant comment elle peut s’appliquer aussi aux formules que nous venons de trouver.
La transformation dont il s’agit consiste en ce que, si l’on fait
![{\displaystyle {\begin{aligned}\delta x=&\mu \left({\frac {x}{3(1+m)}}+{\frac {\xi }{1+\mu }}\right)+\delta x',\\\delta y=&\mu \left({\frac {y}{3(1+m)}}+{\frac {\eta }{1+\mu }}\right)+\delta y',\\\delta z=&\mu \left({\frac {z}{3(1+m)}}+{\frac {\zeta }{1+\mu }}\right)+\delta z',\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19cb77857bfffc1a1d2d17d8e48529a67f558376)
et de plus
![{\displaystyle \mathrm {X} '={\frac {d\mathrm {\left({\cfrac {1}{S}}-{\cfrac {1}{R}}\right)} }{dx}},\quad \mathrm {Y} '={\frac {d\mathrm {\left({\cfrac {1}{S}}-{\cfrac {1}{R}}\right)} }{dy}},\quad \mathrm {Z} '={\frac {d\mathrm {\left({\cfrac {1}{S}}-{\cfrac {1}{R}}\right)} }{dz}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8311286e1fd84ca23a4cce219213d2a627ba27a0)
on aura entre
et
les mêmes équations qu’entre
et
c’est-à-dire, des équations de la mêmes forme que celles du no 22, en y marquant seulement les quantités
chacune d’un trait.
On peut donc appliquer à ces équations les mêmes raisonnements et les mêmes opérations que nous venons de faire dans cette Section sur les équations du no 22, et en tirer des conclusions semblables. Ainsi, si l’on dénote par
des quantités analogues aux quantités
on aura des formules semblables à celles des nos 31 et 37 ci-dessus, en y marquant d’un trait les quantités
Par les premières, on aura les valeurs de
en
et par les autres les valeurs de
en
Supposons maintenant qu’on substitue dans les formules du no 31 les valeurs précédentes de
il est aisé de voir (à cause que ces quantités n’entrent dans les mêmes formules que sous une forme linéaire) que les valeurs des quantités
, deviendront
![{\displaystyle {\begin{alignedat}{2}\delta h=&\mu \mathrm {H} &+&\delta h',\\\delta a=&\mu \mathrm {A} &+&\delta a',\\\delta f=&\mu \mathrm {F} &+&\delta f',\\\delta g=&\mu \mathrm {G} &+&\delta g',\\\delta b=&\mu \mathrm {B} &+&\delta b',\\\delta c=&\mu \mathrm {C} &+&\delta c',\\\delta i=&\mu \mathrm {I} &+&\delta i',\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/625ca570a0ccdc2ee16174cf8dbc2ec65b96ac5b)
en dénotant par
les valeurs de
, provenant de la simple substitution de
à la place de
de
à la place de
et de
à la place de ![{\displaystyle \delta z.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ae10d4362d6fec6d0d35a94049e7ec3e02055b7f)
De sorte qu’en faisant
![{\displaystyle z=0,\quad {\frac {dz}{dt}}=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94d7ecd1388f50b37bada41f21e26b4446f8b8dc)
et mettant à la place de
sa valeur
(nos 29, 31), on aura
![{\displaystyle {\begin{aligned}\mathrm {H} =&{\frac {4h}{3(1+m)}}+{\frac {2h}{1+\mu }}{\frac {\xi dy-\eta dx+xd\eta -yd\xi }{dt{\sqrt {2h(1+m)}}}},\\\mathrm {A} =&{\frac {a^{2}}{3(1+m)}}\left({\frac {1}{r}}+{\frac {dx^{2}+dy^{2}}{(1+m)dt^{2}}}\right)+{\frac {a^{2}}{1+\mu }}\left({\frac {x\xi +y\eta }{r^{3}}}+{\frac {dxd\xi +dyd\eta }{(1+m)dt^{2}}}\right),\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4615307801d144152703ca7734687d5aba097429)
![{\displaystyle {\begin{aligned}\mathrm {F} =&-{\frac {f+{\cfrac {x}{r}}}{3(1+m)}}+{\frac {{\cfrac {y}{r^{3}}}(x\eta -y\xi )}{1+\mu }}\\&-{\frac {\left(f+{\cfrac {x}{r}}\right)(\xi dy-yd\xi )+{\cfrac {y}{r}}(\eta dy-yd\eta )}{(1+\mu )dt{\sqrt {2h(1+m)}}}}+{\frac {2\mathrm {H} dy}{dt{\sqrt {2h(1+m)}}}},\\\\\mathrm {G} =&-{\frac {\cfrac {y}{r}}{3(1+m)}}-{\frac {{\cfrac {x}{r^{3}}}(x\eta -y\xi )}{1+\mu }}\\&+{\frac {\left(f+{\cfrac {x}{r}}\right)(\xi dx-xd\xi )+{\cfrac {y}{r}}(\eta dx-xd\eta )}{(1+\mu )dt{\sqrt {2h(1+m)}}}}-{\frac {2\mathrm {H} dx}{dt{\sqrt {2h(1+m)}}}},\\\\\mathrm {B} =&{\frac {\zeta dy-yd\zeta }{(1+\mu )dt{\sqrt {2h(1+m)}}}},\\\\\mathrm {C} =&{\frac {xd\zeta -\zeta dx}{(1+\mu )dt{\sqrt {2h(1+m)}}}},\\\\\mathrm {I} =&{\frac {3t}{a}}{\sqrt {\frac {2(1+m)}{a^{3}}}}\mathrm {A} +2{\frac {\left({\cfrac {x\eta -y\xi }{1+\mu }}-{\cfrac {r^{2}\mathrm {G} }{f}}\right)}{\sqrt {a^{3}h}}}+{\frac {y(fx-4h)\mathrm {F} }{2{\sqrt {a^{3}h}}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/934987264f9acb65a3b078da02452ca67bcad5a0)
De plus on aura, par les formules du no 37, en y marquant d’un trait les quantités
![{\displaystyle {\begin{aligned}d\,\delta h'=&-\mu {\frac {\sqrt {2h}}{1+\mu }}(x\mathrm {Y} '-y\mathrm {X} ')dt,\\d\,\delta a'=&-{\frac {\mu }{1+m}}a^{2}(\mathrm {X} 'dx+\mathrm {Y} 'dy),\\d\,\delta f'=&-{\frac {\mu ydt}{\sqrt {2h(1+m)}}}\left[\left(f+{\frac {x}{r}}\right)\mathrm {X} '+{\frac {y}{r}}\mathrm {Y} '\right]+{\frac {2dyd\,\delta h'}{dt{\sqrt {2h(1+m)}}}},\\d\,\delta g'=&\quad {\frac {\mu xdt}{\sqrt {2h(1+m)}}}\left[\left(f+{\frac {x}{r}}\right)\mathrm {X} '+{\frac {y}{r}}\mathrm {Y} '\right]-{\frac {2dxd\,\delta h'}{dt{\sqrt {2h(1+m)}}}},\\d\,\delta b'=&\quad {\frac {\mu y}{\sqrt {2h(1+m)}}}\mathrm {Z} 'dt,\\d\,\delta c'=&-{\frac {\mu x}{\sqrt {2h(1+m)}}}\mathrm {Z} 'dt,\\d\,\delta i'=&3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {d\,\delta a'}{a}}-{\frac {2r^{2}d\,\delta g'}{f{\sqrt {a^{3}h}}}}+{\frac {y(fx-4h)d\,\delta f'}{2{\sqrt {ah^{3}}}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/092b7b04d2c98f32e541b74f0b84ce3d0748348d)
Donc, si l’on différentie les valeurs de
données ci-dessus, et
qu’on y substitue ensuite les valeurs précédentes de
![{\displaystyle d\,\delta h',d\,\delta a',\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7d1b6607d826778a1a95e5613351384b2f744ec6)
il viendra
![{\displaystyle {\begin{aligned}d\,\delta h=&\mu d\mathrm {H} -\mu {\frac {\sqrt {2h}}{1+\mu }}(x\mathrm {Y} '-y\mathrm {X} ')dt,\\d\,\delta a=&\mu d\mathrm {A} -{\frac {\mu }{1+m}}a^{2}(\mathrm {X} 'dx+\mathrm {Y} 'dy),\\d\,\delta f=&\mu d\mathrm {F} -{\frac {\mu ydt}{\sqrt {2h(1+m)}}}\left[\left(f+{\frac {x}{r}}\right)\mathrm {X} '+{\frac {y}{r}}\mathrm {Y} '\right]+{\frac {2dyd\,\delta h'}{dt{\sqrt {2h(1+m)}}}},\\d\,\delta g=&\mu d\mathrm {G} +{\frac {\mu xdt}{\sqrt {2h(1+m)}}}\left[\left(f+{\frac {x}{r}}\right)\mathrm {X} '+{\frac {y}{r}}\mathrm {Y} '\right]-{\frac {2dxd\,\delta h'}{dt{\sqrt {2h(1+m)}}}},\\d\,\delta b=&\mu d\mathrm {B} +{\frac {\mu y}{\sqrt {2h(1+m)}}}\mathrm {Z} 'dt,\\d\,\delta c=&\mu d\mathrm {C} -{\frac {\mu x}{\sqrt {2h(1+m)}}}\mathrm {Z} 'dt,\\d\,\delta i=&\mu d\mathrm {I} +3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {d\,\delta a'}{a}}-{\frac {2r^{2}d\,\delta g'}{f{\sqrt {a^{3}h}}}}+{\frac {y(fx-4h)d\,\delta f'}{2{\sqrt {ah^{3}}}}},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b35cf31069962dafcb91f14051ef144dbc174d9c)
formules qu’on pourra employer à la place de celles du no 37, avec lesquelles elles sont identiques dans le fond.
40. En comparant les formules précédentes avec celles du no 37, il est aisé d’en tirer cette conclusion générale, qu’il est permis de changer dans ces dernières les quantités
en
ourvu qu’on ajoute en même temps aux valeurs de
les quantités
De là il s’ensuit que, soit, par exemple,
la valeur de
dans les formules du no 37, on aura, en intégrant,
![{\displaystyle \delta h=\int \mu \Pi dt,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/280d22b4c9e82aa119f804841639b7d9f7dffa13)
cette intégrale étant supposée commencer au point où
Supposons maintenant qu’à commencer d’un point donné de l’orbite on veuille employer les quantités
à la place des
et qu’on dénote
par
![{\displaystyle \delta h'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4d9e5f130959203286b79f3424aecfbdbe5b941e)
la valeur de
![{\displaystyle \delta h}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d71ab197adb11465d6f6af2d885b8b4b11313646)
pour ce point, c’est-à-dire, la valeur de l’intégrale
![{\displaystyle \int \mu \Pi dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b30190506e7ac7915ee56f69752df0d9466366d)
étendue jusqu’à ce point ; soit
![{\displaystyle \Pi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/18754ba38b5204cfb3921f909b25d1c42b642613)
ce que devient
![{\displaystyle \Pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/eed3e3db6cc2028a183af948212ed2551d25c954)
en y changeant
![{\displaystyle \mathrm {X,Y,Z} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e9852b5e32e3746828ae95b7c8664ea44d19655)
en
![{\displaystyle \mathrm {X',Y',Z'} :}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a28f97efd7a4d7ed37cbbf6cb38ccb6569ac2dda)
on aura, en général, par les formules du numéro précédent,
![{\displaystyle d\,\delta h=\mu d\mathrm {H} +\mu \Pi 'dt\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3f3edee76bb341b7c2bab32ed53a03dd3c4546a2)
donc, intégrant,
![{\displaystyle \delta h=\mu \mathrm {H} +\int \mu \Pi 'dt+\mathrm {const} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/60b7986387efc68a3ad9b4de048c2506ab34afe7)
Soit
la valeur de
dans le même point de l’orbite, et supposons que l’intégrale
commence aussi à ce point dans lequel on a supposé que finit l’intégrale
on aura donc dans ce point
![{\displaystyle \delta h'=\mu \mathrm {H} '+\mathrm {const} .,\quad {\text{donc}}\quad \mathrm {const} .=\delta h'-\mu \mathrm {H} '\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/934c21b4dbdbf1794c1e9f67c29d52fe4f8d931f)
donc on aura, en général,
![{\displaystyle \delta h=\mu \mathrm {H} -\mu \mathrm {H} '+\delta h'+\int \mu \Pi 'dt,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a0f7f1f7311cd2d09d561318d9e9982544533fb3)
savoir
![{\displaystyle \delta h=\mu \mathrm {H} -\mu \mathrm {H} '+\int \mu \Pi dt+\int \mu \Pi 'dt.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/76ed8be9bc593d168fc7a3f1eb399a2a0c842907)
Supposons ensuite que, dans un autre point quelconque de l’orbite, on veuille changer de nouveau les quantités
en
et soient dénotées par
et par
les valeurs de
et de
pour ce second point ; on aura donc dans ce point
![{\displaystyle \delta h''=\mu \mathrm {H} ''-\mu \mathrm {H} '+\int \mu \Pi dt+\int \mu \Pi 'dt,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2c9f01d50cd05c7aabd51a14281c8466cd2429e8)
l’intégrale
étant supposée étendue jusqu’à ce second point. Or, lorsqu’on emploie les quantités
on a, en général,
![{\displaystyle d\,\delta h=\mu \Pi dt,\quad {\text{donc}}\quad \delta h=\int \mu \Pi dt+\mathrm {const} .\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5be7e876ace23a0a164ec0cf397fcbc902715fcc)
supposons que l’intégrale
![{\displaystyle \int \mu \Pi dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b30190506e7ac7915ee56f69752df0d9466366d)
commence à ce second point dans lequel
![{\displaystyle \delta h}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d71ab197adb11465d6f6af2d885b8b4b11313646)
devient
![{\displaystyle \delta h'',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8bdb9007caace7813eaddac3c05d2b863fc50a72)
et l’on aura
![{\displaystyle \delta h''=\mathrm {const} .\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b15ef556659c31726e873974ed87a622ebba0c46)
donc, en général,
![{\displaystyle \delta h=\int \mu \Pi dt+\delta h'',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/30704439bd73e3e7203c6349f789681cf9a7785a)
et, substituant la valeur de ![{\displaystyle \delta h'',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8bdb9007caace7813eaddac3c05d2b863fc50a72)
![{\displaystyle \delta h=\mu \mathrm {H} ''-\mu \mathrm {H} '+\int \mu \Pi dt+\int \mu \Pi 'dt+\int \mu \Pi dt\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/633ca27da3b869cfd7904e6781a4febfe4cbfa29)
dans cette formule, la première intégrale
est supposée commencer au point de l’orbite où
est nul et s’étendre seulement jusqu’au point où les quantités
se changent en
la seconde intégrale
est supposée commencer à ce point et s’étendre jusqu’à l’autre point où les quantités
redeviennent
enfin la troisième intégrale
commence à ce dernier point et s’étend indéfiniment de sorte que ces différentes intégrales ne forment proprement qu’une seule intégrale, qui commence au point où
est nul et qui s’étend indéfiniment, mais avec cette condition que la quantité
se change en
dans une certaine étendue.
On voit par là que, dans l’intégration de la valeur de
du no 37, on peut changer à volonté les quantités
en leurs analogues
et rétablir ensuite celles-là à la place de celles-ci, pourvu qu’on ajoute en même temps à la valeur finie de
la quantité
qui est la différence des deux valeurs de
dont l’une
se rapporte au point où
se changent en
et dont l’autre
se rapporte au point où
redeviennent
On fera le même raisonnement sur chacune des autres formules du no 37, et l’on tirera des conclusions semblables. Ainsi, dans l’intégration de la valeur de
on pourra, pour un certain espace à volonté, changer
en
pourvu qu’on ajoute ensuite à la valeur finie de
l’excès de la valeur de
qui répond à la fin de cet espace sur la valeur de
qui répond au commencement du même espace, etc.
Et, si l’on voulait substituer à plusieurs reprises les quantités
à la place de
on ferait la même opération pour chaque nouvelle substitution.
41. Une des déterminations les plus importantes de la Théorie des perturbations des comètes est celle de l’altération du temps périodique. Rien n’est plus facile que de trouver cette altération par le moyen de la formule que nous avons donnée (no 38) pour l’anomalie moyenne dans l’orbite troublée. En effet,
exprimant, en général, l’anomalie moyenne dans l’orbite non altérée, et
l’anomalie moyenne qui a lieu en même temps dans l’orbite troublée, on aura, pour l’instant du périhélie dans l’orbite troublée,
d’où
ou (ce qui revient au même),
D’où l’on voit que, lorsque la comète passera au périhélie dans son orbite troublée, une comète fictive, qu’on supposerait se mouvoir dans l’orbite non altérée, serait encore éloignée de son périhélie de la quantité qui répond à l’anomalie moyenne
dans cette même orbite. Donc, comme on a, en général (no 20),
![{\displaystyle \theta =2t{\sqrt {\frac {2(1+m)}{a^{3}}}}+i,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/13f19a475d7cc8a171d97b39620145327eb2a1fd)
étant une constante dans l’orbite non altérée, si l’on dénote par
le temps qui répond à l’anomalie
dans cette orbite, on aura
![{\displaystyle \delta \theta =2\delta t{\sqrt {\frac {2(1+m)}{a^{3}}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f91851efae0315b40eb1402a6a058762a64a387c)
donc
![{\displaystyle \delta t={\sqrt {\frac {a^{3}}{8(1+m)}}}\delta \theta \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/628972f37bcae7bb8c1f5df582d5860ce906aeb0)
c’est le temps dont le passage au périhélie de l’orbite troublée précédera le passage au périhélie de l’orbite non altérée, ce temps étant exprimé par le mouvement moyen du Soleil qui y répond (numéro cité).
Dénotons par
et
les valeurs de
qui répondent à deux périhélies consécutifs, et par
les valeurs correspondantes de
en sorte que l’on ait
![{\displaystyle \delta t'={\sqrt {\frac {a^{3}}{8(1+m)}}}\delta \theta ',\quad \delta t''={\sqrt {\frac {a^{3}}{8(1+m)}}}\delta \theta ''\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f35f9ecfab72b3d1c4d82748376dcf378e9df98b)
soient de plus
et
les temps des passages par les deux périhélies consécutifs dans l’orbite non altérée on aura, pour les temps de ces passages dans l’orbite troublée,
donc la différence de ces temps, c’est-à-dire, l’intervalle de temps entre deux passages consécutifs au périhélie de l’orbite troublée, sera
où
est le même intervalle pour l’orbite non altérée. D’où il s’ensuit que la durée de la révolution anomalistique dans l’orbite troublée surpassera la même durée, dans l’orbite non altérée, du temps exprimé par
![{\displaystyle \delta t'-\delta t'',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/63371522510183a2be9a4e2c8ddbc86d514c4566)
ou par
![{\displaystyle {\sqrt {\frac {a^{3}}{8(1+m)}}}(\delta \theta '-\delta \theta '')\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d9604ef8e04ec7008ae405755306d16e80734c9c)
c’est l’altération produite par les perturbations.
Il faut remarquer que, pour avoir les valeurs de
et
il faudrait à la rigueur supposer, dans
mais, comme nous négligeons les carrés et les produits des forces perturbatrices, et par conséquent aussi de toutes les quantités résultant de ces forces, il suffira d’y faire
et
Nous venons de déterminer l’altération de la révolution anomalistique de la comète ; si l’on voulait avoir l’altération de sa révolution périodique, il faudrait défalquer de l’altération précédente le temps dû au changement du périhélie. Or nous avons vu (no 38) que le périhélie de l’orbite troublée est plus avancé que celui de l’orbite non altérée de l’angle
donc, si l’on dénote par
et
les valeurs de
qui répondent à
et
on aura
pour l’angle dont le périhélie de l’orbite troublée aura avancé pendant une révolution ; ainsi la quantité à défalquer de l’altération de la révolution anomalistique, pour avoir celle de la révolution périodique, sera le temps qui répond à l’angle ou à l’anomalie vraie
Pour trouver ce temps, on pourra employer la formule différentielle (no 21)
![{\displaystyle dt={\frac {r^{2}d\varphi }{\sqrt {2h(1+m)}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8d166480bb0bc4e2fd6d7580e8275a1638072b5f)
en faisant
et
égal à la distance périhélie dans l’orbite non altérée, laquelle est
![{\displaystyle {\frac {a-ae}{2}}={\frac {a}{2}}(1-e)={\frac {a\left(1-e^{2}\right)}{2(1+e)}}={\frac {2h}{1+e}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dca38057b200afe692f28e4526a281ebdb179faf)
de sorte qu’on aura, pour le temps cherché, la quantité
![{\displaystyle \left({\frac {2h}{1+e}}\right)^{2}{\frac {\delta \varepsilon ''-\delta \varepsilon '}{\sqrt {2h(1+m)}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/86097efb3f88776dfd1428f9afbb475d5927ee3a)
Donc la durée de la révolution périodique de la comète dans l’orbite troublée, c’est-à-dire, le temps qu’elle mettra à faire une révolution entière depuis son départ du périhélie jusqu’à ce qu’elle revienne sur la ligne du même périhélie, surpassera le temps de la révolution entière, dans l’orbite non altérée, de la quantité
![{\displaystyle {\sqrt {\frac {a^{3}}{8(1+m)}}}(\delta \theta '-\delta \theta '')-\left({\frac {2h}{1+e}}\right)^{2}{\frac {\delta \varepsilon ''-\delta \varepsilon '}{\sqrt {2h(1+m)}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/99e3d12cd8f0be265b9cb0441aa6326991826831)
laquelle, en substituant pour
et
leurs valeurs déduites de la formule du no 38, et dénotant par
et par
les valeurs de
qui répondent à
et
se réduit à celle-ci
![{\displaystyle {\frac {3(t''\delta a''-t'\delta a')}{2a}}-{\sqrt {\frac {a^{3}}{8(1+m)}}}(\delta i''-\delta i')-\left({\frac {2h}{1+e}}\right)^{2}{\frac {\delta \varepsilon ''-\delta \varepsilon '}{\sqrt {2h(1+m)}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/88ce2aa0a0b4e5608a8a0d080c78e6787fd213b1)
section quatrième.
application des théories précédentes au calcul des perturbations des comètes, et en particulier au calcul des perturbations de la comète de 1532 et de 1661.
42. Cette application se présente d’elle-même ; il ne s’agit que de trouver les valeurs des quantités
par l’intégration des formules du no 37, et l’on aura immédiatement les altérations des éléments de l’orbite de là comète dues aux perturbations (no 38) ; mais la grande difficulté consiste dans ces intégrations, lesquelles, à cause de la grande excentricité de l’orbite des comètes, ne peuvent s’exécuter, en général, par aucune méthode connue et demandent nécessairement des quadratures de courbes mécaniques.
Nous allons proposer les moyens qui nous paraissent les plus propres pour arriver à ce but.
Je commence par substituer, dans l’es équations du no 37, les valeurs de
(no 22), lesquelles, en effectuant les différentiations indiquées, deviennent
![{\displaystyle \mathrm {X} ={\frac {\xi }{\rho ^{3}}}+{\frac {x-\xi }{\mathrm {R} ^{3}}},\quad \mathrm {Y} ={\frac {\eta }{\rho ^{3}}}+{\frac {y-\eta }{\mathrm {R} ^{3}}},\quad \mathrm {Z} ={\frac {\zeta }{\rho ^{3}}}+{\frac {z-\zeta }{\mathrm {R} ^{3}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e6c30e9fc8fa1a6d2552ce2e614451dd596cbd24)
je substitue de plus, à la place des quantités
leurs valeurs exprimées par l’anomalie excentrique
parce que l’emploi de cette anomalie rend tout à la fois les formules plus simples et plus faciles à calculer ; ces valeurs sont (en faisant
par l’hypothèse du no 25)
![{\displaystyle x={\frac {a}{2}}(\cos u-f),\quad y={\sqrt {ah}}\sin u,\quad z=0\quad }](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e1aef60676ba17033bae58a9fdee5d090985fd8)
(
no 26),
![{\displaystyle r={\frac {a}{2}}(1-f\cos u),\quad dt={\sqrt {\frac {a}{2(1+m)}}}rdu\quad }](https://wikimedia.org/api/rest_v1/media/math/render/svg/b37681807e6286c15575ff7b67275b848fbf601b)
(n
os 20, 21).
Ces substitutions faites, si l’on suppose, pour plus de simplicité,
![{\displaystyle w\mathrm {P} ={\frac {1}{\rho ^{3}}}-{\frac {1}{\mathrm {R} ^{3}}},\quad \varpi ={\frac {1}{\mathrm {R} ^{3}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f7a1b3791b3cc1cc1e3c880e5e19902f8582dcec)
on aura des équations de la forme suivante
![{\displaystyle {\begin{aligned}d\,\delta h=&{\frac {\mu }{1+m}}\left[(\mathrm {H} )\Pi +(h)\varpi \right]du,\\d\,\delta a=&{\frac {\mu }{1+m}}\left[(\mathrm {A} )\Pi +(a)\varpi \right]du,\\d\,\delta f=&{\frac {\mu }{1+m}}\left[(\mathrm {F} )\Pi +(f)\varpi \right]du,\\d\,\delta g=&{\frac {\mu }{1+m}}\left[(\mathrm {G} )\Pi +(g)\varpi \right]du,\\d\,\delta b=&{\frac {\mu }{1+m}}\left[(\mathrm {B} )\Pi +(b)\varpi \right]du,\\d\,\delta c=&{\frac {\mu }{1+m}}\left[(\mathrm {C} )\Pi +(c)\varpi \right]du,\\d\,\delta i=&{\frac {\mu }{1+m}}\left[(\mathrm {I} )\ \,\Pi +(i)\varpi \right]du,\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/456c688a373324a90b34e08ac0d6dd1d8fbde5bd)
dans lesquelles on aura les valeurs suivantes des quantités ![{\displaystyle (\mathrm {H} ),(h),(\mathrm {A} ),(a),\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/31cd66e17333a606cab9dc3e739460f3a8219575)
![{\displaystyle {\begin{aligned}(\mathrm {H} )=&{\sqrt {ah}}(y\xi -x\eta )r,\quad (h)=0,\\(\mathrm {A} )=&a^{2}\left({\frac {a}{2}}\xi \sin u-{\sqrt {ah}}\,\eta \cos u\right),\quad (a)=-{\frac {a^{2}f}{2}}r\sin u,\\(\mathrm {F} )=&-{\sqrt {\frac {a}{4h}}}\left[(fr+x)\xi +y\eta \right]y+2{\sqrt {ah}}(y\xi -x\eta )\cos u,\ \ (f)=-{\sqrt {ah}}\,ry,\\(\mathrm {G} )=&{\sqrt {\frac {a}{4h}}}\left[(fr+x)\xi +y\eta \right]x+a(y\xi -x\eta )\sin u,\quad (g)={\sqrt {ah}}rx,\\(\mathrm {B} )=&{\sqrt {\frac {a}{4h}}}ry\zeta ,\quad (b)=0,\\(\mathrm {C} )=&-{\sqrt {\frac {a}{4h}}}rx\zeta ,\quad (c)=0,\\(\mathrm {I} )=&3t{\sqrt {\frac {2(1+m)}{a^{5}}}}(\mathrm {A} )-{\frac {2r^{2}(\mathrm {G} )}{f{\sqrt {a^{3}h}}}}+{\frac {y(fx-4h)}{2{\sqrt {ah^{3}}}}}(\mathrm {F} ),\\(i)=&3t{\sqrt {\frac {2(1+m)}{a^{5}}}}(a)\ \,-{\frac {2r^{2}(g)}{f{\sqrt {a^{3}h}}}}+{\frac {y(fx-4h)}{2{\sqrt {ah^{3}}}}}(f).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6d2080b8fefaa79c9049672b5fcbb5d0a747dd36)
Dans ces expressions, j’ai conservé, pour plus de simplicité, les lettres
![{\displaystyle x,y,r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aaf32d0fbb52e6eaae8d394fda047120f8464cde)
à la place de leurs valeurs en
![{\displaystyle \sin u}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a12840715005f2cf490e983fe03dc56e6aa26df)
et
![{\displaystyle \cos u\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3966d243b2bf2d06b1f1fcce9aeed25176a484db)
il est facile de les y substituer si on le juge à propos.
43. Il est visible, par les formules précédentes, que les quantités
sont toutes exprimées par des fonctions rationnelles et entières de
de sorte que, si l’on pouvait exprimer de même les quantitéss
et
par des fonctions rationnelles et entières de
et
l’intégration des équations différentielles dont il s’agit n’aurait aucune difficulté. Voyons quels sont les obstacles qui s’opposent à cette réduction dans la Théorie des comètes.
On se rappellera d’abord que les quantités
sont les trois coordonnées rectangles du lieu de la planète perturbatrice, dont la masse est
que
est son rayon vecteur, et
la distance rectiligne entre le lieu de la planète et le lieu de la comète dans l’orbite non altérée (nos 2, 7) ; on se rappellera ensuite que nous prenons pour le plan de projection celui de l’orbite non altérée de la comète, et pour l’axe des abscisses la ligne du périhélie de cette orbite (no 25).
Nommons
l’inclinaison du plan de l’orbite de la planète sur le plan de l’orbite non altérée de la comète, et
la longitude du nœud ascendant de l’orbite de la planète, comptée sur le plan de l’orbite de la comète depuis le périhélie de cette orbite.
Soit, de plus,
l’argument de latitude de la planète, c’est-à-dire, la longitude dans son orbite, moins la longitude de son nœud avec l’orbite de la comète.
Il est facile de comprendre que l’on aura pour
des expressions semblables à celles de
du no 19, en y changeant
en
en
en
et
en
on aura donc ainsi
![{\displaystyle {\begin{aligned}\xi =&\rho (\cos \Omega \cos \lambda -\sin \Omega \cos \Psi \sin \lambda ),\\\eta =&\rho (\sin \Omega \cos \lambda +\cos \Omega \cos \Psi \sin \lambda ),\\\zeta =&\rho \sin \Psi \sin \lambda .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/edc13e95670afec939dbffe329e757103a96cd70)
Or on sait que dans les orbites des planètes, à cause de la petitesse de leur excentricité, on peut exprimer tant l’équation du centre que le rayon vecteur par des suites très-convergentes, qui procèdent suivant les sinus et cosinus de l’anomalie moyenne et de ses multiples (on trouve ces suites développées d’après les principales Tables astronomiques, dans le premier volume du Recueil des Tables, publié par l’Académie de Berlin) ; on pourra donc représenter par de semblables séries les valeurs de
et de
pour chaque planète, et il n’y aura plus qu’à exprimer l’anomalie moyenne de la planète par l’anomalie excentrique
de la comète.
Pour faire cette réduction, soient
le grand axe de l’orbite de la planète, et
son anomalie moyenne comptée à l’ordinaire depuis l’aphélie soit, de plus,
la valeur de l’anomalie moyenne
de la comète pour l’instant du passage de la planète par l’aphélie ; il est visible que
et
seront les anomalies contemporaines de la planète et de la comète, lesquelles doivent être entre elles en raison réciproque de la durée de leurs révolutions, et par conséquent, par les Théorèmes connus, en raison de
d’où il suit qu’on aura
![{\displaystyle \mathrm {M} =(\theta -\mathrm {T} ){\sqrt {\frac {a^{3}}{\alpha ^{3}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ef32f2417163d3909083bc7b8097478432124dd1)
où il n’y aura plus qu’à substituer pour
sa valeur
(no 20).
Comme, dans l’orbite des comètes, l’excentricité e est peu différente de l’unité, il est clair que les sinus et cosinus de
et de ses multiples ne sauraient s’exprimer par de simples sinus et cosinus de
et de ses multiples par conséquent, il est impossible d’exprimer, en général,
et
par des fonctions rationnelles et entières de
et de
C’est la première difficulté qui s’oppose à l’intégration des équations du numéro précédent.
La seconde difficulté vient du dénominateur irrationnel
en effet il est d’abord impossible, par la raison précédente, de réduire l’expression rationnelle de
laquelle est (no 2),
étant
![{\displaystyle \mathrm {R} ^{2}=r^{2}-2(x\xi +y\eta )+\rho ^{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a031d1275d306522fcd83b8b9b09b4dddbabf81f)
à une fonction rationnelle de
et
à plus forte raison le sera-t-il d’y réduire la quantité irrationnelle et rompue ![{\displaystyle {\frac {1}{\mathrm {R} ^{3}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/651f3ef9773bc955093cfc3e6bf3c2b3f58bfce7)
44. On est donc forcé, dans la Théorie des comètes, de renoncer à l’avantage de parvenir à des formules analytiques qui expriment les inégalités de leur mouvement pour un temps quelconque, telles que celles que l’on trouve pour les inégalités des planètes, et la seule ressource qui reste est de déterminer ces inégalités par parties, en partageant l’orbite de la comète en différentes portions, et calculant séparément l’effet des perturbations pour chacune de ces portions.
En effet, tant que l’angle
ne sera pas trop grand, on pourra exprimer son sinus et son cosinus par les séries connues qui procèdent suivant les puissances de l’arc, et par là on remédiera au premier inconvénient.
Ensuite on observera que, tant que le rayon
de la comète sera beaucoup moindre que le rayon
de la planète perturbatrice, et que, par conséquent,
et
seront moindres que
on pourra réduire la quantité
en une série convergente, en prenant
pour le premier terme.
De cette manière, on pourra donc intégrer les valeurs de
du no 42, depuis le périhélie de l’orbite de la comète jusqu’à un point de cette orbite dans lequel
et
soient des quantités encore assez petites.
Soit maintenant
l’anomalie excentrique qui répond à ce point ; on fera, en général,
et, tant que l’angle
sera assez petit, on pourra mettre les quantités à intégrer sous la forme rationnelle
![{\displaystyle \mathrm {\left(L+M\upsilon +N\upsilon ^{2}+\ldots \right)} d\upsilon \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f810ef70220c9b1deaa7ed0759413db938d88d0)
on intégrera donc derechef depuis
![{\displaystyle u=u'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d538ee4c6a02dbab7952a3e8985fad1dc74eff07)
jusqu’à
![{\displaystyle u=u'',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d83ba8ea8ff6b1560835f3e6135df4f4f1975d2f)
en supposant l’arc
![{\displaystyle u''-u'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f96fdbd6511e471a44d1de136592c4eefca10f96)
assez petit, et ainsi de suite.
45. On peut faciliter beaucoup ce calcul par la méthode connue des courbes paraboliques ; mais, pour pouvoir employer cette méthode en toute sûreté, il faut que les quantités qu’on veut exprimer par des formules paraboliques ne souffrent pas de trop grandes ni de trop fréquentes irrégularités ; autrement il arriverait que, parmi les coefficients de la série parabolique, il s’en trouverait de très-grands, ce qui diminuerait la convergence de la série et obligerait à la pousser à un grand nombre de termes. Il est donc nécessaire d’examiner à priori la nature des quantités auxquelles on veut appliquer la méthode des courbes paraboliques.
De ce que nous avons dit dans le no 43, il s’ensuit que les différentes quantités
ainsi que les quantités
et
peuvent être exprimées par des fonctions rationnelles et entières de sinus et de cosinus des angles
et
c’est-à-dire, de l’anomalie excentrique de la comète et du mouvement moyen correspondant de la planète ; donc, si l’on suppose que ces deux angles varient en même temps des angles contemporains
et
chacune des quantités dont il s’agit pourra être représentée, pendant ces variations, par une formule algébrique de la forme
![{\displaystyle \mathrm {L+M\beta +N\gamma +O\beta ^{2}+P\beta \gamma +Q\gamma ^{2}} +\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c7c1dfcb866da84b77aa88016c6c4e65ddb76266)
dans laquelle les quantités
seront toutes aussi des fonctions rationnelles et entières de
Or
donc, faisant croître
de
et
de
on aura
![{\displaystyle \gamma ={\sqrt {\frac {a^{3}}{\alpha ^{3}}}}\left[(1-e\cos u)\beta +{\frac {e\sin u}{2}}\beta ^{2}+\ldots \right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8c33dcbeee80b24d07f8f46dd0994d7e7cc6abcf)
Si donc on substitue cette valeur de
dans la formule précédente, elle
prendra cette forme plus simple
![{\displaystyle \mathrm {L+M\beta +N\beta ^{2}} +\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe5fa6ff720a96ab5d925db07ac4c8b6faccdbb1)
dans laquelle les quantités
seront pareillement des fonctions toutes rationnelles et entières de
en sorte que ces quantités ne pourront jamais augmenter au delà d’un certain terme. Et il est clair que la formule précédente, n’étant poussée que jusqu’au second degré, sera exacte, aux quantités près des ordres de
et de ![{\displaystyle \gamma ^{3}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b5238c7bb8f079834834e085b290ae3a1256ffb3)
Il semble qu’il faudrait faire une exception à l’égard des quantités
et
qui contiennent des termes multipliés par
et qui par conséquent ne sont pas uniquement des fonctions de sinus et cosinus de
et de
mais renferment aussi l’angle même
mais il est facile de se convaincre que cette circonstance ne peut apporter aucun changement à la conclusion précédente.
Si donc on dénote, en général, par
une quelconque des quantités dont il s’agit, et que
soient les valeurs de
qui répondent à
![{\displaystyle u=u_{0},\quad u=u_{1}=u_{0}+\beta ,\quad u=u_{2}=u_{1}+\beta ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/394f7d457585daf3e1ebe5975388df87490a8478)
il résulte de ce que nous venons de démontrer que, pour
(
étant un nombre quelconque compris entre zéro et
), on aura, aux quantités près des ordres de
et de ![{\displaystyle \gamma ^{3},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/078b55c4c5ecf64736277e8fcf708b769ca108cb)
![{\displaystyle \mathrm {V} =\mathrm {V} _{1}+\mathrm {V} '_{1}n+\mathrm {V} ''_{1}n^{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47e745441bc108e67c7a4b6a2eaa432647e76155)
formule qui pourra servir aussi par la même raison, en faisant
négatif depuis zéro jusqu’à ![{\displaystyle -1.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2c1ae9e73ea72a95921a7fbeba221311687f1367)
Or, comme
lorsque
et
lorsque
on aura
![{\displaystyle \mathrm {V_{0}=V_{1}-V'_{1}+V''_{1},\quad V_{2}=V_{1}+V'_{1}+V''_{1}} \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c659cc016a93576c3a1a168307a570fc8e7bdab4)
d’où l’on tire
![{\displaystyle \mathrm {V'_{1}={\frac {V_{2}-V_{0}}{2}},\quad V''_{1}={\frac {V_{2}-2V_{1}+V_{0}}{2}}} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad93a9f06e4cb8bd2822d130419cbe514840e397)
46. Cela posé, séparons, dans les équations différentielles du no 42, les termes divisés par
des autres, et représentons, en général, chacune de ces équations par
![{\displaystyle d\Delta ={\frac {\mu }{1+m}}\left(\mathrm {V} +{\frac {\mathrm {U} }{R^{\frac {3}{2}}}}\right)du,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/518e93901624b8b55358a6323bd61e73335f0fe1)
étant égal à
étant une des quantités
étant respectivement et
et
étant ![{\displaystyle (h)-(\mathrm {H} ),\ (a)-(\mathrm {A} ),\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a700f5f2ac6e484c2186904ace5b755d205c1b57)
Qu’on calcule les valeurs des quantités
et
pour trois anomalies excentriques
dont la commune différence soit
et qu’on marque ces valeurs respectivement
qu’on en déduise ensuite, par les dernières formules du numéro précédent, les valeurs de
ainsi que celles de
et qu’on substitue partout dans l’équation précédente
à la place de
on aura donc, en regardant maintenant
comme variable, la transformée
![{\displaystyle d\Delta ={\frac {\mu \beta }{1+m}}\left(\mathrm {V} _{1}+\mathrm {V} '_{1}n+\mathrm {V} ''_{1}n^{2}+{\frac {\mathrm {U} _{1}+\mathrm {U} '_{1}n+\mathrm {U} ''_{1}n^{2}}{\left(R_{1}+R'_{1}n+R''_{1}n^{2}\right)^{\frac {3}{2}}}}\right)dn,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/180d046e31efc6c97d6e8f3f6c226eaf6bd1373e)
qui, étant intégrée depuis
jusqu’à
donnera, aux quantités près de l’ordre de
et de
la valeur de
ou plutôt l’accroissement de
depuis l’anomalie excentrique
jusqu’à l’anomalie
en sorte que, désignant par
et
les valeurs de
qui répondent à ces deux anomalies, on aura
égale à l’intégrale du second membre de cette équation, prise depuis
jusqu’à ![{\displaystyle n=1.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/04d74ade48a04cf5d7a4d8a0f0a94a0bf6050973)
L’intégration de la partie
![{\displaystyle \left(\mathrm {V} _{1}+\mathrm {V} '_{1}n+\mathrm {V} ''_{1}n^{2}\right)dn}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de88502a9c44189af27861bd90c09a778d47b4db)
n’a aucune difficulté, et l’on trouve sur-le-champ pour l’intégrale totale
![{\displaystyle \mathrm {2V_{1}+{\frac {2}{3}}V''_{1}} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e3cd84a9ceb1807f998ae6b0734782f6c4ea630)
À l’égard de l’autre partie
![{\displaystyle {\frac {\mathrm {U} _{1}+\mathrm {U} '_{1}n+\mathrm {U} ''_{1}n^{2}}{\left(R_{1}+R'_{1}n+R''_{1}n^{2}\right)^{\frac {3}{2}}}}dn,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95666ab559a6f0a76849b709d08efbc2c2d199e5)
elle dépend de la quadrature de l’hyperbole ou du cercle, suivant que
![{\displaystyle R''_{1},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/869116a65058d8785592746c5b5efcce21f2d6d1)
est une quantité positive ou négative.
Pour en trouver l’intégrale, on supposera cette différentielle égale à
![{\displaystyle d{\frac {\mathrm {K+L} n}{\sqrt {R_{1}+R'_{1}n+R''_{1}n^{2}}}}+{\frac {\mathrm {M} dn}{\sqrt {R_{1}+R'_{1}n+R''_{1}n^{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a46d7d2c1f7a4543c0006417829122bfac46026)
et l’on trouvera par la comparaison des termes, après avoir réduit au même dénominateur,
![{\displaystyle {\begin{aligned}\mathrm {K} =&{\frac {{\frac {1}{2}}\mathrm {U} _{1}R'_{1}-\mathrm {U} '_{1}R_{1}+{\cfrac {\mathrm {U} ''_{1}R_{1}R'_{1}}{2R''}}}{R_{1}R''_{1}-{\frac {1}{4}}R_{1}^{'2}}},\\\mathrm {L} =&{\frac {\mathrm {U} _{1}R''_{1}-{\frac {1}{2}}\mathrm {U} '_{1}R'_{1}-\mathrm {U} ''_{1}\left(R_{1}-{\cfrac {R_{1}^{'2}}{2R''_{1}}}\right)}{R_{1}R''_{1}-{\frac {1}{4}}R_{1}^{'2}}},\\\mathrm {M} =&{\frac {\mathrm {U} ''_{1}}{R''_{1}}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aef5c12e9533a8d49a7a6c1c0fda5449c5c84eaa)
or l’intégrale de la première partie est évidemment
![{\displaystyle {\frac {\mathrm {K+L} n}{\sqrt {R_{1}+R'_{1}n+R''_{1}n^{2}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6bd80445f67bf13f45e06f61e65f83b3933ee109)
et celle de la seconde est, en faisant, pour abréger, ![{\displaystyle {\frac {\sqrt {R_{1}+R'_{1}n+R''_{1}n^{2}}}{{\frac {1}{2}}R'_{1}+R''_{1}n}}=\mathrm {N} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc3d118825ecb9fa07f6a924461e7aec0fd29df0)
![{\displaystyle {\frac {\mathrm {M} }{2{\sqrt {R''_{1}}}}}\log {\frac {1+\mathrm {N} {\sqrt {R''_{1}}}}{1-\mathrm {N} {\sqrt {R''_{1}}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53c951d7e373cf1514c4cd95b4effd4df4ff5497)
si
est positif ; mais si
est négatif, cette intégrale devient
![{\displaystyle {\frac {\mathrm {M} }{2{\sqrt {-R''_{1}}}}}\operatorname {arc} \,\operatorname {tang} \mathrm {N} {\sqrt {-R''_{1}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ba91ef1fdd7aefd806abfc067d36aaad6eca9911)
On fera maintenant dans ces formules
et
et l’on retranchera la seconde valeur de la première pour avoir l’intégrale complète ; or, en faisant
la quantité sous le signe
devient
et, en faisant
elle devient
Donc la valeur complète de l’intégrale de la différentielle dont il s’agit sera représentée par
![{\displaystyle {\frac {\mathrm {K+L} }{\sqrt {R_{2}}}}-{\frac {\mathrm {K-L} }{\sqrt {R_{0}}}}+{\frac {\mathrm {M} }{2{\sqrt {\pm R''_{1}}}}}\mathrm {P} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c39e27ff1d8b73d6319d56a5f0c85d64bb032bee)
en faisant
![{\displaystyle \mathrm {P} =\log {\frac {{\frac {1}{2}}R'_{1}+R''_{1}+{\sqrt {R_{2}R''_{1}}}}{{\frac {1}{2}}R'_{1}+R''_{1}-{\sqrt {R_{2}R''_{1}}}}}-\log {\frac {{\frac {1}{2}}R'_{1}+R''_{1}+{\sqrt {R_{0}R''_{1}}}}{{\frac {1}{2}}R'_{1}+R''_{1}-{\sqrt {R_{0}R''_{1}}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0c1f3e9f6d37756acaac377a7256b50c6ca96929)
si
est positif, ou bien
![{\displaystyle \mathrm {P} =\operatorname {arc} \,\operatorname {tang} {\frac {\sqrt {-R_{2}R''_{1}}}{{\frac {1}{2}}R_{1}+R''_{1}}}-\operatorname {arc} \,\operatorname {tang} {\frac {\sqrt {-R_{0}R''_{1}}}{{\frac {1}{2}}R_{1}+R''_{1}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a3e99271f8c160e81aedcaa42d312b363251ea7d)
si
est négatif.
Donc enfin on aura, aux quantités près des ordres de
et
![{\displaystyle \Delta _{2}-\Delta _{0}={\frac {\mu \beta }{1+m}}\left(2V'_{1}+{\frac {2}{3}}V''_{1}+{\frac {\mathrm {K+L} }{\sqrt {R_{2}}}}-{\frac {\mathrm {K-L} }{\sqrt {R_{0}}}}+{\frac {\mathrm {MP} }{2{\sqrt {\pm R''_{1}}}}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d5c1b9dc4057653e4ce70ce707d7453206a3a95)
47. Il n’y a que deux cas où la formule précédente ne puisse pas servir : l’un est celui de
et l’autre celui de
Soit : 1o
on aura à intégrer cette différentielle
![{\displaystyle {\frac {\mathrm {U} _{1}+\mathrm {U} '_{1}n+\mathrm {U} ''_{1}n^{2}}{(R_{1}+R'_{1}n)^{\frac {3}{2}}}}dn,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3964daffdf0f305acf7a70be136db229ed60aa5d)
et, supposant son intégrale de la forme
![{\displaystyle {\frac {\mathrm {K} +\mathrm {L} n+\mathrm {M} n^{2}}{\sqrt {R_{1}+R'_{1}n}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a6d5766001002cd6c4674e05813b4b6616ec1049)
on trouvera par la différentiation et par la comparaison des termes
![{\displaystyle {\begin{aligned}\mathrm {K} =&-{\frac {2\mathrm {U} _{1}}{R'_{1}}}+{\frac {4\mathrm {U} '_{1}R_{1}}{R_{1}^{'2}}}-{\frac {16\mathrm {U} ''_{1}R_{1}^{2}}{3R_{1}^{'3}}},\\\mathrm {L} =&{\frac {2\mathrm {U} _{1}}{R'_{1}}}-{\frac {8\mathrm {U} ''_{1}R_{1}}{3R_{1}^{'2}}},\\\mathrm {M} =&{\frac {2\mathrm {U} ''_{1}}{3R'_{1}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/99f2abbc564d43b0083a194d8c373a4652f50b73)
Complétant donc cette intégrale de la manière que nous l’avons dit, on aura, à la place de la dernière équation du numéro précédent, celle-ci
![{\displaystyle \Delta _{2}-\Delta _{0}={\frac {\mu \beta }{1+m}}\left(2\mathrm {V} _{1}+{\frac {2}{3}}\mathrm {V} ''_{1}+{\frac {\mathrm {K+L+M} }{\sqrt {R_{2}}}}-{\frac {\mathrm {K-L+M} }{\sqrt {R_{0}}}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0dad28988ac61a4c1f916bb354479ae2b9741253)
Soit : 2o
dans ce cas, la quantité
deviendra
et l’on aura à intégrer cette différentielle rationnelle
![{\displaystyle {\frac {\left(\mathrm {U} _{1}+\mathrm {U} '_{1}n+\mathrm {U} ''_{1}n^{2}\right)R_{1}^{\frac {3}{2}}}{\left(R_{1}+{\frac {1}{2}}R'_{1}n\right)^{3}}}dn,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc09320e356e19ac5984d0db3b54298a99199fa3)
qu’on supposera égale à
![{\displaystyle R_{1}^{\frac {3}{2}}\left(d{\frac {\mathrm {K+L} n}{\left(R_{1}+{\frac {1}{2}}R'_{1}n\right)^{2}}}+{\frac {\mathrm {M} dn}{R_{1}+{\frac {1}{2}}R'_{1}n}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f78fc9b00daf4057d501f74633178d4b8bba71c1)
ce qui donnera, en réduisant au même dénominateur et comparant les termes,
![{\displaystyle {\begin{aligned}\mathrm {K} =&-{\frac {\mathrm {U} _{1}}{R'_{1}}}-{\frac {2\mathrm {U} _{1}R_{1}}{R_{1}^{'3}}}+{\frac {12\mathrm {U} _{1}^{''2}R_{1}}{R_{1}^{'2}}},\\\mathrm {L} =&-{\frac {2\mathrm {U} '_{1}}{R'_{1}}}+{\frac {8\mathrm {U} ''_{1}R_{1}}{R_{1}^{'2}}},\\\mathrm {M} =&{\frac {4\mathrm {U} ''_{1}}{R_{1}^{'2}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3f14ba2bc6e3298f1af4253cdf8d52f4bb46c064)
Intégrant donc et complétant dûment l’intégrale, on trouvera, pour le cas dont il s’agit, l’équation
![{\displaystyle \Delta _{2}-\Delta _{0}={\frac {\mu \beta }{1+m}}\left[2\mathrm {V} _{1}+{\frac {2}{3}}\mathrm {V} ''_{1}+\left({\frac {\mathrm {K+L} }{R_{2}}}-{\frac {\mathrm {K-L} }{R_{0}}}\right){\sqrt {R_{1}}}+{\frac {\mathrm {M} R_{1}^{\frac {3}{2}}}{R'_{1}}}\log {\frac {R_{2}}{R_{0}}}\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10465dc5eff2e3223cd00f6cf715285bd5d280d8)
48. Ayant trouvé ainsi la valeur de
pour une portion d’anomalie excentrique
on trouvera de même la valeur de
pour une portion suivante d’anomalie
et ainsi de suite ; et ces différentes valeurs seront exactes, aux quantités près de l’ordre de
et
étant
et
étant la partie correspondante de l’anomalie moyenne de la planète. Ajoutant donc successivement ces valeurs ensemble, on aura la valeur totale de
répondant à une anomalie excentrique quelconque
et, faisant
on aura la valeur de
c’est-à-dire, l’accroissement de la quantité
pour une révolution entière de la comète.
Au reste il est bon de remarquer que les formules précédentes ne doivent proprement être employées que pour les parties de l’anomalie excentrique relativement auxquelles la quantité
sera assez petite et du même ordre que les différences finies
ce qui arrivera vers les minimum de distance entre la comète et la planète ; dans ces cas, les formules dont il s’agit ne sont sujettes à aucun inconvénient, et résolvent le Problème avec toute l’exactitude qu’on peut désirer ; au lieu que la méthode ordinaire des quadratures par les lignes paraboliques serait trop inexacte, à cause que les valeurs de
seront fort grandes et que leurs différences seront fort inégales.
Dans tout autre cas, c’est-à-dire, lorsque la distance entre la comète et la planète sera assez grande et que les variations de cette distance seront fort régulières, on emploiera avec succès la méthode ordinaire, tant pour intégrer la partie
que pour intégrer l’autre partie
et, comme cette méthode est très-connue et très en usage parmi les Géomètres, nous ne croyons pas devoir nous arrêter ici à l’expliquer ; les Ouvrages de Cotes et de Stirling renferment tout ce que l’on peut désirer sur ce sujet.
49. Quoiqu’on puisse, au moyen de ces différentes méthodes, calculer les variations des quantités
pour telles portions de l’orbite qu’on voudra, il ne sera cependant nécessaire de les employer que pour la partie inférieure de l’orbite, dans laquelle la distance de la comète au Soleil sera moindre ou ne sera pas beaucoup plus grande que la distance de la planète au Soleil ; car pour la partie supérieure de l’orbite, dans laquelle la distance de la comète au Soleil surpassera de beaucoup la distance de la planète au Soleil, il sera-bien plus avantageux d’employer la méthode des nos 39 et suivants, laquelle abrège et simplifie considérablement le calcul des perturbations dans cette partie.
Pour faire usage de cette méthode, il ne s’agit que de substituer dans les équations du no 37, à la place des valeurs de
qu’on a employées dans le no 42, celles de
(no 39) ; or nous avons déjà remarqué dans le no 13 que la quantité
n’est autre chose que les deux premiers termes de la quantité
réduite en série ascendante par rapport aux quantités
donc, comme (no 2)
![{\displaystyle \mathrm {R} ^{2}=r^{2}-2(x\xi +y\eta +z\zeta )+\rho ^{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6aad5c620835d3dbc322fe7774834eb5d5961513)
étant égal à
on aura, par la formule connue,
![{\displaystyle {\frac {1}{\mathrm {R} }}={\frac {1}{r}}+{\frac {x\xi +y\eta +z\zeta }{r^{3}}}-{\frac {\rho ^{2}}{2r^{3}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a030f2497cc374116596d3f618a4e88bfa554ba)
![{\displaystyle +{\frac {3(x\xi +y\eta +z\zeta )^{2}}{2r^{5}}}-{\frac {3(x\xi +y\eta +z\zeta )\rho ^{2}}{2r^{5}}}+{\frac {5(x\xi +y\eta +z\zeta )^{3}}{2r^{7}}}+\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/771c204856a844403a0eb77a04c9d89e15936830)
donc
![{\displaystyle {\frac {1}{\mathrm {S} }}={\frac {1}{r}}+{\frac {x\xi +y\eta +z\zeta }{r^{3}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b6aaa86c750246ba01a03746331e187675a560d7)
et par conséquent
![{\displaystyle \mathrm {{\frac {1}{S}}-{\frac {1}{R}}} ={\frac {\rho ^{2}}{2r^{3}}}-{\frac {3(x\xi +y\eta +z\zeta )^{2}}{2r^{5}}}+{\frac {3(x\xi +y\eta +z\zeta )\rho ^{2}}{2r^{5}}}-{\frac {5(x\xi +y\eta +z\zeta )^{3}}{2r^{7}}}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f0d1abd01e3733754123593afc5e0e2d52efcc14)
On différentiera maintenant cette quantité en faisant varier seulement
et les coefficients de
seront les valeurs de
donc, en supposant, pour abréger,
![{\displaystyle \varpi '=-{\frac {3\rho ^{2}}{2r^{5}}}+{\frac {15(x\xi +y\eta +z\zeta )^{2}}{2r^{7}}}-{\frac {15(x\xi +y\eta +z\zeta )\rho ^{2}}{2r^{7}}}+{\frac {35(x\xi +y\eta +z\zeta )^{3}}{2r^{9}}}-\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/90ce36dbaeac533b2f2b5cff650783e75134cf31)
![{\displaystyle \Pi '=-{\frac {3(x\xi +y\eta +z\zeta )}{r^{5}}}+{\frac {3\rho ^{2}}{2r^{5}}}-{\frac {15(x\xi +y\eta +z\zeta )^{2}}{2r^{7}}}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e8db3bd17c3a1267d30e4dcca2b93543d403a9ee)
on trouvera
![{\displaystyle \mathrm {X} '=\Pi '\xi +\varpi 'x,\quad \mathrm {Y} '=\Pi '\eta +\varpi 'y,\quad \mathrm {Z} '=\Pi '\zeta +\varpi 'z.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc0eacd89e2b03661b049d8447ac9a80df4ca328)
En comparant ces expressions de
![{\displaystyle \mathrm {X',Y',Z'} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/bad02fb7b26abf95b10f50f95e63d4bbbc52b753)
avec celles de
![{\displaystyle \mathrm {X,Y,Z} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e9852b5e32e3746828ae95b7c8664ea44d19655)
du
no 42, il est visible qu’elles n’en diffèrent qu’en ce que les quantités
![{\displaystyle \Pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/eed3e3db6cc2028a183af948212ed2551d25c954)
et
![{\displaystyle \varpi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e50d258418b5fa150a86b58f8d5eb40613e3ebf7)
se trouvent changées en
![{\displaystyle \Pi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/18754ba38b5204cfb3921f909b25d1c42b642613)
et
![{\displaystyle \varpi '.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bcdd04b113a9f9a144ebfb86e6f71265453cafd8)
D’où il est aisé de conclure que, par la substitution dont il s’agit, on aura les mêmes équations différentielles que dans le
no 42, en y changeant seulement
![{\displaystyle \Pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/eed3e3db6cc2028a183af948212ed2551d25c954)
et
![{\displaystyle \varpi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e50d258418b5fa150a86b58f8d5eb40613e3ebf7)
en
![{\displaystyle \Pi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/18754ba38b5204cfb3921f909b25d1c42b642613)
et
Il n’y aura donc qu’à employer dans les équations du no 42, à la place de
et
les quantités
et
et l’on pourra continuer à les employer pour telle portion de l’orbite qu’on voudra et reprendre ensuite les premières quantités, pourvu qu’on ajoute aux valeurs totales de
les quantités respectives
étant les valeurs de
du no 39 qui répondent au point de l’orbite où l’on change
en
et
étant les valeurs des mêmes quantités pour le point où l’on reprendra
et
à la place de
et
(no 40).
50. Le grand avantage de la transformation précédente consiste en ce que les quantités
et
qu’on substitue à la place de
et
deviennent très-petites lorsque la distance
de la comète au Soleil est beaucoup plus grande que la distance
de la planète au Soleil, ce qui est visible par les expressions des quantités
et
(numéro précédent) ; tandis que la valeur de
(no 42) demeure toujours finie, quel que soit l’éloignement de la comète, à cause du terme
qui ne dépend que de la distance de la planète au Soleil, et qui est l’effet de l’action de la planète sur le Soleil.
Or, si l’on considère que l’on a, en général,
![{\displaystyle {\begin{aligned}&\left(x^{2}+y^{2}+z^{2}\right)\left(\xi ^{2}+\eta ^{2}+\zeta ^{2}\right)\\&\quad =(x\xi +y\eta +z\zeta )^{2}+(x\eta -y\xi )^{2}+(x\zeta -z\xi )^{2}+(y\zeta -z\eta )^{2},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a484035b81172ad40951e29ffd700d4021f2e6e4)
et que, par conséquent,
est toujours nécessairement renfermé entre
et
on verra que le premier terme de la quantité
sera de l’ordre de
et les deux autres de l’ordre de
et que les deux premiers termes de
seront de l’ordre de
et les deux sui-
vants de l’ordre de
![{\displaystyle {\frac {\rho ^{3}}{r^{6}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc4fd0710647b48c1298bd98a58874a4132f0da2)
et ainsi de suite. Donc, lorsque
![{\displaystyle r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0d1ecb613aa2984f0576f70f86650b7c2a132538)
est assez grand vis-à-vis de
![{\displaystyle \rho ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f12bed314e4bed19299ed16afd79f67ea5c4593c)
en sorte que
![{\displaystyle {\frac {1}{\mathrm {R} ^{3}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a3aff51e7b95e437e2a6f52b8414ce57227191d)
diffère peu de
![{\displaystyle {\frac {1}{r^{3}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d58fac6aeeaecee09f3df8b9d972854caedea3f5)
le rapport de
![{\displaystyle \Pi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/18754ba38b5204cfb3921f909b25d1c42b642613)
à
![{\displaystyle \Pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/eed3e3db6cc2028a183af948212ed2551d25c954)
sera de l’ordre de
![{\displaystyle {\frac {\rho ^{4}}{r^{4}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ac6a17956057d52c00bbf200d8b0c58f235896cc)
et celui de
![{\displaystyle \varpi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ea30ab430af5278b4172b9ca39ff39a64f8a24b2)
à
![{\displaystyle \varpi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e50d258418b5fa150a86b58f8d5eb40613e3ebf7)
de l’ordre de
![{\displaystyle {\frac {\rho ^{2}}{r^{2}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f179131f7b96f663e79269a311e2d04e30a04dd)
la quantité
![{\displaystyle \varpi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e50d258418b5fa150a86b58f8d5eb40613e3ebf7)
étant déjà elle-même très-petite de l’ordre de
Donc, lorsque
sera devenu
ou
on pourra, du moins dans la première approximation, négliger les quantités
et
comme nulles ; ou, si l’on veut absolument y avoir égard, il suffira d’y tenir compte des premiers termes. Dans ce cas, on pourra en toute sûreté employer la méthode ordinaire des quadratures mécaniques pour intégrer les quantités
mais on pourra aussi les intégrer analytiquement, du moins par approximation ; c’est ce que nous allons faire voir.
51. Pour cet effet, on commencera par remettre dans les expressions des quantités
du no 42, à la place de
et
leurs valeurs en
et
savoir
![{\displaystyle \cos u={\frac {2x}{a}}+f,\quad \sin u={\frac {y}{\sqrt {ah}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8090e4bb807c3fe40d3f0070c351486881a48edb)
moyennant quoi ces quantités deviendront des fonctions rationnelles et entières de
et de
dans lesquelles les quantités
ne passeront pas la seconde dimension, excepté les expressions de
et de
où ces quantités monteront à la quatrième dimension ; mais je remarque, à l’égard de l’expression, de
qu’on y peut réduire les dimensions de
à la troisième. En effet il est visible que les termes qui, dans cette expression, peuvent donner des dimensions de
plus hautes que la troisième, sont ceux-ci
![{\displaystyle -{\frac {2r^{2}}{f{\sqrt {a^{3}h}}}}(\mathrm {G} )+{\frac {fyx}{2{\sqrt {ah^{3}}}}}(\mathrm {F} ),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8e77c57fb8060acbae941f4690d8ddd305a2d22a)
autant que les valeurs de
et
contiennent
élevées à la
seconde dimension. Or, en faisant, pour un moment,
![{\displaystyle {\sqrt {\frac {a}{4h}}}\left[(fr+x)\xi +y\eta \right]=\Xi ,\quad y\xi -x\eta =\Upsilon ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/abfdfa43d39f2e2804dea66697f912e469a13b99)
on a
![{\displaystyle (\mathrm {F} )=-\Xi y+\left(4x{\sqrt {\frac {h}{a}}}+2f{\sqrt {ah}}\right)\Upsilon ,\quad (\mathrm {G} )=\Xi y+{\sqrt {\frac {a}{h}}}\,y\Upsilon \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/22d2099cf11f0c4ba24ab3a7d4240ea1f9c7efc4)
donc les termes en question seront
![{\displaystyle -\left({\frac {2r^{2}x}{f{\sqrt {a^{3}h}}}}+{\frac {fy^{2}x}{2{\sqrt {ah^{3}}}}}\right)\Xi +\left(-{\frac {2r^{2}y}{fah}}+{\frac {2fyx^{2}}{ah}}\right)\Upsilon .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/98a7559de6969b505ff6f6d581337cb4c3775fc8)
Maintenant, à cause que nous prenons le grand axe de l’orbite pour celui des abscisses
et que
(no 25), on aura (no 18)
![{\displaystyle x={\frac {2h-r}{f}},\quad y={\frac {2{\sqrt {h}}}{f}}{\sqrt {r-{\frac {r^{2}}{a}}-h}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7fa8839afe2b8c04adcb6a4c53e29d5ec8684765)
donc, substituant cette valeur de
dans le coefficient de
il deviendra
![{\displaystyle {\frac {2r^{2}x}{f{\sqrt {a^{3}h}}}}+{\frac {2x}{f{\sqrt {ah}}}}\left(r-{\frac {r^{2}}{a}}-h\right)={\frac {2x(r-h)}{f{\sqrt {ah}}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b72ba51ac429003ee1b3c7d3c455ad753bc5ff36)
et, substituant la valeur de
dans le coefficient de
il deviendra
![{\displaystyle -{\frac {2r^{2}y}{fah}}+{\frac {2r^{2}y}{fah}}+{\frac {8(h-r)y}{fa}}={\frac {8(h-r)y}{fa}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0b140753ed66981850aa040259eaad3302b5ae64)
Donc, puisque
et
ne contiennent que la première dimension de
il s’ensuit que les termes dont il s’agit de l’expression de
lesquels paraissent, au premier aspect, devoir contenir la quatrième dimension de ces quantités, n’en contiendront réellement que la troisième.
Cela supposé, on mettra, tant dans les expressions de
du numéro cité que dans celles de
et
du no 48, à la place de
et
les valeurs
(no 18),
étant l’anomalie vraie de la comète dans son orbite non altérée ; on verra :
1o Que les expressions de
deviendront des fonctions rationnelles et entières de
et de
et
dans lesquelles
ne montera au plus qu’au second degré, à l’exception des quantités
et
dont la première contiendra
et dont la seconde contiendra
2o Que les expressions de
et
deviendront, à cause de
![{\displaystyle {\begin{aligned}\Pi '=&{\frac {3(\xi \cos \varphi +\eta \sin \varphi )}{r^{4}}}+{\frac {3\rho ^{2}-2(\xi \cos \varphi +\eta \sin \varphi )^{2}}{r^{5}}}+\ldots ,\\\varpi '=&{\frac {-3\rho ^{2}+15(\xi \cos \varphi +\eta \sin \varphi )^{2}}{r^{5}}}\\&+{\frac {-15(\xi \cos \varphi +\eta \sin \varphi )\rho ^{2}+35(\xi \cos \varphi +\eta \sin \varphi )^{3}}{r^{6}}}+\ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f0dee8f4b015ca1dbd0da4ca5762185a4201822b)
On substituera maintenant ces valeurs de
dans les expressions des différentielles
du no 42, et l’on y mettra, à la place de
et
les valeurs précédentes de
et
enfin on mettra pour
sa valeur
déduite de l’équation (no 21)
![{\displaystyle dt={\frac {r^{2}d\varphi }{\sqrt {2h(1+m)}}}={\sqrt {\frac {a}{2(1+m)}}}\,rdu.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/00bf1e342b57d011bef323366ac26d288c7df3e8)
52. Il est aisé de voir que, par ces différentes substitutions, les valeurs des différentielles
du no 42 se trouveront composées de différents termes de la forme
![{\displaystyle {\frac {\Sigma \cos ^{m}\varphi \sin ^{n}\varphi d\varphi }{r^{p}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/32b14231b7d1b579b1400316c09fe5e5820b648a)
étant des nombres entiers positifs ou zéro, et
étant une fonction rationnelle et entière de
j’en excepte seulement les termes de la valeur de
qui seront multipliés par l’angle
et que nous examinerons plus bas. Et il n’est pas difficile de prouver que
ne sera pas
pour les premiers termes de
et
ni
pour les termes suivants, et ainsi du reste.
Or (no 18)
![{\displaystyle r={\frac {2h}{1+f\cos \varphi }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/317d8f907fa04fd8106ee3a8c6bec0e14e6e5fee)
à cause de
donc, si l’on substitue cette valeur dans la formule
précédente, on n’aura dans les valeurs de
![{\displaystyle d\,\delta h,d\,\delta a,d\,\delta f,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/63fd9bea767a8d37409d69140a69aef0b26b40e6)
que des termes de cette forme
![{\displaystyle \mu }](https://wikimedia.org/api/rest_v1/media/math/render/svg/9fd47b2a39f7a7856952afec1f1db72c67af6161)
et
![{\displaystyle \nu }](https://wikimedia.org/api/rest_v1/media/math/render/svg/c15bbbb971240cf328aba572178f091684585468)
étant des nombres entiers positifs, tels que
![{\displaystyle \mu +\nu >5}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78e916f7cfb90347fac9bcb2966f3566ff12d85f)
pour les premiers termes de
![{\displaystyle \Pi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/18754ba38b5204cfb3921f909b25d1c42b642613)
et
![{\displaystyle \varpi ',}](https://wikimedia.org/api/rest_v1/media/math/render/svg/187ed0ba0eee302580758b17d5b6476b645edc64)
ni
![{\displaystyle >7}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5bd6067b9cbad266a62c396ada9da8bd8d26bfa)
pour les termes suivants ; j’excepte toujours les termes affectés de
![{\displaystyle t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/65658b7b223af9e1acc877d848888ecdb4466560)
dans la valeur de
![{\displaystyle d\,\delta i}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9d39f24980e08d7eb2f648a8b4ddbacbb946b91a)
(
voir ci-après le
no 56).
Qu’on substitue maintenant dans
à la place de
leurs valeurs en sinus et cosinus de
(no 43) ; et pour cela on remarquera que, à cause de la petitesse des quantités
et
on peut sans scrupule négliger l’effet de l’excentricité de la planète, et faire simplement
![{\displaystyle \rho ={\frac {\alpha }{2}},\quad \lambda =\mathrm {M} -\Lambda =(\theta -\mathrm {T} ){\sqrt {\frac {a^{3}}{\alpha ^{3}}}}-\Lambda ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7b72b3769f45672a2134bea967aabed475b31cb0)
en dénotant par
l’anomalie vraie de la planète qui répond au nœud ascendant de son orbite sur l’orbite non altérée de la comète ; mais, si l’on voulait absolument avoir égard à l’excentricité de l’orbite de la planète, il n’y aurait qu’à ajouter aux valeurs moyennes de
et de
les inégalités du rayon vecteur et de la longitude de la planète, inégalités dont les premières sont représentées par une suite très-convergente de termes qui procèdent suivant les cosinus de
et dont les autres sont représentées par une semblable suite, mais qui procède suivant les sinus des mêmes angles. (Voir les pages 6 et 8 des Tables astronomiques de Berlin, où
dénote l’anomalie moyenne que nous désignons ici par
)
Ces substitutions rendront la quantité
de la forme
![{\displaystyle \mathrm {A+B\sin \upsilon +C\cos \upsilon +D\sin 2\upsilon +\ldots } ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3397e4665dd588eafe262db95e0c52a698f26f59)
les coefficients
étant constants, et l’angle
étant ![{\displaystyle =\theta {\sqrt {\frac {a^{3}}{\alpha ^{3}}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f3f3715b1cd21171b7e0ae5183868207545fd4f4)
Ainsi les valeurs des différentielles
se trouveront composées de deux sortes de termes : les uns indépendants de l’angle
c’est-à-dire, du mouvement moyen de la planète, les autres affectés des sinus ou cosinus de cet angle ou de ses multiples.
53. À l’égard des termes de la première espèce, il est clair qu’ils seront de la forme
![{\displaystyle \cos ^{\mu }\varphi \sin ^{\nu }\varphi d\varphi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5b9e0039d7d91482b68ac4dfba345e095f407f5c)
et par conséquent tous intégrables,
et
étant, par l’hypothèse, des nombres entiers positifs.
Quant à ceux de l’autre espèce, ils seront évidemment de la forme
![{\displaystyle \cos ^{\mu }\varphi \sin ^{\nu }\varphi \sin \mathrm {N} \upsilon d\varphi ,\quad {\text{ou}}\quad \cos ^{\mu }\varphi \sin ^{\nu }\varphi \cos \mathrm {N} \upsilon d\varphi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3de23c0c61d3fcea02c50534d11254f81cfdcf18)
étant un nombre entier. Ces termes ne sont intégrables par aucune méthode connue ; mais nous allons faire voir que, dans la partie supérieure de l’orbite de la comète, à laquelle est destinée la méthode que nous exposons, ces termes seront considérablement plus petits que les précédents ; en sorte qu’on pourra le plus souvent les négliger sans scrupule.
Pour cet effet, je remarque que
![{\displaystyle d\upsilon =d\theta {\sqrt {\frac {a^{3}}{\alpha ^{3}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be1f7ad394a4d3540302b5fe572d8001a732a0e0)
mais (no 20)
![{\displaystyle d\theta ={\sqrt {\frac {8(1+m)}{a^{2}}}}dt,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e603df1d7e9051d8cc7d13e5e70e1ccbb7ee976)
et (no 21)
![{\displaystyle dt={\frac {r^{2}d\varphi }{\sqrt {2h(1+m)}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ebf3093aca7e6fd8a706b2aef839d24679b6ac93)
donc
![{\displaystyle d\upsilon ={\frac {2r^{2}d\varphi }{\sqrt {\alpha ^{3}h}}},\quad \mathrm {et\ de\ l{\grave {a}}} \quad d\varphi ={\frac {d\upsilon {\sqrt {\alpha ^{3}h}}}{2r^{2}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e36574ba00308c7bf4c469157a064cc621f04f6c)
Si l’on substitue cette valeur de
dans les termes dont il s’agit, et qu’on fasse, pour abréger,
![{\displaystyle {\sqrt {\alpha ^{3}h}}{\frac {\cos ^{\mu }\varphi \sin ^{\nu }\varphi }{2\mathrm {N} r^{2}}}=\Phi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a810483413f1fefd554ac7b52bfcfa9d66e4991b)
ils deviendront
![{\displaystyle -\Phi d\cos \mathrm {N} \upsilon ,\quad \Phi d\sin \mathrm {N} \upsilon ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47e8c5dbb015ca09e38307caeeb33385a30f0eaa)
dont l’intégrale est
![{\displaystyle -\Phi \cos \mathrm {N} \upsilon +\int \cos \mathrm {N} \upsilon d\Phi ,\quad \Phi \sin \mathrm {N} \upsilon -\int \sin \mathrm {N} \upsilon d\Phi .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f45f0449ea3a2a14bb5f3dc1ab6ab9ddbee6002b)
Les expressions
et
représentent, comme l’on voit, les aires des courbes qui auraient
pour abscisse et
ou
pour ordonnée ; et il est facile de concevoir que l’aire totale de chacune de ces courbes sera toujours moindre (abstraction faite du signe) que le produit de l’abscisse totale par lâ plus grande ordonnée, laquelle est égale à
de sorte que, dénotant par
cette abscisse totale, on aura
pour les deux limites entre lesquelles seront nécessairement renfermées les aires
et
Or, dans la partie supérieure de l’orbite, la distance
de la comète au Soleil est supposée beaucoup plus grande que la distance moyenne
de la planète au Soleil ; de plus la distance périhélie
égale à
à très-peu près, est dans la plupart des comètes, et surtout dans celles dont on attend le retour, moindre que l’unité, distance moyenne de la Terre au Soleil ; de sorte que la quantité
sera nécessairement fort petite. Par conséquent les quantités
et
seront beaucoup plus petites, généralement parlant, que la valeur de
Il faut remarquer au reste que, pour avoir la valeur de
pour toute la partie supérieure de l’orbite, c’est-à-dire, la valeur totale de l’intégrale de
pour cet espace, il faut prendre les éléments
toujours avec le même signe. Si donc dans tout cet espace la quantité
n’a ni maximum ni minimum, on prendra l’intégrale à la manière ordinaire, et l’on aura pour
la différence entre les deux valeurs extrêmes de
Mais si entre ces valeurs extrêmes il se trouve des maximum et des minimum, alors la valeur exacte de
sera égale au double de la différence entre la somme de toutes les plus grandes valeurs de
et la somme de toutes les plus petites, en regardant les maximum négatifs comme des minimum et les minimum négatifs comme des maximum, et comptant les deux valeurs extrêmes de
parmi les maximum ou minimum, suivant que
va en diminuant ou en augmentant, mais en ne prenant que la moitié de chacune de ces valeurs. C’est de quoi l’on peut se convaincre aisément par l’inspection d’une figure parabolique quelconque qui aurait différents maximum et minimum.
Or je dis que, si
la quantité
n’aura ni maximum ni minimum lorsque
sera impair, et qu’elle aura un seul minimum au périhélie où
lorsque
sera pair. En effet, à cause de
![{\displaystyle r={\frac {2h}{1+f\cos \varphi }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/317d8f907fa04fd8106ee3a8c6bec0e14e6e5fee)
on aura
![{\displaystyle \cos \varphi ={\frac {1}{f}}\left({\frac {2h}{r}}-1\right)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2161261710432a10301b2ad26ab1a28af460f833)
en sorte que, si
sera négatif et ne changera point de signe, mais
sera positif en deçà de l’aphélie et deviendra négatif au delà. Or
![{\displaystyle \Phi ={\sqrt {\alpha ^{3}h}}{\frac {\cos ^{\mu }\varphi \sin ^{\nu }\varphi }{2\mathrm {N} r^{2}}}={\frac {\sqrt {\alpha ^{3}h}}{(-f)^{\mu }2\mathrm {N} }}\left(1-{\frac {2h}{r}}\right)^{\mu }{\frac {\sin ^{\nu }\varphi }{r^{2}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/851d94e78d1bff71113f3c863a14a59cb0c8275f)
mais la quantité
![{\displaystyle {\frac {1}{r^{2}}}\left(1-{\frac {2h}{r}}\right)^{\mu }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d5e026f254126df236e21d3173df93e706d16c3)
diminue à mesure que
augmente, et vice versâ, du moins tant que
puisque sa différentielle est
![{\displaystyle -2\left(1-{\frac {2h}{r}}\right)^{\mu -1}\left[1-{\frac {(\mu +2)h}{r}}\right]{\frac {dr}{r^{3}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2cf4c6672ba17ebd3c06054e969185f25afcc261)
donc, si
est impair, la quantité
![{\displaystyle \left(1-{\frac {2h}{r}}\right)^{\mu }{\frac {\sin ^{\nu }\varphi }{r^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a64f5072b9c05699f807b041b1512587e4122212)
ira en diminuant jusqu’à l’aphélie où elle sera nulle, et continuera à diminuer au delà de l’aphélie où elle sera négative ; mais, si
est pair, la même quantité, après avoir diminué jusqu’à l’aphélie, augmentera de nouveau au delà de l’aphélie, en demeurant toujours positive.
Donc, si l’on suppose que la partie supérieure de l’orbite commence au point où
et finisse au point semblablement situé au delà de l’aphélie où
et
on aura, pourvu que
![{\displaystyle (\Phi )={\sqrt {\alpha ^{3}h}}{\frac {\cos ^{\mu }\varphi '\sin ^{\nu }\varphi '}{\mathrm {N} r'^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eec0ddb65c373c1d442765930879b8a633d1d2d4)
si
est impair, et
![{\displaystyle (\Phi )={\sqrt {\alpha ^{3}h}}\left[{\frac {\cos ^{\mu }\varphi '\sin ^{\nu }\varphi '}{\mathrm {N} r'^{2}}}-{\frac {\cos ^{\mu }180^{\circ }\sin ^{\nu }180^{\circ }}{\mathrm {N} \left({\cfrac {a(1+e)}{2}}\right)^{2}}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a5a35732c90cb38a8cb59fa4769bcd1490286ed)
si
est pair,
étant la valeur de
dans l’aphélie.
À l’égard de la condition de
comme nous avons vu (no 51) que
ne peut être
pour les premiers termes de
et
auxquels il suffire le plus souvent d’avoir égard, il est clair que cette condition aura toujours lieu dans la partie supérieure de l’orbite où l’on suppose
beaucoup plus grand que
puisque pour Jupiter et Saturne, qui sont les seules planètes qu’on ait à considérer dans la Théorie des perturbations des comètes, on a à peu près
ou
54. Si les limites
n’étaient pas assez petites, en sorte qu’on ne crût pas pouvoir négliger les quantités renfermées entre ces limites, on pourrait les resserrer davantage de la manière suivante.
Les deux différentielles
![{\displaystyle \cos \mathrm {N} \upsilon d\Phi ,\quad \sin \mathrm {N} \upsilon d\Phi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c0a6ada2cc006ded51d694d4a00ea41df43fb39)
étant mises sous la forme
![{\displaystyle {\frac {d\Phi }{d\varphi }}\cos \mathrm {N} \upsilon d\varphi ,\quad {\frac {d\Phi }{d\varphi }}\sin \mathrm {N} \upsilon d\varphi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/764997c2a8f51bd133e0cc07138ebc4ab9cbe390)
se changent, par la substitution de
à
et par la supposition de
en celles-ci
![{\displaystyle -\Phi 'd\cos \mathrm {N} \upsilon ,\quad \Phi 'd\sin \mathrm {N} \upsilon ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cbdfbed968b3d7b0087ee384fc34c031b76ac0ef)
dont l’intégrale est
![{\displaystyle -\Phi '\cos \mathrm {N} \upsilon +\int \cos \mathrm {N} \upsilon d\Phi ',\quad \Phi '\sin \mathrm {N} \upsilon -\int \sin \mathrm {N} \upsilon d\Phi '\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e07d167b1647cba6dc1a7ab706b90c918c8361f)
et l’on pourra appliquer aux quantités
![{\displaystyle \int \cos \mathrm {N} \upsilon d\Phi ',\quad \int \sin \mathrm {N} \upsilon d\Phi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2515ed9073b3ab21569f4f4d690837ff9450e3de)
les mêmes raisonnements que nous avons faits dans le numéro précédent ainsi, dénotant par
la valeur totale de l’intégrale de
prise comme nous l’avons dit dans ce numéro, on aura de nouveau
pour les limites entre lesquelles seront renfermées les valeurs des quantités dont il s’agit.
Or il est facile de se convaincre que la quantité
est nécessairement beaucoup plus petite que la quantité
lorsque
est assez grand vis-à-vis de
ainsi, en négligeant les intégrales renfermées entre ces dernières limites, on commettra une erreur bien plus petite que celle qui pourrait résulter de l’omission des intégrales renfermées dans les limites du numéro précédent.
On voit par là comment on pourrait s’y prendre pour pousser cette approximation plus loin, et diminuer à volonté l’erreur résultant des intégrales qu’on négligerait ; mais il suffira, dans la plupart des cas, de s’en tenir à l’approximation du numéro précédent.
55. Il nous reste encore à examiner les termes multipliés par l’angle
dans la différentielle
termes que nous avons expressément exceptés (no 51). Or on voit, par la valeur générale de
du no 37 (Section précédente), que les termes dont il s’agit ne peuvent venir que du terme
![{\displaystyle 3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {d\,\delta a}{a}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c1c4912a44517cdd3523fb41f2f44196a37df07e)
il suffit donc de considérer ce terme et d’en chercher l’intégrale, en supposant que l’on mette dans la valeur de
les quantités
à la place des quantités
(no 48).
Je reprends pour cela l’expression générale de la différentielle
du même no 37, laquelle est
![{\displaystyle d\,\delta a=-{\frac {\mu }{1+m}}a^{2}(\mathrm {X} dx+\mathrm {Y} dy),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a2219591ff66fbab2396ec176f2bf6de98154830)
et, pour embrasser en même temps toute la généralité possible, je remarque que, si l’on n’avait pas supposé
et
et qu’on eût par conséquent employé dans les calculs de ce numéro la valeur complète de
du no 30 à la place de celle du no 31, on eût trouvé cette expression plus générale de
savoir
![{\displaystyle d\,\delta a=-{\frac {\mu }{1+m}}a^{2}(\mathrm {X} dx+\mathrm {Y} dy+\mathrm {Z} dz).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d3905f7d09baa52a8453c3505bade69fe5551632)
Qu’on change maintenant, dans cette expression, les quantités
en
et qu’on y substitue ensuite, à la place de ces dernières quantités, leurs valeurs, lesquelles
en faisant, pour abréger,
sont exprimées ainsi (no 39)
il est visible que la différentielle
![{\displaystyle \mathrm {X} 'dx+\mathrm {Y} 'dy+\mathrm {Z} 'dz}](https://wikimedia.org/api/rest_v1/media/math/render/svg/39d76038ff4e80ff3317813935a426686da7be5d)
ne sera autre chose que la différence de
prise en faisant varier seulement les quantités
qui appartiennent à l’orbite de la comète, et en regardant comme constantes les coordonnées
de l’orbite de la planète.
De sorte que, si l’on désigne par la caractéristique
cette différence partielle, on aura, en général,
![{\displaystyle d\,\delta a=-{\frac {\mu }{1+m}}a^{2}\operatorname {D} \mathrm {P} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a34dba347c114ba276dad7d5e2e73f978552429)
Or on a, par le
no 48,
![{\displaystyle \mathrm {P} ={\frac {\rho ^{2}}{2r^{3}}}-{\frac {3(x\xi +y\eta +z\zeta )^{2}}{2r^{5}}}+{\frac {3(x\xi +y\eta +z\zeta )\rho ^{2}}{2r^{5}}}-{\frac {5(x\xi +y\eta +z\zeta )^{3}}{2r^{7}}}+\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/229e2d449b2adf8cfba681c6e1e6f9d7ac4a0a40)
Et si l’on substitue, dans cette expression de
les valeurs de
en sinus et cosinus de
(no 51), il est visible qu’elle deviendra de cette forme
![{\displaystyle R+\mathrm {R\sin \upsilon +S\cos \upsilon +T\sin 2\upsilon +V} \cos 2\upsilon +\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/31674110965ec41bdb298caaefec7ac88e315090)
dans laquelle
seront des fonctions rationnelles et entières de
et dont chaque terme sera de plus divisé par une puissance de
dont l’exposant surpassera de trois unités ou davantage la somme des dimensions de
dans le numérateur.
Or, comme l’angle
dépend uniquement des quantités
qui doivent être regardées comme constantes dans la différence partielle
et qu’au contraire les quantités
dépendent uniquement des quantités
qui sont les seules variables dans cette différentielle, il est clair qu’on aura
![{\displaystyle \operatorname {D} \mathrm {P} =dR+\sin \upsilon \mathrm {dR+\cos \upsilon dS+\sin 2\upsilon dT+\cos 2\upsilon dV} +\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80afe88f1b9b37ccb438265ecccf8bfac43c16a2)
étant des différences ordinaires et totales des quantités ![{\displaystyle R,\mathrm {R,S} ,\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d825cde732bf7df070b4890e9855a770e1211df)
Donc on aura, en général,
![{\displaystyle d\,\delta a=-{\frac {\mu a^{2}}{1+m}}\left(dR+\sin \upsilon \mathrm {dR+\cos \upsilon dS+\sin 2\upsilon dT+\cos 2\upsilon dV} +\ldots \right)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2cc8a058b2e7c97d0dcbde5021182cf8b64ca982)
et cette valeur de
en y faisant
et
deviendra identique avec celle du no 50, mais elle sera toujours d’une forme plus simple et plus commode pour l’intégration.
56. En effet on voit d’abord, par l’expression précédente de
que la partie indépendante de l’angle
est intégrable, son intégrale étant
où
est la partie indépendante de
dans la valeur de
laquelle sera par conséquent une fonction rationnelle et entière de
et
en faisant
![{\displaystyle x=r\cos \varphi ,\quad y=r\sin \varphi ,\quad z=0,\quad r={\frac {2h}{1+f\cos \varphi }}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ea959bc488ebcbaefa321b374ac7c3965e8220a)
De là on tire cette conclusion importante, que la valeur
c’est-à-dire, l’altération du grand axe de l’orbite de la comète, en tant qu’elle vient des perturbations de la partie supérieure de l’orbite, ne contient aucun terme proportionnel à l’angle, et qui puisse par conséquent augmenter continuellement.
À l’égard des autres termes de la valeur de
il est clair qu’après la substitution des valeurs de
en
ils deviendront de la forme
![{\displaystyle \cos ^{\mu }\varphi \sin ^{\nu }\varphi \sin \mathrm {N} \upsilon d\varphi ,\quad {\text{ou}}\quad \cos ^{\mu }\varphi \sin ^{\nu }\varphi \cos \mathrm {N} \upsilon d\varphi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3de23c0c61d3fcea02c50534d11254f81cfdcf18)
et pourront être traités par la méthode des nos 52 et suivants.
57. Venons maintenant au terme
![{\displaystyle 3t{\sqrt {\frac {2(1+m)}{a^{3}}}}{\frac {d\,\delta a}{a}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/02f9b18447a73cd34b3df173fc594616e5ae7a7e)
de la valeur de
En y substituant d’abord pour
la quantité
indépendante de
on aura la différentielle
![{\displaystyle -3{\sqrt {\frac {2}{a(1+m)}}}tdR,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/93b5e964c168069369c1c50877478a6d7a5a2881)
dont l’intégrale est
![{\displaystyle -3{\sqrt {\frac {2}{a(1+m)}}}\left(tR-\int Rdt\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4bdbcd3028aee57bfcce97535645279971baebac)
Or on a (no 21),
![{\displaystyle dt={\frac {r^{2}d\varphi }{\sqrt {2h(1+m)}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ebf3093aca7e6fd8a706b2aef839d24679b6ac93)
donc, comme dans l’expression de
les exposants négatifs de
surpassent de trois unités ou davantage la somme des exposants positifs de
(no 54), il est visible qu’en mettant dans
pour
et
![{\displaystyle r\cos \varphi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/7d08a45e0026eac9b3f951e32219fb96768a4119)
et
![{\displaystyle r\sin \varphi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a8cbf964f90ca39007d19f233b82a0e217f27d9)
(
![{\displaystyle z}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf368e72c009decd9b6686ee84a375632e11de98)
étant égal à zéro), la quantité
![{\displaystyle r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0d1ecb613aa2984f0576f70f86650b7c2a132538)
ne s’y trouvera encore qu’au dénominateur ; en sorte que, substituant ensuite
![{\displaystyle {\frac {2h}{1+f\cos \varphi }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9fe917c16396e56dd2b52d81f6ba24530f8e3cca)
pour
![{\displaystyle r,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/250644a0f511e9078be6f89ba78a606a0e08c0a0)
la quantité
![{\displaystyle Rr^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/733fc8497e0a17e872f2733fa897cb2965c47f89)
deviendra une fonction rationnelle et entière de
![{\displaystyle \sin \varphi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5476adfa219ffe9cc3186d8226214b0f4734dce)
et
![{\displaystyle \cos \varphi \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a48d5a7833ac32b5adf33c3bbbbe66f49d205a2a)
d’où il s’ensuit que
![{\displaystyle Rdt={\frac {Rr^{2}d\varphi }{\sqrt {2h(1+m)}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6ff032026b6a966b24e5069d5073490dfa99659c)
sera tout à fait intégrable.
Quant à l’autre partie de la valeur de
elle sera composée, comme nous l’avons vu ci-dessus, de termes de la forme
![{\displaystyle -{\frac {\mu a^{2}}{1+m}}\cos ^{\mu }\varphi \sin ^{\nu }\varphi \sin \mathrm {N} \upsilon d\varphi ,\ \ {\text{ou}}\ \ -{\frac {\mu a^{2}}{1+m}}\cos ^{\mu }\varphi \sin ^{\nu }\varphi \cos \mathrm {N} \upsilon d\varphi \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2689dc9da28cafeca491f672514b61089d10afc6)
donc les termes qui en résulteront dans la valeur de
seront de la forme
![{\displaystyle -3\mu {\sqrt {\frac {2}{a(1+m)}}}\,t\cos ^{\mu }\varphi \sin ^{\nu }\varphi \sin \mathrm {N} \upsilon d\varphi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1d4293d66bdde5d8d83db11b912bfe0c226acece)
![{\displaystyle {\text{ou}}\quad -3\mu {\sqrt {\frac {2}{a(1+m)}}}\,t\cos ^{\mu }\varphi \sin ^{\nu }\varphi \cos \mathrm {N} \upsilon d\varphi .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aa87587e8f60236fbc01dbc84a8b170f5e2defac)
Ainsi il suffira de considérer les différentielles
![{\displaystyle t\cos ^{\mu }\varphi \sin ^{\nu }\varphi \sin \mathrm {N} \upsilon d\varphi ,\quad {\text{et}}\quad t\cos ^{\mu }\varphi \sin ^{\nu }\varphi \cos \mathrm {N} \upsilon d\varphi .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a2d5886ad7038be09151ba1285469ed1382433c)
À l’imitation de ce qu’on a fait plus haut (no 52), on substituera, dans ces différentielles,
au lieu de
et faisant, pour abréger
à cause de
![{\displaystyle {\begin{aligned}\mathrm {V} =&\int t\sin \mathrm {N\upsilon .N} d\upsilon =-t\cos \mathrm {N} \upsilon +{\sqrt {\frac {\alpha ^{3}}{8(1+m)}}}\mathrm {\frac {\sin N\upsilon }{N}} ,\\\mathrm {W} =&\int t\cos \mathrm {N\upsilon .N} d\upsilon =\quad t\sin \mathrm {N} \upsilon -{\sqrt {\frac {\alpha ^{3}}{8(1+m)}}}\mathrm {\frac {\cos N\upsilon }{N}} ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d359bdaf0bb8ff5055e17f5d12a27cc9d4acc1f)
on aura ces transformées
et
en conservant la valeur de
du numéro cité.
Intégrant par parties, on aura
![{\displaystyle \Phi \mathrm {V} -\int \mathrm {V} d\Phi \quad {\text{et}}\quad \Phi \mathrm {W} -\int \mathrm {W} d\Phi ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d083f4be3782165bbdf7b67081f9061f7631813b)
et l’on démontrera, par un raisonnement analogue à celui de ce numéro, que les valeurs des intégrales
et
seront renfermées entre les limites
et
en désignant par
et
les plus grandes valeurs de
et
dans la partie supérieure de l’orbite, et conservant la valeur de
de l’endroit cité. Or les maximum de
et
ayant lieu lorsque
ou
c’est-à-dire, lorsque
ou
il s’ensuit que les plus grandes valeurs des quantités
et
seront égales à
(abstraction faite du signe). Si donc on désigne par
la valeur de
qui répond à toute la partie supérieure de l’orbite, c’est-à-dire, la valeur de
pour le point où finit cette partie de l’orbite, on aura
et
et les valeurs des intégrales
pour toute la partie supérieure de l’orbite, seront renfermées entre ces limites ![{\displaystyle \pm (t)(\Phi ).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b08085cec2a5775a69eb2f2448b0ef81671016d5)
58. Si l’on ne jugeait pas les limites assez approchées, surtout à cause que la valeur de
peut être assez considérable, on pourrait les resserrer davantage par une méthode analogue à celle du no 53.
En effet, en conservant la valeur de
de ce numéro, et faisant, pour abréger,
![{\displaystyle {\begin{aligned}\mathrm {V} _{1}=&\int \mathrm {VN} d\upsilon =-\mathrm {W} -{\sqrt {\frac {\alpha ^{3}}{8(1+m)}}}\mathrm {\frac {\cos N\upsilon }{N}} ,\\\mathrm {W} _{1}=&\int \mathrm {WN} d\upsilon =\quad \mathrm {V} -{\sqrt {\frac {\alpha ^{3}}{8(1+m)}}}\mathrm {\frac {\sin N\upsilon }{N}} ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7311dbb796447b734efbdfa943a4b19395608764)
on transformera les différentielles
et
en celles-ci
et
dont l’intégrale, prise par parties, sera
![{\displaystyle \Phi '\mathrm {V} _{1}-\int \mathrm {V} _{1}d\Phi '\quad {\text{et}}\quad \Phi '\mathrm {W} _{1}-\int \mathrm {W} _{1}d\Phi '\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e6328db9dc44eb1fc87b4a2f92ec71956cec47d3)
et l’on démontrera de la même manière que, si
![{\displaystyle (\mathrm {V} _{1})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b481d60a259e7f24f09a39930c27fccfedcc2dcc)
et
![{\displaystyle (\mathrm {W} _{1})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c78afb15f16241465d176932a60a4e801c7af7af)
sont les plus grandes valeurs de
![{\displaystyle \mathrm {V} _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a529b23683bf150c1675f07d7fa23345e7c2d9f)
et de
![{\displaystyle \mathrm {W} _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78c9a2942bba712cd13edccee7704d2f4c9fb835)
dans la partie supérieure de l’orbite, on aura, pour les valeurs des intégrales
![{\displaystyle \int \mathrm {V} _{1}d\Phi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/45ad73c3ac178a4bfa369fb90efb441770f6e2d3)
et
![{\displaystyle \int \mathrm {W} _{1}d\Phi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf66d16b9c4a3474101b9eef66a9b030c97ee01a)
dans cette partie, les limites
![{\displaystyle \pm (\mathrm {V} _{1})(\Phi ')}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d82f64991b7df0f41ac43bfbc4b2ef2748d07244)
et
![{\displaystyle \pm (\mathrm {W} _{1})(\Phi ').}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7a997e01f1e555991f40e03827a0aa813c054941)
Or, sans chercherles valeurs exactes de
![{\displaystyle (\mathrm {V} _{1})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b481d60a259e7f24f09a39930c27fccfedcc2dcc)
et
![{\displaystyle (\mathrm {W} _{1}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78f5a717fed85be74534c2593737cf071a732b1c)
il suffira de considérer que, les plus grandes valeurs de
![{\displaystyle \mathrm {V} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/664209da7650f00b3507efe25f89aeff9783146c)
et
![{\displaystyle \mathrm {W} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe06244bc27928da09da0f0d8202aed23a10f480)
étant égales à
![{\displaystyle t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/65658b7b223af9e1acc877d848888ecdb4466560)
(abstraction faite du signe), et les plus grandes valeurs de
![{\displaystyle \sin \mathrm {N} \upsilon }](https://wikimedia.org/api/rest_v1/media/math/render/svg/41185e59e3beb80f3768e1627f7eadd0d11b1f7c)
et
![{\displaystyle \cos \mathrm {N} \upsilon }](https://wikimedia.org/api/rest_v1/media/math/render/svg/630f44758b4d8285e46011ba3e62092e1ee479c7)
étant
![{\displaystyle 1,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9cc5fd8163a83100c5330622e9e317fa4e872403)
les plus grandes valeurs de
![{\displaystyle \mathrm {V} _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a529b23683bf150c1675f07d7fa23345e7c2d9f)
et de
![{\displaystyle \mathrm {W} _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78c9a2942bba712cd13edccee7704d2f4c9fb835)
ne pourront jamais être plus grandes que
![{\displaystyle t+{\frac {1}{\mathrm {N} }}{\sqrt {\frac {\alpha ^{3}}{8(1+m)}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1ce8c1e6c16ae00d2a9fffabddb8add0692e20f1)
en sorte qu’on pourra prendre dans les limites précédentes
![{\displaystyle (\mathrm {V} _{1})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b481d60a259e7f24f09a39930c27fccfedcc2dcc)
et
![{\displaystyle (\mathrm {W} _{1})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c78afb15f16241465d176932a60a4e801c7af7af)
égales à
![{\displaystyle (t)+{\frac {1}{\mathrm {N} }}{\sqrt {\frac {\alpha ^{3}}{8(1+m)}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/912426ef3089f81edf0e16b1e4d1dd54d354f3ea)
Et comme la quantité
![{\displaystyle \Phi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bb68da4ee8f61f705569b558e842156ea7548bca)
est nécessairement moindre que
![{\displaystyle \Phi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/aed80a2011a3912b028ba32a52dfa57165455f24)
dans la partie supérieure de l’orbite, il est clair que ces nouvelles limites seront plus approchées que les premières, et qu’ainsi l’erreur qu’on commettrait en négligeant les intégrales
![{\displaystyle \int \mathrm {V} _{1}d\Phi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/45ad73c3ac178a4bfa369fb90efb441770f6e2d3)
et
![{\displaystyle \int \mathrm {W} _{1}d\Phi '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf66d16b9c4a3474101b9eef66a9b030c97ee01a)
sera beaucoup moindre que celle qui résulterait de l’omission des premières intégrales
![{\displaystyle \int \mathrm {V} d\Phi ,\int \mathrm {W} d\Phi \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/09f5102e7273e98cf105e2d6c18f401e6bb361d7)
et ainsi de suite.
59. De ce que nous venons de démontrer depuis le no 48 jusqu’ici, il est aisé de conclure que les perturbations que la comète doit éprouver dans la partie supérieure de son orbite peuvent être déterminées analytiquement sans avoir recours aux quadratures mécaniques, sinon par des formules rigoureuses, du moins par des formules très-approchées et dont on peut pousser l’approximation aussi loin que l’on veut. Si ce Mémoire n’était peut-être pas déjà trop long, je présenterais ici ces formules toutes développées, en sorte qu’il n’y eût plus que les substitutions numériques à faire ; mais, comme cela ne demande plus qu’un travail mécanique et de calcul, nous croyons pouvoir nous en dispenser et nous contenter d’avoir exposé les principes de cette analyse avec tout le détail et la clarté nécessaires.
Nous allons donner maintenant une idée de la manière dont on doit faire usage des Théories précédentes, en montrant comment on doit les appliquer à la comète des années 1532 et 1661, que les Astronomes attendent vers 1789 ou 1790.
60. La comète de l’année 1532 a été observée par Appien, et calculée par Halley ; ses éléments sont (la distance moyenne du
étant
![{\displaystyle {\begin{array}{ll}\mathrm {Temps\ moyen\ du\ p{\acute {e}}rih{\acute {e}}lie,\ {\grave {a}}\ Paris,\ 19\ octobre} &\quad 22^{\text{h}}29^{\text{m}}0^{\text{s}}\\\mathrm {Longitude\ du\ p{\acute {e}}rih{\acute {e}}lie} \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\,3^{\text{s}}21^{\circ }7'0''\\\mathrm {Distance\ p{\acute {e}}rih{\acute {e}}lie} \,\ \ \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\,0{,}50910\\\mathrm {Longitude\ du\ noeud\ ascendant} \,\ \ \ldots \ldots \ldots \ldots \ldots &\,2^{\text{s}}20^{\circ }37'0''\\\mathrm {Inclinaison\ de\ l'orbite} \ \ \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\quad 32^{\circ }36'0''\\\mathrm {Sens\ du\ mouvement} \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\mathrm {direct} .\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d0b5800cb966517f9023f01b8ffe2faf61a25255)
Celle de l’année 1661 a été observée par Hevelius et calculée par Halley ; ses éléments sont
![{\displaystyle {\begin{array}{ll}\mathrm {Temps\ moyen\ du\ p{\acute {e}}rih{\acute {e}}lie,\ {\grave {a}}\ Paris,\ 26\ janvier} &\quad 23^{\text{h}}42^{\text{m}}0^{\text{s}}\\\mathrm {Longitude\ du\ p{\acute {e}}rih{\acute {e}}lie} \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\,3^{\text{s}}25^{\circ }58'40''\\\mathrm {Distance\ p{\acute {e}}rih{\acute {e}}lie} \,\ \ \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\,0{,}44851\\\mathrm {Longitude\ du\ noeud\ ascendant} \,\ \ \ldots \ldots \ldots \ldots \ldots &\,2^{\text{s}}22^{\circ }30'30''\\\mathrm {Inclinaison\ de\ l'orbite} \ \ \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\quad 32^{\circ }35'50''\\\mathrm {Sens\ du\ mouvement} \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\mathrm {direct} .\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/12e0ba421d85a92615ad840e2221ea4b5dec1df9)
Comme les éléments de ces deux comètes sont à très-peu près les mêmes, on est fondé à prendre ces astres pour une même comète, dont la révolution serait d’environ
ans, et qui devrait, par conséquent, reparaître en 1789. Dans cette hypothèse, on peut attribuer les différences qui se trouvent entre les él’éments de 1532 et de 1661 en partie à l’inexactitude des observations, du moins de celles de 1532, et en partie à l’effet des perturbations que la comète a dû éprouver pendant la révolution de 1532 à 1661 par l’action des planètes ; et l’on ne saurait fixer au juste le retour de cette comète qu’en calculant d’avance l’effet des perturbations qu’elle doit éprouver dans la révolution de 1661 à 1789.
61. Comme les observations de 1661 ont été faites par Hevelius, on peut les prendre pour exactes, ainsi que les éléments que Halley en a déduits. On peut supposer de plus que ces éléments soient ceux de l’orbite non altérée, puisque, en faisant abstraction des perturbations qui ont précédé et suivi l’apparition de cette comète en 1661, elle aurait dû se mouvoir toujours dans le même plan et avoir le périhélie placé dans le même lieu du ciel.
Nous prendrons donc, pour plus de simplicité, le temps du passage par le périhélie en 1661, c’est-à-dire, le 21 janvier
temps moyen de Paris, pour l’époque du temps
en supposant
positif après cette époque et négatif avant elle, et en se souvenant que
exprime l’angle du mouvement moyen du Soleil (no 20).
Nous prendrons de plus le plan qui coupe l’écliptique
et sous un angle égal à
(cet angle doit être du côté du Nord et sur la partie de l’écliptique comprise entre
et
pour le plan fixe des coordonnées
et
et nous prendrons l’axe des
dans la ligne menée du Soleil au point de ce plan qui répond à
de longitude comptée depuis le lieu de l’équinoxe en 1661. Nous rapporterons ensuite à ce même plan et à ce même axe les lieux des planètes perturbatrices au moyen des coordonnées rectangles
Ces déterminations s’accordent avec les suppositions des nos 2 et 25.
62. Cela posé, soient, comme dans la Section deuxième,
le grand axe de l’orbite non altérée,
le paramètre de ce grand axe, et
l’excentricité de l’orbite, égale à
la distance périhélie sera
![{\displaystyle {\frac {a-ae}{2}}={\frac {a}{2}}(1-e)={\frac {a\left(1-e^{2}\right)}{2(1+e)}}={\frac {2h}{1+e}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1b47bcbdd2348b29a77e82048672f3b82d6fb0c4)
Or, par les observations de 1661, on a conclu la distance périhélie
(la distance moyenne du Soleil à la Terre étant
) ; on aura donc
![{\displaystyle {\frac {2h}{1+e}}=0{,}44851.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/23c4283c8ac1f3dc3dd016bf456984d880be74a8)
Mais j’observe que, comme les éléments de 1661 ont été calculés dans l’hypothèse de l’orbite parabolique, il paraît naturel d’adopter aussi cette hypothèse dans la détermination de
or, dans la parabole, on a
donc
par conséquent
Au reste, nous verrons ci-après que cette valeur de
ne serait tout au plus diminuée que d’un centième, si l’on voulait la déterminer dans l’hypothèse elliptique (no 64).
63. Il faut déterminer maintenant le grand axe
ce qu’on fera par le principe connu que, dans les orbites elliptiques décrites par une même force tendant au foyer et variant dans la raison inverse des carrés des distances, les temps des révolutions sont comme les racines carrées des cubes des moyennes distances. Or l’intervalle entre le passage par le périhélie en 1532 et le passage par le périhélie en 1661 est de
années
mais il faut remarquer que, comme en 1582 on a retranché
jours, il faut aussi les retrancher du nombre précédent, ce, qui réduira le vrai intervalle entre les deux passages par le périhélie à
années
parmi lesquelles années il y en a
de bissextiles.
Réduisant ce temps en jours et en décimales de jour, on aura donc, pour l’intervalle dont il s’agit,
Si les deux périhélies étaient placés dans le même lieu du ciel, il est clair que l’intervalle qu’on vient de trouver serait en même temps la durée de la révolution de la comète ; mais, comme les lieux des deux périhélies diffèrent\nu peu entre eux, il faut défalquer de l’intervalle trouvé le temps que la comète a mis à aller du lieu du périhélie de 1532 à celui du périhélie de 1661.
Pour cela, j’observe que le périhélie de 1661 est plus avancé en longitude que celui de 1532 de
mais l’équinoxe a reculé, dans l’espace de
de
donc, retranchant cette quantité de la précédente, on aura
pour le vrai espace dont le périhélie de 1661 était plus avancé par rapport aux étoiles fixes que celui de 1532 ; donc la comète, après avoir atteint en 1532 son périhélie, a dû parcourir encore autour du Soleil un angle de
pour arriver au lieu du périhélie de 1661, parce que le mouvement de cette comète se fait suivant l’ordre des signes ; ce serait le contraire si la comète était rétrograde.
Il faut donc chercher le temps qui répond à l’anomalie vraie
dans une parabole dont la distance périhélie est
Or, dans la Table générale du mouvement des comètes (cette Table, calculée d’abord par Halley, rendue ensuite plus commode par l’abbé de la Caille, a été étendue davantage dans le Recueil des Tables publié par l’Académie de Berlin, t. III, p. 2 et suiv.), on trouve que pour l’orbite, dont la distance périhélie serait
ce temps serait de
il faut donc multiplier ce nombre par la racine carrée du cube de
et l’on aura
pour le nombre qu’il faudra retrancher du nombre de jours trouvé plus haut, pour avoir la durée de la révolution de la comète par rapport aux étoiles fixes, laquelle durée sera donc
Or la durée de la révolution périodique de la Terre, c’est-à-dire, l’année sidérale, est de
ou, en décimales de jour, de
Donc, puisque la distance moyenne de la Terre au Soleil est prise pour l’unité, et que la distance moyenne de la comète est
on fera cette proportion
![{\displaystyle 365{,}25636:46840{,}35515=1:\left({\frac {a}{2}}\right)^{\frac {3}{2}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/21f07a4633c689896a94b8a89d9c97914f77f1cf)
d’où l’on tire
![{\displaystyle {\frac {a}{2}}=\left({\frac {46840{,}35515}{365{,}25636}}\right)^{\frac {3}{2}}=25{,}43013.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b37fedcab35ddf3ef6f1daa63aa07cc2525fe80)
C’est la distance moyenne ou le demi-grand axe de l’orbite elliptique de la comète.
64. Il est aisé de conclure de l’équation précédente que, si le temps périodique de la comète était plus long ou plus court d’un petit nombre
d’années sidérales, la distance moyenne
\frac{a}{2}
serait augmentée ou diminuée à très-peu près de la quantité
![{\displaystyle {\frac {2n}{3{\sqrt {\cfrac {a}{2}}}}}=0{,}13220n\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d9c9db0d95829d74f218a6c67d21e68d0ba4c576)
de sorte qu’il faudrait que
![{\displaystyle n}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a601995d55609f2d9f5e233e36fbe9ea26011b3b)
fût plus grand que
![{\displaystyle 7,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a1341efd8b44bd78556844e2ad48200b6da40db0)
pour que la distance moyenne fût changée d’une unité.
Or, quoique les observations de 1532 faites par Appien ne soient peut-être pas tout à fait exactes, cependant, comme la comète dont il s’agit n’a été observée que pendant un mois, dans lequel temps elle a passé par le périhélie, il est visible qu’on ne saurait admettre une erreur de
jours dans le passage au périhélie, et qu’ainsi, à cet égard, on est assuré que la valeur de
est exacte à
près.
Mais il y a une autre source d’erreur qui est bien plus considérable je veux parler de l’effet des perturbations que la comète a dû éprouver dans la période de 1532 à 1661, et qui ont pu allonger ou diminuer cette période de quelques années. En effet, la détermination précédente du grand axe étant fondée sur l’hypothèse que la comète a décrit une orbite régulière autour du Soleil en vertu de la seule attraction de cet astre, cette détermination cesse d’être exacte dès qu’on admet l’action des planètes sur la comète ; dans ce cas, il est clair qu’on ne saurait chercher le grand axe de la véritable orbite décrite par la comète, puisque cette orbite n’est plus une ellipse ; mais on doit chercher plutôt le grand axe de l’orbite que la comète aurait décrite sans les perturbations, et que nous avons nommée, dans le cours de ce Mémoire, orbite non altérée ; et pour cela il est visible qu’on ne doit pas employer la durée observée de la révolution, mais cette durée corrigée de l’effet des perturbations.
Supposons que cet effet consiste à allonger ou à raccourcir le temps périodique de l’orbite non altérée d’un petit nombre
d’années sidérales pendant la période de 1532 à 1661 il est clair que ce temps périodique sera plus court ou plus long de
années que la durée observée de la révolution de la comète ; par conséquent la distance moyenne de l’orbite non altérée sera à très-peu près
Or on ne peut déterminer la valeur de
que par le calcul même des perturbations, calcul dans lequel la quantité
entre comme élément ; mais, comme la valeur de
ne peut être que de quelques unités, il sera permis de prendre pour la valeur de
la quantité
qui aurait lieu sans les perturbations, du moins dans le calcul de ces perturbations. L’erreur qu’on pourra commettre par cette supposition ne sera que de l’ordre des carrés des forces perturbatrices, quantités que nous avons toujours supposées, qu’on néglige dans la Théorie des perturbations des comètes.
65. En faisant donc
et prenant pour
la valeur déterminée ci-dessus (no 62), savoir
on trouvera d’abord l’excentricité (no 25)
![{\displaystyle e={\sqrt {1-{\frac {4h}{a}}}}=0{,}98222=f.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6e3d35da2a95ce6a770d8f7b5e1b11aeb03e2fdc)
Employant cette valeur de
dans l’équation
![{\displaystyle {\frac {2h}{1+e}}=0{,}44851}](https://wikimedia.org/api/rest_v1/media/math/render/svg/64f2e401e263902c2a70c17450160866882286fa)
du numéro cité, on trouvera
c’est la valeur de
dans la supposition que la distance périhélie, déduite des observations, soit la véritable distance périhélie dans l’ellipse ; et l’on voit que cette valeur diffère à peine de
de celle que donne la supposition de l’orbite parabolique ; c’est pourquoi on pourra sans crainte employer la première valeur de
dans le calcul des perturbations.
Comme le demi-petit axe de l’ellipse est
on trouvera pour ce demi-petit axe
Et, si l’on cherche l’angle dont le cosinus sera
on trouvera
c’est la valeur de l’anomalie excentrique qui répond à
degrés d’anomalie vraie, à compter du périhélie, et c’est aussi la valeur de l’anomalie vraie comptée de l’aphélie pour les points de la distance moyenne.
Enfin, comme nous prenons le périhélie de 1661 pour l’époque d’où l’on doit compter le temps
et que nous supposons que, dans ce périhélie, l’orbite troublée coïncide avec l’orbite non troublée (no 60), il s’ensuit qu’on aura non-seulement
mais aussi
et de plus
dans le même périhélie (no 38) ; donc les cinq variables
devront être nulles lorsque
de sorte que, en faisant commencer dans ce point les intégrations des différentielles
et
(no 42), il n’y aura point de constantes à y ajouter.
Nous remarquerions encore que, à cause de
on aura simplement (no 20)
![{\displaystyle \theta =2t{\sqrt {\frac {2(1+m)}{a^{3}}}}\quad \mathrm {et\ par\ cons{\acute {e}}quent} \quad t=\theta {\sqrt {\frac {8(1+m)}{a^{3}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/020a50a1e2f5175645a5cdbc83f50d44edc61559)
et les angles
seront tous nuls à la fois dans le périhélie de 1661 ; ils seront positifs après ce périhélie, et négatifs avant.
À l’égard de la quantité
elle devrait aussi être nulle lorsque
si la valeur de
était exactement égale au grand axe de l’orbite non altérée de la comète ; mais, ayant supposé ci-dessus
on aura (no 63), pour la vraie distance moyenne de l’orbite non altérée,
![{\displaystyle 25{,}43013\mp 0{,}13220n,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/68a9913aedf7b2c32a420a091ab497c24b6f137f)
laquelle doit être, par l’hypothèse, la même que celle de l’orbite troublée dans le périhélie de 1661, où
Or le grand axe de l’orbite troublée étant, en général,
(no 38), on aura donc, lorsque ![{\displaystyle t=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55ff4c2b109c38fe7038da6238ae875f4d37e643)
![{\displaystyle {\frac {a+\delta a}{2}}=25{,}43013\mp 0{,}13220n\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d792fec99324720948cc9d89560ee97689ba2acd)
donc
![{\displaystyle {\frac {\delta a}{2}}=\mp 0{,}13220n.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2432b96031bccf3851a25d3fc3e8e588deff8260)
D’où l’on voit que cette valeur de
dépend du nombre
qui est l’effet des perturbations dans la période précédente.
66. Considérons maintenant le retour de la comète au périhélie ; il est visible que, en faisant abstraction des perturbations, il n’y aura qu’à ajouter l’intervalle trouvé ci-dessus (no 62) de
à l’époque du passage par le périhélie de 1661, pour avoir le temps du retour au périhélie ; et il viendra (en se souvenant que l’année 1700 n’a pas été bissextile) le 26 avril 1789
temps moyen au méridien de Paris. Mais, pour avoir exactement le temps du retour de la comète au périhélie dans l’orbite elliptique dont le demi-grand axe serait
tel qu’on l’a déterminé dans le numéro cité, il faudra retrancher du temps qu’on vient de trouver
c’est-à-dire,
par la raison expliquée dans ce même numéro, ce qui donnera le 25 avril 1789
Cette détermination serait entièrement exacte, même en ayant égard aux perturbations, si les deux révolutions consécutives de 1532 à 1661, et de 1661 à 1789, étaient parfaitement égales, et par conséquent si l’effet des perturbations était le même dans ces deux périodes. Donc, si l’altération de ces deux périodes n’est pas la même, il est clair qu’il ne faudra qu’ajouter, au temps déterminé ci-dessus, l’excès de l’altération de la seconde période sur l’altération de la première.
Or nous avons donné, dans le no 41, la formule qui exprime, en général, l’altération de la révolution périodique de la comète ; appliquant donc ici cette formule, et marquant par un, deux, trois traits les quantités qui répondent aux trois périhélies consécutifs de 1532, 1661, 1789, on aura cette quantité, dans laquelle j’ai substitué au lieu
sa valeur
(no 38),
![{\displaystyle {\frac {3\left(t'''\delta a'''-2t''\delta a''+t'\delta a'\right)}{2a}}-{\sqrt {\frac {a^{3}}{8(1+m)}}}\left(\delta i'''-2\delta i''+\delta i'\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6990f9f198677b44045861aa28f88bfd9dd10195)
![{\displaystyle -\left({\frac {2h}{1+c}}\right)^{2}{\frac {\left(\delta g'''-2\delta g''+\delta g'\right)}{c{\sqrt {2h(1+m)}}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9d200859f1b3c7c385c08686118dc1cee40d3fda)
c’est la correction du temps, c’est-à-dire, le temps qu’il faudra ajouter au 25 avril 1789
pour avoir l’instant du passage de la comète par la ligne du périhélie de 1661, dont la longitude était alors de
et sera, en 1789 (à cause de la précession des équinoxes), de ![{\displaystyle 3^{\text{s}}27^{\circ }46'16''.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/391f28384b4184f9a5bee49b307df715b0ed28ea)
Il est bon de remarquer que la dernière partie de la quantité précédente, celle qui contient les quantités
dépend uniquement du déplacement du périhélie, comme on peut le voir par le no 41 ; de sorte qu’en rejetant cette partie la quantité restante sera la correction du temps pour le passage de la comète par le vrai périhélie de 1789 ; mais il faudra alors ajouter cette correction au 26 avril 1789
temps du passage par le périhélie, dans le cas où la révolution anomalistique de 1661 à 1789 serait égale à celle de 1532 à 1661.
Pour réduire ces quantités en temps, on se souviendra que nous exprimons le temps par le mouvement moyen du Soleil (no 20) ; de sorte que, si la quantité à réduire en temps est exprimée en angles, en la divisant par
degrés, on aura des années sidérales de
et, si elle est exprimée en nombres absolus (la moyenne distance du Soleil étant l’unité), il faudra la diviser par le rapport de la circonférence au rayon, c’est-à-dire, par
pour la réduire de même en années sidérales.
67. La formule précédente est générale ; mais, dans notre cas, on aura, par ce qu’on a établi dans le no 64,
![{\displaystyle t''=0,\quad \delta i''=0,\quad \delta g''=0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ef642172c53d2965e9c82617a2301b66225f570d)
de plus, à cause de
(
étant l’anomalie moyenne de la comète), on aura
![{\displaystyle t'''=360^{\circ }{\sqrt {\frac {a^{3}}{8(1+m)}}},\quad t'=-360^{\circ }{\sqrt {\frac {a^{3}}{8(1+m)}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b2d4276ad205928dd0960e58249108175dcdb53a)
de sorte que la première partie de la formule dont il s’agit se réduira à
![{\displaystyle {\frac {3.360^{\circ }}{4}}{\sqrt {\frac {a}{2(1+m)}}}\left(\delta a'''-\delta a'\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/75dba2b09035f77b236f8b22257911958984d8d3)
où
sera l’intégrale totale de la différentielle
pour la période entière de 1661 à 1789, en faisant
positif, et où
sera de même l’intégrale totale de
pour la période de 1661 à 1532, en faisant
négatif de sorte que, comme il n’y a que la différence de ces deux intégrales qui entre en ligne de compte, il n’y aura point de constantes à y ajou-
ter, et l’on pourra prendre chaque intégrale en sorte qu’elle commence au périhélie de 1661, où
![{\displaystyle t=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/43469ec032d858feae5aa87029e22eaaf0109e9c)
.
Les deux autres parties de la même formule deviendront
![{\displaystyle -{\sqrt {\frac {a^{3}}{8(1+m)}}}\left(\delta i'''+\delta i'\right)-\left({\frac {2h}{1+e}}\right)^{2}{\frac {\delta g'''+\delta g'}{e{\sqrt {2h(1+m)}}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/06d8ef34464148b26b2ab91d8dcb43b8c4c522ec)
où (à cause de
et
) il faudra prendre pour
et
les intégrales des quantités
et
depuis le périhélie de 1661 jusqu’à celui de 1789, et pour
et
les intégrales des mêmes quantités, mais depuis le périhélie de 1661 jusqu’à celui de 1532, en supposant
négatif.
L’altération du temps périodique est celle qu’il est le plus important de déterminer dans la Théorie des perturbations des comètes.
Quant aux altérations des autres éléments de l’orbite, on les déterminera directement par l’intégration des quantités
(nos 38, 42 et suivants), en faisant commencer les intégrales au périhélie de 1661, et les étendant jusqu’au périhélie de 1789 ou de 1532, suivant qu’on voudra déterminer ces altérations pour la dernière période de la comète ou pour la période précédente.
68. Voilà toutes les données et les formules nécessaires pour calculer les perturbations causées à l’orbite de la comète dont il s’agit par l’action des planètes. Or, parmi toutes les planètes, il n’y a que Jupiter et Saturne dont l’action sur la comète puisse être sensible, tant parce que les masses des autres planètes sont trop petites, que parce qu’elles sont trop proches du Soleil. Ainsi l’on prendra successivement Jupiter et Saturne pour la planète perturbatrice dont on a supposé la masse
et dont le rayon vecteur est
et les coordonnées rectangles
Les Tables astronomiques de Halley donneront toutes les valeurs des quantités qui dépendent des lieux de ces planètes dans un temps quelconque ; et nous ne croyons pas qu’il soit nécessaire d’entrer là-dessus dans aucun détail.
Comme la distance de Jupiter au Soleil est environ
et celle de Saturne environ
il est clair que, si l’on fait commencer la partie supérieure de l’orbite de la comète aux points de la moyenne-distance, c’est-à-dire, aux extrémités du petit axe, alors
sera toujours beaucoup plus grand que
et la méthode des nos 48 et suivants aura toute l’exactitude qu’on peut désirer, en négligeant même tout à fait les termes qui dépendent du mouvement moyen de la planète perturbatrice dans les formules des quantités différentielles
et en ne tenant compte que des termes indépendants de l’angle
que nous avons vus être toujours intégrables (nos 52 et 56).
Il y aura encore un autre avantage à prendre ainsi la moitié supérieure de l’orbite pour ce que nous appelons la partie, supérieure, car on aura alors, pour le commencement de cette partie,
et, pour la fin,
de sorte que, à cause de
![{\displaystyle x={\frac {a}{2}}(\cos u-e),\ \ y={\sqrt {ah}}\sin u,\ \ r={\frac {a}{2}}(1-e\cos u),\ \ dt={\sqrt {\frac {a}{2(1+m)}}}rdu}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b26fd893851b028dd5d3b132cdd3a657b0c5ee89)
(
étant l’anomalie excentrique comptée depuis le périhélie), on aura, pour le commencement de la partie supérieure,
![{\displaystyle x=-{\frac {ae}{2}},\quad y=\pm {\sqrt {ah}},\quad r={\frac {a}{2}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69d86e8b3810b330858d0262d9e327ea2dcb247b)
![{\displaystyle {\frac {dx}{dt}}=\mp {\sqrt {\frac {2(1+m)}{a}}},\quad {\frac {dy}{dt}}=0,\quad {\frac {dr}{dt}}=\pm e{\sqrt {\frac {2(1+m)}{a}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1ef3748a168f3aa3dca8af386530ace1885b3a6f)
et pour la fin
![{\displaystyle x=-{\frac {ae}{2}},\quad y=\mp {\sqrt {ah}},\quad r={\frac {a}{2}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a0ee7475b3023cc9566727e3007baf9a24af361e)
![{\displaystyle {\frac {dx}{dt}}=\pm {\sqrt {\frac {2(1+m)}{a}}},\quad {\frac {dy}{dt}}=0,\quad {\frac {dr}{dt}}=\mp e{\sqrt {\frac {2(1+m)}{a}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d3d5deec4f2a1dc65eda1ec8ea250e0ef0fe24b0)
(les signes supérieurs étant pour le cas où l’on prendra les angles
et
positifs, c’est-à-dire, pour la période qui suit le périhélie de 1661, et les signes inférieurs étant pour le cas où
et
seront négatifs, c’est-à-dire, pour la période qui précède ce périhélie), ce qui simplifiera beaucoup
les valeurs des quantités
![{\displaystyle \mathrm {H',H'',A',A''} ,\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc84668d123ee2f53482e7f6dd908baba167cbd4)
c’est-à-dire, les valeurs de
![{\displaystyle \mathrm {H,A} ,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f484afb77ca23a546992784adf8780c118e344e)
pour le commencement et pour la fin de la partie supérieure de l’orbite (
no 48).
Quant à la masse
de la comète, que nous avons conservée, pour plus d’exactitude et de généralité, dans les formules de ce Mémoire, elle est inconnue, et rien ne saurait conduire à la déterminer ; mais il est naturel de la supposer très-petite vis-à-vis de la masse du Soleil ; de sorte qu’on pourra négliger partout la quantité
vis-à-vis de l’unité.