Mécanique analytique/Partie 2/Section 7

La bibliothèque libre.
Gauthier-Villars et Fils (Œuvres de Lagrange. Tome XIIp. 1-166).

SECONDE PARTIE.

LA DYNAMIQUE (Suite).

Séparateur


SECTION SEPTIÈME.

Sur le mouvement d’un système de corps libres, regardés comme des points, et animés par des forces d’attraction..


On peut ranger en trois classes tous les systèmes de corps qui agissent les uns sur les autres, et dont on peut déterminer le mouvement par les lois de la Mécanique ; car leur action mutuelle ne peut s’exercer que de trois manières différentes qui nous soient connues ou par des forces d’attraction, lorsque les corps sont isolés, ou par des liens qui les unissent, ou enfin par la collision immédiate. Notre système planétaire appartient à la première classe, et par cette raison les problèmes qui s’y rapportent doivent tenir le premier rang parmi tous les problèmes de la Dynamique. Nous allons en faire l’objet de cette Section.

Quoique dans les systèmes de cette classe, où les corps sont supposés se mouvoir librement, il soit très facile de trouver les équations de leur mouvement, puisqu’il ne s’agit que de réduire toutes les forces à trois directions perpendiculaires entre elles et d’égaler, par le principe des forces accélératrices, la force suivant chacune de ces directions à l’élément de la vitesse relative à la même direction, divisé par l’élément du temps ; néanmoins, l’usage des formules données dans la Section IV est toujours préférable, parce qu’elles fournissent directement, et sans aucune décomposition préalable de forces, les équations différentielles les plus simples, quelles que soient les coordonnées qu’on emploie pour déterminer la position des corps, même lorsque les corps, au lieu d’être tout à fait libres, sont contraints de se mouvoir sur des surfaces ou des lignes données.

Nous commencerons par rappeler les formules dont nous ferons usage.


1. Soient les masses des différents corps regardés comme des points ; les coordonnées rectangles du corps celles du corps et ainsi de suite, ces coordonnées étant toutes rapportées aux mêmes axes fixes de l’espace. On fera

Et si, à la place des coordonnées rectangles on veut employer d’autres coordonnées quelconques il n’y aura qu’à substituer les valeurs de en dans la formule de même, on substituera dans les valeurs de en si l’on veut transformer les coordonnées rectangles en et ainsi de suite. De cette manière, la quantité deviendra une fonction des variables et de leurs différences premières.

Soient maintenant les forces avec lesquelles chaque point de la masse tend vers des centres fixes ou non, dont les distances soient lesquelles, étant données en deviendront aussi des fonctions de on fera

soit que soit une différentielle complète ou non ; et, dénotant par les mêmes lettres marquées d’un trait, de deux traits, les quantités analogues relatives aux corps on fera de plus

Si, outre ces forces dirigées vers des centres donnés, il y avait des forces d’attraction mutuelle entre toutes les molécules des corps et en nommant la distance de ces corps, regardés comme des points, et la force d’attraction, dépendante de la distance ou non, il faudrait ajouter à le terme et ainsi pour tous les autres corps qui s’attireraient mutuellement.

Or, les corps étant supposés libres, les coordonnées qui déterminent leur position dans l’espace sont indépendantes, et chacune d’elles, comme, donnera une équation de la forme


2. Lorsque les quantités sont des différentielles complètes, ce qui a toujours lieu dans le cas où les forces sont proportionnelles à des fonctions quelconques de leurs distances aux centres d’attraction, lequel est celui de la nature, il sera plus simple de prendre d’abord les intégrales lesquelles seront

et la quantité deviendra

laquelle étant réduite en fonction des variables il sera aisé d’en déduire par la différentiation les différences partielles Dans ce cas, si les fonctions et ne renferment point

le temps fini on aura toujours l’intégrale

étant une constante arbitraire, laquelle renferme le principe des forces vives.

CHAPITRE PREMIER.

Du mouvement d’un corps regardé comme un point et attiré vers un centre fixe par des forces proportionnelles à une fonction de la distance, et en particulier du mouvement des planètes et des comètes autour du soleil.

3. Lorsqu’on ne considère que le mouvement d’un corps isolé, on peut supposer sa masse égale à l’unité, et l’on aura simplement

et

Dans ce cas, quelles que soient les trois coordonnées qui déterminent le lieu du corps dans l’espace, comme elles sont indépendantes, elles donneront trois équations différentielles de la forme

auxquelles on pourra joindre l’équation du premier ordre

qui tiendra lieu de l’une d’entre elles.

Si le mouvement se faisait dans un milieu résistant, en désignant la résistance par il n’y aurait qu’à ajouter à la valeur de les termes (Sect. II, art. 8)

mais l’équation n’aurait plus lieu.


4. Supposons que le corps soit attiré vers un centre fixe par une force fonction de la distance du corps au centre, on aura simplement

Prenons la distance pour l’une des coordonnées du corps, et prenons, pour les deux autres, l’angle que le rayon vecteur fait avec le plan des et l’angle que la projection de sur ce plan fait avec l’axe des en plaçant l’origine des coordonnées rectangles dans le centre des rayons de manière que l’on ait

on trouve facilement

et de là

On aura donc ces trois équations différentielles relatives à

lesquelles deviennent

et l’équation

donnera tout de suite cette première intégrale

dans laquelle est une constante arbitraire.


5. La dernière des trois équations différentielles est intégrale d’elle-même ; son intégrale est

étant une constante arbitraire ; et la seconde devient intégrable en y substituant pour sa valeur tirée de celle-ci, et en la multipliant par l’intégrale est

étant une nouvelle constante arbitraire.

Je remarque d’abord sur cette intégrale que, si l’on suppose que et soient nuls à la fois dans un instant, ils seront toujours nécessairement nuls ; car, en faisant pour un instant et la dernière équation donne et elle devient, par la substitution de au lieu de

laquelle ne peut avoir lieu qu’en faisant

La supposition dont il s’agit revient à faire en sorte que le corps se meuve dans un instant dans le plan des ce qui est toujours possible, puisque la position de ce plan est arbitraire ; alors le corps continuera de se mouvoir dans le même plan et décrira nécessairement une orbite plane, c’est-à-dire une ligne à simple courbure. C’est ce

qu’on peut aussi démontrer directement par l’intégration de la même équation. Car, en y substituant pour sa valeur tirée de la première intégrale, elle devient

Soit, lorsque on aura

et la dernière équation se changera en

d’où l’on tire

équation séparée dont l’intégrale est

ou bien

étant la valeur de lorsque

Cette équation fait voir que et sont les deux côtés d’un triangle sphérique rectangle dans lequel est l’angle opposé au côté Ainsi, puisque l’arc est pris sur le plan des et que l’arc est toujours perpendiculaire à ce même plan, il s’ensuit que l’arc qui joint ces deux-ci, et qui forme l’hypoténuse du triangle, fera avec la base l’angle constant par conséquent, cet arc passera par les extrémités de tous les arcs, et tous les rayons se trouveront dans le plan du même arc, lequel sera ainsi le plan de l’orbite du corps, dont l’inclinaison sur le plan des sera l’angle constant et dont l’intersection avec ce même plan fera avec l’axe des l’angle

Si, pour fixer les idées, on prend le plan des pour l’écliptique, sera la longitude sur l’écliptique, la latitude ; la longitude du nœud de l’orbite et son inclinaison.


6. Prenons maintenant l’intégrale

en y substituant pour sa valeur en trouvée ci-dessus, elle devient

laquelle doit être combinée avec l’autre intégrale

Si l’on y substitue la valeur de tirée de celle-ci et qu’on fasse

on aura l’équation

d’où l’on tire

En intégrant cette équation, on aura l’expression de en et réciproquement celle de en


7. On aura ensuite par l’équation

or, comme le plan des angles est arbitraire, si on le fait coïncider avec le plan de l’orbite, en faisant on aura aussi (art. 5), par conséquent et, dans ce cas, l’angle sera celui que le

rayon décrit dans le plan de l’orbite. Donc, si l’on désigne en général cet angle par on aura

et, substituant la valeur de en

équation dont l’intégrale donnera la valeur de en et réciproquement celle de en

Ensuite on aura en par l’équation

laquelle, en substituant pour sa valeur tirée de l’équation

trouvée plus haut, devient

d’où l’on tire par l’intégration

étant une constante arbitraire ; et de là

équation qui indique que est l’hypoténuse du même triangle sphérique rectangle dont la base est et l’angle adjacent (art. 5), et dont le côté opposé à est

On voit par là que est l’angle décrit par le rayon dans le plan de l’orbite, et dont l’origine est à la ligne d’intersection de ce plan avec celui des  ; que est l’angle décrit par la projection de ce rayon sur le même plan, et que est l’inclinaison du plan de l’orbite sur le plan fixe des

8. Le problème est donc résolu, puisqu’il ne dépend plus que de l’intégration des deux équations séparées entre et les six constantes arbitraires nécessaires pour l’intégration complète des trois équations différentielles en et seront et les deux que l’intégration introduira dans les valeurs de et de

Dans la solution que nous venons de donner, nous avons pris pour coordonnées le rayon vecteur avec les deux angles de longitude et de latitude, pour nous conformer à l’usage des astronomes ; aussi cette solution a-t-elle l’avantage d’offrir directement la plupart des théorèmes que l’on ne trouve ordinairement que par la Trigonométrie sphérique. Mais, en l’envisageant du côté analytique, elle est moins simple que si l’on avait conservé les coordonnées rectangles primitives ; c’est ce qu’il est bon de faire voir, d’autant qu’il en résultera de nouvelles formules qui pourront être utiles par la suite.

9. En prenant pour les trois variables indépendantes, les formules générales de l’article 3 donnent tout de suite les trois équations différentielles

et l’équation intégrale

En chassant des trois équations différentielles, on a immédiatement

trois équations intégrables et dont les intégrales sont

étant des constantes arbitraires dont la première est la même que celle de l’équation

de l’article 5, parce qu’en effet celle-ci n’est qu’une transformée de l’équation

par la substitution des valeurs de de l’article 4.

Ces trois intégrales répondent à celles que nous avons données pour un système de corps, dans l’article 9 de la Section III, d’où nous aurions pu les emprunter.

10. En ajoutant ensemble les carrés des trois dernières équations et employant cette réduction connue

on a l’équation

laquelle, en y substituant pour sa valeur tirée de la première intégrale, et faisant, pour abréger,

donne

d’où l’on tire tout de suite

comme dans l’article 6.

Les mêmes équations, étant ajoutées ensemble, après avoir multiplié la première par la deuxième par et la troisième par donnent celle-ci

laquelle est à un plan passant par l’origine des coordonnées et fait voir que l’orbite décrite par le corps est une courbe plane décrite autour du centre des forces.

11. Nommons les coordonnées rectangles de cette courbe, l’axe des étant pris dans la ligne d’intersection du plan de la courbe avec celui des nommons de plus, comme dans l’article 5, l’angle formé par ces deux plans, et l’angle que la même ligne d’intersection fait avec l’axe des ces deux quantités et seront constantes, et, par les formules connues de la transformation des coordonnées, on aura

Ces valeurs, étant substituées dans les mêmes équations, donneront celles-ci :

Ajoutant leurs carrés ensemble et extrayant ensuite la racine, on a (art. 6)

de sorte que les valeurs des constantes seront

Or, désignant par comme dans l’article 7, l’angle que le rayon fait avec la ligne d’intersection du plan de l’orbite et du plan fixe des il est clair qu’on aura

et la dernière des équations précédentes deviendra

laquelle donne le théorème connu des secteurs proportionnels au temps

Substituant la valeur de on aura

comme dans l’article cité.

Ainsi le problème est de nouveau réduit à l’intégration des deux équations séparées en et que nous avions déjà trouvées ci-dessus (art. 6 et 7) ; mais cette intégration dépend de l’expression de la force centrale en fonction du rayon

12. On voit, par ces équations, que ce rayon sera le plus grand ou le plus petit, soit relativement au temps soit relativement à l’angle lorsqu’il sera déterminé par l’équation

Supposons qu’en intégrant ces mêmes équations on prenne les intégrales en de leurs seconds membres, de manière qu’elles commencent au point où est un minimum, et que l’angle commence aussi à ce point ; l’angle sera alors celui que le rayon qui passe par le même point fera avec la ligne d’intersection de l’orbite avec le plan fixe (art. 7) ; et cette constante jointe à celle que l’intégration peut ajouter à et aux constantes ou complétera le nombre des six constantes arbitraires que l’intégration des trois équations différentielles en et doit donner.

13. Si maintenant on fait

il est clair que et seront les coordonnées rectangles de la courbe, placées dans son plan et ayant la même origine que le rayon les abscisses étant dirigées vers le point où est un minimum ; et si l’on substitue ces quantités dans les valeurs de et de l’article 11, on aura

Substituons ces valeurs dans celles de du même article et faisons, pour abréger,

on aura

expressions qui ont cet avantage, que les quantités dépendantes du mouvement dans l’orbite sont séparées des quantités qui dépendent uniquement de la position de l’orbite relativement au plan fixe des

Ces expressions de sont conformes à la théorie générale exposée (Sect. II, art. 10), et l’on aurait pu les en déduire immédiatement.

En effet, en considérant tout de suite le mouvement dans l’orbite, on a les coordonnées la troisième étant nulle, lesquelles, ne renfermant que trois constantes arbitraires, peuvent être regardées comme des valeurs particulières des coordonnées générales ensuite on aura celles-ci, par le moyen des coefficients qui renferment les trois autres constantes.

14. Si, au lieu de considérer le mouvement dans l’orbite propre du corps, on rapportait ce mouvement à un plan quelconque, par les trois coordonnées lesquelles ne continssent aussi que trois constantes arbitraires, on aurait alors par la même théorie les expressions générales

et comme on a trouvé (Sect. III, art. 10)

on aurait

Ces valeurs de renfermant les trois arbitraires satisfont d’une manière générale aux six équations de condition données (Part. I, Sect. III, art. 10)

Après avoir donné les formules générales pour le mouvement d’un corps attiré vers un point fixe, il ne reste qu’à les appliquer au mouvement des planètes et des comètes ; c’est l’objet des paragraphes suivants.

§ 1. — Du mouvement des planètes et des comètes autour du Soleil
supposé fixe.

15. Dans le système du monde, la force attractive étant en raison inverse du carré des distances, on fera étant la force attractive d’une planète vers le Soleil, à la distance ce qui donnera

Substituant cette valeur dans l’équation entre et (art. 11), on voit que la quantité sous le signe devient

laquelle peut se mettre sous la forme

alors le second membre de l’équation exprimera la différentielle de l’angle ayant pour cosinus la quantité

de sorte que, intégrant, ajoutant à la constante arbitraire et passant des arcs à leurs cosinus, on aura

On voit que la plus petite valeur de aura lieu lorsque l’angle est nul de sorte que, comme nous avons supposé (art. 12) que l’angle commence au point qui répond au minimum de on aura

Donc, en faisant, pour abréger,

on aura

équation polaire d’une section conique dont est le paramètre, l’excentricité, c’est-à-dire le rapport de la distance des foyers au grand axe, le rayon vecteur partant d’un des foyers, et l’angle qu’il fait avec la partie du grand axe qui répond au sommet le plus proche de ce foyer.

La plus grande et la plus petite valeur de étant et leur demi-somme sera c’est la distance moyenne que nous désignerons par de sorte qu’on aura

et si l’on substitue ici pour et leurs valeurs en et on aura

d’où l’on voit que la constante doit être négative pour que l’orbite soit elliptique ; si elle était nulle, l’axe serait infini, et l’orbite deviendrait parabolique ; mais, si elle était positive, l’axe serait négatif[1], et l’orbite serait hyperbolique. Dans le premier cas, la valeur de l’excentricité sera moindre que l’unité ; elle sera dans le deuxième cas et plus grande que dans le troisième.

Il y a encore une autre hypothèse d’attraction qui donne aussi une orbite elliptique, c’est l’attraction en raison directe des distances ; mais, comme elle n’est point applicable aux planètes, nous ne nous y arrêterons pas ici. On peut voir les Principes de Newton et les Ouvrages où l’on a traduit ses théories en Analyse.

16. Revenons maintenant à l’équation qui donne en (art. 10), et substituons-y à la place de à la place de et à la place de elle deviendra

Faisons ce qui donne

on aura

et, intégrant avec une constante arbitraire

Cette équation donnera en et comme on a en on aura, par la substitution, en

Si l’on fait la même substitution dans l’équation entre et de l’article 11, on aura celle-ci

dont l’intégrale est

Mais on peut avoir la valeur de en sans une nouvelle intégration, par la simple comparaison des valeurs de laquelle donne l’équation

d’où l’on tire, à cause de

et, de là,

On voit, par ces formules, que, lorsque l’angle est augmenté de le rayon revient le même, et que l’angle est aussi augmenté de Ainsi la planète revient au même point, après avoir fait une révolution entière. Or, l’angle augmentant de le temps se trouve augmenté de c’est le temps que la planète emploie pour revenir au même point de son orbite, et qu’on nomme le temps périodique. Ainsi ce temps ne dépend que du grand axe et il est le même que si la planète décrivait un cercle ayant pour rayon la distance moyenne Dans ce cas, on aurait

ainsi le temps serait proportionnel aux angles parcourus. Et si l’on suppose et qu’on prenne la distance moyenne de la Terre pour l’unité des distances, les temps seront représentés par les angles mêmes que la Terre décrirait si elle se mouvait dans un cercle dont la distance moyenne serait le rayon, avec une vitesse égale à l’unité. Le mouvement dans ce cercle est celui que les astronomes appellent mouvement moyen de la Terre ou du Soleil, et auquel ils rapportent communément les mouvements des autres planètes.

17. Lorsque l’orbite est hyperbolique, le grand axe devient négatif et l’angle imaginaire. Pour appliquer les formules précédentes à ce cas, faisons

on aura par les formules connues, étant le nombre dont le logarithme hyperbolique est

et les équations de l’article précédent deviendront

à cause de

18. L’équation

trouvée dans l’article 15, donne, en substituant pour (art. 13),

Substituant pour sa valeur en on aura

et, comme on trouvera

expressions fort simples qu’on pourra substituer dans les expressions générales de du même article.

Ainsi il ne s’agira plus que de substituer la valeur de en tirée de l’équation donnée dans l’article 16, pour avoir les trois coordonnées en fonction du temps.

19. L’angle que nous venons d’introduire à la place de est ce qu’on appelle en Astronomie anomalie excentrique, et qui répond à l’anomalie moyenne et à l’anomalie vraie mais les astronomes ont coutume de compter ces angles depuis le sommet de l’ellipse le plus éloigné du foyer où le Soleil est supposé placé, et qu’on nomme aphélie ou apside supérieure, au lieu que, dans les formules précédentes, ils sont supposés comptés depuis le sommet le plus proche du même foyer, qu’on nomme périhélie ou apside inférieure. Pour les rapporter à l’aphélie, il n’y aurait qu’à y ajouter l’angle de ou, ce qui revient au même, changer le signe de la quantité mais, en prenant l’origine des anomalies au périhélie, on a l’avantage d’avoir des formules également applicables aux planètes, dont l’excentricité est assez petite, et aux comètes, dont l’excentricité est presque égale à l’unité, leur grand axe étant très grand tandis que le paramètre conserve une valeur finie.

20. Il nous reste à déterminer en c’est-à-dire l’anomalie excentrique par l’anomalie moyenne ; c’est le problème connu sous le nom de problème de Kepler, parce qu’il est le premier qui l’ait proposé et qui en ait cherché la solution. Comme l’équation entre et est transcendante, il est impossible d’avoir, en général, la valeur de en par une expression finie ; mais, en supposant l’excentricité fort petite, on peut l’avoir par une série plus ou moins convergente. Pour y parvenir de la manière la plus simple, nous ferons usage de la formule générale que nous avons démontrée ailleurs[2], pour la résolution en série d’une équation quelconque.

Soit une équation de la forme

