Annales de mathématiques pures et appliquées/Tome 03/Analise transcendante, article 4
# Annales de mathématiques pures et appliquées/Tome 03/Analise transcendante, article 4

## ANALISE TRANSCENDANTE.

À M. le Rédacteur des Annales,

On sait qu’en procédant à l’intégration des équations linéaires, à coefficiens constans, la substitution de eᵐˣ au lieu de y, semble en défaut, lorsque deux ou un plus grand nombre de racines de l’équation en m sont égales entre elles ; et que d’Alembert publia en 1748, une méthode très-ingénieuse pour écarter cette difficulté. Cependant, quel que soit le mérite de cette méthode, adoptée par Euler, dans son calcul intégral, et depuis par les auteurs de tous les traités sur cette matière, il m’a semblé qu’elle laissait à désirer un procédé plus rigoureux.

Je n’ignore pas qu’au fond le moyen employé par d’Alembert, et par les autres géomètres après lui, peut être entièrement justifié, soit par des considérations tirées de la théorie des limites, soit en faisant adroitement disparaître dans les termes à conserver (par un calcul un peu long quand il y a plus de deux racines égales) la quantité infiniment petite dont on a supposé que les racines venaient à différer. Mais cette petite différence k que, dans tous les traités que je connais, l’on annulle, sans que les quantités ck, ck², ck³, …

deviennent nulles en même temps, occasione toujours de l’embarras aux commençans, qui ne peuvent pas encore saisir le véritable esprit de la démonstration.

Je pense donc que la méthode suivante, qui n’est point sujette aux mêmes difficultés, et qui a l’avantage de donner immédiatement l’expression générale de l’intégrale, quels que soient l’ordre de l’équation et le nombre des racines égales, pourrait être introduite, avec avantage, dans les élémens ; et c’est pour lui donner la publicité nécessaire que je me suis déterminé, Monsieur, à vous l’adresser.

Soit l’équation

dⁿy+Adⁿ⁻¹y.dx+Bdⁿ⁻²y.dx²+…+Nydxⁿ=0;(P)

on sait que son intégrale complette est

y=a′eᵐ′ˣ+a″eᵐ″ˣ+a‴eᵐ‴ˣ+…+a⁽ⁿ⁾eᵐ⁽ⁿ⁾ˣ;

m′, m″, m‴, …m⁽ⁿ⁾ étant les n racines de l’équation

mⁿ+Amⁿ⁻¹+Bmⁿ⁻²+…+N=0, (Q)

qui provient de la substitution de eᵐˣ, au lieu de y, dans l’équation (P).

Quand toutes les racines de l’équation (Q) sont égales entre elles et à m, l’intégrale assignée se réduit à

y=(a′+a″+a‴+…+a⁽ⁿ⁾)eᵐˣ=α.eᵐˣ;

mais cette valeur de y n’est plus qu’une intégrale particulière, puisqu’elle ne renferme qu’une seule constante arbitraire α.

La simple substitution de eᵐˣ, au lieu de y, parait donc être ici en défaut, et ne pouvoir faire connaître la véritable intégrale de la proposée (P).

Cependant, puisque cette substitution satisfait toujours à l’équation différentielle, et puisque le défaut apparent de la méthode dépend d’une certaine relation existante entre les coefficiens constans A, B, …N; supposons

y=u.eᵐˣ;

u étant une fonction de x qu’on peut espérer de déterminer en telle sorte que l’intégrale renferme le nombre de constantes arbitraires nécessaire à la question.

Remarquons auparavant que, dans le cas où l’équation (Q) a toutes ses racines égales, comme elle est équivalente à (m−m)ⁿ=0, on a

A=−n/1.m, B=+n/1.n−1/2.m², C=−n/1.n−1/2.n−2/3.m³, ………………, N=±mⁿ;

le signe supérieur ou le signe inférieur devant être pris, dans la valeur de N, suivant que n est pair ou impair.

En conséquence, la proposée devient

dⁿy−n/1 dⁿ⁻¹y.mdx+n/1.n−1/2 dⁿ⁻²y.m²dx²−…±mⁿdxⁿ=0;

or, on sait que, lorsqu’une fonction y est de la forme u.t, u et t étant des fonctions de x, on a

dⁿy=dⁿ(u.t)=(du+dt)ⁿ,

pourvu qu’on écrive le développement avec les précautions convenables ; c’est-à-dire, qu’on a

dⁿy=tdⁿu+n/1 dⁿ⁻¹udt+n/1.n−1/2 dⁿ⁻²u.d²t+…+u.dⁿt;

