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

## ANALISE TRANSCENDANTE.

### Application du calcul aux différences partielles à la résolution de quelques problèmes d’analise ;

Le calcul aux différences partielles, qui doit son origine à des questions de géométrie et de mécanique, a été postérieurement appliqué, d’une manière très-heureuse, par des géomètres du premier ordre, à des questions de pure analise. Ce sont quelques essais de ce genre d’application que nous nous proposons de présenter ici, en employant successivement cette branche de calcul et au développement des fonctions polynomiales en séries, et au problème du retour des suites.

M. Paoli s’est déjà proposé de déduire des seuls principes du calcul différentiel tout ce qui est nécessaire pour parvenir au développement en séries des fonctions polynomiales : c’est du même sujet que nous nous proposons de nous occuper ici. Notre méthode étant un peu plus simple que celle de l’illustre Italien, nos résultats doivent aussi être moins compliqués que les siens.

Soit u une fonction de y qu’il soit question de développer suivant les puissances de x ;y étant donnée par l’équation

y=a+a₁x+a₂x²+a₃x³+…

On aura d’abord cette suite d’équations

d⁡u/d⁡x=d⁡u/d⁡y.d⁡y/d⁡x=d⁡u/d⁡y (a₁+2a₂x+3a₃x²+4a₄x³+…)

d⁡u/d⁡a=d⁡u/d⁡y.d⁡y/d⁡a =d⁡u/d⁡y, d⁡u/d⁡a₁=d⁡u/d⁡y.d⁡y/d⁡a₁ =x d⁡u/d⁡y, d⁡u/d⁡a₂=d⁡u/d⁡y.d⁡y/d⁡a₂ =x² d⁡u/d⁡y, …………… d⁡u/d⁡aₙ=d⁡u/d⁡y.d⁡y/d⁡aₙ =xⁿ d⁡u/d⁡y ;}(1)

dans lesquelles a, a₁, a₂, …aₙ sont traitées comme indépendantes.

Éliminant d⁡u/d⁡y entre chacune des équations (1) et la dernière, on trouve

d⁡u/d⁡aₙ=xⁿ d⁡u/d⁡a, d⁡u/d⁡aₙ=xⁿ⁻¹ d⁡u/d⁡a₁, …………… d⁡u/d⁡aₙ=xⁿ⁻ᵏ d⁡u/d⁡aₖ.}(2)

Mettant ensuite dans d⁡u/d⁡x les valeurs de d⁡u/d⁡y, x d⁡u/d⁡y, x² d⁡u/d⁡y, … données par les équations (1), on trouve

d⁡u/d⁡x=a₁ d⁡u/d⁡a+2a₂ d⁡u/d⁡a₁+3a₃ d⁡u/d⁡a₂+4a₄ d⁡u/d⁡a₃+…

Posons, pour abréger,

u=A, d⁡u/d⁡x=A₁, 1/2! d²⁡u/d⁡x²=A₂, 1/3! d³⁡u/d⁡x³=A₃, … 1/n! dⁿ⁡u/d⁡xⁿ=Aₙ.…

La première des équations (2) donne

d.dⁿ⁻ᵅ⁡u/d⁡xⁿ⁻ᵅ /d⁡aₙ=dⁿ⁻ᵅ.(xⁿ d⁡u/d⁡a )/d⁡xⁿ⁻ᵅ=0;

pourvu qu’après la différentiation on suppose x=0. Partant

d⁡A_{n−α}/d⁡aₙ=0.

On tirerait de même de la dernière équation (2)

(d⁡Aₙ₊ₚ)/d⁡aₙ=(d⁡Aₖ₊ₚ)/d⁡aₖ ;

et enfin de l’équation (3)

nAₙ=a₁ d⁡Aₙ₋₁/d⁡a+2a₂ d⁡Aₙ₋₁/d⁡a₁+3a₃ d⁡Aₙ₋₁/d⁡a₂+…

équation qu’en vertu de l’équation de condition (5) qu’on peut changer en celle-ci,

nAₙ=a₁ d⁡Aₙ₋₁/d⁡a+2a₂ d⁡Aₙ₋₂/d⁡a+3a₃ d⁡Aₙ₋₃/d⁡a+…;

ou encore

