# Annales de mathématiques pures et appliquées/Tome 16/Analise transcendante, article 3
## ANALISE TRANSCENDANTE.
### Exposition des principes du calcul des variations ;

I. Dans le calcul différentiel proprement dit, on ne fait varier, dans les fonctions qu’on y considère, que les seules variables $x,y,z\,;$ c’est-à-dire, les coordonnées de certaines courbes ou surfaces déterminées, ce qui répond à des points déterminés de ces mêmes surfaces, pour lesquels on enseigne ce que signifient les différentielles de ces coordonnées ; et on y enseigne également à déterminer les différentielles de toutes sortes de fonctions.
Mais, dans une fonction telle, par exemple, que
$y=\operatorname {f} (x,a,b,c,\ldots ),$
rien n’empêche de faire varier les paramètres $a,b,c,\ldots \,;$ et dès lors on sort de la courbe où l’on était d’abord, pour passer à une autre courbe voisine de celle-là. On peut obtenir ainsi une infinité de nouvelles courbes, suivant les valeurs diverses qu’on attribuera à $a,b,c,\ldots .$ On rencontre déjà des exemples de pareils changemens dans la recherche des solutions particulières des équations différentielles, où l’on différentie les équations des courbes par rapport à leurs paramètres, afin que les nouveaux coefficiens différentiels ainsi obtenus donnent, en les égalant à zéro, les solutions particulières demandées.
Si on ne fait varier que les paramètres, dans une courbe
$y=\operatorname {f} (x,a,b,c,\ldots )\,;$
cette courbe se trouvera transportée, suivant ses ordonnées, de $\mathrm {AB}$ en $\mathrm {A'B'}$ (fig. 1). On peut alors chercher l’augmentation de l’ordonnée ou de $y,$ qui est $\mathrm {MM'} \,;$ On peut aussi demander l’augmentation de d$dy$ qui de $\mathrm {TQ}$ est devenue $\mathrm {T'Q',}$ l’augmentation du rayon de courbure, etc. On pourrait bien se contenter, pour les questions de géométrie, de faire varier les courbes de la sorte, en suivant les ordonnées ; mais ce point de vue n’est pas assez général pour les questions de mécanique. On suppose alors que tous les points d’une courbe se transportent sur une autre courbe, en décrivant eux-mêmes des courbes intermédiaires $\mathrm {MM} ',mm'$ (fig. 2). Alors, en même temps qu’on fait varier les paramètres, on fait aussi varier l’abcisse. Ce sont les accroissemens que l’on donne, soit aux paramètres soit aux variables dépendantes, que l’on appelle leurs variations, et que l’on désigne par la caractéristique $\delta .$ Par analogie avec les dénominations admises dans le calcul différentiel, nous appellerons variation de la fonction la quantité à laquelle se réduit son accroissement total, lorsqu’on en supprime les termes affectés des puissances supérieures de $\delta x,\delta a,\delta b,\ldots .$ D’après cette définition, $u$ étant une fonction quelconque de $x,a,b,\ldots$ de même qu’on a
$\operatorname {d} u={\frac {\operatorname {d} u}{\operatorname {d} x}}\operatorname {d} x+{\frac {\operatorname {d} u}{\operatorname {d} a}}\operatorname {d} a+{\frac {\operatorname {d} u}{\operatorname {d} b}}\operatorname {d} b+\ldots ,$
on aura aussi
$\delta u={\frac {\operatorname {d} u}{\operatorname {d} x}}\delta x+{\frac {\operatorname {d} u}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} u}{\operatorname {d} b}}\delta b+\ldots .$
Il est inutile de mettre ${\frac {\delta u}{\delta x}},{\frac {\delta u}{\delta a}},\ldots$ au lieu de ${\frac {\operatorname {d} u}{\operatorname {d} x}},{\frac {\operatorname {d} u}{\operatorname {d} a}},\ldots \,;$ parce que ces derniers rapports, étant indépendans des accroissemens donnés à $x,a,\ldots ,$ restent les mêmes, quels que soient ces accroissemens.
La courbe $\mathrm {A'B',}$ sur laquelle on transporte la courbe $\mathrm {AB,}$ n’a pas besoin d’être de même forme que celle-ci, pour qu’on puisse prendre sur elle les variations comme nous venons de le dire. En effet, on peut développer en séries les ordonnées de ces deux courbes, comme on le voit ici
$\mathrm {AB} ,\qquad y=b\,+{\frac {b'}{1}}(x-\alpha )\ +{\frac {b''}{1.2}}(x-\alpha )^{2}+\ldots ,$
$\mathrm {A'B'} ,\quad y'=B+{\frac {B'}{1}}(x-\alpha )+{\frac {B''}{1.2}}(x-\alpha )^{2}+\ldots .$
Elles sont ainsi présentées sous la même forme, ayant un nombre infini de paramètres $b,b',b'',\ldots ,B,B',B'',\ldots .$ Or, ce que nous avons dit tout à l’heure, étant tout-à-fait indépendant du nombre des paramètres, rien n’empêche de supposer $B=b+\delta b,\ B'=b'+\delta b',\ B''=b''+\delta b'',\ldots ,$ et par conséquent de regarder la courbe $A'B'$ comme étant la variation de la première $\mathrm {AB}$.
Soient $\mathrm {PP} '=\delta x,\ \mathrm {P} p=\operatorname {d} x\,;$ on aura (fig. 2)
$\operatorname {d} \delta x=\mathrm {P} 'p'-\mathrm {P} p=(pp'+p\mathrm {P} ')-(\mathrm {PP} '+p\mathrm {P} ')=pp'-\mathrm {PP} ',$
$\delta \operatorname {d} x=pp'-\mathrm {PP} '$
d’où l’on conclura
$\operatorname {d} \delta x=\delta \operatorname {d} x\,;$
c’est-à-dire que la différentielle de la variation est égale à la variation de la différentielle.
Le même théorème a lieu pour les fonctions. En effet, soit
$y=\operatorname {f} (x,a,b,c,\ldots ),$
qui donne
$\operatorname {d} y={\frac {\operatorname {d} y}{\operatorname {d} x}}\operatorname {d} x\,;$
Si nous appliquons à cette équation une variation, en nous rappelant que, d’après la définition de la variation, semblable à celle de la différentielle,
$\delta .uv=v\delta u+u\delta v,$
nous aurons
$\delta \operatorname {d} y={\frac {\operatorname {d} y}{\operatorname {d} x}}\delta \operatorname {d} x+\operatorname {d} x\delta {\frac {\operatorname {d} y}{\operatorname {d} x}},$
c’est-à-dire,
$\delta \operatorname {d} y={\frac {\operatorname {d} y}{\operatorname {d} x}}\delta \operatorname {d} x+\operatorname {d} x\left({\frac {\operatorname {d} ^{2}y}{\operatorname {d} x^{2}}}\delta x+{\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} a}}\delta a+{\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} b}}\delta b+\ldots \right).$
Mais, d’un autre côté
$\delta y={\frac {\operatorname {d} y}{\operatorname {d} x}}\delta x+{\frac {\operatorname {d} y}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} y}{\operatorname {d} b}}\delta b+\ldots$
En différentiant cette dernière équation et observant que, puisque la différentielle se prend le long de la courbe, en ne faisant varier que $x$ et $y,$ il s’ensuit que $\delta a,\ \delta b,$ doivent ici être considérées comme des constantes, on aura
$\operatorname {d} \delta x={\frac {\operatorname {d} y}{\operatorname {d} x}}\operatorname {d} \delta x+\operatorname {d} x\left({\frac {\operatorname {d} ^{2}y}{\operatorname {d} x^{2}}}\operatorname {d} x+{\frac {\operatorname {d} ^{2}y}{\operatorname {d} a\operatorname {d} x}}\delta a+{\frac {\operatorname {d} ^{2}y}{\operatorname {d} b\operatorname {d} x}}\delta b+\ldots \right).$
En comparant entre eux les développemens de $\delta \operatorname {d} y$ et de $\operatorname {d} \delta y,$ en observant que, d’après ce qui a été prouvé ci-dessus,
${\frac {\operatorname {d} y}{\operatorname {d} x}}\delta \operatorname {d} x={\frac {\operatorname {d} y}{\operatorname {d} x}}\operatorname {d} \delta x$
et que d’ailleurs, parce que l’ordre des différentiations est indifférent,
${\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} a}}={\frac {\operatorname {d} ^{2}y}{\operatorname {d} a\operatorname {d} x}},\qquad {\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} b}}={\frac {\operatorname {d} ^{2}y}{\operatorname {d} b\operatorname {d} x}},\ldots$
on verra que ces deux développemens sont égaux, et qu’ainsi
$\operatorname {d} \delta y=\delta \operatorname {d} y.$
On peut tirer de ceci une conclusion semblable, par rapport au signe d’intégration. Soit, en effet
$\int v\operatorname {d} x=u,\quad$
$\quad v\operatorname {d} x=\operatorname {d} u,$
donc
$\delta .v\operatorname {d} x=\delta \operatorname {d} u=\operatorname {d} \delta u\,;$
d’où, en intégrant
$\int \delta .v\operatorname {d} x=\delta u=\delta \int v\operatorname {d} x.$
Ainsi, on peut intervertir l’ordre des signes $\int$ et $\delta .$
La courbe $\mathrm {AB}$ ayant été transportée sur $\mathrm {A'B'}$ (fig. 3), de manière que les points $\mathrm {M}$ et $m$ soient venus en $\mathrm {M} '$ et $m'\,;$ si nous menons les tangentes $\mathrm {MT}$ et $\mathrm {MS,}$ aux courbes $\mathrm {MB}$ et $\mathrm {MM'} \,;$ il est clair, d’après la définition que nous avons donnée de $\delta y,$ qui est, du reste,
$\delta y={\frac {\operatorname {d} y}{\operatorname {d} x}}\delta x+{\frac {\operatorname {d} y}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} y}{\operatorname {d} b}}\delta b+\ldots \,;$
il est clair, disons nous, que $\delta y=\mathrm {SQ} \,;$ puisque ce doit être la partie de $\mathrm {M'Q}$ qui ne contient pas les termes infiniment petits par rapport aux premiers. Mais, si l’on suppose que la position $\mathrm {A'B'}$ soit infiniment voisine de la position $\mathrm {AB} ,$ on aura, à la limite,
$\delta y=\mathrm {M'Q} ,\quad {\frac {\operatorname {d} y}{\operatorname {d} x}}\delta x\quad$
$\quad p\delta x=\mathrm {HQ} ,$
donc
$\delta y-p\delta x=\mathrm {M'Q-HQ} .$
De même
${\frac {\operatorname {d} y}{\operatorname {d} x}}\delta x=p\delta x=\mathrm {TQ}$
sera, à la limite, la distance des deux courbes. Or, on sait déjà que
$\delta y-p\delta x={\frac {\operatorname {d} y}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} y}{\operatorname {d} b}}\delta b+\ldots$
Nous représenterons désormais par $\omega$ cette quantité, dépendant seulement de la variation des paramètres ; quantité qui, lorsque les variations deviennent infiniment petites, représente, ainsi que nous venons de le voir, la distance entre les deux courbes. Nous aurons donc ainsi
$\delta y=p\delta x+\omega \,;$
c’est cet $\omega ,$ qui reste seul, quand on suppose $\delta x=0,$ et que les paramètres seuls varient.
Comme $p$ est une fonction de $x,$ de même qu’on a établi
$\delta y=p\delta x+\omega ,$
on peut poser aussi
${\begin{aligned}&\delta p=q\delta x+\omega ',\\&\delta q=r\delta x+\omega '',\\&\ldots \ldots \ldots \ldots \ldots \end{aligned}}$
Enfin si $\mathrm {V}$ est une fonction de $x,y,p,q,\ldots \,;$ comme on peut toujours exprimer $y,p,q,\ldots$ en fonction de $x$ seul, et par conséquent $\mathrm {V}$ aussi ; on aura de même
$\delta \mathrm {V} =\left[{\frac {\operatorname {d} \mathrm {V} }{\operatorname {d} x}}\right]\delta x+\Omega \,;$
où $\left[{\frac {\operatorname {d} \mathrm {V} }{\operatorname {d} x}}\right]$ représente la différentielle de $\mathrm {V}$ par rapport aux $x$ non seulement en évidence, mais même contenus dans $y,p,q,\ldots .$
Il existe des relations remarquables entre ces $\omega ,\omega ',\omega '',\ldots \Omega .$
Soit, en effet, comme ci-dessus
$y=\operatorname {f} (x,a,b,c,\ldots )\quad$
$\quad {\frac {\operatorname {d} y}{\operatorname {d} x}}=p,$
d’où
$\delta y={\frac {\operatorname {d} y}{\operatorname {d} x}}\delta x+{\frac {\operatorname {d} y}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} y}{\operatorname {d} b}}\delta b+\ldots$
et
$\omega =\delta y-p\delta x={\frac {\operatorname {d} y}{\operatorname {d} a}}\delta a+{\frac {\operatorname {d} y}{\operatorname {d} b}}\delta b+\ldots$
On entend par ${\frac {\operatorname {d} \omega }{\operatorname {d} x}}$ la dérivée de cette fonction de $x,a,b,\ldots \delta a,\delta b,\ldots$ en n’y faisant varier que $x,$ parce que
$\delta a=a'-a,\qquad \delta b=b'-b,\ldots$
sont constantes, quand $a,b,\ldots a',b',\ldots$ le sont. Ainsi
${\frac {\operatorname {d} \omega }{\operatorname {d} x}}={\frac {\operatorname {d} ^{2}y}{\operatorname {d} a\operatorname {d} x}}\delta a+{\frac {\operatorname {d} ^{2}y}{\operatorname {d} b\operatorname {d} x}}\delta b+\ldots .$
D’un autre côté, $p={\frac {\operatorname {d} y}{\operatorname {d} x}}$ étant une fonction de $x,a,b,\ldots$ on a
$\delta p={\frac {\operatorname {d} ^{2}y}{\operatorname {d} x^{2}}}\delta x+{\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} a}}\delta a+{\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} b}}\delta b+\ldots \,;$
en se rappelant donc que
${\frac {\operatorname {d} ^{2}y}{\operatorname {d} x^{2}}}={\frac {\operatorname {d} p}{\operatorname {d} x}}=q,$
et qu’en outre
${\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} a}}={\frac {\operatorname {d} ^{2}y}{\operatorname {d} a\operatorname {d} x}},\qquad {\frac {\operatorname {d} ^{2}y}{\operatorname {d} x\operatorname {d} b}}={\frac {\operatorname {d} ^{2}y}{\operatorname {d} b\operatorname {d} x}},\ldots$
on trouve
$\delta p=q\delta x+{\frac {\operatorname {d} \omega }{\operatorname {d} x}}.$
On trouvera de même
$\delta q=r\delta x+{\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}}$
d’où l’on voit que
$\omega '={\frac {\operatorname {d} \omega }{\operatorname {d} x}},\qquad \omega ''={\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}},\ldots$
Pour $\mathrm {V} =F(x,y,p,q,\ldots ),$ où les paramètres $a,b,c,\ldots$ entrent implicitement, on a
$\delta \mathrm {V} =\left[{\frac {\operatorname {d} \mathrm {V} }{\operatorname {d} x}}\right]\delta x+\left[{\frac {\operatorname {d} \mathrm {V} }{\operatorname {d} a}}\right]\delta a+\left[{\frac {\operatorname {d} \mathrm {V} }{\operatorname {d} b}}\right]\delta b+\ldots \,;$
mais, en désignant par $\mathrm {M,N,P,Q} ,\ldots$ les coefficiens différentiels de $\mathrm {V}$ par rapport aux $x,y,p,q,\ldots$ en évidence, on a
$\operatorname {d} V=M\operatorname {d} x+N\operatorname {d} y+P\operatorname {d} p+Q\operatorname {d} q+\ldots ,$
$\left[{\frac {\operatorname {d} \mathrm {V} }{\operatorname {d} x}}\right]=M+Np+Pq+Qr+\ldots \,;$
et il en résulte
$\delta V=M\delta x+N\delta y+P\delta p+Q\delta q+\ldots .$
Si nous substituons pour $\delta y,\delta p,\delta q,\ldots$ les valeurs précédemment obtenues, il en résultera
$\delta V=(M+Np+Pq+Qr+\ldots )\delta x+N\omega +P{\frac {\operatorname {d} \omega }{\operatorname {d} x}}+Q{\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}}+\ldots \,;$
valeur qui, comparée à
$\delta V=\left[{\frac {\operatorname {d} V}{\operatorname {d} x}}\right]\delta x+\left[{\frac {\operatorname {d} V}{\operatorname {d} a}}\right]\delta a+\left[{\frac {\operatorname {d} V}{\operatorname {d} b}}\right]\delta b+\ldots =\left[{\frac {\operatorname {d} V}{\operatorname {d} x}}\right]\delta x+\Omega ,$
fait voir que
$\Omega =N\omega +P{\frac {\operatorname {d} \omega }{\operatorname {d} x}}+Q{\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}}+\ldots$
La variation de $V,$ savoir :
$\delta V=(M+Np+Pq+Qr+\ldots )\delta x+N\omega +P{\frac {\operatorname {d} \omega }{\operatorname {d} x}}+Q{\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}}+\ldots$
une fois trouvée, celle de l’intégrale $\int V\operatorname {d} x$ est facile à en déduire. On sait en effet que
$\delta \int V\operatorname {d} x=\int \delta (V\operatorname {d} x)=\int (\delta V.\operatorname {d} x+V\delta \operatorname {d} x)=\int (V\operatorname {d} \delta x+\delta V.\operatorname {d} x),$
substituant pour $\delta V$ la valeur que nous venons d’obtenir, on aura
$\delta \int V\operatorname {d} x=$
$\int \left\{V\operatorname {d} \delta x+\left[{\frac {\operatorname {d} V}{\operatorname {d} x}}\right]\operatorname {d} x\delta x+N\omega \operatorname {d} x+P{\frac {\operatorname {d} \omega }{\operatorname {d} x}}\operatorname {d} x+Q{\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} ^{2}x}}+\ldots \right\}.$
Or, en premier lieu,
$\int \left(V\operatorname {d} \delta x+\left[{\frac {\operatorname {d} V}{\operatorname {d} x}}\right]\operatorname {d} x\delta x\right)=V\delta x\,;$
on trouve ensuite, en intégrant par parties,
${\begin{aligned}\int P{\frac {\operatorname {d} \omega }{\operatorname {d} x}}\operatorname {d} x&=P\omega -\int {\frac {\operatorname {d} P}{\operatorname {d} x}}\omega \operatorname {d} x,\\\\\int Q{\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}}\operatorname {d} x&=Q{\frac {\operatorname {d} \omega }{\operatorname {d} x}}-{\frac {\operatorname {d} Q}{\operatorname {d} x}}\omega +\int {\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}\omega \operatorname {d} x,\\\\\int R{\frac {\operatorname {d} ^{3}\omega }{\operatorname {d} x^{3}}}\operatorname {d} x&=R{\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}}-{\frac {\operatorname {d} R}{\operatorname {d} x}}.{\frac {\operatorname {d} \omega }{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}\omega -\int {\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}\omega \operatorname {d} x,\end{aligned}}$
réunissant alors tous ces termes, on aura finalement
$\delta \int V\operatorname {d} x=V\delta x+\left(P-{\frac {\operatorname {d} Q}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\ldots \right)\omega +\left(Q-{\frac {\operatorname {d} R}{\operatorname {d} x}}+\ldots \right){\frac {\operatorname {d} \omega }{\operatorname {d} x}}+\ldots$
$+\int \left(N-{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}-{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}+\ldots \right)\omega \operatorname {d} x.$
La dernière formule que nous venons d’obtenir sert à résoudre une des questions les plus importantes du calcul des variations, celle des maxima et minima. Une partie de cette question ressort du calcul différentiel ; mais cette partie est peu étendue ; elle ne comprend que le cas où l’on demande le maximum ou le minimum d’une fonction dans laquelle on ne fait varier que les variables courantes $x,y,\ldots .$ Ainsi, ce sera la détermination de l’ordonnée maximum ou minimum d’une certaine courbe ou d’une certaine surface. Mais si, dans une fonction, on fait varier les constantes ou paramètres, on peut se demander quelles sont les valeurs à leur attribuer, pour que la fonction jouisse d’une certaine propriété à un degré maximum ou minimum. Cette importante question a beaucoup occupé les géomètres du dernier siècle. Euler l’a résolue, en partie, dans un ouvrage ayant pour titre : Methodus inveniendi, etc.
Or, la série de Taylor donne, comme on sait, dans le cas de plusieurs variables, et conséquemment pour
$u=\operatorname {f} (x,y,p,q,\ldots a,b,c,\ldots ),$
$\operatorname {f} \left(x+\delta x,\ y+\delta y,\ p+\delta p,\ q+\delta q,\ldots a+\delta a,\ b+\delta b,\ c+\delta c,\ldots \right)$
$=u+{\frac {\operatorname {d} u}{\operatorname {d} x}}\delta x+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} x^{2}}}{\frac {\delta x^{2}}{1.2}}+\ldots$
$+{\frac {\operatorname {d} u}{\operatorname {d} y}}\delta y+2{\frac {\operatorname {d} ^{2}u}{\operatorname {d} x\operatorname {d} y}}.{\frac {\delta x\delta y}{1.2}}+\ldots$
${\begin{aligned}&+{\frac {\operatorname {d} u}{\operatorname {d} p}}\delta p+{\frac {\operatorname {d} ^{2}u}{\operatorname {d} y^{2}}}.{\frac {\delta y^{2}}{1.2}}+\ldots \\\\&+\ldots \ldots +2{\frac {\operatorname {d} ^{2}u}{\operatorname {d} x\operatorname {d} p}}.{\frac {\delta p\delta x}{1.2}}+\ldots \\\\&+{\frac {\operatorname {d} u}{\operatorname {d} a}}\delta a+\ldots \ldots \ldots +\ldots \\\\&+{\frac {\operatorname {d} u}{\operatorname {d} b}}\delta b+2{\frac {\operatorname {d} ^{2}u}{\operatorname {d} x\operatorname {d} a}}.{\frac {\delta x\delta a}{1.2}}+\ldots \\&+\ldots \ldots +\ldots \ldots \ldots +\ldots \end{aligned}}$
ou bien
$U=u+{\frac {\delta u}{1}}+{\frac {\delta ^{2}u}{1.2}}+{\frac {\delta ^{3}u}{1.2.3}}+\ldots ,$
en désignant par $\delta ^{2}u$ la variation de $\delta u,$ par $\delta ^{3}u$ celle de $\delta ^{2}u$ et ainsi de suite.
Or, si $\mathrm {U}$ est un maximum ou un minimum, on sait que $\delta u$ doit être nul, puisqu’il ne contient que les premières puissances des variations $\delta x,\ \delta y,\ldots ,$ et change conséquemment de signe avec elles. Au contraire $\delta ^{2}u$ n’en devra pas changer, puisqu’il contient ces variations à la seconde puissance.
Quand on demande le maximum ou le minimum d’une fonction, on ne peut souvent exprimer cette fonction que par des intégrales définies, entre les limites où la fonction doit avoir lieu. Il faut alors considérer le développement de la série de Taylor relatif à $\int V\operatorname {d} x,$ savoir :
$\int V\operatorname {d} x+{\frac {\delta \int V\operatorname {d} x}{1}}+{\frac {\delta ^{2}\int V\operatorname {d} x}{1.2}}+{\frac {\delta ^{3}\int V\operatorname {d} x}{1.2.3}}+\ldots$
S’il y a maximum ou minimum, il faut que $\delta \int V\operatorname {d} x=0.$ Les valeurs de $\delta x,\delta y,\ldots$ qui satisfont à cette condition, devront rendre, en outre, $\delta ^{2}\int V\operatorname {d} x$ négatif dans le cas du maximum, et positif dans le cas du minimum.
Il est encore à remarquer que, si le maximum ou le minimum doit avoir lieu entre certaines limites, non seulement $\int V\operatorname {d} x$ mais aussi
$\delta \int V\operatorname {d} x,\qquad \delta ^{2}\int V\operatorname {d} x,\qquad \delta ^{3}\int V\operatorname {d} x,\ldots$
devront être pris entre ces limites. Or nous avons trouvé ci-dessus
$\delta \int v\operatorname {d} x=V\delta x+\left(P-{\frac {\operatorname {d} Q}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\ldots \right)\omega +\left(Q-{\frac {\operatorname {d} R}{\operatorname {d} x}}+\ldots \right){\frac {\operatorname {d} \omega }{\operatorname {d} x}}+\ldots$
$+\int \left(N-{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}-\ldots \right)\omega \operatorname {d} x\,;$
et si nous nous rappelons que
$\omega =\delta y-p\delta x,\quad {\frac {\operatorname {d} \omega }{\operatorname {d} x}}=\delta p-q\delta x,\quad {\frac {\operatorname {d} ^{2}\omega }{\operatorname {d} x^{2}}}=\delta q-r\delta x,\ldots \,;$
il viendra, en substituant,
$0=\delta \int V\operatorname {d} x=$
$\left\{V-P\left(P-{\frac {\operatorname {d} Q}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\ldots \right)-q\left(Q-{\frac {\operatorname {d} R}{\operatorname {d} x}}+\ldots \right)-r(R-\ldots )-\ldots \right\}\delta x$
$+\left(P-{\frac {\operatorname {d} Q}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\ldots \right)\delta y+\left(Q-{\frac {\operatorname {d} R}{\operatorname {d} x}}+\ldots \right)\delta p+(R-\ldots )+\ldots$
$+\int \left(N-{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}-{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}+\ldots \right)\omega \operatorname {d} x,$
qu’il faudra prendre entre les limites données.
Il se présente ici plusieurs cas à examiner.
1.$^o$ Supposons d’abord que l’on donne, pour les deux limites, les valeurs $x',x'',y',y'',p',p'',\ldots$ de $x,y,p,\ldots .$ Alors, ces valeurs extrêmes étant fixes, leurs variations seront nulles, de sorte que l’équation à laquelle nous venons de parvenir se réduira simplement à
$\int \left(N-{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}-{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}+\ldots \right)\omega \operatorname {d} x=0.$
Or, on ne peut supposer $\omega =0,$ parce que, sous le signe $\int ,$ c’est le $\omega$ courant, et qu’il ne peut être nul qu’aux limites ; il faut donc que
$N-{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}-{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}+\ldots =0\,;\qquad$
de sorte que cette équation sera l’équation différentielle du maximum ou du minimum cherché. Mais si
$V=\operatorname {f} (x,y,p,q,\ldots )$
contient un coefficient différentiel de l’ordre $n,$ il est clair que cette équation sera de l’ordre $2n,$ de sorte qu’en l’intégrant elle donnera $2n$ constantes arbitraires ; or ; comme aux deux limites $x',x'',$ il faut qu’elle reproduise $y',y'',p',p'',q',q'',\ldots$ en nombre $2n,$ on aura $2n$ relations pour déterminer les $2n$ constantes arbitraires introduites par l’intégration ; et on le pourra aisément.
2.$^o$ Supposons qu’on ne donne que les deux limites $x',x'',y',y'',$ sans assigner, pour ces limites, les valeurs de $p',p'',q',q'',\ldots .$ Il faudra toujours, dans le développement de $\delta \int V\operatorname {d} x,$ égaler séparément à zéro la partie contenue sous le signe $\int$ et la partie en avant de ce signe. En effet, supposons que nous ayons trouvé la courbe qui satisfait au maximum ou au minimum, entre les deux points donnés, et prenons, pour les points limites, les valeurs de $p',p'',q',q'',\ldots \,;$ cherchons ensuite la courbe qui satisfait au maximum ou au minimum, pour ces valeurs fixes de $p,q,\ldots \,;$ il est manifeste que nous devrons retomber de nouveau sur la courbe dont il s’agit ; or, en fixant $p',p'',q',q'',\ldots$ tout aussi bien que $x',x'',y',y'',$ il faut, comme nous l’avons vu ci-dessus, que l’on ait l’équation (A) ; donc, cette équation devra encore avoir lieu, lorsqu’on ne fixera pas $p',p'',q',q'',\ldots$. Comme ces quantités sont indépendantes entre elles, pour que la première partie de $\delta \int V\operatorname {d} x$ soit nulle, il faudra égaler séparément à zéro les coefficiens
$Q'-{\frac {\operatorname {d} R'}{\operatorname {d} x'}}+{\frac {\operatorname {d} ^{2}S'}{\operatorname {d} x'^{2}}}-\ldots$
$\delta p',\quad Q''-{\frac {\operatorname {d} R''}{\operatorname {d} x''}}+{\frac {\operatorname {d} ^{2}S''}{\operatorname {d} x''^{2}}}-\ldots$
$\delta p'',$
$R'-{\frac {\operatorname {d} S'}{\operatorname {d} x'}}+\ldots \ldots \ldots \ldots$
$\delta q',\quad R''-{\frac {\operatorname {d} S''}{\operatorname {d} x''}}+\ldots \ldots \ldots \ldots$
$\delta q'',$
$S'-\ldots \ldots \ldots \ldots \ldots \ldots \ \$
$\delta r',\quad S''-\ldots \ldots \ldots \ldots \ldots \ldots \quad$
$\delta r'',$
Les équations résultantes ne seront qu’au nombre de $2n-2,$ parce que le $\delta$ du dernier coefficient différentiel manque dans la première partie de $\delta \int V\operatorname {d} x,$ comme le fait voir bien évidemment la forme des coefficiens de $\delta p,\delta q,\ldots .$ Par exemple, s’il n’y a que cinq coefficiens différentiels $p,q,r,s,t,$ la première partie de $\delta \int V\operatorname {d} x,$ finira par $T\delta s,$ et ne contiendra point $\delta t.$ Mais il faut qu’aux deux limites l’intégrale de l’équation donne pour $x=x'$ ou $x'',\ y=y'$ ou $y''.$ Cela fournira deux nouvelles relations qui, jointes aux $2n-2$ dont il vient d’être question, compléteront le nombre total de celles qui seront nécessaires pour la détermination des $2n$ constantes arbitraires introduites par l’intégration de (A).
Il est presque superflu d’observer que, si quelques-unes seulement des quantités $p',p'',q',q'',$ étaient données, il n’y aurait à égaler à zéro que les coefficiens des $\delta p',\ \delta p'',\ \delta q',\ \delta q'',\ldots$ non donnés. Les autres relations nécessaires pour déterminer les $2n$ constantes arbitraires de l’intégrale de (A), seraient fournies par celles des quantités $p',p'',q',q'',\ldots$ dont les valeurs seraient connues.
Nous remettons à traiter plus tard du cas où les limites ne seraient pas invariablement fixées, mais assujetties seulement à la condition de se trouver sur deux lieux géométriques donnés ; et nous nous occuperons, pour le présent, des cas où l’équation (A) est susceptible de quelques simplifications.
Il est de cas où l’équation
$N-{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}-{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}+\ldots =0\qquad$
peut être simplifiée.
1.$^o$ Quand $V$ ne contient pas $y,$ on a $N=0,$ et par suite
${\frac {\operatorname {d} P}{\operatorname {d} x}}-{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}+{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}-\ldots =0\,;$
ce qui donne, en multipliant par $\operatorname {d} x$ et intégrant,
$P-{\frac {\operatorname {d} Q}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\ldots =C\,;\qquad$
équation de l’ordre $2n-1.$ Si en outre $p$ manquait dans $V,$ l’équation pourrait, par un procédé analogue, être amenée à l’ordre $2n-2,$ et ainsi de suite.
2.$^o$ Quand $V$ ne contient pas $x,$ et que l’on a conséquemment
$V=\operatorname {f} (y,p,q,\ldots ),$
et par suite
$\left[{\frac {\operatorname {d} V}{\operatorname {d} x}}\right]=Np+Pq+Qr+\ldots .$
Si nous substituons ici la valeur de $N$ tirée de l’équation (A), qui est
$N={\frac {\operatorname {d} P}{\operatorname {d} x}}-{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}+{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}-\ldots ,$
nous aurons
${\frac {\operatorname {d} V}{\operatorname {d} x}}=Pq+Qr+Rs+\ldots +p{\frac {\operatorname {d} P}{\operatorname {d} x}}-p{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}+p{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}-\ldots \,;$
multipliant par $\operatorname {d} x,$ et observant que
$q\operatorname {d} x=\operatorname {d} p,\qquad r\operatorname {d} x=\operatorname {d} q,\qquad s\operatorname {d} x=\operatorname {d} r,\ldots ,$
il viendra
${\frac {\operatorname {d} V}{\operatorname {d} x}}\operatorname {d} x=P\operatorname {d} p+Q\operatorname {d} q+R\operatorname {d} r+\ldots +p{\frac {\operatorname {d} P}{\operatorname {d} x}}\operatorname {d} x-p{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}\operatorname {d} x+p{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}\operatorname {d} x-\ldots$
et par conséquent
$V=\int P\operatorname {d} p+\int Q\operatorname {d} q+\int R\operatorname {d} r+\ldots$
$+\int p{\frac {\operatorname {d} P}{\operatorname {d} x}}\operatorname {d} x-\int p{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}\operatorname {d} x+\int p{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}\operatorname {d} x-\ldots$
Or, en intégrant par parties, on trouve
${\begin{aligned}&\int P\operatorname {d} p=Pp-\int p{\frac {\operatorname {d} P}{\operatorname {d} x}}\operatorname {d} x,\\\\&\int Q\operatorname {d} q=Qq-p{\frac {\operatorname {d} Q}{\operatorname {d} x}}+\int p{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}\operatorname {d} x,\\\\&\int R\operatorname {d} r=Rr-q{\frac {\operatorname {d} R}{\operatorname {d} x}}+p{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\int p{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}\operatorname {d} x,\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \,;\end{aligned}}$
ce qui donne, en subtituant,
$V=p\left(P-{\frac {\operatorname {d} Q}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\ldots \right)+q\left(Q-{\frac {\operatorname {d} R}{\operatorname {d} x}}+\ldots \right)$
$+r(R-\ldots )+\ldots +C\,;\quad$
qui n’est plus que de l’ordre $2n-1.$
3.$^o$ Quand $V$ ne contient ni $x$ ni $y,$ on peut abaisser l’équation (A), qu’on appelle quelquefois l’équation indéfinie, parce que son intégrale donne l’équation du maximum ou du minimum, avec des constantes arbitraires. On a alors, à la fois,
$M=0,\qquad N=0,$
${\frac {\operatorname {d} V}{\operatorname {d} x}}Pq+Qr+Rs+\ldots \,;$
l’équation indéfinie se réduit donc à
${\frac {\operatorname {d} P}{\operatorname {d} x}}-{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}+{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}-\ldots =0,$
et donne, en intégrant
$P={\frac {\operatorname {d} Q}{\operatorname {d} x}}-{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}+\ldots +C\,;$
puis, en substituant
${\frac {\operatorname {d} V}{\operatorname {d} x}}=Qr+Rs+\ldots +q{\frac {\operatorname {d} Q}{\operatorname {d} x}}-q{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}+\ldots +Cq\,;$
et ensuite, en multipliant par $\operatorname {d} x$ et intégrant,
$V=\int Q\operatorname {d} q+\int R\operatorname {d} r+\ldots +\int q{\frac {\operatorname {d} Q}{\operatorname {d} x}}\operatorname {d} x-\int q{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}\operatorname {d} x+\ldots +Cp+C'\,;$
mais, nous avons trouvé ci-dessus, en intégrant par parties
${\begin{aligned}&\int Q\operatorname {d} q=Qq-p{\frac {\operatorname {d} Q}{\operatorname {d} x}}+\int p{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}\operatorname {d} x,\\\\&\int R\operatorname {d} r=Rr-q{\frac {\operatorname {d} R}{\operatorname {d} x}}+p{\frac {\operatorname {d} ^{2}R}{\operatorname {d} x^{2}}}-\int p{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}\operatorname {d} x,\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \,;\end{aligned}}$
il viendra donc, en substituant,
$V=q\left(Q-{\frac {\operatorname {d} R}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}S}{\operatorname {d} x^{2}}}-\ldots \right)+r\left(R-{\frac {\operatorname {d} S}{\operatorname {d} x}}+\ldots \right)$
$+s(S-\ldots )+\ldots +Cp+C'\,;\qquad$
équation qui ne s’élève qu’à l’ordre $n-2.$ seulement, et qu’on pourrait abaisser encore si $p$ était nul.
Soit proposé d’assigner la courbe, joignant deux points donnés, dont la révolution autour de l’axe des $x$ engendre une surface minimum ?
La surface de révolution considérée comme indéfinie, étant exprimée par $2\varpi \int y\operatorname {d} x{\sqrt {1+p^{2}}},$ il faudra qu’on ait
$\delta \int y\operatorname {d} x{\sqrt {1+p^{2}}}=0\,;$
on a donc ici
${\begin{aligned}V&=y{\sqrt {1+p^{2}}},\\M&={\frac {\operatorname {d} V}{\operatorname {d} x}}=0,\\\\N&={\frac {\operatorname {d} V}{\operatorname {d} y}}={\sqrt {1+p^{2}}},\\\\P&={\frac {\operatorname {d} V}{\operatorname {d} p}}={\frac {py}{\sqrt {1+p^{2}}}},\end{aligned}}$
et les quantités $Q,R,S,\ldots$ sont toutes nulles. Puis donc que $x$ n’entre pas dans $V,$ on peut faire usage de l’équation (D), qui est simplement ici
$V=pP+c,$
et qui devient, par les substitutions,
$y{\sqrt {1+p^{2}}}={\frac {p^{2}y}{\sqrt {1+p^{2}}}}+C,$
ou encore
${\frac {y}{\sqrt {1+p^{2}}}}=c,$
d’où on tire
$p={\frac {\operatorname {d} y}{\operatorname {d} x}}={\frac {c}{\sqrt {y^{2}-c^{2}}}},$
et par suite
${\frac {\operatorname {d} x}{c}}={\frac {\operatorname {d} y}{\sqrt {y^{2}-c^{2}}}}\,;$
ce qui donne, en intégrant
${\frac {x}{c}}=\operatorname {Log} .\left(y+{\sqrt {y^{2}-c^{2}}}+c'\right).$
Il resterait à déterminer $c$ et $c'$ par la condition que la courbe passe par les deux points donnés ; mais leur détermination ne peut avoir lieu, parce qu’elle entraînerait la résolution d’une équation exponentielle.
Toutefois, on peut reconnaître quelle est la nature de cette courbe. En effet, il est visible que $y$ ne saurait être moindre que $c\,;$ de manière que $c$ est l’ordonnée minimum. Si donc nous prenons cette ordonnée pour axe des $y,$ il faudra qu’en posant $x=0,$ on trouve $y=c,$ d’où $c'=-\operatorname {Log} .c\,;$ donc
$\operatorname {Log} .{\frac {y+{\sqrt {y^{2}-c^{2}}}}{c}}={\frac {x}{c}},$
donc aussi
$\operatorname {Log} .{\frac {c}{y+{\sqrt {y^{2}-c^{2}}}}}\quad$
$\quad \operatorname {Log} .{\frac {y-{\sqrt {y^{2}-c^{2}}}}{c}}=-{\frac {x}{c}}\,:$
cela donne
$y+{\sqrt {y^{2}-c^{2}}}=ce^{+{\frac {x}{c}}},$
$y-{\sqrt {y^{2}-c^{2}}}=ce^{-{\frac {x}{c}}}\,;$
d’où en ajoutant,
$2y=c\left(e^{+{\frac {x}{c}}}+e^{-{\frac {x}{c}}}\right)\,;$
équation que l’on reconnaît pour celle d’une chaînette dont l’axe principal est celui autour duquel la révolution s’exécute.
Soit proposé, en second lieu, de trouver la courbe pour laquelle l’aire comprise entre un arc, les deux normales extrêmes et l’arc correspondant de la développée soit un maximum ?
Soit $u$ l’aire dont il s’agit, et $r$ le rayon de courbure ; on aura
$\operatorname {d} u={\frac {1}{2}}r{\sqrt {\operatorname {d} x^{2}+\operatorname {d} y^{2}}}={\frac {\left(1+p^{2}\right)^{\frac {3}{2}}{\sqrt {\operatorname {d} x^{2}+\operatorname {d} y^{2}}}}{2q}}={\frac {\left(1+p^{2}\right)^{2}\operatorname {d} x}{2q}},$
d’où
$u=\int {\frac {\left(1+p^{2}\right)^{2}\operatorname {d} x}{2q}},$
de sorte qu’il faudra poser
$\delta \int {\frac {\left(1+p^{2}\right)^{2}\operatorname {d} x}{2q}}=0,$
on aura donc ici
${\begin{aligned}V&={\frac {(1+p^{2})^{2}}{2q}},\\M&={\frac {\operatorname {d} V}{\operatorname {d} x}}=0,\\\\N&={\frac {\operatorname {d} V}{\operatorname {d} y}}=0,\\\\P&={\frac {\operatorname {d} V}{\operatorname {d} p}}={\frac {2p(1+p^{2})}{q}},\\\\Q&=-{\frac {(1+p^{2})^{2}}{q}},\end{aligned}}$
$R,S\ldots$ seront nuls ; et comme ni $x$ ni $y$ n’entrent dans $V,$ nous pourrons faire usage de l’équation (E) qui se réduira à
$V=qQ+Cp+C',$
et deviendra, par les substitutions,
${\frac {(1+p^{2})^{2}}{q}}=Cp+C'.$
Il faudra intégrer cette dernière équation, et déterminer les quatre constantes arbitraires, en exprimant que la courbe passe par les deux points limites, et qu’en ces deux points $Q'$ et $Q''$ sont nuls.
Mais si l’on veut seulement savoir quelle est la nature de la courbe, on remarquera que, $r$ désignant toujours le rayon de courbure, son équation différentielle revient à
$r{\sqrt {1+p^{2}}}=Cp+C',$
ce qui donne successivement
${\begin{aligned}r&={\frac {Cp+C'}{\sqrt {1+p^{2}}}},\\\\\operatorname {d} r&={\frac {(C-C'p)\operatorname {d} p}{(1+p^{2})^{\frac {3}{2}}}},\\\\r\operatorname {d} v={\frac {(1+p^{2})^{2}}{q}}.{\frac {(C-C'p)\operatorname {d} p}{(1+p^{2})^{\frac {3}{2}}}}&={\frac {(C-C'p)\operatorname {d} p}{q}}=(C-C'p)\operatorname {d} x,\end{aligned}}$
c’est-à-dire,
$r\operatorname {d} r=C\operatorname {d} x-C'\operatorname {d} y$
d’où, en intégrant et transformant les constantes,
$r^{2}=Cx+C'y+C''.\qquad$
On peut faire disparaître ces constantes, en choisissant les axes d’une manière convenable. On peut d’abord porter l’origine au point pour lequel le rayon de courbure est nul. En désignant alors par $\alpha ,\beta$ les coordonnées de ce point, on aura
$C\alpha +C'\beta +C''=0.\qquad$
Nous pourrons ensuite faire tourner le système des axes autour de la nouvelle origine $(\alpha ,\beta ).$ Les formules générales de cette transformation étant, comme l’on sait
$x=\alpha +t\operatorname {Cos} .m-u\operatorname {Sin} .m,\qquad y=\beta +t\operatorname {Sin} .m+u\operatorname {Cos} .m$
on trouvera, en substituant dans (1) et ayant égard à la relation (2)
$r^{2}=t(C\operatorname {Cos} .m+C'\operatorname {Sin} .m)+u(C'\operatorname {Cos} .m-C\operatorname {Sin} .m).$
L’angle $m$ étant arbitraire, on pourra en disposer de manière que
$C\operatorname {Cos} .m+C'\operatorname {Sin} .m=0,$
il en résultera
$\operatorname {Cos} .m={\frac {C'}{\sqrt {C^{2}+C'^{2}}}},\qquad \operatorname {Sin} .m=-{\frac {C}{\sqrt {C^{2}+C'^{2}}}}\,;$
on aura conséquemment
$r^{2}=u(C'\operatorname {Cos} .m-C\operatorname {Sin} .m)=u{\frac {C^{2}+C'^{2}}{\sqrt {C^{2}+C'^{2}}}}=u{\sqrt {C^{2}+C'^{2}}}\,;$
ou, en transformant les constantes
$r^{2}=au\,;$
Soient faits
${\frac {\operatorname {d} u}{\operatorname {d} t}}=p_{1},\quad {\frac {\operatorname {d} ^{2}u}{\operatorname {d} t^{2}}}={\frac {\operatorname {d} p_{1}}{\operatorname {d} t}}=q_{1},\quad$
$r={\frac {\left(1+p_{1}^{2}\right)^{\frac {2}{3}}}{q_{1}}}\,;$
il en résultera
${\frac {\frac {\operatorname {d} p_{1}}{\operatorname {d} t}}{\left(1+p_{1}^{2}\right)^{\frac {3}{2}}}}={\frac {1}{r}}={\frac {1}{2{\sqrt {a}}}}.{\frac {1}{\sqrt {u}}},\quad$
$\quad {\frac {2p_{1}\operatorname {d} p_{1}}{\left(1+p_{1}^{2}\right)^{\frac {3}{2}}}}={\frac {1}{\sqrt {a}}}.{\frac {p_{1}\operatorname {d} t}{\sqrt {u}}}={\frac {1}{\sqrt {a}}}.{\frac {\operatorname {d} u}{\sqrt {u}}}$
d’où, en intégrant,
$-{\frac {1}{\sqrt {1+p_{1}^{2}}}}={\sqrt {\frac {u}{a}}},\quad 1+p_{1}^{2}={\frac {a}{u}},\quad p_{1}={\frac {\operatorname {d} u}{\operatorname {d} t}}=\pm {\sqrt {{\frac {a}{u}}-1}}\,;$
équation qui appartient à une cycloïde engendrée par un cercle dont le diamètre est égal à $a.$
Souvent la courbe qui doit donner le maximum ou le minimum est assujettie à certaines conditions de constance, comme d’être d’une longueur donnée, de produire une aire de révolution de grandeur donnée, etc. On exprime ces conditions en écrivant que les variations des expressions des fonctions qui doivent demeurer constantes sont nulles. Ainsi, soient les conditions
$\int V'\operatorname {d} x=l',\quad \int V''\operatorname {d} x=l'',\quad \int V'''\operatorname {d} x=l''',\ldots$
en nombre quelconque, on écrira
$\delta \int V'\operatorname {d} x=0,\quad \delta \int V''\operatorname {d} x=0,\quad \delta \int V'''\operatorname {d} x=0,\ldots$
Soit, du reste,
$\delta \int V\operatorname {d} x=0$
la condition du maximum ou du minimum.
Il est clair que, si une fonction a un maximum ou un minimum, elle en aura encore un si on lui ajoute des quantités constantes ; puis donc que $\int V'\operatorname {d} x,\ \int V''\operatorname {d} x,\ \int V'''\operatorname {d} x,\ldots$ sont constantes, il en sera de même de $c'\int V'\operatorname {d} x,\ c''\int V''\operatorname {d} x,\ c'''\int V'''\operatorname {d} x,\ldots ,$ $c',c'',c''',\ldots$ représentant des constantes indéterminées ; donc on devra avoir
$\delta \left(\int V\operatorname {d} x+c'\int V'\operatorname {d} x+c''\int V''\operatorname {d} x+\ldots \right)$
$=\delta \int \left(V+c'V'+c''V''+\ldots \right)\operatorname {d} x=0.$
Nous pourrons traiter $\delta \int \left(V+c'V'+c''V''+\ldots \right)\operatorname {d} x,$ comme nous avons traité $\delta \int V\operatorname {d} x,$ dans ce qui précède, et nous trouverons ainsi le maximum ou le minimum cherché, avec les conditions auxquelles il doit satisfaire. Mais, une fois parvenus à l’équation indéfinie, il faudra la déterminer de manière à satisfaire non seulement aux conditions que nous avons indiquées, en traitant du maximum et du minimum de $\int V\operatorname {d} x,$ mais encore aux conditions $\int V'\operatorname {d} x=l',$ $\int V''\operatorname {d} x=l'',\ldots \,;$ et c’est à quoi serviront les constantes $c',c'',c''',\ldots .$ Il était nécessaire de les introduire, pour n’avoir pas moins de quantités à déterminer que de conditions à remplir. Du reste, nous aurons ici autant de conditions qu’il en faudra pour en assigner les valeurs.
Soit, par exemple, à trouver, parmi toutes les courbes isopérimètres, celle qui à la plus grande aire.
On a, en général, pour la courbe de plus grande aire.
$\delta \int y\operatorname {d} x=0.$
Si $l$ est la longueur donnée, on aura, en outre
$\int \operatorname {d} x{\sqrt {1+p^{2}}}=l\,;$
nous poserons donc
$\delta \int \left(y+c{\sqrt {1+p^{2}}}\right)\operatorname {d} x=0\,;$
nous aurons ainsi
${\begin{aligned}V&=y+c{\sqrt {1+p^{2}}},\\M&={\frac {\operatorname {d} V}{\operatorname {d} x}}=0,\\\\N&={\frac {\operatorname {d} V}{\operatorname {d} y}}=1,\\\\P&={\frac {\operatorname {d} V}{\operatorname {d} p}}={\frac {cp}{\sqrt {1+p^{2}}}},\end{aligned}}$
et tous les coefficiens ultérieurs seront nuls. L’équation (D) qui est celle qui convient au cas actuel, et qui se réduit à
$V=pP+C,$
deviendra, en substituant
$y+c{\sqrt {1+p^{2}}}={\frac {cp^{2}}{\sqrt {1+p^{2}}}}+C\,;$
ou bien
$y+{\frac {c}{\sqrt {1+p^{2}}}}=C,\qquad$
$\quad p={\frac {\operatorname {d} y}{\operatorname {d} x}}={\frac {\sqrt {c^{2}-(y-C)^{2}}}{y-C}}\,;$
et, par suite
${\frac {(y-C)\operatorname {d} y}{\sqrt {c^{2}-(y-C)^{2}}}}=\operatorname {d} x\,;$
ce qui donne, en intégrant
${\sqrt {c^{2}-(y-C)^{2}}}=-x+C',$
c’est-à-dire
$(y-C)^{2}+(x-C')^{2}=c^{2}\,;$
équation du cercle. Il restera à déterminer les constantes $c,C,C',$ par les conditions que l’arc se termine à deux points donnés, et que sa longueur entre ces deux points soit égale à $l.$
Il est à remarquer qu’ayant dû poser ici
$\delta \int y\operatorname {d} x=0\,;\qquad \delta \int \operatorname {d} x{\sqrt {1+p^{2}}}=0\,;$
ce sont aussi les mêmes équations que nous aurions posées s’il avait été question d’assigner, parmi les arcs de courbes qui comprennent une aire donnée, celui de moindre longueur ; de sorte que le cercle doit résoudre, à la fois, les deux problèmes. Il est clair, en effet, que si, à longueur égale, le cercle embrasse la plus grande aire, il s’ensuit qu’il contient le plus de surface sous le moindre développement possible ; et, par suite, à aire égale, il aura la moindre longueur.
Proposons-nous encore d’assigner, parmi les courbes isopérimètres, celle qui engendre l’aire de révolution minimum ?
Nous aurons ici l’équation
$\delta \int (y+c){\sqrt {1+p^{2}}}\operatorname {d} x=0.$
En faisant
$y+c=y',\quad$
${\frac {\operatorname {d} y}{\operatorname {d} x}}={\frac {\operatorname {d} y'}{\operatorname {d} x}},\qquad$
$p=p',$
notre équation deviendra
$\delta \int y'{\sqrt {1+p'^{2}}}.\operatorname {d} x=0\,;$
équation que nous avons déjà traitée ci-dessus, et qui nous a conduit à l’équation de la chainette ; c’est donc encore cette courbe qui résout ce dernier problème. Il arrive seulement ici que son axe principal se trouve à la distance $c$ de l’axe de révolution. Du reste, la constante $c,$ ainsi que les autres qu’introduira l’intégration de l’équation indéfinie, provenant de $\delta \int y'{\sqrt {1+p'^{2}}}.\operatorname {d} x=0,$ se détermineront comme dans l’exemple qui précède celui-ci.
Si, dans l’équation
$\delta \int V\operatorname {d} x=0,$
$V$ contenait, outre $x$ et $y,$ une troisième variable, $z$ liée aux deux premières par une équation de relation, on l’exprimerait au moyen de celles-ci, et on ramènerait $V$ à ne plus contenir ainsi que $x,y,p,q,\ldots .$
Qu’il soit question, par exemple, de trouver la courbe la moins longue que l’on puisse tracer sur un cylindre, entre deux points de sa surface. En prenant l’axe du cylindre pour axe des $y$ et désignant par $a$ son rayon, l’équation de sa surface sera
$z={\sqrt {a^{2}-x^{2}}}\quad$
$\quad {\frac {\operatorname {d} y}{\operatorname {d} x}}=-{\frac {x}{\sqrt {a^{2}-x^{2}}}}.$
On aura ensuite, pour la longueur de l’arc
$s=\int {\sqrt {\operatorname {d} x^{2}+\operatorname {d} y^{2}+\operatorname {d} z^{2}}}=\int \operatorname {d} x{\sqrt {1+p^{2}+{\frac {a^{2}}{a^{2}-x^{2}}}}}$
$=\int \operatorname {d} x{\sqrt {{\frac {a^{2}}{a^{2}-x^{2}}}+p^{2}}}\,;$
de sorte que la condition du minimum sera
$\delta \int \operatorname {d} x{\sqrt {{\frac {a^{2}}{a^{2}-x^{2}}}+p^{2}}}=0.$
On aura donc ici
$V={\sqrt {{\frac {a^{2}}{a^{2}-x^{2}}}+p^{2}}},\qquad P={\frac {p}{\sqrt {{\frac {a^{2}}{a^{2}-x^{2}}}+p^{2}}}}.$
Comme $y$ manque dans $V,$ nous pourrons prendre pour équation indéfinie l’équation (C) qui, dans le cas actuel, est simplement $P=C\,;$ c’est-à-dire, en substituant,
${\frac {p}{\sqrt {{\frac {a^{2}}{a^{2}-x^{2}}}+p^{2}}}}=C,\qquad$
$\quad p={\frac {\operatorname {d} y}{\operatorname {d} x}}={\frac {C}{\sqrt {1-C^{2}}}}.{\frac {a}{\sqrt {a^{2}-x^{2}}}}$
ce qui donne, en intégrant,
$y={\frac {C}{\sqrt {1-C^{2}}}}.a\operatorname {Arc} \left(\operatorname {Tang} .={\frac {x}{a}}\right)+C'.$
Si l’on prend l’une des extrémités de l’arc dans le plan des $yz,$ à une distance $b$ de l’axe des $z,$ on devra avoir, en même temps, $x=0,\ y=b\,;$ en mettant donc simplement $C$ pour ${\frac {C}{\sqrt {1-C^{2}}}},$ on aura
$y=C.\operatorname {Arc} \left(\operatorname {Tang} .={\frac {x}{a}}\right)+b\,;$
dans laquelle $C$ se déterminera en fixant la situation de l’autre extrémité de l’arc demandé.
L’équation que nous venons d’obtenir est celle d’une hélice ; et c’est ce qu’il était facile de prévoir à l’avance. La plus courte ligne qu’on puisse tracer sur un cylindre, entre deux points de sa surface, doit être telle, en effet, qu’elle demeure encore la plus courte, en développant la surface de ce cylindre sur un plan ; il faut donc que, par l’effet du développement de cette surface, elle devienne une ligne droite, ce qui est la propriété caractéristique de l’hélice.
Nous avons vu, dans ce qui précède, comment, quand on fixe les deux limites entre lesquelles le maximum ou le minimum doit avoir lieu, on peut déterminer les $2n$ constantes qu’introduit, en général, l’intégration de l’équation indéfinie
$0=N-{\frac {\operatorname {d} P}{\operatorname {d} x}}+{\frac {\operatorname {d} ^{2}Q}{\operatorname {d} x^{2}}}-{\frac {\operatorname {d} ^{3}R}{\operatorname {d} x^{3}}}+\ldots$
On pourrait ne pas donner proprement les deux limites, mais seulement les assujettir à la condition de se trouver sur deux courbes données. Il y aura alors $2n+4$ inconnues à déterminer, savoir, les $2n$ constantes $c,c',c'',\ldots ,$ introduites par l’intégration de l’équation indéfinie, et les quatre coordonnées $x',y',x'',y'',$ des deux extrémités de la courbe cherchée. Or ici les $\delta x',\delta y',$ $\delta x'',\delta y'',$ ne sont plus nulles ; mais, puisqu’elles ont lieu le long des courbes données, ce sont les différentielles des coordonnées de ces courbes. En prenant donc les différentielles de leurs équations, nous pourrons exprimer $\delta y'$ et $\delta y''$ respectivement, en fonction de $\delta x'$ et $\delta x''.$ En faisant la substitution de leurs valeurs dans la première partie du développement de $\delta \int V\operatorname {d} x,$ il s’y trouvera $2n$ variations $\delta y',\delta y'',\delta p',\delta p'',\ldots$ dont nous pourrons séparément égaler les coefficiens à zéro. Ensuite, $x',y',x'',y'',$ devront satisfaire aux deux courbes données, ce qui fournira deux relations ; ils devront enfin satisfaire à la courbe trouvée, ce qui en fournira deux autres ; de sorte qu’on aura en tout $2n+4$ conditions pour déterminer les $2n+4$ inconnues.
Supposons, par exemple, que l’on demande une courbe, qui se terminant à deux hyperboles équilatères
$y'={\frac {a^{2}}{x'}},\qquad y''=-{\frac {b^{2}}{x''}},$
produise une aire de révolution maximum autour de l’axe des $x.$ On posera toujours
$\delta \int y\operatorname {d} x={\sqrt {1+p^{2}}}=0\,;$
ce qui donnera, en opérant comme ci-dessus
$y={\frac {1}{2c'}}\left(e^{+{\frac {x}{c}}}+c^{2}e^{-{\frac {x}{c}}}\right).$
Pour déterminer les deux constantes $c,c',$ ainsi que les points extrêmes de la courbe, c’est-à-dire, pour déterminer les six inconnues $c,c',x',y',x'',y'',$ nous écrirons les premiers termes du développement de $\delta \int y\operatorname {d} x={\sqrt {1+p^{2}}},$ qui sont ici
${\begin{aligned}&{\frac {y'}{\sqrt {1+p'^{2}}}}\delta x'+{\frac {p'y'}{\sqrt {1+p'^{2}}}}\delta y'\\\\-&{\frac {y''}{\sqrt {1+p''^{2}}}}\delta x''-{\frac {p''y''}{\sqrt {1+p''^{2}}}}\delta y''.\end{aligned}}$
Or, la différentiation des équations des deux hyperboles donne
$\delta y'=-{\frac {a^{2}}{x'^{2}}}\delta x',\qquad \delta y''=-{\frac {a^{2}}{x''^{2}}}\delta x'',$
en substituant donc, il viendra
${\frac {y'}{\sqrt {1+p'^{2}}}}\left(1-{\frac {a^{2}p'}{x'^{2}}}\right)\delta x'-{\frac {y''}{\sqrt {1+p''^{2}}}}\left(1+{\frac {b^{2}p''}{x''^{2}}}\right)\delta x''.$
Il faudra donc poser
$1-{\frac {a^{2}p'}{x'^{2}}}=0,\qquad 1+{\frac {b^{2}p''}{x''^{2}}}=0\,;$
et comme d’ailleurs
$-{\frac {a^{2}}{x'^{2}}}={\frac {\operatorname {d} y'}{\operatorname {d} x'}},\qquad {\frac {b^{2}}{x''^{2}}}={\frac {\operatorname {d} y''}{\operatorname {d} x''}}\,;$
cela reviendra à
$1+p'{\frac {\operatorname {d} y'}{\operatorname {d} x'}}=0,\qquad 1+p''{\frac {\operatorname {d} y''}{\operatorname {d} x''}}=0\,;$
équations qui expriment qu’aux points limites la chainette génératrice doit être normale aux deux hyperboles données.
En remplaçant, dans ces équations $p'$ et $p''$ par leurs valeurs, déduites de l’équation de la chainette ; savoir
$p'={\frac {1}{2cc'}}\left(e^{+{\frac {x'}{c}}}-c^{2}e^{-{\frac {x'}{c}}}\right),\qquad p''={\frac {1}{2cc'}}\left(e^{+{\frac {x''}{c}}}-c^{2}e^{-{\frac {x''}{c}}}\right),$
elles deviennent
$2cc'x'^{2}-a^{2}\left(e^{+{\frac {x'}{c}}}-c^{2}e^{-{\frac {x'}{c}}}\right),\qquad 2cc'x''^{2}+b\left(e^{+{\frac {x''}{c}}}-c^{2}e^{-{\frac {x''}{c}}}\right)=0.$
On a d’ailleurs
$y'=+{\frac {a^{2}}{x'}},\qquad y''=-{\frac {b^{2}}{x''}},$
$y'={\frac {1}{2c'}}\left(e^{+{\frac {x'}{c}}}+c^{2}e^{-{\frac {x'}{c}}}\right),\qquad y''={\frac {1}{2c'}}\left(e^{+{\frac {x''}{c}}}+c^{2}e^{-{\frac {x''}{c}}}\right)\,;$
on se trouve donc avoir six équations, pour déterminer $c,c',x',y',x'',y''.$
Si la courbe qui doit satisfaire au maximum ou au minimum était assujettie à des conditions de constance, comme d’avoir une aire, une longueur, etc., déterminée, il serait facile d’opérer comme il a déjà été dit.
Si cependant ces conditions, au lieu d’être exprimées par des intégrales, l’étaient par des équations algébriques ; si, par exemple, au lieu de donner $\int V\operatorname {d} x=l,$ on donnait $\varphi (x',y',x'',y'',l)=0\,;$ on remplacerait d’abord $y'$ et $y''$ par leurs valeurs en $x'$ et $x'',$ ce qui donnerait $\psi (x',x'',l)=0,$ d’où $\delta x''=\delta x'\Phi (x',l).$ Remplaçant donc $\delta x''$ par cette valeur, on n’aurait plus que le coefficient de $\delta x'$ à égaler à zéro, ce qui ne laisserait plus que $2n+3$ relations ; mais celle qui aurait disparu se trouverait alors suppléée par l’équation donnée $\varphi (x',y',x'',y'',l)=0.$ C’est ainsi, par exemple, qu’on en agirait, si la chaînette qui doit se terminer aux deux hyperboles était assujettie à avoir une corde d’une longueur donnée. L’équation de condition serait ainsi $(x'-x'')^{2}+(y'-y'')^{2}=l^{2}.$
* ↑ Cette exposition a été rédigée par l’auteur, pour son cours d’analise à l’École polytechnique. On peut aussi consulter, sur le même sujet, un mémoire inséré au commencement du XIII.e volume du présent recueil. J. D. G.
* ↑ Il est bien essentiel de remarquer que les courbes pour lesquelles il s’agit de prouver que l’une peut être considérée comme la variation de l’autre, sont des courbes individuelles, c’est-à-dire, à paramètres invariables ou numériques. Il est bien clair alors que les coefficiens du développement, par la suite de Maclaurin, doivent être regardés comme indépendans les uns des autres, puisque ces $b,b',b'',\ldots B,B',B'',\ldots$ sont numériques. Si l’on pouvait donner diverses valeurs aux paramètres des deux courbes que nous considérons, les coefficiens du développement dépendraient les uns des autres, et on ne pourrait plus regarder les coefficiens du second développement comme résultant de la variation de ceux du premier.
* ↑ Il est essentiel de remarquer que $a',b',\ldots ,$ étant de nouvelles valeurs de $a,b,\ldots ,$ sont constantes, dans les différentiations relatives à $\operatorname {d} \,;$ et que $\delta a=a'-a,$$\delta b=b'-b,\ldots \,;$ parce que, dans les différentiations relatives à $\delta ,$ les quantités $a,b,\ldots$ sont, dans la fonction, $y=\operatorname {f} (x,a,b,\ldots )$ et dans ses dérivées, des variables indépendantes, pour lesquelles les différentielles marquées par $\delta a,\delta b,\ldots$ sont la même chose que les différences entières représentées par $a'-a,\ b'-b,\ldots$
* ↑ On ne donne que $2n-2$ coefficiens différentiels $p',p'',q',q'',\ldots$ s’il doit y en avoir $n$ dans l’équation indéfinie ; parce que les conditions de $y=y'$ et $y=y'',$ lorsque $x=x'$ et $x=x''$ fournissent des relations restantes, le cas où l’on donnerait $2n$ coefficiens différentiels ou plus n’a pas encore été traité.
* ↑ M. Ampère est ici le premier qui ait expliqué nettement pourquoi, dans tous les cas, la quantité sous le signe $\int$ que contient l’équation commune au maximum et au minimum doit être séparément nulle. C’est le seul point sur lequel notre mémoire du tom. XIII, déjà cité, nous avait semblé vulnérable.J. D. G.
* ↑ Voyez Annales, tom. 1, pag. 58. J. D. G.