3-37-1. Simon Newcomb to H. Poincaré

Washington le  Juin, 1890

À Monsieur H. Poincaré, Membre de l’Institut

Cher Monsieur,

Permetter me de vous indiquer un problème qui me paraît digne de votre admirable génie : c’est celui des inégalités qui peuvent être produites dans le mouvement de la lune par l’action des planètes. Les observations de la lune depuis 1670 montrent qu’il y a des inégalités dans ce mouvement notablement dans la longitude moyenne, qui ne sont pas produites par l’action du Soleil. Si elles sont réelles l’action des planètes est la seule cause à laquelle on peut les attribuer.

Je pense depuis vingt ans à une méthode d’aborder le problème que j’ai indiquée dans le Journal de M. Liouville pour 1870. Mais j’ai trouvé que l’application pratique de cette méthode demande des calculs analytiques et numériques tellement longues et compliqués que je n’ai pu les achever.

J’ai cherché vainement à trouver quelque transformation des équations du problème qui peut rendre plus facile les calculs nécessaires à sa solution. Voilà pourquoi je prends la liberté de attirer votre attention au sujet en espérant que, sans bien de difficulté, vous pouvez trouver une telle transformation ou au moins montrer qu’elle n’est pas possible.

Je pose

$X$, $Y$, $Z$, co-ordonnées rectangulaires du Soleil rapportées au centre de gravité de la terre et de la lune;

$x$, $y$, $z$, celles de la lune rapportées au centre de la terre;

$m_1$, $m_2$, $m_3$, $m_4$, les masses du Soleil, de la terre, de la lune, et d’une planète;

$x_4$, $y_4$, $z_4$, co-ordonnées de la planète rapportées au Soleil.

${\mu }_1$ $={\displaystyle \frac{m_1(m_2+m_3)}{m_1+m_2+m_3}}$

${\mu }_2$ $={\displaystyle \frac{m_2{m}_3}{m_2+m_3}}$

$\mathrm{\Delta }$ $=x(x_4+X)+y(y_4+Y)+z(z_4+Z)$

$\rho$ $={\left\{{(x_4+X)}^2+{(y_4+Y)}^2+{(z_4+Z)}^2\right\}}^{\frac{1}{2}}$

$R^2$ $=X^2+Y^2+Z^2$

$r^2$ $=x^2+y^2+z^2$

$r_4^2$ $=x_4^2+y_4^2+z_4^2$

$W$ $=\text{angle entre}R\text{et}r.$

La part de la fonction de force qui concerne le mouvement relatif des trois corps Soleil, Terre, Lune peut être réduite à la forme suivante:

$\mathrm{\Omega }$ $={\mathrm{\Omega }}_0+{\mathrm{\Omega }}_1+{\mathrm{\Omega }}_2$

$\mathrm{\Omega }$ $={\displaystyle \frac{m_1(m_2+m_3)}{R}}+{\displaystyle \frac{m_2{m}_3}{r}}$

${\mathrm{\Omega }}_1$ $={\displaystyle \frac{1}{2}}{m}_1{\mu }_2{\displaystyle \frac{r^2}{R^3}}(3{\cos }^{2}W-1)$

${\mathrm{\Omega }}_2$ $=m_4{\mu }_1\left({\displaystyle \frac{1}{\rho }}+{\displaystyle \frac{x_{4}X+y_{4}Y+z_{4}Z}{r_4^3}}\right)+m_4{\mu }_2\left({\displaystyle \frac{3}{2}}{\displaystyle \frac{{\mathrm{\Delta }}^2}{{\rho }^5}}-{\displaystyle \frac{1}{2}}{\displaystyle \frac{r^2}{{\rho }^3}}\right).$

Alors les six équations différentielles du mouvement relatif sont

${\mu }_1{\displaystyle \frac{d^{2}X}{d{t}^2}}$ $={\displaystyle \frac{\partial \mathrm{\Omega }}{\partial X}}$

et mêmes pour Y et Z,

${\mu }_2{\displaystyle \frac{d^{2}x}{d{t}^2}}$ $={\displaystyle \frac{\partial \mathrm{\Omega }}{\partial x}}$

et mêmes pour $y$ et $z$.

La solution demande trois approximations.

Première Approximation

On néglige ${\mathrm{\Omega }}_1$ et ${\mathrm{\Omega }}_2$. La solution est celle du mouvement elliptique : $x$, $y$, $z$ sont exprimés en fonction de six constantes arbitraires et $X$, $Y$, $Z$ en fonction de six autres constantes arbitraires.

Deuxième approximation

On pose $\mathrm{\Omega }={\mathrm{\Omega }}_0+{\mathrm{\Omega }}_1$; la solution est de la forme