dénotant une fonction quelconque de on aura réciproquement

En général, si l’on demande la valeur d’une fonction quelconque de désignée par on fera

et l’on aura

21. Pour appliquer cette formule à l’équation de l’article 16, on fera

on aura immédiatement

où il n’y aura plus qu’à exécuter les différentiations indiquées ; mais, pour avoir des expressions plus simples, il conviendra de développer auparavant les puissances des sinus en sinus et cosinus d’angles multiples de

On aura de même

On aura ainsi, par les formules des articles 16 et 17,

22. On pourrait tirer de là la valeur de l’angle par la série qui donne l’angle par la tangente, mais on aurait difficilement, de cette manière, une série dont on pût connaître la loi. Pour obtenir une telle série, il faudra tirer d’abord la valeur de l’angle de l’équation

ce qu’on peut faire d’une manière élégante, en employant les exponentielles imaginaires. On aura ainsi cette transformée, en prenant pour le nombre dont le logarithme hyperbolique est l’unité,

laquelle se réduit à celle-ci

d’où l’on tire, en faisant

ou bien, en supposant

Prenons maintenant les logarithmes des deux membres, on aura, en divisant par

réduisant les logarithmes du second membre en série, et substituant ensuite, à la place des exponentielles imaginaires, les sinus réels qui y

répondent, on aura enfin la série[3]

Il ne s’agira donc plus que de substituer pour sa valeur en Si donc on fait, pour abréger,

on aura

On peut réduire la valeur de à une forme finie, et l’on trouve

Ces formules ont l’avantage de donner la loi des séries, qui n’était pas connue auparavant.

23. Puisqu’en prenant le plan des pour celui de l’écliptique supposé fixe, et supposant l’axe des dirigé vers le premier point d’Aries, l’angle est ce qu’on appelle la longitude de la planète, l’angle est la longitude du nœud, l’angle est la latitude, il est clair que l’angle dont est la projection sur l’écliptique, sera la longitude dans l’orbite comptée du nœud, ou ce qu’on appelle l’argument de la latitude ; et l’équation (art. 7)

qui donne l’angle par pourra, lorsque l’inclinaison est assez petite, se résoudre en série par la méthode des exponentielles imaginaires employée ci-dessus. Il n’y aura qu’à mettre, dans l’expression de en à la place de à la place de et à la place de ce qui donnera

et l’on aura

L’équation qui donne en (art. 5),

pourrait aussi se résoudre de la même manière, mais il en résulterait une série moins élégante. On aurait d’abord l’équation en exponentielles imaginaires, étant le nombre dont le logarithme hyperbolique est

d’où l’on tirerait


et prenant les logarithmes

enfin, en développant les puissances des exponentielles imaginaires et y substituant les sinus qui y répondent, on aurait

Les séries que nous venons de donner ne sont convergentes qu’à raison de la petitesse de l’excentricité [4] ou de l’inclinaison et ne sont, par conséquent, applicables qu’aux orbites elliptiques peu différentes du cercle et peu inclinées, telles que celles des planètes et de leurs satellites ; il n’y aurait d’exception que pour Pallas, une des quatre nouvelles petites planètes, dont l’inclinaison sur l’écliptique est d’environ ce qui donne pour une fraction encore assez petite ; de sorte que la série de la valeur de en sera très convergente, mais la série de en le sera beaucoup moins.

24. Après le cas où l’excentricité est très petite, le problème de Kepler est encore résoluble analytiquement dans le cas où l’excentricité est peu différente de l’unité, et qui est celui des orbites presque paraboliques, comme celles des comètes. Dans ce cas, le demi grand axe a devient très grand, et l’équation de l’article 15

dans laquelle est le demi-paramètre, donne

L’équation entre et (art. 16), étant mise sous la forme

fait voir que, lorsque est très grand, devient très petit, de sorte qu’on peut developper en

En faisant ces substitutions dans l’équation précédente, on aura

où l’on voit que la quantité est de l’ordre de Si donc on fait

et qu’on ne pousse l’approximation que jusqu’aux termes de l’ordre de on aura

On trouvera par les mêmes réductions

Soit la valeur de lorsque ce qui est le cas de la paraboles ; on a, pour déterminer en l’équation du troisième degré

laquelle donne

et si l’on fait

on aura

et, de là,

Mais l’irrationnalité de l’expression de empêchera toujours que ces formules ne soient d’un grand usage dans le calcul analytique des orbites paraboliques ou presque paraboliques.

25. Il est bon de remarquer, relativement au mouvement parabolique, qu’on peut déterminer le temps employé à parcourir un arc quelconque de la parabole par une formule assez simple, qui ne renferme que la somme des rayons vecteurs qui répondentaux deux extrémités de l’arc et la corde qui sous-tend cet arc.

En faisant infini et les formules précédentes donnent

Marquons par un trait les mêmes quantités rapportées à un autre point de la parabole ; la différence ou le temps employé à parcourir un arc de parabole contenu entre deux points donnés, sera exprimé par la formule

Or, on a

et, si l’on nomme la corde qui joint les deux extrémités des rayons et on aura

Soit, pour abréger,

on a

équation d’où il s’agit de tirer la valeur de

Faisant disparaître les radicaux et ordonnant les termes par rapport à on a

d’où l’on tire

ou bien, en multipliant le haut et le bas par

Maintenant on a

donc


donc

substituant la valeur de cette quantité devient

Donc enfin, en remettant pour sa valeur, et faisant on aura

expression qui peut se mettre sous la forme suivante, plus simple,

comme on peut s’en assurer en prenant les carrés.

26. Cette formule élégante a été donnée d’abord par Euler, dans le septième Volume des Miscellanea Berolinensia. On pourraitla déduire du lemme X du troisième Livre des Principes mathématiques, en traduisant en analyse la construction par laquelle Newton détermine la vitesse qui ferait parcourir uniformément la corde d’un arc de parabole, dans le même temps que l’arc serait parcouru par une comète, et en observant que, dans la parabole, la demi-somme des rayons vecteurs qui aboutissent aux extrémités d’un arc quelconque est toujours égale au rayon vecteur qui aboutit au sommet du diamètre mené par le milieu de la corde parallèlement à l’axe, plus à la partie de ce diamètre interceptée entre l’arc et la corde ; d’où et du lemme IX on tire la valeur de e dernier rayon, exprimée par la corde et par la somme des rayons vecteurs qui répondent à ses deux extrémités.

On verra plus bas comment on peut étendre la même formule au mouvement elliptique ou hyperbolique[5].

27. Enfin l’équation entre et est toujours résoluble par approximation, lorsqu’on suppose le temps très petit ; on a alors, pour et, par conséquent, pour toutes les variables qui en dépendent, des séries ordonnées suivant les puissances de et qui seront d’autant plus convergentes que la valeur de sera plus petite. Mais, dans ce cas, il est plus simple d’en tirer la solution directement des équations différentielles en et de l’article 9, en y faisant

En regardant les variables comme des fonctions de et supposant qu’elles deviennent lorsque devient on a, en général, par le théorème connu,

et il ne s’agira que d’y substituer les valeurs des différentielles de déduites des trois équations

auxquelles on pourra joindre, pour simplifier le calcul, l’équation en de l’article 10

laquelle, étant différentiée et divisée par donne

d’où, en différentiant de nouveau et faisant, pour abréger,

on a celle-ci

laquelle est tout à fait semblable aux précédentes.

On aura ainsi, par des différentiations et des substitutions successives,

et ainsi de suite.

On aura de pareilles expressions pour les différentielles de et en changeant seulement en et

28. On fera donc ces substitutions, et, comme les quantités et leurs différentielles se rapportent, dans ces formules, au commencement du temps si l’on y change en et qu’on désigne par les valeurs de qui répondent à et qu’on suppose, pour abréger,

on aura ces expressions

À l’égard des constantes et que renferment ces expressions, il est bon de remarquer qu’elles se réduisent immédiatement aux constantes et d’où dépendent les éléments de l’orbite elliptique, comme nous l’avons vu dans l’article 8. Car, en rapportant au commencement du temps les deux équations en de l’article précédent, on a

savoir

et substituant pour et leurs valeurs et (art. 15), on aura

d’où l’on tire

On voit par là que les quantités et ne dépendent que de la figure de l’orbite, et nullement de la position de son plan.

29. Comme la quantité ou est déterminée par une équation différentielle semblable à celle qui détermine on aura aussi pour cette quantité une expression semblable, en changeant seulement et en et On aura ainsi

De là, en intégrant et ajoutant la constante

où les intégrales doivent être prises de manière qu’elles soient nulles

lorsque On aura ainsi, en substituant les valeurs de et et ordonnant les termes par rapport aux puissances de

Cette expression de doit être identique avec celle que donneraient les valeurs de car, puisque on aura aussi

Or

et

(art. 9) (art. 27)

de sorte qu’on aura

valeur qui coïncide avec la précédente.

§ II. — Détermination des éléments du mouvement elliptigue
ou parabolique.

30. Dans la théorie des planètes, on nomme éléments les six quantités constantes qui servent à déterminer la figure de l’orbite, sa position par rapport à un plan fixe, qu’on prend pour celui de l’écliptique, et l’époque ou le moment du passage par l’aphélie ou par le périhélie.

Soient, comme dans le paragraphe précédent, le demi grand axe ou la distance moyenne et le demi-paramètre ; ces deux éléments déterminent la figure de l’orbite ; et si l’on nomme l’excentricité, ou plutôt le rapport de la distance des deux foyers au grand axe, on a

et par conséquent,

Soit, de plus, le temps qui répond au passage de la planète par le périhélie ; cet élément, avec les deux précédents, servira à déterminer le mouvement elliptique, indépendamment de la position de l’orbite dans l’espace.

Pour déterminer cette position, soit la longitude du périhélie comptée depuis la ligne des nœuds, c’est-à-dire l’angle que la partie du grand axe qui répond au périhélie fait avec la ligne d’intersection du plan de l’orbite avec un plan fixe ; cet élément détermine la position de l’ellipse sur le plan de l’orbite.

Soit enfin l’inclinaison de ce plan sur le plan fixe auquel on le rapporte et qu’en Astronomie on prend ordinairement pour l’écliptique (nous le prenons dans nos formules pour celui des coordonnées ), et soit la longitude du nœud, c’est-à-dire l’angle que l’intersection des deux plans fait avec une ligne fixe, que les astronomes supposent dirigée vers le premier point d’Aries et que nous prenons pour l’axe des

Ces six quantités sont les éléments qu’il s’agit de déterminer, d’après quelques circonstances du mouvement elliptique donné.

31. Le cas le plus simple de ce problème est celui où l’on connaît la position du mobile, sa vinesse et sa direction dans un instant quelconque donné. Dans ce cas, les données sont les valeurs de pour un instant donné, valeurs que nous désignerons par les lettres romaines et il s’agira d’exprimer par ces six quantités les six éléments

L’article 9 donne d’abord, en mettant à la place de et changeant en

et les articles 11 et 15 donnent

On aura ainsi immédiatement, par ces formules, les valeurs du demi-axe du demi-paramètre d’où l’on tire l’excentricité et les angles et et il ne restera qu’à connaître les quantités et

32. Il est bon de remarquer que la valeur de et celle de peuvent se réduire à une forme plus simple. En effet, il est clair que est le carré de la vitesse initiale, laquelle étant nommée on aura

d’où l’on voit que le grand axe de la section conique et, par conséquent aussi, le temps périodique (art. 16) ne dépendent que de la distance primitive du corps au foyer attractif et de la vitesse de projection.

À l’égard du paramètre on a réduit, dans l’article 11, la quantité à la forme est l’angle décrit par le rayon dans l’instant de sorte que est le petit arc décrit par le même rayon ; par conséquent, est la vitesse perpendiculaire à ce rayon et que le corps a pour tourner autour du foyer.

Si l’on désignait cette vitesse de rotation par on aurait

et par conséquent

Ainsi le paramètre ne dépend que du rayon et de la partie de la vitesse par laquelle le corps tend à tourner autour du foyer vers lequel il est attiré.

33. Pour trouver la valeur de l’élément qui détermine le temps du passage par le périhélie, on remarquera que cette constante n’est entrée dans le calcul que par l’intégration qui a donné la valeur de en (art. 16).

Donc, si l’on dénote par la valeur de qui répond à on aura par les formules de l’article cité, en y faisant ce qui change en et en

Ainsi on aura par l’élimination de la valeur de en puisque et sont déjà connues.

Enfin, pour déterminer le dernier élément qui est aussi entré par l’intégration de l’équation entre et (art. 15), on remarquera d’abord que l’on a (art. 4), en changeant en et rapportant l’angle au commencement de

Ensuite l’article 7 donne

de sorte que, et étant déjà connus, on aura, par l’angle intermédiaire l’angle en et et de là on aura par l’équation de

l’article 15 rapportée à l’instant où

34. Si l’on connaissait deux lieux du mobile dans son orbite, avec le temps écoulé entre les instants où il a occupé ces deux lieux, on aurait aussi six données par les coordonnées qui répondent aux deux points donnés de l’orbite, et les six éléments seraient aussi déterminés par les valeurs de ces coordonnées ; mais l’expression transcendante du temps empêcherait de donner une solution générale et algébrique du problèmé. On pourra seulement le résoudre par approximation, si l’intervalle de temps entre les deux lieux est assez petit, en faisant usage des formules de l’article 28.

Soient les trois coordonnées du premier lieu dans l’orbite et celles du second lieu ; en prenant pour le temps écoulé entre les passages du mobile par ces deux lieux, on aura, en général (art. 28),

Supposons qu’on ne veuille porter la précision que jusqu’aux troisièmes puissances de on aura

Comme l’expression de renferme la constante

on commencera par la déterminer en ajoutant ensemble les trois équations précédentes, après avoir multiplié la première par la deuxième par et la troisième par on aura ainsi l’équation

d’où l’on tirera la valeur de qu’on substituera ensuite dans l’expression de

Les valeurs de et de étant connues, les mêmes équations donneront les valeurs des différentielles ainsi le problème sera réduit au cas précédent.

35. Enfin, si l’on ne connaissait que trois rayons vecteurs avec les temps et écoulés entre les passages par et et par et on pourrait encore déterminer l’orbite par les formules de l’article 29, en supposant les temps et assez petits.

Car, en faisant dans la valeur de et ne poussant les séries pour les valeurs de et que jusqu’aux et on aura

équations d’où l’on tirera les valeurs et Ces deux dernières donnent tout de suite, par les formules de l’article 19,

ensuite on aura l’angle compris entre le rayon et celui du périhélie, par la formule (art. 15)


étant égal

Si l’orbite était une parabole, on aurait et, par conséquent,

alors il suffirait de connaître deux distances et la première donnerait la valeur de et la seconde la valeur de par l’équation

36. Les éléments des planètes sont assez connus ; c’est par des observations de longitudes et de latitudes qu’on les a déterminés, et la petitesse de leurs excentricités et de leurs inclinaisons sur le plan de l’écliptique a beaucoup contribué à faciliter ces déterminations.

En prenant ce plan pour celui des les angles et (art. 4) représentent, l’un la longitude du corps, et l’autre sa latitude ; et nous avons donné dans les articles 22 et 23, pour les valeurs de et en des séries qui sont d’autant plus convergentes que l’excentricité et l’inclinaison sont plus petites. En prenant six observations, trois de longitude et trois de latitude correspondante ou, en général, de longitude ou de latitude, à des instants donnés, on aura six équations par lesquelles on pourra déterminer les six éléments, du moins pour le Soleil et la Lune, qui tournent immédiatement autour de la Terre.

Pour les autres planètes qui tournent autour du Soleil, le calcul est un peu plus compliqué, parce que l’observation ne donne immédiatement que les longitudes et latitudes vues de la Terre, qu’on nomme géocentriques ; mais, en supposant le mouvement du Soleil connu, on peut toujours déduire de chaque observation une équation ; de sorte que six observations suffiront, à la rigueur, pour la détermination des six éléments.

Ce problème est surtout important pour les comètes, dont les éléments, lorsqu’elles paraissent, sont tout à fait inconnus ; aussi, depuis Newton, qui a le premier tenté de le résoudre, il y a peu de géomètres et d’astronomes qui ne s’en soient occupés. Ne pouvant établir l’approximation sur la petitesse de l’excentricité et de l’inclinaison, comme pour les planètes, ils ont tous supposé que les intervalles de temps entre les observations sont très petits, et ils ont donné des méthodes plus ou moins approchées pour déduire les éléments des comètes de trois longitudes et d’autant de latitudes observées. Comme celle que j’ai proposée, dans les Mémoires de Berlin[6] pour 1783, me paraît offrir la solution la plus directe et la plus générale du problème des comètes, je crois pouvoir la donner ici, mais un peu simplifiée et accompagnée de remarques nouvelles ; elle fournira une application importante des principales formules que nous avons développées dans le paragraphe précédent.

§ III. — Sur la détermination des orbites des comètes.

37. Soient, dans un instant quelconque, la distance de la comète à la Terre et les cosinus des angles que la ligne ou le rayon visuel fait avec trois axes perpendiculaires entre eux et supposés fixes dans l’espace ; on aura pour les trois coordonnées rectangles de la comète, parallèles à ces axes et ayant leur origine dans le centre de la Terre. La quantité sera l’inconnue, mais les trois quantités seront connues par l’observation de la comète et devront être telles que l’on ait la condition

parce que, par l’hypothèse, on doit avoir

Soient de même les quantités correspondantes relativement au Soleil, en sorte que soient les coordonnées rectangles du lieu du Soleil par rapport à la Terre, et parallèles aux mêmes axes ; ces quantités doivent être censées connues par le calcul du lieu du Soleil dans le même instant de l’observation de la comète, et l’on aura aussi la condition

Enfin soient les coordonnées rectangles du lieu de la comète par rapport au Soleil, parallèles aux mêmes axes, et le rayon vecteur de son orbite autour du Soleil ; il est visible qu’on aura ces trois équations

et, comme on aura

Or on sait que l’expression est celle du cosinus de l’angle formé entre les deux rayons et partant du centre commun de la Terre et dirigés, l’un à la comète, l’autre au Soleil ; de sorte que, si l’on désigne cet angle par on aura

Si donc on a trois observations de la même comète, faites à des intervalles de temps connus, on aura trois systèmes pareils d’équations qui contiendront chacun une nouvelle inconnue et les propriétés de la parabole[7] donneront trois autres équations.

38. Ce qui se présente de plus simple pour cet objet est d’employer la formule donnée dans l’article 25, par laquelle on a le temps que la comète emploie à décrire un arc quelconque, exprimé par la corde de l’arc et par la somme des rayons vecteurs qui aboutissent à ses deux extrémités, et dégagée de tous les éléments de l’orbite ; car les trois intervalles de temps entre les trois observations, prises deux à deux, donneront les trois équations demandées.

Nous marquerons par un trait les lettres qui désignent les quantités analogues dans la seconde observation ; nous aurons ainsi

Pour la corde de l’arc parcouru par la comète, dans l’intervalle des deux observations, il est clair qu’on aura

En substituant pour et pour leurs valeurs, on aura

Or, par les théorèmes connus, l’expression doit repré-

senter le cosinus de l’angle compris entre les deux rayons et partant du centre de la Terre et dirigés aux deux lieux de la comète dans les deux observations ; de même, sera le cosinus de l’angle formé au centre de la Terre par les deux rayons dirigés aux deux lieux du Soleil, et ainsi des autres expressions semblables.

Donc si, pour plus de clarté, on imagine que les deux lieux apparents de la comète soient marqués sur la surface de la sphère par les lettres et de même les lieux apparents du Soleil par les lettres et qu’on joigne par des arcs de grands cercles les quatre points il est évident que les arcs représenteront les angles que nous avons dénotés par et que les arcs et représenteront les angles dont les cosinus sont

et qu’enfin les arcs et représenteront les angles dont les cosinus sont

Ainsi, en considérant le quadrilatère sphérique qui est censé donné par les deux observations de la comète et par les deux lieux calculés du Soleil, on aura

