Théorie du mouvement des corps célestes/L1S3

Carl Friedrich Gauss

Théorie du mouvement des corps célestes

Traduction par Edmond Dubois.
1864 (p. 105-164).

TROISIÈME SECTION.

RELATIONS ENTRE PLUSIEURS POSITIONS DANS L’ORBITE.

78

La comparaison de deux ou plusieurs positions d’un corps céleste, soit dans son orbite, soit dans l’espace, fournit une telle abondance de propositions élégantes, qu’elles rempliraient facilement le volume entier. Notre but ne tend pas assurément à épuiser ce sujet fécond, mais principalement à en retirer des ressources importantes pour la solution du grand problème de la détermination des orbites inconnues, d’après les observations ; c’est pourquoi, en négligeant les questions qui seraient trop étrangères à notre but, nous développerons le plus soigneusement toutes celles qui, de quelque manière, peuvent y conduire. Nous faisons précéder ces recherches de quelques propositions trigonométriques auxquelles il faudra très-souvent avoir recours, puisqu’elles sont plus habituellement employées.

I. En désignant par $\mathrm {A} ,\,\mathrm {B} ,\,\mathrm {C}$ des angles quelconques, on a

${\begin{aligned}\sin \mathrm {A} \sin(\mathrm {C} -\mathrm {B} )+\,\sin \mathrm {B} \sin(\mathrm {A} -\mathrm {C} )+\sin \mathrm {C} \sin(\mathrm {B} -\mathrm {A} )&=0,\\\cos \mathrm {A} \sin(\mathrm {C} -\mathrm {B} )+\cos \mathrm {B} \sin(\mathrm {A} -\mathrm {C} )+\cos \mathrm {C} \sin(\mathrm {B} -\mathrm {A} )&=0.\end{aligned}}$

II. Si deux quantités $p,\,\mathrm {P}$ doivent être déterminées par des équations telles que

${\begin{aligned}p\sin(\mathrm {A} -\mathrm {P} )&=a,\\p\sin(\mathrm {B} -\mathrm {P} )&=b,\end{aligned}}$

ceci s’obtiendra généralement par le secours des formules

${\begin{aligned}p\sin(\mathrm {B} -\mathrm {A} )\sin(\mathrm {H} -\mathrm {P} )&=b\sin(\mathrm {H} -\mathrm {A} )-a\sin(\mathrm {H} -\mathrm {B} ),\\p\sin(\mathrm {B} -\mathrm {A} )\cos(\mathrm {H} -\mathrm {P} )&=b\cos(\mathrm {H} -\mathrm {A} )-a\cos(\mathrm {H} -\mathrm {B} ),\end{aligned}}$

dans lesquelles $\mathrm {H}$ est un angle arbitraire. De là se déduiront (art. 14, II) l’angle $(\mathrm {H} -\mathrm {P} ),$ et $p\sin(\mathrm {B} -\mathrm {A} )\,;$ et de là $\mathrm {P}$ et $p$. La plupart du temps il est ajouté la condition que $p$ doit être une quantité positive, d’où se trouve écartée l’ambiguïté qui résulte de la détermination de l’angle $\mathrm {H} -\mathrm {P}$ par sa tangente ; mais en l’absence de cette condition, l’incertitude sera levée arbitrairement. Pour la plus grande facilité du calcul, il sera convenable de prendre l’angle auxiliaire $\mathrm {H} =\mathrm {A}$ ou $=\mathrm {B}$ ou $={\frac {1}{2}}(\mathrm {A} +\mathrm {B} ).$ Dans le premier cas, les équations, pour la détermination de $\mathrm {P}$ et $p,$ seront

${\begin{aligned}p\sin(\mathrm {A} -\mathrm {P} )&=a,\\p\cos(\mathrm {A} -\mathrm {P} )&={\frac {b-a\cos(\mathrm {B} -\mathrm {A} )}{\sin(\mathrm {B} -\mathrm {A} )}}.\end{aligned}}$

Dans le second cas, les équations seront entièrement analogues ; mais dans le troisième on aura

${\begin{aligned}p\sin \left({\frac {1}{2}}\mathrm {A} +{\frac {1}{2}}\mathrm {B} -\mathrm {P} \right)&={\frac {b+a}{2\cos {\dfrac {1}{2}}(\mathrm {B} -\mathrm {A} )}},\\p\cos \left({\frac {1}{2}}\mathrm {A} +{\frac {1}{2}}\mathrm {B} -\mathrm {P} \right)&={\frac {b-a}{2\sin {\dfrac {1}{2}}(\mathrm {B} -\mathrm {A} )}}.\end{aligned}}$

Et alors, si l’on introduit l’angle auxiliaire $\zeta$, dont la tangente $={\frac {a}{b}}$, on trouvera $\mathrm {P}$ par la formule

$\operatorname {tang} \left({\frac {1}{2}}\mathrm {A} +{\frac {1}{2}}\mathrm {B} -\mathrm {P} \right)=\operatorname {tang} (45^{\circ }\!+\zeta )\operatorname {tang} {\frac {1}{2}}(\mathrm {B} -\mathrm {A} )$,

et ensuite $p$ par l’une des formules précédentes, dans lesquelles

${\begin{aligned}{\frac {1}{2}}(b+a)&=\sin(45^{\circ }\!+\zeta ){\sqrt {\frac {ab}{\sin 2\zeta }}}\\&={\frac {a\sin(45^{\circ }\!+\zeta )}{\sin \zeta {\sqrt {2}}}}={\frac {b\sin(45^{\circ }\!+\zeta )}{\cos \zeta {\sqrt {2}}}},\\[1ex]{\frac {1}{2}}(b-a)&=\cos(45^{\circ }\!+\zeta ){\sqrt {\frac {ab}{\sin 2\zeta }}}\\&={\frac {a\cos(45^{\circ }\!+\zeta )}{\sin \zeta {\sqrt {2}}}}={\frac {b\cos(45^{\circ }\!+\zeta )}{\cos \zeta {\sqrt {2}}}}.\end{aligned}}$

III. Si $p$ et $\mathrm {P}$ doivent être déterminées d’après les équations

${\begin{aligned}p\cos(\mathrm {A} -\mathrm {P} )&=a,\\p\cos(\mathrm {B} -\mathrm {P} )&=b,\end{aligned}}$

tout ce qui est exposé dans II peut immédiatement s’appliquer, pourvu qu’à la place de $\mathrm {A}$ et $\mathrm {B}$ on mette partout $90^{\circ }\!+\mathrm {A} ,$ $90^{\circ }\!+\mathrm {B} \,;$ mais pour que l’application en soit plus facile, nous ne craignons pas d’ajouter les formules développées. Les formules générales seront

${\begin{alignedat}{5}p\sin(\mathrm {B} -\mathrm {A} )\sin(\mathrm {H} -\mathrm {P} )&=-&&b\cos(\mathrm {H} -\mathrm {A} )&{}+{}&a\cos(\mathrm {H} -\mathrm {B} ),\\p\sin(\mathrm {B} -\mathrm {A} )\cos(\mathrm {H} -\mathrm {P} )&=&&b\sin(\mathrm {H} -\mathrm {A} )&{}-{}&a\sin(\mathrm {H} -\mathrm {B} ),\end{alignedat}}$

de manière que pour $\mathrm {H} =\mathrm {A}$ elles se changent en

${\begin{aligned}p\sin(\mathrm {A} -\mathrm {P} )&={\frac {a\cos(\mathrm {B} -\mathrm {A} )-b}{\sin(\mathrm {B} -\mathrm {A} )}},\\p\cos(\mathrm {A} -\mathrm {P} )&=a.\end{aligned}}$

Pour $\mathrm {H} =\mathrm {B}$ elles prennent une forme semblable ; mais pour $\mathrm {H} ={\frac {1}{2}}(\mathrm {A} +\mathrm {B} )$ elles deviennent

${\begin{aligned}p\sin \left({\frac {1}{2}}\mathrm {A} +{\frac {1}{2}}\mathrm {B} -\mathrm {P} \right)&={\frac {a-b}{2\sin {\dfrac {1}{2}}(\mathrm {B} -\mathrm {A} )}},\\p\cos \left({\frac {1}{2}}\mathrm {A} +{\frac {1}{2}}\mathrm {B} -\mathrm {P} \right)&={\frac {a+b}{2\cos {\dfrac {1}{2}}(\mathrm {B} -\mathrm {A} )}},\end{aligned}}$

de sorte que par l’introduction de l’angle auxiliaire $\zeta$, dont la tangente $={\frac {a}{b}},$ il vient

$\operatorname {tang} \left({\frac {1}{2}}\mathrm {A} +{\frac {1}{2}}\mathrm {B} -\mathrm {P} \right)=\operatorname {tang} (\zeta -45^{\circ }\!)\operatorname {cotang} {\frac {1}{2}}(\mathrm {B} -\mathrm {A} )$.

Enfin, si nous désirons déterminer $p$ immédiatement au moyen de $a$ et de $b,$ sans faire le calcul préalable de l’angle $\mathrm {P}$, nous avons la formule

$p\sin(\mathrm {B} -\mathrm {A} )={\sqrt {\left[a^{2}+b^{2}-2ab\cos(\mathrm {B} -\mathrm {A} )\right]}}$,

aussi bien dans le problème actuel que dans II.

79

Pour la détermination complète d’une section conique dans son plan, trois quantités sont demandées : la position du périhélie, l’excentricité et le demi-paramètre. Si ces quantités doivent être déterminées d’après les quantités données qui en dépendent, il faut qu’il y ait assez de données pour pouvoir former trois équations indépendantes les unes des autres. Tout rayon vecteur donné en grandeur et en position fournit une équation ; c’est pourquoi trois rayons donnés en grandeur et en position sont nécessaires pour la détermination d’une orbite. Mais si l’on en a deux seulement, un élément même doit être déjà donné, ou au moins quelque autre quantité à l’aide de laquelle il soit permis d’établir une troisième équation. De là surgit une variété de problèmes que nous traiterons maintenant successivement.