nAₙ=a₁ (d⁡A_{n+α})/(d⁡a_{α+1})+2a₂ (d⁡A_{n+α})/(d⁡a_{α+2})+3a₃ (d⁡A_{n+α})/(d⁡a_{α+3})+…

L’équation (6) exprime suivant quelle loi chacune des quantités A₁, A₂, A₃, … dérive de celle qui la précède immédiatement.

La formule (7), ne renfermant de différentiations que par rapport à a, permet d’employer les valeurs particulières ou numériques de a₁, a₂, a₃, … et de faire ainsi les réductions à mesure qu’elles se présenteront ; ce qui, dans bien de cas, la rendra préférable.

Enfin, la formule (8) donne le moyen de revenir de l’une quelconque des quantités A₂, A₃, A₄, … à celle qui la précède. Elle devient illusoire lorsque n=0 ; ce qu’il n’était pas difficile de prévoir.

Au moyen de la relation que l’équation (8) établit entre les divers coefficiens de x, dans le développement de u, on peut donner une infinité de formes différentes à la formule (6). Nous avons rapporté seulement les plus remarquables ; mais son emploi peut devenir plus intéressant. En effet ; les équations (4) et (8) peuvent se concentrer en celle-ci

d⁡Aₘ/d⁡aₙ=d⁡A_{m−α}/d⁡a_{n−α},

de laquelle on tirera facilement

d²⁡Aₘ/d⁡aₙ²=d²⁡A_{m−α}/d⁡a_{n−α}.d⁡aₙ=d²⁡A_{m−2α}/(d⁡a_{n−α})² ;

et, en général,

dᵖ⁡Aₘ/d⁡aₙᵖ=dᵖ⁡A_{m−pα}/(d⁡a_{n−α})ᵖ.

Maintenant, on a

Aₘ=Aₘ+aₙ d⁡Aₘ/d⁡aₙ+aₙ²/2! d²⁡Aₘ/d⁡aₙ²+aₙ³/3! d³⁡Aₘ/d⁡aₙ³+…

pourvu qu’après les différenciations on fasse aₙ=0, dans les quantités Aₘ, d⁡Aₘ/d⁡aₙ, d²⁡Aₘ/d⁡aₙ², … du second membre. L’on aura donc, en vertu de l’équation (9)

Aₘ=Aₘ+aₙ d⁡A_{m−α}/d⁡a_{n−α}+aₙ²/2!.d²⁡A_{m−2α}/(d⁡a_{n−α})²+aₙ³/3!.d³⁡A_{m−3α}/(d⁡a_{n−α})³+…;

moyennant les mêmes restrictions ; et comme, dans la dernière équation, les différentiations ne sont plus relatives à aₙ, on trouvera la valeur de Aₘ, dans la supposition de aₙ quelconque, si l’on a les valeurs de Aₘ, A_{m−α}, A_{m−2α}, … dans la supposition de aₙ=0, sans que pour cela on soit obligé de recommencer les calculs.

Par le moyen des formules précédentes, on pourra s’élever, de proche en proche, à la valeur de Aₙ, en fonction de a₁, a₂, a₃, …A, et des coefficiens différentiels de cette dernière quantité ; mais cette marche, d’ailleurs très-laborieuse, ne serait fondée que sur l’analogie. Voici, pour le même objet, une méthode en même temps plus expéditive et plus rigoureuse.

Si l’on fait, pour abréger,

t=a₁+a₂x+a₃x²+a₄x³+…;

on aura

y=a+tx;

u=f⁡(a+tx);

ou bien, en développant,

u=A+d⁡A/d⁡a.tx+d²⁡A/d⁡a².t²x²/2!+d³⁡A/d⁡a³.t³x³/3!+…

On a d’ailleurs

tⁿ⁻¹/(n−1)!=1/n! d⁡tⁿ/d⁡a₁, tⁿ⁻²/(n−2)!=1/n! d².tⁿ/d⁡a₁², …………, tⁿ⁻ᵏ/(n−k)!=1/n! dᵏ.tⁿ/d⁡a₁ᵏ ;

ainsi

u=A+1/n! (x d⁡A/d⁡a.dⁿ⁻¹.tⁿ/d⁡a₁ⁿ⁻¹+x² d²⁡A/d⁡a².dⁿ⁻².tⁿ/d⁡a₁ⁿ⁻²+…+xⁿ dⁿ⁡A/d⁡aⁿ.tⁿ)