donc, comme la différence des temps et qui répondent aux deux observations, c’est-à-dire leur intervalle en temps, est censée donnée, on aura, par la formule de l’article cité, l’équation

dans laquelle il n’y aura d’inconnues que les deux distances et

Si l’on a une troisième observation, pour laquelle les quantités analogues soient désignées par les mêmes lettres marquées de deux traits, on aura, par la comparaison de la première observation avec celle-ci, une deuxième équation tout à fait semblable, dans laquelle les lettres marquées d’un trait dans l’équation précédente le seront par deux traits, et qui ne contiendra que les deux inconnues et

On aura de même une troisième équation semblable, par la comparaison de la deuxième observation avec la troisième, en ne faisant que marquer d’un trait, dans la première équation, toutes les lettres qui n’ont point de trait, et de deux traits toutes celles qui en ont un. Ainsi cette troisième équation ne contiendra que les mêmes inconnues et de sorte que les trois équations ne contiendront que les trois inconnues et suffiront pour les déterminer. Mais, quoique ces équations se présentent sous une forme assez simple, leur résolution offre des difficultés presque insurmontables, parce que les inconnues y sont mêlées entre elles et renfermées dans différents radicaux.

39. Au reste, si l’on pouvait parvenir d’une manière quelconque à trouver les valeurs des rayons on aurait tout de suite le demi-paramètre qui, dans la parabole, est égal au double de la distance périhélie, par la formule (art. 25)

et comme on a en général (art. 10) l’équation

dans laquelle

(art. 11), on aura ces deux-ci

d’où l’on tirera facilement les valeurs de et ce qui don-

nera la position du plan de la parabole relativement au plan qu’on aura choisi pour les axes des et

On peut même remarquer que, par le moyen de ces équations, qui dépendent de ce que l’orbite de la comète est supposée dans un plan passant par le Soleil, on peut d’abord réduire les trois inconnues à deux seulement.

En effet, si l’on fait, pour abréger,

on aura, en substituant les valeurs de l’équation

d’où l’on tire

et l’on aura de même les expressions de et de en marquant d’un trait et de deux traits les lettres, à l’exception de et qui sont les mêmes pour toutes les observations.

De cette manière, les trois inconnues seront réduites aux deux et de sorte qu’il ne faudra employer que deux équations pour leur détermination, ce qui simplifie un peu la solution du problème.

40. Pour la simplifier davantage, il ne paraît pas qu’il y ait d’autre moyen que de supposer les intervalles de temps entre les observations assez petits pour qu’on puisse négliger plusieurs termes comme insensibles, ce qui ne donnera d’abord qu’une solution approchée qu’on pourra rendre plus exacte ensuite par de nouvelles corrections. C’est aussi ce qu’on a fait jusqu’ici dans toutes les solutions qu’on a données de ce problème.

En appliquant cette hypothèse à la solution précédente, la corde deviendra très petite, et, en ne retenant que les deux premiers termes des radicaux qui entrent dans l’expression du temps écoulé entre les deux premières observations, on aura

ce qui donne

et il n’y aura plus d’autres radicaux dans cette équation que ceux qui entrent dans les expressions de et mais les équations entre les trois inconnues ou entre les deux seront encore trop compliquées pour qu’on puisse les employer avec succès.

On peut conclure de là que ces inconnues ne sont pas celles dont l’emploi est le plus avantageux dans la question présente ; et, lorsqu’on ne demande d’abord qu’une solution approchée, il est beaucoup plus simple de faire usage des formules que nous avons données dans l’article 28, pour le cas où l’on suppose le temps très petit.

41. Pour appliquer ces formules à la détermination de l’orbite des comètes, il n’y aura qu’à y substituer à la place de les expressions données dans l’article 37 ; on aura ainsi, en général,

où les quantités répondent au commencement du temps et sont regardées comme constantes, et où et sont des fonctions rationnelles de et des constantes

Comme le commencementdu temps est arbitraire, on peut le fixer au moment de la première observation ; or, en faisant on a

donc on aura pour la première observation ce premier système d’équations

Pour la deuxième observation, distante de la première du temps on aura, en marquant d’un trait les lettres ce second système d’équations

On aura des équations pareilles pour la troisième observation, distante de la première du temps en marquant de deux traits les lettres marquées d’un trait dans les dernières équations, et d’un trait les lettres et pour indiquer que le dont elles sont fonctions doit être changé en on aura ainsi ce troisième système d’équations

On peut éliminer des premières équations de chacun des trois systèmes les deux constantes et et faisant, pour abréger,

on aura

Éliminant de même les deux constantes des secondes équations

des mêmes systèmes, on aura

et l’élimination des constantes des dernières équations de ces systèmes donnera pareillement

De ces trois équations on tire

en supposant, pour abréger,

et en dénotant par ce que devient lorsqu’on y change respectivement en en et en en dénotant de même par et par ce que devient en faisant subir les mêmes changements aux quantités ainsi qu’aux quantités analogues

Maintenant les trois observations donnent aussi (art. 37, 38) les équations

Donc, si l’on y substitue les valeurs précédentes de on aura trois équations finales qui ne contiendront que des quantités connues, avec les quantités et qui sont données en fonctions du temps et des trois constantes d’où dépendent les éléments

de l’orbite (art. 28, 29) ; de sorte qu’on pourra déterminer ces trois constantes.

42. En ne poussant l’approximation que jusqu’aux quatrièmes puissances de on a

et de même,

et comme on aura

En faisant ces substitutions dans les valeurs de et supposant le coefficient étant donné par le rapport des deux intervalles entre les trois observations, il est clair que la quatrième dimension de disparaîtra par la division, et qu’ainsi il suffira d’avoir égard à la troisième dans les valeurs de et

Or on a, en général, aux près,

mais nous avons supposé que la première observation répond à et que les deux suivantes répondent aux temps et ainsi l’on aura

On fera donc ces substitutions dans les trois dernières équations de l’article précédent, et, rejetant les termes qui contiendraient des puissances de supérieures à la troisième, on aura trois équations entre les trois inconnues et dont les deux dernières n’y paraîtront que sous la forme linéaire, de sorte qu’il sera très facile de les éliminer et de réduire le problème à une seule équation en C’est en quoi consiste le principal avantage de la méthode que nous proposons.

Si l’on voulait pousser l’approximation plus loin et avoir égard à un plus grand nombre de termes dans les valeurs de on aurait des équations où les inconnues et ne seraient plus linéaires, mais monteraient successivement à des dimensions plus hautes, ce qui rendrait leur élimination plus difficile et l’équation finale encore plus compliquée.

43. Pour donner là-dessus un essai de calcul, nous nous contenterons d’avoir égard, dans les valeurs de aux troisièmes dimensions de et de ce qui fera disparaître les termes affectés de l’inconnue nous ferons, pour plus de simplicité, en prenant la distance moyenne de la Terre au Soleil pour l’unité des distances, et représentant les temps par les mouvements moyens du Soleil (art. 23) ; et, supposant nous aurons

Les valeurs de deviendront ainsi de la forme

en supposant, pour abréger,

et dénotant par et ce que et deviennent en y changeant en et en

Ces valeurs de distances de la comète à la Terre dans les trois observations, ne contiennent, comme l’on voit, que la seule inconnue rayon vecteur de la comète dans la première observation. Si donc on substitue la valeur de dans l’équation (art. 25)

on aura une équation finale en laquelle montera au huitième degré, et le problème sera réduit à la résolution de cette équation.

Ayant trouvé la valeur de on aura par les formules précédentes celles de et de là on aura, par les formules de l’article 42, les valeurs des trois rayons vecteurs ainsi que celles des coordonnées et de leurs différentielles et l’on pourra déterminer l’orbite par les formules du § II, ou, si l’on aime mieux, par les formules trigonométriques connues, d’après les trois distances de la comète à la Terre.

44. Les expressions des distances peuvent être simplifiées par la considération suivante comme la Terre et la comète se meuvent autour du Soleil par la même force attractive de cet astre, si l’on nomme les coordonnées rectangles de la Terre autour du Soleil lorsque et qu’on désigne par ce que deviennent les fonctions et lorsqu’on y change les éléments de l’orbite de la comète en ceux de la Terre, on aura, comme dans l’article 28, les trois équations

parce que, ayant dénoté (art. 24) par les coordonnées rec-

tangles du lieu du Soleil par rapport à la Terre, on aura pour celles de la Terre par rapport au Soleil.

Comme ces équations ne diffèrent de celles de l’article 41 que parce que sont changées en et que y est nul, il est clair qu’on aura des résultats analogues en faisant ces mêmes changements dans ceux que nous venons de trouver dans l’article précédent. Ainsi, puisque les expressions de données à la fin de cet article ne contiennent d’autre quantité dépendante des éléments de l’orbite que le rayon vecteur si l’on change en rayon vecteur de l’orbite de la Terre, on aura

d’où l’on tire

Ces valeurs étant maintenant substituées dans les mêmes expressions de et négligeant, dans le dénominateur, le terme très petit du second ordre en vis-à-vis du terme fini en on aura ces expressions plus simples

Si donc on substitue la valeur de dans l’équation

et qu’on fasse, pour abréger,

quantité toute-connue par les observations, en multipliant par on

aura l’équation

où l’inconnue montera au huitième degré, mais qui, étant divisible par ne sera, après la division, que du septième degré.

Cet abaissement de l’équation en est dû à ce que nous avons représenté le mouvement de la Terre, comme celui de la comète, par des formules approchées où l’on a négligé les puissances de supérieures à la troisième ; il n’aurait pas lieu en employant la valeur de de l’\piticle précédent, dans laquelle les lieux du Soleil sont supposés exacts, étant déterminés d’après les Tables.

45. On peut ramener l’équation précédente à une construction assez simple. Ayant mené d’un point donné deux droites qui fassent entre elles un angle égal à l’arc distance apparente de la comète au Soleil dans la première observation, et dont la première soit égale à et la seconde égale à il s’agira de trouver dans la première un point tel que la partie comprise entre ce point et l’extrémité de la même droite soit à la droite entière comme le cube de la seconde droite est au cube de la droite qui joindra l’extrémité de celle-ci et le point cherché alors cette dernière droite sera égale à et la partie de la première interceptée entre le point donné et le point cherché sera égale à Car, par cette construction, on aura la proportion

laquelle donne

et ensuite

d’où résulte l’équation ci-dessus en

Lambert, est, je crois, le premier qui ait réduit le problème des comètes, envisagé d’une manière approchée, mais exacte[8], à une équation unique à une seule inconnue. Il y est parvenu par une considération ingénieuse, fondée sur ce que le lieu apparent de la comète, dans la deuxième observation, s’écarte du grand Cercle mené par les lieux apparents dans la première et dans la troisième observation ; et la détermination de cet écartement l’a conduit directement à une construction analogue à celle que nous venons de donner, et qui se réduit à une équation en du septième degré. (Voir les Mémoires de l’Académie de Berlin pour l’année 1771.)

Connaissant ainsi les valeurs de et on aura

et les deux équations (art. 40 et 41)

donneront les valeurs des constantes et et de là celles des éléments et de l’orbite, par les formules de l’article 28, étant le grand axe et le paramètre.

46. Si l’on suppose l’orbite parabolique, on aura infini, ce qui donne Dans ce cas, les deux dernières équations ne contiendront plus que l’inconnue laquelle étant éliminée, on aura une nouvelle équation en qui devra avoir une racine commune avec celle qu’on a déjà trouvée, ce qui servira à faciliter la recherche de cette racine.

En adoptant d’abord pour les comètes l’hypothèse de la parabole, il sera préférable de faire dépendre la solution uniquement de cette dernière équation, parce qu’elle a l’avantage d’être exempte de la quantité qui est très petite du troisième ordre lorsque les intervalles de temps et ou sont très petits du premier, comme on le verra plus bas, de sorte que les erreurs des observations d’où cette quantité dépend peuvent y avoir une influence très grande.

En faisant, pour abréger,

et négligeant les termes affectés de dans les coefficients de l’élimination de cette quantité donnera l’équation en

qui, étant combinée avec l’équation

donnera, par l’élimination de une équation en du sixième degré ; et si, dans la combinaison des deux équations, on néglige le carré du terme qui serait du quatrième ordre, l’équation finale ne montera plus qu’au cinquième degré. On pourrait même, dans la première approximation, négliger ce terme, qui n’est que du second ordre ; alors l’équation finale ne serait plus que du quatrième degré et pourrait se résoudre directement par les méthodes connues.

La valeur de donnera celles de et de là celle de par les formules de l’article précédent ; et, comme on suppose infini, on aura

où le demi-paramètre devient double de la distance périhélie de la comète.

47. Après avoir réduit le problème des comètes à des équations finales à une seule inconnue, il reste à examiner les quantités qui doivent être supposées connues ; ces quantités sont :

1o Les trois rayons qui représentent les distances du Soleil à la Terre dans les trois observations, et qui doivent être calculés par les Tables du Soleil ;

2o Les quantités d’où dépendent les valeurs de (art. 41 et 43). Celles-ci doivent être déterminées par les trois observations de la comète et par le calcul des lieux du Soleil ; mais on peut les ramener à des expressions plus simples, qui en rendront la détermination beaucoup plus facile.

Commençons par la quantité dont les autres ne sont que des dérivées on a (art. 41)

le carré de cette expression peut se mettre sous la forme

comme on peut s’en convaincre par le développement. Or, par la nature des quantités (art. 37),

Donc, faisant, pour abréger,

on aura

Or nous avons déjà remarqué (art. 38) que la quantité est égale au cosinus de l’angle compris entre les deux rayons et dirigés vers la comète dans les deux premières observations, angle que nous avons désigné par le côté du triangle sphérique supposé tracé sur la sphère en joignant par des arcs de grands cercles les trois lieux apparents de la comète dans les trois observations. Ce triangle est entièrement donné par les observations de la comète, de quelque manière qu’elles aient été faites ; et nous pouvons regarder comme connus ses trois côtés ainsi que les angles qui sont respectivement opposés aux côtés et

On aura donc

et de même

et l’expression de la quantité deviendra

Cette expression de peut encore se réduire à une forme plus simple ; car il est facile de se convaincre, par le développement des termes, qu’elle est la même chose que celle-ci

laquelle, par les transformations connues, devient celle-ci

formule très commode pour le calcul logarithmique.

Si l’on veut employer les angles du même triangle, on peut avoir encore une expression plus simple de la quantité car on a, par les formules connues,

si l’on fait cette substitution dans la première expression de on aura, après les réductions,

et, par conséquent, en tirant la racine carrée,

On peut trouver de la même manière

Il est facile de prouver que la quantité n’est autre chose que la solidité, prise six fois, de la pyramide triangulaire qui a le sommet au centre de la sphère dont le rayon est supposé égal à l’unité, et qui s’appuie sur le triangle sphérique c’est-à-dire qui a pour base le triangle rectiligne formé par les cordes des trois arcs car, si l’on considère une des faces triangulaires de cette pyramide, celle, par exemple, qui a pour base la corde de l’arc on aura pour l’aire de ce triangle isoscèle. Ensuite, si l’on considère la-face ad\sqrt{a}cette qui a pour base la corde de l’arc il est visible que l’inclinaison mutuelle de ces deux faces sera égale à l’angle du triangle sphérique ; par conséquent, la perpendiculaire menée de l’angle sur la première face sera égale à Cette perpendiculaire devient la hauteur de la pyramide, en la supposant couchée sur la première face égale à donc la solidité de la pyramide sera égale à

et, par conséquent, égale à

48. Nous dénoterons, en général, par le symbole la fonction des côtés et des angles de tout triangle sphérique par laquelle nous avons exprimé la quantité

Ainsi, ayant marqué sur un globe les trois lieux apparents de la comète donnés par les trois observations, et formé le triangle sphérique on aura tout de suite

Si ensuite on place sur le même globe les trois lieux du Soleil dans les trois observations, et qu’en joignant ces lieux et ceux de la comète par des arcs de grands cercles on forme différents triangles sphériques il est facile de voir, par ce que nous avons dit dans l’article 40, relativement aux quantités que les trois premières seront données par des fonctions semblables des triangles que les trois autres seront données par de pareilles fonctions des triangles et que les trois dernières le seront par de semblables fonctions des triangles On aura donc ainsi, d’après la même notation,

Ces quantités ne dépendent, comme l’on voit, que de la position mutuelle des lieux apparents de la comète et du Soleil, et, comme elles sont les seules qui entrent dans les équations qui déterminent les étéments absolus de l’orbite, notre analyse a l’avantage de séparer la détermination de ces éléments de celle des autres éléments qu’on peut appeler relatifs, parce qu’ils se rapportent à la position de l’orbite dans l’espace.

49. On peut remarquer encore que les expressions que nous venons de donner ont lieu quelle que soit la position des lieux apparents de la comète et du Soleil ; mais lorsque, comme nous l’avons supposé, les lieux de la comète sont peu distants entre eux, les arcs seront très petits, et l’angle compris entre ces arcs sera peu différent de deux droits ; il serait égal à deux droits si la Terre et la comète décrivaient, dans l’intervalle de la première à la troisième observation, des lignes droites, parce qu’alors les trois lieux apparents de la comète seraient dans un même grand cercle. Les sinus de et de seront donc très petits, et la quantité

sera très petite du troisième ordre ; mais les quantités

ne seront que du premier ; et, comme d’ailleurs il n’entre dans la valeur de que des quantités dépendantes des lieux apparents de la comète, au lieu que les quantités dépendent en partie des lieux du Soleil, qui, étant donnés par les Tables, peuvent être regardés comme exacts, il s’ensuit que la valeur de la quantité sera toujours beaucoup plus sujette à erreur que celles des quantités et qu’ainsi il conviendra, autant qu’il est possible, de l’éviter, comme nous l’avons montré dans l’article 46.

50. Nous remarquerons enfin que, comme l’observation d’une \inftymète donne ordinairement son ascension droite et sa déclinaison, si l’on veut employer immédiatement ces données dans nos formules, il n’y aura qu’à supposer que les trois axes auxquels nous avons rapporté les rayons dirigés vers la comète, et les rayons dirigés vers le Soleil, soient dirigés, le premier vers l’équinoxe du printemps, le deuxième à angle droit sur le plan de l’équateur et suivant l’ordre des sjgnes, et le troisième vers le pôle boréal de l’équateur alors, nommant l’ascension droite de la comète, sa déclinaison dans la première observation, et de même l’ascension droite du Soleil, sa déclinaison au même instant, il est facile de voir qu’on aura

De là on aura (article cité)

et, pareillement,

en marquant, comme nous l’avons fait, par un trait et par deux traits les quantités analogues qui se rapportent à la deuxième et à la troisième observation.

On aura de la même manière

et ainsi des autres cosinus.

Si ensuite on substitue ces mêmes valeurs de dans l’expression de on aura

et l’on en déduira les valeurs de en changeant et en et en et en et et celles de en faisant les mêmes changements sur et et celles de en faisant ces mêmes changements sur et On aura ainsi

et pour avoir les valeurs de et de il n’y aura qu’à marquer, dans les expressions de les lettres et d’un trait et de deux traits.

Il est inutile d’observer que si, au lieu des ascensions droites et des déclinaisons, on avait pour données les longitudes et les latitudes, il n’y aurait qu’à substituer ces données à la place de celles-là dans les mêmes formules ; l’orbite se trouverait alors rapportée à l’écliptique, au lieu de l’être à l’équateur.

51. Après avoir calculé ces valeurs, on calculera celles des quantités par la formule de l’article 42 ; et si l’on veut employer la méthode de l’article 44, comme la plus courte, on aura tout de suite l’équation finale en dont la résolution ne sera pas difficile, en la réduisant pour la première approximation au quatrième degré.

Si les intervalles entre les observations étaient égaux, on aurait et, par conséquent, ce qui donnerait