Soient $r$, $r'$ deux rayons vecteurs qui font avec une droite arbitraire menée par le Soleil, dans le plan de l’orbite, les angles $\mathrm {N}$, $\mathrm {N} '$ selon la direction du mouvement ; soit ensuite $\Pi$ l’angle que fait, avec la même droite, le rayon vecteur mené au périhélie, de telle sorte que les anomalies vraies $\mathrm {N} -\Pi$, $\mathrm {N} '-\Pi$ répondent aux rayons vecteurs $r$, $r'$ ; soient enfin $e$ l’excentricité, $p$ le demi-paramètre. On obtient alors les équations

${\begin{aligned}{\frac {p}{r}}&=1+e\cos(\mathrm {N} -\Pi ),\\{\frac {p}{r'}}&=1+e\cos(\mathrm {N} '-\Pi ),\end{aligned}}$

desquelles, si l’une des quantités $p$, $e$, $\Pi$ est en outre donnée, on pourra déterminer les deux autres.

Supposons d’abord que le demi-paramètre $p$ soit donné, et il est évident que la détermination des quantités $e,$ $\Pi ,$ au moyen des équations

${\begin{aligned}e\cos(\mathrm {N} -\Pi )&={\frac {p}{r}}-1,\\e\cos(\mathrm {N} '-\Pi )&={\frac {p}{r'}}-1,\end{aligned}}$

peut s’effectuer d’après la règle du lemme III de l’article précédent. Nous avons donc

${\begin{aligned}\operatorname {tang} (\mathrm {N} -\Pi )=\operatorname {cotang} (\mathrm {N} '-\mathrm {N} )-{\frac {r(p-r')}{r'(p-r)\sin(\mathrm {N} '-\mathrm {N} )}},\\\operatorname {tang} \left({\frac {1}{2}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} '-\Pi \right)={\frac {(r'-r)\operatorname {cotang} {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} )}{r'+r-{\dfrac {2rr'}{p}}}}.\end{aligned}}$

80

Si l’angle $\Pi$ est donné, $p$ et $e$ seront déterminés au moyen des équations

${\begin{aligned}p&={\frac {rr'\left[\cos(\mathrm {N} -\Pi )-\cos(\mathrm {N} '-\Pi )\right]}{r\cos(\mathrm {N} -\Pi )-r'\cos(\mathrm {N} '-\Pi )}},\\e&={\frac {r'-r}{r\cos(\mathrm {N} -\Pi )-r'\cos(\mathrm {N} '-\Pi )}},\end{aligned}}$

On peut réduire le dénominateur commun dans ces formules, à la forme $a\cos(\mathrm {A} -\Pi ),$ de telle sorte que $a$ et $\mathrm {A}$ soient indépendantes de $\Pi .$ En désignant en effet par $\mathrm {H}$ un angle arbitraire, nous avons

${\begin{alignedat}{5}r\cos(\mathrm {N} -\Pi )-&{}r'\cos(\mathrm {N} '-\Pi )&=&\\{}&{\big [}r\cos(\mathrm {N} -\mathrm {H} )&{}-{}&r'\cos(\mathrm {N} '-\mathrm {H} ){\big ]}&&\cos(\mathrm {H} -\Pi )\\-{}&{\big [}r\sin(\mathrm {N} -\mathrm {H} )&{}-{}&r'\sin(\mathrm {N} '-\mathrm {H} ){\big ]}&&\sin(\mathrm {H} -\Pi ),\end{alignedat}}$

et par suite,

$=a\cos(\mathrm {A} -\Pi )$,

si $a$ et $\mathrm {A}$ sont déterminées par les équations

${\begin{aligned}r\cos(\mathrm {N} -\mathrm {H} )-r'\cos(\mathrm {N} '-\mathrm {H} )&=a\cos(\mathrm {A} -\mathrm {H} ),\\r\sin(\mathrm {N} -\mathrm {H} )-r'\sin(\mathrm {N} '-\mathrm {H} )&=a\sin(\mathrm {A} -\mathrm {H} ),\end{aligned}}$

De cette manière il vient.