si donc l’on suppose

tⁿ=B+B₁x+B₂x²+B₃x³+…

on aura

Aₙ=1/n! { d⁡A/d⁡a.dⁿ⁻¹.Bₙ₋₁/d⁡a₁ⁿ⁻¹+d²⁡A/d⁡a².dⁿ⁻².Bₙ₋₂/d⁡a₁ⁿ⁻²+…+dⁿ⁡A/d⁡aⁿ.B};

ou, en renversant l’ordre des termes,

Aₙ=1/n! { dⁿ⁡A/d⁡aⁿ.B+dⁿ⁻¹⁡A/d⁡aⁿ⁻¹.d⁡B₁/d⁡a₁+… d⁡A/d⁡a.dⁿ⁻¹.Bₙ₋₁/d⁡a₁ⁿ⁻¹ }.

Mais, en suivant la marche qui nous a conduit aux formules (6) et (7), nous trouverions

kBₖ=a₂ d⁡Bₖ₋₁/d⁡a₁+2a₃ d⁡Bₖ₋₁/d⁡a₂+3a₄ d⁡Bₖ₋₁/d⁡a₃+…

Partant, si l’on fait

B=C, d⁡B₁/d⁡a₁=C₁, d²⁡B₂/d⁡a₁²=C₂, … dᵏ⁡Bₖ/d⁡a₁ᵏ=Cₖ,

on aura

kCₖ=a₂ d²⁡Cₖ₋₁/d⁡a₁.d⁡a₁+2a₃ d²⁡Cₖ₋₁/d⁡a₁.d⁡a₂+3a₄ d²⁡Cₖ₋₁/d⁡a₁.d⁡a₃+…;

et

kCₖ=a₂ d²⁡Cₖ₋₁/d⁡a₁²+2a₃ d³⁡Cₖ₋₂/d⁡a₁³+3a₄ d⁴⁡Cₖ₋₃/d⁡a₁⁴+…;

et par conséquent

Aₙ=1/n! (C dⁿ⁡A/d⁡aⁿ+C₁ dⁿ⁻¹⁡A/d⁡aⁿ⁻¹+C₂ dⁿ⁻²⁡A/d⁡aⁿ⁻²+…+Cₙ₋₁ d⁡A/d⁡a )

La réunion de la formule (12) avec une des formules (10) et (11) remplit l’objet proposé.

Voici enfin une méthode très-simple, pour construire tout d’un coup l’entier développement de u. Soit fait

1/1−htx=1+htx+h²t²x²+h³t³x³+…=H;

on aura alors

u=AH+1/1!1!.d⁡A/d⁡a.d⁡H/d⁡h+1/2!2!.d²⁡A/d⁡a².d²⁡H/d⁡h²+1/3!3!.d³⁡A/d⁡a³ d³⁡H/d⁡h³+…

pourvu que, dans le second membre, on fasse h=0, après les differentiations. D’où l’on voit qu’étant donné le développement de H, on aura très-facilement celui de u. Supposons donc

H=1+E₁x+E₂x²+E₃x³+…;

on trouvera facilement l’équation

Eₙ=h(a₁Eₙ₋₁+a₂Eₙ₋₂+a₃Eₙ₋₃+…aₙ₋₁E₁+aₙ),

pour déterminer chaque terme de H, en fonction de ceux qui le précèdent.

Dans ce qui précède, nous avons supposé que u n’était fonction que d’une seule quantité y, qui n’était elle-même fonction que de x ; mais à présent supposons

u=f⁡(y, y′),

y=a_{0, 0}+a_{1, 0}x+a_{2, 0}x²+…a_{0, 1}z+a_{1, 1}xz+…a_{0, 2}z²+…y′=b_{0, 0}+b_{1, 0}x+b_{2, 0}x²+…b_{0, 1}z+b_{1, 1}xz+…b_{0, 2}z²+…

y et y′ étant la somme de tous les termes de la forme a_{m, n}xᵐzⁿ, b_{m, n}xᵐzⁿ que l’on peut faire, en prenant pour m, n tous les nombres entiers positifs, zéro compris, pourront se mettre sous cette forme