en désignant par la caractéristique la différence seconde des quantités dans lesquelles il n’y a que les quantités relatives au Soleil qui varient. Or, comme on suppose les observations peu distantes entre elles, les différences de ces quantités seront très petites, par conséquent la différence seconde sera très petite du second ordre et pourra être négligée vis-à-vis de la quantité finie ce qui réduira la valeur de à cette seule quantité, et l’on pourra faire les mêmes réductions sur les quantités analogues de sorte qu’on aura simplement

ce qui abrégera encore le calcul de la première approximation.

À l’égard de la mesure du temps, comme ce temps doit être représenté par le mouvement moyen du Soleil, si on veut l’exprimer en jours moyens, il suffira de multiplier le nombre des jours et des décimales de jour par l’angle du mouvement moyen du Soleil dans un jour, réduit en parties du rayon. Cet angle est de et donne en parties du rayon le nombre par lequel il faudra donc multiplier les intervalles de temps réduits en jours moyens.

CHAPITRE DEUXIÈME.

Sur la variation des éléments des orbites elliptiques produite par une force d’impulsion ou par des forces accélératrices.

52. Un des premiers et des plus beaux résultats de la théorie de Newton sur le système du monde consiste en ce que toutes les orbites des corps célestes sont de même nature et ne diffèrent entre elles qu’à raison de la force de projection que ces corps peuvent être supposés avoir reçue dans l’origine des choses. Il suit de là que, si une planète ou une comète venait à recevoir une impulsion étrangère quelconque, son orbite en serait dérangée ; mais il n’y aurait que les éléments, qui sont les constantes arbitraires de l’équation, qui pourraient changer c’est ainsi que l’orbite circulaire ou elliptique d’une planète pourrait devenir parabolique ou même hyperbolique, ce qui transformerait la planète en comète.

Il en est de même de tous les problèmes de Mécanique. Comme les constantes arbitraires introduites par les intégrations dépendent de l’état initial du système, qui peut être placé dans un instant quelconque, si l’on suppose que les corps viennent à recevoir pendant leur mouvement des impulsions quelconques, les vitesses produites par ces impulsions, étant composées avec les vitesses acquises par les corps, pourront être regardées comme des vitesses initiales et ne feront que changer les valeurs des constantes.

Et si, au lieu d’impulsions finies, qui n’agissent que dans un instant, on suppose des impulsions infiniment petites, mais dont l’action soit continuelle, les mêmes constantes deviendront tout à fait variables et serviront à déterminer l’effet de ces sortes de forces, qu’il faudra regarder comme des forces perturbatrices. On aura alors le problème dont nous avons donné une solution générale dans la Section et que nous appliquerons ici aux orbites des planètes.

§ I. — Du changement produit dans les éléments de l’orbite d’une planète,
lorsqu’elle est supposée recevoir une impulsion quelconque
.

53. Nous avons vu, dans le § II du Chapitre précédent, comment on peut exprimer tous les éléments du mouvement elliptique d’une planète par des fonctions des coordonnées et de leurs différentielles qui expriment les vitesses suivant les directions de ces coordonnées. Si donc on suppose qu’une planète, pendant qu’elle se meut, reçoive dans un lieu quelconque de son orbite une impulsion qui lui communique les vitesses suivant les mêmes coordonnées et tendantes à les augmenter, il n’y aura qu’à mettre dans les mêmes fonctions

à la place de et l’on aura les éléments de la nouvelle orbite que la planète décrira après l’impulsion.

Si, à la place des coordonnées rectangles on prend, comme dans l’article 5, le rayon vecteur avec les angles et , dont le premier soit l’inclinaison de sur le plan fixe des et dont l’autre soit l’angle de la projection de sur ce plan avec l’axe fixe des les expressions de l’orbite deviennent plus simples.

En effet, en substituant et à la place de on trouve, pour les éléments

Dans ces formules, les expressions différentielles et représentent les vitesses dans la direction du rayon dans une direction perpendiculaire à ce rayon et parallèle au plan de projection, et dans une direction normale au plan des deux autres composantes.

54. Prenons, pour plus de simplicité, le plan de projection dans le plan même de l’orbite, et supposions que la vitesse reçue par l’impulsion soit décomposée en trois, l’une suivant le rayon l’autre perpendiculaire à ce rayon dans le plan de l’orbite, et la troisième perpendiculaire à ce plan. Si l’on désigne la première par la deuxième par et la troisième par on aura les éléments de la nouvelle orbite après l’impulsion, en mettant dans les expressions précédentes à la place de et faisant alors la position de la nouvelle orbite se trouvera rapportée au plan de l’orbite primitive.

Soient ce que les éléments deviennent pour la nouvelle orbite ; on aura

donc

en effet, il est clair que le nœud de la nouvelle orbite avec l’orbite primitive doit être dans le lieu où se fait l’impulsion.

Si l’on fait aussi et dans les expressions des éléments primitifs et on a

et de là on tire

En substituant ces valeurs, on aura les éléments de la nouvelle orbite exprimés par ceux de l’orbite primitive et par les vitesses produites par l’impulsion.

55. Supposons maintenant qu’on demande l’impulsion nécessaire pour changer les éléments primitifs en et pour rendre la nouvelle orbite inclinée à la première avec l’angle il ne s’agira que d’avoir les expressions en et Les formules que nous venons de trouver donnent

Soit la vitesse imprimée par l’impulsion, et soient les angles que la direction de l’impulsion fait avec trois axes dont l’un soit le rayon prolongé, l’autre perpendiculaire à ce rayon dans le plan de l’orbite primitive et dans le sens du mouvement de la planète, et le troisième perpendiculaire au même plan ; on aura, par le principe de la décomposition, pour les trois vitesses suivant ces axes, lesquelles sont aussi celles que nous avons désignées par On aura donc

d’où l’on tire, à cause de

Donc, si l’on fait, pour abréger,

on aura

Mais, si l’on voulait rapporter la direction de l’impulsion à deux autres axes placés dans le plan de l’orbite primitive, dont l’un serait perpendiculaire et l’autre tangent à cette orbite, alors, en nommant l’angle que la perpendiculaire à l’orbite fait avec le rayon vectcur et dont la tangente est exprimée par les vitesses imprimées suivant ces deux axes seront

la vitesse suivant le troisième axe perpendiculaire au plan de l’orbite demeurant la même. Si donc on désigne par et par les angles que la direction de l’impulsion fait avec ces nouveaux axes, on aura

Or on a

d’où l’on tire, en substituant la valeur de

et de là on aura

où l’on remarquera que est la vitesse dans l’orbite primitive.

À l’égard des signes ambigus des radicaux qui entrent dans ces formules, on remarquera :

1o Que étant la valeur de exprime la vitesse suivant le rayon dans l’orbite primitive, et que exprimera la vitesse suivant ce rayon dans la nouvelle orbite ; ainsi il faudra prendre ces quantités positivement ou négativement, suivant que les vitesses qu’elles représentent tendront à augmenter ou à diminuer le rayon c’est-à-dire à éloigner ou à rapprocher le corps du foyer ;

2o Que étant égal à représente la vitesse circulatoire autour du foyer dans l’orbite primitive, et que, de même, représentera la vitesse circulatoire dans la nouvelle orbite, et sera cette vitesse circulatoire rapportée au plan de l’orbite primitive. Ainsi, en prenant positivement, il faudra prendre l’autre radical positivement ou négativement, suivant que la nouvelle orbite sera, par rapport au plan de l’orbite primitive, dans le même sens que dans cette orbite ou en sens contraire, c’est-à-dire suivant que le mouvement dans la nouvelle orbite sera direct ou rétrograde, relativement au mouvement dans l’orbite primitive.

56. Lorsqu’on voudra appliquer ces formules aux planètes et aux comètes, on fera en prenant la distance moyenne de la Terre au Soleil pour l’unité des distances et la vitesse moyenne de la Terre dans son orbite pour l’unité des vitesses. Cette vitesse est à peu près de lieues, de au degré, par seconde. La vitesse d’un boulet de au sortir du canon, est d’environ pieds, ou toises par seconde, laquelle est aussi à peu près celle d’un point de l’équateur dans le mouvement diurne de la Terre, celle-ci étant de toises par seconde. Donc si, pour rendre nos estimations plus sensibles, nous prenons pour unité cette vitesse d’un boulet de laquelle est à peu près d’un dixième de lieue, la vitesse de la Terre dans son orbite sera exprimée par le nombre par conséquent, il faudra multiplier par la valeur de la vitesse d’impulsion.

Voyons quelle peut être la plus grande valeur de

En nommant l’excentricité de l’orbite primitive et l’anomalie vraie qui répond au rayon on a (art. 15)

donc

Ainsi la plus petite valeur de sera et de même la plus petite valeur de sera en nommant l’excentricité de la nouvelle orbite. Donc la plus grande valeur des termes sera et cette expression aura lieu aussi pour les orbites hyperboliques où et surpasseraient l’unité.

Par les mêmes formules, on a dont la plus grande valeur est la plus grande valeur de sera de même donc la plus grande valeur de sera mais il est facile de prouver que

car la différence de leurs carrés est donc on aura toujours

Il faut encore chercher les plus grandes valeurs de et Or, les plus petites valeurs de et de étant la plus grande valeur de sera et de même la plus grande valeur de sera

Donc, puisque dans les expressions de et peuvent avoir les signes ou en prenant positivement les termes qui contiennent ces radicaux, et donnant aussi à sa plus grande valeur on aura

Cette limite se réduira à lorsque l’orbite primitive sera circulaire ou presque circulaire comme celle des planètes, et que la nouvelle sera parabolique comme celles des comètes.

57. Les principales circonstances du mouvement des planètes autour du Soleil nous portent à croire qu’elles ont eu une origine commune c’est le contraire pour les comètes ; elles n’ont de commun entre elles que le mouvement dans une parabole, ou, en général, dans une section conique, et elles paraissent avoir été jetées au hasard dans l’espace.

Ne peut-on pas supposer que la cause qui a produit nos planètes en a produit en même temps un plus grand nombre d’autres placées au delà de Saturne, et décrivant des orbites semblables, comme Uranus, mais dont plusieurs seront devenues ensuite comètes, en éclatant par une explosion interne ? Car, une planète étant brisée en deux ou plusieurs morceaux, par la force de l’explosion, chacun de ces morceaux recevra une impulsion qui lui fera décrire une orbite différente de celle de la planète ; et, pour que cette orbite soit parabolique, il suffira que la vitesse imprimée par l’explosion n’excède pas fois celle d’un boulet de canon[9]. Pour Saturne on a et pour Uranus en supposant il suffira d’une vitesse moindre que fois celle d’un boulet, qui n’est produite que par une poignée de poudre.

L’hypothèse d’une planète brisée par une explosion interne a déjà été proposée par M. Olbers, pour expliquer la presque égalité des éléments des quatre nouvelles planètes ; et ce qui pourrait la confirmer, ce sont les variations de lumière qu’on observe dans ces planètes, et qui, en indiquant un mouvement de rotation, indiquent en même temps que leur figure n’est pas de révolution comme celles des autres planètes ; que, par conséquent, elles ne pouvaient pas être fluides, mais qu’elles devaient être déjà durcies lorsqu’elles sont devenues planètes comme elles le sont dans l’état actuel.

Si l’on suppose l’orbite primitive circulaire, et l’orbite changée par l’explosion elliptique, mais peu différente d’un cercle et peu inclinée au plan de l’orbite primitive, et qu’on n’ait égard qu’aux premières dimensions de l’excentricité et du sinus de l’inclinaison on a

l’angle étant celui que le rayon fait avec le rayon du périhélie,

Ainsi, puisque les excentricités et les inclinaisons des planètes ne gardent entre elles aucune loi et n’ont de commun que leur petitesse, on pourrait supposer que les orbites des planètes ont été circuiaires dans leur formation, et qu’elles sont devenues ensuite elliptiques et inclinées par l’effet de petites explosions internes. En effet, si un petit morceau de la masse d’une planète en avait été détaché et lancé avec une vitesse capable d’en faire une comète, la planète n’aurait reçu en sens contraire qu’une petite vitesse qui aurait pu changer son orbite circulaire en elliptique et inclinée, comme celles de nos planètes, et la même impulsion aurait pu produire aussi quelque changement sur sa rotation, comme nous le verrons plus bas.

§ II. — Variations des éléments des planètes produites
par des forces perturbatrices
.

58. Supposons maintenant que les impulsions qui changent les constantes arbitraires soient infiniment petites et continuelles ; ces constantes deviennent variables, et l’on pourra, de cette manière, réduire l’effet des forces perturbatrices des planètes aux variations des éléments de leurs orbites.

Soient les forces perturbatrices décomposées suivant les directions des coordonnées rectangles et tendantes à augmenter ces coordonnées ; ces forces engendreront pendant l’instant les petites vitesses qu’il faudra ajouter aux vitesses dans l’expression de chacun des éléments comme dans l’article 52. Mais, comme ces vitesses additionnelles sont ici infiniment petites, elles ne produiront dans les éléments que des variations infiniment petites, qu’on pourra déterminer par le Calcul différentiel.

Faisons, pour abréger,

chacun des éléments sera exprimé par une fonction donnée de Soit un quelconque de ces éléments ; on aura sa variation en augmentant des quantités infiniment petites on aura ainsi

et l’on aura de pareilles équations pour les autres éléments de l’orbite

Pour faire usage de ces équations, il faudra substituer à la place des variables leurs valeurs en et en données par les formules trouvées dans le premier Chapitre ; on aura ainsi autant d’équations du premier ordre entre le temps et les éléments devenus variables qu’il y a de ces éléments, et il ne s’agira plus que de les intégrer.

Si l’on voulait introduire directement les forces perturbatrices dans les équations de l’orbite primitive (art. 4), il n’y aurait qu’à ajouter respectivement les quantités aux termes de ces équations. Ainsi l’on peut regarder les équations précédentes entre les nouvelles variables comme des transformées des équations en mais ces transformations seraient peu utiles pour la solution générale du problème. Leur grande utilité est lorsque la solution rigoureuse est impossible, et que les forces perturbatrices sont très petites ; elles fournissent alors un moyen d’approximation que nous avons exposé d’une manière générale dans la Section V.

59. Cette approximation, fondée sur la variation des éléments, est surtout applicable aux orbites elliptiques des planètes, en tant qu’elles sont dérangées par l’action des autres planètes, et les géomètres l’ont souvent employée dans la théorie des planètes et des comètes ; on peut dire que ce sont les observations elles-mêmes qui l’ont fait connaître avant qu’on y eût été conduit par le calcul ; elle a l’avantage de conserver la forme elliptique des orbites, et même de supposer l’ellipse invariable pendant un temps infiniment petit, de manière que non seulement le lieu de la planète, mais aussi sa vitesse et sa direction[10] ne soient point affectés de la variation instantanée des éléments.

En effet, en regardant les coordonnées comme des fonctions du temps et des éléments devenus variables, on a, par la différentiation,

et il est facile de prouver que la partie qui contient les variations devient nulle par la substitution de la valeur de donnée ci-

dessus, et des valeurs semblables de car, en faisant ces substitutions dans les termes et ordonnant par rapport aux quantités on aura

Mais, en regardant d’abord comme fonctions de et ensuite comme fonctions de on a

Substituant dans l’expression de ces valeurs de on doit avoir des équations identiques ; par conséquent, il faudra que les termes affectés de dans les expressions de deviennent nuls ; ce qui donnera, par rapport à les équations identiques

On aura donc simplement

et l’on trouvera de la même manière

comme si les constantes ne variaient point.

60. Lorsque les forces perturbatrices viennent des attractions d’autres corps fixes ou mobiles, et que ces attractions sont proportionnelles à des fonctions des distances, alors, si l’on désigne, comme dans l’article 8 de la Section V, par la somme des intégrales de chaque force multipliée par l’élément de sa distance au centre d’attraction, et qu’on regarde la quantité comme fonction de les forces sont de la forme

je donne ici le signe à la quantité parce que j’ai supposé que les forces tendent à augmenter les distances au lieu que, dans les fonctions les forces perturbatrices, dirigées sur des centres, sont supposées tendre à diminuer les distances des corps à ces centres.

Dans ce cas, qui est celui de la nature, les variations des éléments peuvent s’exprimer d’une manière plus simple, en employant, au lieu des différences partielles de relatives à ses différences partielles relatives à après la substitution des valeurs de en et c’est cette considération qui a fait naître la nouvelle théorie de la variation des constantes arbitraires.

Si l’on regarde comme fonctions de on aura

et, ces valeurs étant substituées dans l’expression de de l’article 58, à la place de elle deviendra

On peut faire disparaître de cette expression les termes multipliés par par la considération que, ne contenantpointles variables on a

Donc, si l’on soustrait de la valeur de la quantité


Cette expression de est, en apparence, plus compliquée que la formule primitive d’où nous sommes partis ; mais elle a, d’un autre côté, le grand avantage que les coefficients des différences partielles deviennent indépendants du temps après la substitution des valeurs de en et données par le mouvement elliptique de la planète, comme on peut s’en assurer par la différentiation, en faisant varier le temps dans les coefficients.

61. En effet, puisque est censé fonction de et que varient aussi avec de manière que

et

par les équations différentielles du problème (art. 4), il s’ensuit qu’on aura, en différentiant par rapport à

Mais, étant une des constantes arbitraires introduites par l’intégration des mêmes équations, sa différentielle relative à doit devenir identiquement nulle par les mêmes valeurs de on aura donc

équation identique qui subsistera, par conséquent, en faisant varier séparément

Faisons varier on aura donc aussi

Donc la valeur de la différentielle de se réduira à

On trouve de la même manière

On aura ensuite

Mais, en faisant varier dans l’équation identique et observant que la fonction est supposée ne pas contenir les variables on a

donc on aura simplement

et l’on trouvera de la même manière

On aura des expressions semblables pour les différentielles en changeant seulement en et ainsi pour les autres quantités semblables.

Si maintenant on différentie le coefficient de dans l’expression de de l’article 60, et qu’on y substitue les valeurs qu’on vient de trouver pour les différentielles de on verra d’abord que les termes qui contiendront les différentielles de se détruiront mutuellement, et que tes termes qui contiendront les différences de étant ordonnés par rapport aux différences partielles de se détruiront aussi mutuellement dans chacun des coefficients de ces différences partielles.

D’où l’on peut conclure que le coefficient de dans l’expression de sera constant à l’égard du temps et ne pourra être qu’une fonction de après la substitution des valeurs de en et de sorte que la variable disparaîtra d’elle-même, et qu’il suffira d’y substituer les valeurs de qui répondent ou à ou à une valeur quelconque de

On prouvera de la même manière que disparaîtra des autres coefficients des différences partielles de dans la même expression de . Ainsi la variation de sera représentée par une formule qui ne contiendra que les différences partielles de par rapport à multipliées chacune par une fonction de sans Et la même chose aura lieu à l’égard des variations des autres constantes arbitraires

Ce résultat important, que nous venons de trouver a posteriori, n’est qu’un cas particulier de la théorie générale de la variation des constantes arbitraires, que nous avons exposée dans le § II de la Section V ; et nous aurions pu le déduire immédiatement de cette théorie ; mais nous avons cru qu’il n’était pas inutile de montrer comment on peut y arriver en partant des formules qui donnent directement les variations des éléments dues aux forces perturbatrices, et surtout comment ces variations acquièrent une forme simple et élégante par la réduction des forces perturbatrices aux différences partielles d’une même fonction, relatives à ces mêmes éléments regardés comme variables.

62. Nous avons supposé, dans l’article 60, que les forces pouvaient s’exprimer par les différences partielles d’une même fonction relatives à Cette hypothèse simplifie le calcul, mais n’est pas absolument nécessaire pour son exactitude, puisque les équations différentielles sont toujours indépendantes de la nature des forces accélératrices du mobile ; il s’agit seulement de savoir ce qu’on doit substituer à la place des différences partielles de relatives aux constantes arbitraires Or ces constantes n’entrent dans la fonction que parce qu’elles entrent dans les expressions des variables dont est supposé fonction ; ainsi l’on aura

et, remettants à la place de on aura