mais, quand t=eᵐˣ, on a, par la nature de la fonction eᵐˣ,

dⁿt=t.(mdx)ⁿ;

donc,

dⁿy=t(dⁿu+n/1 dⁿ⁻¹u.mdx+n/1.n−1/2 dⁿ⁻²u.m²dx²+…+u.mⁿdxⁿ); (R)

mais la forme de la proposée, dans le cas actuel, est

dⁿy=n/1 dⁿ⁻¹y.mdx−n/1.n−1/2 dⁿ⁻²y.m²dx²+…∓mⁿdxⁿ;

en mettant donc, dans le second membre de cette dernière équation, au lieu de dⁿ⁻¹y, dⁿ⁻²y, .…, leurs valeurs en t, u, x, tirées de la forme générale (R), égalant ensuite entre elles les deux valeurs de dⁿy, et divisant de part et d’autre par t, on aura l’équation identique

dⁿu+n/1 dⁿ⁻¹u.mdx+n/1.n−1/2 dⁿ⁻²u.m²dx²+…−n/1 dⁿ⁻¹u.mdx−2 n/1.n−1/2 dⁿ⁻²u.m²dx²−…+n/1.n−1/2 dⁿ⁻²u.m²dx²+…−…|=0.

ce qui donne dⁿu=0; et, en intégrant,

u=a′+a″x+a‴x²+…+a⁽ⁿ⁾.xⁿ⁻¹;

donc

y=u.eᵐˣ=(a′+a″x+a‴x²+…+a⁽ⁿ⁾.xⁿ⁻¹).eᵐˣ;

valeur qui, renfermant n constantes arbitraires, est l’intégrale complette de la proposée (P).

Lorsqu’il n’y a que p racines de l’équation (Q) qui soient égales entre elles et à m ; n étant supposé égal à p+q l’intégrale se réduit à

y=(a′+a″+…+a⁽ᵖ⁾).eᵐˣ+a⁽ᵖ⁺¹⁾.eᵐ⁽ᵖ⁺¹⁾ˣ+…+a⁽ᵖ⁺𐞥⁾.eᵐ⁽ᵖ⁺𐞥⁾ˣ;

ou

y=α.eᵐˣ+a⁽ᵖ⁺¹⁾.eᵐ⁽ᵖ⁺¹⁾ˣ+…+a⁽ᵖ⁺𐞥⁾.eᵐ⁽ᵖ⁺𐞥⁾ˣ;

intégrale qui n’est que particulière, puisqu’au lieu de n ou p+q constantes arbitraires, elle n’en renferme que 1+q.

Dans ce cas, l’équation (Q) revient à

(m−mᵖ).(m−m⁽ᵖ⁺¹⁾)…(m−m⁽ᵖ⁺𐞥⁾)=0.

Considérons séparément le premier facteur, et posons l’équation

mᵖ+A′mᵖ⁻¹+B′mᵖ⁻²+…+H′=0.(Q′)

Il est évident, par le cas général que nous venons de traiter, que cette équation se rapporte à l’équation différentielle

dᵖy+A′dᵖ⁻¹y.dx+B′dᵖ⁻²y.dx²+…+H′y.dxᵖ=0.(P′)

dont toutes les solutions y=eᵐˣ seraient égales entre elles ; en sorte que son intégrale se présenterait sous la forme particulière y=αeᵐˣ.

Si donc, en raisonnant comme dans le cas général, nous instituons les mêmes calculs, nous trouverons, pour l’intégrale complette de cette équation (P′),

y=u.eᵐˣ=(a′+a″x+a‴x²+…+a⁽ᵖ⁾.xᵖ⁻¹).eᵐˣ;

valeur qui renfermera p constantes arbitraires.

Mais, d’après la propriété des équations différentielles linéaires, l’on sait que, si l’on a n valeurs particulières de y, leur somme donne immédiatement l’expression générale de cette fonction.

Donc, en réunissant la valeur précédente de y aux solutions fournies par les q facteurs inégaux de l’équation (Q), lesquelles renferment chacune une constante arbitraire, nous aurons enfin pour intégrale complette de la proposée (P)

y=(a′+a″x+a‴x²+…+a⁽ᵖ⁾.xᵖ⁻¹).eᵐˣ+a⁽ᵖ⁺¹⁾.eᵐ⁽ᵖ⁺¹⁾ˣ+…+a⁽ᵖ⁺𐞥⁾.eᵐ⁽ᵖ⁺𐞥⁾ˣ.