$$x,y,z,X,Y,Z=f(a_1,a_2,\mathrm{\dots },a_6,a_1^{\prime },\mathrm{\dots },a_6^{\prime },t).$$ (1)

On peut choisir les douze constantes $a_1,\mathrm{\dots },a_6,a_1^{\prime },\mathrm{\dots },a_6^{\prime }$ de telle sorte que, quand ${\mathrm{\Omega }}_1$ s’évanouit, $a_1,\mathrm{\dots },a_6$ deviennent ceux du mouvement elliptique de $x$, $y$, $z$, et $a_1^{\prime },\mathrm{\dots },a_6^{\prime }$ celles du mouvement elliptique de $X$, $Y$, $Z$.

La partie la plus difficile de cette approximation est achevée par notre excellent ami feu M. Delaunay dans sa “Théorie du mouvement de la Lune.”

Troisième Approximation

On tient compte de ${\mathrm{\Omega }}_2$ en faisant varier les douze éléments (1) par la méthode de Lagrange. Maintenant, nous arrivons au point délicat du problème. Le travail de Delaunay nous donne les valeurs de $x$, $y$ et $z$ dans la forme

$$x,y,z=f(a_1,\mathrm{\dots },a_6,a_1^{\prime },\mathrm{\dots },a_6^{\prime },t).$$ (2)

La détermination des valeurs variables de $a_1^{\prime },\mathrm{\dots },a_6^{\prime }$ dues à l’action de $m_4$ n’offre aucune difficulté, et comme on peut censer ces variations comme actuellement existantes, on est porté à supposer que la méthode de solution la plus simple doit consister en regardant les douze éléments comme simultanément variables. Telle est la méthode que j’ai imaginée. Pour l’appliquer il faut former les valeurs de

${\displaystyle \frac{\partial {\mathrm{\Omega }}_2}{\partial a}}$ $={\displaystyle \frac{\partial {\mathrm{\Omega }}_2}{\partial x}}{\displaystyle \frac{\partial x}{\partial a}}+{\displaystyle \frac{\partial {\mathrm{\Omega }}_2}{\partial y}}{\displaystyle \frac{\partial y}{\partial a}}+{\displaystyle \frac{\partial {\mathrm{\Omega }}_2}{\partial z}}{\displaystyle \frac{\partial z}{\partial a}}$ (3)

$+{\displaystyle \frac{\partial {\mathrm{\Omega }}_2}{\partial X}}{\displaystyle \frac{\partial X}{\partial a}}+{\displaystyle \frac{\partial {\mathrm{\Omega }}_2}{\partial Y}}{\displaystyle \frac{\partial Y}{\partial a}}+{\displaystyle \frac{\partial {\mathrm{\Omega }}_2}{\partial Z}}{\displaystyle \frac{\partial Z}{\partial a}}.$ (4)

Maintenant, voilà la difficulté : Quoique les valeurs de $\frac{\partial X}{\partial a}$ soient moindres que celles $\frac{\partial x}{\partial a}$ dans le rapport

$$\frac{m_3}{m_2}\frac{r}{R}=1÷\mathrm{25\hspace{0.17em}000\hspace{0.17em}000},$$

cependant, $\partial {\mathrm{\Omega }}_2/\partial X$ est plus grand comme $\partial {\mathrm{\Omega }}_2/\partial x$ dans le même rapport; ainsi (3) et (4) sont du même ordre de grandeur. Aussi, le calcul analytique de $X$, $Y$, $Z$ en fonction de $a_1,\mathrm{\dots },a_6$ paraît demander un travail comparable avec celui de Delaunay; et je désire naturellement éviter ce travail.

Maintenant : voilà comment on peut éviter ce travail par ce qu’on peut appeler la méthode ordinaire, comme Delaunay l’a appliquée dans son calcul des inégalités de longue période produit par l’action de Vénus sur la lune (Appendix à la Conn. du Temps 1862?) On regarde $a_1^{\prime },\mathrm{\dots },a_6^{\prime }$ dans (2) comme des constantes numériquement données, $a_1,\mathrm{\dots },a_6$ comme seules variables et, en ${\mathrm{\Omega }}_2$, $X$, $Y$, $Z$ comme fonctions du temps numériquement données. Ainsi $\partial X/\partial a$ disparaît des équations et son calcul n’est pas nécessaire.

Les deux questions qui se présentent sont :

1. Cette seconde méthode est-elle bien rigoreuse?

2. Y a-t-il un rapprochement possible des deux méthodes qui combine la rigueur de l’une avec la facilité de l’autre?

ADft fragment, 13p. Simon Newcomb Papers 35, Library of Congress, USA.

Time-stamp: "28.01.2016 23:00"