Mémoires extraits des recueils de l’Académie royale de Berlin/Sur les intégrales particulières des équations différentielles

Joseph-Louis Lagrange

Sur les intégrales particulières des équations différentielles

Mémoires extraits des recueils de l’Académie royale de Berlin, Gauthier-Villars, 1868, Œuvres de Lagrange. Tome IV (p. 5-108).

SUR LES
INTÉGRALES PARTICULIÈRES
DES
ÉQUATIONS DIFFÉRENTIELLES^[1].

---

(Nouveaux Mémoires de l’Académie royale des Sciences et Belles-Lettres
de Berlin, année 1774.)

Dans un Mémoire de feu M. Clairaut, imprimé parmi ceux de l’Académie des Sciences de Paris pour l’année 1734, on trouve cette remarque singulière, qu’il y a des équations différentielles qu’on peut intégrer par la différentiation, et que les intégrales trouvées de la sorte ne sont jamais comprises dans les intégrales complètes que donnent les règles ordinaires de l’intégration, quoique d’ailleurs ces mêmes intégrales satisfassent aux équations différentielles proposées et résolvent très-bien les Problèmes géométriques qui conduisent à ces équations (voyez les Mémoires de 1734, pages 209 et suivantes).

M. Euler a mis ensuite ces deux espèces de paradoxes dans un plus grand jour, et il les a confirmés par différents Exemples tirés de la Géométrie c’est le sujet d’un Mémoire donné à cette Académie et imprimé dans le volume de 1756 sous le titre d’Exposition de quelques paradoxes dans le Calcul intégral. Ce grand Géomètre avait aussi déjà remarqué, dans sa Mécanique, qu’il y a souvent des solutions particulières qui échappent à la solution générale, et il avait même donné une formule pour trouver ces solutions particulières dans un grand nombre de cas (voyez Mechanica, tome II, Articles 268, 303, 335) ; mais ni M. Clairaut ni M. Euler n’avaient encore cherché les moyens de reconnaître à priori si une équation finie qui satisfait à une équation différentielle donnée est comprise ou non dans l’intégrale complète de cette équation différentielle, sans connaître cette intégrale.

Ce Problème, qui est d’une grande importance dans la Théorie du Calcul intégral, a depuis été résolu par M. Euler dans le premier volume de son Calcul intégral. M. d’Alembert s’en est occupé aussi et en a rendu la solution plus rigoureuse et plus générale dans les Mémoire de l’Académie des Sciences de Paris, année 1769 (voyez pages 84 et suivantes). On trouve de plus quelques principes généraux sur le même sujet dans les Ouvrages de M. le Marquis de Condorcet (voyez son Calcul intégral, page 67, les Mémoires de Turin, tome IV, pages 7 et suivantes). Enfin je viens de lire un Mémoire sur les solutions particulières des équations différentielles, que M. de Laplace a donné depuis peu à l’Académie des Sciences, et qui doit paraître dans le volume de 1772, mais dont l’Auteur a bien voulu m’envoyer d’avance un exemplaire imprimé. Dans ce Mémoire, M. de Laplace perfectionne et étend plus loin la théorie déjà connue des solutions particulières, et, ce que personne n’avait encore fait, il donne des méthodes pour trouver directement toutes les solutions particulières qui peuvent satisfaire à une équation différentielle donnée, et qui ne seraient point comprises dans la solution générale de cette équation.

Cette lecture a réveillé d’anciennes idées que j’avais sur la même matière et a occasionné les recherches suivantes, dans lesquelles je me flatte de pouvoir présenter aux Géomètres une Théorie nouvelle et complète sur le point d’Analyse dont il s’agit.

Article I^er. — Des intégrales particulières des équations différentielles du premier ordre à deux variables, et de la manière de les déduire des intégrales complètes.

1. J’entends, en général, par intégrale complète d’une équation différentielle du premier ordre, une équation finie qui satisfait à cette équation différentielle et qui renferme une constante arbitraire ; si l’on donne à cette constante une valeur déterminée, l’intégrale devient alors incomplète, parce qu’elle ne renferme plus de constante arbitraire ; mais elle sera toujours comprise dans l’intégrale complète. Les intégrales particulières dont nous allons traiter ici sont celles qui, ne renfermant point de constante arbitraire, ne sont pas non plus comprises dans l’intégrale complète, et par conséquent échappent à la méthode ordinaire d’intégration. Par exemple, l’intégrale complète de l’équation

$xdx+ydy=dy{\sqrt {x^{2}+y^{2}-b^{2}}}$

est

${\sqrt {x^{2}+y^{2}-b^{2}}}=y+a,$

ou bien

$x^{2}-2ay-a^{2}-b^{2}=0,$

$a$ étant une constante arbitraire. Si l’on donnait à $a$ une valeur déterminée quelconque, comme si l’on faisait $a=3b,$ on aurait

$x^{2}-6by-10b^{2}=0,$

qui serait une intégrale incomplète ; mais l’équation précédente, malgré la constante arbitraire $a,$ n’est pas la seule équation finie qui satisfasse à l’équation différentielle proposée ; car il est aisé de voir que l’équation

$x^{2}+y^{2}-b^{2}=0$

y satisfait aussi, équation qu’on voit bien n’être pas comprise dans celle-là, puisque l’une est à un cercle dont le rayon est $b,$ et l’autre est à une parabole ayant $2a$ pour paramètre.

L’équation

$x^{2}+y^{2}-b^{2}=0$

sera donc une intégrale particulière de l’équation différentielle

$xdx+ydy=dy{\sqrt {x^{2}+y^{2}-b^{2}}}\,;$

on trouvera de même que l’équation

$x^{2}+y^{2}-b^{2}=0$

est l’intégrale particulière de l’équation différentielle

$ydx-xdy=b{\sqrt {dx^{2}+dy^{2}}},$

dont l’intégrale complète est

$y-ax-b{\sqrt {1+a^{2}}},$

$a$ étant la constante arbitraire. Il en est de même d’une infinité d’autres équations différentielles qui admettent des intégrales particulières, lesquelles ne sauraient être comprises dans les intégrales complètes de ces équations.

2. Après nous être assurés à posteriori de l’existence des intégrales particulières, cherchons, maintenant à priori, et d’après les seuls principes du Calcul intégral, quelle est l’origine de ces sortes d’intégrales.

Pour considérer les choses d’une manière générale, soit

$\mathrm {Z} =0$

une équation différentielle quelconque, $\mathrm {Z}$ étant une fonction de $x,y$ et de ${\frac {dy}{dx}}\,;$ supposons que l’intégrale complète de cette équation soit

$\mathrm {V} =0,$

$\mathrm {V}$ étant une fonction de $x,y$ et d’une constante arbitraire $a$ qui n’entre point dans la fonction $\mathrm {Z} ,$ et voyons comment cette équation $\mathrm {V} =0$ satisfait à l’équation différentielle $\mathrm {Z} =0.$

L’équation $\mathrm {V} =0,$ étant différentiée, donne celle-ci

$dy=pdx,\quad {\text{ou bien}}\quad {\frac {dy}{dx}}=p,$

où $p$ est une fonction finie de $x,y$ et $a\,;$ puis donc que ces deux équations

$\mathrm {V} =0\quad {\text{et}}\quad {\frac {dy}{dx}}-p=0$

ont lieu en même temps, on peut en éliminer la quantité $a,$ et l’on aura une équation entre $x,y$ et ${\frac {dy}{dx}}$ où $a$ n’entrera plus, et qui aura donc lieu en même temps que l’équation finie $\mathrm {V} =0\,;$ ce sera donc nécessairement l’équation $\mathrm {Z} =0.$

D’où l’on peut conclure que, si $\mathrm {V} =0$ est l’intégrale complète d’une équation différentielle du premier ordre $\mathrm {Z} =0,$ celle-ci ne peut être autre chose que le résultat de l’élimination de la constante arbitraire $a,$ à l’aide des deux équations $\mathrm {V} =0$ et $dy=pdx.$

Ainsi, l’équation finie

$x^{2}-2ay-a^{2}-b^{2}=0$

donne par la différentiation

${\frac {dy}{dx}}={\frac {x}{a}}\,;$

éliminant $a,$ on aura

${\frac {dy}{dx}}={\frac {x}{-y+{\sqrt {x^{2}+y^{2}-b^{2}}}}}\,;$

ce qui se réduit à

$xdx+ydy=dy{\sqrt {x^{2}+y^{2}-b^{2}}}\,;$

par conséquent, l’équation ci-dessus sera l’intégrale complète de cette équation différentielle, $a$ étant la constante arbitraire (1).

3. Maintenant, puisque l’équation différentielle $\mathrm {Z} =0$ résulte des deux équations $\mathrm {V} =0$ et ${\frac {dy}{dx}}=p,$ en éliminant la constante $a,$ il est visible que l’on aura toujours la même équation $\mathrm {Z} =0,$ quelle que soit la valeur de la quantité $a$ qu’on doit éliminer ; ainsi, quand même cette quantité ne serait pas constante, comme on l’a supposé jusqu’ici, il est clair que l’équation finie $\mathrm {V} =0$ satisfera toujours à l’équation différentielle $\mathrm {Z} =0,$ pourvu que, par la différentiation de l’équation $\mathrm {V} =0,$ on ait également, dans le cas de $a$ variable, $dy=pdx.$ Or, en faisant varier dans l’équation $\mathrm {V} =0$ les quantités $x,y$ et $a$ à la fois, il est clair qu’on aura cette équation différentielle

$dy=pdx+qda,$

$p$ et $q$ étant des fonctions de $x,y$ et $a\,;$ et l’on voit que, pour que cette équation se réduise à $ay=pdx,$ comme dans le cas de $a$ constante, il n’y a qu’à supposer la quantité $q$ égale à zéro. Faisant donc

$q=0,$

on aura une équation par laquelle on pourra déterminer la valeur de $a$ en $x$ et $y,$ et cette valeur de $a$ étant ensuite substituée dans l’équation finie $\mathrm {V} =0,$ cette équation satisfera encore à l’équation différentielle $\mathrm {Z} =0,$ et en sera une intégrale particulière.

4. Comme en faisant varier à la fois les quantités $x,y$ et $a$ dans l’équation $\mathrm {V} =0$ on a

$dy=pdx+qda,$

on aura, en faisant varier seulement $y$ et $a,$

$dy=qda,$

et par conséquent

$q={\frac {dy}{da}}\,;$

donc si $\mathrm {V} =0$ est l’intégrale complète d’une équation différentielle du premier ordre $\mathrm {Z} =0,$ et que $a$ soit la constante arbitraire introduite par l’intégration ; qu’on fasse varier dans l’équation $\mathrm {V} =0,$ $y$ et $a,$ et qu’on suppose ensuite ${\frac {dy}{da}}=0\,;$ qu’on détermine $a$ par le moyen de cette équa-

tion, et qu’on substitue sa valeur dans l’intégrale complète $\mathrm {V} =0,$ ou bien, ce qui revient au même, qu’on élimine $a$ par le moyen des deux équations ${\frac {dy}{da}}=0$ et $\mathrm {V} =0,$ on aura une intégrale particulière de la même équation différentielle $\mathrm {Z} =0.$ Si l’équation ${\frac {dy}{da}}=0$ renferme la quantité $a$ mêlée avec les variables $x$ et $y,$ alors la valeur de $a$ tirée de cette équation sera une fonction des mêmes variables ; par conséquent l’intégrale particulière, qui résultera de la substitution de cette valeur de $a$ dans l’équation $\mathrm {V} =0,$ sera nécessairement différente de l’intégrale complète $\mathrm {V} =0,$ dans laquelle $a$ est supposée constante.

Mais il peut arriver, ou que l’équation ${\frac {dy}{da}}=0$ ne renferme que la quantité $a$ avec des constantes, mais sans $x$ ni $y\,;$ ou que cette équation ne renferme que $x$ et $y$ sans $a.$ Dans le premier cas, la valeur de $a$ sera nécessairement constante ; ainsi il n’y aura point alors d’intégrale particulière proprement dite.

Dans le second cas, l’équation ${\frac {dy}{da}}=0$ sera elle-même une intégrale de la proposée ; mais, pour pouvoir juger si c’est une intégrale particulière ou non, il faudra combiner cette équation avec l’équation $\mathrm {V} =0$ en éliminant $x$ ou $y,$ et voir si la résultante donne a variable ou constante. S’il arrivait que la valeur de $a$ demeurât indéterminée ou $={\frac {0}{0}},$ ce serait une marque que l’équation ${\frac {dy}{da}}=0$ est un facteur de l’équation $\mathrm {V} =0,$ indépendant de la constante arbitraire $a,$ et par conséquent étranger à l’équation différentielle $\mathrm {Z} =0.$

Au reste, si l’équation ${\frac {dy}{da}}=0$ avait des facteurs, il faudrait appliquer à chacune des équations qui en résulteraient ce que nous venons de dire en général sur l’équation ${\frac {dy}{da}}=0.$

5. Si, au lieu de faire varier $y$ et $a$ dans l’équation $\mathrm {V} =0,$ on y fait varier $x$ et $a,$ et qu’on suppose ${\frac {dx}{da}}=0,$ cette équation, traitée comme l’équation ${\frac {dy}{da}}=0,$ servira aussi à déterminer les intégrales particulières de la même équation différentielle $\mathrm {Z} =0\,;$ ce qui est aisé à démontrer par les mêmes principes que nous avons établis ci-dessus. L’équation ${\frac {dx}{da}}=0$ donnera le plus souvent les mêmes résultats que l’équation ${\frac {dy}{da}}=0\,;$ mais il y a des cas où ces équations donnent des résultats différents il faudra donc avoir égard à ces deux équations, pour pouvoir trouver toutes les intégrales particulières de l’équation $\mathrm {Z} =0,$ et il est facile de démontrer qu’il n’y a pas d’autres combinaisons possibles qui puissent fournir des intégrales de cette espèce non comprises dans l’intégrale complète $\mathrm {V} =0.$

6. Pour éclaircir la théorie précédente par quelques exemples, je prends d’abord l’équation différentielle

$xdx+ydy=dy{\sqrt {x^{2}+y^{2}-b^{2}}},$

dont nous avons déjà vu que l’intégrale complète est

$x^{2}-2ay-a^{2}-b^{2}=0,$

où $a$ est la constante arbitraire. Faisant donc varier d’abord $y$ et $a,$ on aura

${\frac {dy}{da}}=-{\frac {a+y}{a}},$

et faisant varier $x$ et $a,$ on a

${\frac {dx}{da}}={\frac {a+y}{x}}\,;$

les deux équations.

${\frac {dy}{da}}=0,\quad {\frac {dx}{da}}=0$

donnent également

$a+y=0,\quad {\text{d’où}}\quad a=-y\,;$

ce qui étant substitué dans l’intégrale complète, on a

$x^{2}+y^{2}-b^{2}=0,$

intégrale particulière de l’équation différentielle proposée, et la seule qui ait lieu.

Soit ensuite l’équation différentielle

$ydx-xdy=b{\sqrt {dx^{2}+dy^{2}}},$

dont l’intégrale complète est.

$y-ax-b{\sqrt {1+a^{2}}}=0.$

En faisant varier $y$ et $a,$ on en tire

${\frac {dy}{da}}=x+{\frac {ba}{\sqrt {1+a^{2}}}}\,;$

donc

$x+{\frac {ba}{\sqrt {1+a^{2}}}}=0,\quad {\text{ou bien}}\quad x=-{\frac {ba}{\sqrt {1+a^{2}}}}\,;$

et combinant cette équation avec la précédente, on aura

$y={\frac {b}{\sqrt {1+a^{2}}}}\,;$

donc

$x^{2}+y^{2}=b^{2},$

intégrale particulière de la proposée. Si l’on faisait varier $x$ et $a,$ on aurait

${\frac {dx}{da}}=-{\frac {x}{a}}-{\frac {b}{\sqrt {1+a^{2}}}}\,;$

de sorte que l’équation ${\frac {dx}{da}}=0$ donnera le même resultat que l’équation ${\frac {dy}{da}}=0\,;$ et qu’ainsi il n’y aura d’autre intégrale particulière possible que celle qu’on a trouvée.

Soit de plus l’équation différentielle

$ydx-xdy+bdy={\sqrt {c^{2}\left(dx^{2}+dy^{2}\right)-b^{2}dx^{2}}},$

dont on trouve que l’intégrale complète est

$y-a(x-b)={\sqrt {c^{2}\left(1+a^{2}\right)-b^{2}}},$

$a$ étant la constante arbitraire.

Faisant done varier $y$ et $a,$ on aura

${\frac {dy}{da}}=x-b+{\frac {c^{2}a}{\sqrt {c^{2}\left(1+a^{2}\right)-b^{2}}}}.$

donc

$x-b+{\frac {c^{2}a}{\sqrt {c^{2}\left(1+a^{2}\right)-b^{2}}}}=0,$

par conséquent

$x-b=-{\frac {c^{2}a}{\sqrt {c^{2}\left(1+a^{2}\right)-b^{2}}}}=0,$

et combinant cette équation avec la précédente, on aura

$y={\frac {c^{2}-b^{2}}{\sqrt {c^{2}\left(1+a^{2}\right)-b^{2}}}},$

de sorte qu’on aura, en éliminant $a,$

${\frac {(x-b)^{2}}{c^{2}}}+{\frac {y^{2}}{c^{2}-b^{2}}}=1,$

intégrale particulière de la proposée, et qui n’est point comprise dans l’intégrale complète.

Si l’on fait varier $x$ et $a,$ on aura

${\frac {dx}{da}}=-{\frac {x-b}{a}}-{\frac {c^{2}}{\sqrt {c^{2}\left(1+a^{2}\right)-b^{2}}}},$

et l’équation ${\frac {dx}{da}}$ redonnera le même résultat que nous venons de trouver d’après l’équation ${\frac {dy}{da}}=0.$ Ainsi la proposée n’admet point d’autre intégrale particulière que la précédente. Voyez au reste, sur l’intégration de ces deux équations, les n^os 17 et 19 ci-après.

Considérions enfin l’équation différentielle séparée

${\frac {dx}{\sqrt {\mathrm {X} }}}={\frac {dy}{\sqrt {\mathrm {Y} }}},$

dans laquelle

${\begin{aligned}\mathrm {X} &=\mathrm {A} +\mathrm {B} x+\mathrm {C} x^{2}+\mathrm {D} x^{3}+\mathrm {E} x^{4},\\\mathrm {Y} &=\mathrm {A} +\mathrm {B} y+\mathrm {C} y^{2}+\mathrm {D} y^{3}+\mathrm {E} y^{4}\,;\end{aligned}}$

l’intégrale complète de cette équation est, comme j’ai fait voir ailleurs^[2],

${\sqrt {\mathrm {X} }}+{\sqrt {\mathrm {Y} }}=(x-y){\sqrt {a+\mathrm {U} }},$

où $a$ est la constante arbitraire, et

$\mathrm {U} =\mathrm {D} (x+y)+\mathrm {E} (x+y)^{2}.$

Faisant d’abord varier $y$ et $a,$ et ensuite $x$ et $a$ à la fois, et supposant, pour plus de simplicité,

${\frac {d\mathrm {X} }{dx}}=\mathrm {X} ',\quad {\frac {d\mathrm {Y} }{dy}}=\mathrm {Y} ',\quad {\frac {d\mathrm {U} }{dx}}={\frac {d\mathrm {U} }{dy}}=\mathrm {U} ',$

on aura

${\begin{aligned}{\frac {dy}{da}}=&{\frac {(x-y){\sqrt {\mathrm {Y} }}}{\mathrm {Y} '{\sqrt {a+\mathrm {U} }}-\left[\mathrm {U} '(x-y)+2(a+\mathrm {U} )\right]{\sqrt {\mathrm {Y} }}}}\\{\frac {dx}{da}}=&{\frac {(x-y){\sqrt {\mathrm {X} }}}{\mathrm {X} '{\sqrt {a+\mathrm {U} }}-\left[\mathrm {U} '(x-y)+2(a+\mathrm {U} )\right]{\sqrt {\mathrm {X} }}}}.\end{aligned}}$

Les équations ${\frac {dy}{da}}=0$ et ${\frac {dx}{da}}=0$ donnent d’abord la même équation

$x-y=0,$

laquelle, ne contenant point la quantité $a,$ peut être ou n’être pas une intégrale particulière. Pour pouvoir en juger, je reprends l’intégrale complète, et j’en tire

$a={\frac {\mathrm {\left({\sqrt {X}}+{\sqrt {Y}}\right)} ^{2}}{(x-y)^{2}}}-\mathrm {U} ,$

où l’on voit qu’en faisant $x=y,$ $a$ devient $=\infty$ par conséquent, l’équation $x-y=0$ n’est point une intégrale particulière, mais un cas de l’intégrale complète.

Rejetant donc le facteur $x-y,$ l’équation ${\frac {dy}{da}}=0$ donne encore $\mathrm {Y} =0\,;$ mais, comme le dénominateur de la fraction se réduit alors à $\mathrm {Y} '{\sqrt {a+\mathrm {U} }},$ qui devient aussi nul lorsque $\mathrm {Y} '$ est en même temps égal à zéro, il s’ensuit que l’équation $\mathrm {Y} =0$ peut être une intégrale particulière, pourvu que $\mathrm {Y} '$ ne soit pas à la fois égal à zéro. De même, l’équation ${\frac {dx}{da}}=0$ donnera $\mathrm {X} =0,$ pourvu que $\mathrm {X} '$ ne soit pas en même temps égal à zéro. Or il est clair, par l’expression de $a$ trouvée ci-dessus, que la valeur de $a$ ne devient point constante par la supposition de $\mathrm {Y} =0,$ ni par celle de $\mathrm {X} =0.$ Donc on peut conclure que les équations $\mathrm {Y} =0$ et $\mathrm {X} =0$ seront deux intégrales particulières de la proposée, pourvu que l’on n’ait pas en même temps $\mathrm {Y} '=0$ et $\mathrm {X} '=0\,;$ ainsi donc les intégrales particulières de l’équation dont il s’agit seront toutes comprises sous cette forme

$x=u\quad {\text{ou}}\quad y=u,$

en prenant pour $u$ une des racines simples quelconques de l’équation

$\mathrm {A} +\mathrm {B} u+\mathrm {C} u^{2}+\mathrm {D} u^{3}+\mathrm {E} u^{4}=0.$

Si l’équation proposée était

${\frac {dx}{\sqrt {\mathrm {X} }}}+{\frac {dy}{\sqrt {\mathrm {Y} }}}=0,$

l’intégrale complète serait

${\sqrt {\mathrm {X} }}-{\sqrt {\mathrm {Y} }}=(x-y){\sqrt {a+\mathrm {U} }}\,;$

d’où l’on tirerait les mêmes valeurs de ${\frac {dy}{da}}$ et ${\frac {dx}{da}}$ que ci-dessus, à l’exception que dans la première le radical ${\sqrt {\mathrm {Y} }}$ y serait avec un signe différent.

On aurait donc d’abord l’équation

$x-y=0\,;$

mais comme la valeur de $a,$ qui dans ce cas est

$\left({\frac {{\sqrt {\mathrm {X} }}-{\sqrt {\mathrm {Y} }}}{x-y}}\right)^{2}-\mathrm {U} ,$

devient, par la supposition de $x=y,$ égale à ${\frac {0}{0}},$ il s’ensuit que l’équation $x-y=0$ doit être rejetée comme étrangère à l’équation différentielle

${\frac {dx}{\sqrt {\mathrm {X} }}}+{\frac {dy}{\sqrt {\mathrm {Y} }}}=0,$

quoiqu’elle soit contenue (4) dans l’intégrale

${\sqrt {\mathrm {X} }}-{\sqrt {\mathrm {Y} }}=(x-y){\sqrt {a+\mathrm {U} }}.$

Ensuite on trouvera, comme ous haut, les intégrales particulières

$x=u\quad {\text{ou}}\quad y=u,$

$u$ étant une des racines simples de l’équation

$\mathrm {A} +\mathrm {B} u+\mathrm {C} u^{2}+\mathrm {D} u^{3}+\mathrm {E} u^{4}=0.$

7. On voit donc par ce que nous venons de démontrer comment, lorsqu’on a trouvé l’intégrale complète d’une équation différentielle du premier ordre, on en peut aisément déduire les intégrales particulières qui satisfont à la même équation ; on voit aussi que si ces intégrales particulières ne sont pas comprises dans l’intégrale complète, ce n’est nullement une imperfection du Calcul intégral, comme on pourrait le croire, faute de donner à ce Calcul toute la généralité dont il est susceptible. Ainsi l’on doit regarder la théorie que nous venons de donner, moins comme une exception que comme un supplément nécessaire à la règle générale du Calcul intégral.

Article II. — De l’étendue des intégrales particulières des équations différentielles du premier ordre, et de la manière de trouver ces intégrales sans connaître les intégrales complète.

8. L’équation finie

$\mathrm {V} =0,$