y=Σ(a_{m, n}xᵐzⁿ), y′=Σ(b_{m, n}xᵐzⁿ).

Cela posé, regardant comme entièrement indépendantes les diverses quantités x, z, a_{m, n}, b_{m, n}, et différentiant dans cette vue, on trouvera facilement

d⁡u/d⁡x=Σ{mxᵐ⁻¹zⁿ(a_{m, n} d⁡u/d⁡y+b_{m, n} d⁡u/d⁡y′ )},

d⁡u/d⁡z=Σ{nxᵐzⁿ⁻¹(a_{m, n} d⁡u/d⁡y+b_{m, n} d⁡u/d⁡y′ )},

d⁡u/d⁡a_{p, q}=xᵖz𐞥 d⁡u/d⁡y,

d⁡u/d⁡b_{p, q}=xᵖz𐞥 d⁡u/d⁡y′.

En vertu de ces deux dernières équations, on peut changer (1) et (2) en celles-ci

d⁡u/d⁡x=Σ{m(a_{m, n} d⁡u/d⁡a_{m−1, n}+b_{m, n} d⁡u/d⁡b_{m−1, n} )},

d⁡u/d⁡z=Σ{n(a_{m, n} d⁡u/d⁡a_{m, n−1}+b_{m, n} d⁡u/d⁡b_{m, n−1} )}.

Mettant, dans l’équation (3), p+h et q+k, au lieu de p et q, elle deviendra

d⁡u/(d⁡a_{p+h, q+k})=xᵖ⁺ʰz𐞥⁺ᵏ d⁡u/d⁡y,

qui, comparé avec l’équation (3), donne

d⁡u/(d⁡a_{p+h, q+k})=xʰzᵏ d⁡u/d⁡a_{p, q} ;

et on tirera de même de l’équation (4)

d⁡u/(d⁡b_{p+h, q+k})=xʰzᵏ d⁡u/d⁡a_{p, q}.

Maintenant, représentons en général par A_{r, s} la valeur de

1/r!s! (dʳ⁺ˢ⁡u)/d⁡xʳd⁡zˢ,

lorsqu’après les différentiations on y suppose x et z nuls. On aura d’abord

u=A_{0, 0}+A_{1, 0}x+A_{2, 0}x²+…+A_{0, 1}z+A_{1, 1}zx+…+A_{0, 2}z²+…

et ensuite, on tirera des équations (5), (6), (7), (8) les formule suivantes

rA_{r, s}=Σ{m(a_{m, n} d⁡A_{r−1, s}/d⁡a_{m−1, n}+b_{m, n} d⁡A_{r−1, s}/d⁡b_{m−1, n} )},

sA_{r, s}=Σ{n(a_{m, n} d⁡A_{r, s−1}/d⁡a_{m, n−1}+b_{m, n} d⁡A_{r, s−1}/d⁡b_{m, n−1} )},

(d⁡A_{f+h, g+k})/(d⁡a_{p+h, q+k})=d⁡A_{f, g}/d⁡a_{p, q},

(d⁡A_{f+h, g+k})/(d⁡b_{p+h, q+k})=d⁡A_{f, g}/d⁡b_{p, q} ;

au moyen des relations données par les équations (11), (12) on pourra donner une infinité de formes différentes aux formules (9), (10) : nous nous contenterons de rapporter les principales, qui sont

rA_{r, s}=Σ{m(a_{m, n} d⁡A_{r−m, s−n}/d⁡a_{0, 0}+b_{m, n} d⁡A_{r−m, s−n}/d⁡b_{0, 0} )},

sA_{r, s}=Σ{n(a_{m, n} d⁡A_{r−m, s−n}/d⁡a_{0, 0}+b_{m, n} d⁡A_{r−m, s−n}/d⁡b_{0, 0} )},