quelles que soient les valeurs de Il en sera de même à l’égard de en changeant en

En général, si l’on dénote par la caractéristique les variations de relatives aux constantes arbitraires on aura

et si l’on suppose que les forces perturbatrices soient tendantes à des centres dont les distances respectives soient ce qui donne

on aura aussi, relativement aux constantes arbitraires,

Je donne ici à le signe parce que les forces sont supposées tendre à diminuer les distances au lieu que les

forces sont supposées tendre à augmenter les lignes comme nous l’avons déjà observé dans l’article [11].

63. Pour appliquer les formules générales de l’article 18 de la Section citée aux éléments d’une planète, il n’y a qu’à considérer que les coordonnées étant indépendantes, doivent être prises pour les variables et, comme il n’y a qu’un seul corps mobile dont la masse peut être supposée égale à l’unité, on aura simplement, comme dans l’article_3,

donc

Ainsi les constantes et qui représentent les valeurs de et de lorsque (Sect. V, art. 12), seront ici (art. 31), et les variations des éléments deviendront de la forme

les coefficients représentés par les symboles étant exprimés ainsi :

On voit que ces expressions de coïncident avec celles que nous avons trouvées ci-dessus (art. 60), si ce n’est qu’à la place des lettres il y a les lettres qui représentent les valeurs de lorsque est égal à zéro, ou à une valeur quelconque, puisque le commencement du temps est arbitraire ; ce qui revient au même, parce que, les coefficients devant être indépendants de les quantités doivent être les mêmes fonctions de que de

64. Comme les quantités sont aussi des constantes arbitraires, on peut les prendre à la place des six constantes Changeant donc en en en en on aura

et tous les autres coefficients deviendront nuls ; de sorte que les variations de seront représentées par ces formules très simples

lesquelles résultent aussi de celles auxquelles nous sommes parvenus directement dans l’article 14 de la Section V ; ainsi il y aurait toujours de l’avantage à employer ces constantes à la place des autres constantes

Mais, quelles que soient les constantes elles ne peuvent être que des fonctions des constantes donc, réciproquement, on peut regarder celles-ci comme fonctions de celles-là. On aura ainsi

donc, substituant les valeurs de du présent article, on

aura

Or étant fonctions de on a

Substituant ces valeurs et ordonnant les termes par rapport aux variations on aura

où les symboles sont exprimés par ces formules

On aura de même, à cause de

et ainsi de suite, en changeant simplement les quantités entre elles, prises deux à deux, et en observant que l’on a, en général,

de sorte que la valeur des symboles ne fait que changer de signe par la permutation des deux quantités qu’ils contiennent.

Si l’on compare les valeurs de ces symboles marqués par des crochets carrés avec celles des symboles analogues marqués par des crochets ronds (art. 63), on y remarque une analogie singulière, qui consiste en ce qu’elles sont exprimées de la même manière en différences partielles de relatives à ou de relatives à

65. Ces dernières formules sont celles que j’avais trouvées directement, dans mon premier Mémoire sur la variation des constantes arbitraires[12], et elles résultent aussi immédiatement de la formule de l’article 12 de la Section V, laquelle, en faisant les substitutions indiquées ci-dessus (art. 61), se réduit à

Dans cette formule, les différences marquées par doivent se rapporter aux variations de toutes les constantes mais les différences marquées par peuvent se rapporter à la variation de chacune de ces constantes en particulier (Section citée, art. 10). Ainsi l’on aura, en rapportant la caractéristique successivement à

et de même en changeant en

Mais on a

et de même pour en faisant ces substitutions, on a pour les mêmes formules trouvées ci-dessus.

Mais une conséquence importante qui résulte de ces formules, c’est que la variation de la fonction en tant qu’elle dépend de celle des éléments est toujours nulle. En effet, si, dans la différentielle

on substitue les valeurs de en on trouve que tous les termes se détruisent, ce qui est un résultat très remarquable.

66. Comme la solution du problème principal dans lequel on n’a point égard aux forces perturbatrices doit donner les valeurs des variables en avec les constantes arbitraires il n’y a qu’à faire d’abord dans ces valeurs et dans celles de leurs différentielles relatives à et prendre ensuite leurs différences partielles relatives à On a ainsi facilement les coefficients des différences dans les valeurs de et il ne s’agit plus que de chercher ces différences mêmes par des éliminations linéaires, comme vje l’ai pratiqué dans le Mémoire cité, à l’égard des éléments des planètes.

À cet égard, les formules de l’article 63 paraissent avoir de l’avantage, en ce qu’elles donnent directement les mêmes différences ; mais elles demandent d’abord qu’on ait les expressions des constantes arbitraires par les variables et leurs différentielles, ce qui, dans plusieurs cas, ne peut s’obtenir que par des éliminations d’un genre supérieur aux linéaires ; ensuite, après avoir pris leurs différences partielles relatives à il faut y remettre les valeurs de ces variables en puisqu’en dernière analyse les coefficients doivent devenir des fonctions de sans ce qui constitue l’essence et la force de cette analyse.

Quoi qu’il en soit, ayant donné, dans le § I, des expressions fort simples des coordonnées en et nous y appliquerons les formules du dernier article, pour en déduire les variations des éléments comme nous l’avons pratiqué dans le Mémoire cité, parce que le calcul par ces formules acquiert une simplicité et une élégance qu’il n’aurait pas, à beaucoup près, par les autres formules.

67. Reprenons les expressions de données dans l’article 13,

dans lesquelles (art.  17)

l’angle étant déterminé par l’équation (art.  16)

Ces formules ont l’avantage que les trois éléments de l’orbite ne se trouvent que dans les quantités variables et sont, par conséquent, séparés des trois éléments qui dépendent de la position de l’orbite et dont les coefficients sont fonctions (art.  13).

Considérons d’abord la formule

et substituons-y les valeurs de données ci-dessus ; en faisant

on aura pour les mêmes expressions, où les quantités seront marquées d’un trait ; et comme les constantes n’entrent que dans et , on aura

et changeant en et on aura les valeurs de

Ces différentes valeurs étant substituées dans la formule précédente, et ayant égard aux équations de condition

qui ont lieu entre les coefficients (art. 14), cette formule se réduira à la forme

où l’on voit que les quantités, qui dépendent de la position de l’orbite, ont disparu.

On aura un pareil résultat par rapport aux différences partielles relatives à et il n’y aura qu’à changer dans la formule précédente et en

Si donc on substitue dans leurs valeurs en qu’ensuite on fasse égal à zéro ou à une quantité quelconque déterminée, et qu’on désigne par ce que deviennent, on aura (art. 64)

et l’on aura de même les valeurs de en changeant en et en dans les différences partielles.

68. Or on a

donc, puisque on aura

Mais l’équation

donne, par la différentiation,

donc on aura

Il faut maintenant différentier ces formules en faisant varier les trois constantes nous dénoterons par la caractéristique les variations relatives à ces constantes ; ainsi l’on aura d’abord

ensuite

On peut faire ici mais il est plus simple de faire ce qui donne aussi ainsi l’on aura, en changeant en

69. Ici nous avons conservé la quantité qui est l’excentricité ; mais si, à sa place, on veut employer le paramètre on aura, par la différentiation,

et les expressions de et de qui contiennent deviendront

De là nous aurons les différences partielles

et l’on trouvera, par la substitution de ces valeurs dans les expressions des symboles de l’article 67,

On aurait encore les mêmes résultats en changeant en si l’on voulait conserver l’excentricité à la place du paramètre.

70. Considérons ensuite la formule

Comme la quantité ne se trouve que dans les coefficients, qui ne contiennentpoint on a

et changeant en et en on aura les valeurs de À l’égard des différences partielles relatives à elles seront les mêmes que dans l’article précédent.

En faisant ces substitutions, on remarquera qu’en vertu des mêmes équations de condition différentiées, on aura

de sorte qu’en faisant, pour abréger,

[j’emploie l’expression différentielle [13], quoique la valeur de ne soit pas une différentielle complète], la formule

se réduira à la forme

et la formule

à cette forme semblable

Donc, en retranchant la seconde de ces quantités de la première, et observant que

on aura pour la transformée de la formule dont il s’agit, contenant les différences partielles relatives à et

et il en sera de même en changeant en et et en et

71. Il nous reste à considérer les formules où il n’y a que des différences partielles relatives à de sorte que, comme ces quantités n’entrent que dans les coefficients il n’y aura aussi que ces coefficients qui deviendront variables.

Les différentielles de ces coefficients se réduisent à une forme très simple, en employant les coefficients analogues et en ayant égard aux équations de condition entre ces différents coefficients (art. 14)).

En effet, si l’on suppose

les trois équations

étant ajoutées ensemble, après avoir multiplié la première par la deuxième par et la troisième par on a

et si on les multiplie par et par qu’on les ajoute ensuite, on a pareillement

De même les trois équations

donneront

Enfin les trois équations

donneront pareillement

72. Par le moyen de ces formules, on aura

et, affectant les quantités d’un ou de deux traits, on aura les valeurs de pur avoir celles de il n’y aura qu’à


, affecter d’un trait les quantités et Il en sera de même des différences partielles relatives à et en ne faisant que changer en et

En faisant ces substitutions et ayant toujours égard aux mêmes équations de condition, la formule

se réduira à celle-ci

Il en sera de même des formules semblables, en changeant en et

Comme les coefficients sont fonctions des trois éléments (art. 13 et 14), les trois quantités que nous avons introduites dans les formules précédentes doivent être aussi fonctions des mêmes éléments ; et si, dans les valeurs de ces trois quantités, on substitue les expressions de données dans les articles cités, on trouve, après quelques réductions fort simples,

Mais ces quantités ne servent pas seulement à simplifier le calcul. ; elles représentent d’une manière fort simple les variations instantanées de la position de l’orbite. En effet, comme le plan des auquel nous avons rapporté l’inclinaison et la longitude du nœud, est arbitraire, on peut le faire coïncider dans un instant avec celui de l’orbite, en faisant on aura alors

Dans ce cas, sera l’angle que le grand axe de l’ellipse fait avec une ligne fixe ; par conséquent, ou sera la rotation élémentaire du grand axe de l’orbite sur son plan.

L’angle élémentaire sera l’inclinaison comprise entre deux positions successives du plan de l’orbite devenue mobile, et l’angle sera la longitude du nœud formé par ces deux positions, comptée sur le même plan ; de sorte qu’en désignant par et ces deux éléments, on aura

Ainsi la variation instantanée de la position de l’orbite est déterminée par les trois éléments d’une manière indépendante de tout plan de projection.

73. Il est maintenant très facile de trouver les valeurs des autres coefficients représentés par les symboles il n’ya qu’à substituer pour c’est-à-dire pour sa valeur, qui est égale à (art. 11) ou à (art. 15), et pour leurs valeurs en de l’article précédent ; mais, à la place de l’élément nous retiendrons l’élément [14], qui exprime l’angle que le grand axe de l’orbite parcourt en tournant sur son plan mobile c’est proprement le mouvement de l’aphélie ou du périhélie sur le plan même de l’orbite. Nous aurons ainsi

donc

ensuite

et toutes les autres différences partielles seraient nulles, ce qui donne

De là on aura

74. Ces valeurs, jointes à celles que nous avons déjà trouvées (art. 69), donneront enfin

d’où résultent ces expressions très simples des variations des éléments elliptiques

75. On aurait des formules un peu moins simples, si, à la place du paramètre on voulait conserver l’excentricité Alors, à cause de on aurait

ce qui donnerait

et les valeurs de deviendraient

d’où l’on tire, en substituant pour sa valeur donnée ci-dessus,

qu’on substituera à la place des valeurs de de l’article précédent, les autres valeurs demeurant les mêmes.

Par ces formules, on peut donc avoir l’effet des forces perturbatrices sur le mouvement d’une planète, en rendant variables les quantités qui, sans ces forces, seraient constantes ; mais, quoiqu’on puisse, de cette manière, déterminer toutes les inégalités dues aux perturbations, c’est surtout pour les inégalités qu’on nomme séculaires que les formules que nous venons de donner sont utiles, parce que ces inégalités, étant indépendantes des périodes relatives aux mouvements des planètes, affectent essentiellement leurs éléments et y produisent des variations, ou croissantes avec le temps, ou périodiques, mais avec des périodes propres et d’une longue durée.

76. Pour déterminer les variations séculaires, il n’y aura qu’à substituer pour la partie non périodique de cette fonction, c’est-à-dire le premier terme du développement de en séries du sinus et du cosinus d’angles dépendants des moyens mouvements de la planète troublée et des planètes perturbatrices. Car, n’étant fonction que des coordonnées elliptiques de ces planètes, lesquelles peuvent toujours, du moins tant que les excentricités et les inclinaisons sont peu considérables, se réduire en séries du sinus et cosinus d’angles proportionnels aux anomalies et aux longitudes moyennes, on pourra aussi développer la fonction dans une série du même genre, et le premier terme sans sinus et cosinus sera le seul qui puisse donner des équations séculaires.

Désignons par ce premier terme de lequel sera une simple fonction des éléments de la planète troublée et des éléments semblables des planètes perturbatrices ; il est clair que l’élément qui est joint au temps ne s’y trouvera pas ainsi, en substituant à la place de on aura, pour les variations séculaires, les formules

77. L’équation

fait voir que le demi grand axe, ou la distance moyenne n’est sujette à aucune variation séculaire, ce cas particulier du

théorème général que nous avons démontré dans l’article 23 de la Section V ; car la quantité de cet article est la même que la quantité des articles 3 et suivants de la Section précédente, et l’on voit, par l’article 15, que l’on a

Ainsi il faut appliquer à la distance moyenne des planètes les résultats que nous avons trouvés sur la valeur de la force vive d’un système quelconque (Sect. V, § III).

La variation produit une altération dans le mouvement moyen ; car, étant l’anomalie moyenne, c’est-à-dire l’angle du mouvement moyen compté depuis le périhélie (art. 19), cette anomalie sera sujette à une variation exprimée par à cause de et si l’on ajoute la variation du lieu du périhélie dans l’orbite, on aura pour la variation séculaire de la longitude moyenne, que nous désignerons par On aura ainsi

à cause de

78. Lorsque l’excentricité est fort petite, les valeurs de et ont l’inconvénient d’avoir au dénominateur la quantité très petite Mais il est facile d’y remédier en substituant à la place de et les transformées

En effet, si l’on fait

on aura

donc, substituant les valeurs de et ,

Or, en regardant comme fonction de et comme fonction de on a

équation identique qui, par la substitution des valeurs de et donne ces deux-ci :

donc, faisant ces substitutions, on aura les équations

qu’on pourra employer à la place de celles qui donnent les valeurs de et (art. 75).

On peut faire des transformations analogues sur les dernières équations qui donnent les valeurs de et

Soit, pour cela,

on trouvera par un procédé analogue

79. Les forces perturbatricesque l’on considère dans la théorie des planètes viennent de l’attraction des autres planètes, et nous donnerons plus bas la valeur de qui résulte de cette attraction ; mais on pourrait aussi regarder comme force perturbatrice la résistance qu’elles éprouveraient de la part d’un fluide très subtil, dans lequel on les supposerait nager. Dans ce cas, en prenant pour la résistance, on ferait, comme on l’a vu dans l’article 8 de la Section II,

le fluide résistant étant supposé en repos.

Il en résultera ainsi, dans la valeur de les termes (art. 62)

On suppose ordinairement la résistance proportionnelle au carré de la vitesse, laquelle est représentée par et à la densité du milieu, que nous désignerons par ainsi les termes dus à la résistance, dans l’expression de seront

Pour évaluer la quantité il n’y a qu’à employer les formules des articles 67 et 70, en observant que la caractéristique se rapporte au temps qui n’entre que dans les valeurs de et et que la caractéristique doit se rapporter aux constantes arbitraires qui entrent dans et et dans les coefficients

On aura ainsi, en changeant en dans les expressions de

De là, en ayant égard aux équations de condition entre les coefficients (art. 14), on aura

et si l’on y substitue pour et leurs valeurs (art. 13), on aura

Donc les termes à ajouter à à raison de la résistance du milieu, seront représentés par

où il n’y aura plus qu’à substituer pour et leurs valeurs en données par les formules des articles 21, 22, en faisant attention que la caractéristique se rapporte à la variable et la aux constantes arbitraires.

CHAPITRE III.

Sur le mouvement d’un corps attiré vers deux centres fixes par des forces réciproquement proportionnelles aux carrés des distances.

80. Quoique ce problème ne puisse avoir aucune application au système du monde, où tous les centres d’attraction sont en mouvement, il est néanmoins assez intéressant du côté analytique pour mériter d’être traité en particulier avec quelque détail.

Supposons qu’un corps isolé soit attiré à la fois vers deux centres fixes par des forces proportionnelles à des fonctions quelconques des distances.

Soient, comme dans l’article 4, l’un des centres à l’origine des coordonnées, et la force attractive ; et, pour l’autre centre, supposons que sa position soit déterminée par les coordonnées parallèles aux soient, de plus, sa force attractive et la distance du corps à ce centre ; il est clair qu’on aura

c’est-à-dire, en substituant pour leurs valeurs en (art. 4),

en faisant distance des deux centres.

Il est clair que la valeur de sera la même que dans le problème du Chapitre I, mais la valeur de se trouvera augmentée du terme et comme est fonction de et fonction de ce terme donnera, dans les différentielles les termes suivants, savoir qu’il faudra, par conséquent, ajouter respectivement aux premiers membres des équations différentielles de l’article cité.

On aura donc, pour le mouvement du corps attiré vers deux centres par les forces et les trois équations suivantes :

(1)
(2)
(3)

Et si le corps était attiré en même temps vers d’autres centres, il n’y aurait qu’à ajouter à ces équations des termes semblables pour chacun de ces centres.

L’équation donnera cette quatrième équation, qui est une intégrale des précédentes,

et il est visible, en effet, que les trois équations précédentes, étant mul-

tipliées respectivement par et ajoutées ensemble, donnent une équation intégrable, et dont l’intégrale est celle que nous venons de présenter.

On tire de cette équation

valeur qui, étant substituée dans la première équation multipliée par la réduit à

Or, puisque

on aura, en faisant varier

donc, substituant cette valeur de on aura enfin

Cette équation a l’avantage de ne contenir que les deux variables et et elle indique en même temps qu’il doit y avoir une pareille équation entre et en changeant simplement et ainsi que et entre elles ; car il est indifférent de rapporter le mouvement du corps à l’un ou à l’autre des deux centres fixes, et il est clair qu’en le rapportant au centre de la force on trouverait, par une analyse semblable à la précédente,

ainsi l’on pourra, par ces deux équations, déterminer directement les deux rayons et

Je remarque maintenant qu’on peut, sans rien ôter à la généralité, supposer les deux coordonnées et du centre des forces nulles, ce qui revient à placer l’axe des coordonnées dans la ligne qui joint les deux centres. Par cette supposition, on aura et la quantité deviendra

laquelle ne contenant plus on aura donc

Par conséquent, la troisième équation différentielle se réduira à

dont l’intégrale est

étant une constante arbitraire ; d’où l’on tire

Mais on a

donc

par conséquent, en substituant cette valeur, on aura

de sorte que, connaissant et en on aura aussi en

Or, puisque et sont déjà données en et il est clair qu’on peut réduire la quatrième équation à ne contenir que et et alors elle sera nécessairement, à raison de la constante arbitraire une intégrale complète des deux équations ci-dessus en et En effet, on aura

ajoutant et réduisant, il viendra

De plus, on aura

donc, faisant ces substitutions dans la quatrième équation et ôtant le dénominateur, on aura cette intégrale

(a)

Et il est facile de voir maintenant, d’après la forme de cette équation, qu’elle résulte des deux équations en et multipliées respectivement par

ajoutées ensemble et intégrées ensuite ; mais il aurait été assez difficile de découvrir cette intégrale a priori.

81. Pour achever la solution, il faut avoir encore une autre intégrale des mêmes équations ; mais on ne saurait y parvenir que pour des valeurs particulières de et