dans laquelle $\mathrm {V}$ est une fonction des variables $x,y$ et d’une arbitraire $a,$ donne par la différentiation, et en faisant varier à la fois $x,y$ et, $a,$

$dy=pdx+qda,$

$p,q$ étant des fonctions finies de $x,y$ et $a\,;$ différentiant $p$ et substituant pour $dy$ sa valeur, on aura

$dp=p'dx+q'da,$

et différentiant de même $p',$ on aura

$dp'=p''dx+q''da,$

et ainsi de suite. Maintenant si l’on regarde $a$ comme constante, on a pour le premier ordre

$dy=pdx\,;$

et toute équation différentielle du premier ordre, telle que

$\mathrm {Z} =0,$

à laquelle satisfera l’équation finie $\mathrm {V} =0,$ la constante $a$ demeurant arbitraire, sera nécessairement produite par la combinaison deux équations

$V=0,\quad {\frac {dy}{dx}}=p,$

de manière que $a$ s’évanouisse.

En regardant toujours $a$ comme constante, on a

${\frac {dy}{dx}}=p\quad {\text{et}}\quad dp=p'dx\,;$

donc (en prenant $dx$ pour constante)

${\frac {d^{2}y}{dx^{2}}}=p'\,;$

et toute équation différentielle du second ordre, telle que

$\mathrm {Z} '=0,$

à laquelle satisfera l’équation finie

$\mathrm {V} =0,$

la constante $a$ demeurant arbitraire, sera nécessairement formée par la combinaison des équations

$\mathrm {V} =0,\quad {\frac {dy}{dx}}=p,\quad {\frac {d^{2}y}{dx^{2}}}=p',$

en sorte que $a$ disparaisse.

En continuant ainsi, dans l’hypothèse de $a$ constante, on aura

$dp'=p''dx,$

par conséquent

${\frac {d^{3}y}{dx^{3}}}=p''\,;$

et toute équation différentielle du troisième ordre, telle que

$\mathrm {Z} ''=0,$

à laquelle satisfera l’équation finie

$\mathrm {V} =0,$

$a$ demeurant arbitraire, sera formée par la combinaison des équations

$\mathrm {V} =0,\quad {\frac {dy}{dx}}=p,\quad {\frac {d^{2}y}{dx^{2}}}=p',\quad {\frac {d^{3}y}{dx^{3}}}=p'',$

de manière que $a$ disparaisse ; et ainsi de suite.

9. Voyons maintenant dans quels cas l’équation $\mathrm {V} =0$ pourra satisfaire aux mêmes équations