rA_{r, s}=Σ{m(a_{m, n} (d⁡A_{r+h, s+k})/(d⁡a_{m+h, n+k})+b_{m, n} (d⁡A_{r+h, s+k})/(d⁡b_{m+h, n+k}) )}, (15)sA_{r, s}=Σ{n(a_{m, n} (d⁡A_{r+h, s+k})/(d⁡a_{m+h, n+k})+b_{m, n} (d⁡A_{r+h, s+k})/(d⁡b_{m+h, n+k}) )}, (16)rBᵣ=Σ′{m(aₘ′ d⁡Bᵣ₋₁/d⁡a_{m−1, 0}+bₘ′ d⁡Bᵣ₋₁/d⁡b_{m−1, 0} )}, (17)rBᵣ=Σ′{m(aₘ′ d⁡Bᵣ₋ₘ/d⁡a_{0, 0}+bₘ′ d⁡Bᵣ₋ₘ/d⁡b_{0, 0} )}, (18)rBᵣ=Σ′{m(aₘ′ (d⁡Bᵣ₊ₕ)/(d⁡a_{m+h, 0})+bₘ′ (d⁡Bᵣ₊ₕ)/(d⁡b_{m+h, 0}) )}.(19)

Dans ces trois dernières formules, le signe Σ′ ne se rapporte qu’à m ; et nous avons fait, pour abréger,

Bᵣ=A_{r, 0}+z/x A_{r−1, 1}+z²/x² A_{r−2, 2}+… zʳ/xʳ A_{0, r}, aₘ′=a_{m, 0}+z/x a_{m−1, 1}+z²/x² a_{m−2, 2}+… zᵐ/xᵐ a_{0, m}, bₘ′=b_{m, 0}+z/x b_{m−1, 1}+z²/x² b_{m−2, 2}+… zᵐ/xᵐ b_{0, m}.

Quant à la manière d’y parvenir, voici pour cela une méthode que je crois plus simple que l’emploi des formules (9), (10), (11), (12).

Les équations (1), (2) donnent

d⁡u/d⁡x+z/x d⁡u/d⁡z=Σ{(m+n)(xᵐ⁻¹zⁿ)(a_{m, n} d⁡u/d⁡y+b_{m, n} d⁡u/d⁡y′ )}

ou, en faisant m+n=α,

d⁡u/d⁡x+z/x d⁡u/d⁡z=Σ{α zⁿ/xⁿ xᵅ⁻¹(a_{α−n, n} d⁡u/d⁡y+b_{α−n, n} d⁡u/d⁡y′ )},

ou encore

d⁡u/d⁡x+z/x d⁡u/d⁡z=Σ′{α(a_{α}′ d⁡u/d⁡a_{α−1, 0}+b_{α}′ d⁡u/d⁡b_{α−1, 0} )};

en mettant respectivement

d⁡u/d⁡a_{m−1, 0}

d⁡u/d⁡b_{m−1, 0},

à la place de xᵐ⁻¹ d⁡u/d⁡y et xᵐ⁻¹ d⁡u/d⁡y′, ce qui est permis, en vertu des équations (3) et (4) ; et le signe Σ′ ne se rapportant qu’à α. Enfin, on tirera de l’équation (20) une formule qui ne différera de la formule (17) qu’en ce que α y sera au lieu de m. Cette formule une fois trouvée, on en déduira facilement (18) et (19), au moyen des relations (11) et (12).

Enfin, l’on sait que le développement de u est égal à une suite de termes de la forme

(dᵐ⁺ⁿ⁡A_{0, 0})/d⁡a_{0, 0}ᵐd⁡b_{0, 0}ⁿ.(y−a_{0, 0})ᵐ/m!.(y′−b_{0, 0})ⁿ/n! ;

faisant ensuite

(y−h)ᵖ(y′−k)ᵖ=H,

on aura

(y−a_{0, 0})ᵖ⁻ʳ(y′−b_{0, 0})ᵖ⁻ˢ/(p−r)!(p−s)!=1/p!p! (dʳ⁺ˢ⁡H)/d⁡a_{0, 0}ʳd⁡b_{0, 0}ˢ ;

pourvu qu’après les différentiations l’on fasse h=a_{0, 0}, k=b_{0, 0} dans le développement de H. On déduira donc par là, du développement de H, tous les termes de celui de u, dans lesquels la somme des exposans de x et z ne sera pas plus grande que p.

Soit encore

H=1/{1−h(y−a_{0, 0})}{1−k(y′−b_{0, 0})},

on aura

(y−a_{0, 0})ᵐ(y′−b_{0, 0})ⁿ=1/m!n! (dᵐ⁺ⁿ.H)/d⁡hᵐ.d⁡kⁿ ;