Si l’on suppose, ce qui est le cas de la nature,

on trouve alors que ces équations, multipliées l’une par et l’autre

par donnent une somme intégrable et dont l’intégrale est
(b)

étant une nouvelle constante arbitraire.

Cette équation, étant multipliée par et ajoutée à l’intégrale (a) trouvée précédemment, donne, dans l’hypothèse présente, une réduite de la forme

(c)

Et la même équation, étant multipliée par et ensuite ajoutée à celle-ci, ou retranchée, donnera cette double équation

(d)

de sorte qu’en faisant on aura ces deux-ci

(e)

d’où l’on tire d’abord cette équation, où les variables sont séparées,

(f)

ensuite

(g)
Les mêmes substitutions étant employées dans l’expression de trouvée ci-dessus, on aura

et, substituant la valeur de

(h)

Si l’on pouvait intégrer chacune de ces différentielles, on aurait d’abord une équation entre et ensuite on aurait et en fonction de et donc on aurait et, de là, et en fonction de et comme

on aurait aussi en Mais ces différentielles se rapportant à la rectification des sections coniques, on ne saurait les intégrer que par approximation, et la meilleure méthode pour cela me paraît celle que j’ai donnée ailleurs[15] pour l’intégration de toutes les différentielles qui renferment un radical carré où la variable monte à la quatrième dimension sous le signe.

82. Si, outre les deux forces et qui attirent le corps vers les deux centres fixes, il y avait une troisième force proportionnelle à la distance, qui l’attirât vers le point placé au milieu de la ligne qui joint les deux centres, il est visible que cette force pourrait se décomposer en deux tendantes aux mêmes points, et proportionnelles aussi aux distances. Dans ce cas donc, on aurait

et l’on trouverait que l’intégrale (b) aurait aussi lieu dans ce cas ; seulement, il faudrait ajouter à son premier membre les termes

ensuite il y aurait à ajouter au premier membre de l’équation les termes

et, par conséquent, au premier membre de l’équation les termes

de sorte qu’il n’y aura qu’à augmenter les polynômes en et sous le signe, dans les équations des termes respectifs

ce qui ne rend guère la solution plus compliquée.

83. Quoiqu’il soit impossible d’intégrer en général l’équation trouvée entre et et d’avoir, par conséquent, une relation finie entre ces deux variables, on peut néanmoins en avoir deux intégrales particulières représentées par const et const.

En effet, si l’on représente en général cette équation par

il est clair qu’elle aura aussi lieu en faisant ou nuls, pourvu que les dénominateurs ou soient aussi nuls en même temps, et du même ordre.

Pour déterminer les conditions nécessaires dans ce cas, on fera

étant une constante, et une quantité infiniment petite, et désignant par ce que devient lorsqu’on change en le membre deviendra

il faudra donc, pour qu’il y ait le même nombre de dimensions de en haut et en bas, que l’on ait

alors, à cause de infiniment petit, la différentielle dont il s’agit se réduira à

dont l’intégrale est

étant une constante arbitraire. Si donc on fait et qu’on prenne en même temps aussi la valeur de deviendra indéterminée, et l’équation pourra toujours subsister, quelque valeur que puisse avoir l’autre membre Or on sait, et il est visible par soi-même, que

sont les conditions qui rendent une racine double de l’équation D’où[16] il suit, en général, que, si le polynôme a une ou plu-

sieurs racines doubles, chacune de ces racines fournira une valeur particulière de il en sera de même pour le polynôme

Maintenant il est clair que l’équation ou représente une ellipse dont les deux foyers sont aux deux centres des rayons et et dont le grand axe est égal à De même, l’équation ou représente une hyperbole dont les foyers sont aux mêmes centres, et dont le premier axe est

Ainsi les solutions particulières dont nous venons de parler donnent des ellipses ou des hyperboles décrites autour des centres des forces pris pour foyers. Et comme les polynômes et contiennent les trois constantes arbitraires dépendant de la direction et de la vitesse initiales du corps, il est visible qu’on pourra toujours prendre ces éléments tels que le corps décrive une ellipse ou une hyperbole donnée autour des foyers donnés. Ainsi la même section conique qui peut être décrite en vertu d’une force tendante à l’un des foyers et agissant en raison inverse des carrés des distances, ou tendante au centre et agissant en raison directe des distances, peut l’être encore en vertu de trois forces pareilles tendantes aux deux foyers et au centre, ce qui est très remarquable[17].

84. S’il n’y avait qu’un seul centre vers lequel le corps fût attiré par la force on aurait le cas de l’orbite elliptique que nous avons résolu dans le Chapitre I. Dans ce cas, on aurait et les deux polynômes et deviendraient semblables et ne passeraient pas le quatrième degré ; les équations de l’article 81 seraient alors intégrables par les méthodes connues, et le mouvement du corps serait déterminé par des formules en et c’est-à-dire par les distances aux deux centres, dont l’un, celui dont l’attraction est nulle, pourrait être placé où l’on voudrait ces formules ne seraient donc que de pure curiosité ; mais il y a un cas où elles se simplifient et donnent un résultat remarquable, c’est celui où le centre d’attraction nulle est placé sur le périmètre de l’ellipse.

Pour obtenir ce cas, on déterminera les constantes et de manière que, le rayon étant nul, l’autre rayon soit égal à distance entre les deux autres ; par conséquent, il faudra que les variables et deviennent à la fois égales à Les équations de l’article 81 sont très propres à cette détermination.

Faisant la première de ces équations donne

ensuite, la différence de ces équations étant divisée par si l’on y fait on a, à cause de

d’où l’on tire

Par les substitutions de ces valeurs, le polynôme

devient

ce qui se réduit à la forme

il en sera de même du polynôme en

Or, par l’article 15, on a, dans ce cas,

étant le demi grand axe de l’ellipse ; donc les équations et deviendront

et si de cette dernière on retranche la première, multipliée par et qu’ensuite on divise les numérateurs et les dénominateurs respectivement par on aura

expression qui a l’avantage de ne contenir d’autre élément que le grand axe

85. Si l’on fait

l’intégrale étant prise de manière qu’elle commence lorsque a une valeur quelconque donnée, et qu’on remette pour et leurs valeurs on aura, en intégrant,

où l’on voit que lorsque de quelque manière que l’intégrale soit prise.

Or, puisque est le rayon vecteur qui part du foyer, est le rayon qui part de l’autre centre, qui est pris dans un point de l’ellipse, et dont la distance au foyer est il est clair que et seront deux rayons vecteurs, et que sera la corde de l’arc intercepté entre ces deux rayons ; par conséquent, l’expression précédente de sera le temps employé par le mobile à décrire cet arc dans l’ellipse, lequel sera donné ainsi par la somme des rayons vecteurs par la corde et par le grand axe

L’intégrale que nous avons désignée par fonction dépend des arcs de cercle ou des logarithmes, suivant que est positif ou négatif ; mais, lorsque l’axe est très grand, cette fonction se réduit à une série très convergente : on a alors

Le premier terme donne l’expression du temps dans la parabole, et l’on a

laquelle coïncide avec celle que nous avons trouvée dans l’article 25. Le reste de la série donne la différence des temps employés à parcourir un arc de parabole et un arc d’ellipse ou d’hyperbole ayant la même corde et la même somme des rayons vecteurs.

Cette belle propriété du mouvement dans les sections coniques a été trouvée par Lambert[18], qui en a donné une démonstration ingénieuse dans son Traité intitulé : Insigniores orbitæ cometarum proprietates. Voir aussi les Mémoires de l’Académie de Berlin pour l’année 1778[19].

Le problème que nous venons de résoudre l’a été d’abord par Euler, dans le cas où il n’y a que deux centres fixes qui attirent en raison inverse des carrés des distances, et où le corps se meut dans un plan passant par les deux centres (Mémoires de Berlin de 1760) ; sa solution est surtout remarquable par l’art avec lequel il a su employer différentes substitutions, pour ramener au premier ordre et aux quadratures des équations différentielles qui, par leur complication, se refusaient à toutes les méthodes connues.

En donnant une autre forme à ces équations, je suis parvenu directement aux mêmes résultats, et j’ai même pu les étendre au cas où la courbe n’est pas dans un même plan, et où il y a, de plus, une force proportionnelle à la distance et tendante à un centre fixe placé au milieu des deux autres centres. Voir le quatrième Volume des anciens Mémoires de Turin[20] d’où l’analyse précédente est tirée, et dans lequel on trouvera aussi l’examen du cas où l’un des centres s’éloignant à l’infini, la force tendante à ce centre deviendrait uniforme et agirait suivant des lignes parallèles et il est remarquable que, dans ce cas, la solution ne se simplifie guère seulement les radicaux qui forment les dénominateurs des équations séparées, au lieu de contenir les quatrièmes puissances des variables, ne contiennent que les troisièmes, ce qui fait également dépendre leur intégration de la rectification des sections coniques.

CHAPITRE IV.

Du mouvement de deux ou plusieurs corps libres qui s’attirent mutuellement, et en particulier du mouvement des planètes autour du soleil, et des variations séculaires de leurs éléments.

86. Lorsque plusieurs corps s’attirent réciproquement avec des forces proportionnelles aux masses et à des fonctions des distances, on a, pour leurs mouvements, les formules générales des articles 1 et 2, en prenant les corps mêmes pour les centres d’attraction.

Soient les masses des corps, leurs coordonnées rectangles rapportées à des axes fixes dans l’espace ; la quantité sera, comme dans l’article 1,

Soient les distances des corps au corps et les fonctions de ces distances auxquelles les attractions entre ces corps sont proportionnelles.

Soient aussi les distances des corps au corps et les fonctions de ces distances proportionnelles aux attractions.

Soient de même les distances des corps au corps et les fonctions de ces distances proportionnelles aux attractions ; et ainsi de suite.

On aura

et la quantité (art. 2) sera

Or, quelles que soient les coordonnées indépendantes qu’on voudra adopter, on aura toujours, par rapport à chacune d’elles, comme une équation de la forme canonique

Et comme, dans le système que nous considérons, il n’y a aucun point fixe, on pourra prendre l’origine des coordonnées partout où l’on voudra, et l’on aura toujours, comme on l’a vu dans la Section III, les trois intégrales finies relatives au centre de gravité, ainsi que les trois intégrales du premier ordre relatives aux aires, et enfin l’intégrale des forces vives

On aura de cette manière le mouvement absolu des corps dans l’espace mais, comme la solution de ce problème n’est importante qu’à l’égard des planètes, et qu’il n’y a que leurs mouvements relatifs par rapport au Soleil, regardé comme immobile, qui intéressent l’Astronomie, il nous reste à voir comment on peut transporter aux mouvements relatifs l’équation générale des mouvements absolus des corps du système.


§ I. — Équations générales pour le mouvement relatif des corps
qui s’attirent mutuellement.

87. Supposons qu’on demande les mouvements relatifs des corps par rapport au corps désignons par les coordonnées rectangles du corps rapporté au corps en prenant ce dernier corps pour l’origine des coordonnées ; soient de même les coordonnées rectangles du corps par rapport au même corps et ainsi de suite ; la question consistera à trouver une formule générale qui ne contienne que ces coordonnées.

Il est d’abord évident qu’on aura

et la quantité deviendra

Comme les variables après ces substitutions, n’entrent plus dans la quantité et que ces variables n’entrent point dans sous la forme finie, on aura, relativement à ces mêmes variables, les équations

ce qui donne

étant des constantes arbitraires.

Ainsi, on aura les trois équations

les quantités étant des constantes.

Si maintenant on substitue dans l’expression précédente de les valeurs de tirées de ces équations, et qu’on fasse, pour abréger,

on aura

88. Les variables étant indépendantes, et la quantité ne contenant point ces variables sous la forme finie, on aura tout de suite, par rapport à chacune d’elles, les équations

Si l’on ajoute ensemble les premières équations relatives aux variables on a

ce qui donne

et l’on trouvera de même, par l’addition des deuxièmes équations et par celle des troisièmes,

valeurs qu’on pourra substituer dans les équations précédentes.

On aura ainsi, pour le mouvement du corps autour de les trois équations

et l’on aura de pareilles équations pour le mouvement des corps autour du corps en changeant seulement entre elles les quantités affectées de deux traits, ou de trois, etc.

Il n’y aura donc qu’à substituer la valeur de et prendre ses différences partielles relatives aux différentes variables ; mais cette substitution peut se simplifier par la considération suivante.

89. Dénotons par la somme de tous les termes de la quantité qui contiennent les distances et remarquons que les expressions de ces distances sont telles, qu’elles demeurent les mêmes en augmentant les coordonnées qui y entrent d’une même quantité quelconque ; d’où il suit qu’en faisant varier ces mêmes coordonnées d’une même quantité infiniment petite, la variation de sera nulle, ce qui donnera l’équation

On trouvera de la même manière, parce que la même propriété a lieu par rapport aux coordonnées et aux coordonnées

Donc, puisque

ne contenant que ne contenant que et ainsi de suite, la première équation deviendra par ces substitutions, en la divisant par

Or, dans la quantité il n’y a que les termes qui contiennent qui dépendent des variables (art. 86) ; ainsi l’on peut réduire la valeur de à

substituant la valeur de dans l’équation précédente, elle deviendra

et l’on aura de la même manière

90. On peut ramener ces équations à la forme générale, qui a l’avantage de s’appliquer également à des coordonnées quelconques.

Si l’on multiplie la première par la deuxième par la troisième par et qu’on les ajoute ensemble, on aura d’abord la partie différentielle

laquelle, en transformant les coordonnées en d’autres coordonnées indépendantes donnera pour les termes multipliés par la formule (Sect. IV, art. 7)

en faisant

À l’égard des termes qui contiennent les forces il est facile de voir qu’en changeant la caractéristique en tous ces termes sont intégrables par rapport aux variables l’intégrale contiendra d’abord les termes


ensuite elle contiendra les termes

Or on a

et comme

On aura

et de même

et ainsi des autres expressions semblables. Donc, nommant l’intégrale totale, on aura

Ainsi, après la transformation des coordonnées, les termes multipliés par se réduiront à et, comme on suppose les nouvelles coordonnées indépendantes, chacune d’elles, comme donnera une équation de la forme

91. S’il n’y a que deux corps, et l’expression de devient

ainsi les valeurs de et de sont les mêmes que pour un corps attiré vers un centre fixe avec une force proportionnelle à une fonction de la distance (art. 4). Donc le mouvement relatif du corps autour du corps sera le même que si celui-ci était fixe et que la masse attirante fût la somme des deux masses, ce qui est connu depuis Newton.

Lorsque la masse du corps autour duquel les autres sont censés se mouvoir est beaucoup plus grande que la somme des masses ce qui est le cas du Soleil par rapport aux planètes, on a, à très peu près,

Le mouvement du corps autour du corps sera donc, dans ce cas, à très peu près le même que si celui-ci était fixe et que la somme des masses y fût réunie ; et en regardant les autres forces comme des forces perturbatrices, on pourra employer la théorie de la variation des constantes arbitraires pour déterminer l’effet de ces forces il ne s’agira que de prendre, conformément à l’article 9 de la Section V, égale à la somme de tous les autres termes de la valeur de donnée ci-dessus. On fera ainsi, en accentuant la lettre pour la rapporter à la planète

et l’on aura, par les formules générales de l’article 14 de la même Section, les variations des éléments du mouvement du corps autour du corps regardé comme fixe.

§ II. — Formules générales pour les variations séculaires des éléments des orbites des planètes autour du Soleil.

92. Pour appliquer ces formules au mouvement des planètes autour du Soleil, on prendra la masse pour celle du Soleil, la masse pour celle de la planète dont on cherche les perturbations, et les masses pour les masses des planètes perturbatrices, et l’on fera

On substituera ensuite dans la fonction au lieu des coordonnées de ces différents corps autour de leurs valeurs exprimées en fonction de conformément aux formules que nous avons données dans le Chapitre Ier pour les expressions des coordonnées en y faisant

ou simplement pour le rapporter au corps et l’on aura, par les articles 69 et suivants, les variations des six éléments de l’orbite de la planète autour du Soleil.

Nous nous contenterons ici de chercher les variations séculaires de ces éléments qui sont les plus importantes, et qui ne dépendent que du premier terme tout constant du développement de

L’expression de deviendra

93. Commençons par développer la quantité

on y mettra d’abord pour les expressions de de l’article 13, en marquant par un trait, ou par deux, les quantités qui se

rapportent aux masses On aura

Si l’on fait les multiplications, qu’on développe les produits des sinus et cosinus et qu’on fasse, pour abréger,

on aura

Les quantités sont fonctions des éléments de l’orbite de la planète donnée par les formules de l’article 13, en marquant toutes les lettres d’un trait, et les quantités sont fonctions semblables des éléments de l’orbite de la planète en marquant les lettres de deux traits ainsi les quantités sont fonctions de ces mêmes éléments ; mais, par la considération suivante, on peut voir ce qu’elles expriment.

94. Les coordonnées primitives rapportées à un plan donné étant pour les transformer dans les coordonnées rapportées au plan de l’orbite de on a, par les formules générales de l’article 14,

où les coefficients dépendent des constantes qui déterminent la position des nouveaux axes par rapport aux primitifs, étant l’inclinaison des deux plans.

De même, si l’on voulait transformer les mêmes coordonnées dans les coordonnées rapportées au plan de l’orbite de on aurait

où les coefficients seraient des fonctions semblables des constantes qui déterminent la position de ce nouveau plan patrapport au même plan primitif, et où serait l’inclinaison de ces plans.

Si l’on compare maintenant ces expressions, on aura

Comme les coefficients sont assujettis aux mêmes équations de condition que les coefficients de l’article 14, si l’on ajoute ensemble les trois équations précédentes, après les avoir multipliées respectivement par par et par on aura, en vertu de ces équations,

en faisant, pour abréger,

Il est évident que, par ces formules, les coordonnées sont transformées dans les coordonnées ainsi les coefficients seront exprimés d’une manière semblable aux coefficients analogues et, prenant les constantes à la place des on aura, par les formules générales de l’article 13,

La constante représentera l’angle d’inclinaison des deux plans où sont placées les orbites des planètes et et nous la désignerons par pour indiquer qu’elle se rapporte aux orbites de et et si, dans l’expression de en on substitue les valeurs de ces coefficients en (art. 13), on a

On voit que les quantités que nous venons de désigner par sont les mêmes fonctions de que celles que nous avons désignées par les mêmes lettres dans l’article 93 ; ainsi l’on aura dans les formules de cet article, en y substituant pour ces quantités les valeurs que nous venons de trouver,

95. Donc, faisant ces substitutions dans l’expression de de l’article 93, on aura

Faisons, pour un moment,

on aura

c’est la valeur qu’il faudra substituer dans l’expression de de l’article 92 ; et l’on aura de la même manière

En marquant de trois traits les lettres qui ne sont marquées que de deux, on aura les termes multipliés par dans et ainsi de suite.

Il faudra ensuite substituer pour et pour leurs valeurs exprimées par les anomalies moyennes suivant les formules des articles 21 et 22 ; et, dans le développement, nous nous contenterons d’avoir égard aux secondes dimensions des excentricités et des inclinaisons mutuelles des orbites de sur celles de en regardant ces quantités comme très petites du même ordre, et en négligeant les termes où elles formeraient des produits de plus de deux dimensions.

On aura ainsi

96. On sait que les puissances d’une fonction de la forme

peuvent se développer en séries de cosinus d’angles multiples de ainsi l’on peut supposer

sont des fonctions de exprimées en séries ou par des intégrales définies, dans lesquelles les quantités et entrent de la même manière et forment des fonctions homogènes des dimensions ou

Ainsi, en faisant

on aura

où il faudra faire (art. 21 et 22)

et, par conséquent,

étant égal à

Comme nous négligeons les quantités au-dessus du second ordre dans les termes multipliés par on pourra mettre tout de suite et au lieu de au lieu de et au lieu de et, en développant les produits des cosinus, on verra facilement qu’il n’y aura de terme indépendant des angles et que celui-ci :