$\mathrm {Z=0,\quad Z'=0,\quad Z''} =0,\ldots ,$

en supposant que $a$ soit une quantité variable.

Et d’abord il est clair que cela aura lieu pour l’équation du premier ordre $\mathrm {Z} =0,$ si $q=0\,;$ parce qu’alors on aura également ${\frac {dy}{dx}}=p,$ comme dans le cas de $a$ constante. De là naissent les intégrales particulières, ainsi que nous l’avons vu dans l’Article précédent.

Pour l’équation du second ordre $\mathrm {Z} '=0,$ il faudra que l’on ait de plus $q'=0,$ afin que l’on ait aussi ${\frac {d^{2}y}{dx^{2}}}=p',$ comme dans l’hypothèse de $a$ constante.

De même pour l’équation du troisième ordre $\mathrm {Z} ''=0,$ il faudra que l’on ait encore $q''=0$ pour que la valeur de ${\frac {d^{3}y}{dx^{3}}}$ soit également $=p''\,;$ et ainsi de suite.

Donc, en général, l’équation finie $\mathrm {V} =0$ sera une intégrale particulière de l’équation du premier ordre $\mathrm {Z} =0,$ si $a$ est une quantité telle, que l’on ait $q=0$. Elle sera une intégrale particulière de l’équation du second ordre $\mathrm {Z} '=0,$ si l’on a à la fois $q=0$ et $q'=0.$ Elle sera une intégrale particulière de l’équation du troisième ordre $\mathrm {Z} ''=0,$ si l’on a en même temps $q=0,\ q'=0,\ q''=0\,;$ et ainsi de suite.

10. En regardant $y$ comme une fonction de $x$ et $a$ donnée par l’équation $\mathrm {V} =0,$ on a, suivant la notation reçue (8).

$p={\frac {dy}{dx}},\quad q={\frac {dy}{da}},$

ensuite

$p'={\frac {dp}{dx}},\quad q'={\frac {dp}{da}},\quad p''={\frac {dp'}{dx}},\quad q''={\frac {dp'}{da}},\ldots ,$

donc

$q={\frac {dy}{da}},\quad q'={\frac {d^{2}y}{dxda}},\quad q''={\frac {d^{3}y}{dx^{2}da}},\ldots .$

Ainsi l’on aura pour les équations différentielles du premier ordre la condition

${\frac {dy}{da}}=0,$

pour celles du second ordre les deux conditions

${\frac {dy}{da}}=0,\quad {\frac {d^{2}y}{dxda}}=0,$

pour celles du troisième ordre les trois conditions

${\frac {dy}{da}}=0,\quad {\frac {d^{2}y}{dxda}}=0,\quad {\frac {d^{3}y}{dx^{2}da}}=0\,;$

et ainsi de suite.

Et comme on peut échangera $y$ en $x,$ en regardant $x$ comme une fonction de $y$ et $a,$ on aura de même ($dy$ étant pris pour constante)

${\frac {dx}{da}}=0$

pour le premier ordre,

${\frac {dx}{da}}=0,\quad {\frac {d^{2}x}{dyda}}=0$

pour le second ordre,

${\frac {dx}{da}}=0,\quad {\frac {d^{2}x}{dyda}}=0,\quad {\frac {d^{3}x}{dy^{2}da}}=0$

pour le troisième ordre ; et ainsi de suite.

11. De là il s’ensuit que si $\mathrm {Z} '=0$ est une équation différentielle du second ordre dont l’intégrale aux premières différences soit l’équation différentielle $\mathrm {Z} =0,$ l’intégrale particulière de cette dernière, trouvée d’après la condition de

${\frac {dy}{da}}=0\quad {\text{ou}}\quad {\frac {dx}{da}}=0,$

ne satisfera pas, en général, à l’équation proposée $\mathrm {Z} '=0,$ à moins que l’on n’ait à la fois

${\frac {dy}{da}}=0\quad {\frac {d^{2}y}{dxda}}=0,\quad {\text{ou}}\quad {\frac {dx}{da}}=0,\quad {\frac {d^{2}x}{dyda}}=0.$

De même, si $\mathrm {Z} ''=0$ est une équation différentielles du troisième ordre, dont l’intégrale aux premières différences soit l’équation $\mathrm {Z} =0,$ l’intégrale particulière de cette dernière équation, déduite de la condition

${\frac {dy}{da}}=0\quad {\text{ou}}\quad {\frac {dx}{da}}=0,$

ne satisfera pas à la proposée $\mathrm {Z} ''=0,$ à moins que l’on n’ait à la fois

${\frac {dy}{da}}=0,\quad {\frac {d^{2}y}{dxda}}=0,\quad {\frac {d^{3}y}{dx^{2}da}}=0,$

ou bien

${\frac {dx}{da}}=0,\quad {\frac {d^{2}x}{dyda}}=0,\quad {\frac {d^{3}x}{dy^{2}da}}=0\,;$

et ainsi de suite.

Donc, en général, si dans la solution d’un Problème on a été conduit directement à une équation différentielle d’un ordre supérieur au premier, et qu’on ait déjà ramené cette équation au premier ordre à l’aide d’une ou de plusieurs intégrations, l’intégrale particulière de cette équation du premier ordre ne résoudra pas le Problème, à moins que toutes les conditions relatives à l’ordre de l’équation différentielle primitive ne se trouvent remplies.

Mais, si l’équation primitive du Problème n’est que du premier ordre, l’intégrale particulière de cette équation résoudra la question tout aussi bien que l’intégrale complète.

12. L’équation

$xd^{2}y-dydx=0$

a pour intégrale complète du premier ordre

$xdx+ydy-dy{\sqrt {x^{2}+y^{2}-b^{2}}}=0\,;$

et celle-ci a pour intégrale complète finie

$x^{2}-2ay-a^{2}-b^{2}=0,$

d’où l’on tire

${\frac {dy}{da}}=-{\frac {a+y}{a}},$

ce qui, étant fait égal à zéro, donne, pour l’intégrale particulière,

$a+y=0,\quad \quad {\text{d’où}}\quad a=-y,$

et par conséquent

$x^{2}+y^{2}-b^{2}=0,$

comme on l’a déjà vu (6).

Maintenant, pour que cette intégrale particulière satisfasse aussi à l’équation différentio-diflérentielle, il faudra que l’on ait en même temps

${\frac {dy}{da}}=0,\quad {\text{et}}\quad {\frac {d^{2}y}{dxda}}=0\,;$

différentiant donc la valeur trouvée de ${\frac {dy}{da}},$ on aura

${\frac {d^{2}y}{dxda}}=-{\frac {1}{a}}{\frac {dy}{dx}}\,;$

mais de l’équation

$x^{2}-2ay-a^{2}-b^{2}=0$

on tire

${\frac {dy}{dx}}={\frac {x}{a}}\,;\quad {\text{donc}}\quad {\frac {d^{2}y}{dxda}}=-{\frac {x}{a^{2}}},$

ce qui ne peut pas être égal à zéro, en général. D’où il faut conclure que, quoique l’équation

$x^{2}+y^{2}-b^{2}=0$

satisfasse à l’équation différentielle du premier ordre

$xdx+ydy-dy{\sqrt {x^{2}+y^{2}-b^{2}}}=0,$

elle ne satisfera cependant pas à l’équation différentio-différentielle

$xd^{2}y-dydx=0,$

qui en est dérivée. En effet, on trouve

${\frac {dy}{dx}}={\frac {-x}{\sqrt {b^{2}-x^{2}}}},\quad {\frac {d^{2}y}{dx^{2}}}={\frac {-b^{2}}{\left(b^{2}-x^{2}\right)^{\frac {3}{2}}}},$

ce qui, comme l’on voit, ne satisfait pas à l’équation dont il s’agit.

Prenons maintenant l’équation différentio-difféuentielle

$\left(d^{2}y\right)^{2}+4xd^{2}ydx^{2}-4dydx^{3}=0,$

dont l’intégrale du premier ordre est

$x^{2}dx+dy=\left(b^{3}-x^{3}-3y\right)^{\frac {2}{3}}dx,$

laquelle a pour intégrale complète

$y-a^{2}x+ax^{2}+{\frac {a^{3}-b^{3}}{3}}=0.$

On aura donc

${\frac {dy}{da}}=2ax-x^{2}-a^{2}=-(a-x)^{2}\,;$

ce qui, étant fait égal à zéro, donne $a-x=0,$ par conséquent $a=x,$ et

$y+{\frac {a^{3}-b^{3}}{3}}=0$

pour l’intégrale particulière. Pour que cette intégrale satisfasse donc aussi à l’équation différentio-différentielle, il faudra que ${\frac {d^{2}y}{dxda}}=2(a-x)$ soit

nul en même temps ; ce qui est en effet ; donc, etc : On peut s’assurer à posteriori que l’équation

$y={\frac {b^{3}-x^{3}}{3}}$

satisfait à la proposée : car on a

${\frac {dy}{dx}}=-x^{2}\quad {\text{et}}\quad {\frac {d^{2}y}{dx^{2}}}=-2x\,;$

ce qui étant substitué dans les termes

$\left(d^{2}y\right)^{2}+4xd^{2}ydx^{2}-4dydx^{3},$

tout se détruit de soi-même.

13. S’il arrivait que l’on eût en même temps

${\frac {dy}{da}}=0,\quad {\frac {d^{2}y}{dxda}}=0,\quad {\frac {d^{3}y}{dx^{2}da}}=0,\ldots ,$

ou bien

${\frac {dx}{da}}=0,\quad {\frac {d^{2}x}{dyda}}=0,\quad {\frac {d^{3}x}{dy^{2}da}}=0,\ldots ,$

et ainsi de suite à l’infini ; alors l’intégrale particulière satisferait nonseulement à l’équation différentielle du premier ordre, mais aussi à toutes les équations différentielles des ordres ultérieurs qui en seraient dérivées. Cette intégrale aurait donc les mêmes propriétés que l’intégrale complète ; et nous allons prouver qu’elle sera alors effectivement comprise dans celle-ci ; de sorte qu’elle cessera d’être une intégrale particulière, et devra être rangée dans la classe des intégrales incomplètes. En effet, puisqu’en regardant la quantité $y$ comme une fonction de $x$ et de $a$ donnée par l’équation $\mathrm {V} =0,$ on a

${\frac {dy}{da}}=0,\quad {\frac {d^{2}y}{dadx}}=0,\quad {\frac {d^{3}y}{dadx^{2}}}=0,\ldots ,$ à l’infini,

il est visible que la quantité ${\frac {dy}{da}}$ ne doit pas contenir $x,$ et ne peut être,

par conséquent, qu’une fonction de $a$ mêlée avec des constantes. Ainsi, l’équation ${\frac {dy}{da}}=0$ donnera $a$ égal à une constante ; donc, etc. Ce sera la même chose si, en regardant $x$ comme une fonction de $y$ et de $a,$ on a en même temps

${\frac {dx}{da}}=0,\quad {\frac {d^{2}x}{dady}}=0,\quad {\frac {d^{3}x}{dady^{2}}}=0,\ldots ,$ à l’infini.

14. Cette considération nous conduit à une méthode directe pour trouver l’intégrale particulière d’une équation différentielle du premier ordre sans en connaître l’intégrale complète. Soit $\mathrm {Z} =0$ l’équation différentielle du premier ordre dont on cherche l’intégrale particulière, et dont l’intégrale complète est $\mathrm {V} =0,$ $\mathrm {Z}$ étant une fonction de $x,y,{\frac {dy}{dx}},$ et $\mathrm {V}$ une fonction de $x,y$ et de l’arbitraire $a.$ Puisque l’équation $\mathrm {Z} =0$ est indépendante de la quantité $a,$ il s’ensuit qu’on aura également ${\frac {d\mathrm {Z} }{da}}=0,$ en regardant $y$ comme une fonction de $x$ et $a,$ ou $x$ comme une fonction de $y$ et $a,$ donnée par l’équation $\mathrm {V} =0\,;$ et pour avoir la valeur de ${\frac {d\mathrm {Z} }{da}}$ il faudra différentier $\mathrm {Z}$ en faisant varier $y$ seul dans le premier cas, ou $x$ seul dans le second.

Supposons, pour plus de simplicité, que l’équation $\mathrm {Z} =0$ ne renferme point de fonctions transcendantes, et imaginons, ce qui est toujours possible et ne change point la nature de l’équation, qu’elle soit délivrée des fractions et des radicaux, en sorte que $\mathrm {Z}$ soit une fonction entière et rationnelle de $x,y$ et ${\frac {dy}{dx}}\,;$ on aura, en général, par la différentiation,

$d\mathrm {Z} =\mathrm {A} d{\frac {dy}{dx}}+\mathrm {B} dy+\mathrm {C} dx,$

$\mathrm {A,B}$ et $\mathrm {C}$ étant aussi des fonctions rationnelles et entières des mêmes quantités ; donc, en regardant $y$ comme une fonction de $x$ et $a,$ on aura

${\frac {d\mathrm {Z} }{da}}=\mathrm {A} {\frac {d^{2}y}{dxda}}+\mathrm {B} {\frac {dy}{da}}=0\,;$

or dans le cas de l’intégrale particulière on a ${\frac {dy}{da}}=0\,;$ donc, puisque $\mathrm {B}$ ne peut devenir infini étant une fonction sans dénominateur, le terme $\mathrm {B} {\frac {dy}{da}}$ deviendra nul, et il faudra qu’on ait

$\mathrm {A} {\frac {d^{2}y}{dxda}}=0\,;$

donc, si ${\frac {d^{2}y}{dxda}}$ n’est pas nul, il faudra que $\mathrm {A} =0.$

Si ${\frac {d^{2}y}{dxda}}=0$ en même temps que ${\frac {dy}{da}}=0,$ l’équation

$\mathrm {A} {\frac {d^{2}y}{dxda}}+\mathrm {B} {\frac {dy}{da}}=0\,;$

aura lieu d’elle-même ; mais en prenant la différentielle de cette équation, $x$ et $y$ variant à la fois et $a$ demeurant constante, j’aurai

$\mathrm {A} {\frac {d^{3}y}{dx^{2}da}}+\left({\frac {d\mathrm {A} }{dx}}+\mathrm {B} \right){\frac {d^{2}y}{dxda}}+{\frac {d\mathrm {B} }{dx}}{\frac {dy}{da}}=0\,;$

or on a ${\frac {dy}{da}}$ et ${\frac {d^{2}y}{dxda}}$ nuls à la fois par l’hypothèse ; donc, puisque les coefiicients de ces quantités ne sauraient devenir infinis, étant des fonctions sans dénominateur, l’équation précédente se réduira à

$\mathrm {A} {\frac {d^{3}y}{dx^{2}da}}=0,$

laquelle, si ${\frac {d^{3}y}{dx^{2}da}}$ n’est pas nul, donne de nouveau $\mathrm {A} =0.$

Si ${\frac {d^{3}y}{dx^{2}da}}$ est nul aussi, on trouvera, par une nouvelle différentiation, que l’on aura nécessairement $\mathrm {A} =0,$ à moins que ${\frac {d^{4}y}{dx^{3}da}}$ ne soit nul ; et ainsi de suite à l’infini.

Mais nous avons vu ci-dessus que pour que l’équation ${\frac {dy}{da}}=0$ donne une intégrale particulière, il faut que les quantités

${\frac {d^{2}y}{dxda}},\quad {\frac {d^{3}y}{dx^{2}da}},\ldots$ à l’infini

ne soient pas toutes nulles à la fois ; donc il faudra que quelqu’une de ces quantités ne soit pas nulle ; par conséquent il faudra nécessairement qu’on ait

$\mathrm {A} =0.$

Or l’équation différentielle proposée $\mathrm {Z} =0$ donne par la différentiation

$d\mathrm {Z} =\mathrm {A} d{\frac {dy}{dx}}+\mathrm {B} dy+\mathrm {C} dx=0\,;$

donc, puisque dans le cas de l’intégrale particulière $\mathrm {A}$ doit être nul, cette équation se réduira à

$\mathrm {B} dy+\mathrm {C} dx=0\quad {\text{ou bien}}\quad \mathrm {B} {\frac {dy}{dx}}+\mathrm {C} =0,$

laquelle devra s’accorder avec l’équation $\mathrm {A} =0,$ après avoir chassé la valeur de ${\frac {dy}{dx}},$ au moyen de l’équation proposée $\mathrm {Z} =0.$

15. Donc, puisque la valeur de ${\frac {d^{2}y}{dx^{2}}}$ tirée de l’équation différentielle $\mathrm {Z} =0$ au moyen de la différentiation est exprimée, en général, par

$-{\frac {\mathrm {B} {\cfrac {dy}{dx}}+\mathrm {C} }{\mathrm {A} }},$

il s’ensuit de ce que nous venons de démontrer que cette valeur deviendra égale à ${\frac {0}{0}}$ dans le cas de l’intégrale particulière tirée de la condition ${\frac {dy}{da}}=0\,;$ et l’on prouvera de même que la condition ${\frac {dx}{da}}=0$ rendra la valeur de ${\frac {d^{2}x}{dy^{2}}}$ égale à ${\frac {0}{0}},$ en prenant ici $dy$ pour constante au lieu de $dx.$

Et quoique la démonstration précédente soit fondée sur l’hypothèse que l’équation proposée ne renferme aucune fonction transcendante, il n’est cependant pas difficile de se convaincre que la même conclusion aura lieu quelles que soient la nature et la forme de cette équation.

16. En supposant que l’équation différentielle proposée $\mathrm {Z} =0$ donne par la différentiation

$d\mathrm {Z} =\mathrm {A} d{\frac {dy}{dx}}+\mathrm {B} dy+\mathrm {C} dx=0\,;$

s’il arrive que les deux termes $\mathrm {B} dy+\mathrm {C} dx$ se détruisent d’eux-mêmes, on aura

$\mathrm {A} d{\frac {dy}{dx}}=0,\quad$et par conséquent$\quad d{\frac {dy}{dx}}={\frac {0}{\mathrm {A} }}\,;$

de sorte que dans ce cas l’une et l’autre condition ${\frac {d^{2}y}{dx^{2}}}={\frac {0}{0}}$ et ${\frac {d^{2}x}{dy^{2}}}={\frac {0}{0}}$ sera remplie par la condition unique $\mathrm {A} =0$ ; ainsi il n’y aura qu’à éliminer la quantité ${\frac {dy}{dx}}$ au moyen des deux équations $\mathrm {A} =0$ et $\mathrm {Z} =0,$ et l’équation résultante entre $x$ et $y$ sera l’intégrale particulière de la proposée. Quant à l’intégrale complète, elle est facile à déduire de l’équation

$\mathrm {A} d{\frac {dy}{dx}}=0\,;$

car cette équation, lorsque $\mathrm {A}$ n’est point nul, donne

$d{\frac {dy}{dx}}=0,\quad$et par conséquent$\quad {\frac {dy}{dx}}=a,$

$a$ étant une constante arbitraire ; il n’y aura donc qu’à substituer cette valeur de ${\frac {dy}{dx}}$ dans l’équation donnée $\mathrm {Z} =0,$ et l’on aura l’intégrale complète où $a$ sera la constante arbitraire.

Voyons maintenant quels sont les cas où l’on aura

$\mathrm {B} dy+\mathrm {C} dx=0.$

Pour plus de simplicité, je fais ${\frac {dy}{dx}}=p,$ en sorte que $\mathrm {Z}$ soit une fonction de $x,y,p,$ dont la différentielle $d\mathrm {Z} =\mathrm {A} dp+\mathrm {B} dy+\mathrm {C} dx\,;$ puis donc

que $\mathrm {B} dy+\mathrm {C} dx=0,$ on aura

$\mathrm {C} =-\mathrm {B} {\frac {dy}{dx}}=-\mathrm {B} p\,;$

donc

$d\mathrm {Z} =\mathrm {A} dp+\mathrm {B} (dy-pdx)=0\,;$

donc

$dy-pdx+{\frac {\mathrm {A} dp}{\mathrm {B} }}=0,$

et intégrant

$y-px+\int \left(\mathrm {\frac {A}{B}} +x\right)dp=0\,;$

de sorte qu’il faudra que la quantité $\mathrm {\frac {A}{B}} +x$ soit une fonction de $p$ sans $x$ ni $y\,;$ et alors l’équation sera

$y-px+f(p)=0,$

$f(p)$ dénotant une fonction quelconque de $p$ seul.

17. Toute équation donc de la forme

$y-px+f(p)=0,$

$p$ étant ${\frac {dy}{dx}},$ donnera par la différentiation $\left[{\text{en supposant }}f'(p)={\frac {df(p)}{dp}}\right]$ celle-ci

$\left[f'(p)-x\right]dp=0\,;$

en faisant $dp=0,$ on aura

$p=a,$

et

$y-ax+f(a)=0$

sera l’intégrale complète où $a$ est arbitraire ; en faisant

$f'(p)-x=0,$

et éliminant $p$ au moyen de cette équation et de la proposée

$y-px+f(p)=0,$

on aura l’intégrale particulière de cette dernière équation. On voit par là que l’intégrale complète ne donnera jamais autre chose qu’une ligne droite, tandis que l’intégrale particulière donnera toujours une courbe ; nous en donnerons la raison à priori dans l’Article suivant.

Ces sortes d’équations sont celles qui donnent lieu aux paradoxes dont il est question dans les Mémoires de MM. Clairaut et Euler que nous avons cités au commencement de ce Mémoire ; et l’on doit voir maintenant que le vrai dénouement de ces paradoxes tient à la théorie des intégrales particulières que nous venons d’exposer.

18. Reprenons les Exemples que nous avons apportés dans le n^o 6, et voyons si la règle ci-dessus donnera les mêmes intégrales particulières que nous avons trouvées d’après les intégrales complètes.

L’équation

${\frac {dy}{dx}}={\frac {x}{{\sqrt {x^{2}+y^{2}-b^{2}}}-y}}$

donne par la différentiation

${\frac {d^{2}y}{dx^{2}}}={\frac {y^{2}-b^{2}-xy{\cfrac {dy}{dx}}+\left(x{\cfrac {dy}{dx}}-y\right){\sqrt {x^{2}+y^{2}-b^{2}}}}{\left({\sqrt {x^{2}+y^{2}-b^{2}}}-y\right)^{2}{\sqrt {x^{2}+y^{2}-b^{2}}}}}\,;$

faisant cette quantité égale à ${\frac {0}{0}},$ on a les deux équations

$y^{2}-b^{2}-xy{\frac {dy}{dx}}+\left(x{\frac {dy}{dx}}-y\right){\sqrt {x^{2}+y^{2}-b^{2}}}=0,$

$\left({\sqrt {x^{2}+y^{2}-b^{2}}}-y\right)^{2}{\sqrt {x^{2}+y^{2}-b^{2}}}=0\,;$

la seconde donne d’abord,

${\text{ou}}\quad {\sqrt {x^{2}+y^{2}-b^{2}}}-y=0,\quad {\text{ou}}\quad {\sqrt {x^{2}+y^{2}-b^{2}}}=0.$

Dans le premier cas, on aura donc

${\sqrt {x^{2}+y^{2}-b^{2}}}=y,$

et la première équation deviendra par là

$-b^{2}=0,$

ce qui ne donne rien. Dans le second cas, la première équation deviendra

$y^{2}-b^{2}-xy{\frac {dy}{dx}}=0\,;$

mais la proposée donne, en supposant ${\sqrt {x^{2}+y^{2}-b^{2}}}=0,$

${\frac {dy}{dx}}=-{\frac {x}{y}}\,;$

donc l’équation précédente deviendra

$y^{2}+x^{2}-b^{2}=0,$

laquelle s’accorde avec

${\sqrt {x^{2}+y^{2}-b^{2}}}=0\,;$

ainsi cette équation est une intégrale particulière de la proposée.

Si l’on cherche la valeur de ${\frac {d^{2}x}{dy^{2}}}$ en prenant $dy$ pour constante, on aura

${\frac {d^{2}x}{dy^{2}}}={\frac {yx-\left(y^{2}-b^{2}\right){\cfrac {dx}{dy}}+\left(y{\cfrac {dx}{dy}}-x\right){\sqrt {x^{2}+y^{2}-b^{2}}}}{x^{2}{\sqrt {x^{2}+y^{2}-b^{2}}}}},$

d’où l’on tire ces deux équations, en égalant le numérateur et le dénominateur chacun à zéro,

$yx-\left(y^{2}-b^{2}\right){\frac {dx}{dy}}+\left(y{\frac {dx}{dy}}-x\right){\sqrt {x^{2}+y^{2}-b^{2}}}=0,$

$x^{2}{\sqrt {x^{2}+y^{2}-b^{2}}}=0\,;$

la dernière de ces équations donne

${\text{ou}}\quad x=0,\quad {\text{ou}}\quad {\sqrt {x^{2}+y^{2}-b^{2}}}=0\,;$

dans le premier cas, la première équation deviendra

$\left(y{\sqrt {y^{2}-b^{2}}}-y^{2}+b^{2}\right){\frac {dx}{dy}}=0\,;$

mais

${\frac {dx}{dy}}={\frac {{\sqrt {x^{2}+y^{2}-b^{2}}}-y}{x}}=\infty \quad {\text{lorsque}}\quad x=0\,;$

donc $x=0$ n’est pas une intégrale particulière ; reste donc le cas de

${\sqrt {x^{2}+y^{2}-b^{2}}}=0,$

dans lequel la première équation devient

$yx-\left(y^{2}-b^{2}\right){\frac {dx}{dy}}=0\,;$

mais on a, dans ce même cas,

${\frac {dx}{dy}}=-{\frac {y}{x}}\,;$

donc substituant cette valeur et multipliant par ${\frac {x}{y}},$ l’équation précédente deviendra

$x^{2}+y^{2}-b^{2}=0,$

qui s’accorde avec

${\sqrt {x^{2}+y^{2}-b^{2}}}=0,$

en sorte que cette équation sera une intégrale particulière.

Ainsi les deux conditions ${\frac {d^{2}y}{dx^{2}}}={\frac {0}{0}}$ et ${\frac {d^{2}x}{dy^{2}}}={\frac {0}{0}}$ donnent, dans le cas présent, la même intégrale particulière

$x^{2}+y^{2}-b^{2}=0\,;$

ce qui s’accorde avec ce que nous avons trouvé (6), d’où il s’ensuit que

cette équation est l’unique intégrale particulière dont l’équation différentielle proposée soit susceptible.

19. Les deux autres Exemples du n^o 6 appartiennent à la formule

$y-px+f(p)=0$

que nous avons considérée, en général, dans le n^o 17 ci-dessus. En effet, en faisant ${\frac {dy}{dx}}=p,$ les équations différentielles des deux Exemples dont nous parlons se réduisent à ces formes

$y-px-b{\sqrt {1+p^{2}}}=0$

et

$y-px+bp-{\sqrt {c^{2}\left(1+p^{2}\right)-b^{2}}}=0,$

qui sont évidemment des cas particuliers de la forme générale

$y-px+f(p)=0\,;$

or nous avons déjà vu (numéro cité) que cette équation admet toujours une intégrale particulière, laquelle est le résultat de l’élimination de $p$ des deux équations

$y-px+f(p)=0,\quad f'(p)-x=0\,;$

ainsi il ne s’agit que d’examiner si cette intégrale est la même qu’on tirerait de l’intégrale complète par la règle de l’Article I (n^os 4, 5).

L’intégrale dont il s’agit est (17)

$y-ax+f(a)=0,$

d’où l’on tire

${\frac {dy}{da}}=x-f'(x),\quad {\frac {dx}{da}}={\frac {f'(a)-x}{a}}\,;$

ainsi les deux conditions ${\frac {dy}{da}}=0$ et ${\frac {dx}{da}}=0$ donnent également $f'(a)-x=0\,;$

et l’intégrale particulière sera le résultat de l’élimination de $a$ au moyen de cette équation

$f'(a)-x=0,$

et de l’équation

$y-ax+f(a)=0\,;$

or il est visible que ce résultat sera le même que celui de l’élimination de $p$ au moyen des équations

$f'(p)-x=0,\quad y-px+f(p)=0\,;$

donc, etc.

20. Le dernier Exemple du n^o 6 est tiré de l’équation différentielle

${\frac {dx}{\sqrt {\mathrm {X} }}}={\frac {dy}{\sqrt {\mathrm {Y} }}},$

dans laquelle $\mathrm {X}$ est un quinôme en $x,$ et $\mathrm {Y}$ un quinôme semblable en $y\,;$ mais nous supposerons ici que $\mathrm {X}$ soit, en général, un polynôme quelconque en $x,$ et $\mathrm {Y}$ un polynôme quelconque en $y,$ et nous désignerons par $\mathrm {X} '$ et $\mathrm {Y} '$ les valeurs de ${\frac {d\mathrm {X} }{dx}}$ et de ${\frac {d\mathrm {Y} }{dy}},$ lesquelles seront par conséquent aussi des polynômes en $x$ et $y,$ mais d’un degré inférieur d’une unité.

Puis donc que

${\frac {dy}{dx}}=\mathrm {\frac {\sqrt {Y}}{\sqrt {X}}} ,$

on trouvera par la différentiation, après avoir réduit au dénominateur commun,

${\frac {d^{2}y}{dx^{2}}}={\frac {\mathrm {XY} '{\cfrac {dy}{dx}}-\mathrm {YX} '}{2\mathrm {X{\sqrt {XY}}} }},\quad {\frac {d^{2}x}{dy^{2}}}={\frac {\mathrm {YX} '{\cfrac {dx}{dy}}-\mathrm {XY} '}{2\mathrm {Y{\sqrt {XY}}} }},$

en prenant dans la première formule $dx$ constant et dans la seconde $dy$ constant.

Supposons d’abord que les quantités $\mathrm {X}$ et $\mathrm {X} '$ n’aient aucun diviseur commun, non plus que les quantités $\mathrm {Y}$ et $\mathrm {Y} ',$ ce qui arrive lorsque les équations $\mathrm {X} =0$ et $\mathrm {Y} =0$ n’ont point de racines égales ; dans ce cas le numérateur et le dénominateur de l’une et de l’autre quantité ${\frac {d^{2}y}{dx^{2}}}$ et ${\frac {d^{2}x}{dy^{2}}}$ n’auront non plus de diviseur commun.

Donc :

1^o En faisant ${\frac {d^{2}y}{dx^{2}}}={\frac {0}{0}}$ on aura les deux équations

$\mathrm {XY} '{\frac {dy}{dx}}-\mathrm {YX} '=0,\quad \mathrm {X{\sqrt {XY}}} =0,$

dont la seconde donne ou $\mathrm {X} =0$ ou $\mathrm {Y} =0\,;$ mais la première donne, par la substitution de la valeur de ${\frac {dy}{dx}},$

$\mathrm {Y'{\sqrt {XY}}-YX'} =0\,;$

faisant $\mathrm {X} =0,$ cette équation se réduit à $\mathrm {YX} '=0,$ laquelle donnerait $\mathrm {X} '=0,$ ce qui est contre l’hypothèse ; faisant $\mathrm {Y} =0,$ l’équation précédente se trouve remplie d’elle-même ; ainsi $\mathrm {Y} =0$ est une intégrale particulière.

2^o Si l’on fait ${\frac {d^{2}x}{dy^{2}}}={\frac {0}{0}}$ on trouvera, par un raisonnement semblable, l’intégrale particulière $\mathrm {X} =0\,;$ de sorte que ces deux intégrales particulières auront lieu en même temps.

Si l’on suppose que $\mathrm {Y}$ et $\mathrm {Y} '$ aient un diviseur commun, alors il est aisé de voir que ce diviseur disparaîtra entièrement par la division du dénominateur de la quantité ${\frac {d^{2}y}{dx^{2}}}\,;$ par conséquent il ne pourra servir à rendre cette quantité égale à ${\frac {0}{0}}\,;$ il en sera de même relativement à la quantité ${\frac {d^{2}x}{dy^{2}}},$ si $\mathrm {X}$ et $\mathrm {X} '$ ont un diviseur commun.

D’où il faut conclure, en général, que l’équation proposée

${\frac {x}{\sqrt {\mathrm {X} }}}={\frac {y}{\sqrt {\mathrm {Y} }}}$

aura pour intégrales particulières tous les facteurs simples des deux

équations $\mathrm {X} =0$ et $\mathrm {Y} =0\,;$ ce qui s’accorde avec ce que nous avons trouvé dans le n^o 6 pour le cas particulier où $\mathrm {X}$ et $\mathrm {Y}$ étaient des quinômes semblables.

Article III. — Dans lequel on déduit la théorie des intégrales
particulières de la considération des courbes.

21. Soit $\mathrm {V} =0$ l’intégrale complète d’une équation différentielle du premier ordre $\mathrm {Z} =0,$ $\mathrm {Z}$ étant une fonction de $x,y$ et ${\frac {dy}{dx}}\,;$ on sait que $\mathrm {V}$ sera une fonction finie de $x,y$ et d’une constante arbitraire $a\,;$ donc si l’on considère la courbe exprimée par l’équation $\mathrm {V} =0,$ en prenant $x$ et $y$ pour les deux coordonnées, cette courbe exprimera aussi l’équation différentielle $\mathrm {Z} =0,$ quelque valeur qu’on donne à la constante $a\,;$ de sorte qu’en donnant successivement à $a$ toutes les valeurs possibles depuis zéro jusqu’à l’infini positif et négatif, on aura un assemblage d’une infinité de courbes toutes de la même famille, et infiniment peu différentes l’une de l’autre, dont chacune représentera également l’équation différentielle $\mathrm {Z} =0.$

Je dis maintenant que la courbe, qui touchera toutes les courbes dont il s’agit, satisfera aussi à la même équation différentielle $\mathrm {Z} =0.$ Car cette équation détermine la valeur de ${\frac {dy}{dx}}$ par une fonction de $x$ et $y\,;$ par conséquent elle détermine la position de la tangente a chaque point par la position de ce point dans le plan des coordonnées $x$ et $y\,;$ donc toute courbe, qui dans un point quelconque aura la même tangente qu’une des courbes dont nous venons de parler, satisfera aussi nécessairement à l’équation $\mathrm {Z} =0\,;$ or il est visible que la courbe, qui touche toutes les courbes données par l’équation $\mathrm {V} =0,$ en faisant varier le paramètre $a,$ a cette propriété ; donc, etc.

22. Si l’on considère deux points infiniment proches de la courbe touchante, il est facile de concevoir que les deux courbes touchées dans ces points doivent nécessairement se couper dans un point intermédiaire ; par conséquent, en faisant coïncider les deux points d’attouchement, le point d’intersection des deux courbes touchées se confondra avec eux ; d’où il suit que la courbe touchante est formée par l’intersection mutuelle et successive des courbes données par l’équation $\mathrm {V} =0,$ en faisant varier le paramètre $a\,;$ donc cette courbe satisfera à l’équation $\mathrm {Z} =0\,;$ ce qui est d’ailleurs évident, puisque suivant ce point de vue la courbe dont nous parlons n’est composée que de portions infiniment petites des courbes représentées par l’équation $\mathrm {V} =0,$ et dont chacune satisfait à la même équation $\mathrm {Z} =0.$

23. Maintenant, si l’on regarde $y$ comme une fonction de $x$ et de $a$ donnée par l’équation $\mathrm {V} =0,$ il est clair que, pour la même abscisse $x,$ les coordonnées qui répondent à deux courbes infiniment proches seront, en général, $y$ et $y+{\frac {dy}{da}}da\,;$ donc, au point d’intersection de ces deux courbes, on aura ${\frac {dy}{da}}=0\,;$ par conséquent, si l’on élimine $a$ au moyen des deux équations $\mathrm {V} =0$ et ${\frac {dy}{da}}=0,$ on aura l’équation de la courbe formée par les intersections continuelles de toutes les courbes contenues dans l’équation $\mathrm {V} =0,$ laquelle sera aussi la courbe qui touchera toutes ces mêmes courbes.

On prouvera de même, en regardant $x$ comme fonction de $y$ et de $a,$ que la condition ${\frac {dx}{da}}=0,$ combinée avec l’équation $\mathrm {V} =0$ en sorte que $a$ disparaisse, donnera aussi la courbe touchante des mêmes courbes.

D’où et de ce que nous avons démontré plus haut (4 et 5) on doit conclure que l’intégrale particulière d’une équation différentielle du premier ordre est représentée par la courbe qui touche toutes les différentes courbes représentées par l’intégrale complète de cette équation, en faisant varier la constante arbitraire, c’est-à-dire toutes les différentes courbes qui peuvent être représentées à la fois par la même équation différentielle.

Ainsi les caustiques par réflexion et par réfraction ne sont autre chose que les courbes représentées par l’intégrale particulière de l’équation différentielle qui exprime à la fois toutes les lignes droites suivant lesquelles les rayons sont réfléchis ou réfractés.

Et les développées ne sont que les courbes représentées par l’intégrale particulière de l’équation différentielle qui exprime toutes les lignes droites qui coupent la développante à angles droits ; et ainsi du reste.

24. Toute équation différentielle du premier ordre représente donc premièrement une infinité de courbes de la même famille, qui ne diffèrent entre elles que par la valeur de la constante arbitraire, laquelle tient lieu de paramètre ; en second lieu, cette équation représente aussi la courbe qui touche toutes ces mêmes courbes ; en sorte qu’on peut regarder en quelque façon tant les courbes touchées que la courbe touchante comme une seule courbe ayant une infinité de branches liées entre elles par la même équation. Ainsi, à chaque point de la courbe touchante il y aura deux branches qui se rencontrent dans ce point et qui ont une tangente commune ; l’une c’est la courbe touchante même, et l’autre c’est la courbe qu’elle touche dans ce même point ; donc à chaque valeur de ${\frac {dy}{dx}}$ il devra répondre une valeur double de ${\frac {d^{2}y}{dx^{2}}}\,;$ par conséquent l’expression de la quantité ${\frac {d^{2}y}{dx^{2}}},$ tirée de l’équation différentielle proposée au moyen de la différentiation, devra devenir égale à ${\frac {0}{0}}$ pour tous les points de la courbe touchante, par une raison semblable à celle par laquelle on prouve que la valeur de ${\frac {dy}{dx}}$ devient égale à ${\frac {0}{0}}$ dans les points doubles des courbes ; et l’on dira la même chose à l’égard de la quantité ${\frac {d^{2}x}{dy^{2}}},$ en supposant $dy$ constant au lieu de $dx.$ De cette manière on pourra donc déduire de l’équation différentielle même celle de la courbe qui toucherait toutes les différentes courbes représentées par cette équation différentielle ; et comme l’équation de la courbe touchante n’est autre chose que l’intégrale particulière de l’équation différentielle dont nous parlons, ainsi qu’on l’a démontré ci-dessus, il résulte de là la même règle pour trouver ces sortes d’intégrales, que nous avons donnée dans le n^o 15, d’après d’autres principes.

25. Pour jeter un plus grand jour sur la théorie précédente et rendre bien sensible la liaison qu’il y a entre les intégrales complètes et les intégrales particulières, nous allons apporter quelques Exemples tirés de la Géométrie dans lesquels l’application de cette théorie se présente naturellement.

Supposons qu’on demande une courbe telle, que toutes les perpendiculaires menées d’un point donné sur les tangentes de cette courbe soient d’une grandeur donnée.

Il est visible que le cercle résout d’abord la question, pourvu qu’on place le centre dans le point donné et qu’on fasse le rayon égal à la grandeur donnée ; mais comme le Problème conduit naturellement à une équation différentielle du premier ordre, il s’ensuit que la solution complète doit renfermer une constante arbitraire ; par conséquent, puisque le cercle qui résout le Problème est nécessairement donné de grandeur et de position, on ne peut pas avoir par son moyen une solution complète, mais seulement une solution particulière.

Or si l’on considère que la ligne droite satisfait aussi au même Problème, et que pour cela il suffit que la perpendiculaire menée du point donné sur cette ligne soit donnée, on verra qu’il y a une infinité de droites qui résolvent le Problème ; de sorte que ces droites en donnent la véritable solution complète, puisque dans l’équation qui les représente toutes il entre nécessairement une constante arbitraire ; et l’on verra de plus que toutes ces droites sont nécessairement les tangentes du cercle qui donne la solution particulière.

Pour confirmer par le calcul ce que nous venons de trouver synthétiquement, soient $x,y$ les coordonnées de la courbe cherchée, et $t,u$ les coordonnées d’une quelconque de ses tangentes considérée comme une ligne droite ; on aura donc, en général, entre $t$ et $u$ l’équation

$u=mt+n,$

$m$ et $n$ étant constantes pour la même tangente, mais variables d’une tangente à l’autre. Or, comme la droite et la courbe doivent d’abord se rencontrer dans un point, on aura dans ce point $u=y,$ $t=x\,;$ donc

$y=mx+n\,;$

ensuite, comme elles doivent de plus se toucher dans le même point, on aura encore ${\frac {du}{dt}}={\frac {dy}{dx}}\,;$ mais ${\frac {du}{dt}}=m\,;$ donc $m={\frac {dy}{dx}}\,;$ et par conséquent $n=y-x{\frac {dy}{dx}}\,;$ donc l’équation à la tangente sera

$u=t{\frac {dy}{dx}}+y-x{\frac {dy}{dx}}.$

Prenons l’origine des coordonnées pour le point donné, et nommant $\rho$ une ligne menée de ce point à la tangente, on aura

$\rho ^{2}=t^{2}+u^{2}\,;$

donc, pour que cette ligne soit perpendiculaire, il faudra que $d\rho =0,$ ce qui donne

$tdt+udu=0\quad {\text{et}}\quad {\frac {du}{dt}}=-{\frac {t}{u}}\,;$

donc

${\frac {t}{u}}=-{\frac {dy}{dx}},\quad t=-u{\frac {dy}{dx}}\,;$

donc, substituant cette valeur de $t$ dans l’équation ci-dessus, on en tire

$u={\frac {(ydx-xdy)dx}{dx^{2}+dy^{2}}}\,;$

donc

$\rho ^{2}=u^{2}\left(1+{\frac {dy^{2}}{dx^{2}}}\right)={\frac {(ydx-xdy)^{2}}{dx^{2}+dy^{2}}}\,;$

c’est le carré de la perpendiculaire menée du point donné sur la tan-

gente ; nommant donc cette perpendiculaire $b,$ on aura l’équation

$ydx-xdy=b{\sqrt {dx^{2}+dy^{2}}},$

qui servira à résoudre le Problème. Or cette équation a déja été examinée dans le n^o 6, et nous avons vu qu’elle donne l’intégrale complète

$y-ax-b{\sqrt {1+a^{2}}}=0,$

et ensuite l’intégrale particulière

$x^{2}+y^{2}=b^{2}\,;$

ce qui s’accorde avec les résultats trouvés plus haut.

26. Ayant tiré d’un point donné acne perpendiculaire à la tangente d’une courbe, et menant du point où cette perpendiculaire rencontre la tangente à un autre point donné une droite, on demande quelle doit être la nature de la courbe pour que cette droite comprise entre les deux points dont il s’agit soit d’une grandeur donnée.

Par les propriétés connues des sections coniques il est facile de voir que si l’on décrit une section conique qui ait le premier des deux points donnés pour l’un des foyers, l’autre point pour centre, et la grandeur donnée pour demi-axe, cette section conique résoudra le Problème ; mais la section étant entièrement déterminée par ces données, elle ne pourra pas fournir une solution complète du Problème, lequel conduit naturellement à une équation différentielle du premier ordre, et par conséquent indéterminée.

Outre la section conique, on voit aisément qu’il y a une infinité de droites qui peuvent aussi résoudre la question ; car si l’on décrit autour du second point donné un cercle dont le rayon soit égal à la grandeur donnée, toute ligne droite qui coupera ce cercle en un point quelconque, de manière qu’elle fasse un angle droit avec la droite menée de ce point d’intersection au premier point donné, aura évidemment les propriétés requises.

L’équation générale de toutes ces lignes droites renfermant donc une constante arbitraire ; elle donnera nécessairement la solution complète du Problème ; et il est facile de prouver, par les propriétés connues des sections coniques, que toutes ces droites seront tangentes à la section conique que nous avons vu résoudre aussi le Problème ; de sorte que la solution par une section conique ne sera qu’une solution particulière.

En effet, pour réduire le Problème en équation, on remarquera que, si de l’origine des coordonnées on mène une perpendiculaire à une tangente quelconque d’une courbe dont les coordonnées soient $x$ et $y,$ et qu’on nomme $t$ et $u$ les coordonnées qui se rapportent au point de la tangente sur laquelle tombe la perpendiculaire, on remarquera, dis-je, que les formules trouvées dans le n^o 25 ci-dessus donneront

$u={\frac {(ydx-xdy)dx}{dx^{2}+dy^{2}}},\quad t=-u{\frac {dy}{dx}}=-{\frac {(ydx-xdy)dy}{dx^{2}+dy^{2}}}\,;$

maintenant si l’on fait passer l’axe des abscisses par les deux points donnés, qu’on prenne le premier de ces deux points pour l’origine, et qu’on nomme $b$ la distance entre les deux points et $c$ la grandeur donnée, il est aisé de concevoir qu’on aura

$c^{2}=(b-t)^{2}+u^{2}\,;$

donc

$t^{2}+u^{2}-2bt=c^{2}-b^{2},$

et substituant pour $t$ et $u$ les valeurs ci-dessus,

${\frac {(ydx-xdy)^{2}+2b(ydx-xdy)dy}{dx^{2}+dy^{2}}}=c^{2}-b^{2}\,;$

d’où, en multipliant par $dx^{2}+dy^{2},$ et extrayant la racine carrée après avoir ajouté de part et d’autre $b^{2}dy^{2},$ on aura l’équation du Problème

$ydx-xdy+bdy={\sqrt {c^{2}\left(dx^{2}+dy^{2}\right)-b^{2}dx^{2}}}.$

Nous avons déjà traité cette équation dans le n^o 6, et nous avons vu que

son intégrale complète est

$y-a(x-b)={\sqrt {c^{2}\left(1+a^{2}\right)-b^{2}}},$

ce qui donne différentes lignes droites suivant la valeur de la constante arbitraire $a\,;$ nous avons vu ensuite que cette même équation est susceptible d’une intégrale particulière, laquelle est

${\frac {(x-b)^{2}}{c^{2}}}+{\frac {y^{2}}{c^{2}-b^{2}}}=1,$

et représente par conséquent une ellipse dans laquelle les abscisses sont prisses depuis l’un des foyers, et où $b$ est l’excentricité etc le demi grand axe ; de sorte que cette ellipse est la même que celle dont nous avons parlé ci-dessus.

Article IV. — Des intégrales particulières des équations différentielles du second ordre et des ordres plus élevés.

27. Soit

$\mathrm {Z} '=0$

une équation différentielle du second ordre, $\mathrm {Z} '$ étant une fonction de $x,y,$${\frac {dy}{dx}},d{\frac {dy}{dx}}\,;$ et soit

$\mathrm {V} =0$

l’intégrale finie et complète de cette équation : $\mathrm {V}$ sera, dans ce cas, une fonction de $x,y$ et de deux constantes arbitraires $a$ et $b.$ Or, puisque $a$ et $b$ sont arbitraires, on peut supposer, en général, que $b$ soit une fonction quelconque de $a\,;$ alors $\mathrm {V}$ sera une fonction de $x,y$ et $a\,;$ et, de ce que nous avons démontré dans l’Article II, il s’ensuit que l’équation $\mathrm {V} =0$ satisfera également à l’équation $\mathrm {Z} '=0,$ en supposant $a$ variable, pourvu que l’on ait

${\frac {dy}{da}}=0,\quad {\frac {d^{2}y}{dadx}}=0,$

ou

${\frac {dx}{da}}=0,\quad {\frac {d^{2}x}{dady}}=0\,;$

et dans ce cas l’équation $\mathrm {V} =0$ deviendraune intégrale particulière (10).

28. Considérons les deux conditions

${\frac {dy}{da}}=0,\quad {\frac {d^{2}y}{dxda}}=0\,;$

et supposant

$b=f(a),\quad db=f'(a)da,$

il est clair que si l’on regarde $y$ comme une fonction de $a$ et de $b,$ la valeur complète de ${\frac {dy}{da}}$ sera représente par ${\frac {dy}{da}}+{\frac {dy}{db}}f'(a),$ de sorte que les deux conditions dont il s’agit seront exprimées ainsi

${\frac {dy}{da}}+{\frac {dy}{db}}f'(a)=0,\quad {\frac {d^{2}y}{dxda}}+{\frac {d^{2}y}{dxdb}}f'(a)=0.$

Au moyen de ces deux équations on déterminera les valeurs de $a$ et de $f(a)$ ou $b$ en $x$ et $y,$ et on les substituera ensuite dans l’équation $\mathrm {V} =0,$ ou, ce qui revient au même, on éliminera $a$ et $f(a)$ au moyen des trois équations dont il s’agit ; et l’équation résultante sera l’intégrale particulière de l’équation différentio-différentielle $\mathrm {Z} '=0.$

Si l’on remet $db$ à la place de $f'(a)da,$ on aura les deux équations de condition

${\frac {dy}{da}}+{\frac {dy}{db}}{\frac {db}{da}}=0,\quad {\frac {d^{2}y}{dxda}}+{\frac {d^{2}y}{dxdb}}{\frac {db}{da}}=0,$

au moyen desquelles et de l’équation $\mathrm {V} =0$ il faudra éliminer $a$ et $b.$.

En éliminant d’abord les différentiels ${\frac {db}{da}},$ on aura l’équation

${\frac {dy}{da}}{\frac {d^{2}y}{dxdb}}-{\frac {dy}{db}}{\frac {d^{2}y}{dxda}}=0,$

laquelle, étant combinée avec l’équation $\mathrm {V} =0,$ servira à déterminer $a$

et $b$ en $x$ et $y\,;$ et il n’y aura plus qu’à substituer ces valeurs de $a$ et $b$ dans l’équation

${\frac {dy}{da}}da+{\frac {dy}{db}}db=0\,;$

ce qui donnera une équation différentielle du premier ordre en $x$ et $y,$ laquelle sera par conséquent l’intégrale particulière cherchée de l’équation du second ordre

$\mathrm {Z} '=0.$

Mais puisque l’équation $\mathrm {V} =0$ donne par la différentiation, en faisant varier à la fois $x,y,a,b,$

$dy=pdx+{\frac {dy}{da}}da+{\frac {dy}{db}}db=0,$

on aura

${\frac {dy}{da}}da+{\frac {dy}{db}}db=dy-pdx\,;$

par conséquent l’équation

${\frac {dy}{da}}da+{\frac {dy}{db}}db=0$

sera équivalente à celle-ci

$dy-pdx=0\,;$

ainsi il n’y aura qu’à substituer les valeurs de $a$ et de $b$ dans cette dernière équation, ou, ce qui revient au même, éliminer les valeurs de $a$ et de $b$ au moyen des équations

$\mathrm {V} =0,\quad {\frac {dy}{da}}{\frac {d^{2}y}{dxdb}}-{\frac {dy}{db}}{\frac {d^{2}y}{dxda}}=0,\quad dy-pdx=0\,;$

l’équation du premier ordre qui en résultera sera l’intégrale particulière dont il s’agit.

29. De là je conclus, en général, que pour trouver l’intégrale particulière de l’équation différentio-différentielle $\mathrm {Z} '=0,$ dont l’intégrale finie et complète est $\mathrm {V} =0,$ il n’y a qu’à éliminer les quantités $a,b$ et ${\frac {db}{da}}$ au moyen des équations

$\mathrm {V} =0,\quad {\frac {dy}{dx}}-p=0,\quad {\frac {dy}{da}}+{\frac {dy}{db}}{\frac {db}{da}}=0,\quad {\frac {d^{2}y}{dxda}}+{\frac {d^{2}y}{dxdb}}{\frac {db}{da}}=0\,;$

et comme au lieu de regarder $y$ comme une fonction de $x,a,b,$ on peut vice versâ regarder $x$ comme une fonction de $y,a,b,$ on pourra aussi, à la place des deux dernières équations, substituer ces deux-ci

${\frac {dx}{da}}+{\frac {dx}{db}}{\frac {db}{da}}=0,\quad {\frac {d^{2}x}{dyda}}+{\frac {d^{2}x}{dydb}}{\frac {db}{da}}=0.$

30. Soit l’équation du second ordre

$y-x{\frac {dy}{dx}}+{\frac {x^{2}}{2}}{\frac {d^{2}y}{dx^{2}}}=\left({\frac {dy}{dx}}-x{\frac {d^{2}y}{dx^{2}}}\right)^{2}+\left({\frac {d^{2}y}{dx^{2}}}\right)^{2},$

dont l’intégrale finie et complète est

$y={\frac {ax^{2}}{2}}+bx+a^{2}+b^{2},$

$a$ et $b$ étant les deux constantes arbitraires. On tire par la différentiation

${\frac {dy}{dx}}=ax+b,\quad {\frac {dy}{da}}={\frac {x^{2}}{2}}+2a,\quad {\frac {dy}{db}}=x+2b,\quad {\frac {d^{2}y}{dxda}}=x,\quad {\frac {d^{2}y}{dxdb}}=1\,;$

on aura donc ces quatre équations

${\begin{array}{ll}y={\dfrac {ax^{2}}{2}}+bx+a^{2}+b^{2},&{\dfrac {dy}{dx}}=ax+b,\\{\dfrac {x^{2}}{2}}+2a+(x+2b){\dfrac {db}{da}}=0,&x+{\dfrac {db}{da}}=0\,;\end{array}}$

au moyen desquelles, éliminant les quantités $a,b,{\frac {db}{da}},$ on aura pour résultante l’intégrale particulière de la proposée.

Les trois dernières donnent

$a={\frac {{\cfrac {x^{2}}{4}}+x{\cfrac {dy}{dx}}}{1+x^{2}}},\quad b={\frac {{\cfrac {dy}{dx}}-{\cfrac {x^{2}}{4}}}{1+x^{2}}},$

et ces valeurs étant substituées dans la première, on aura

$y={\frac {-x^{4}+\left(8x^{3}+16x\right){\cfrac {dy}{dx}}+16{\cfrac {dy^{2}}{dx^{2}}}}{16\left(1+x^{2}\right)}}\,;$

c’est l’intégrale particulière aux premières différences de l’équation différentio-différentielle dont il s’agit.

Si l’on intègre cette équation, on aura alors l’intégrale particulière finie de la proposée. Pour cela, je tire par l’extraction de la racine carrée la valeur de ${\frac {dy}{dx}},$ j’ai

$4{\frac {dy}{dx}}+2x+x^{2}={\sqrt {1+x^{2}}}{\sqrt {16y+4x^{2}+x^{4}}}\,;$

donc, divisant par ${\sqrt {16y+4x^{2}+x^{4}}}$ et multipliant par $2,$ on aura

${\frac {8dy+4xdx+2x^{3}dx}{\sqrt {16y+4x^{2}+x^{4}}}}=2dx{\sqrt {1+x^{2}}},$

équation intégrable, et dont l’intégrale est, en ajoutant une constante arbitraire $\alpha ,$

${\sqrt {16y+4x^{2}+x^{4}}}=x{\sqrt {1+x^{2}}}-\log \left({\sqrt {1+x^{2}}}-x\right)+\alpha .$

Il est remarquable que tandis que l’intégrale complète de la proposée est algébrique, l’intégrale particulière en est transcendante.

31. L’intégrale particulière aux différences premières que nous avons trouvée ci-dessus admet, outre l’intégrale complète précédente, encore une intégrale particulière finie, qu’on peut trouver par les méthodes des Articles I et II. Déduisons-la de l’intégrale complète au moyen de la condition ${\frac {dy}{d\alpha }}=0\,;$ la différentiation de la dernière équation donne

${\frac {dy}{d\alpha }}={\frac {\sqrt {16y+4x^{2}+x^{4}}}{8}}\,;$

ainsi l’on aura

${\sqrt {16y+4x^{2}+x^{4}}}=0\,;$

comme cette équation ne renferme point la quantité $\alpha ,$ il faut la combiner avec l’intégrale complète en éliminant l’une des variables $x$ ou $y,$ pour voir si la valeur résultante de $\alpha$ est constante ou variable (4) ; or l’équation que nous venons de trouver donne

$y=-{\frac {x^{2}}{4}}-{\frac {x^{4}}{16}}\,;$

et cette valeur étant substituée dans l’intégrale complète (numéro précédent), on a

$x{\sqrt {1+x^{2}}}-\log \left({\sqrt {1+x^{2}}}-x\right)+\alpha =0\,;$

d’où l’on voit que $\alpha$ est déterminée par une fonction de $x\,;$ par conséquent l’équation dont il s’agit est une intégrale particulière.

Mais cette intégrale particulière, quoiqu’elle satisfasse à l’intégrale particulière aux premières différences de l’équation proposée, il ne s’ensuit pas qu’elle doive aussi satisfaire à cette dernière équation ; au contraire, elle n’y satisfera pas à moins que l’on n’ait en même temps

${\frac {dy}{d\alpha }}=0,\quad {\frac {d^{2}y}{dxd\alpha }}=0,$

à cause qu’il s’agit d’une équation différentielle du second ordre (11) ; or ayant

${\frac {dy}{d\alpha }}={\frac {\sqrt {16y+4x^{2}+x^{4}}}{8}},$

on aura

${\frac {d^{2}y}{dxd\alpha }}={\frac {8{\cfrac {dy}{dx}}+4x+2x^{3}}{8{\sqrt {16y+4x^{2}+x^{4}}}}},$

et mettant pour ${\frac {dy}{dx}}$ sa valeur

$-{\frac {x}{2}}-{\frac {x^{3}}{4}}-{\frac {1}{4}}{\sqrt {1+x^{2}}}{\sqrt {16y+4x^{2}+x^{4}}},$

il viendra

${\frac {d^{2}y}{dxd\alpha }}={\frac {\sqrt {1+x^{2}}}{4}},$

ce qui n’est pas nul ; d’où il s’ensuit que l’équation dont il s’agit, savoir

$y=-{\frac {x^{2}}{4}}-{\frac {x^{4}}{16}},$

ne satisfait pas à la proposée du second ordre, comme on peut aisément s’en assurer. Cet Exemple peut servir de confirmation à la théorie donnée dans l’Article II,

Au reste, il est bon de remarquer que cette intégrale particulière

$y=-{\frac {x^{2}}{4}}-{\frac {x^{4}}{16}}$

peut aussi se déduire immédiatement de l’intégrale finie et complète

$y={\frac {ax^{2}}{2}}+bx+a^{2}+b^{2},$

en faisant en même temps

${\frac {dy}{da}}=0\quad {\text{et}}\quad {\frac {dy}{db}}=0\,;$

ce qui donne les deux équations

${\frac {x^{2}}{2}}+2a=0,\quad x+2b=0,$

d’où l’on tire

$a=-{\frac {x^{2}}{4}},\quad b=-{\frac {x}{2}},$

ce qui, étant substitué dans l’intégrale complète, donne

$y=-{\frac {x^{4}}{16}}-{\frac {x^{2}}{4}}.$

Et cette règle est générale pour toutes les équations différentio-différentielles dont on connaît l’intégrale finie et complète.

32. Si, au moyen de l’équation finie $\mathrm {V} =0$ et de l’équation aux premières différences ${\frac {dy}{dx}}-p=0$ qui en est dérivée par la différentiation, on élimine l’une des deux constantes $a,b,$ on a une équation différentielle du premier ordre $\mathrm {Z} =0,$ qui sera l’intégrale complète aux différences premièresde l’équation différentio-différentielle $\mathrm {Z} '=0\,;$ et, comme on peut éliminer à volonté l’une ou l’autre des deux constantes arbitraires $a,b,$ on aura ainsi deux intégrales aux premières différences ; ce qui est connu des Géomètres.

Supposons maintenant qu’on ait éliminé $b,$ en sorte que dans l’équation $\mathrm {Z} =0$ la quantité $\mathrm {Z}$ soit une fonction de $x,y,{\frac {dy}{dx}}$ et $a\,;$ si l’on différentie cette équation en faisant varier $x,y$ et $a,$ et qu’on suppose, en général,

$d\mathrm {Z} =\mathrm {A} d{\frac {dy}{dx}}+\mathrm {B} dy+\mathrm {C} dx+\mathrm {E} da,$

on aura

$d{\frac {dy}{dx}}=-{\frac {\mathrm {B} dy+\mathrm {C} dx+\mathrm {E} da}{\mathrm {A} }}\,;$

donc, en faisant varier $a$ seul, on aura

${\frac {d^{2}y}{dxda}}=-\mathrm {\frac {B}{A}} {\frac {dy}{da}}-\mathrm {\frac {E}{A}} ,$

et, faisant varier $b$ seul,

${\frac {d^{2}y}{dxdb}}=-\mathrm {\frac {B}{A}} {\frac {dy}{db}}\,;$

ces valeurs étant substituées dans l’équation de condition

${\frac {d^{2}y}{dxda}}+{\frac {d^{2}y}{dxdb}}{\frac {db}{da}}=0,$

du n^o 29, on aura

$-\mathrm {\frac {B}{A}} \left({\frac {dy}{da}}+{\frac {dy}{db}}{\frac {db}{da}}\right)-\mathrm {\frac {E}{A}} =0\,;$

mais on doit avoir aussi (numéro cité)

${\frac {dy}{da}}+{\frac {dy}{db}}{\frac {db}{da}}=0\,;$

donc on aura

$-\mathrm {\frac {E}{A}} =0\,;$

or si dans l’équation $\mathrm {Z} =0$ on fait varier uniquement ${\frac {dy}{dx}}$ et $a,$ en regardant $y$ et $x$ comme constantes, on a

$\mathrm {A} d{\frac {dy}{dx}}+\mathrm {E} da=0,$

d’où

$\mathrm {\frac {E}{A}} =-{\frac {d{\cfrac {dy}{dx}}}{da}}=-{\frac {d^{2}y}{dxda}}\,;$

ainsi l’équation de condition se réduira à

${\frac {d^{2}y}{dxda}}=0\,;$

laquelle, étant combinée avec l’équation $\mathrm {Z} =0,$ donnera par l’élimination de $a$ la même équation qu’on eût obtenue d’après les quatre équations du n^o 29.

Et si au lieu d’éliminer $b$ on eût éliminé $a,$ en sorte que $\mathrm {Z}$ fût une fonction de $x,y,{\frac {dy}{dx}}$ et $b,$ alors on aurait l’équation de condition

${\frac {d^{2}y}{dxdb}}=0,$

laquelle donnerait encore le même résultat en éliminant $b$ au moyen de l’équation $\mathrm {Z} =0.$

33. Il s’ensuit de là que si l’on ne connait pas l’intégrale finie et complète $\mathrm {V} =0$ de l’équation différentio-différentielle $\mathrm {Z} '=0,$ mais seulement une des deux intégrales aux premières différences de cette équation, telle que $\mathrm {Z} =0,$ $\mathrm {Z}$ étant une fonction de $x,y,{\frac {dy}{dx}}$ et d’une constante arbitraire $a,$ on pourra également trouver l’intégrale particulière de la même équation $\mathrm {Z} '=0\,;$ pour cela il n’y aura qu’à faire varier dans l’équation $\mathrm {Z} =0$ les deux quantités ${\frac {dy}{dx}}$ et $a,$ et à supposer ensuite

${\frac {d^{2}y}{dxda}}=0\,;$

cette équation, étant combinée avec l’équation $\mathrm {Z} =0,$ en éliminant la quantité $a,$ donnera l’intégrale cherchée.

Cette règle peut aussi se démontrer directement, et indépendamment de la considération de l’intégrale finie et complète $\mathrm {V} =0.$ En effet, puisque l’équation du premier ordre $\mathrm {Z} =0$ satisfait à l’équation du second ordre $\mathrm {Z} '=0,$ quelle que soit la valeur de la constante $a$ contenue danse, il s’ensuit que cette équation $\mathrm {Z} '=0$ ne peut être que le résultat de l’élimination de $a$ au moyen de l’équation $\mathrm {Z} =0$ et de l’équation $d{\frac {dy}{dx}}=p'dx$ déduite de celle-là au moyen d’une différentiation. Or il est clair que ce résultat sera toujours le même, quelle que soit la quantité à éliminer $a,$ constante ou non, pourvu que les deux équations

$\mathrm {Z} =0,\quad d{\frac {dy}{dx}}=p'dx$

soient les mêmes ; mais en regardant $a$ comme variable, on a

$d{\frac {dy}{dx}}=p'dx+q'da\,;$

équation qui se réduira à la forme précédente en faisant $q'=0\,;$ donc si l’on détermine $a$ en sorte que $q',$ ou ce qui est la même chose ${\frac {d{\cfrac {dy}{dx}}}{da}},$ ou bien ${\frac {d^{2}y}{dxda}},$ soit nul, l’équation $\mathrm {Z} =0$ satisfera encore à l’équation $\mathrm {Z} '=0\,;$ et, comme $a$ devient dans ce cas égal à une quantité variable, l’équation $\mathrm {Z} =0$ ne sera plus qu’une intégrale particulière de la même équation.

34. Pour confirmer cette règle par un Exemple, reprenons l’équation différentio-différentielledu n^o 30, et nous trouverons aisément, d’après l’intégrale finie et complète qu’on connaît déjà, ces deux intégrales aux premières différences

${\begin{aligned}y=&\left(-{\frac {a}{2}}+a^{2}\right)x^{2}+(1-2a)x{\frac {dy}{dx}}+a^{2}+{\frac {d^{2}y}{dx^{2}}},\\y=&{\frac {\left(b+{\cfrac {dy}{dx}}\right)x}{2}}+{\frac {\left(b-{\cfrac {dy}{dx}}\right)^{2}}{x^{2}}}+b^{2},\end{aligned}}$

$a$ et $b$ étant les constantes arbitraires.

Faisant varier dans la première les quantités ${\frac {dy}{dx}}$ et $a,$ on en tire

${\frac {d^{2}y}{dxda}}={\frac {\left({\cfrac {1}{2}}-2a\right)x^{2}+2x{\cfrac {dy}{dx}}-2a}{(1-2a)x+2{\cfrac {dy}{dx}}}},$

et supposant cette quantité égale à zéro, on aura l’équation

$\left({\frac {1}{2}}-2a\right)x^{2}+2x{\frac {dy}{dx}}-2a=0,$

d’où résulte

$a={\frac {{\cfrac {x^{2}}{2}}+2x{\cfrac {dy}{dx}}}{2\left(1+x^{2}\right)}}\,;$

et cette valeur étant substituée dans l’équation ci-dessus, il viendra,

après les réductions,

$y=x{\frac {dy}{dx}}+{\frac {d^{2}y}{dx^{2}}}-{\frac {\left(4x{\cfrac {dy}{dx}}+x^{2}\right)^{2}}{16\left(1+x^{2}\right)}},$

ou bien

$y={\frac {-x^{4}+\left(8x^{3}+16x\right){\cfrac {dy}{dx}}+16{\cfrac {d^{2}y}{dx^{2}}}}{16\left(1+x^{2}\right)}}\,;$

c’est, comme l’on voit, la même équation qu’on a trouvée dans le n^o 30.

On trouvera encore le même résultat si, dans la seconde équation ci-dessus, on fait varier les quantités ${\frac {dy}{dx}}$ et $b,$ et qu’on suppose ensuite ${\frac {d^{2}y}{dxdb}}=0\,;$ on aura en effet, par la différentiation de cette équation,

${\frac {d^{2}y}{dxdb}}={\frac {{\cfrac {x^{3}}{2}}+2bx^{2}+2\left(b-{\cfrac {dy}{dx}}\right)}{2\left(b-{\cfrac {dy}{dx}}\right)-{\cfrac {x^{3}}{2}}}},$

ce qui étant supposé égal à zéro donne

$b={\frac {2{\cfrac {dy}{dx}}-{\cfrac {x^{3}}{2}}}{2\left(1+x^{2}\right)}},$

et cette valeur étant substituée à la place de $b,$ on aura, après les réductions,

$y={\frac {x}{2}}{\frac {dy}{dx}}+{\frac {1}{x^{2}}}{\frac {y^{2}}{x^{2}}}-{\frac {\left(4{\cfrac {dy}{dx}}-x^{3}\right)^{2}}{16x^{2}\left(1+x^{2}\right)}},$

ou bien

$y={\frac {16x{\cfrac {dy}{dx}}+8x^{3}{\cfrac {dy}{dx}}+16{\cfrac {dy^{2}}{dx^{2}}}-x^{4}}{16\left(1+x^{2}\right)}}.$

Ainsi les deux intégrales aux premières différences, quoique très-différentes entre elles, donnent cependant la même intégrale particulière ; la raison en est claire par l’analyse du n^o 32.

35. La méthode du n^o 33 pour trouver l’intégrale particulière d’une équation différentio-différelltielle lorsqu’on connaît seulement une de ses intégrales complètes aux premières différences est, comme l’on voit, absolument analogue à celle du n^o 4 pour trouver l’intégrale particulière d’une équation différentielle du premierordre au moyen de son intégrale finie et complète ; l’une et l’autre sont fondées sur les mêmes principes, et doivent par conséquent donner lieu à des conséquences semblables. Ainsi tout ce qu’on a dit dans l’Article II relativement à l’étendue des intégrales particulières des équations différentielles du premier ordre, pourra s’appliquer aussi aux intégrales particulières des équations différentio-différentielles.

Donc, si l’équation différentio-différentielle $\mathrm {Z} '=0$ est elle-même l’intégrale d’une équation différentielle du troisième ordre $\mathrm {Z} ''=0,$ l’intégrale particulière de l’équation $\mathrm {Z} '=0,$ déduite de la condition ${\frac {d^{2}y}{dxda}}=0,$ ne satisfera pas à l’équation $\mathrm {Z} ''=0,$ à moins que l’on n’ait en même temps ${\frac {d^{2}y}{dxda}}=0$ et ${\frac {d^{3}y}{dx^{2}da}}=0\,;$ et ainsi de suite.

Si l’on a

${\frac {d^{2}y}{dxda}}=0,\quad {\frac {d^{3}y}{dx^{2}da}}=0,\ldots$ à l’infini

on prouvera, comme dans le n^o 13, que la condition ${\frac {d^{2}y}{dxda}}=0$ ne donnera plus une intégrale particulière ; et de là, par un raisonnement semblable à celui du n^o 14, on déduira une règle pour trouver immédiatement l’intégrale particulière d’une équation différentio-différentielle $\mathrm {Z} '=0$ sans connaître aucune de ses intégrales complètes.

Cette règle consiste à supposer égale à ${\frac {0}{0}}$ la valeur de ${\frac {d^{3}y}{dx^{3}}}$ tirée de l’équation proposée $\mathrm {Z} '=0$ au moyen de la différentiation ; on aura ainsi deux équations en $x,y,{\frac {dy}{dx}},{\frac {d^{2}y}{dx^{2}}},$ d’où éliminant ${\frac {d^{2}y}{dx^{2}}}$ à l’aide de la même équation $\mathrm {Z} '=0,$ on aura deux équations en $x,y,{\frac {dy}{dx}}$ qui devront s’accorder entre elles et se réduire à une même équation, si la proposée est susceptible d’une intégrale particulière ; et alors cette équation sera l’intégrale particulière cherchée.

36. Si l’équation différentio-différentielle $\mathrm {Z} '=0$ était telle, que l’on eût

$d\mathrm {Z} '=\mathrm {A} 'd{\frac {d^{2}y}{dx^{2}}},$

alors la condition de ${\frac {d^{3}y}{dx^{3}}}={\frac {0}{0}}$ donnerait cette équation unique $\mathrm {A} '=0\,;$ par conséquent, en chassant la quantité ${\frac {d^{2}y}{dx^{2}}}$ au moyen des deux équations

$\mathrm {A} '=0,\quad \mathrm {Z} '=0,$

la résultante sera toujours une intégrale particulière de la proposée $\mathrm {Z} '=0.$

Or dans ce cas on peut aussi trouver aisément l’intégrale finie et complète de la même équation. En effet, puisque $\mathrm {Z} '=0,$ on aura aussi, en différentiant,

$d\mathrm {Z} '=\mathrm {A} 'd{\frac {d^{2}y}{dx^{2}}}=0\,;$

donc ou $\mathrm {A} '=0,$ ce qui, comme nous venons de le voir, donne l’intégrale particulière, ou $d{\frac {d^{2}y}{dx^{2}}}=0,$ et par conséquent

${\frac {d^{2}y}{dx^{2}}}=a,\quad {\frac {dy}{dx}}=ax+b,\quad y={\frac {ax^{2}}{2}}+bx+c,$

$a,b,c$ étant trois constantes arbitraires ; or comme l’équation $\mathrm {Z} '=0$ n’est (hypothèse) que du second ordre, il s’ensuit que son intégrale finie et complète ne peut renfermer que deux constantes arbitraires ; ainsi, si l’on y substitue les valeurs précédentes de $y,{\frac {dy}{dx}}$ et ${\frac {d^{2}y}{dx^{2}}},$ il viendra nécessairement une équation entre les constantes $a,b,c,$ sans $x$ ni $y,$ par laquelle il faudra déterminer l’une de ces constantes par les deux autres, qui resteront par conséquent arbitraires.

De là on peut déduire la forme générale de ces sortes d’équations ; car soit, en général, $c$ une fonction de $a$ et $b$ représentée par $f{a,b},$ et puisqu’on a

$a={\frac {d^{2}y}{dx^{2}}},\quad b={\frac {dy}{dx}}-ax={\frac {dy}{dx}}-x{\frac {d^{2}y}{dx^{2}}},$

on aura

$c=f\left({\frac {d^{2}y}{dx^{2}}},{\frac {dy}{dx}}-x{\frac {d^{2}y}{dx^{2}}}\right)\,;$

donc, substituant ces valeurs dans l’équation

$y={\frac {ax^{2}}{2}}+bx+c,$

on aura celle-ci

$y=x{\frac {dy}{dx}}-{\frac {x^{2}}{2}}{\frac {d^{2}y}{dx^{2}}}+f\left({\frac {d^{2}y}{dx^{2}}},{\frac {dy}{dx}}-x{\frac {d^{2}y}{dx^{2}}}\right),$

ou bien, si Pon fait pour plus de simplicité ${\frac {dy}{dx}}=p,\ {\frac {dp}{dx}}=p',$

$y=px-{\frac {p'x^{2}}{2}}+f(p',p-xp')\,;$

toute équation donc qui sera réductible à cette forme aura, comme celle du n^o 17, la propriété de pouvoir être facilement intégrée au moyen d’une nouvelle différentiation, et son intégrale finie et complète sera

$y={\frac {ax^{2}}{2}}+bx+f(a,b),$

laquelle représente toujours une parabole.

De plus, l’équation précédente aura la propriété d’admettre toujours une intégrale particulière, qu’on trouvera en éliminant $p'$ au moyen de l’équation

$-{\frac {x^{2}}{2}}+{\frac {df(p',p-p'x)}{dp'}}=0,$

et qui pourra représenter différentes courbes.

L’équation qui a servi d’exemple dans le n^o 30 est comprise sous la forme précédente.

37. La théorie que nous venons de donner sur les intégrales particulières des équations différentielles du premier et du second ordre peut s’appliquer aisément aux équations d’un ordre quelconque plus élevé.

Soit par exemple $\mathrm {Z} ''=0$ une équation différentielle du troisième ordre, $\mathrm {Z} ''$ étant une fonction de $x,y,{\frac {dy}{dx}},{\frac {d^{2}y}{dx^{2}}},{\frac {d^{3}y}{dx^{3}}},$ si l’on connait son intégrale tinie et complète $\mathrm {V} =0,$ $\mathrm {V}$ étant une fonction de $x,y$ et de trois constantes arbitraires $a,b,c,$ on déterminera l’intégrale particulière de l’équation $\mathrm {Z} ''=0$ en éliminant les trois quantités $a,b,c$ et les deux ${\frac {db}{da}},{\frac {dc}{da}}$ de ces six équations

${\begin{aligned}&\mathrm {V} =0,\quad {\frac {dy}{dx}}-p=0,\quad {\frac {d^{2}y}{dx^{2}}}-p'=0,\\&{\frac {dy}{da}}+{\frac {dy}{db}}{\frac {db}{da}}+{\frac {dy}{dc}}{\frac {dc}{da}}=0,\\&{\frac {d^{2}y}{dadx}}\ \ +{\frac {d^{2}y}{dbdx}}{\frac {db}{da}}\ \,+{\frac {d^{2}y}{dcdx}}{\frac {dc}{da}}\ \,=0,\\&{\frac {d^{3}y}{dadx^{2}}}+{\frac {d^{3}y}{dbdx^{2}}}{\frac {db}{da}}+{\frac {d^{3}y}{dcdx^{2}}}{\frac {dc}{da}}=0.\end{aligned}}$

Si l’on connaît seulement une des intégrales complètes aux premières différences telle que $\mathrm {Z} =0,$ $\mathrm {Z}$ étant une fonction de $x,y,{\frac {dy}{dx}}$ et de deux constantes arbitraires $a$ et $b,$ on déterminera l’intégrale particulière en éliminant les quantités $a,b,{\frac {db}{da}},$ au moyen de ces quatre équations

${\begin{aligned}&\mathrm {Z} =0,\quad {\frac {d^{2}y}{dx^{2}}}-p'=0,\\&{\frac {d^{2}y}{dxda}}+{\frac {d^{2}y}{dxdb}}{\frac {db}{da}}=0,\\&{\frac {d^{3}y}{dx^{2}da}}+{\frac {d^{3}y}{dx^{2}db}}{\frac {db}{da}}=0.\end{aligned}}$

Si l’on ne connaît qu’une intégrale complète aux secondes différences, telle que $\mathrm {Z} '=0,$ $\mathrm {Z} '$ étant une fonction de $x,y,{\frac {dy}{dx}},{\frac {d^{2}y}{dx^{2}}}$ et d’une constante arbitraire $a,$ on pourra déterminer l’intégrale particulière en éliminant $a$ au moyen des deux équations

$\mathrm {Z} '=0,\quad {\frac {d^{3}y}{dx^{2}da}}=0.$

Enfin on pourra aussi déterminer cette intégrale d’après la seule équation différentielle $\mathrm {Z} ''=0\,;$ pour cela il n’y aura qu’à chercher par la différentiation la valeur de ${\frac {d^{4}y}{dx^{4}}}$ et la supposer égale à ${\frac {0}{0}}\,;$ les deux équations qu’on aura de cette manière devront revenir à la même, après l’élimination de ${\frac {d^{3}y}{dx^{3}}}$ faite par le moyen de la proposée $\mathrm {Z} ''=0,$ si celle-ci est susceptible d’une intégrale particulière, et alors l’équation résultante de l’élimination dont il s’agit sera l’intégrale particulière en question.

On pourra faire au reste sur les intégrales particulières des équations différentielles du troisième ordre des remarques analogues à celles qu’on a faites plus haut sur les intégrales particulières des équations du second ordre ; c’est sur quoi nous ne croyons pas qu’il soit nécessaire de nous étendre davantage.

38. Nous terminerons cet Article par faire remarquer qu’il y $a,$ dans chaque ordre, des équations différentielles qui ont des propriétés analogues à celles des équations des n^os 17 et 36.

Soit $\mathrm {Z} ''=0$ la forme générale de ces sortes d’équations pour le troisième ordre, on aura par la différentiation

$d\mathrm {Z} ''=\mathrm {A} ''d{\frac {d^{3}y}{dx^{3}}}=0\,;$

d’où l’on tire ou $\mathrm {A} ''=0,$ ce qui donnera une intégrale particulière après l’élimination de la quantité ${\frac {d^{3}y}{dx^{3}}},$ parce qu’alors ${\frac {d^{4}y}{dx^{4}}}={\frac {0}{0}}\,;$ ou bien $d{\frac {d^{3}y}{dx^{3}}}=0,$ ce qui donnera l’intégrale complète

${\frac {d^{3}y}{dx^{3}}}=a,$

d’où

${\frac {d^{2}y}{dx^{2}}}=ax+b,\quad {\frac {dy}{dx}}={\frac {ax^{2}}{2}}+bx+c,\quad y={\frac {ax^{3}}{2.3}}+{\frac {bx^{2}}{2}}+cx+e\,;$

or comme l’équation $\mathrm {Z} ''=0$ n’est que du troisième ordre, son intégrale finie et complète ne peut renfermer que trois constantes arbitraires ; par conséquent, si l’on substitue dans cette équation les valeurs précédentes de $y,{\frac {dy}{dx}},$${\frac {d^{2}y}{dx^{2}}},$${\frac {d^{3}y}{dx^{3}}},$ il arrivera nécessairement qu’on aura une équation entre les constantes $a,b,c,e$ sans $x$ ni $y,$ par laquelle il faudra déterminer une de ces constantes par les trois autres.

Supposons donc, en général, que l’on ait

$e=f(a,b,c),$

et si l’on fait pour plus de simplicité

${\frac {dy}{dx}}=p,\quad {\frac {dp}{dx}}=p',\quad {\frac {dp'}{dx}}=p'',$

on aura

$a=p'',\quad b=p'-xp'',\quad c=p-xp'+{\frac {x^{2}p''}{2}}\,;$

donc substituant ces valeurs dans l’équation

$y={\frac {ax^{3}}{2.3}}+{\frac {bx^{2}}{2}}+cx+e,$

on aura pour la forme générale des équations dont il s’agit

$y={\frac {x^{3}p''}{2.3}}-{\frac {x^{2}p'}{2}}+xp+f\left(p'',p'-xp'',p-xp'+{\frac {x^{2}p''}{2}}\right).$

L’intégrale finie et complète sera

$y={\frac {ax^{3}}{2.3}}+{\frac {bx^{2}}{2}}+cx+f(a,b,c)\,;$

et l’intégrale particulière se trouvera en éliminant $p''$ au moyen de l’équation

${\frac {x^{3}}{2.3}}+{\frac {df\left(p'',p'-xp'',p-xp'+{\frac {x^{2}p''}{2}}\right)}{dp''}}=0.$

Article V. — Des intégrales particulière des équations aux différences partielles, avec des remarques nouvelles sur la nature et sur l’intégration de ces sortes d’équations.

39. Soit

$\mathrm {V} =0$

une équation entre trois variables $x,y,z$ et deux constantes $a,b\,;$ je dis qu’on en peut déduire une équation à différences partielles du premier ordre dans laquelle les constantes $a$ et $b$ ne se trouvent plus. En effet, supposons qu’en faisant varier à la fois, $x,y$ et $z$ on ait

$dz=pdx+qdy,$

$p$ et $q$ étant des fonctions connues de $x,y,z,a$ et $b\,;$ donc en regardant $z$ comme une fonction de $x$ et $y,$ on aura, suivant la notation ordinaire des différences partielles,

${\frac {dz}{dx}}=p,\quad {\frac {dz}{dy}}=q\,;$

qu’on élimine les deux quantités $a$ et $b$ dans les trois équations

$V=0,\quad {\frac {dz}{dx}}-p=0,\quad {\frac {dz}{dy}}-q=0,$

et l’on aura pour résultante une équation entre $x,y,z$ et ${\frac {dz}{dx}},{\frac {dz}{dy}},$ où les constantes $a$ et $b$ ne se trouveront plus, et que nous représenterons, en général, par

$\mathrm {Z} =0.$

On peut donc regarder l’équation finie $\mathrm {V} =0$ comme l’intégrale complète de l’équation aux différences partielles du premier ordre $\mathrm {Z} =0\,;$ et comme les deux constantes $a$ et $b$ demeurent arbitraires dans l’équation $\mathrm {V} =0,$ il s’ensuit que l’intégrale complète de toute équation aux différences partielles du premier ordre entre trois variables doit nécessairement renfermer deux constantes arbitraires.

Soit par exemple l’équation finie

$z=a+bx+mby\,;$

on a par la différentiation

$dz=bdx+mbdy\,;$

donc

${\frac {dz}{dx}}=b,\quad {\frac {dz}{dy}}=mb,$

et éliminant $a$ et $b,$ on aura

$m{\frac {dz}{dx}}={\frac {dz}{dy}},$

équation dont l’intégrale complète sera donc

$z=a+b(x+my),$

$a$ et $b$ étant arbitraires.

Soit l’équation

$z=a+bx+cy,$

dans laquelle $c$ soit une fonction quelconque de $a$ et $b,$ que nous désignerons par $f(a,b)\,;$ on aura par la différentiation

${\frac {dz}{dx}}=b,\quad {\frac {dz}{dy}}=c=f(a,b)\,;$

donc

$z=a+x{\frac {dz}{dx}}+y{\frac {dz}{dy}},$

d’où

$a=z-{\frac {dz}{dx}}-y{\frac {dz}{dy}}\,;$

par conséquent l’équation différentielle sera

${\frac {dz}{dy}}=f\left(z-{\frac {dz}{dx}}-y{\frac {dz}{dy}},{\frac {dz}{dx}}\right),$

laquelle aura pour intégrale complète

$z=a+bx+f(a,b)y.$

Soit encore l’équation

$z=c+ax+by\quad {\text{et}}\quad c=f(a,b),$

on aura par la différentiation

${\frac {dz}{dx}}=a,\quad {\frac {dz}{dy}}=b\,;$

donc l’équation différentielle sera

$z=x{\frac {dz}{dx}}+y{\frac {dz}{dy}}+f\left({\frac {dz}{dx}},{\frac {dz}{dy}}\right),$

et son intégrale complète sera

$z=f(a,b)+ax+by.$

40. Nous avons supposé les quantités $a$ et $b$ constantes ; mais si elles étaient variables on parviendrait toujours à la même équation différentielle $\mathrm {Z} =0,$ pourvu que l’on eût également

$dz=pdx+qdy,$

comme dans le cas où ces quantités seraient constantes ; car il est clair que le résultat de l’élimination de $a$ et $b,$ dans les équations :

$\mathrm {V} =0,\quad {\frac {dz}{dx}}-p=0,\quad {\frac {dz}{dy}}=0,$

sera toujours le même, quelles que soient les valeurs de $a$ et $b.$ Or en faisant varier à la fois les quantités $x,y,a$ et $b,$ on aura

$dz=pdx+qdy+rda+sdb\,;$

donc la condition dont il s’agit aura lieu si $rda+sdb=0,$ par consé-

quent si l’on détermine les quantités $a$ et $b$ en sorte que l’on ait

$rda+sdb=0,$

et qu’on substitue ensuite leurs valeurs dans l’équation finie $\mathrm {V} =0,$ on aura une nouvelle intégrale de l’équation proposée $\mathrm {Z} =0.$

41. La manière la plus simple de satisfaire à l’équation

$rda+sdb=0$

est de supposer séparément

$r=0,\quad s=0,$

ce qui donne deux équations qui serviront à déterminer $a$ et $b.$ Or, en regardant $z$ comme fonction de $a$ et $b,$ il est visible que

$r={\frac {dz}{da}},\quad s={\frac {dz}{db}}\,;$

donc les deux conditions dont il s’agit seront représentées par

${\frac {dz}{da}}=0,\quad {\frac {dz}{db}}=0\,;$

lesquelles étant analogues à la condition ${\frac {dy}{da}}=0$ que nous avons trouvée dans l’Article I pour la détermination des intégrales particulières des équations à deux variables, on pourra regarder aussi les intégrales provenantes de ces conditions comme des intégrales particulières des équations à différences partielles entre trois variables.

42. Prenons par exemple l’équation

$z=x{\frac {dz}{dx}}+y{\frac {dz}{dy}}+f\left({\frac {dz}{dx}},{\frac {dz}{dy}}\right),$

dont nous avons déjà vu que l’intégrale complète est (39)

$z=f(a,b)+ax+by.$

En faisant varier successivement $a$ et $b,$ on aura

${\frac {dz}{da}}={\frac {df(a,b)}{da}}+x,\quad {\frac {dz}{db}}={\frac {df(a,b)}{db}}+y\,;$

ainsi on aura, pour la détermination de l’intégrale particulière, les deux équations

$x+{\frac {df(a,b)}{da}}=0,\quad y+{\frac {df(a,b)}{db}}=0,$

au moyen desquelles on éliminera $a$ et $b$ de l’équation

$z=f(a,b)+ax+by\,;$

et la résultante sera l’intégrale particulière qu’on cherche.

Supposons que l’équation proposée soit

$z=y+{\frac {dz}{dy}}+x+{\frac {dz}{dx}}+h{\sqrt {1+\left({\frac {dz}{dx}}\right)^{2}+\left({\frac {dz}{dy}}\right)^{2}}},$

on aura dans ce cas

$f\left({\frac {dz}{dx}},{\frac {dz}{dy}}\right)=h{\sqrt {1+\left({\frac {dz}{dx}}\right)^{2}+\left({\frac {dz}{dy}}\right)^{2}}}\,;$

donc

$f(a,b)=h{\sqrt {1+a^{2}+b^{2}}},$

${\frac {df(a,b)}{da}}={\frac {ha}{\sqrt {1+a^{2}+b^{2}}}},\quad {\frac {df(a,b)}{db}}={\frac {hb}{\sqrt {1+a^{2}+b^{2}}}}\,;$

ainsi on aura d’abord cette intégrale complète

$z=h{\sqrt {1+a^{2}+b^{2}}}+ax+by\,;$

ensuite, pour avoir l’intégrale particulière, il n’y aura qu’à éliminer dans cette équation les quantités $a$ et $b$ au moyen de ces deux-ci

$x+{\frac {ha}{\sqrt {1+a^{2}+b^{2}}}}=0,\quad y+{\frac {hb}{\sqrt {1+a^{2}+b^{2}}}}=0,$

lesquelles donnent

$a=-{\frac {x}{\sqrt {h^{2}-x^{2}-y^{2}}}},\quad b=-{\frac {y}{\sqrt {h^{2}-x^{2}-y^{2}}}},$

${\sqrt {1+a^{2}+b^{2}}}={\frac {h}{\sqrt {h^{2}-x^{2}-y^{2}}}}\,;$

de sorte que l’intégrale particulière sera

$z={\sqrt {h^{2}-x^{2}-y^{2}}}.$

43. Pour rendre plus sensible l’analogie qu’il y a entre les intégrales particulières des équations aux différences partielles et celles des équations différentiellesà deux variables, on remarquera que, si l’on exprime les variables $x,y,z$ par les coordonnées rectangles d’une surface courbe, l’équation $\mathrm {V} =0$ pourra représenter une infinité de surfaces courbes, en donnant aux arbitraires $a$ et $b$ toutes les valeurs possibles, et chacune de ces différentes surfaces satisfera également à l’équation du premier ordre $\mathrm {Z} =0\,;$ ensuite on prouvera, par un raisonnement semblable à celui des n^os 21 et suivants, que la surface qui touchera toutes celles-ci satisfera aussi à la même équation différentielle $\mathrm {Z} =0\,;$ enfin on démontrera aisément que, pour avoir l’équation de la surface touchante dont il s’agit, il n’y aura qu’à éliminer $a$ et $b$ de l’équation $\mathrm {V} =0$ au moyen des deux équations

${\frac {dz}{da}}=0,\quad {\frac {dz}{db}}=0\,;$

d’où il s’ensuit que cette surface touchante exprimera l’intégrale particulière de l’équation $\mathrm {Z} =0\,;$ ce qui est conforme à la théorie donnée dans l’Article relativement aux lignes courbes.

44. Pour confirmer cette théorie par un exemple, supposons qu’on demande la surface courbe qui aura cette propriété, qu’en menant d’un point donné sur un quelconque des plans touchants de cette surface une perpendiculaire, elle soit toujours d’une même grandeur donnée.

Il est d’abord visible qu’une sphère décrite autour du point donné avec un rayon égal à la grandeur donnée satisfera à la question ; mais comme dans cette solution tout est déterminé et que le Problème conduit naturellement à une équation aux différences partielles du premier ordre, comme on le verra ci-après, il s’ensuit qu’on ne peut avoir de cette manière qu’une solution particulière du Problème.

De plus, il est clair qu’il y a une infinité de plans qui résolvent ce Problème ; car il suffit pour cela que la position du plan soit telle, que la perpendiculaire qu’on y abaisserait du point donné soit égale à la grandeur donnée ; et si l’on cherche l’équation générale de tous les, plans qui ont cette propriété, on verra sans peine que cette équation renfermera deux constantes arbitraires ; de sorte qu’on pourra la regarder comme l’intégrale complète de l’équation différentielle du Problème.

Enfin il est aisé de se convaincre que tous les différents plans dont il s’agit toucheront une surface sphérique décrite autour du point donné avec un rayon égal à la valeur donnée de la perpendiculaire ; c’est-à-dire la même surface qui nous a déjà donné une solution particulière du Problème.

Appliquons maintenant le calcul à cette question, et nommons $s,t,$ il les trois coordonnées rectangles qui répondent à un point quelconque d’un des plans touchants de la surface cherchée dont les coordonnées rectangles et parallèles à celles-là sont $x,y,z\,;$ on aura, par la nature du plan, l’équation

$u=l+ms+nt,$

$l,m,n$ étant des constantes. Or, puisque le plan et la surface passent par un même point, on aura dans ce point $s=x,\ t=y,\ u=z,$ donc

$z=l+mx+ny,\quad {\text{d’où}}\quad l=z-mx-ny\,;$

ensuite, puisque dans le même point le plan et la surface se touchent, on aura aussi

${\frac {du}{ds}}={\frac {dz}{dx}}\quad {\text{et}}\quad {\frac {du}{dt}}={\frac {dz}{dy}}\,;$

mais

$u{\frac {ds}{d=m,\quad }}u{\frac {dt}{d=n}}\,;$

donc

$m={\frac {dz}{dx}},\quad n={\frac {dz}{dy}},$

et par conséquent

$l=z-x{\frac {dz}{dx}}-y{\frac {dz}{dy}}\,;$

donc, substituant ces valeurs dans l’équation du plan touchant, elle deviendra

$u=s{\frac {dz}{dx}}+t{\frac {dz}{dy}}+z-x{\frac {dz}{dx}}-y{\frac {dz}{dy}}.$

Supposons, pour plus de simplicité, que le point donné soit l’origine des coordonnées, et cherchons l’expression générale de la perpendiculaire menée de ce point sur le plan dont l’équation est

$u=l+ms+nt\,;$

soit $\rho$ la valeur d’une ligne menée du point dont il s’agit à un point quelconque de ce plan, on aura, en général,

$\rho ^{2}=s^{2}+t^{2}+u^{2},$

et, dans le cas où cette ligne devient perpendiculaire, on aura $d\rho =0$ et par conséquent

$sds+tdt+udu=0\,;$

mais l’équation au plan donne

$du=mds+ndt,$

donc

$(s+mu)ds+(t+nu)dt=0,$

et par conséquent

$s+mu=0,\quad l+nu=0,$

d’où l’on tire

$s=-mu,\quad t=-nu\,;$

donc

$u=l-m^{2}u-n^{2}u\,;$

donc

${\begin{aligned}u={\frac {l}{1+m^{2}+n^{2}}},\quad \rho ^{2}=s^{2}+t^{2}+u^{2}&=\left(1+m^{2}+n^{2}\right)u^{2}\\&={\frac {l^{2}}{1+m^{2}+n^{2}}}\,;\end{aligned}}$

donc

$\rho ={\frac {l}{\sqrt {1+m^{2}+n^{2}}}}\,;$

donc, substituant pour $l,m,n$ les valeurs trouvées ci-dessus, on aura enfin pour l’expression générale de la perpendiculaire $\rho$

$\rho ={\frac {z-x{\cfrac {dz}{dx}}-y{\cfrac {dz}{dy}}}{\sqrt {1+\left({\cfrac {dz}{dx}}\right)^{2}+\left({\cfrac {dz}{dy}}\right)^{2}}}}.$

Supposant donc cette perpendiculaire égale à une constante donnée $h,$ on aura enfin

$z=x{\frac {dz}{dx}}+y{\cfrac {dz}{dy}}+h{\sqrt {1+\left({\cfrac {dz}{dx}}\right)^{2}+\left({\cfrac {dz}{dy}}\right)^{2}}}$

pour l’équation du Problème.

Cette équation est la même que nous avons déjà traitée plus haut (42), et dont nous avons trouvé que l’intégrale complète est

$z=ax+by+h{\sqrt {1+a^{2}+b^{2}}},$

et que l’intégrale particulière est

$z={\sqrt {h^{2}-x^{2}-y^{2}}}\,;$

ce qui s’accorde avec les conclusions trouvées ci-dessus.

45. Si l’on rapproche la théorie que nous venons de donner sur les intégrales particulières des équations aux différences partielles de celle que nous avons donnée plus haut sur les intégrales particulières des équations différentielles à deux variables, on en pourra déduire une règle analogue à celle des n^os 15 et 24 pour trouver les intégrales particulières sans connaître les intégrales complètes ; car on prouvera aisément, par des principes analogues à ceux qu’on a employés dans les endroits cités, que, dans le cas de l’intégrale particulière, les différences des quantités ${\frac {dz}{dy}},\ {\frac {dz}{dx}}$ déduites de l’équation différentielle proposée $\mathrm {Z} =0,$ au moyen d’une nouvelle différentiation, devront rester indéterminées.

Ainsi, si après avoir différentié l’équation $\mathrm {Z} =0,$ et avoir fait dispapaître les fractions on a une équation de la forme

$\mathrm {M} d{\frac {dz}{dx}}+\mathrm {N} d{\frac {dz}{dy}}+\mathrm {P} dx+\mathrm {Q} dy=0,$

$\mathrm {M,N,P,Q,R}$ étant des fonctions connues et entières de $x,y,z,{\frac {dz}{dx}},{\frac {dz}{dy}},$ il faudra, pour obtenir l’intégrale particulière de l’équation dont il s’agit, faire séparément les quantités $\mathrm {M,N,P,Q}$ chacune égale à zéro ; ce qui donnera, comme l’on voit, quatre équations, lesquelles étant combinées avec l’équation $\mathrm {Z} =0$ donneront, par l’élimination des deux quantités ${\frac {dz}{dx}},{\frac {dz}{dy}},$ trois équations finales en $x,y,z,$ qui devront avoir lieu en même temps. Par conséquent, si ces équations ont un facteur commun, ce facteur sera l’intégrale particulière cherchée sinon la proposée n’admettra point d’intégrale particulière.

46. Si l’équation $\mathrm {Z} =0$ était telle, que l’on eût par la différentiation

$d\mathrm {Z} =\mathrm {A} d{\frac {dz}{dx}}+\mathrm {B} d{\frac {dz}{dy}},$

alors on aurait l’équation différentielle

$\mathrm {A} d{\frac {dz}{dx}}+\mathrm {B} d{\frac {dz}{dy}}=0,$

et les conditions de l’intégrale particulière seraient remplies en faisant $\mathrm {A} =0$ et $\mathrm {B} =0\,;$ on n’aurait donc, dans ce cas, que deux équations de condition, lesquelles serviraient a éliminer les deux quantités ${\frac {dz}{dx}},\ {\frac {dz}{dy}}$

dans l’équation $\mathrm {Z} =0\,;$ et l’équation résultante serait toujours l’intégrale particulière de cette même équation $\mathrm {Z} =0.$

Au reste on peut aussi trouver son intégrale complète en remarquant que l’équation

$d\mathrm {Z} =\mathrm {A} d{\frac {dz}{dx}}+\mathrm {B} d{\frac {dz}{dy}}=0$

donne aussi

$d{\frac {dz}{dx}}=0,\quad d{\frac {dz}{dy}}=0,$

d’où

${\frac {dz}{dx}}=a,\quad {\text{et}}\quad z=ax+\mathrm {Y} ,$

$a$ étant une constante et $\mathrm {Y}$ une fonction de $y$ sans $x\,;$ mais, puisqu’on doit avoir en même temps $d{\frac {dz}{dy}}=0,$ il faudra que $d{\frac {d\mathrm {Y} }{dy}}=0,$ donc

${\frac {d\mathrm {Y} }{dy}}=b,\quad {\text{et}}\quad \mathrm {Y} =by+c,$

$b$ et $c$ étant des constantes ; donc

$z=ax+by+c\,;$

si l’on substitue cette valeur de $z$ dans l’équation différentielle $\mathrm {Z} =0,$ il arrivera nécessairement que les quantités $x$ et $y$ s’en iront et que l’on aura une équation entre les constantes $a,b,c,$ par laquelle il faudra en déterminer une par les deux autres.

Soit donc, en général, $c=f(a,b),$ l’intégrale complète sera alors

$z=ax+by+f(a,b),$

et l’équation différentielles sera, comme on l’a déjà vu (39),

$z=x{\frac {dz}{dx}}+y{\frac {dz}{dy}}+f\left({\frac {dz}{dx}},{\frac {dz}{dy}}\right)\,;$

c’est la forme générale des équations différentielles qui peuvent avoir la propriété dont il s’agit.

En effet, si l’on différentie cette équation, on aura, à cause de

$dz={\frac {dz}{dx}}dx+{\frac {dz}{dy}}dy,$

celle-ci, où je mets, pour plus de simplicité, $p$ et $q$ à la place de ${\frac {dz}{dx}}$ et ${\frac {dz}{dy}},$

$\left[x+{\frac {df(p,q)}{dp}}\right]dp+\left[y+{\frac {df(p,q)}{dq}}\right]dq=0\,;$

ainsi l’on aura pour l’intégrale particulière les deux équations

$x+{\frac {df(p,q)}{dp}}=0,\quad y+{\frac {df(p,q)}{dq}}=0,$

à l’aide desquelles il faudra éliminer les quantités $p$ et $q$ dans la proposée

$z=px+qy+f(p,q).$

Et il est visible qu’on aura de cette manière le même résultat que par la méthode du n^o 42.

Les équations à différences partielles de la forme dont il s’agit répondent, comme l’on voit, à celles qu’on a considérées plus haut (17).

47. Nous avons vu ci-dessus (40) que, pour que l’équation finie $\mathrm {V} =0$ satisfasse à l’équation aux différences partielles du premier ordre $\mathrm {Z} =0,$ sans supposer que les arbitraires $a$ et $b$ soient constantes, il suffit que ces quantités soient telles, que l’on ait la condition

$rda+sdb=0,$

c’est-à-dire, à cause de

$r={\frac {dz}{da}},\quad s={\frac {dz}{db}}$

suivant la notation ordinaire,

${\frac {dz}{da}}da+{\frac {dz}{db}}db=0.$

Dans les n^os 41 et suivants nous avons satisfait à cette condition en faisant séparément ${\frac {dz}{da}}=0$ et ${\frac {dz}{db}}=0,$ ce qui nous a donné l’intégrale particulière de l’équation $\mathrm {Z} =0\,;$ mais il est clair que cette supposition est trop limitée et qu’on peut remplir la condition dont il s’agit d’une manière plus générale. En effet, comme il y a deux indéterminées $a$ et $b,$ on peut supposer que l’une soit une fonction quelconque de l’autre, par exemple,

$b=\varphi (a)$

(en prenant la caractéristique $\varphi$ pour dénoter une fonction indéterminée) ; substituant cette valeur de $b$ en $a$ et faisant

$d\varphi (a)=\varphi '(a)da,$

on aura l’équation de condition

${\frac {dz}{da}}+{\frac {dz}{db}}\varphi '(a)=0,$

au moyen de laquelle on pourra éliminer $a$ dans l’équation $\mathrm {V} =0.$ L’équation résultante de cette élimination satisfera également à l’équation différentielle $\mathrm {Z} =0,$ et comme elle renferme une fonction arbitraire, elle sera beaucoup plus générale que l’intégrale complète $\mathrm {V} =0\,;$ c’est pourquoi, et pour la distinguer de celle-ci, nous la nommerons l’intégrale générale de l’équation $\mathrm {Z} =0.$

48. Ainsi donc, connaissant l’intégrale complète d’une équation à différences partielles du premier ordre, on pourra toujours, par la méthode précédente, en déduire l’intégrale générale, laquelle résoudra la même équation dans toute son étendue. Il n’y aura pour cela qu’à supposer que dans l’intégrale complète l’une des deux arbitraires soit une fonction quelconque de l’autre, différentier ensuite cette intégrale en faisant varier uniquement l’arbitraire restante, et éliminer cette arbitraire. Appliquons cette méthode à quelques Exemples.

Soit proposée d’abord l’équation

${\frac {dz}{dy}}=m{\frac {dz}{dx}},$

dont nous avons vu ci-dessus que l’intégrale complète est

$z=a+b(x+my).$

Je fais $a=\varphi (b),$ j’ai

$z=\varphi (b)+b(x+my)\,;$

je différentie en faisant varier $b$ seul et divisant par $db,$ j’ai

$\varphi '(b)+x+my=0\,;$

au moyen de ces deux équations on éliminera $b$ et l’on aura l’intégrale générale. Or comme

$-\varphi '(b)=x+my,$

on aura $b$ égale à une fonction de $x+my\,;$ par conséquent

$\varphi (b)+b(x+my)$

sera aussi nécessairement une fonction de $x+my,$ fonction qui restera indéterminée, à cause que $\varphi (b)$ est une fonction indéterminée de $b.$ Ainsi l’on aura pour l’intégrale générale de la proposée

$z=\Phi (x+my),$

la caractéristique $\Phi$ dénotant une fonction quelconque.

49. Soit l’équation

$z=x{\frac {dz}{dx}}+y{\frac {dz}{dy}}+f\left({\frac {dz}{dx}},{\frac {dz}{dy}}\right),$

dont l’intégrale complète est (39)

$z=f(a,b)+ax+by\,;$

faisons $b=\varphi (a),$ on aura

$z=f\left[a,\varphi (a)\right]+ax+\varphi (a)y,$

et faisant varier $a$ seul, on aura, après avoir divisé par $da,$

$f'\left[a,\varphi (a)\right]+x+\varphi '(a)y=0\,;$

et il n’y aura plus qu’à éliminer $a$ au moyen de ces deux équations.

Si l’équation était par exemple

$z=x{\frac {dz}{dx}}+y{\frac {dz}{dy}}+h{\sqrt {1+\left({\frac {dz}{dx}}\right)^{2}+\left({\frac {dz}{dy}}\right)^{2}}},$

comme dans le Problème du n^o 44, on aurait

$f(a,b)=h{\sqrt {1+a^{2}+b^{2}}}\,;$

faisant donc $b=\varphi (a),$ on aurait ces deux équations

${\begin{aligned}&z=ax+\varphi (a)y+h{\sqrt {1+a^{2}+\left[\varphi (a)\right]^{2}}},\\&x+\varphi '(a)+h{\frac {a+\varphi (a)\varphi '(a)}{\sqrt {1+a^{2}+\left[\varphi (a)\right]^{2}}}}=0,\end{aligned}}$

d’où il faudrait éliminer $a\,;$ on aurait ainsi la solution générale du Problème dont il s’agit.

Cette élimination est impossible, en général, c’est-à-dire tant que la fonction $\varphi (a)$ est indéterminée ; ainsi nous nous contenterons d’examiner quelques cas particuliers.

Supposons, ce qui est le cas le plus simple,

$\varphi (a)=m+na,$

$m$ et $n$ étant des coefficients constants quelconques, on aura

$\varphi '(a)=n,$

et les deux équations précédentes deviendront

${\begin{aligned}&z=a(x+ny)+my+h{\sqrt {1+a^{2}+(m+na)^{2}}},\\&x+ny+h{\frac {\left(1+n^{2}\right)a+mn}{\sqrt {1+a^{2}+(m+na)^{2}}}}=0.\end{aligned}}$

Pour chasser $a$ de ces deux équations, je commence par tirer de la seconde la valeur du radical, j’ai

${\sqrt {1+a^{2}+(m+na)^{2}}}=-h{\frac {\left(1+n^{2}\right)a+mn}{x+ny}},$

ce qui, étant substitué dans la première, donne

$z=a\left[x+ny-{\frac {\left(1+n^{2}\right)h^{2}}{x+ny}}\right]+my-{\frac {mnh^{2}}{x+ny}}.$

Maintenant je carre l’équation précédente et je la réduis à celle-ci

$\left(1+n^{2}\right)a^{2}+2mna+{\frac {m^{2}n^{2}h^{2}-\left(1+m^{2}\right)(x+ny)^{2}}{\left(1+n^{2}\right)h^{2}-(x+ny)^{2}}}=0,$

d’où je tire

$a=-{\frac {mn}{1+n^{2}}}+{\frac {\sqrt {1+m^{2}+n^{2}}}{1+n^{2}}}{\frac {x+ny}{\sqrt {\left(1+n^{2}\right)h^{2}-(x+ny)^{2}}}}\,;$

cette valeur étant enfin substituée dans l’équation ci-dessus, on aura celle-ci

$z={\frac {m(y-nx)}{1+n^{2}}}\pm {\frac {\sqrt {1+m^{2}+n^{2}}}{1+n^{2}}}{\sqrt {\left(1+n^{2}\right)h^{2}-(x+ny)^{2}}}.$

Cette équation est celle d’un cylindre droit qui a le rayon de la base égal à $h\,;$ en effet, si l’on change les deux coordonnées rectangles $x,y$ en deux autres aussi rectangles $x',y',$ telles, que

${\frac {x+ny}{\sqrt {1+n^{2}}}}=x',\quad {\frac {y-nx}{\sqrt {1+n^{2}}}}=y',$

on aura

$z={\frac {my'}{\sqrt {1+n^{2}}}}\pm {\frac {\sqrt {1+m^{2}+n^{2}}}{1+n^{2}}}{\sqrt {h^{2}-x'^{2}}},$

et si l’on change encore les deux coordonnées $z,y'$ en deux autres $z',y''$ telles, que

${\frac {z{\sqrt {1+n^{2}}}-my'}{\sqrt {1+m^{2}+n^{2}}}}=z',\quad {\frac {y'{\sqrt {1+n^{2}}}-mz}{\sqrt {1+m^{2}+n^{2}}}}=y'',$

on aura cette équation-ci

$z'={\sqrt {h^{2}-x'^{2}}},$

laquelle est évidemment celle d’un cylindre droit dont l’axe coïncide avec l’axe des coordonnées $y'',$ et dont la base a le rayon égal à $h.$

Or, comme en changeant les coordonnées nous n’avons fait que changer la position du cylindre relativement aux coordonnées primitives $x,y,z,$ il s’ensuit que tout cylindre droit dont l’axe passera par le point donné qui a été pris pour l’origine des coordonnées et dont la base aura la quantité $h$ pour rayon, satisfera au Problème du n^o 44 ; ce qui est évident par soi-même.

Supposons

${\sqrt {1+a^{2}+\left[\varphi (a)\right]^{2}}}=k,$

ce qui donne

$\varphi (a)={\sqrt {k^{2}-1-a^{2}}},$

on aura alors ces deux équations

${\begin{aligned}z&=ax+y{\sqrt {k^{2}-1-a^{2}}}+hk,\\x&-{\frac {ay}{\sqrt {k^{2}-1-a^{2}}}}=0,\end{aligned}}$

d’où, éliminant $a,$ il viendra

$z=hk+{\sqrt {k^{2}-1^{2}}}{\sqrt {x^{2}+y^{2}}},$

équation à un cône droit dont l’axe coïncide avec l’axe des ordonnées $z,$ le sommet est distant du plan des ordonnées $x,y$ de la quantité $hk,$ et la base prise dans ce plan est un cercle dont le rayon est ${\frac {hk}{\sqrt {k^{2}-1^{2}}}}\,;$ de sorte que si du centre de la base on mène une perpendiculaire à la surface du

cône, cette perpendiculaire sera égale à $h\,;$ ce qui est la condition du Problème.

On trouverait un cône semblable, mais dans une situation oblique à l’axe des ordonnées $z,$ si l’on prenait la quantité $\varphi (a)$ telle, que

${\sqrt {1+a^{2}+\left[\varphi (a)\right]^{2}}}=k+ma+n\varphi (a)\,;$

le calcul en étant un peu long, nous ne nous y arrêterons pas.

Comme la forme de la fonction $\varphi (a)$ est arbitraire, on pourra trouver une infinité de solutions différentes du Problème proposé ;` mais il est remarquable que la solution qui donne une sphère, et qui est en quelque façon la plus simple, n’est point comprise dans cette infinité de solutions qui résultent de l’hypothèse que $b$ est une fonction quelconque de $a.$

En effet, si l’on reprend les équations générales et qu’on cherche à déterminer $\varphi (a)$ en sorte qu’il en résulte l’équation à la sphère

$h^{2}=z^{2}+y^{2}+x^{2},$

il faudra qu’en substituant pour $z$ sa valeur ${\sqrt {h^{2}-x^{2}-y^{2}}}$ et éliminant ensuite l’une des deux variables $x$ ou $y$ l’autre disparaisse aussi, de manière que l’équation résultante soit uniquement entre les quantités $a$ et $\varphi (a)\,;$ or en éliminant par exemple $x,$ on aura une équation où $y$ montera au second degré, en sorte qu’il faudrait faire évanouir séparément les coefficients des trois termes de cette équation ordonnée par rapport à $y\,;$ ce qui donnerait trois équations au lieu d’une ; et l’on se convaincra aisément par le calcul qu’il est impossible de satisfaire à ces trois équations à la fois.

Au reste ce résultat ne doit pas paraître surprenantquand on considère que la sphère est donnée par une intégrale particulière (44) ; et nous n’avons fait cette remarque, qui peut d’ailleurs être appliquée à toutes les intégrales particulières des équations à différences partielles, que pour montrer la nécessité d’avoir égard à ces sortes d’intégrales pour avoir toutes les solutions possibles des équations de l’espèce dont il s’agit.

Une autre chose digne d’être remarquée, c’est que les surfaces, qui représentent l’intégrale particulière et l’intégrale générale d’une même équation à différences partielles du premier ordre, touchent dans chaque point une des surfaces exprimées par l’intégrale complète de la même équation, mais avec cette différence que la surface représentée par l’intégrale particulière touche absolument toutes les surfaces possibles que donne l’intégrale complète (43), au lieu que la surface représentée par l’intégrale générale ne touche que celles de ces surfaces qui se rapportent à une certaine espèce dépendante du rapport qu’on établit entre les deux quantités $a$ et $b,$ qui sont les constantes arbitraires de l’intégrale complète ; c’est de quoi on peut se convaincre, en général, par les principes de la méthode du n^o 48, et dont l’Exemple précédent fournit des preuves particulières, puisqu’il est visible que les cylindres et les cônes, que nous avons déduits de l’intégrale générale, sont touchés partout par desplans exprimés par l’intégrale complète.

50. Soit l’équation

$f\left({\frac {dz}{dx}},x\right)=\mathrm {F} \left({\frac {dz}{dy}},y\right),$

les caractéristiques $f$ et $\mathrm {F}$ dénotant des fonctions quelconques données de deux quantités.

Je suppose

$f\left({\frac {dz}{dx}},x\right)=a,$

$a$ étant une constantes ; je tire de cette équation la valeur de ${\frac {dz}{dx}},$ qui sera exprimée en $x$ et $a\,;$ et après avoir multiplié par $dx$ j’intègre en faisant varier $x$ seul ; j’aurai

$z=\mathrm {X} +\Psi ,$

$\mathrm {X}$ étant une fonction connue de $x$ où entrera aussi la constante $a,$ et $\Psi$ une fonction quelconque indéterminée de $y\,;$ j’aurai ensuite

$\mathrm {F} \left({\frac {dz}{dy}},y\right)=a,$

et j’en tirerai de même

$z=\mathrm {Y} +\Xi ,$

$\mathrm {Y}$ étant une fonction connue de $y$ où entrera aussi $a$ comme constante, et $\Xi$ une fonction quelconque de $x\,;$ donc, puisque ces deux valeurs de $z$ doivent être identiques, il faudra que

$\Psi =\mathrm {Y} \quad {\text{et}}\quad \Xi =\mathrm {X} \,;$

par conséquent la valeur de $z$ sera $\mathrm {X+Y} ,$ et il est visible qu’on peut ajouter à cette valeur une constante quelconque, puisque dans l’équation différentielle la quantité finie $z$ ne se trouve pas.

On aura donc

$z=\mathrm {X+Y} +b,$

intégrale complète de la proposée, puisqu’elle renferme deux constantes arbitraires $a$ et $b.$

Pour en tirer l’intégrale générale, on fera $b=\varphi (a),$ ensuite on différentiera en faisant varier $a$ seul ; on aura ainsi les deux équations

$z=\mathrm {X+Y} +\varphi (a),\quad {\frac {d\mathrm {X} }{da}}+{\frac {d\mathrm {Y} }{da}}+\varphi '(a)=0,$

au moyen desquelles on éliminera $a,$ et la résultante sera l’intégrale cherchée.

51. Soit l’équation

$f\left({\frac {dz}{dx}},{\frac {dz}{dy}},z\right)=0.$

Je fais

${\frac {dz}{dx}}=\mathrm {Z} ,\quad {\frac {dz}{dy}}=a\mathrm {Z} \,;$

substituant ces valeurs, j’aurai une équation en $z$ et $\mathrm {Z} ,$ d’où je tirerai $\mathrm {Z}$ exprimé par une fonction de $z$ seul dans laquelle $a$ entrera aussi comme constante ; or l’équation

${\frac {dz}{dx}}=\mathrm {Z}$

donnera, en intégrant,

$x=\int {\frac {dz}{\mathrm {Z} }}+\mathrm {Y} ,$

$\mathrm {Y}$ étant une fonction de $y$ seul ; de même l’équation

${\frac {dz}{dy}}=a\mathrm {Z}$

donnera

$ay=\int {\frac {dz}{\mathrm {Z} }}+\mathrm {X} ,$

$\mathrm {X}$ étant une fonction de $x$ seul ; donc, pour que ces deux équations deviennent la même chose, il faudra que l’on ait

$\mathrm {X} =-x,\quad \mathrm {Y} =-ay,$

et alors on aura, en ajoutant une constante $b$ à l’intégrale $\int {\frac {dz}{\mathrm {Z} }},$ cette intégrale complète

$x+ay=\int {\frac {dz}{\mathrm {Z} }}+b.$

Je fais maintenant $b=\varphi (a),$ et je différentie en faisant varier $a$ seul j’aurai les deux équations

$x+ay=\int {\frac {dz}{\mathrm {Z} }}+\varphi (a),\quad y={\frac {d\int {\cfrac {dz}{\mathrm {Z} }}}{da}}+\varphi '(a),$

au moyen desquelles il faudra éliminer $a$ pour avoir l’intégrale générale de la proposée.

52. Soit l’équation

${\frac {dz}{dx}}=\mathrm {V} {\frac {dz}{dy}}+\mathrm {Z} ,$

$\mathrm {V}$ étant une fonction quelconque de $x$ et $y,$ et $\mathrm {Z}$ une fonction quelconque de $x,y,z.$

Je multiplie l’équation par $dx,$ et j’ajoute ensuite à l’un et à l’autre membre la quantité ${\frac {dz}{dy}}dy,$ j’ai à cause de

$dz={\frac {dz}{dx}}dx+{\frac {dz}{dy}}dy,$

cette équation-ci

$dz=(\mathrm {V} dx+dy){\frac {dz}{dy}}+\mathrm {Z} dx\,;$

je suppose pour un moment

$\mathrm {V} dx+dy=0,$

j’ai une équation entre deux variables $x$ et $y,$ que j’intègre en $y$ ajoutant la constante arbitraire $\alpha \,;$ je regarde maintenant $\alpha$ comme une fonction de $x$ et $y$ déterminée par cette même équation ; j’aurai par la différentiation

$\mathrm {V} dx+dy=\mathrm {A} d\alpha ,$

$\mathrm {A}$ étant une fonction connue de $x,y$ et $\alpha \,;$ ainsi substituant cette valeur dans l’équation précédente, elle deviendra

$dz=\mathrm {A} {\frac {dz}{dy}}d\alpha +\mathrm {Z} dx.$

Or, si l’on substitue partout dans cette équation à la place de $y$ sa valeur en $x$ et $\alpha ,$ on aura une équation entre les trois variables $x,z,\alpha \,;$ et, supposant $\alpha$ constante, on aura l’équation

$dz=\mathrm {Z} dx$

entre les deux variables $x$ et $z,$ laquelle étant intégrée en y ajoutant une constante arbitraire, qui pourra être une fonction quelconque indéterminée de $\alpha ,$ donnera sur-le-champ l’intégrale générale de la proposée ; car il n’y aura plus qu’à y remettre à la place de $\alpha$ sa valeur en $x$ et $y.$

53. Soit encore l’équation

$x=\mathrm {V} y+\mathrm {Z} ,$

$\mathrm {V}$ étant une fonction quelconque de ${\frac {dz}{dx}}$ et ${\frac {dz}{dy}},$ et $\mathrm {Z}$ une fonction quelconque de ${\frac {dz}{dx}},{\frac {dz}{dy}},z-x{\frac {dz}{dx}}-y{\frac {dz}{dy}}.$

Je fais pour plus de simplicité

${\frac {dz}{dx}}=p,\quad {\frac {dz}{dy}}=q,\quad z-x{\frac {dz}{dx}}-y{\frac {dz}{dy}}=r,$

en sorte que $\mathrm {V}$ soit une fonction quelconque de $p$ et $q,$ et $\mathrm {Z}$ une fonction quelconque de $p,q,r\,;$ je multiplie toute l’équation par $dp$ et j’ajoute aux deux membres la quantité $ydq,$ j’ai

$xdp+ydq=y(\mathrm {V} dp+dq)+\mathrm {Z} dp\,;$

or, puisque

$r=z-xp-yq,$

et que

$dz={\frac {dz}{dx}}dx+{\frac {dz}{dy}}dy=pdx+qdy,$

on aura

$dr=-xdp-ydy\,;$

donc l’équation précédente deviendra

$-dr=y(\mathrm {V} dp+dq)+\mathrm {Z} dp.$

Je suppose

$\mathrm {V} dp+dq=0,$

ce qui fait une équation entre $p$ et $q$ que j’intègre en y ajoutant une constante arbitraire $\alpha \,;$ de sorte que j’ai une équation finie entre $p,q$ et $\alpha ,$ dans laquelle je puis regarder $\alpha$ comme une fonction de $p$ et $q,$ et qui donnera par la différentiation

$\mathrm {V} dp+dq=\mathrm {A} d\alpha ,$

$\mathrm {A}$ étant une fonction connue de $p,q$ et $\alpha \,;$ ainsi l’équation précédente deviendra

$-dr=y\mathrm {A} d\alpha +\mathrm {Z} dp.$

Qu’on substitue partout dans cette équation à la place de $q$ sa valeur en $p$ et $\alpha ,$ on aura une équation entre les trois variables $r,p,\alpha ,$ laquelle, en supposant $\alpha$ constante, sera

$-dr=\mathrm {Z} dp\,;$

qu’on intègre donc cette équation entre les deux variables $r$ et $p,$ et soit $\mathrm {R} =0$ l’intégrale, dans laquelle on pourra supposer que la constante arbitraire soit une fonction quelconque indéterminée de $\alpha \,;$ faisant ensuite varier dans l’équation $\mathrm {R} =0$ les trois quantités $r,p$ et $\alpha$ à la fois, il viendra

$dr=-\mathrm {Z} dp+\mathrm {Q} d\alpha ,$

$\mathrm {Q}$ étant une fonction connue de $r,p,\alpha \,;$ donc, substituant cette valeur dans l’équation ci-dessus, on aura

$-\mathrm {Q} d\alpha =y\mathrm {A} d\alpha ,\quad {\text{d’où}}\quad \mathrm {A} y+\mathrm {Q} =0.$

On a donc ainsi deux équations finies

$\mathrm {R} =0,\quad \mathrm {A} y+\mathrm {Q} =0$

entre les quantités $p,q,\alpha ,y\,;$ et, comme $\alpha$ est donnée par une fonction connue de $p$ et $q,$ substituant cette valeur de $\alpha ,$ on aura deux équations entre les trois quantités $p,q$ et $y,$ à l’aide desquelles on pourra éliminer $p$ et $q,$ c’est-à-dire ${\frac {dz}{dx}}$ et ${\frac {dz}{dy}},$ dans l’équation proposée

$x=\mathrm {V} y+\mathrm {Z} \,;$

et l’équation résultante sera l’intégrale générale de cette même équation, puisqu’elle contient déjà une fonction indéterminée.

54. Soit enfin l’équation

${\frac {dz}{dx}}=\mathrm {V} y+\mathrm {Z} ,$

$\mathrm {V}$ étant une fonction quelconque de $x$ et ${\frac {dz}{dy}},$ et $\mathrm {Z}$ une fonction quelconque de $x,{\frac {dz}{dy}},z-y{\frac {dz}{dy}}.$

Je fais

${\frac {dz}{dx}}=p,\quad {\frac {dz}{dy}}=q,\quad z-yq=r,$

j’aurai l’équation

$p=\mathrm {V} y+\mathrm {Z} ,$

dans laquelle $\mathrm {V}$ sera une fonction quelconque de $x,q,$ et où $\mathrm {Z}$ sera une fonction quelconque de $x,q,r.$ Je multiplie cette équation par $dx,$ et j’ajoute aux deux membres la quantité $-ydq,$ j’aurai, à cause de

$dr=dz-d(yq)=pdx+qdy-d(yq)=pdx-ydq,$

j’aurai, dis-je, l’équation

$dr=y(\mathrm {V} dx-dq)+\mathrm {Z} dx.$

Je suppose $\mathrm {V} dx-dq=0,$ j’ai une équation entre $x$ et $q$ que j’intègre en y ajoutant une constante arbitraire $\alpha \,;$ ensuite, regardant $\alpha$ comme variable, je différentie de nouveau, j’ai

$\mathrm {V} dx-dq=\mathrm {A} d\alpha ,$

$\mathrm {A}$ étant une fonction connue de $x,q,\alpha \,;$ cette substitution ainsi que celle de la valeur de $q$ en $x$ et $\alpha$ étant faites dans l’équation ci-dessus, elle deviendra

$dr=y\mathrm {A} d\alpha +\mathrm {Z} dx,$

où $\mathrm {A}$ et $\mathrm {Z}$ seront maintenant des fonctions connues de $r,x$ et $\alpha .$ Regardant donc $\alpha$ comme constante, on aura l’équation

$dr=\mathrm {Z} dx$

entre les variables $r$ et $x,$ dont l’intégrale pourra contenir comme constante une fonction quelconque indéterminée de $\alpha .$ Soit

$\mathrm {R} =0$

cette intégrale ; en y faisant varier à la fois $r,x$ et $\alpha ,$ on en tirera

$dr=\mathrm {Z} dx+\mathrm {Q} d\alpha \,;$

mais on a

$dr=y\mathrm {A} d\alpha +\mathrm {Z} dx\,;$

donc

$\mathrm {Q} =y\mathrm {A} .$

Cette équation étant combinée avec l’équation

$\mathrm {R} =0,$

on pourra éliminer $q,$ après avoir remis pour $r$ sa valeur $z-yq,$ et pour $\alpha$ sa valeur en $x$ et $q\,;$ l’équation résultante ne contiendra que $x,y,z,$ et sera l’intégrale générale de la proposée à cause de la fonction indéterminée qui s’y trouvera.

On pourra, par un procédé semblable, trouver l’intégrale générale de toute équation de la forme

${\frac {dz}{dy}}=\mathrm {V} x+\mathrm {Z} ,$

$\mathrm {V}$ étant une fonction de $y,{\frac {dz}{dx}},$ et $\mathrm {Z}$ une fonction de $y,{\frac {dz}{dx}},z-x{\frac {dz}{dx}}.$

Il n’y aura pour cela qu’à changer dans l’analyse précédente $x$ en $y$ et vice versâ.

On voit par là que les équations de la forme

$z=y{\frac {dz}{dy}}+f\left(x,{\frac {dz}{dx}},{\frac {dz}{dy}}\right),$

ou

$z=x{\frac {dz}{dx}}+f\left(y,{\frac {dz}{dy}},{\frac {dz}{dx}}\right)$

(la caractéristique $f$ dénotant une fonction quelconque donnée de trois quantités), sont intégrables en général car si l’on tire de la première la

valeur de ${\frac {dz}{dx}},$ ou de la seconde celle de ${\frac {dz}{dy}},$ on aura des équations de la forme

${\frac {dz}{dx}}=\mathrm {Z} ,\quad {\text{ou}}\quad {\frac {dz}{dy}}=\mathrm {Z} ,$

et par conséquent intégrables par la méthode ci-dessus.

55. Les Exemples précédents renferment d’une manière générale tous les cas connus d’intégration des équations de différences partielles du premier ordre entre trois variables ; et c’est pour cette raison que nous avons ajouté les trois derniers Exemples, quoique la méthode qu’on y a suivie n’ait pas un rapport immédiat avec la méthode générale du n^o  48 ; nous avions déjà donné ailleurs l’intégration de l’équation du n^o  52 [voyez les Mémoires pour 1772^[3]] ; mais celle des équations des n^os 53 et 54 est, si je ne me trompe, entièrement nouvelle.

56. On a vu ci-dessus que l’intégrale particulière n’est renfermée ni dans l’intégrale complète ni dans l’intégrale générale ; mais il n’en est pas de même de l’intégrale complète par rapport à l’intégrale générale ; car il est facile de se convaincre, soit d’après notre théorie de la formation des intégrales générales, soit d’après la seule considération de la nature de ces intégrales, laquelle consiste en ce qu’elles doivent renfermer une fonction arbitraire, il est aisé, dis-je, de se convaincre que ces sortes d’intégrales doivent toujours renfermer les intégrales complètes comme des cas particuliers. En effet, si dans l’intégrale générale on donne à la fonction indéterminée une valeur particulière dans laquelle il y ait des coefficients arbitraires, en sorte qu’il se trouve deux de ces coefficients dans l’intégrale, cette intégrale sera alors une intégrale complète, et conduira nécessairement par la différentiation à la même équation aux différences partielles du premier ordre que l’intégrale générale dont elle est dérivée. On voit par là qu’on peut donner différentes formes aux intégrales complètes, mais que ces formes différentes sont néanmoins liées entre elles, en sorte que dès qu’on en connaît une on peut en déduire toutes les autres, puisqu’il n’y a qu’à chercher d’abord, par notre méthode, l’intégrale générale, et ensuite donner à la fonction indéterminée des valeurs particulières.

57. Nous allons maintenant considérer les équations à différences partielles du second ordre ; et nous remarquerons d’abord que si $\mathrm {V} =0$ est une équation finie entre les trois variables $x,y,z$ et cinq constantes $a,b,c,h,g,$ on pourra en déduire une équation aux différences partielles du second ordre dans laquelle ces constantes ne se trouveront plus. Car, en regardant $z$ comme une fonction de $x$ et $y$ et faisant varier successivement ces deux quantités dans l’équation donnée, on en tirera par une double différentiation ces cinq équations-ci

${\frac {dz}{dx}}-p=0,\quad {\frac {dz}{dy}}-q=0,$

${\frac {d^{2}z}{dx^{2}}}-r=0,\quad {\frac {d^{2}z}{dxdy}}-s=0,\quad {\frac {d^{2}z}{dy^{2}}}-t=0,$

$p,q$ étant des fonctions connues de $x,y,z$, et $r,s,t$ des fonctions aussi connues de $x,y,z,{\frac {dz}{dx}},{\frac {dz}{dy}}.$ Si donc au moyen de ces cinq équations différentielles et de l’équation finie $\mathrm {V} =0$ on élimine les cinq constantes $a,b,c,$$h,g,$ on aura une équation finale entre les quantités

$x,\ \ y,\ \ z,\ \ {\frac {dz}{dx}},\ \ {\frac {dz}{dy}},\ \ {\frac {d^{2}z}{dx^{2}}},\ \ {\frac {d^{2}z}{dxdy}},\ \ {\frac {d^{2}z}{dy^{2}}},$

dans laquelle les constantes dont il s’agit ne se trouveront plus, et que nous représenterons, en général, par $\mathrm {Z} '=0.$

On pourra donc regarder l’équation $\mathrm {V} =0$ comme l’intégrale finie et complète de l’équation aux différences partielles du second ordre $\mathrm {Z} '=0,$ et comme les cinq constantes $a,b,c,g,h$ demeurent arbitraires, il s’ensuit que l’intégrale finie et complète de toute équation aux différences partielles du second ordre doit renfermer cinq constantes arbitraires.

Soit par exemple l’équation finie

$z=a+bx+cy+hx^{2}+gxy+mhy^{2}\,;$

on aura par la différentiation

${\frac {dz}{dx}}=b+2hx+gy,\quad {\frac {dz}{dy}}=c+gx+2mhy,$

${\frac {d^{2}z}{dx^{2}}}=2h,\quad {\frac {d^{2}z}{dxdy}}=g,\quad {\frac {d^{2}z}{dy^{2}}}=2mh\,;$

donc, éliminant les constantes $a,b,c,h,g,$ on aura cette équation finale

${\frac {d^{2}z}{dy^{2}}}=m{\frac {d^{2}z}{dx^{2}}},$

dont l’équation précédente sera par conséquent l’intégrale complète.

58. Soit encore l’équation

$z=k+ax+by+cx^{2}+gxy+hy^{2},$

dans laquelle $k$ soit une fonction donnée de $a,b,c,g,h\,;$ on aura par la différentiation

${\frac {dz}{dx}}=a+2cx+gy,\quad {\frac {dz}{dy}}=b+2hy+gx,$

${\frac {d^{2}z}{dx^{2}}}=2c,\quad {\frac {d^{2}z}{dxdy}}=g,\quad {\frac {d^{2}z}{dy^{2}}}=2h\,;$

donc

$h={\frac {1}{2}}{\frac {d^{2}z}{dy^{2}}},\quad g={\frac {d^{2}z}{dxdy}},\quad c={\frac {1}{2}}{\frac {d^{2}z}{dx^{2}}},$

$a={\frac {dz}{dx}}-x{\frac {d^{2}z}{dx^{2}}}-y{\frac {d^{2}z}{dxdy}},\quad b={\frac {dz}{dy}}-y{\frac {d^{2}z}{dy^{2}}}-x{\frac {d^{2}z}{dxdy}}\,;$

donc, substituant ces valeurs, on aura l’équation différentielle du second ordre

${\begin{aligned}&z=x{\frac {dz}{dx}}+y{\frac {dz}{dy}}-{\frac {x^{2}}{2}}{\frac {d^{2}z}{dx^{2}}}-xy{\frac {d^{2}z}{dxdy}}-{\frac {y^{2}}{2}}{\frac {d^{2}z}{dy^{2}}}\\&+f\left({\frac {dz}{dx}}-x{\frac {d^{2}z}{dx^{2}}}-y{\frac {d^{2}z}{dxdy}},\ {\frac {dz}{dy}}-y{\frac {d^{2}z}{dy^{2}}}-x{\frac {d^{2}z}{dxdy}},\ {\frac {1}{2}}{\frac {d^{2}z}{dx^{2}}},\ {\frac {d^{2}z}{dxdy}},\ {\frac {1}{2}}{\frac {d^{2}z}{dy^{2}}},\right),\end{aligned}}$

dont l’intégrale complète sera

$z=ax+by+cx^{2}+gxy+hy^{2}+f(a,b,c,g,h),$

la caractéristique $f$ dénotant une fonction quelconque de cinq quantités.

59. En suivant les principes établis dans ce Mémoire, il est facile de démontrer que pour avoir l’intégrale particulière d’une équation à différences partielles du second ordre $\mathrm {Z} '=0,$ dont on connaît l’intégrale finie complète $\mathrm {V} =0,$ il n’y aura qu’à faire varier, tant dans la quantité $\mathrm {V}$ que dans les deux quantités $p$ et $q$ des équations

${\frac {dz}{dx}}-p=0,\quad {\frac {dz}{dy}}-q=0,$

les cinq constantes arbitraires $a,b,c,g,h,$ et supposer les différences de ces trois quantités $\mathrm {V} ,p,q$ nulles ; ce qui donnera trois équations qui contiendront les cinq différentielles $da,db,dc,dg,dh$ sous une forme linéaire ; on éliminera deux de ces différentielles, et l’on fera ensuite évanouir séparément dans l’équation résultante les coefficients des trois différentielles restantes ; on aura ainsi trois équations de condition qui étant combinées avec les trois équations

$\mathrm {V} =0,\quad {\frac {dz}{dx}}-p=0,\quad {\frac {dz}{dy}}-q=0$

donneront, par l’élimination des cinq quantités $a,b,c,g,h,$ une équation à différences partielles du premier ordre, laquelle sera l’intégrale particulière aux premières différences de la proposée $\mathrm {Z} '=0$ du second ordre.

Par exemple, l’équation différentio-différentielle du numéro précédent donne

$\mathrm {V} =-z+ax+by+cx^{2}+gxy+hy^{2}+f(a,b,c,g,h),$

$p=a+2cx+gy,\quad q=b+2hy+gx\,;$

donc les équations de condition de l’intégrale particulière seront, en fai-

sant, pour abréger, $f(a,b,c,g,h)=\mathrm {F} ,$

$\left(x+{\frac {d\mathrm {F} }{da}}\right)da+\left(y+{\frac {d\mathrm {F} }{db}}\right)db$

$+\left(x^{2}+{\frac {d\mathrm {F} }{dc}}\right)dc+\left(xy+{\frac {d\mathrm {F} }{dg}}\right)dg+\left(y^{2}+{\frac {d\mathrm {F} }{dh}}\right)dh=0.$

$da+2xdc+ydg=0,\quad db+2ydh+xdg=0.$

Mettant dans la première équation les valeurs de $da$ et $db$ tirées de ces deux dernières, et égalant ensuite à zéro les coefficients des différences $dc,dh,dg,$ on aura ces trois équations

${\frac {d\mathrm {F} }{dc}}-2x{\frac {d\mathrm {F} }{da}}-x^{2}=0,\quad {\frac {d\mathrm {F} }{dh}}-2y{\frac {d\mathrm {F} }{db}}-y^{2}=0,$

${\frac {d\mathrm {F} }{dg}}-x{\frac {d\mathrm {F} }{db}}-y{\frac {d\mathrm {F} }{da}}-xy=0,$

à l’aide desquelles et des équations

$\mathrm {V} =0,\quad {\frac {dz}{dx}}-p=0,\quad {\frac {dz}{dy}}-q=0$

on éliminera les cinq quantités $a,b,c,g,h\,;$ et la résultante sera l’intégrale particulière aux premières différences de l’équation dont il s’agit.

60. Soit à présent $\mathrm {Z} =0$ une équation entre les variables finies $x,y,z$ et les différences partielles du premier ordre ${\frac {dz}{dx}},{\frac {dz}{dy}}$ dans laquelle entrent deux constantes arbitraires $a$ et $b\,;$ on pourra en déduire par la différentiation une équation aux différences partielles du second ordre dans laquelle les deux constantes $a$ et $b$ ne se trouveront plus. Car, si l’on fait varier successivement $x$ et $y,$ on aura les deux équations

${\frac {d\mathrm {Z} }{dx}}=0,\quad {\frac {d\mathrm {Z} }{dy}}=0,$

à l’aide desquelles on pourra éliminer dans la proposée $\mathrm {Z} =0$ les deux constantes $a$ et $b\,;$ et l’équation résultante sera entre les variables $x,y,z$ et leurs différences partielles ${\frac {dz}{dx}},\ {\frac {dz}{dy}},\ {\frac {d^{2}z}{dx^{2}}},\ {\frac {d^{2}z}{dxdy}},\ {\frac {d^{2}z}{dy^{2}}}\,;$ nous la designerons, en général, par $\mathrm {Z} '=0.$

Ainsi l’on pourra regarder l’équation $\mathrm {Z} =0$ comme l’intégrale complète du premier ordre de l’équation du second ordre $\mathrm {Z} '=0,$ et, comme les deux constantes $a$ et $b$ restent arbitraires, on peut conclure, en général, que toute intégrale complète du premier ordre d’une équation aux différences partielles du second ordre doit contenir deux constantes arbitraires.

Soit par exemple l’équation aux différences premières

${\frac {dz}{dy}}-m{\frac {dz}{dx}}=a+bx+nby,$

$a$ et $b$ étant deux constantes qui doivent rester indéterminées, et $m,n$ deux coefficients donnés.

Faisant varier successivement $x$ et $y,$ on aura les deux équations

${\frac {d^{2}z}{dy^{2}}}-m{\frac {d^{2}z}{dxdy}}=nb,\quad {\frac {d^{2}z}{dxdy}}-m{\frac {d^{2}z}{dx^{2}}}=b\,;$

d’où chassant $b,$ il viendra l’équation

${\frac {d^{2}z}{dy^{2}}}-(m+n){\frac {d^{2}z}{dxdy}}+mn{\frac {d^{2}z}{dx^{2}}}=0,$

dont l’équation ci-dessus sera par conséquent l’intégrale complète du premier ordre.

Ainsi, ayant une équation du second ordre de la forme

${\frac {d^{2}z}{dy^{2}}}-\mathrm {A} {\frac {d^{2}z}{dxdy}}+\mathrm {B} {\frac {d^{2}z}{dx^{2}}}=0,$

il n’y aura qu’à chercher les racines de l’équation

$s^{2}-\mathrm {A} s+\mathrm {B} =0,$

et nommant ces racines $m,n,$ on aura l’intégrale complète du premier ordre

${\frac {dz}{dy}}-m{\frac {dz}{dx}}=a+b(x+ny)\,;$

et, comme on peut échanger entre elles les racines $m$ et $n,$ on aura aussi cette autre intégrale du premier ordre

${\frac {dz}{dy}}-n{\frac {dz}{dx}}=h+g(x+my),$

$h$ et $g$ étant d’autres constantes arbitraires.

Au moyen de ces deux intégrales du premier ordre on pourra, si l’on veut, trouver l’intégrale complète finie ; car en chassant par exemple ${\frac {dz}{dy}},$ on aura l’équation

$(n-m){\frac {dz}{dx}}=a-h+(b-g)x+(nb-mg)y,$

qui peut être intégrée en faisant varier $x$ seul ; et l’on aura

$(n-m)z=\mathrm {Y} +(a-h)x+(b-g){\frac {x^{2}}{2}}+(nb-mg)xy,$

$\mathrm {Y}$ étant une fonction quelconque de $y.$

Or, si des mêmes équations on chasse ${\frac {dz}{dx}},$ on a

$(n-m){\frac {dz}{dy}}=na-mh+(nb-mg)x+\left(n^{2}b-m^{2}g\right)y,$

et, intégrant par rapport à $y$ seul, on aura

$(n-m)z=\mathrm {X} +(na-mh)y+(nb-mg)xy+\left(n^{2}b-m^{2}g\right){\frac {y^{2}}{2}},$

et, comme ces deux équations doivent être la même chose, il faudra faire

${\begin{aligned}\mathrm {Y} =&c+(na-mh)y+\left(n^{2}b-m^{2}g\right){\frac {y^{2}}{2}},\\\mathrm {X} =&c+(a-h)x+(b-g){\frac {x^{2}}{2}}\,;\end{aligned}}$

$c$ étant une constante arbitraire.

Ainsi l’intégrale complète finie de la proposée sera

${\begin{aligned}(n-m)z=c&+(a-h)x+(na-mh)y+(b-g){\frac {x^{2}}{2}}\\&+(nb-mg)xy+(n^{2}b-m^{2}g){\frac {y^{2}}{2}},\end{aligned}}$

où $a,b,c,g,h$ sont les cinq arbitraires.

Au reste il est facile de voir que cette méthode de trouver l’intégrale complète finie d’une équation du second ordre, lorsqu’on connaît deux différentes intégrales complètes du premier ordre de la même équation, est générale et réussira toujours, quelle que soit la forme de ces intégrales complètes ; car en éliminant ${\frac {dz}{dx}}$ on aura toujours une équation où $y$ pourra être traitée comme constante, et en éliminant ${\frac {dz}{dy}}$ on en aura une où $x$ pourra être regardée elle-même comme constante ; et, l’intégration introduisant nécessairement une nouvelle constante arbitraire, on aura dans l’intégrale finie le nombre de cinq arbitraires ; ce qui est le caractère des intégrales complètes (57).

61. Lorsqu’on connaît l’intégrale complète du premier ordre $\mathrm {Z} =0$ d’une équation du second ordre $\mathrm {Z} '=0,$ on en peut déduire sans peine son intégrale particulière.

Car il ne faudra que faire varier dans l’équation $\mathrm {Z} =0$ les deux constantes arbitraires $a$ et $b,$ et supposer ensuite les coefficients des deux différentielles $da$ et $db$ chacun égal à zéro ; on aura ainsi deux équations qui serviront à éliminer les quantités $a$ et $b$ dans l’équation $\mathrm {Z} =0,$ et la résultante sera l’intégrale particulière aux premières différences de la proposée $\mathrm {Z} '=0$.

62. Enfin, si l’on ne connaît aucune des intégrales complètes de l’équation $\mathrm {Z} '=0$ du second ordre, on pourra néanmoins trouver son intégrale particulière aux premières différences.

Il n’y aura pour cela qu’à différentier l’équation proposée, et, ayant fait disparaître les fractions pour avoir une équation de la forme

$\mathrm {M} d{\frac {d^{2}z}{dx^{2}}}+\mathrm {N} d{\frac {d^{2}z}{dy^{2}}}+\mathrm {P} d{\frac {d^{2}z}{dxdy}}+\mathrm {Q} dx+\mathrm {R} dy=0.$

$\mathrm {M,N,P,Q,R}$ étant des fonctions connues et entières de

$x,\ \ y,\ \ z,\ \ {\frac {dz}{dx}},\ \ {\frac {dz}{dy}},\ \ {\frac {d^{2}z}{dx^{2}}},\ \ {\frac {d^{2}z}{dxdy}},\ \ {\frac {d^{2}z}{dy^{2}}},$

on supposera séparément égale à zéro chacune des cinq quantités $\mathrm {M,N,P,Q,R} ,$ ce qui donnera cinq équations, lesquelles étant combinées avec l’équation $\mathrm {Z} '=0,$ en sorte que les trois quantités ${\frac {d^{2}z}{dx^{2}}},\ {\frac {d^{2}z}{dy^{2}}},\ {\frac {d^{2}z}{dxdy}}$ disparaissent, il résultera trois équations finales en $x,y,z,{\frac {dz}{dx}},{\frac {dz}{dy}}\,;$ il faudra donc que ces trois équations puissent avoir lieu en même temps, c’est-à-dire qu’elles aient un facteur commun, pour que la proposée soit susceptible d’une intégrale particulière ; et ce facteur, s’il y en a un, sera l’intégrale cherchée.

La démonstration de cette méthode et de celle du numéro précédent est facile à déduire des principes exposés dans ce Mémoire, et nous ne croyons pas devoir nous y arrêter.

63. Si l’équation $\mathrm {Z} '=0$ est telle, que

$d\mathrm {Z} '=\mathrm {A} d{\frac {d^{2}z}{dx^{2}}}+\mathrm {B} d{\frac {d^{2}z}{dxdy}}+\mathrm {C} d{\frac {d^{2}z}{dy^{2}}},$

alors on n’aura pour la détermination de l’intégrale particulière que les trois équations de condition

$\mathrm {A=0,\quad B=0,\quad C=0} ,$

lesquelles serviront à éliminer les trois quantités ${\frac {d^{2}z}{dx^{2}}},\ {\frac {d^{2}z}{dxdy}},\ {\frac {d^{2}z}{dy^{2}}}$ dans l’équation $\mathrm {Z} '=0\,;$ et la résultante sera l’intégrale particulière de cette même équation

64. Pour trouver la forme générale des équations qui ont cette propriété, il n’y a qu’à remarquer qu’on peut satisfaire à l’équation $d\mathrm {Z} '=0$ en faisant séparément

$d{\frac {d^{2}z}{dx^{2}}}=0,\quad d{\frac {d^{2}z}{dxdy}}=0,\quad d{\frac {d^{2}z}{dy^{2}}}=0,$

ce qui donne, en prenant des constantes arbitraires,

${\frac {d^{2}z}{dx^{2}}}=2c,\quad {\frac {d^{2}z}{dxdy}}=g,\quad {\frac {d^{2}z}{dy^{2}}}=2h,$

et intégrant de nouveau

${\frac {dz}{dx}}=a+2cx+gy,\quad {\frac {dz}{dy}}=b+2hy+gx,$

et enfin

$z=k+ax+by+cx^{2}+gxy+hy^{2}\,;$

où il se trouve, comme l’on voit, six constantes indéterminées. Or, si l’on substitue ces valeurs dans l’équation $\mathrm {Z} '=0,$ il arrivera nécessairement que les quantités $x$ et $y$ s’en iront d’elles-mêmes, en sorte qu’il ne restera qu’une équation entre les constantes $k,a,b,\ldots ,$ par laquelle il faudra en déterminer une par les autres.

Supposons donc qu’on ait déterminé $k,$ en sorte que l’on ait, en général,

$k=f(a,b,c,g,h),$

alors on aura l’équation finie

$z=ax+by+cx^{2}+gxy+hy^{2}+f(a,b,c,g,h)$

qui, contenant cinq constantes arbitraires, sera l’intégrale complète de la proposée $\mathrm {Z} '=0\,;$ et d’où l’on pourra par conséquent déduire la forme générale de cette même équation, ainsi que nous l’avons déjà fait plus haut (58).

L’équation différentio-différentielle, que nous avons trouvée dans le n^o 58, sera donc la formule générale de touttes les équations qui peuvent avoir la propriété en question ; et si l’on différentie cette équation, qu’ensuite on suppose égal à zéro chacun des coefficients des trois différences $d{\frac {d^{2}z}{dx^{2}}},\ d{\frac {d^{2}z}{dxdy}},\ d{\frac {d^{2}z}{dy^{2}}},$ on aura trois équations, qui étant combinées avec l’équation proposée donneront, par l’élimination des quantités ${\frac {d^{2}z}{dx^{2}}},\ {\frac {d^{2}z}{dy^{2}}},$${\frac {d^{2}z}{dxdy}},$ le même résultat que l’on aura par la méthode du n^o 62.

65. Après avoir vu comment on peut déduire les intégrales particulières des intégrales complètes, voyons comment on peut en déduire aussi les intégrales générales. Et d’abord il est facile de prouver d’après les principes exposés ci-dessus (57) que pour que l’équation $\mathrm {V} =0,$ qu’on suppose être l’intégrale complète finie de l’équation du second ordre $\mathrm {Z} '=0,$ satisfasse à cette même équation, en y regardant les cinq arbitraires $a,b,c,g,h$ comme variables, il suffira que ces quantités soient telles, qu’elles satisfassent aux trois équations différentielles qu’on aura en égalant à zéro les différences des quantités $\mathrm {V} ,p,q,$ dans lesquelles on n’aura fait varier que les quantités $a,b,c,g,h.$

Or, comme de cette manière on n’a que trois équations pour la détermination des cinq variables $a,b,c,g,h,$ il est clair qu’il y en aura deux à volonté qu’on pourra supposer être des fonctions quelconques indéterminées des trois autres ; il s’agira seulement de faire en sorte qu’on puisse déterminer ensuite les valeurs des cinq variables dont il s’agit d’une manière finie, en fonctions de $x,y,z\,;$ alors il n’y aura plus qu’à substituer ces valeurs dans l’équation $\mathrm {V} =0,$ et l’on aura l’intégrale générale de la proposée $\mathrm {Z} '=0$ du second ordre, laquelle intégrale contiendra deux fonctions indéterminées.

Nous avons vu ci-dessus (57) que l’équation du second ordre

${\frac {d^{2}z}{dy^{2}}}=m{\frac {d^{2}z}{dx^{2}}}$

a pour intégrale finie complète

$z=a+bx+cy+hx^{2}+gxy+mhy^{2}.$

De là on a

$\mathrm {V} =-z+a+bx+cy+hx^{2}+gxy+mhy^{2},$

$p={\frac {dz}{dx}}=b+2hx+gy,\quad q={\frac {dz}{dy}}=c+gx+2mhy\,;$

donc, faisant varier les quantités $a,b,c,g,h,$ on a les équations de condition

${\begin{aligned}&da+xdb+ydc+x^{2}dh+xydg+my^{2}dh=0,\\&db+2xdh+ydg=0,\\&dc+xdg+2mydh=0.\end{aligned}}$

J’ajoute ces deux dernières équations ensemble apyès en avoir multiplié une par un coefficient constant arbitraire $\mu ,$ ce qui donne

$\mu db+dc+x(dg+2\mu dh)+\mu y\left(dg+{\frac {2mdh}{\mu }}\right)=0,$

où l’on voit que, si l’on fait $\mu ={\frac {m}{\mu }},$ ce qui donne $\mu ={\sqrt {m}},$ on aura

$dc+db{\sqrt {m}}+\left(x+y{\sqrt {m}}\right)\left(dg+2df{\sqrt {m}}\right)=0\,;$

cette équation ne contenant plus que deux variables $c+b{\sqrt {m}}$ et $g+2h{\sqrt {m}},$ il n’y aura qu’à supposer, à l’imitation de ce que nous avons pratiqué plus haut (47),

$c+b{\sqrt {m}}=\varphi (g+2h{\sqrt {m}}),$

ce qui donnera

$dc+db{\sqrt {m}}=\varphi '\left(g+2h{\sqrt {m}}\right)\left(dg+2dh{\sqrt {m}}\right),$

et l’on aura, après la substitution et la division par $dg+2dh{\sqrt {m}},$

$\varphi '\left(g+2h{\sqrt {m}}\right)+x+y{\sqrt {m}}=0\,;$

d’où l’on tirera la valeur de $g+2h{\sqrt {m}},$ en fonction de $x+y{\sqrt {m}}\,;$ or, comme le radical ${\sqrt {m}}$ peut être pris également en moins, on aura de

même l’équation

$dc-db{\sqrt {m}}+\left(x-y{\sqrt {m}}\right)\left(dg-2dh{\sqrt {m}}\right)=0,$

laquelle, en faisant

$c-b{\sqrt {m}}=\psi \left(g-2h{\sqrt {m}}\right),$

donnera

$\psi '\left(g-2h{\sqrt {m}}\right)+x-y{\sqrt {m}}=0\,;$

moyennant quoi on connaîtra séparément $g$ et $h$ en $x$ et $y.$

Soit, pour plus de simplicité,

$x+y{\sqrt {m}}=t,\quad x-y{\sqrt {m}}=u,$

et soient $\mathrm {T,U}$ deux fonctions de $t$ et de $u\,;$ il est clair qu’on aura

$g+2h{\sqrt {m}}=\mathrm {T} ,\quad g-2h{\sqrt {m}}=\mathrm {U} \,;$

donc

$dc+db{\sqrt {m}}+td\mathrm {T} =0,\quad dc-db{\sqrt {m}}+ud\mathrm {U} =0,$

par conséquent

$c+b{\sqrt {m}}=-\int td\mathrm {T} ,\quad c-b{\sqrt {m}}=-\int ud\mathrm {U} .$

On aura donc ainsi

${\begin{alignedat}{2}c=&-{\frac {\int td\mathrm {T} +\int ud\mathrm {U} }{2}},\qquad &b=&-{\frac {\int td\mathrm {T} +\int ud\mathrm {U} }{2{\sqrt {m}}}},\\g=&{\frac {\mathrm {T+U} }{2}},&h=&{\frac {\mathrm {T-U} }{4{\sqrt {m}}}},\end{alignedat}}$

et, substituant ces valeurs dans la première des trois équations de condition, on aura, après les réductions,

$da-{\frac {t^{2}d\mathrm {T} -u^{2}d\mathrm {U} }{4{\sqrt {m}}}}=0,$

d’où l’on tire

$a={\frac {\int t^{2}d\mathrm {T} -\int u^{2}d\mathrm {U} }{4{\sqrt {m}}}}.$

Connaissant ainsi les valeurs des cinq quantités $a,b,c,g,h,$ il n’y aura plus qu’à les substituer dans l’équation $\mathrm {V} =0,$ et l’on aura

$z={\frac {\int t^{2}d\mathrm {T} -2t\int td\mathrm {T} +t^{2}\mathrm {T} }{4{\sqrt {m}}}}-{\frac {\int u^{2}d\mathrm {U} -2u\int ud\mathrm {U} +u^{2}\mathrm {U} }{4{\sqrt {m}}}},$

ce qui se réduit à cette forme plus simple

$z={\frac {\int dt\int \mathrm {T} dt-\int du\int \mathrm {U} du}{2{\sqrt {m}}}}\,;$

or, $\mathrm {T}$ étant une fonction quelconque de $t$ ou $x+y{\sqrt {m}},$ et $\mathrm {U}$ une fonction quelconque de $u$ ou $x-y{\sqrt {m}},$ il est visible que $\int dt\int \mathrm {T} dt$ et $\int du\int \mathrm {U} du$ seront aussi des fonctions quelconques des mêmes quantités. De sorte qu’on aura, en général,

$z=\Phi \left(x+y{\sqrt {m}}\right)+\Psi \left(x-y{\sqrt {m}}\right),$

ce qui s’accorde avec ce que l’on sait depuis longtemps.

Au reste on voit par cet Exemple, qui est d’ailleurs un des plus simples, que la méthode dont il s’agit, quoique directe et générale, est en quelque façon plus curieuse qu’utile, à cause des difficultés qui peuvent se rencontrer dans l’intégration des équations de condition ; c’est pourquoi nous ne nous arrêterons pas davantage sur ce sujet.

66. Il est beaucoup plus aisé de tirer l’intégrale générale aux premières différences de l’intégrale complète du premier ordre. Car nous avons vu (60) que, si $\mathrm {Z} =0$ est l’intégrale complète du premier ordre de l’équation du second ordre $\mathrm {Z} '=0,$ il suffit qu’il y ait dans l’équation $\mathrm {Z} =0$ deux constantes arbitraires $a$ et $b\,;$ et il est facile de prouver que cette équation satisfera également à l’équation $\mathrm {Z} '=0$ en y regardant $a$ et $b$ comme des variables, pourvu que l’on ait

${\frac {d\mathrm {Z} }{da}}da+{\frac {d\mathrm {Z} }{db}}db=0\,;$

de sorte qu’en faisant, en général,

$b=\varphi (a),$

on n’aura que cette seule équation de condition

${\frac {d\mathrm {Z} }{da}}+{\frac {d\mathrm {Z} }{db}}\varphi '(a)=0,$

laquelle servira à éliminer $a$ dans l’équation $\mathrm {Z} =0$ ; et la résultante contenant une fonction arbitraire sera nécessairement l’intégrale générale aux différences premières de la proposée $\mathrm {Z} '=0\,;$ il n’y a, comme l’on voit, aucune difficulté dans l’application de cette méthode.

L’équation

${\frac {d^{2}z}{dy^{2}}}-\mathrm {A} {\frac {d^{2}z}{dxdy}}+\mathrm {B} {\frac {d^{2}z}{dx^{2}}}=0$

a pour intégrale complète du premier ordre

${\frac {dz}{dy}}-m{\frac {dz}{dx}}=a+b(x+ny),$

comme on l’a vu plus haut (60). Pour en déduire l’intégrale générale on fera donc varier $a$ et $b,$ ce qui donnera

$da+(x+ny)db=0\,;$

et faisant $b=\varphi (a)$ on aura

$1+(x+ny)\varphi '(a)=0,$

d’où l’on voit que

$\varphi '(a)=-{\frac {1}{x+ny}},$

donc $a$ est égale à une fonction de $x+ny,$ et par conséquent $b$ est aussi

égal à une fonction de $x+ny\,;$ donc aussi $a+(x+ny)b$ est égal à une fonction de $x+ny,$ que je dénoterai par $\Phi (x+ny)\,;$ et il est visible que cette fonction pourra être quelconque. Ainsi l’intégrale générale du premier ordre sera

${\frac {dz}{dy}}-m{\frac {dz}{dx}}=\Phi (x+ny).$

De plus on a vu que la même équation a aussi pour intégrale complet du premier ordre

${\frac {dz}{dy}}-n{\frac {dz}{dx}}=h+g(x+my)\,;$

d’où l’on tirera, par un procédé semblable au précédent, la nouvelle intégrale générale

${\frac {dz}{dy}}-n{\frac {dz}{dx}}=\Psi (x+ny),$

la caractéristique $\Psi$ dénotant aussi une fonction quelconque.

Maintenant de même que dans l’endroit cité on a tiré l’intégrale finie complète des deux intégrales complètes du premier ordre, de même aussi pourra-t-on déduire des deux intégrales générales du premier ordre qu’on vient de trouver l’intégrale générale finie de l’équation proposée du second ordre.

En effet on a d’abord, en éliminant ${\frac {dz}{dy}},$

$(n-m){\frac {dz}{dx}}=\Phi (x+ny)-\Psi (x+my)\,;$

multipliant par $dx$ et intégrant relativement à $x,$ on aura

$(n-m)z=^{\text{‵}}\!\Phi (x+ny)-^{\text{‵}}\!\Psi (x+my)+\mathrm {Y} .$

$\mathrm {Y}$ dénote une fonction quelconque de $y$, et les caractéristiques $^{\text{‵}}\Phi ,^{\text{‵}}\!\Psi$ dénotent les intégrales des fonctions exprimées par les $\Phi ,\Psi .$ En éliminant de même la quantité ${\frac {dz}{dx}},$ on a

$(n-m){\frac {dz}{dy}}=n\Phi (x+ny)-m\Psi (x+my),$

et intégrant relativement à $y$ seul, on aura

$(n-m)z=^{\text{‵}}\!\Phi (x+ny)-^{\text{‵}}\!\Psi (x+my)+\mathrm {X} ,$

$\mathrm {X}$ étant une fonction quelconque de $x\,;$ or, comme ces deux valeurs de $z$ doivent être identiques, il faudra que

$\mathrm {X=Y} ,$

ce qui ne peut avoir lieu, à moins que les deux quantités $\mathrm {X}$ et $\mathrm {Y}$ ne soient constantes ; mais, comme les caractéristiques $^{\text{‵}}\Phi ,^{\text{‵}}\!\Psi$ expriment des fonctions quelconques, il serait superflu d’ajouter à ces fonctions une constante quelconque. Ainsi l’intégrale générale finie de la proposée sera

$(n-m)z=^{\text{‵}}\!\Phi (x+ny)-^{\text{‵}}\!\Psi (x+my).$

Cette méthode est générale, et l’on pourra toujours par son moyen trouver l’intégrale générale finie de toute équation à différences partielles du second ordre, dont on connaîtra deux intégrales complètes du premier ordre.

67. On pourrait croire qu’il suffit de connaître l’intégrale complète finie d’une équation à différences partielles du second ordre pour pouvoir trouver deux intégrales complètes du premier ordre de la même équation, comme cela a lieu pour les équations différentielles à deux variables (32) ; mais il n’en est pas ainsi des équations à différences partielles. En effet, dans ces sortes d’équations, lorsqu’il n’y a que trois variables, toute différentiation simple fournit deux équations, mais une différentiation double en fournit cinq, et une différentiation triple en fournira neuf, et ainsi de suite. De là vient que toute intégrale complète simple ou du premier ordre doit contenir deux constantes arbitraires, toute intégrale double ou du second ordre doit contenir cinq arbitraires, et ainsi du reste. Ainsi, pour pouvoir déduire une intégrale complète du premier ordre d’une intégrale complète du second, il faudrait pouvoir éliminer, par une simple différentiation de celle-ci, trois constantes vphitraires à la fois ; ou bien, si l’on n’en élimine que deux, il faudrait que la résultante fût telle, que les trois arbitraires restantes pussent être éliminées à la fois par une nouvelle différentiation simple ; or c’est ce qui est impossible, en général, et ne peut guère avoir lieu que dans des cas particuliers.

À plus forte raison sera-t-il impossible de déduire, en général, d’tune intégrale complète du troisième ordre une intégrale complète du second ou du premier, et ainsi du reste.

Nous allons rendre tout cela sensible par quelques Exemples.

68. Soit prise d’abord l’équation du n^o 58 ci-dessus, dont l’intégrale complète finie est

$z=k+ax+by+cx^{2}+gxy+hy^{2},$

$k$ étant $=f(a,b,c,g,h).$

Une différentiation simple donne les deux équations

${\frac {dz}{dx}}=a+2cx+gy,\quad {\frac {dz}{dy}}=b+gx+2hy\,;$

il faut donc voir si en combinant ces deux équations avec la précédente on peut chasser à la fois trois des cinq arbitraires $a,b,c,g,h.$ Je fais d’abord cette combinaison

$x{\frac {dz}{dx}}+y{\frac {dz}{dy}}=ax+by+2\left(cx^{2}+gxy+hy^{2}\right),$

et je remarque que si je retranche cette équation de celle ci-dessus multipliée par $2,$ j’en aurai une du premier ordre qui ne contiendra plus que les deux arbitraires $a$ et $b,$ pourvu qu’on suppose que $k$ soit une fonction de $a$ et $b$ seulement.

J’aurai donc dans ce cas l’équation

$2z-x{\frac {dz}{dx}}-y{\frac {dz}{dy}}=2f(a,b)+ax+by,$

qui sera l’intégrale complète du premier ordre de l’équation du second

ordre

${\begin{aligned}z=&x{\frac {dz}{dx}}+y{\frac {dz}{dy}}-{\frac {x^{2}}{2}}{\frac {d^{2}z}{dx^{2}}}-xy{\frac {d^{2}z}{dxdy}}-{\frac {y^{2}}{2}}{\frac {d^{2}z}{dy^{2}}}\\&+f\left({\frac {dz}{dx}}-x{\frac {d^{2}z}{dx^{2}}}-y{\frac {d^{2}z}{dxdy}},\ {\frac {dz}{dy}}-y{\frac {d^{2}z}{dy^{2}}}-x{\frac {d^{2}z}{dxdy}}\right).\end{aligned}}$

De là on pourra donc trouver, si l’on veut, l’intégrale générale du premier ordre de cette même équation (66).

Car en faisant varier $a$ et $b,$ et supposant

$b=\varphi (a),\quad db=\varphi '(a)da,$

on aura l’équation

$x+2{\frac {df(a,b)}{da}}+\left[y+2{\frac {df(a,b)}{db}}\right]\varphi '(a)=0,$

par laquelle on éliminera $a$ dans l’intégrale complète après avoir mis partout $\varphi (a)$ à la place de $b\,;$ la résultante sera l’intégrale générale cherchée.

Dans cet Exemple nous avons pu éliminer à la fois les trois arbitraires $c,g,h,$ ce qui nous a donné une intégrale complète du premier ordre ; mais pour avoir une autre intégrale complète il faudrait pouvoir éliminer à la fois trois autres des cinq arbitraires $a,b,c,g,h,$ ce qui n’est guère possible, comme on peut s’en assurer aisément par le calcul.

Au reste on voit aussi par l’Exemple précédent que, si la quantité était une fonction qui contînt les arbitraires $c,g,h,$ il ne serait plus possible de parvenir à une intégrale complète du premier ordre, par l’élimination simultanée de trois arbitraires.

69. Soit maintenant l’équation

${\frac {d^{2}z}{dx^{2}}}+{\frac {2y}{x}}{\frac {d^{2}z}{dxdy}}+{\frac {y^{2}}{x^{2}}}{\frac {d^{2}z}{dy^{2}}}=0,$

dont l’intégrale complète est

$z=a+bx+cy+g{\frac {y}{x}}+h{\frac {y^{2}}{x}},$

comme il est facile de s’en assurer par la substitution de cette valeur de $z.$

J’aurai par la différentiation

${\frac {dz}{dx}}=b-g{\frac {y}{x^{2}}}-h{\frac {y^{2}}{x^{2}}},\quad {\frac {dz}{dy}}=c+{\frac {g}{x}}+{\frac {2hy}{x}},$

d’où l’on tire cette combinaison

$x{\frac {dz}{dx}}+y{\frac {dz}{dy}}=bx+cy+h{\frac {y^{2}}{x}},$

et, retranchant cette équation de la précédente, j’aurai celle-ci

$z-x{\frac {dz}{dx}}-y{\frac {dz}{dy}}=a+g{\frac {y}{x}},$

laquelle ne contenant plus que deux arbitraires sera par conséquent l’intégrale complète du premier ordre de la proposée.

On pourra donc tirer de celle-ci l’intégrale générale du premier ordre, laquelle sera

$z-x{\frac {dz}{dx}}-y{\frac {dz}{dy}}=\varphi \left({\frac {y}{x}}\right).$

Je remarque de plus que l’on a

${\frac {dz}{dx}}+{\frac {y}{x}}{\frac {dz}{dy}}=b+c{\frac {y}{x}}+h{\frac {y^{2}}{x^{2}}},$

et que cette équation est telle, que les trois arbitraires qu’elle renferme peuvent être éliminées à la fois au moyen de ses deux différentielles

${\begin{aligned}&{\frac {d^{2}z}{dx^{2}}}+{\frac {y}{x}}{\frac {d^{2}z}{dxdy}}-{\frac {y}{x^{2}}}{\frac {dz}{dy}}=-{\frac {cy}{x^{2}}}-{\frac {2hy^{2}}{x^{3}}},\\&{\frac {d^{2}z}{dxdy}}+{\frac {y}{x}}{\frac {d^{2}z}{dy^{2}}}+{\frac {1}{x}}{\frac {dz}{dy}}\ \ \,={\frac {c}{x}}+{\frac {2hy}{x^{2}}}\,;\end{aligned}}$

car, en ajoutant ces deux équations-ci ensemble après avoir multiplié la

seconde par ${\frac {y}{z}},$ on a

${\frac {d^{2}z}{dx^{2}}}+{\frac {2y}{x}}{\frac {d^{2}z}{dxdy}}+{\frac {y^{2}}{x^{2}}}{\frac {d^{2}z}{dy^{2}}}=0,$

qui est l’équation même proposée.

De là il s’ensuit donc que l’équation dont il s’agit sera aussi une intégrale complète du premier ordre, et même on y pourra pour plus de simplicité supposer égale à zéro une quelconque des trois arbitraires qu’elle renferme.

On aura donc de cette manière cette nouvelle intégrale complète

${\frac {dz}{dx}}+{\frac {y}{x}}{\frac {dz}{dy}}=b+c{\frac {y}{x}},$

d’où l’on tirera aussi une nouvelle intégrale générale du premier ordre, laquelle sera

${\frac {dz}{dx}}+{\frac {y}{x}}{\frac {dz}{dy}}=\psi \left({\frac {y}{x}}\right).$

Au moyen de cette intégrale et de la précédente on pourra avoir sur-le-champ l’intégrale générale finie ; car multipliant la dernière par $x$ et l’ajoutant à la première on aura

$z=\varphi \left({\frac {y}{x}}\right)+x\psi \left({\frac {y}{x}}\right),$

$\varphi \left({\frac {y}{x}}\right)$ et $\psi \left({\frac {y}{x}}\right)$ étant deux fonctions arbitraires de ${\frac {y}{x}}.$

70. Nous nous dispenserons d’examiner les équations à différences partielles des ordres plus élevés, comme aussi celles où il y aurait plus de trois variables ; l’application des principes que nous venons d’établir à ces sortes d’équations ne doit pas être difficile à présent, et d’ailleurs ce Mémoire n’est déjà que trop long.

1. Ce Mémoire a été lu dans les Assemblées du 12 octobre et du 9 novembre 1775.
2. Œuvres de Lagrange, t. II, p. 17.
3. Œuvres de Lagrange, t. III, p. 549.