pourvu que, dans le développement de H on suppose h et k nuls après les différentiations ; ce qui fournit une autre méthode pour trouver le développement de u.

Nous pourrions donner, pour parvenir a ce même développement, une infinité d’autres méthodes plus ou moins compliquées ; mais nous nous sommes contentés de rapporter les plus simples. Nous pourrions encore nous occuper du cas où u est fonction de plus de deux fonctions y, y′…; du cas où y, y′…, seraient elles-mêmes fonctions de plus de deux variables indépendantes x, z; mais ces divers cas ne présenteront aucune difficulté sérieuse à ceux qui auront bien saisi l’esprit de notre méthode. Au surplus, de quelque nombre d’application que nos formules puissent être susceptibles, nous n’avons pas pensé qu’il dût être nécessaire d’en faire comprendre l’usage par des exemples qui n’auraient fait que donner à ce mémoire un surcroit d’étendue que nous avons sur-tout cherché à éviter. Nous terminerons sur ce sujet en observant que, bien que nous ayons supposé que u était fonction de y, y′, sans a_{m, n}, b_{m, n}, on peut cependant étendre notre méthode à ce cas, et cela, par un artifice ingénieux du à l’illustre auteur de la Mécanique céleste ; il consiste à remplacer momentanément ces quantités par d’autres, que l’on regardera comme constantes dans les différentiations. On peut faire une semblable remarque pour le cas où u n’est fonction que de y seulement.

Dans son Traité de calcul différentiel et intégral, (2.ᵉ édit., tom. I, pag. 298), M. Lacroix observe que, quelque élégant que soit l’emploi du Théorème de Lagrange, dans l’opération du retour des suites, il ne saurait être pourtant regardé comme indiquant la loi des formules auxquelles il conduit.

Frappé de cette remarque, qui nous a paru très-fondée, nous avons cherché une solution du problème qui ne fût pas sujette à cet inconvénient ; et la suivante nous paraît remplir le but. À la vérité, elle sera jugée peut-être moins générale et moins élégante que celle de Lagrange ; mais aussi est-il bien loin de notre pensée de prétendre lutter contre cet illustre géomètre.

Soit l’équation

a=a₁z+a₂z²+a₃z³+a₄z⁴+…

et proposons-nous d’en tirer la valeur de z, ordonnée suivant les puissances de a. Cette valeur, quelle qu’elle soit, est évidemment une fonction des quantités a, a₁, a₂, a₃, …, que nous regarderons comme indépendantes. Ainsi, différentiant successivement l’équation (1), par rapport à chacune de ces quantités, et posant, pour abréger,

k=a₁+2a₂z+3a₃z²+4a₄z³+…

nous aurons

k d⁡z/d⁡a−1=0, k d⁡z/d⁡a₁+z=0, k d⁡z/d⁡a₂+z²=0, …………k d⁡z/d⁡aₙ+zⁿ=0.}

Éliminant k entre la première des équations (2) et chacune des suivantes, on trouve

d⁡z/d⁡a₁+z d⁡z/d⁡a =0, d⁡z/d⁡a₂+z² d⁡z/d⁡a =0, ………… d⁡z/d⁡aₙ+zⁿ d⁡z/d⁡a =0.}

D’ailleurs, en remettant pour k sa valeur dans la première des équations (2), elle deviendra

1=a₁ d⁡z/d⁡a+2a₂z d⁡z/d⁡a+3a₃z² d⁡z/d⁡a+4a₄z³ d⁡z/d⁡a+…

ou bien, en vertu des équations (3),

a₁ d⁡z/d⁡a=1+2a₂ d⁡z/d⁡a₁+3a₃ d⁡z/d⁡a₂+4a₄ d⁡z/d⁡a₃+…

Supposons

z=A₁a+A₂a²+A₃a³+A₄a⁴+…;

en substituant cette valeur dans l’équation (4), la comparaison des termes semblables dans les deux membres donnera

a₁=1, 2a₁=2a₂ d⁡A₁/d⁡a₁+3a₃ d⁡A₁/d⁡a₂+4a₄ d⁡A₁/d⁡a₃+…, 3a₁=2a₂ d⁡A₂/d⁡a₁+3a₃ d⁡A₂/d⁡a₂+4a₄ d⁡A₂/d⁡a₃+…, …………………na₁=2a₂ d⁡Aₙ₋₁/d⁡a₁+3a₃ d⁡Aₙ₋₁/d⁡a₂+4a₄ d⁡Aⁿ⁻¹/d⁡a₃+…;