Considérons maintenant les fonctions en y faisant les substitutions précédentes à la place de et de et conservant les secondes dimensions de et de on aura

et il en sera de même des autres fonctions semblables. On aura pareillement

on développera de la même manière les cosinus des angles multiplcs de et l’on multipliera les expressions de ces cosinus par celles des coefficients

Comme nous ne cherchons que les termes indépendants de toute période, il faudra rejeter tous ceux qui se trouveront multipliés par des cosinus d’angles multiples de et de qui sont les angles des mouvements moyens de et de Ainsi le terme de l’expression ne donnera que ceux-ci :

Le terme donnera ceux-ci :

Le terme et les suivants ne donneront que des termes où formeraient plus de deux dimensions, et que nous négligeons.

Il reste à développer la quantité (art. 93, 94)

On fera ici pour les mêmes substitutions que ci-dessus ; on aura d’abord

mais, dans le terme qui contient et qui est déjà du second ordre, il suffira de mettre à la place de et comme

il est clair que ce terme ne donnera aucune quantité constante.

En multipliant l’expression précédente de par celle de de l’article précédent, et ne retenant que les termes constants où et ne passent pas la seconde dimension, on trouve aisément ceux-ci :

de sorte que, dans la quantité dont il s’agit, les termes constants se détruisent.

97. La somme de tous les termes que nous venons de trouver, étant multipliée par sera la partie constante de la fonction due à l’action de la planète et l’on aura une expression semblable pour la partie due à l’action de la planète en rapportant à celle-ci les quantités relatives à la planète

Nous avons désigné par cette partie non périodique de la fonction si donc on fait, pour abréger,

et par conséquent, puisque et entrent de la même manière dans la fonction

et, de plus

on aura

Cette valeur est exacte aux quantités du troisième ordre près, en regardant les excentricités et des orbites de et de ainsi que leur inclinaison mutuelle comme de très petites quantités du premier ordre, quelles que soient d’ailleurs les inclinaisons de ces orbites sur le plan fixe de projection.

98. On peut simplifier beaucoup les expressions des fonctions et par les propriétés connues des coefficients des séries en En effet, si l’on différentie logarithmiquement par rapport à et ensuite par rapport à l’équation identique

et qu’après avoir multiplié en croix, on compare les termes multipliés par les mêmes cosinus, on aura d’abord

ensuite les différentielles relatives à et donneront

Substituant ces valeurs, on aura

Mais on peut avoir des expressions plus simples de ces fonctions, en employant les coefficients de la série

car, en différentiantlogarithmiquement et multipliant ensuite en croix, on trouve d’abord, comme ci-dessus,

substituant cette valeur de et comparant la série multipliée par avec la série

avec laquelle elle doit devenir identique, il est facile de déduire ces relations

et, par la substitution de ces valeurs, on aura celles-ci

qu’on substituera dans l’expression de de l’article précédent.

À l’égard de la valeur des coefficients en fonction de on peut les trouver par le développement des radicaux en puissances de et par le développement de ces puissances en cosinus d’angles multiples de comme Euler l’a fait le premier dans ses Recherches sur Jupiter et Saturne ; mais j’ai trouvé, depuis longtemps, qu’on pouvait les avoir d’une manière plus simple en décomposant le trinôme en ses deux facteurs imaginaires

et en développantpar la formule du binôme les puissances de chacun de ces facteurs.

Soit, pour abréger,

on aura, en général

et si l’on multiplie ensemble les deux séries qui répondent à et à et qu’on repasse des exponentielles imaginaires aux cosinus des angles multiples, on aura

en faisant

Ces séries sont toujours convergentes, lorsque mais si l’on avait il n’y aurait qu’à changer en et en puisque, dans la fonction non développée, les quantités et entrent également.

Une conséquence qui résulte de la forme de ces séries est que, tant que est un nombre positif, tous les coefficients ont toujours des valeurs positives.

Si l’on fait ces coefficients deviendront et, si l’on fait ils deviendront

99. Il nous reste encore à déterminer l’angle Comme nous négligeons les quantités du troisième ordre etque, dans l’expression de est déjà multiplié par on pourra, dans la détermination de l’angle faire abstraction des quantités très petites du premier ordre, et par conséquent y supposer Or (art. 96)

et, faisant dans les formules de l’article 94, on a

on a aussi, par les formules de cet article,

Différentions cette valeur de faisons varier les quantités substituons à la place de leurs différentielles les expressions données dans l’article 71, en accentuant les quantités respectives, et remettons les quantités à la place de leurs valeurs en on trouvera facilement

Mais

donc, en divisant par on aura

et, intégrant,

où il n’est pas nécessaire d’ajouter des constantes, puisque l’origine des angles est arbitraire. L’angle [21] est en général celui que l’orbite décrit en tournant dans son plan, et que nous avons substitué à la place de la longitude du périhélie (art. 68).

100. La fonction est maintenant réduite à la forme la plus simple et la plus propre pour le calcul des variations séculaires ; il n’y aura qu’à la substituer dans les formules de l’article 76, en marquant d’un trait les lettres de ces formules, pour les rapporter à la planète dont on cherche les variations ; et, en changeant simplement entre elles les lettres marquées d’un trait et de deux, on aura des formules semblables pour les variations de la planète et ainsi des autres.

On voit que cette fonction est composée de deux fonctions distinctes entre elles, dont l’une ne renferme que les excentricités et les lieux des aphélies dans les orbites, et dont l’autre ne renferme que les inclinaisons des orbites sur un plan fixe avec les lieux de leurs nœuds. Si l’on désigne la première par et la seconde par en sorte que l’on ait

on aura

les angles étant les inclinaisons de l’orbite de la planète sur celles des planètes

On substituera ainsi au lieu de dans les équations des variations séculaires (art. 76), et l’on accentuera les lettres, pour les rapporter à la planète dont on cherche les variations ; on aura, en négligeant les et mettant simplement au lieu de dans les coefficients des fonctions et qui sont déjà du second ordre,

On aura de pareilles équations pour les variations des éléments de la planète dans son orbite autour de et il n’y aura pour cela qu’à marquer de deux traits les lettres qui ne sont marquées que d’un trait, et au contraire ne marquer que d’un trait les lettres qui le sont de deux.

Ainsi, en observant que les fonctions de et représentées par des parenthèses ne changent pas en échangeant entre elles les quantités on aura

et les équations des variations seront

et ainsi de suite pour les variations des éléments des orbites de autour de

101. Mais nous remarquerons qu’on peut réduire à une seule fonction les différentes fonctions ainsi que les fonctions ce qui mettra plus de simplicité et d’uniformité dans les formules des variations : en effet, si l’on fait

en faisant toutes les combinaisons deux à deux des masses

et des fonctions qui y sont relatives, il est facile de voir que, dans les différences partielles de on pourra changer ces fonctions en pourvu qu’on divise par les différences partielles relatives à et par les différences partielles relatives à et ainsi de suite.

De sorte que les équations des variations des excentricités et des aphélies deviendront

Ces équations donnent

donc, comme est fonction des variables sans on aura

et par conséquent égale à une constante. C’est une relation générale entre les excentricités et les lieux des aphélies des planètes, qui doit toujours subsister, quelques variations que les excentricités et les lieux des aphélies subissent à la longue, pourvu qu’elles soient très petites.

102. Mais la nature de la fonction donne encore naissance à d’autres relations générales entre ces mêmes éléments.

En effet, il est facile de voir qu’on a l’équation

et si l’on substitue à la place de ces différences partielles leurs valeurs données par les équations de l’ar-

ticle précédent, on aura, en intégrant relativement à l’équation finie

étant une constante égale à la valeur du premier membre de cette équation dans un instant quelconque.

Cette équation fait voir que les excentricités ont nécessairement des limites qu’elles ne peuvent passer ; car, comme elles sont nécessairement réelles, tant que les orbites sont des sections coniques, chaque terme, comme sera toujours positif, et son maximum sera la constante

Il suit de là que, si les excentricités des orbites qui appartiennent à des masses très grandes sont une fois très petites, elles le seront toujours, ce qui est le cas de Jupiter et Saturne ; mais celles qui appartiennent à des masses fort petites pourront croître jusqu’à l’unité et au delà, et l’on ne pourra déterminer leurs véritables limites que par l’intégration des équations différentielles, comme on le verra ci-après[22].

De plus, comme la quantité regardée comme fonction de est une fonction homogène de deux dimensions, on aura, par la propriété connue de ces fonctions,

Substituant dans cette équation les valeurs des différences partielles de relatives à tirées des mêmes équations de l’article précédent, on aura

étant la valeur de dans un instant quelconque.

Dans cette équation, les quantités expriment les vitesses angulaires des mouvements des aphélies, et par conséquent elle donne une relation invariable entre ces vitesses, par laquelle on voit qu’elles ont aussi nécessairement des limites, tant qu’elles sont toutes de même signe.

103. Si l’on emploie les transformations de l’article 78, en faisant

on aura

et les équations des variations seront

Si, dans ces équations, on substitue la valeur de et qu’on exécute les différentiations partielles, on a des équations linéaires en faciles à intégrer, et ces équations seront entièrement identiques avec celles que j’avais trouvées par une autre voie, dans les Mémoires de Berlin de 1781, page 262[23], comme il est facile de s’en assurer en comparant entre elles les dénominations différentes des mêmes quantités.

Dans les Mémoires de l’année 1782[24], j’ai appliqué ces équations aux six planètes principales, en adoptant pour leurs masses les valeurs les plus probables, et j’en ai tiré, par l’intégration, des formules générales pour les variations de leurs excentricités et des lieux de leurs aphélies, lesquelles donnent les valeurs de ces éléments, tant pour la Terre que pour les autres planètes, au bout d’un temps quelconque indéfini, soit avant ou après l’époque de 1700 ; et comme, par ces formules, les excentricités demeurent toujours très petites, ainsi qu’on l’a supposé dans le calcul, on est assuré de leur exactitude dans tous les temps passés et à venir. On trouve ensuite, dans le Volume des mêmes Mémoires pour 1786-87, imprimé en 1792, un supplément[25] à cette théorie, relatif à la nouvelle planète d’Herschel, où l’on détermine de la même manière, et par des formules aussi générales, les variations séculaires de l’excentricité et du lieu de l’aphélie de cette planète, produites par les actions de Jupiter et de Saturne ; on a seulement négligé l’effet de l’action d’Herschel sur ces deux planètes, ainsi que sur les autres planètes inférieures, à cause de la petitesse de sa masse et de sa grande distance.

104. On peut réduire de la même manière à une forme plus simple les équations de la variation des nœuds et des inclinaisons. Soit

en faisant aussi toutes les combinaisons deux à deux des masses qui sont supposées agir les unes sur les autres, avec les fonctions correspondantes de et des inclinaisons

qui sont déterminées en général par la formule

on aura, en substituant à la place de


Ces équations donnent aussi

et, par conséquent, puisque est une fonction de sans aucune autre variable,

et égal à une constante.

De plus, il est visible, par la forme de la fonction que l’on a cette équation

Substituant, pour les différentielles de relatives à leurs valeurs données par les équations précédentes, on aura une équation différentielle en dont l’intégrale sera

et qu’on pourra mettre aussi sous la forme

étant la valeur du premier membre dans un instant quelconque. On peut tirer de cette équation, relativement aux liniités des quantités

des conséquences analogues à celles que nous avons déduites d’une équation semblable en dans l’article 101.

105. Dans le cas où l’on ne considère que l’action de deux planètes et l’expression de se réduit au seul terme multiplié par et l’inclinaison mutuelle des deux orbites devient alors constante ; c’est à très peu près le cas de Jupiter et Saturne.

Dans ce cas, nous remarquerons encore que, si l’on suppose que le plan de la planète perturbatrice coïncide dans un instant avec le plan fixe, on aura et, par conséquent, ce qui donnera

et de là

c’est l’expression de la vitesse du mouvement rétrograde du nœud de l’orbite de sur le plan de l’orbite de tandis que leur inclinaison mutuelle deméure constante ; d’où l’on voit que l’action de la planète sur la planète pour faire varier la position de son orbite, se réduit à donner au nœud de son orbite, sur celle de la planète perturbatrice un mouvement rétrograde instantané exprimé par

sans affecter l’inclinaison mutuelle des deux orbites.

De la même manière, l’action de la planète sur la planète pour changer la position de son orbite, fait rétrograder le nœud de cette planète sur le plan de l’orbite de d’un mouvement instantané

et ainsi des autres planètes.

En combinant ainsi deux à deux toutes les planètes, on pourra trouver les variations de leurs nœuds et de leurs inclinaisons réciproques, puisque, par la nature du Calcul différentiel, la somme des valeurs particulières d’une différentielle en forme la valeur complète. C’est ainsi qu’on avait trouvé les changements annuels des nœuds et des inclinaisons des planètes, produits par leurs attractions mutuelles, avant qu’on eût une théorie directe et générale des variations séculaires.

106. Pour donner un exemple de cette méthode, considérons trois planètes dont les orbites s’entrecoupent, étant les inclinaisons de la seconde et de la troisième sur la première, et étant l’inclinaison de la seconde sur la troisième ; il est facile de voir qu’elles formeront sur la sphère un triangle sphérique dont les trois angles, en supposant les inclinaisons de et de du même côté, seront et que nous désignerons, pour plus de simplicité, par

La planète fera rétrograder sur son orbite le nœud de la planète de la quantité élémentaire

et la même planète fera rétrograder en même temps sur son orbite le nœud de l’orbite de la planète de la quantité élémentaire

tandis que les inclinaisons et demeureront constantes.

Ainsi, dans le triangle formé par l’intersection des trois orbites, la portion de l’orbite de interceptée entre les orbites de et de c’est-à-dire le côté adjacent aux angles et , croîtra de la quantité faisant, pour abréger,

les angles et demeurent constantes.

Or, dans un triangle sphérique dont les angles sont et dont le côté adjacent à et et par conséquent opposé à est on a

Donc, faisant varier de on aura

Mais la même équation donne

d’où l’on tire

Faisons, pour abréger,

on aura

On trouvera de la même manière, en considérant la rétrogradation des orbites de et de sur celle de laquelle augmente le côté adjacent aux angles et par conséquent opposé à l’angle de la quantité élémentaire en faisant

tandis que les angles et demeurent constants,

parce que la quantité est une fonction symétrique des trois cosinus.

Enfin la rétrogradation des orbites de et sur celle de donnera aussi

en faisant

Dans ces équations, les trois coefficients sont constants ; ainsi il n’y a de variable que la quantité qui est fonction des trois cosinus de c’est-à-dire des inclinaisons respectives des orbites ainsi on pourra déterminer leurs valeurs en fonction de

Si l’on ajoute ensemble ces trois équations, après avoir multiplié la première par la seconde par la troisième par on a

Or, en substituant les valeurs de on voit que le second membre se réduit à zéro, par la destruction mutuelle de tous les termes, et, le premier membre étant intégrable, on a, en remettant pour leurs valeurs

équation qui s’accorde, dans le cas de trois orbites, avec l’intégrales trouvée plus haut (art. 104).

Si l’on fait, pour plus de simplicité,

on a les trois équations

La première, combinée avec la seconde et avec la troisième, donne, par l’élimination de

et de là

Substituant ces valeurs dans l’expression de la variable montera au troisième degré sous le radical, et l’équation donnera

équation où les variables sont séparées, mais dont le second membre ne sera intégrable que par la rectification des sections coniques.

Mais, comme les inclinaisons mutuelles des orbites doivent être supposées très petites, si l’on fait

ce qui donne

les quantités devront être très petites, et l’on pourra négliger, dans l’expression de leurs troisièmes dimensions vis-à-vis des secondes. On aura ainsi

et

Si l’on différentie cette valeur de et qu’après la substitution des valeurs de on divise l’équation par qu’ensuite on la redifférentie et qu’on y fasse encore les mêmes substitutions, on aura, étant constant,

équation intégrable par des exponentielles ou des sinus, suivant que le coefficient de sera positif ou négatif ; mais, comme

il est évident que la valeur de en ne peut pas contenir d’exponentielles désignant donc par le coefficient de dans l’équation

précédente, on aura

et étant deux constantes arbitraires qu’il faudra déterminer par l’état initial ; et, comme et sont, par hypothèse, des quantités très petites, aura aussi une valeur très petite.

De là on aura, par l’intégration, les valeurs de qui ne contiendront que dans et qui, étant une fois très petites, le seront toujours nécessairement, de sorte que la solution sera toujours bonne.

On connaîtra donc ainsi les inclinaisons réciproques des orbites pour un temps quelconque, mais on ne connaîtra pas encore par là leurs positions absolues dans l’espace, lesquelles dépendent des angles c’est pourquoi il est plus simple de chercher directement ces angles par l’intégration des formules de l’article 104.

107. Mais, au lieu d’employer ces équations sous la forme où elles se présentent, il sera plus avantageux de les transformer par les substitutions de l’article 78, en faisant

on aura ainsi, en accentuant les lettres pour les rapporter successivement aux planètes et mettant la fonction \frac{\Psi}{\mathrm m} à la place de (art. 104),

La fonction sera, comme dans l’article cité,

mais les valeurs deviendront, par les mêmes substitutions,

Faisant ces substitutions et exécutant les différentiations relatives à

· · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · ·

108. Ces équations ont lieu, quelles que soient les valeurs des variables parce que nos formules ne supposent point que les inclinaisons des orbites sur le plan fixe soient très petites, comme on l’a vu jusqu’ici dans toutes les formules que l’on a données pour le mouvement des nœuds et les variations des inclinaisons, mais elles supposent seulement la petitesse des inclinaisons mutuelles des orbites.

Quant à leur intégration, elle paraît très difficile en général, et il n’y a peut-être que le cas de deux orbites où elle réussisse.

Faisons dans ce cas, pour abréger,

on aura les équations

Elles donnent d’abord

d’où l’on tire

et étant des constantes.

Maintenant, si l’on différentie l’équation et qu’on y substitue les valeurs de on a, après la division par

on trouve de même

et, de là,

étant une constante.

Les intégrales que nous venons de trouver donnent

ces valeurs étant substituées dans les trois équations

on obtient les équations suivantes

qui, étant linéaires et à coefficients constants, sont toujours intégrables.

On aura de pareilles équations en changeant en mais, lorsqu’on connaîtra les trois premières de ces quantités, on aura les trois dernières par les trois intégrales précédentes.

Les expressions de en contiendront trois constantes arbitraires, et, comme les constantes sont aussi arbitraires, on aura en tout six constantes arbitraires, mais qui se réduiront à quatre, pour satisfaire aux équations supposées

on aura ainsi les valeurs complètes des quatre variables qui donnent la position des deux orbites dans l’espace.

Mais notre analyse est fondée sur la supposition que l’inclinaison mutuelle des deux orbites soit très petite ; le cosinus de cette inclinaison est exprimé par la formule (art. 107)

dont la différentielle devient égale à zéro, d’après les équations différentielles ci-dessus ; cette quantité sera donc égale à une constante,

comme nous l’avons déjà trouvé (art. 105), et il faudra que cette constante soit supposée fort peu différente de pour l’exactitude de la solution précédente.

Il serait difficile, peut-être même impossible, de résoudre de la même manière le cas de trois ou d’un plus grand nombre d’orbites mais nous observerons que, comme la position du plan de projection est arbitraire, on peut toujours le prendre de manière que les inclinaisons des orbites sur ce plan soient très petites, puisque leurs inclinaisons mutuelles doivent être très petites ; et si, les inclinaisons étant très petites dans un instant quelconque, elles demeurent toujours très petites, la solution fondée sur cette hypothèse sera légitime.

109. En supposant les inclinaisons des orbites sur le plan fixe très petites, les variables seront aussi très petites, et l’on pourra, dans les équations de l’article 107, entre ces variables, mettre simplement à la place des radicaux

ce qui les réduira à la forme linéaire, dont l’intégration est facile.