${\begin{aligned}p&={\frac {2rr'\sin {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} )\sin \left({\dfrac {1}{2}}\mathrm {N} +{\dfrac {1}{2}}\mathrm {N} '-\Pi \right)}{a\cos(\mathrm {A} -\mathrm {B} )}},\\e&={\frac {r'-r}{a\cos(\mathrm {A} -\mathrm {B} )}}.\end{aligned}}$

Ces formules sont principalement commodes toutes les fois que $p$ et $e$ sont calculés pour plusieurs valeurs de $\Pi \,;$ $r,$ $r',$ $\mathrm {N} ,$ $\mathrm {N} '$ ne changeant pas. — Comme pour le calcul des quantités auxiliaires $a,$ $\mathrm {A}$ il est permis de prendre l’angle $\mathrm {H}$ arbitrairement, on pourra alors poser $\mathrm {H} ={\frac {1}{2}}(\mathrm {N} +\mathrm {N} '),$ d’après quoi les formules se changent en celles-ci :

${\begin{aligned}(r'-r)\cos {\frac {1}{2}}(\mathrm {N} '-\mathrm {N} )&=-a\cos \left(\mathrm {A} -{\frac {1}{2}}\mathrm {N} -{\frac {1}{2}}\mathrm {N} '\right),\\(r+r')\sin {\frac {1}{2}}(\mathrm {N} '-\mathrm {N} )&=-a\sin \left(\mathrm {A} -{\frac {1}{2}}\mathrm {N} -{\frac {1}{2}}\mathrm {N} '\right).\end{aligned}}$

L’angle $\mathrm {A}$ étant donc déterminé par l’équation

$\operatorname {tang} \left(\mathrm {A} -{\frac {1}{2}}\mathrm {N} -{\frac {1}{2}}\mathrm {N} '\right)={\frac {r'+r}{r'-r}}\operatorname {tang} {\frac {1}{2}}(\mathrm {N} '-\mathrm {N} )$,

on a aussitôt

$e=-{\frac {\cos \left(\mathrm {A} -{\dfrac {1}{2}}\mathrm {N} -{\dfrac {1}{2}}\mathrm {N} '\right)}{\cos {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} )\cos(\mathrm {A} -\Pi )}}$.

Le calcul du logarithme de l’expression ${\frac {r'+r}{r'-r}}$ pourra se simplifier par une méthode déjà souvent expliquée.

81

Si l’excentricité est donnée, on trouve l’angle $\Pi$ par l’équation

$\cos(\mathrm {A} -\Pi )=-{\frac {\cos \left(\mathrm {A} -{\dfrac {1}{2}}\mathrm {N} -{\dfrac {1}{2}}\mathrm {N} '\right)}{e\cos {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} )}}$,

après que l’angle auxiliaire $\mathrm {A}$ a été déterminé par l’équation

$\operatorname {tang} \left(\mathrm {A} -{\frac {1}{2}}\mathrm {N} -{\frac {1}{2}}\mathrm {N} '\right)={\frac {r'+r}{r'-r}}\operatorname {tang} {\frac {1}{2}}(\mathrm {N} '-\mathrm {N} )$.

L’ambiguïté qui existe dans la détermination de l’angle $\mathrm {A} -\Pi$, par son cosinus, résulte de la nature du problème, de manière qu’on peut satisfaire à la question par deux solutions différentes ; il faudra décider, d’autre part, laquelle devra être adoptée et laquelle rejetée ; dans ce but, la valeur au moins approchée de $\Pi$ devra déjà être connue. — Une fois $\Pi$ trouvé, on calculera $p$ par les formules

$p=r\left[1+e\cos(\mathrm {N} -\Pi )\right]=r'\left[1+e\cos(\mathrm {N} '-\Pi )\right]$,

ou par celle-ci :

$p={\frac {2rr'e\sin {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} )\sin \left({\dfrac {1}{2}}\mathrm {N} '+{\dfrac {1}{2}}\mathrm {N} -\Pi \right)}{r'-r}}$

82

Supposons enfin, que nous connaissions les trois rayons vecteurs $r$, $r',$ $r''$ qui font, avec la droite menée arbitrairement par le Soleil dans le plan de l’orbite, les angles $\mathrm {N} ,$ $\mathrm {N} ',$ $\mathrm {N} ''.$ Nous aurons alors, en conservant les mêmes notations, les équations (I) :

${\begin{aligned}{\frac {p}{r}}&=1+e\cos(\mathrm {N} -\Pi ),\\{\frac {p}{r'}}&=1+e\cos(\mathrm {N} '-\Pi ),\\{\frac {p}{r''}}&=1+e\cos(\mathrm {N} ''-\Pi ),\end{aligned}}$

qui permettront d’obtenir $p$, $\Pi$, $e$ de plusieurs manières. Si l’on veut calculer la quantité $p$ avant les autres, on doit multiplier les trois équations (I) respectivement par $\sin(\mathrm {N} ''-\mathrm {N} '),$ $-\sin(\mathrm {N} ''-\mathrm {N} ),$ $\sin(\mathrm {N} '-\mathrm {N} ),$ et l’on aura, d’après le lemme 1, art. 78, en ajoutant les produits,

$p={\frac {\sin(\mathrm {N} ''-\mathrm {N} ')-\sin(\mathrm {N} ''-\mathrm {N} )+\sin(\mathrm {N} '-\mathrm {N} )}{{\dfrac {1}{r}}\sin(\mathrm {N} ''-\mathrm {N} ')-{\dfrac {1}{r'}}\sin(\mathrm {N} ''-\mathrm {N} )+{\dfrac {1}{r''}}\sin(\mathrm {N} '-\mathrm {N} )}}$.

Cette expression mérite d’être considérée plus attentivement. Le numérateur devient évidemment,

${\begin{aligned}2\sin {\frac {1}{2}}(\mathrm {N} ''&-\mathrm {N} ')\cos {\frac {1}{2}}(\mathrm {N} ''-\mathrm {N} ')\\&\quad -2\sin {\frac {1}{2}}(\mathrm {N} ''-\mathrm {N} ')\cos \left({\frac {1}{2}}\mathrm {N} ''+{\frac {1}{2}}\mathrm {N} '-\mathrm {N} \right)\\&=4\sin {\frac {1}{2}}(\mathrm {N} ''-\mathrm {N} ')\sin {\frac {1}{2}}(\mathrm {N} ''-\mathrm {N} )\sin {\frac {1}{2}}(\mathrm {N} '-\mathrm {N} ).\end{aligned}}$

En posant, ensuite,

${\begin{alignedat}{5}r'r''&\sin(\mathrm {N} ''&{}-{}&\mathrm {N} ')&{}={}&n,\\rr''&\sin(\mathrm {N} ''&{}-{}&\mathrm {N} )&{}={}&n',\\rr'&\sin(\mathrm {N} '&{}-{}&\mathrm {N} )&{}={}&n'',\end{alignedat}}$

il est évident que ${\frac {1}{2}}n$, ${\frac {1}{2}}n'$, ${\frac {1}{2}}n''$ sont les aires des triangles compris entre le second rayon vecteur et le troisième, entre le premier et le troisième, entre le premier et le second. De là on apercevra facilement, dans la nouvelle formule

$p={\frac {4\sin {\dfrac {1}{2}}(\mathrm {N} ''-\mathrm {N} ')\sin {\dfrac {1}{2}}(\mathrm {N} ''-\mathrm {N} )\sin {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} ).rr'r''}{n-n'+n''}}$

que le dénominateur est le double de l’aire du triangle compris entre les extrémités des trois rayons vecteurs, c’est-à-dire entre les trois positions du corps céleste dans l’espace. Toutes les fois que ces positions sont peu écartées les unes des autres, cette aire est toujours une quantité très-petite et certainement du troisième ordre si $\mathrm {N} '-\mathrm {N} ,$ $\mathrm {N} ''-\mathrm {N} '$ sont considérées comme de petites quantités du premier ordre. On conclut de là en même temps, que si une ou plusieurs des quantités $r,$ $r',$ $r'',$ $\mathrm {N} ,$ $\mathrm {N} ',$ $\mathrm {N} ''$ sont affectées d’erreurs même légères, il pourra en résulter, sur la détermination de $p$, une erreur très-grande ; c’est pourquoi cette méthode-ci n’admet jamais une grande précision, à moins que les trois lieux héliocentriques ne soient distants l’un de l’autre d’un intervalle considérable.

Enfin, du moment que le demi-paramètre $p$ sera trouvé, $e$ et $\Pi$ se détermineront par la combinaison de deux quelconques des équations (I), d’après la méthode de l’art. 79.

83

Si nous aimons mieux commencer la résolution du même problème par le calcul de l’angle $\Pi$, nous emploierons la méthode suivante. Nous retranchons, dans les équations (I) la troisième de la seconde, la troisième de la première, la seconde de la première ; nous obtenons ainsi les trois nouvelles équations (II) :

${\begin{aligned}{\frac {{\dfrac {1}{r'}}-{\dfrac {1}{r''}}}{2\sin {\dfrac {1}{2}}(\mathrm {N} ''-\mathrm {N} ')}}={\frac {e}{p}}\sin \left({\frac {1}{2}}\mathrm {N} '+{\frac {1}{2}}\mathrm {N} ''-\Pi \right),\\{\frac {{\dfrac {1}{r}}-{\dfrac {1}{r''}}}{2\sin {\dfrac {1}{2}}(\mathrm {N} ''-\mathrm {N} )}}={\frac {e}{p}}\sin \left({\frac {1}{2}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} ''-\Pi \right),\\{\frac {{\dfrac {1}{r}}-{\dfrac {1}{r'}}}{2\sin {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} )}}={\frac {e}{p}}\sin \left({\frac {1}{2}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} '-\Pi \right).\end{aligned}}$

Deux quelconques de ces équations donneront, d’après le lemme II, art. 78, $\Pi$ et ${\frac {e}{p}}$, d’où l’on obtiendra $e$ et $p$ par l’une ou l’autre des équations (I). Si nous adoptons la troisième solution enseignée dans l’art. 78, II, la combinaison de la première équation avec la troisième produit l’algorithme suivant : l’angle auxiliaire $\zeta$ sera déterminé par l’équation

$\operatorname {tang} \zeta ={\frac {{\dfrac {r'}{r}}-1}{1-{\dfrac {r'}{r''}}}}.{\frac {\sin {\dfrac {1}{2}}(\mathrm {N} ''-\mathrm {N} ')}{\sin {\dfrac {1}{2}}(\mathrm {N} '-\mathrm {N} )}}$,

et l’on aura

$\operatorname {tang} \left({\frac {1}{4}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} '+{\frac {1}{4}}\mathrm {N} ''-\Pi \right)=\operatorname {tang} (45^{\circ }\!+\zeta )\operatorname {tang} {\frac {1}{4}}(\mathrm {N} ''-\mathrm {N} ).$

En changeant le second lieu avec le premier ou avec le troisième, on obtiendra deux autres solutions entièrement analogues à celle-ci.

Comme en employant cette méthode les formules relatives à ${\frac {e}{p}}$ se résolvent moins promptement, il vaudra mieux déduire $e$ et $p$ des deux équations (I) par la méthode de l’art. 80. Enfin, l’ambiguïté dans la détermination de $\Pi$ par la tangente de l’angle ${\frac {1}{4}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} '+{\frac {1}{4}}\mathrm {N} ''-\Pi$ devra être ici écartée par la considération que $e$ est une quantité positive ; il est en effet évident, que $e$ acquière des valeurs opposées si l’on prend pour $\Pi$ des valeurs différant entre elles de 180°. Mais le signe de $p$ ne dépend pas de cette ambiguïté et sa valeur ne peut être prise négativement, à moins que les trois points donnés n’appartiennent à un arc d’hyperbole opposé au Soleil, cas contraire aux lois de la nature, et auquel nous n’avons pas ici égard.

Les quantités que l’on obtiendrait, d’après l’application de la première méthode de l’art. 78, II, après de pénibles substitutions, peuvent être, dans le cas actuel, déterminées plus commodément de la manière suivante : Que l’on multiplie la première des équations II par $\cos {\frac {1}{2}}(\mathrm {N} ''-\mathrm {N} '),$ la troisième par $\cos {\frac {1}{2}}(\mathrm {N} '-\mathrm {N} ),$ et que l’on retranche le dernier produit du premier. Alors, par l’application exacte

du lemme I de l’art. 78^[1], on obtiendra l’équation

${\begin{aligned}{\frac {1}{2}}\left({\frac {1}{r'}}\!-\!{\frac {1}{r''}}\right)\operatorname {cotang} {\frac {1}{2}}(\mathrm {N} ''\!-\!\mathrm {N} ')-{\frac {1}{2}}\left({\frac {1}{r}}\!-\!{\frac {1}{r'}}\right)\operatorname {cotang} {\frac {1}{2}}(\mathrm {N} '\!-\!\mathrm {N} )\\={\frac {e}{p}}\sin {\frac {1}{2}}(\mathrm {N} ''\!-\!\mathrm {N} )\cos \left({\frac {1}{2}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} ''-\Pi \right).\end{aligned}}$

En combinant cette équation avec la seconde des équations II, $\Pi$ et ${\frac {e}{p}}$ seront déterminés ; et alors $\Pi$, par la formule.

${\begin{aligned}&\operatorname {tang} \left({\frac {1}{2}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} ''-\Pi \right)\\[0ex]&\qquad ={\frac {{\dfrac {r'}{r}}-{\dfrac {r'}{r''}}}{\left(1\!-\!{\dfrac {r'}{r''}}\right)\operatorname {cotang} {\dfrac {1}{2}}(\mathrm {N} ''\!-\!\mathrm {N} )-\left({\dfrac {r'}{r}}\!-\!1\right)\operatorname {cotang} {\dfrac {1}{2}}(\mathrm {N} '\!-\!\mathrm {N} )}}.\end{aligned}}$

De là aussi, dérivent deux autres formules entièrement analogues en changeant le second lieu avec le premier ou avec le troisième.

84

Puisqu’on peut déterminer l’orbite entière, au moyen de deux rayons vecteurs donnés de grandeur et de position, avec un élément de l’orbite, on pourra aussi, avec ces données, déterminer le temps pendant lequel le corps céleste se meut d’un rayon vecteur à l’autre, si à la vérité nous négligeons la masse de l’astre ou si nous la regardons comme connue : nous nous en tiendrons à la première hypothèse, à laquelle la seconde est facilement réduite. De là, réciproquement, il est évident que deux rayons vecteurs donnés de grandeur et de position avec le temps pendant lequel le corps céleste décrit l’espace compris, déterminent l’orbite entière. Mais ce problème, qui doit être considéré comme très-important dans la théorie du mouvement des corps célestes, n’est pas si facilement résolu, puisque l’expression du temps en fonction des éléments est transcendante et de plus assez compliquée. Il est, pour cela, le plus digne d’être traité avec tout le soin possible ; c’est pourquoi, nous espérons qu’il ne sera pas désagréable au lecteur que, outre la solution enseignée plus loin, solution qui parait ne rien laisser à désirer, nous ayons considéré comme devant aussi être arrachée à l’oubli, celle dont nous nous sommes fréquemment servis, avant que la première ne se soit offerte à nous. Il est toujours favorable d’attaquer par plusieurs voies un problème plus difficile, et de ne pas mépriser la bonne quoiqu’on préfère la meilleure. Nous entamons la question par l’exposition de cette méthode plus ancienne.

85

Nous conserverons aux lettres $r$, $r'$, $\mathrm {N}$, $\mathrm {N} '$, $p$, $e$, $\Pi$ la même signification que précédemment. Nous désignerons par $\Delta$ la différence $\mathrm {N} '-\mathrm {N} ,$ et par $t,$ le temps pendant lequel le corps céleste se transporte de la première position à la seconde. Il est maintenant évident, que si l’on connaît la valeur approchée de l’une des quantités $p,$ $e,$ $\Pi ,$ on pourra aussi en déduire la valeur des deux autres, et ensuite, par les méthodes développées dans la première section, le temps correspondant au mouvement de l’astre pour aller de la première position à la seconde. Si ce temps se trouve égal à l’intervalle proposé $t,$ la valeur supposée de $p,$ $e$ ou $\Pi$ est exacte et l’orbite est déjà trouvée : si cela n’est pas, le calcul recommencé avec une autre valeur un peu différente de la première, montrera quelle variation, dans la valeur du temps, correspond à une petite variation dans la valeur de $p,$ $e,$ ou $\Pi \,;$ de là, par une simple interpolation, on déterminera la valeur corrigée. Si le calcul est de nouveau recommencé avec cette valeur, le temps qu’on en déduira s’accordera entièrement avec le proposé, ou au moins, en différera d’une très-petite quantité, de telle sorte certainement, qu’on pourra atteindre, par de nouvelles corrections, un accord aussi parfait que le permettent les tables logarithmiques et trigonométriques.

Le problème se réduit donc à ceci ; — que, pour le cas où l’orbite est entièrement inconnue, nous sachions déterminer une valeur au moins approchée de l’une quelconque des quantités $p,$ $e,$ $\Pi .$ Nous donnerons maintenant une méthode par laquelle la valeur de $p$ est obtenue avec tant de précision que, pour de petites valeurs de $\Delta ,$ elle n’exige aucune nouvelle correction ; et par suite, l’orbite entière est déterminée, par un premier calcul, avec toute la précision que permettent les tables vulgaires. Mais autrement, il ne faudra presque jamais recourir à cette méthode, si ce n’est pour des valeurs médiocres de $\Delta ,$ parce qu’on ne peut guère entreprendre la détermination d’une orbite entièrement inconnue, à cause de la complication trop embarrassée du problème, qu’avec des observations peu écartées les unes des autres, ou plutôt telles qu’elles répondent à un mouvement héliocentrique peu considérable.

86

En désignant par $\rho$ le rayon vecteur indéterminé ou variable qui répond à l’anomalie vraie $v-\Pi ,$ l’aire du secteur décrit par le corps céleste dans le temps $t$, sera ${\frac {1}{2}}\int \rho ^{2}dv$, cette intégrale étant prise depuis $v=\mathrm {N}$ jusqu’à $v=\mathrm {N} ',$ et par suite, en prenant $k$ avec sa signification de l’art. 6, on a $kt{\sqrt {\overset {}{p}}}=\int \rho ^{2}dv.$ Il est actuellement évident, d’après les formules développées par Cotes, que si $\varphi x$ exprime une fonction quelconque de $x,$ des valeurs de plus en plus approchées de l’intégrale $\int \varphi x.dx,$ prise depuis $x=u$ jusqu’à $x=u+\Delta$ s’obtiendront par les formules

${\begin{aligned}&{\frac {1}{2}}\Delta \left[\varphi u+\varphi (u\!+\!\Delta )\right],\\&{\frac {1}{6}}\Delta \left[\varphi u+4\varphi \left(u\!+\!{\frac {1}{2}}\Delta \right)+\varphi (u\!+\!\Delta )\right],\\&{\frac {1}{8}}\Delta \left[\varphi u+3\varphi \left(u\!+\!{\frac {1}{3}}\Delta \right)+3\varphi \left(u\!+\!{\frac {2}{3}}\Delta \right)+\varphi (u+\Delta )\right]+{\text{etc.}},\,{\text{etc.}}\end{aligned}}$

Il suffira, pour notre but, de s’arrêter aux deux premières formules.

Par la première, nous avons donc, dans notre problème,

$\int \rho ^{2}dv={\frac {1}{2}}\Delta (r^{2}+r'^{2})={\frac {\Delta rr'}{\cos 2\omega }}$,

si nous posons

${\frac {r'}{r}}=\operatorname {tang} (45^{\circ }\!+\omega )$.

C’est pourquoi, la première valeur approchée de ${\sqrt {\overset {}{p}}}$ sera ${\frac {\Delta rr'}{kt\cos \omega }}\,;$ que nous posons $=3\alpha$.

Par la seconde formule, nous avons plus exactement

$\int \rho ^{2}dv={\frac {1}{6}}\Delta (r^{2}+r'^{2}+4\mathrm {R} ^{2})$,

en désignant par $\mathrm {R}$ le rayon vecteur qui correspond à l’anomalie intermédiaire ${\frac {1}{2}}\mathrm {N} +{\frac {1}{2}}\mathrm {N} '-\Pi .$

En exprimant maintenant $p,$ en fonction de $r,$ $\mathrm {R} ,$ $r',$ $\mathrm {N} ,$ $\mathrm {N} +{\frac {1}{2}}\Delta ,$ $\mathrm {N} +\Delta ,$ d’après les formules données dans l’article 82, nous trouvons

$p={\frac {4\sin ^{2}{\dfrac {1}{4}}\Delta \sin {\dfrac {1}{2}}\Delta }{\left({\dfrac {1}{r}}+{\dfrac {1}{r'}}\right)\sin {\dfrac {1}{2}}\Delta -{\dfrac {1}{\mathrm {R} }}\sin \Delta }},$

et de là

${\begin{aligned}{\frac {\cos {\dfrac {1}{2}}\Delta }{\mathrm {R} }}&={\frac {1}{2}}\left({\frac {1}{r}}+{\frac {1}{r'}}\right)-{\frac {2\sin ^{2}{\dfrac {1}{4}}\Delta }{p}}\\&={\frac {\cos \omega }{\sqrt {rr'\cos 2\omega }}}-{\frac {2\sin ^{2}{\dfrac {1}{4}}\Delta }{p}}.\end{aligned}}$

En posant donc,

${\frac {2\sin ^{2}{\dfrac {1}{4}}\Delta {\sqrt {rr'\cos 2\omega }}}{\cos \omega }}=\delta ,$

il vient

$\mathrm {R} ={\frac {\cos {\dfrac {1}{2}}\Delta {\sqrt {rr'\cos 2\omega }}}{\cos \omega \left(1-{\dfrac {\delta }{p}}\right)}},$

d’où l’on obtient la seconde valeur approchée de ${\sqrt {\overset {}{p}}}$,

${\sqrt {\overset {}{p}}}=\alpha +{\frac {2\alpha \cos ^{2}{\dfrac {1}{2}}\Delta \cos ^{2}2\omega }{\cos ^{2}\omega \left(1-{\dfrac {\delta }{p}}\right)^{2}}}=\alpha +{\frac {\varepsilon }{\left(1-{\dfrac {\delta }{p}}\right)^{2}}}$,

si nous posons

$2\alpha \left({\frac {\cos {\dfrac {1}{2}}\Delta \cos 2\omega }{\cos \omega }}\right)^{2}=\varepsilon$.

C’est pourquoi, en écrivant $\pi$ à la place de ${\sqrt {\overset {}{p}}},$ on déterminera $\pi$ par l’équation

$(\pi -\alpha )\left(1-{\frac {\delta }{\pi ^{2}}}\right)^{2}=\varepsilon$,

qui convenablement développée s’élèverait au cinquième degré. Posons $\pi =q+\mu ,$ de telle sorte que $q$ soit une valeur approchée de $\pi ,$ et $\mu$ une très-petite quantité dont on peut négliger le carré et les puissances supérieures ; on trouve par cette substitution,

$(q-\alpha )\left(1-{\frac {\delta }{q^{2}}}\right)^{2}+\mu \left[\left(1-{\frac {\delta }{q^{2}}}\right)^{2}+{\frac {4\delta (q-\alpha )}{q^{3}}}\left(1-{\frac {\delta }{q^{2}}}\right)\right]=\varepsilon$

ou

$\mu ={\frac {\varepsilon q^{5}-(q^{2}-\alpha q)(q^{2}-\delta )^{2}}{(q^{2}-\delta )(q^{3}+3\delta q-4\alpha \delta )}},$

et par suite,

$\pi ={\frac {\varepsilon q^{5}-(q^{2}-\delta )(\alpha q^{2}+4\delta q-5\alpha \delta )q}{(q^{2}-\delta )(q^{3}+3\delta q-4\alpha \delta )}}.$

Nous avons maintenant, dans notre problème, la valeur approximative de $\pi ,$ à savoir $3\alpha ,$ qui étant substituée à la place de $q$ dans la formule précédente, donne la valeur corrigée

$\pi ={\frac {243\alpha ^{4}\varepsilon +3\alpha (9\alpha ^{2}-\delta )(9\alpha ^{2}+7\delta )}{(9\alpha ^{2}-\delta )(27\alpha ^{2}+5\delta )}}.$

En posant donc,

${\frac {\delta }{27\alpha ^{2}}}=\beta ,$  ${\frac {\varepsilon }{(1-3\beta )\alpha }}=\gamma ,$

la formule prend cette forme-ci :

$\pi ={\frac {\alpha (1+\gamma +21\beta )}{1+5\beta }},$

et toutes les opérations nécessaires à la solution du problème sont contenues dans les cinq formules :

I.

${\frac {r'}{r}}=\operatorname {tang} (45^{\circ }\!+\omega ).$

II.

${\frac {\Delta rr'}{3kt\cos 2\omega }}=\alpha .$

III.

${\frac {2\sin ^{2}{\dfrac {1}{4}}\Delta {\sqrt {rr'\cos 2\omega }}}{27\alpha ^{2}\cos \omega }}=\beta .$

IV.

${\frac {2\cos ^{2}{\dfrac {1}{2}}\Delta \cos ^{2}2\omega }{(1-3\beta )\cos ^{2}\omega }}=\gamma .$

V.

${\frac {\alpha (1+\gamma +21\beta )}{1+5\beta }}={\sqrt {\overset {}{p}}}.$

S’il plaît de diminuer quelque chose de la précision de ces formules, on pourra obtenir des expressions encore plus simples. En faisant, en effet, $\cos \omega$ et $\cos 2\omega =1,$ et en développant la valeur de ${\sqrt {\overset {}{p}}}$ en série suivant les puissances croissantes de $\Delta ,$ on obtient, en négligeant les quatrièmes puissances et celles plus élevées,

${\begin{aligned}{\sqrt {\overset {}{p}}}=\alpha \left(3-{\frac {1}{2}}\Delta ^{2}+{\frac {\Delta ^{2}{\sqrt {rr'}}}{18\alpha ^{2}}}\right).\end{aligned}}$

dans lesquelles $\Delta$ doit être exprimé en parties du rayon. C’est pourquoi, en faisant ${\frac {\Delta rr'}{kt}}={\sqrt {p'}},$ on a

VI.

$p=p'\left(1-{\frac {1}{3}}\Delta ^{2}+{\frac {\Delta ^{2}{\sqrt {rr'}}}{3p'}}\right).$

De la même manière, en développant ${\sqrt {\overset {}{p}}}$ en série, suivant les puissances croissantes de $\sin \Delta ,$ on obtient, en posant ${\frac {rr'\sin \Delta }{kt}}={\sqrt {p''}},$

VII.

${\sqrt {\overset {}{p}}}=\left(1+{\frac {\sin ^{2}\Delta {\sqrt {rr'}}}{6p''}}\right){\sqrt {p''}},$

ou

VIII.

$p=p''+{\frac {1}{3}}\sin ^{2}\Delta {\sqrt {rr'}}.$

Les formules VII et VIII s’accordent avec celles que l’illustre Euler a développées dans le « Theoria motus planetarum et cometarum, » et la formule VI avec celle qui a été employée dans les « Recherches et calculs sur la vraie orbite elliptique de la comète de 1769, » p. 80.

87

Les exemples suivants éclairciront la pratique des méthodes précédentes, et permettront en même temps d’estimer le degré de précision.

I. Soient $\log r=0,3307640,$ $\log r'=0.3222239,$ $\Delta =7^{\circ }34'53''\!73=27293''\!,73,$ $t=21,93391$ jours. On trouve ici $\omega =-33'47''\!,90,$ d’où le calcul ultérieur se fait de la manière suivante :

	$\log \Delta .........$	4,4360629		${\tfrac {1}{2}}$	$rr'\cos 2\omega .....$	0,3264519
	$\log rr'........$	0,6529879		$2$	$\log \sin {\tfrac {1}{4}}\Delta ....$	7,0389972
$\mathrm {c^{t}}$	$\log 3k\,........$	5,9728722			$\log {\tfrac {2}{27}}........$	8,8696662
$\mathrm {c^{t}}$	$\log t..........$	8,6588840		$\mathrm {c^{t}}$	$\log \alpha ^{2}........$	0,5582180
$\mathrm {c^{t}}$	$\log \cos 2\omega .....$	0,0000840		$\mathrm {c^{t}}$	$\log \cos \omega ......$	0,0000210
	$\log \alpha .........$	9,7208910			$\log \beta .........$	6,7933543
$\beta =...$ 0,0006213757
	$\log 2\,.........$	0,3010300
$2$	$\log \cos {\tfrac {1}{2}}\Delta ....$	9,9980976		$1+\gamma +21\beta =$ 3,0074471
$2$	$\log \cos 2\omega .....$	9,9998320			$\log \;..........$	0,4781980
$\mathrm {c^{t}}$	$\log(1-3\beta )...$	0,0008103			$\log \alpha .........$	9,7208910
$2\mathrm {c^{t}}$	$\log \cos 2\omega .....$	0,0000420		$\mathrm {c^{t}}$	$\log(1+5\beta )...$	9,9986528
	$\log \gamma ..........$	0,2998119			$\log {\sqrt {\overset {}{p}}}.......$	0,1977418
$\gamma ={}$ 1,9943982					$\log p.........$	0,3954836
$21\beta ={}$ 0,0130489

Cette valeur de $\log p$ diffère de la vraie valeur, à peine d’une unité du septième ordre : la formule VI, dans cet exemple, donne $\log p=0,3954822\,;$ la formule VII donne $3954780\,;$ enfin, la formule VIII, $0,3954754.$

II. Soient $\log r=0,4282792,$ $\log r'=0,4062033,$ $\Delta =62^{\circ }55'16''\!,64,$ $t=259,88477$ jours. On déduit de là $\omega =-1^{\circ }27'20''\!,14,$ $\log \alpha =9,7482348,$ $\beta =0,04535216,$ $\gamma =1,681127,$ $\log {\sqrt {\overset {}{p}}}=0,2198027,$ $\log p=0,4396054,$ qui est inférieur à la vraie valeur, de 183 unités du septième ordre.

La vraie valeur est en effet, dans cet exemple, $0,4396237\,;$ par la formule VI on trouve $0,4368730\,;$ par la formule VII, $0,4159824\,;$ enfin, par la formule VIII, $0,4051103\,;$ les deux dernières valeurs diffèrent tellement de la vérité qu’elles ne peuvent même tenir lieu d’approximations.

88

L’exposition de la seconde méthode fournira l’occasion de développer un grand nombre de relations nouvelles et élégantes ; comme elles prennent différentes formes suivant les diverses espèces de sections coniques, il sera convenable de traiter chacune d’elles en particulier ; nous commencerons par l’ELLIPSE.

Qu’à deux positions de l’astre correspondent les anomalies vraies $v,\,v'$ (dont $v$ répond à l’époque antérieure), les anomalies excentriques $\mathrm {E} ,$ $\mathrm {E} '$ et les rayons vecteurs $r,$ $r'\,;$ soient ensuite $p$ le demi-paramètre, $e=\sin \varphi$ l’excentricité, $a$ le demi grand axe, $t$ le temps pendant lequel l’astre passe de la première position à la seconde ; posons enfin, $v'-v=2f,$ $v'+v=2\mathrm {F} ,$ $\mathrm {E} '\!-\!\mathrm {E} =2g,$ $\mathrm {E} '\!+\!\mathrm {E} =2\mathrm {G} ,$ $a\cos \varphi ={\frac {p}{\cos \varphi }}=b.$ Ces conventions faites, nous déduisons facilement, par la combinaison des formules V, VI, art. 8, les équations suivantes :

[1]

$b\sin g=\sin f{\sqrt {rr'}},$

[2]

$b\sin \mathrm {G} =\sin \mathrm {F} {\sqrt {rr'}},$

$p\cos g=\left[\cos {\frac {1}{2}}v\cos {\frac {1}{2}}v'.(1\!+\!e)+\sin {\frac {1}{2}}v\sin {\frac {1}{2}}v'.(1\!-\!e)\right]{\sqrt {rr'}}$

ou

[3]

$p\cos g=(\cos f+e\cos \mathrm {F} ){\sqrt {rr'}},$

et de même

[4]

$p\cos \mathrm {G} =(\cos \mathrm {F} +e\cos f){\sqrt {rr'}}.$

De la combinaison des équations [3], [4], on obtient ensuite

[5]

$\cos f{\sqrt {rr'}}=(\cos g-e\cos \mathrm {G} )a,$

[6]

$\cos \mathrm {F} {\sqrt {rr'}}=(\cos \mathrm {G} -e\cos g)a.$

Au moyen de la formule III, art. 8, nous obtenons

[7]

$r'-r=2ae\sin g\sin \mathrm {G} .$

$r'+r=2a-2ae\cos g\cos \mathrm {G} =2a\sin ^{2}g+2\cos f\cos g{\sqrt {rr'}},$

d’où

[8]

$a={\frac {r'+r-2\cos f\cos g{\sqrt {rr'}}}{2\sin ^{2}g}}.$

Posons

[9]

${\frac {{\sqrt {\dfrac {r'}{r}}}+{\sqrt {\dfrac {r}{r'}}}}{2\cos f}}=1+2l,$

et l’on aura

[10]

$a={\frac {2\left(l+\sin ^{2}{\dfrac {1}{2}}g\right)\cos f{\sqrt {rr'}}}{\sin ^{2}g}}.$

et aussi

${\sqrt {\overset {}{a}}}=\pm {\frac {\sqrt {2\left(l+\sin ^{2}{\dfrac {1}{2}}g\right)\cos f{\sqrt {rr'}}}}{\sin g}},$

dans laquelle il faudra prendre le signe supérieur ou le signe inférieur, selon que $\sin g$ sera positif ou négatif. La formule XII de l’art. 8 nous donne l’équation

${\begin{aligned}{\frac {kt}{a^{\frac {3}{2}}}}&=\mathrm {E} '-e\sin \mathrm {E} '-\mathrm {E} +e\sin \mathrm {E} \\&=2g-2e\sin g\cos \mathrm {G} \\&=2g-\sin 2g+2\cos f\sin g{\frac {\sqrt {rr'}}{a}}.\end{aligned}}$

Si maintenant on substitue dans cette équation à la place de $a$ sa valeur [10], et qu’on pose par abréviation

[11]

${\frac {kt}{2^{\frac {3}{2}}\cos ^{\frac {3}{2}}\!f\,(rr')^{\frac {3}{2}}}}=m$,

on trouve, toutes réductions convenablement faites,

[12]

$\pm m=\left(l+\sin ^{2}{\frac {1}{2}}g\right)^{\frac {1}{2}}\!+\left(l+\sin ^{2}{\frac {1}{2}}g\right)^{\frac {3}{2}}\left({\frac {2g-\sin 2g}{\sin ^{3}g}}\right)$

dans laquelle le signe supérieur ou le signe inférieur doit être placé devant $m$, suivant que $\sin g$ est positif ou négatif.

Toutes les fois que le mouvement héliocentrique est compris entre 180° et 360°, ou plus généralement, toutes les fois que $\cos f$ est négatif, la quantité $m$ déterminée par la formule [11] devient imaginaire. Pour éviter cela, nous adopterons dans ce cas là, à la place des équations [9] et [11], les suivantes :

[9∗]

${\frac {{\sqrt {\dfrac {r'}{r}}}+{\sqrt {\dfrac {r}{r'}}}}{2\cos f}}=1-2\mathrm {L}$,

[11∗]

${\frac {kt}{2^{\frac {3}{2}}(-\cos f)^{\frac {3}{2}}(rr')^{\frac {3}{2}}}}=\mathrm {M}$,

d’où nous obtiendrons, à la place des formules [10] et (12]

[10∗]

$a={\dfrac {-2\left(\mathrm {L} -\sin ^{2}{\frac {1}{2}}g\right)\cos f{\sqrt {rr'}}}{\sin ^{2}g}}.$

[12∗]

$\;\;\pm \mathrm {M} =-\left(\mathrm {L} -\sin ^{2}{\frac {1}{2}}g\right)^{\frac {1}{2}}\!\!+\!\left(\mathrm {L} -\sin ^{2}{\frac {1}{2}}g\right)^{\frac {3}{2}}\!\left(\!{\frac {2g-\sin 2g}{\sin ^{3}g}}\!\right)$

dans laquelle, le signe à prendre est déterminé de la même manière que précédemment.

89

Un double travail se présente maintenant à nous : premièrement, que de l’équation transcendante [12], puisqu’elle n’admet pas de solution directe, nous déterminions l’inconnue $g$ le plus commodément ; secondement, que de l’angle $g$ trouvé nous déduisions les éléments eux-mêmes. Avant d’entreprendre ces questions, effectuons une transformation particulière au moyen de laquelle le calcul des quantités auxiliaires $l$ ou $\mathrm {L}$ est promptement achevé, et en outre plusieurs formules développées plus loin sont réduites à une forme plus élégante.

En introduisant en effet l’angle auxiliaire $\omega ,$ devant être déterminé par la relation

${\sqrt[{4}]{\frac {r'}{r}}}=\operatorname {tang} (45^{\circ }\!+\omega ),$

il vient

${\begin{aligned}{\sqrt {\frac {r'}{r}}}+{\sqrt {\frac {r}{r'}}}&=2+\left[\operatorname {tang} (45^{\circ }\!+\omega )-\operatorname {cotang} (45^{\circ }\!+\omega )\right]^{2}\\&=2+4\operatorname {tang} ^{2}{2\omega }\,;\end{aligned}}$

d’où l’on a

${\begin{aligned}l&={\frac {\sin ^{2}{\dfrac {1}{2}}f}{\cos f}}+{\frac {\operatorname {tang} ^{2}2\omega }{\cos f}}&\mathrm {L} &=-{\frac {\sin ^{2}{\dfrac {1}{2}}f}{\cos f}}-{\frac {\operatorname {tang} ^{2}2\omega }{\cos f}}\end{aligned}}.$

90

Nous considérerons d’abord le cas où, par la résolution de l’équation [12], on obtient pour $g$ une valeur qui n’est pas trop grande, de telle sorte que ${\frac {2g-\sin 2g}{\sin ^{3}g}}$ puisse être développé en série suivant les puissances croissantes de $\sin {\frac {1}{2}}g.$ Le numérateur de cette expression que nous désignons par $\mathrm {X} ,$ devient

$={\frac {32}{3}}\sin ^{3}{\frac {1}{2}}g-{\frac {16}{5}}\sin ^{5}{\frac {1}{2}}g-{\frac {4}{7}}\sin ^{7}{\frac {1}{2}}g-\dots \mathrm {etc.}$

et le dénominateur

$=8\sin ^{3}{\frac {1}{2}}g-12\sin ^{5}{\frac {1}{2}}g+3\sin ^{7}{\frac {1}{2}}g+\dots \mathrm {etc.}$

d’où $\mathrm {X}$ prend la forme

${\frac {4}{3}}+{\frac {8}{5}}\sin ^{2}{\frac {1}{2}}g+{\frac {64}{35}}\sin ^{4}{\frac {1}{2}}g+\dots \mathrm {etc.}$

Mais pour mettre en évidence la loi des coefficients de cette série, différentions l’équation

$\mathrm {X} \sin ^{3}g=2g-\sin 2g\,;$

d’où l’on déduit

$3\mathrm {X} \cos g\sin ^{2}g+\sin ^{3}g{\frac {d\mathrm {X} }{dg}}=2-2\cos 2g=4\sin ^{2}g.$

En posant en outre

$\sin ^{2}{\frac {1}{2}}g=x,$

on a

${\frac {dx}{dg}}={\frac {1}{2}}\sin g,$

d’où l’on conclut

${\frac {d\mathrm {X} }{dx}}={\frac {8-6\mathrm {X} \cos g}{\sin ^{2}g}}={\frac {4-3\mathrm {X} (1-2x)}{2x(1-x)}},$

et par conséquent

$(2x-2x^{2}){\frac {d\mathrm {X} }{dx}}=4-(3-6x)\mathrm {X} .$

Si donc, nous posons

$\mathrm {X} ={\frac {4}{3}}(1+\alpha x+\beta x^{2}+\gamma x^{3}+\delta x^{4}+\dots \mathrm {etc.} ),$

nous obtenons l’équation

${\begin{aligned}&{\frac {8}{3}}\left[\alpha x+(2\beta -\alpha )x^{2}+(3\gamma -2\beta )x^{3}+(4\delta -3\gamma )x^{4}+\,\mathrm {etc.} \right]\\={}&(8-4\alpha )x+(8\alpha -4\beta )x^{2}+(8\beta -4\gamma )x^{3}+(8\gamma -4\delta )x^{4}+\,\mathrm {etc.} \end{aligned}}$

qui doit être identique.

De là, nous concluons

$\alpha ={\frac {6}{5}},$  $\beta ={\frac {8}{7}}\alpha ,$  $\gamma ={\frac {10}{9}}\beta ,$  $\delta ={\frac {12}{11}}\gamma \dots \mathrm {etc.} ,$

d’où l’on déduit la loi de la progression. Nous avons donc

$\mathrm {X} ={\frac {4}{3}}+{\frac {4.6}{3.5}}x+{\frac {4.6.8}{3.5.7}}x^{2}+{\frac {4.6.8.10}{3.5.7.9}}x^{3}+{\frac {4.6.8.10.12}{3.5.7.9.11}}x^{4}+\mathrm {etc.}$

On peut transformer cette série en la fraction continue suivante :

$\mathrm {X} ={\cfrac {\cfrac {4}{3}}{1-{\cfrac {{\cfrac {6}{5}}x}{1+{\cfrac {{\cfrac {2}{5.7}}x}{1-{\cfrac {{\cfrac {5.8}{7.9}}x}{1-{\cfrac {{\cfrac {1.4}{9.11}}x}{1-{\cfrac {{\cfrac {7.10}{11.13}}x}{1-{\cfrac {{\cfrac {3.6}{13.15}}x}{1-{\cfrac {{\cfrac {9.12}{15.17}}x}{1-\mathrm {etc.} }}}}}}}}}}}}}}}}$

La loi suivant laquelle se forment les coefficients ${\frac {6}{5}},$ $-{\frac {2}{5.7}},$ ${\frac {5.8}{7.9}},$ ${\frac {1.4}{9.11}},$ etc., est évidente ; le $n$^ième terme de cette série, si $n$ est pair, est $={\frac {n-3n}{(2n+1)(2n+3)}},$ et si $n$ est impair, $={\frac {(n+2)(n+5)}{(2n+1)(2n+3)}}\,;$ un développement plus étendu de cette question serait trop étranger notre sujet.

Si nous posons maintenant

${\cfrac {x}{1+{\cfrac {{\cfrac {2}{5.7}}x}{1-{\cfrac {{\cfrac {5.8}{7.9}}x}{1-{\cfrac {{\cfrac {1.4}{9.11}}x}{1-\,\mathrm {etc.} ,}}}}}}}}=x-\xi$

on a

$\mathrm {X} ={\frac {1}{{\dfrac {3}{4}}-{\dfrac {9}{10}}(x-\xi )}},$  et  $\xi =x-{\frac {5}{6}}+{\frac {10}{9\mathrm {X} }},$

ou

$\xi ={\frac {\sin ^{2}g-{\dfrac {3}{4}}(2g-\sin 2g)\left(1-{\dfrac {6}{5}}\sin ^{2}{\dfrac {1}{2}}g\right)}{{\dfrac {9}{10}}(2g-\sin 2g)}}.$

Le numérateur de cette expression est une quantité du septième ordre, le dénominateur du troisième ordre, et par suite $\xi$ du quatrième, pourvu que $g$ soit considéré comme une quantité du premier ordre, ou $x$ comme une quantité du second ordre. On conclut de là que cette dernière formule n’est pas convenable pour la détermination numérique exacte de $\xi ,$ toutes les fois que $g$ n’exprime pas un angle très-considérable ; mais alors, les formules suivantes sont convenablement employées, formules qui ne diffèrent l’une de l’autre que par le changement des numérateurs dans les coefficients fractions, et dont la première se déduit facilement au moyen de la valeur supposée de $x-\xi$^[2].

[13]

$\xi ={\cfrac {{\cfrac {2}{35}}x^{2}}{1+{\dfrac {2}{35}}x-{\cfrac {{\cfrac {40}{63}}x}{1-{\cfrac {{\cfrac {4}{99}}x}{1-{\cfrac {{\cfrac {70}{143}}x}{1-{\cfrac {{\cfrac {18}{195}}x}{1-{\cfrac {{\cfrac {108}{255}}x}{1-\,{\text{etc.}}}}}}}}}}}}}}$

ou,

$\xi ={\cfrac {{\cfrac {2}{35}}x^{2}}{1-{\dfrac {18}{35}}x-{\cfrac {{\cfrac {4}{63}}x}{1-{\cfrac {{\cfrac {40}{99}}x}{1-{\cfrac {{\cfrac {18}{143}}x}{1-{\cfrac {{\cfrac {70}{195}}x}{1-{\cfrac {{\cfrac {40}{255}}x}{1-\,{\text{etc.}}}}}}}}}}}}}}$

Dans la troisième table annexée à cet ouvrage, se trouvent calculées avec sept décimales, les valeurs correspondantes de $\xi$ pour toutes les valeurs de $x,$ de millième en millième, depuis 0 jusqu’à 0,3. Cette table montre, au premier aspect, l’exiguïté de la valeur de $\xi$ pour des valeurs médiocres de $g\,;$ ainsi, pour $\mathrm {E} '\!-\!\mathrm {E} =10^{\circ }$ ou $g=5^{\circ },$ ce qui donne $x=0,00195,$ on a $\xi =0,0000002$. Il eût été superflu de continuer cette table au delà puisque le dernier terme $x=0,3$ répond à $g=66^{\circ }25',$ ou $\mathrm {E} '\!-\!\mathrm {E} =132^{\circ }50'.$ Enfin, la troisième colonne qui contient les valeurs de $\xi$ correspondant aux valeurs négatives de $x,$ sera expliquée plus loin, en son lieu.

91

L’équation 12, dans laquelle relativement au cas dont il s’agit il convient évidemment d’adopter le signe supérieur, prend, en introduisant la quantité $\xi ,$ la forme

$m=(l+x)^{\frac {1}{2}}+{\frac {(l+x)^{\frac {3}{2}}}{{\dfrac {3}{4}}-{\dfrac {9}{10}}(x-\xi )}}.$

En posant donc, ${\sqrt {(l+x)}}={\frac {m}{y}},$ et

[14]

${\frac {m^{2}}{{\dfrac {5}{6}}+l+\xi }}=h,$

on a, après toutes les réductions convenables,

[15]

$h={\frac {(y-1)y^{2}}{y+{\dfrac {1}{9}}}}.$

C’est pourquoi, si l’on peut regarder $h$ comme une quantité connue, $y$ s’en déduira à l’aide d’une équation du troisième degré, et l’on aura ensuite,

[16]

$x={\frac {m^{2}}{y^{2}}}-l.$

Maintenant, quoique $h$ contienne une quantité $\xi$ encore inconnue, on pourra, dans une première approximation, la négliger et prendre pour valeur de $h,$ ${\frac {m^{2}}{{\frac {5}{6}}+l}}$ puisque certainement, dans le cas que nous considérons, $\xi$ est toujours une quantité extrêmement petite.

De là, par les équations 15 et 16 on obtiendra $y$ et $x\,;$ de $x$ on déduira $\xi ,$ au moyen de la table III, et avec cette quantité on trouvera par la formule 14, la valeur de $h$ corrigée, avec laquelle recommençant le même calcul, on obtiendra les valeurs exactes de $y$ et $x.$ Le plus souvent ces valeurs diffèrent si peu des précédentes que $\xi$ pris de nouveau dans la table, n’est pas différent de la première valeur ; s’il en était autrement, il faudrait recommencer de nouveau le calcul jusqu’à ce qu’il ne déterminât plus aucun changement. Dès que la quantité $x$ sera trouvée, on aura $g$ par la formule,

$\sin ^{2}{\frac {1}{2}}g=x.$

Ces principes se rapportent au premier cas dans lequel $\cos f$ est positif ; dans l’autre cas, où il est négatif, nous posons

${\sqrt {\mathrm {L} -{\overset {}{x}}}}={\frac {\mathrm {M} }{\mathrm {Y} }},$

et

[14∗]

${\frac {\mathrm {M} ^{2}}{\mathrm {L} -{\dfrac {5}{6}}-\xi }}=\mathrm {H} ,$

d’où l’équation 12^∗, convenablement réduite, se change en celle-ci,

[15∗]

$\mathrm {H} ={\frac {(\mathrm {Y} +1)\mathrm {Y} ^{2}}{\mathrm {Y} -{\dfrac {1}{9}}}}.$

On pourra donc déterminer $\mathrm {Y}$ d’après $\mathrm {H} ,$ au moyen de cette équation cubique, d’où l’on déduira ensuite $x$ par l’équation

[16∗]

$x=\mathrm {L} -{\frac {\mathrm {M} ^{2}}{\mathrm {Y} ^{2}}}.$

Dans une première approximation, on pourra prendre pour $\mathrm {H} ,$ la valeur

${\frac {\mathrm {M} ^{2}}{\mathrm {L} -{\dfrac {5}{6}}}}\,;$

avec la valeur de $x$ déduite de là, par le moyen des équations 15^∗ et 16^∗, on extraira $\xi$ de la table III ; de là, par la formule 14^∗ on aura la valeur corrigée de $\mathrm {H}$ avec laquelle on recommencera le calcul de la même manière. Enfin, l’angle $g$ sera déterminé, d’après $x,$ de la même manière que dans le premier cas.

92

Quoique dans certains cas, les équations 15, 15^∗ puissent avoir trois racines réelles, il ne sera néanmoins jamais douteux laquelle dans notre problème devra être adoptée. Puisqu’en effet $h$ est évidemment une quantité positive, on conclut facilement de la théorie des équations, que l’équation 15 a une seule racine positive avec deux imaginaires ou deux négatives. Maintenant, puisque $y={\frac {m}{\sqrt {l+{\overset {}{x}}}}}$ doit être nécessairement une quantité positive, il est évident qu’il ne reste ici aucune incertitude.

Mais, en ce qui regarde l’équation 15^∗, nous observons d’abord que $\mathrm {L}$ est nécessairement plus grand que $1\,;$ ce qui est facilement démontré, si l’équation donnée dans l’art. 89 est mise sous la forme

$\mathrm {L} =1+{\frac {\cos ^{2}{\dfrac {1}{2}}f}{-\cos f}}+{\frac {\operatorname {tang} ^{2}2\omega }{-\cos f}}.$

En substituant ensuite, dans l’équation 12^∗, $\mathrm {Y} {\sqrt {\mathrm {L} -{\overset {}{x}}}}$ à la place de $\mathrm {M} ,$ il vient

$\mathrm {Y} +1=(\mathrm {L} -x)\mathrm {X} ,$

et, par suite,

$\mathrm {Y} \!+\!1>(1\!-\!x)\mathrm {X} >{\frac {4}{3}}+{\frac {4}{3.5}}x+{\frac {4.6}{3.5.7}}x^{2}+{\frac {4.6.8}{3.5.7.9}}x^{3}+...>{\frac {4}{3}},$

et par conséquent, $\mathrm {Y} >{\frac {1}{3}}.$ En posant, donc, $\mathrm {Y} ={\frac {1}{3}}+\mathrm {Y} ',$ $\mathrm {Y} '$ sera nécessairement une quantité positive ; par là aussi, l’équation 15^∗ se change en celle-ci

$\mathrm {Y} '^{3}+2\mathrm {Y} '^{2}+(1-\mathrm {H} )\mathrm {Y} '+{\frac {4}{27}}-{\frac {2}{9}}\mathrm {H} =0,$

qui ne peut avoir plusieurs racines positives, ainsi que le prouve facilement la théorie des équations. On conclut de là que l’équation 15^∗ ne peut avoir qu’une racine plus grande que ${\frac {1}{3}}$^[3] qu’il faudra, les autres étant négligées, adopter dans notre problème.

93

Pour rendre la solution de l’équation 15 la plus facile possible, pour les cas qui se présentent le plus fréquemment dans la pratique, nous ajoutons à la fin de cet ouvrage une table particulière (la table II) qui donne, pour les valeurs de $h$ comprises depuis $0$ jusqu’à $0,6$, les logarithmes correspondants de $y^{2}$ calculés avec le plus grand soin jusqu’à la septième décimale. L’argument $h$ est donné de dix-millième en dix-millième depuis $0$ jusqu’à $0,04$ ; par ce fait, les différences secondes de $\log y^{2}$ sont rendues insensibles, de sorte que dans cette partie de la table il suffit assurément d’une simple interpolation. Mais comme la table serait devenue beaucoup trop volumineuse, si elle avait eu partout le même développement, on a dû ne la donner depuis $h=0,04$ jusqu’à la fin, que de millième en millième ; c’est pourquoi il faudra, dans cette dernière partie, avoir égard aux différences secondes, si nous désirons à la vérité éviter des erreurs de quelques unités dans la septième décimale. Au reste, les petites valeurs de $h$ sont, dans la pratique, de beaucoup les plus fréquentes.

Toutes les fois que $h$ sort des limites de la table, la solution de l’équation 15, et aussi celle de l’équation 15^∗, pourront être effectuées sans difficulté par la méthode indirecte ou par d’autres méthodes assez connues. Au reste, il ne sera pas étranger à la question d’avertir qu’une petite valeur de $g$ ne peut exister avec une valeur négative de $\cos f$, si ce n’est dans les orbites très-excentriques, ainsi que cela ressortira spontanément de l’équation 20, donnée plus loin dans l’art. 95^[4].

94

La manière de traiter les équations 12, 12^∗, expliquée dans les art. 91, 92, 93, est basée sur la supposition que l’angle $g$ n’est pas trop grand, et certainement inférieur à la limite 66° 25′, au delà de laquelle nous n’avons pas étendu la table III. Toutes les fois que cette hypothèse n’est pas exacte, ces équations ne réclament pas d’artifices aussi grands ; car elles pourront toujours, sans changement de forme, être résolues par tâtonnements en toute sûreté et très-commodément. Sûrement, en effet, puisque la valeur de l’expression

${\frac {2g-\sin 2g}{\sin ^{2}g}}$,

dans laquelle $2g$ doit évidemment être exprimé en parties du rayon, peut être calculée pour les grandes valeurs de $g,$ avec toute précision, par le moyen des tables trigonométriques, ce qui ne peut certainement avoir lieu tant que $g$ est un petit angle ; commodément, puisque les lieux héliocentriques distants l’un de l’autre d’un aussi grand intervalle ne sont presque jamais employés pour la détermination d’une orbite encore entièrement inconnue, tandis que par le moyen des équations 1 ou 3 de l’art. 88, une valeur approchée de $g$ se déduit presque sans travail d’une connaissance quelconque de l’orbite ; enfin, d’une valeur approchée de $g,$ on obtiendra toujours, par un petit nombre d’essais, une valeur corrigée satisfaisant avec toute la précision désirable à l’équation 12 ou 12^∗. Au reste, toutes les fois que deux lieux héliocentriques donnés embrassent plus d’une révolution entière, il importe de se rappeler que pour l’anomalie excentrique tout autant de révolutions complètes auront été achevées, de telle sorte que les angles $\mathrm {E} '-\mathrm {E} ,$ $v'-v$ tomberont tous les deux entre 0 et 360° ou entre les mêmes multiples de la circonférence entière, et par suite $f$ et $g$ entre 0 et 180° ou entre les mêmes multiples de la demi-circonférence. Si enfin, l’orbite était entièrement inconnue, et qu’on ne pût même établir si le corps céleste, en allant du premier rayon vecteur au second, a décrit une partie seulement de la révolution, ou, en outre, une ou plusieurs révolutions, notre problème pourrait admettre quelquefois plusieurs solutions différentes ; toutefois nous ne nous arrêterons pas à ce cas qui ne se rencontre presque jamais dans la pratique.

95

Nous passons à une seconde question, à savoir la détermination des éléments d’après l’angle $g$ trouvé. Le demi-grand axe est ici obtenu immédiatement par les formules 10, 10^∗, à la place desquelles peuvent aussi être employées les suivantes :

[17]

$a={\frac {2m^{2}\cos f{\sqrt {rr'}}}{y^{2}\sin ^{2}\!g}}={\frac {k^{2}t^{2}}{4y^{2}rr'\cos ^{2}\!f\sin ^{2}\!g}},$

[17∗]

$a={\frac {-2\mathrm {M} ^{2}\cos f{\sqrt {rr'}}}{\mathrm {Y} ^{2}\sin ^{2}\!g}}={\frac {k^{2}t^{2}}{4\mathrm {Y} ^{2}rr'\cos ^{2}\!f\sin ^{2}\!g}}.$i

Le demi-petit axe $b={\sqrt {{\overset {}{a}}p}}$ est déterminé au moyen de l’équation 1 qui, étant combinée avec les précédentes, donne

[18]

$p=\left({\frac {yrr'\sin 2f}{kt}}\right)^{2},$

[18∗]

$p=\left({\frac {\mathrm {Y} rr'\sin 2f}{kt}}\right)^{2}.$

Maintenant, le secteur elliptique compris entre deux rayons vecteurs et un arc d’ellipse est ${\frac {1}{2}}kt{\sqrt {\overset {}{p}}},$ mais le triangle compris entre les mêmes rayons vecteurs et la corde $={\frac {1}{2}}rr'\sin 2f\,;$ c’est pourquoi le rapport du secteur au triangle est comme $y:1$ ou comme $\mathrm {Y} :1.$ Cette remarque est d’une très-grande importance et éclaire en même temps très-bien les équations 12 et 12^∗ ; de là, il est en effet évident que dans l’équation 12 les parties $m,$ $(l+x)^{\frac {1}{2}},$ $\mathrm {X} (l+x)^{\frac {3}{2}},$ et dans l’équation 12^∗ les parties $\mathrm {M} ,$ $(\mathrm {L} -x)^{\frac {1}{2}},$ $(\mathrm {L} -x)^{\frac {3}{2}},$ sont respectivement proportionnelles à l’aire du secteur (compris entre les rayons vecteurs et l’arc d’ellipse), à l’aire du triangle (entre les rayons vecteurs et la corde), à l’aire du segment (entre l’arc et la corde), puisque la première aire est évidemment égale à la somme ou à la différence des deux autres, selon que $v'-v$ tombe entre 0 et 180° ou entre 180° et 360°. Dans le cas où $v'-v$ est plus grand que 360°, il faut concevoir l’aire de l’ellipse entière ajoutée à l’aire du secteur et aussi à l’aire du segment autant de fois que le mouvement contient de révolutions entières.

Puisque $b=a\cos \varphi ,$ on trouve ensuite, par la combinaison des équations 1, 10, 10^∗,

[19]

$\cos \varphi ={\frac {\sin g\operatorname {tang} f}{2\left(l+\sin ^{2}\!{\frac {1}{2}}g\right)}},$

[19∗]

$\cos \varphi ={\frac {-\sin g\operatorname {tang} f}{2\left(\mathrm {L} -\sin ^{2}\!{\frac {1}{2}}g\right)}}\,;$

d’où, en substituant à la place de $l$ et $\mathrm {L}$ leurs valeurs de l’art. 89, il vient

[20]

$\cos \varphi ={\frac {\sin f\sin g}{1-\cos f\cos g+2\operatorname {tang} ^{2}\!2\omega }}.$

Cette formule n’est pas propre au calcul exact de l’excentricité toutes les fois que celle-ci a une petite valeur ; mais de cette relation on déduit facilement la formule suivante, qui est plus convenable,

[21]

$\operatorname {tang} ^{2}\!{\frac {1}{2}}\varphi ={\frac {\sin ^{2}\!{\dfrac {1}{2}}(f-g)+\operatorname {tang} ^{2}2\omega }{\sin ^{2}\!{\dfrac {1}{2}}(f+g)+\operatorname {tang} ^{2}2\omega }},$

à laquelle on peut aussi donner la forme suivante (en multipliant le numérateur et le dénominateur par $\cos ^{2}2\omega$),

[22]

$\operatorname {tang} ^{2}\!{\frac {1}{2}}\varphi ={\frac {\sin ^{2}\!{\dfrac {1}{2}}(f-g)+\cos ^{2}\!{\dfrac {1}{2}}(f-g)\sin ^{2}2\omega }{\sin ^{2}\!{\dfrac {1}{2}}(f+g)+\cos ^{2}\!{\dfrac {1}{2}}(f-g)\sin ^{2}2\omega }}.$

On pourra toujours déterminer l’angle $\varphi ,$ avec toute précision, au moyen de l’une ou l’autre formule (en employant si cela convient, les angles auxiliaires dont les tangentes sont ${\frac {\operatorname {tang} 2\omega }{\sin {\frac {1}{2}}(f-g)}},$ ${\frac {\operatorname {tang} 2\omega }{\sin {\frac {1}{2}}(f+g)}},$ pour la première, ou ${\frac {\sin 2\omega }{\operatorname {tang} {\frac {1}{2}}(f-g)}},$ ${\frac {\sin 2\omega }{\operatorname {tang} {\frac {1}{2}}(f+g)}}$ pour la dernière).

Pour la détermination de l’angle $\mathrm {G} ,$ on peut employer la formule suivante qui se déduit naturellement de la combinaison des équations 5, 7 et de celle qui suit non numérotée,

[23]

$\operatorname {tang} \mathrm {G} ={\frac {(r'-r)\sin g}{(r'+r)\cos g-2\cos f{\sqrt {rr'}}}},$

de laquelle, en introduisant $\omega ,$ on trouve facilement

[24]

$\operatorname {tang} \mathrm {G} ={\frac {\sin g\sin 2\omega }{\cos ^{2}\!2\omega \sin {\dfrac {1}{2}}(f\!-\!g)\sin {\dfrac {1}{2}}(f\!+\!g)+\sin ^{2}\!2\omega \cos g}}.$

L’ambiguïté qui se présente ici est facilement éludée par le secours de l’équation 7 qui apprend que l’on doit prendre $\mathrm {G}$ entre 0 et 180°, ou entre 180° et 360°, selon que le numérateur dans ces deux formules est positif ou négatif.

En combinant l’équation 3 avec celles-ci, qui découlent immédiatement de l’équation II art. 8,

${\begin{aligned}{\frac {1}{r}}-{\frac {1}{r'}}&={\frac {2e}{p}}\sin f\sin \mathrm {F} ,\\{\frac {1}{r}}+{\frac {1}{r'}}&={\frac {2}{p}}+{\frac {2e}{p}}\cos f\cos \mathrm {F} ,\end{aligned}}$

la formule suivante se déduira sans peine,

[25]

$\operatorname {tang} \mathrm {F} ={\frac {(r'-r)\sin f}{2\cos g{\sqrt {rr'}}-(r'+r)\cos f}},$

de laquelle, en introduisant l’angle $\omega ,$ il vient

[26]

$\operatorname {tang} \mathrm {F} ={\frac {\sin f\sin 2\omega }{\cos ^{2}\!2\omega \sin \!{\frac {1}{2}}(f\!-\!g)\sin \!{\frac {1}{2}}(f\!+\!g)-\sin ^{2}\!2\omega \cos f}}.$

L’ambiguïté est ici écartée comme précédemment. Aussitôt que les angles $\mathrm {F}$ et $\mathrm {G}$ auront été trouvés, on aura $v=\mathrm {F} -f,$ $v'=\mathrm {F} +f,$ d’où la position du périhélie sera connue ; et aussi $\mathrm {E} =\mathrm {G} -g,$ $\mathrm {E} '=\mathrm {G} +g.$ Enfin, le mouvement moyen pendant le temps $t,$ sera

${\frac {kt}{a^{\frac {3}{2}}}}=2g-2e\cos \mathrm {G} \sin g\,;$

l’accord de ces expressions servira à confirmer le calcul ; l’époque de l’anomalie moyenne correspondant à l’instant compris entre les deux époques proposées sera $\mathrm {G} -e\sin \mathrm {G} \cos g,$ laquelle pourra à volonté, être transportée à tout autre instant.

Il est encore un peu plus commode de calculer les anomalies moyennes pour les deux époques données par les formules $\mathrm {E} -e\sin \mathrm {E} ,$ $\mathrm {E} '-e\sin \mathrm {E} ',$ et d’employer leur différence, en la comparant à la quantité ${\frac {kt}{a^{\frac {3}{2}}}},$ à confirmer l’exactitude du calcul,

96

Les équations développées dans l’article précédent jouissent en vérité de tant de justesse qu’il semble qu’on ne peut rien désirer de plus. Néanmoins, on