formules qui indiquent suivant quelle loi le coefficient d’une puissance quelconque de a dérive de celui de la puissance immédiatement inférieure.

L’équation (3) donne

zⁿ=−n d⁡∫zd⁡a/d⁡aₙ₋₁ ;

ou ; en mettant pour z sa valeur,

zⁿ=−n (d⁡( A₁/2 a²+A₂/3 a³+A₃/4 a⁴+…))/d⁡aₙ₋₁ ;

d’ailleurs z″ est de la forme

Baⁿ+B₁aⁿ⁺¹+B₂aⁿ⁺²+…;

partant

d⁡Aₙ₋₁₋ₖ/d⁡aₙ₋₁=0,

tant que k ne sera pas nul ; et

zⁿ=−( d⁡Aₙ₋₁/d⁡aₙ₋₁ aⁿ+d⁡Aₙ/d⁡aₙ₋₁.(naⁿ⁺¹)/(n+1)+(d⁡Aₙ₊₁)/d⁡aₙ₋₁.(naⁿ⁺²)/(n+2)+…);

qui donnera facilement la valeur de zⁿ, quand celle de z sera connue.

Souvent on a

a=b₁x+b₂x²+b₃x³+b₄x⁴+…;

et alors c’est suivant les puissances de x qu’il faut ordonner le développement de z. Dans ce cas on a

d⁡z/d⁡x=d⁡z/d⁡a d⁡a/d⁡x ;

ou bien

d⁡z/d⁡x=d⁡z/d⁡a (b₁+2b₂x+3b₃x²+4b₄x³+…).

On trouverait de même

d⁡z/d⁡b₁=x d⁡z/d⁡a, d⁡z/d⁡b₂=x² d⁡z/d⁡a, d⁡z/d⁡b₃=x³ d⁡z/d⁡a, ……………, d⁡z/d⁡bₙ=x d⁡z/d⁡a, }

partant

d⁡z/d⁡x=b₁ d⁡z/d⁡a+2b₂ d⁡z/d⁡b₁+3b₃ d⁡z/d⁡b₂+4b₄ d⁡z/d⁡b₃+…

mettant, dans cette dernière équation, au lieu de d⁡z/d⁡a, sa valeur tirée de l’équation (4), et multipliant par a₁, on la changera en celle-ci

a₁ d⁡z/d⁡x=b₁+2a₁b₂ d⁡z/d⁡b₁+3a₁b₃ d⁡z/d⁡b₂+…+2b₁a₂ d⁡z/d⁡a₁+3b₁a₃ d⁡z/d⁡a₂+…

Supposant ensuite

z=C₁x+C₂x²+C₃x³+C₄x⁴+…

et mettant cette valeur dans la formule précédente ; on trouvera d’abord

a₁C₁=b₁,

et ensuite, en général,

na₁Cₙ=2a₁b₂ d⁡Cₙ₋₁/d⁡b₁+3a₁b₃ d⁡Cₙ₋₁/d⁡b₂+4a₁b₄ d⁡Cₙ₋₁/d⁡b₃+…+2b₁a₂ d⁡Cₙ₋₁/d⁡a₁+3b₁a₃ d⁡Cₙ₋₁/d⁡a₂+4b₁a₄ d⁡Cₙ₋₁/d⁡a₃+…}

pour la loi qui lie les coefficiens de deux puissances quelconques consécutives de x, dans le développement de z.

Si dans l’équation

zⁿ=−n d⁡∫zd⁡a/d⁡aₙ₋₁,

on met pour d⁡a sa valeur

d⁡x(b₁+2b₂x+3b₃x²+4b₄x³+…),

et pour z son développement ; on trouvera que le coefficient de zᵐ dans zⁿ est

−n/m . (d(b₁Cₘ₋₁+2b₂Cₘ₋₂+3b₃Cₘ₋₃+⋯))/da.

D’où l’on voit qu’étant donné le développement de z, on en déduira facilement celui de zⁿ.