On aura ainsi des équations tout à fait semblables à celles que j’avais trouvées, par une autre méthode, dans les Mémoires de l’Académie de Berlin pour 1782, et dont j’ai fait aussi l’application aux six planètes principales, en donnant les expressions finies des variables pour un temps indéfini ; et les Mémoires de 1787, de la même Académie, renferment, de plus, ce qui est relatif à l’orbite d’Herschel. Il faut seulement remarquer que, dans les formules de ces Mémoires, les tangentes des inclinaisons tiennent lieu des sinus qui se trouvent dans les valeurs des variables

mais, à cause de la petitesse des inclinaisons, cette différence n’est d’aucune considération.

Lorsque l’on connaît les valeurs de ces variables, on peut déterminer tout de suite les inclinaisons mutuelles des orbites par les formules de l’article 107 ; mais ces formules se simplifient dans le cas où les quantités sont très petites. Dans ce cas, on aura, en négligeant les troisièmes dimensions de ces quantités,

et de là, à cause de

et ainsi des autres.

110. Il reste encore, pour compléter la théorie des variations séculaires, à considérer la variation du mouvement moyen que nous avons désigné par dans l’article 77, et qui, en négligeant le carré de l’excentricité qui est supposée fort petite vis-à-vis de l’unité, et accentuant les lettres pour les rapporter respectivement aux planètes devient

pour la planète on aura de même, pour la planète la variation en ajoutant un trait aux lettres qui n’en ont qu’un, et ainsi de suite.

On substituera donc dans cette formule, au lieu de la fonction la somme suivant l’article 100, et comme la fonction ne contient point les excentricités on aura simplement

et, pour avoir une formule uniforme pour toutes les planètes il n’y aura qu’à mettre, suivant les remarques des articles 101

et 104, à la place de et à la place de ce qui donnera

les fonctions et étant données dans les mêmes articles et étant les mêmes pour toutes les planètes.

Mais si, à la place de ces fonctions exprimées en on veut employer les expressions des articles 103, 107, en il faudra, suivant les formules de l’article 73, changer en On aura ainsi

et, pour avoir il n’y aura qu’à changer en

Dans ces formules, les différentielles relatives à n’affectent que les coefficients des fonctions et à la place desquels il suffira de substituer

pour avoir les valeurs de et, par les formules de l’article 98, on trouvera les valeurs des différences partielles

Il faudra ensuite y substituer les valeurs de en trouvées par l’intégration des équations différentielles des articles 103 et 109, et que nous avons données, pour toutes les planètes, dans les Mémoires cités de l’Académie de Berlin ; et, comme ces valeurs sont exprimées par des suites de sinus et de cosinus, les variations seront intégrables : les termes constants donneront, dans des termes proportionnels à qui se confondront avec les mouvements moyens ; et les termes en sinus et cosinus donneront de pareils termes, qui exprimeront les variations séculaires de ces mouvements.

J’avais trouvé, par une autre méthode, dans les Mémoires cités de l’Académie de Berlin, des formules pour déterminer les variations séculaires des moyens mouvements des planètes, et elles m’avaient donné, pour Jupiter et Saturne, des résultats presque insensibles ; mais les formules précédentes sont peut-être plus rigoureuses, et il sera bon d’en faire l’application aux planètes ; c’est un objet dont je pourrai m’occuper ailleurs ici, je n’ai eu en vue que de montrer l’usage de la nouvelle théorie des variations des constantes arbitraires dans la détermination des variations séculaires des éléments des orbites des planètes.

§ III. — Sur les équations séculaires des éléments des planètes,
produites par la résistance d’un milieu très rare.

§ 111. Pour ne rien laisser à désirer sur les variations séculaires des planètes, nous devons encore considérer l’effet d’un milieu peu résistant dans lequel il est possible qu’elles se meuvent, et où elles devraient nécessairement se mouvoir, si la lumière était due aux oscillations d’un fluide.

Nous avons déjà vu, dans l’article 79, que, pour avoir égard à la résistance, il suffit d’ajouter à la valeur de les termes

étant la densité du milieu, qui peut être une fonction de et qui doit être supposée très petite, et d’y substituer pour et leurs valeurs en données par le mouvement elliptique de la planète, en se souvenant que la caractéristique se rapporte au temps et la caractéristique aux éléments de la planète.

Comme nous ne cherchons ici que les variations séculaires, il faudra, ainsi que nous l’avons pratiqué plus haut, rejeter tous les termes périodiques et ne retenir que les termes constants.

112. En désignant, comme dans l’article 21, l’anomalie moyenne de la planète par

on a vu que et peuvent s’exprimer par des séries dont la première ne contient que des cosinus et la seconde des sinus d’angles multiples de donc ne contiendra que des sinus sans terme constant, et des cosinus de ces mêmes angles ; par conséquent, la quantité ne contiendra que des cosinus, et il en sera de même de la série qui exprimera la valeur de Ainsi, en faisant la quantité sera exprimée par une série de cosinus de multiples de

Maintenant, pour avoir et il faudra faire varier dans les séries de et de les coefficients des cosinus et des sinus qui sont donnés en et et, de plus, l’angle à raison des constantes et qu’il renferme. Désignons par et les parties de et de qui contiennent les variations des coefficients ; on aura

et, de même,

donc

Or il est clair que ne contiendra que des cosinus ; et, comme ne contient que des sinus, ne contiendra aussi que des sinus sans terme constant. De même, ne contiendra que des sinus, et, comme ne contient que des cosinus, ne contiendra aussi que des sinus. D’ailleurs, ne contient que des cosinus ; par conséquent, ne contiendra que des sinus. Donc la quantité

ne contenant que des sinus d’angles multiples de sans aucun terme constant, devrà être négligée.

113. Ainsi, pour les équations séculaires, on aura simplement

Or donc, faisant ces substitutions, on aura, pour les termes à ajouter à à cause de la résistance du milieu,

où il faudra substituer pour et leurs valeurs en ou et ne retenir dans les résultats que les termes indépendants des sinus et cosinus de

Par les propriétés du mouvement elliptique, on a tout de suite

ainsi les termes dont il s’agit se réduiront à

114. Tels sont les termes qu’il faudra substituer à la place de dans les formules générales des variations des éléments des planètes (art. 74), et l’on aura

et les variations des autres éléments, seront nulles ; d’où l’on peut d’abord conclure que le grand axe ou la ligne des apsides, ainsi que le nœud et l’inclinaison, ne seront sujets à aucune variation séculaire par conséquent, la résistance ne déplacera point l’orbite de la planète, mais en changera seulement à la longue le grand axe et l’excentricité, et produira en même temps une équation séculaire dans l’anomalie moyenne et dépendante de la variation de

Si l’on combine ensemble les deux premières équations, on a

donc

divisant par et intégrant, on trouve

et, comme on aura

en supposant que l’intégrale commence au point où ainsi tout dépend de la variation de la distance moyenne

Si, dans la première approximation, on néglige l’excentricité supposée très petite, on a ce qui donne

et, comme la densité du milieu ne peut être qu’une fonction de elle scra une fonction de et la première équation donnera

Supposons constant ; on aura, en intégrant,

étant la valeur de lorsque donc

et la valeur de deviendra

où le coefficient doit être supposé très petit.

115. Pour avoir la variation séculaire de l’excentricité il faudra substituer, dans les expressions irrationnelle et de la valeur de l’expression de en et ne retenir dans le développement que les termes constants. Or, en ne conservant que les secondes dimensions de on a, par l’article 21,

d’où l’on tire

par conséquent, en rejetant les et [26], on aura

Si l’on substitue pour sa valeur et qu’on néglige d’abord les on aura l’équation

dont l’intégrale donne

étant la valeur de lorsque

Mais, comme l’existence d’un milieu résistant et, à plus forte raison, la loi de la densité de ce milieu ne sont qu’hypothétiques, les résultats précédents ne doivent être regardés que comme une application de nos formules générales des variations séculaires.

§ IV. — Du mouvement autour du centre commun de gravité
de plusieurs corps qui s’attirent mutuellement
.

116. Nous avons démontré dans l’article 6 de la troisième Section que, dans tout système libre, les équations du mouvement des corps du système sont les mêmes, soit qu’on les rapporte au centre de gravité du système, ou à un point quelconque fixe hors du système. Ainsi, dans les formules de l’article 86, on pourra établir dans le centre de gravité de tous les corps l’origine de leurs coordonnées et, par les propriétés du centre de gravité, on aura les trois équations

lesquelles donnent tout de suite les valeurs des coordonnées de par celles des autres corps

Considérons en particulier le mouvement du corps autour du centre commun de gravité. Comme ses coordonnées sont indépendantes, on pourra, dans les formules de l’article cité, réduire les quantités et aux termes multipliés par qui sont les seuls qui contiennent les variables et diviser ensuite ces quantités par ainsi, dans l’équation générale, on pourra substituer et à la place de et en faisant

et

et l’on aura, pour chacune des trois coordonnées de l’orbite de autour du centre commun de gravité, une équation de la forme

étant une quelconque de ces coordonnées.

117. S’il n’y avait que deux corps et la valeur de se réduirait au seul terme et l’on aurait étant supposé fonction de

Pour avoir la valeur de relative à il faut différentier la variable

en ne faisant varier que et ensuite y substituer pour leurs valeurs

On aura ainsi

Or, par les mêmes substitutions, on a

Donc, faisant

rayon vecteur de l’orbite de on aura

et, par conséquent,

donc aussi de sorte que la valeur de deviendra étant maintenant une fonction de semblable à la fonction supposée de

Dans le cas de la nature, on a donc la force dirigée vers le centre commun de gravité, sera représentée par ce qui est connu.

118. Considérons maintenant le cas où le système est composé de plus de deux corps, et supposons, pour simplifier la question, que la masse soit beaucoup plus grande que chacune des autres masses ce qui est le cas des planètes à l’égard du Soleil ; les quantités deviendronttrès petites vis-à-vis des quantités dans le rapport des masses à la masse en vertu des équations données dans l’article précédent ; et l’on pourra, dans le développement, s’en tenir aux premières puissances de du moins tant qu’on ne voudra pas avoir égard aux carrés des masses.

Comme est supposé fonction de sera aussi une fonction de que nous dénoterons par et la quantité sera exprimée par suivant la notation des fonctions dérivées. Or on a

et

donc

et différentiant par les quantités demeurant constantes,

où il faudra mettre pour leurs valeurs

119. Supposons que la force d’attraction soit comme la puissance de la distance on aura

et la fonction sera une fonction homogène de du degré de sorte qu’on aura, par la propriété de ces fonctions,

donc, différentiant par et observant que

on aura

Donc, si l’on fait, pour abréger,

on aura, après les substitutions,

et telle sera la valeur de dans la différence (art. 116) ; de sorte qu’on aura, pour le premier terme de

Dans le système du monde, l’attraction des planètes est en raison inverse des carrés des distances ; on a ainsi

et l’on trouve

 

en faisant

Donc, si l’on fait ces substitutions dans la valeur de (art. 116), et qu’on y suppose aussi on aura, dans le cas de la nature, pour le mouvement du corps autour du centre commun de gravité,

Le premier terme de s’il était seul, donnerait, comme on l’a vu dans le Chapitre I, une orbite elliptique dans laquelle et comme les autres termes sont fort petits par rapport à celui-ci, étant multipliés par les masses supposées très petites vis-à-vis de on peut les regarder comme dus à des forces perturbatrices dont l’effet est de faire varier les éléments de l’orbite elliptique. Ainsi en faisant, comme dans l’article 90,

on pourra déterminer, par les formules données dans l’article 74, les variations de ces éléments.

120. Si l’on compare cette valeur de qu’on vient de trouver pour le mouvement du corps autour du centre de gravité du système avec celle de la quantité pour le mouvement du même corps autour du corps on voit qu’elles sont semblables ; les rayons vecteurs des orbites étant représentés dans celle-ci par et dans la précédente par et les quantités étant les mêmes dans l’une et dans l’autre, puisqu’elles représentent les distances rectilignes du corps aux autres corps il n’y a de différence qu’en ce qu’en changeant en il faut échanger entre elles les quantités ainsi que les quantités et ainsi des autres. Or, si l’on ne cherche que les variations séculaires des éléments, comme nous l’avons fait pour l’orbite de autour de il est facile de voir qu’on aura la même valeur de pour les deux orbites et, par conséquent, les mêmes formules pour ces variations, ce qui est assez remarquable.

Au reste, dans la valeur de trouvée ci-dessus, on pourrait effacer les termes parce qu’étant de la même forme que le premier terme de la quantité ils peuvent se joindre à ce terme, lequel deviendrait ainsi

de sorte que le corps décrirait autour du centre de gravité une orbite comme s’il y avait dans ce centre une masse égale à

et les forces perturbatrices de cette orbite seraient par conséquent, toutes choses d’ailleurs égales, moindres que celles de l’orbite du même corps autour du corps le plus gros


Séparateur

  1. Lorsque les formules donnent pour une valeur négative, l’équation apprend que est plus grand que l’unité, car est positif et égal à c’est pour cette raison que la trajectoire est une hyperbole. (J. Bertrand.)
  2. Voir les Mémoires de Berlin, années 1768-69 ; la Théorie des fonctions, chap. XVI, Ire Partie, et le Traité de la Résolution des équations, note 11. (Note de Lagrange.)
  3. Voir dans les Mémoires de l’Académie de Berlin de 1776 plusieurs applications de cette méthode*. (Note de Lagrange.)

    * Le Mémoires auquel renvoie Lagrange a pour titre Solution de quelques problèmes d’Astronomie par le moyen des séries et se trouve inséré au tome IV des Œuvres de Lagrange, p. 275 G. D.

    .
  4. Dans les Mémoires de l’Académie des Sciences pour 1823, Laplace a donné la condition nécessaire pour la convergence des séries précédentes. La même question a été traitée depuis par M. Cauchy dans les Comptes rendus de l’Académie des Sciences, et dans les Exercices d’Analyse et de Physique mathématique de 1841. L’analyse de M. Cauchy a été enfin développée et complétée par M. Puiseux dans le Journal de Mathématiques de M. Liouville, t. XIV ; 1849. Voir une Note à la fin du Volume. Legendre a donné d’ailleurs, dans les Exercices de Calcul intégral (Ve Partie, no 116) une série beaucoup plus convergente pour exprimer en fonction de Cette série s’appliquerait même à la planète Pallas*. (J. Bertrand.)

    * La série de Legendre est la suivante :

    G. D.

    .
  5. La formule relative au temps nécessaire pour parcourir un arc de parabole a été souvent attribuée à Lambert, qui, en effet, y est parvenu en 1761 sans avoir eu connaissance du Mémoire d’Euler où elle se trouve démontrée, et qui date cependant de 1744. Lagrange lui-même a partagé longtemps l’erreur dont nous parlons ; car, dans son premier Mémoire Sur le probléme de la détermination des orbites des comètes d’après trois observations (Œuvres de Lagrange, t. IV, p. 439), il dit, à propos de son théorème : « M. Lambert est parvenu à un des théorèmes les plus élégants et les plus utiles qui aient été trouvés sur ce sujet, et qui a en même temps l’avantage de s’appliquer aux orbites elliptiques. » Cette phrase a été imprimée en 1780, c’est-à-dire trois années avant la mort d’Euler, qui n’a jamais réclamé son droit de priorité. Dans les Mémoires de Berlin pour 1771, Lambert cite le même théorème et s’en attribue la découverte. (J. Bertrand.)
  6. Œuvres de Lagrange, t. IV, p. 439.
  7. C’est-à-dire les propriétés du mouvement parabolique développées plus haut. (J. Bertrand.)
  8. Euler avait donné, en 1744, une solution approchée du problème des comètes mais sa méthode exigeait plusieurs fausses hypothèses et l’emploi d’une quatrième observation. Quant à la solution de Lambert, on doit observer que celle dont Lagrange donne ici une idée succincte est la seconde qui ait été proposée par ce géomètre, qui a même négligé d’en développer les calculs. Le problème avait été résolu par lui d’une manière très différente dans un Ouvrage spécial, Insigniores orbitæ cometarum proprietates, qui date de 1761. (J. Bertrand.)
  9. C’est-à-dire, il n’est pas nécessaire que la vitesse excède fois celle d’un boulet de canon. (J. Bertrand.)
  10. C’est-à-dire la formule qui exprime la vitesse et celles qui font connaître la direction du mouvement. (J. Bertrand.)
  11. L’introduction de la fonction dans un raisonnement relatif au cas où cette fonction n’existe pas donne à ce paragraphe une apparence d’obscurité qu’un peu d’attention fait néanmoins disparaître. (J. Bertrand.)
  12. Voir les Mémoires de la première Classe de l’Institut pour 1808*. (Note de Lagrange.)

    * Œuvres de Lagrange, t. VI, p. 769.    G. D.

    .
  13. Il n’y a dans la question qu’une variable indépendante qui est le temps. Toute expression différentielle peut donc être intégrée, et l’observationde Lagrange semble, par conséquent, inutile. J’ajouterai qu’elle est de nature à embarrasser le lecteur, qui doit voir plus loin (art. 73)) adopter cette lettre comme l’une des constantes variables du problème. (J. Bertrand.)
  14. Il faut avouer que l’application ultérieure des formules de l’article 64 exigerait quelques explications, car ces formules supposaient exprimé en fonction des constantes du mouvement elliptique, et la lettre n’en est pas une. Cette lettre n’a pas même de sens précis, puisqu’elle n’a été définie que par sa différentielle. Voir, à ce sujet, des observations très fondées de M. Binet (Journal de l’École Polytechnique, XXVIIIe Cahier, t. XVII, p. 76). (J. Bertrand.)
  15. Voir le quatrième Volume des anciens Mémoires de Turin*.
    (Note de Lagrange.)

    * Œuvres de Lagrange, t. II, p. 253.    G. D.

  16. Dans sa Thèse présentée à la Faculté des Sciences de Paris, M. Serret a signalé la démonstration précédente comme peu rigoureuse ; il semble incontestable, en effet, que quelques développements sont nécessaires pour la rendre entièrement satisfaisante. (Voir une Note à la fin du Volume.) (J. Bertrand.)
  17. M. Ossian Bonnet a donné (Journal de M. Liouville, t. IX, p. 105) une raison très simple de ce fait qui est susceptible d’une généralisation intéressante. (Voir une Note à la fin du Volume.)(J. Bertrand.)
  18. Voir la note de la page 30.
  19. Œuvres de Lagrange, t. IV, p. 559. Le Mémoire auquel Lagrange fait allusion est intitulé Sur une manière particulière d’exprimer le temps dans les sections coniques décrites par des forces tendantes au foyer et réciproquement proportionnelles aux carrés des distances. Lagrange y indique différentes manières de démontrer le théorème d’Euler et de Lambert. G. D.
  20. Œuvres de Lagrange, t. II, p. 67. Le Mémoire est intitulé Recherches sur le mouvement d’un corps qui est attiré par deux centres fixes. G. D.
  21. Voir la Note relative à cet angle page 94. On peut remarquer, en outre, que l’équation est subordonnée à l’hypothèse que l’on puisse négliger l’inclinaison des deux orbites l’une sur l’autre. Si l’on n’avait pas supposé les formules seraient tout autres, et l’angle ne s’y introduirait pas. (J. Bertrand.)
  22. Voir à ce sujet un Mémoire de Laplace, Mémoires de l’Académie des Sciences de Paris pour 1784, le no 57 du second Livre de la Mécanique céleste, et le Chapitre II du livre XV (cinquième Volume) où Laplace reproche à Lagrange d’avoir énoncé sous une forme dubitative un théorème démontré par lui depuis longtemps. (J. Bertrand.)
  23. Œuvres de Lagrange, t. V, p. 125. G. D.
  24. Œuvres de Lagrange, t. V, p. 211. G. D.
  25. Ce supplément n’est pas de Lagrange. Voir à ce sujet la Note du tome V des Œuvres de Lagrange, p. 489. G. D.
  26. Qui sont des termes périodiques. (J. Bertrand.)