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

## ANALISE TRANSCENDANTE.

### Du calcul des dérivations, ramené à ses véritables principes, ou théorie du développement des fonctions, et du retour des suites ;

Depuis l’invention du Théorème de Taylor, sur le développement des fonctions d’un binôme, et du Théorème de Lagrange, sur le retour des fonctions et des séries, bien des géomètres se sont occupés d’étendre et de généraliser les découvertes de ces deux géomètres célèbres ; et, sous le rapport de la théorie générale, on peut dire que les résultats auxquels ils sont parvenus ne laissent plus rien à désirer ; mais les formules qui les contiennent, quelques précieuses qu’elles soient comme solutions générales, ne font qu’indiquer une suite d’opérations ultérieures, souvent si compliquées qu’elles découragent le calculateur le plus intrépide. Il restait donc à trouver une méthode simple, facile et uniforme, pour exécuter complètement et immédiatement tous ces développemens, tant directs que de retour. Les géomètres allemands sont les premiers qui ont réussi dans cette recherche : leurs travaux ont donné naissance à un nouveau calcul, appelé Analise combinatoire par son inventeur Hindenhurg. Ce calcul résout à la vérité la question, mais d’une manière trop disparate avec les procédés ordinaires de l’analise : il oblige à former d’abord séparément les groupes de lettres, et ensuite leurs coefficiens numériques, pour lesquels on a besoin de tables de combinaisons calculées d’avance. Il était réservé à Arbogast de donner la solution générale, complète et analitique de cette question difficile, dans son Calcul des dérivations. Malheureusement cet ouvrage est entaché de plusieurs défauts très-graves, qui ont dégoûté les géomètres de sa lecture, et ont empêché qu’il ne fût étudié et connu autant qu’il le mérite. Ces défauts sont 1.ᵒ de n’avoir pas assez justifié l’introduction de ses nouvelles notations ; 2.ᵒ de n’avoir pas défini assez nettement ses dérivées et ses dérivations ; 3.ᵒ de déduire sa théorie d’un principe qui n’est ni assez clair ni assez évident (n.ᵒ 6) ; 4.ᵒ de l’exposer d’une manière trop longue et trop embarrassée ; 5.ᵒ enfin d’avoir noyé des résultats vraiment remarquables dans une foule de choses qui sont, pour ainsi dire, hors d’œuvre, et sans liaison avec l’objet principal de son ouvrage ; de sorte que ce qui pouvait être présenté dans quelques feuilles d’impression est devenu un gros in-4.ᵒ .

Je me propose, dans ce petit écrit, de remédier, le mieux que je pourrai, à ces défauts de l’ouvrage d’Arbogast, en déduisant la véritable théorie du calcul des dérivations du seul théorème de Taylor, sans l’emploi d’aucun principe nouveau ; de sorte que ce calcul ne sera, à proprement parler, qu’une extension de ce théorème.

Afin de rendre l’exposition de cette théorie plus rapide, et de présenter de suite aux géomètres toute la partie usuelle de ce calcul, je me contenterai quelquefois de généraliser les résultats par des conclusions d’induction ; sauf à démontrer ces conclusions dans un article à part. Pour la même raison, je réléguerai dans des remarques toutes les observations ; soit sur les notations, soit sur le fond même de la théorie.

### Article premier.

1. On sait, par le théorème de Taylor, que toute fonction d’un binôme α+x (x étant une variable qui reste indéterminée), peut être développée en une série de cette forme (1)

f⁡(α+x)=a+a₁x+a₂x²+a₃x³+…;

les coeffiiens a, a₁, a₂, a₃, … étant déterminés par les valeurs que prennent la fonction f⁡(α+r) et ses coefficiens différentiels (d.f⁡(α+x))/d⁡x, (d².f⁡(α+x))/1.2d⁡x², (d³.f⁡(α+x))/1.2.3d⁡x³, …,

lorsqu’après la différentiation on y fait x=0.

Nous représenterons ces valeurs particulières des coefficiens différentiels par D.f⁡α, 1/1.2 D².f⁡α, 1/1.2.3 D³.f⁡α…; et nous appellerons, avec Arbogast, les quantités D.f⁡α, D².f⁡α, D³.f⁡α, …, dérivées, première, seconde, troisième,… de fα. Ainsi la dérivée Dⁿ.f⁡α, d’un ordre quelconque n, n’est autre chose que ce que devient le coefficient différentiel (dⁿ.f⁡(α+x))/d⁡xⁿ, lorsqu’on y fait x=0. Nous aurons donc

(2){a₁=1/1 Df⁡α=1/1 D⁡a, a₂=1/1.2 D²f⁡α=1/1.2 D²⁡a, a₃=1/1.2.3 D³f⁡α=1/1.2.3 D³⁡a, …………………aₙ=1/1.2.…n Dⁿf⁡α=1/1.2.…n Dⁿ⁡a;}d′ou^{`} l′on tirere^{´}ciproquement{Df⁡a=a₁, D²⁡a=1.2a₂, D³⁡a=1.2.3a₃, ……………Dⁿ⁡a=1.2.…naₙ.}(3)

Les valeurs (2) étant substituées dans l’équation (1) donnent

(4)

f(α+x)=f⁡α+Df⁡α.x+1/1.2 D²f⁡α.x²+1/1.2.3 D³f⁡α.x³+…

=a+D⁡a.x+1/1.2 D²⁡a.x²+1/1.2.3 D³⁡a.x³+…

Mais ici il faut observer que ce développement peut devenir impossible ; ce qui arrive toutes les fois que f⁡α et ses déîivées deviennent infinies : ce cas a lieu pour la fonction Log.x, par exemple.

2. Remarque. La loi des Dérivations, ou de la formation des dérivées successives, est évidemment la même que celle des differentiations, à la seule exception près que nous supprimons les dénominateurs inutiles D⁡α, D²⁡α, D³⁡α, …; ce qui revient à supposer D⁡α=1. En effet, on sait, et il est d’ailleurs évident que le coefficient différentiel d’un ordre quelconque de la fonction d’un binôme α+x est le même, soit qu’on la différencie en regardant x comme variable et α comme constant, soit qu’on la différencie en regardant α comme variable, et x comme constant ; on a donc (dⁿ.f⁡(α+x))/d⁡xⁿ=(dⁿ.f⁡(α+x))/d⁡αⁿ ; or, en faisant x=0, dans (dⁿ.f⁡(α+x))/d⁡αⁿ, on obtient dⁿ.f⁡α/d⁡αⁿ=Dⁿf⁡α. Il est donc démontré que la loi des dérivations est la même que celle des différentiations ; et il s’ensuit aussi que les dérivées d’un ordre quelconque ne sont autre chose que les coefficiens différentiels du même ordre, pris relativement à la quantité constante α, que l’on feint être variable.

Il se présente ici naturellement une objection que l’on a faite, dès l’origine, contre la notation du calcul des dérivations : c’est que les opérations dérivatives étant les mêmes que celles du calcul différentiel, il fallait les indiquer par les mêmes notations. Je réponds d’abord qu’absolument parlant la chose eût été possible ; mais que, pour l’ordre et la précision, et d’après les règles d’une saine logique, des opérations qui, bien qu’identiques pour la forme, diffèrent entièrement pour le fond, doivent être représentées par des signes différens. Je dis que ces opérations diffèrent entièrement pour le fond ; car les signes différentiels se rapportent aux variables et à leurs variabilité de grandeur, tandis, que les signes dérivatifs ne se rapportent qu’aux seules constantes, et que les dérivées successives n’indiquent qu’une dépendance d’ordre et de succession dans les termes d’un développement. Si Lagrange a été autorisé à introduire une notation nouvelle, dans le calcul des variations, pour indiquer une opération entièrement identique avec la différentiation, et qui se rapporte aux variables même, seulement parce qu’elle ne se rapporte pas à leur variabilité de grandeur, mais à leur variabilité de forme ; à plus forte raison sera-t-il permis, ou plutôt nécessaire, de représenter par une notation particulière une opération qui ne se rapporte pas même aux variables, ni à aucune espèce de variabilité.

Une autre raison, qui suffirait à elle seule pour justifier l’introduction d’une notation particulière pour les dérivations, c’est qu’elles peuvent se trouver, et se trouvent réellement souvent combinées avec les différentielles, dans une même formule ; il faut donc qu’on ne puisse pas les confondre : ce qui arriverait infailliblement, si elles étaient représentées par une même notation.

Quant à la suppression des dénominateurs D⁡α, D⁡α², D⁡α³, …, leur inutilité seule suffit pour la justifier. Si des personnes habituées aux considérations d’infiniment petits tiennent à conserver ces dénominateurs, dans le calcul différentiel, où leur considération abrège quelquefois les raisonnemens dans des questions de géométrie et de mécanique, mais où elle peut aussi égarer ; il n’y a pas la moindre raison de les conserver dans l’analise pure, ni, à plus forte raison, dans les dérivations où toute idée d’infiniment petit serait plus que déplacée.

3. Réciproquement, tout polynôme de la forme a+a₁x+a₂x² +a₃x³+… (terminé ou non), peut être représenté par a+D⁡a.x +1/1.2 D²⁡a.x²+1/1.2.3 D³⁡a.x³+… où, entre les coefficiens a, a₁, a₂, a₃, … et les dérivées successives de a, on a les relations (3).

En effet, en supposant, ce qui est toujours permis,

a+a₁x+a₂x²+a₃x³+…=f⁡(α+x),

on a, d’après la définition des dérivées (n.ᵒ 1),

a=f⁡α, a₁=Df⁡α=D⁡a, a₂=1/1.2 D²f⁡α=1/1.2 D²⁡a,

a₃=1/1.2.3 D³f⁡α=1/1.2.3 D³⁡a, …;

d’où l’on tire les relations (3).

Si le polynôme est terminé, et n’a qu’un nombre n de termes, aₙ₋₁ sera le dernier des coefficlens ; et tous les suivans aₙ, aₙ₊₁, aₙ₊₂, …

ainsi que leurs valeurs correspondantes, en dérivées de a, seront égaux à zéro.

4. Remarque. Il est bon d’observer que l’équation (5) peut être satisfaite d’une infinité de manières, et que α est entièrement indéterminé. En supposant α=0, on a

a=f⁡0, a₁=Df⁡0=D⁡a, a₂=1/1.2 D²f⁡0=1/1.2 D²⁡a,

a₃=1/1.2.3 D³f⁡0=1/1.2.3 D³⁡a, …;

ce qui donne encore les relations (3).

5. Proposons-nous maintenant de développer la fonction d’un polynôme quelconque, ordonné selon les puissances ascendantes de la variable, en une série qui procède selon les mêmes puissances ; c’est-à-dire, de déterminer les coefficiens du second membre de l’équation suivante :

ϕ(a+a₁x+a₂x²+a₃x³+…)²=A+A₁x+A₂x²+A₃x³+…

D’après le n.ᵒ 3, le polynôme sous le signe de la fonction peut être représenté par f⁡(α+x) ; le premier membre de cette équation, peut donc être mis sous la forme ϕf⁡(α+x). Il suffit donc de substituer, dans l’équation (4), ϕf au lieu de f ; ce qui donne

ϕ(a+a₁x+a₂x²+a₃x³+…)=ϕf⁡α+D.ϕf⁡α.x+1/2 D².ϕf⁡α.x³+1/6 D³.ϕf⁡α.x³+…=ϕa+D.ϕa.x+1/2 D².ϕa.x²+1/6 D².ϕa.x²+…

en observant cependant que, dans ce cas, chacun des signes de dérivation indique des opérations ultérieures à exécuter, que nous expliquerons dans la remarque suivante : c’est pourquoi nous avons mis un point après chaque signe de dérivation, pour marquer cette différence.

En comparant terme à terme l’équation (7) avec celle (6), on en tire les valeurs suivantes :

A=ϕa, A₁=D.ϕa, A₂=1/2 D².ϕa, A₃=1/6 D³.ϕa, …

et par conséquent

ϕa=A, D.ϕa=A₁, D².ϕa=2A₂, D³.ϕa=6A₃, …,

6. Remarque. Nous avons démontré, dans la remarque du n.ᵒ 2, qu’on a, en général, Dⁿf⁡α=dⁿf⁡α/d⁡αⁿ ; on a donc aussi

D.ϕf⁡α=d⁡(ϕf⁡α/d⁡α, D².ϕf⁡α=d²⁡(ϕf⁡α/d⁡α², D³.ϕf⁡α=d³⁡(ϕf⁡α/d⁡α³, …

D.ϕa=d⁡(ϕa)/d⁡α), D².ϕa=d²⁡(ϕa)/d⁡α², D³.ϕa=d³⁡(ϕa)/d⁡α³, …

Or, d’après les règles ordinaires de la différentiation des fonctions de fonctions, ou des fonctions de variables qui sont elles-mêmes fonctions de la variable principale, on a

d⁡(ϕf⁡α)/d⁡α=d⁡(ϕf⁡α)/df⁡α.df⁡α/d⁡α,

d²⁡(ϕf⁡α)/d⁡α²=d⁡[ d⁡(ϕf⁡α)/df⁡α.df⁡α/d⁡α ]/d⁡α=d⁡(ϕf⁡α)/df⁡α.d².f⁡α/d⁡α²+d²⁡(ϕf⁡α)/(df⁡α)².(df⁡α/d⁡α)²,

ou, en mettant a au lieu de f⁡α,

d⁡(ϕa)/d⁡α=d⁡(ϕa)/d⁡a.d⁡a/d⁡α,

d²⁡(ϕa)/d⁡α²=d⁡[ d⁡(ϕa)/d⁡a.d⁡a/d⁡α ]/d⁡α=d⁡(ϕa)/d⁡α.d²⁡a/d⁡α²+d²⁡(ϕa)/d⁡a².(d⁡a/d⁡α)²,

Passant donc aux notations dérivatives, en supprimant les dénominateurs, nous aurons

D.ϕf⁡α=D⁡ϕf⁡α.D⁡ϕf⁡α,

D²⁡ϕf⁡α=D⁡ϕf⁡α.D²⁡ϕf⁡α+D²⁡ϕf⁡α.(Df⁡α)²,

et

D.ϕa=D⁡ϕa.D⁡a,

D²⁡ϕa=D⁡ϕa.D²⁡a+D²⁡ϕa.D⁡a²,

Le développement de ces deux premières dérivées suffit pour expliquer la différence qui existe entre les dérivées sans point et celles qui sont suivies d’un point ; entre D⁡ϕa, D²⁡ϕa, … et D.ϕa, D².ϕa, … Les premières sont les dérivées de ϕa, en supposant D⁡a=1, D²⁡a=0, D³⁡a=0, …, et les autres sont les dérivées de ϕa, en supposant que a est fonction de α, que ses dérivées successives sont D⁡a, D²⁡a, D³⁡a, …, et qu’elles ne sont pas toutes égales à zéro. En un mot, les dérivées sans point, D⁡ϕa, D²⁡ϕa, … sont les coefficiens du développement de ϕ(a+x) et les dérivées suivies d’un point D.ϕa, D².ϕa, … sont les coefficiens du développement de ϕ(a+D⁡a.x+1/2 D²⁡a.x²+1/6 D³⁡a.x³+…)

7. En exécutant, d’après la remarque précédente, les opérations indiquées par les dérivées suivies d’un point, on obtient pour les six premières, en suivant les règles ordinaires de la différenciation, lorsqu’aucune différentielle n’est constante :

D.ϕa=D⁡ϕa.D⁡a, D².ϕa=D⁡ϕa.D²⁡a+D²⁡ϕa.D⁡a², D².ϕa=D⁡ϕa.D²⁡a+D²⁡ϕa.D⁡a², D³.ϕa=D⁡ϕa.D³⁡a+D²⁡ϕa.3D⁡aD²⁡a+D³⁡ϕaD⁡a³, D⁴.ϕa=D⁡ϕa.D⁴⁡a+D²⁡ϕa[4D⁡aD³⁡a+3(D²⁡a)²]+D³⁡ϕa.6D⁡a²D²⁡a+D⁴⁡ϕa.D⁡a⁴, D⁵.ϕa=D⁡ϕa.D⁵⁡a+D²⁡ϕa(5D⁡aD⁴⁡a+10D²⁡aD³⁡a)+D³⁡ϕa[10D⁡a²D³⁡a+15D⁡a(D²⁡a)²], +D⁴⁡ϕa.10D⁡a³D²⁡a+D⁵⁡ϕa.D⁡a⁵,

D⁶.ϕa=D⁡ϕa.D⁶⁡a+D²⁡ϕa[6D⁡aD⁵⁡a+15D²⁡aD⁴⁡a+10(D³⁡a)²]+D³⁡ϕa[15D⁡a²D⁴⁡a+60D⁡aD²⁡aD³⁡a+15(D²⁡a)³]+D⁴⁡ϕa[20D⁡a³D³⁡a+45D⁡a²(D²⁡a)²+D⁵⁡ϕa.15D⁡a⁴D²⁡a+D⁶⁡ϕa.D⁡a⁶,

et ainsi de suite.

En substituant ces valeurs dans les équations (8), et remplaçant les dérivées de a par leurs valeurs (3), on obtient

(10){A=ϕa, A₁=D.ϕa=D⁡ϕa.a₁, A₂=1/2 D².ϕa=D⁡ϕa.a₂+1/2 D²⁡ϕa.a₁², A₃=1/6 D³.ϕa=D⁡ϕa.a₃+1/2 D²⁡ϕa.2a₁a₂+1/6 D³⁡ϕa.a₁³, A₄=1/24 D⁴.ϕa=D⁡ϕa.a₄+1/2 D²⁡ϕa(2a₁a₃+a₂²)+1/6 D³⁡ϕa.3a₁²a₂+1/24 D⁴⁡ϕa.a₁⁴, A₅=1/120 D⁵.ϕa=D⁡ϕa.a₅+ 1/2 D²⁡ϕa(2a₁a₄+2a₂a₃)+1/6 D³⁡ϕa(3a₁²a₃+3a₁a₂²)+1/24 D⁴⁡ϕa.4a₁³a₂+1/120 D⁵⁡ϕa.a₁⁵, A₆=1/720 D⁶.ϕa=D⁡ϕa.a₆+1/2 D²⁡ϕa(2a₁a₅+a₂a₄+a₃²)+1/6 D³⁡ϕa(3a₁²a₄+2.3a₁a₂a₃+a₂³)+1/24 D⁴⁡ϕa(4a₁³a₄+4.3/2 a₁²a₂²)+1/120 D⁵⁡ϕa.5a₁⁴a₂+1/720 D⁶⁡ϕa.a₁⁶,

et ainsi de suite.

Nous voici donc parvenus au développement complet des sept premiers termes de la série (6), qui est l’objet fondamental du calcul des dérivations, sans supposer autre chose que le théorème de Taylor, et les règles ordinaires de la différentiation. En examinant de près la composition successive de ces termes, on en conclut aisément la règle pratique suivante, pour déduire immédiatement un terme quelconque de celui qui le précède.

8. Pour déduire à développement de Aₙ₊₁ de celui de Aₙ, les lettres étant disposées d’après leur ordre de succession ;

1.ᵒ On ne fera varier, dans chaque terme, que la dernière lettre ou sa puissance, en suivant les règles ordinaires de la différentiation, et en mettant simplement a₁ pour D⁡a, a₂ pour D⁡a₁, a₃ pour D⁡a₂, …, sans autre coefficient que l’unité ;

2.ᵒ On fera varier de plus (d’après les mêmes règles de la différentiation) l’avant-dernière lettre, sa puissance ou sa fonction, si elle se trouve être la lettre qui, dans l’ordre numérique des indices, précède immédiatement la dernière du terme ; et comme la puissance de la dernière lettre augmente alors d’une unité, on divisera par son exposant ainsi augmenté.

Pour faire une application de cette règle, et pour mieux en faire comprendre l’usage, nous allons déduire le développement de A₆ de celui de A₅ [formules (10)].

Le premier terme D⁡ϕa.a₅ donne, d’après la première partie de la règle, qui seule y est applicable, D⁡ϕa.a₆. Le second terme 1/2 D²⁡ϕa(2a₁a₄+2a₂a₃) donne 1/2 D²⁡ϕa(2a₁a₅+a₂a₄+a₃²), dont les deux premiers termes sont dûs à la première partie de la règle, et le dernier à la seconde partie. Le troisième terme 1/6 D³⁡ϕa(3a₁²a₃+3a₁a₂²) donne 1/6 D³⁡ϕa(3a₁²a₄+2.3a₁a₃+a₂³), dont les deux premiers termes d’après la première partie de la règle et le dernier d’après la seconde. Le quatrième terme 1/24 D⁴⁡ϕa.4a₁³a₂ donne 1/24 D⁴⁡ϕa(4a₁³a₄+4.3/2 a₁²a₂²), dont le premier terme d’après la première partie de la règle et le suivant d’après la seconde. Enfin, le terme 1/120 D⁵⁡ϕa.a₁⁵ donne 1/120 D⁵⁡ϕa.5a₁⁴a₂, d’après la première parue de la règle et 1/720 D⁶⁡ϕa.a₁⁶, d’après la seconde. En rassemblant tous ces différens termes, on obtient exactement le développement de A₆, tel que nous l’avons donné [formules (10)].

On voit, d’après cet exemple, que la règle est d’une exécution très-facile, et qu’elle fournit immédiatement les termes successifs du développement, tout ordonnés et réduits à leur plus simple expression, sans qu’on ait, pour ainsi dire, d’autre peine que celle de les écrire. Il est vrai que, jusqu’à présent, cette règle n’est qu’une conclusion d’induction ; mais nous nous proposons de la démontrer dans un des articles suivant. Nous réservons pour le même article la règle pour former immédiatement un terme quelconque du développement, indépendamment de ceux qui le précèdent.

9. Remarque I. On a pu remarquer, en examinant la composition des formules (10), que le signe de fonction, ainsi que les signes de dérivation, n’affectent que la première lettre du polynôme dont il s’agit de développer la fonction. Arbogast appelle cette première lettre origine de dérivations ou premier terme de polynôme, et toutes les suivantes quantités polynômiales. Il résulte de cette observation que la composition des termes successifs du développement, en quantités polynômiales, reste la même, quelle que soit la fonction à développer ; et que toute la différence, dans le développement des diverses sortes de fonctions, consiste dans les valeurs des dérivées D⁡ϕa, D²⁡ϕa, D³⁡ϕa, …, qui n’affectent que la première lettre du polynôme. Ainsi on a, pour (a+a₁x+a₂x²+…)ᵐ

D⁡ϕa=maᵐ⁻¹, D²⁡ϕa=m(m−1)aᵐ⁻²,

D³⁡ϕa=m(m−1)(m−2)aᵐ⁻³, …;

pour e^{(a+a1x+a2x2+…)}

D⁡ϕa=eᵃ, D²⁡ϕa=eᵃ, D³⁡ϕa=eᵃ, …;

pour Cos.(a+a₁x+a₂x²+…)

D⁡ϕa=−Sin.a, D²⁡ϕa=−Cos.a, D³⁡ϕa=+Sin.a, …;

et ainsi de suite, pour d’autres formes de fonctions.

Si le polynôme sous le signe de fonction était terminé, et composé de n termes, on aurait aₙ=0, aₙ₊₁=0, aₙ₊₂=0, …;

il suffirait donc alors de rejeter, dans les formules (10), tous les termes où il entrerait, comme facteur, une des quantités polynômiales dont l’indice serait supérieure à n−1 ; et, dans l’application de la règle du n.ᵒ précédent, on s’arrêterait, dans chaque coefficient au terme affecté de aₙ₋₁.

10. Remarque II. L’inspection des termes successifs (10) du développement de l’équation (6) fait aisément découvrir la loi remarquable qui y règne. Elle consiste en ce que Aₙ est composé de n termes, formés des dérivées successives

D⁡ϕa/1, D²⁡ϕa/1.2, D³⁡ϕa/1.2.3, … Dⁿ⁡ϕa/1.2.3…n,

dont les coefficiens se composent de la manière suivante : 1.ᵒ le coefficient de Dʳ⁡ϕa/1.2.3…r

est composé de tous les produits de r lettres qu’on peut former avec les quantités polynômiales a₁, a₂, a₃, …, aₙ de manière que la somme des indices de chaque produit soit n, chaque lettre étant supposée écrite autant de fois qu’il y a d’unités dans son exposant ; 2.ᵒ les coefficiens numériques de chaque produit indiquent le nombre de permutations dont les lettres de ce produit sont susceptibles. Ainsi, le coefficient de D³⁡ϕa/1.2.3 dans A₆ est composé des trois produits 3a₁a₁a₄+2.3a₁a₂a₃+a₂a₂a₂ qui sont les seuls qu’on peut former avec trois lettres, de manière que la somme des indices dans chacun soit égale à 6 : leurs coefficiens numériques indiquent, comme on voit, le nombre des permutations dont les lettres qui les composent sont susceptibles.

Cette remarque, traduite en deux règles pratiques, l’une relative à la formation des groupes de lettres, et l’autre relative à celle des coefficiens numériques, d’après la théorie des combinaisons, constitue l’Analise combinatoire des géomètres allemands.

11. Si, au lieu de la fonction d’un polynôme ϕ(a+a₁x+a₂x²+…) on avait à développer la fonction ϕf⁡(α+x) d’une fonction de binôme, il suffirait, d’après le n.ᵒ 5, de substituer, dans les équations (10), ou dans les résultats obtenus par la règle du n.ᵒ 8, pour les quantités polynômiales a, a₁, a₂, a₃, …, leurs valeurs (2).

Enfin, si l’on avait à développer la fonction d’une fonction de polynôme ψϕ(a+a₁x+a₂x²+a₃x³+…), on aurait, d’après le même, n.ᵒ 5,

ψϕ(a+a₁x+a₂x²+a₃x³+…)

=ψ(ϕa+D.ϕa.x+1/2 D².ϕa.x²+1/6 D³.ϕa.x³+…);

a, a₁, a₂, a₃, …,

ϕa, D.ϕa,

1/2 D².ϕa, 1/6 D³.ϕa, …, qui doivent être elles-mêmes développées selon les règles du n.ᵒ 8.

Il serait même aisé, dans ce cas, d’obtenir immédiatement le développement de la fonction proposée ; car, en mettant dans l’équation (7) ψϕ à la place de ϕ, on obtiendrait

ψϕ(a+a₁x+a₂x²+a₃x³+…)

=ψϕa+D.ϕa.x+1/2 D².ϕa.x²+…;

et au moyen d’une légère extension donnée à la règle du n.ᵒ 8, on trouverait

D.ψϕa=D⁡ψϕa.D⁡ϕa.a₁, 1/2 D².ψϕa=D⁡ψϕa.(D⁡ϕa.a₂+1/2 D²⁡ϕa.a₁²)+ 1/2D²⁡ψϕa(D⁡ϕa)².a₁², 1/6 D³.ψϕa=D⁡ψϕa.(D⁡ϕa.a₃+1/2 D²⁡ϕa.2a₁a₂+1/6 D³⁡ϕa.a₁³)+1/2 D²⁡ψϕa[(D⁡ϕa)²2.a₁a₂+2D⁡ϕa.1/2 D²⁡ϕa.a₁³]+1/6 D³⁡ψϕa.(D⁡ϕa)³.a₁³,

et ainsi de suite.

Le développement précédent équivaut à celui d’une fonction triple d’un binôme ψϕf⁡(α+x). En effet, pour avoir le développement de cette fonction triple, il suffit de substituer, dans les formules précédentes, pour a, a₁, a₂, a₃, …, les valeurs (2), en fonction de α.

Nous ne pousserons pas plus loin ces observations, sur le développement des fonctions multiples ; ce que nous venons de dire suffit pour faire apercevoir la possibilité de cette extension, au moyen du calcul des dérivations.

12. Jusqu’à présent nous n’avons considéré que les fonctions d’un seul polynôme, et nous avons complètement résolu le problème de leur développement successif : il nous resterait maintenait à résoudre la même question pour les fonctions de plusieurs polynômes, ainsi que pour celles des polynômes à double ou à triple entrée ; mais les limites que nous avons dû prescrire à cet écrit, ne nous permettent pas d’entrer dans ces détails : nous nous contenterons d’ajouter encore le développement des produits de deux polynômes et de deux fonctions de polynômes, parce que nous en aurons besoin pour la théorie du retour des suites.

Proposons-nous d’abord de développer le produit

(a+a₁x+a₂x²+a₃x³+…)×(b+b₁x+b₂x²+b₃x³+…)

en une série de la forme

(A+A₁x+A₂x²+A₃x³+…)

En effectuant [a multiplication, par le procédé ordinaire, on obtient

{A=ab, A₁=ab₁+a₁b, A₂=ab₁+a₁b, A₃=ab₃+a₁b₂+a₂b₁+a₃b, ………………Aₙ=abₙ+a₁bₙ₋₁+a₂bₙ₋₂+…+aₙ₋₂b₂+aₙ₋₁b₁+aₙ;

où la loi est évidente.

D’après le n.ᵒ 3, le problème peut se mettre sous la forme

(13) (a+D⁡a.x+1/2 D²⁡a.x²+1/6 D³⁡a.x³+…)×(b+D⁡b.x+1/2 D²⁡b.x²+1/6 D³⁡b.x³+…)

=A+A₁x+A₂x²+A₃x³+…

En substituant donc, dans les équations (12), pour a₁, a₂, a₃, …, b₁, b₂, b₃, … leurs valeurs en dérivées de a et de b, d’après le n.ᵒ 3, on obtient

{A=ab, A₁=D⁡(ab)=a.D⁡b+D⁡a.b, A₂=1/2 D²⁡(ab)=a.1/2 D²⁡b+D⁡a.D⁡b+1/2 D²⁡a.b, A₃=1/6 D³⁡(ab)=a.1/6 D³⁡b+D⁡a.1/2 D²⁡b+1/2 D²⁡a.D⁡b+1/6 D³⁡a.b, ………………………Aₙ=1/1.2…n Dⁿ⁡(ab)=a.Dⁿ⁡b/1.2…n+D⁡a.Dⁿ⁻¹⁡b/1.2…(n−1)+…+Dⁿ⁻²⁡a/1.2…(n−1) D⁡b+Dⁿ⁡a/1.2…n b;

où a et b (d’après les n.ᵒˢ 3 et 4)sont supposés être respectivement fonctions de deux quantités arbitraires α et β, dont les dérivées sont D⁡α=1, , D²⁡α=0, , D⁡β=1, , D²⁡β=0, …

De la comparaison des formules (12) et (14) résulte la règle pratique suivante ;

Pour déduire Aₙ₊₁ de Aₙ [formules (12)], 1.ᵒ ne faites varier dans chaque terme, que b et ses dérivées, en écrivant b₁ pour D⁡b, b₂ pour D⁡b₁, b₃ pour D⁡b₂…; 2.ᵒ dans le dernier terme seulement, qui contient la plus haute dérivée de a, faites varier cette dérivée, en écrivant aₙ₊₁ pour D⁡aₙ.

13. Soit maintenant à développer le produit de deux fonctions de polynômes

ϕ(a+a₁x+a₂x²+a₃x³+…)×ψ(b+b₁x+b₂x²+b₃x³+…)

=A+A₁x+A₂x²+A₃x³+…

D’après le n.ᵒ 5, cette équation peut être mise sous la forme

(ϕa+D.ϕa.x+1/2 D².ϕa.x²+…)(ψb+D.ψb.x+1/2 D².ψb.x²+…)

=A+D.A.x+1/2 D².A.x²+1/6 D³.A.x³+…

qui, étant comparée à celle (13), fait voir qu’il suffit de remplacer, dans celle-ci, a par ϕa, b par ψb, et les signes de dérivation sans points par des signes de dérivation avec points ; on aura donc, d’après le même n.ᵒ 5, les équations analogues à celles (14), c’est-à-dire,

A=ϕa.ψb, A₁=D⁡(ϕa.ψb)=ϕa.D⁡ψb+D.ϕa.ψb, A₂=1/2 D²⁡(ϕa.ψb)=ϕa.1/2 D².ψb+D.ϕa.D.ψb+1/2 D².ϕa.ψb, A₃=1/6 D³⁡(ϕa.ψb)=ϕa.1/6 D³.ψb+D⁡ϕa.1/2 D².ψb+1/2 D².ϕa.D.ψb+1/6 D³.ϕa.ψb, ………………………Aₙ=1/1.2…n Dⁿ⁡(ϕa.ψb)=ϕa.1/1.2…n Dⁿ.ψb+D.ϕa.1/1.2…(n−1) Dⁿ⁻¹.ψb+…+1/1.2…(n−1) Dⁿ⁻¹.ϕa.D.ψb+1/1.2…n Dⁿ.ϕa.ψb;

où chaque dérivée de ϕa et de ψb doit être développée comme les équations (10) ; c’est-à-dire, d’après la règle du n.ᵒ 8. Les formules précédentes contiennent donc, au fond, tout ce qu’il faut pour la solution complète de la question ; mais, pour ne rien laisser à désirer, nous allons en déduire les moyens d’exécuter immédiatement le développement complet de ces formules.

14. En exécutant les dérivations indiquées, au moyen de la règle 8, effectuant les multiplications, et ordonnant le tout par rapport aux exposans des dérivées, on obtient

(18){A=ϕa.ψbA₁=ϕa.D⁡ψb.b₁+D⁡ϕa.ψb.a₁A₂=ϕa.D⁡ψb.b₂+ϕa.1/2 D²⁡ψb.b₁²+D⁡ϕa.ψb.a₂+D⁡ϕa.D⁡ψb.a₁b₁+1/2 D²⁡ϕa.ψb.a₁²A₃=ϕa.D⁡ψ.b₃+ϕa.1/2 D².ψb.2b₁b₂+ϕa.1/6 D³⁡ψb.b₁³+D⁡ϕa.ψb.a₃+D⁡ϕa.D⁡ψb(a₁b₂+a₂b₁)D⁡ϕa.1/2 D²⁡ψb.a₁b₂²+1/2 D²⁡ϕa.ψb.2a₁a₂+1/2 D²⁡ϕa.D⁡ψb.a₁²b₁+1/6 D³⁡ϕa.ψb.a₁³A₄=ϕa.D⁡ψ.b₄+ϕa.1/2 D².ψb(2b₁b₃+b₂²)+ϕa.1/6 D³⁡ψb.3b₁²b₂+D⁡ϕa.ψb.a₄+D⁡ϕa.D⁡ψb(a₁b₃+a₂b₂+a₃b₁)+D⁡ϕa.1/2 D²⁡ψb.(2a₁b₁b₂+a₂b₁²)+1/2 D²⁡ϕa.ψb(2a₁a₃+a₂²)+1/2 D²⁡ϕa.D⁡ψb.(a₁²b₂+2a₁a₂b₁)+1/6 D³⁡ϕa.ψb.3a₁²a₂+ϕa.1/24 D⁴⁡ψb.b₁⁴+D⁡ϕa.1/6 D³⁡ψb.a₁b₁³+1/2 D²⁡ϕa.1/2 D²⁡ψb.a₁²b₁²+1/6 D³⁡ϕa.D⁡ψb.a₁³b₁+1/24 D⁴⁡ϕa.ψb.a₁⁴A₅=ϕa.D⁡ψ.b₅+ϕa.1/2 D².ψb(2b₁b₄+2b₂b₃)+ϕa.1/6 D³⁡ψb(3b₁²b₃+3b₁b₂²)+D⁡ϕa.ψb.a₅+D⁡ϕa.D⁡ψb(a₁b₄+a₂b₃+a₃b₂+a₄b₁)+D⁡ϕa.1/2 D²⁡ψb(2a²b₁b₃+a₁b₁+2a₂b₁b₂+a₃b₁)+1/2 D²⁡ϕa.ψb(2a₁a₄+2a₂a₃)+1/2 D²⁡ϕa.D⁡ψb.(a₁²b₃+2a₁a₂b₂+2a₁a₃b₁+a₂²b₁)+1/6 D³⁡ϕa.ψb(3a₁²a₃+3a₁a₂²)+ϕa.1/24 D⁴⁡ψb.4b₁³b₂ +ϕa.1/120 D⁵⁡ψb.b⁵+D⁡ϕa.1/6 D⁴⁡ψb(3a₁b₁²b₂+a₂b₁³) +D⁡ϕa.1/24 D⁴⁡ψb.a₁b₁⁴+1/2 D²⁡ϕa.1/2 D²⁡ψb(2a₁²b₁b₂+2a₁a₂b₁²)+ 1/2 D²⁡a.1/6 D³⁡ψb.a₁²b₁³+1/6 D³⁡ϕa.D⁡ψb(a₁³b₂+3a₁²a₂b₁) +1/6 D³⁡ϕa.1/6 D²⁡ψb.a₁³b₁²+1/24 D⁴⁡ϕa.ψb.4a₁³a₂+1/24 D⁴⁡ϕa.D⁡ψb.a₁⁴b₁+1/120 D⁵⁡ϕa.ψb.a₁⁵

et ainsi de suite.

En examinant la composition successive de ces coefficiens, on en conclut la règle pratique suivante, pour déduire immédiatement un coefficient quelconque de celui qui le précède.

15. Pour déduire le développement de Aₙ₊₁ de celui de Aₙ; les dérivées des fonctions étant disposées en colonnes, d’après les dimensions de leurs exposens, et les lettres d’après leur ordre de succession ;

1.ᵒ On ne fera varier, dans chaque terme de chaque colonne, que les coefficiens composés des quantités polynômiales a₁, a₂, a₃, …, b₁, b₂, b₃, …, d’après la règle du n.ᵒ 8 ; en observant, pour ceux qui contiennent à la fois des a et des b, de ne faire varier d’abord que les b, et ensuite les a, mais dans le dernier terme seulement de chaque coefficient.

2.ᵒ On fera varier de plus, mais dans la dernière colonne seulement, la fonction ψb, dans tous les termes, et, comme la puissance de b₁ augmente alors d’une unité, on divisera par son exposant ainsi augmenté ;

3.ᵒ Enfin, on fera encore varier, mais dans le dernier terme de la dernière colonne seulement, la fonction ϕa; et, comme la puissance de a₁ augmente alors d’une unité, on divisera par son exposant ainsi augmenté.

Donnons des exemples de chacune des trois parties de cette règles.

1.ᵒ Le coefficient de 1/2 D²⁡ϕa.D⁡ψb, dans A₄, est a₁²b₂+2a₁a₂b₁. Pour en déduire celui du même terme dans A₅, je fais d’abord varier les b, ce qui donne a₁²b₃+2a₁a₂b₁; faisant ensuite varier les a dans le dernier terme 2a₁a₂b₁, on a, d’après la première partie de la règle du n.ᵒ 8, 2a₁a₃b, et d’après la seconde partie de cette règle a₂²b₁. Rassemblant tous ces termes, on a le coefficient de 1/2 D²⁡ϕa.ψb dans A₅.

2.ᵒ En appliquant la seconde partie de la règle ci-dessus aux cinq termes de la dernière colonne de A₄, on obtient les cinq premiers termes de la dernière colonne de A₅ ;

3.ᵒ Enfin, en appliquant la troisième partie de la règle ci-dessus au dernier terme 1/24 D⁴⁡ϕa.ψb.a₁⁴ de la dernière colonne de A₄, on obtient le dernier terme 1/120 D⁵⁡ϕa.ψb.a₁⁵ de la dernière colonne de A₅.

Cette règle est encore d’une exécution très-facile, et si expéditive qu’on peut écrire de suite, et sans s’arrêter, les termes successifs du développement. Elle n’est, jusqu’à présent, de même que celle du n.ᵒ 8, qu’une conclusion d’induction ; mais nous la démontrerons complètement dans la suite, et nous donnerons aussi une règle très-simple, pour écrire immédiatement un terme quelconque du développement, indépendamment de ceux qui le précèdent.

16. Remarque I. En examinant la composition des termes successifs (18) du développement de l’équation (15), on découvre la loi remarquable suivante qui y règne. Le terme général Aₙ est composé de n colonnes, ordonnées selon les dimensions des exposans des dérivées de ϕa et ψb, de manière que la m.ᵐᵉ colonne contient les m+1 termes

ϕa.Dᵐ⁡ψb/1.2…m, D⁡ϕa.Dᵐ⁻¹⁡ψb/1.2…(m−1), … Dᵐ⁻¹⁡ϕa/1.2…(m−1).D⁡ψb, Dᵐ⁡ϕa/1.2…m ψb.

Chacun de ces termes a pour coefficient une fonction des quantités polynômiales a₁, a₂, a₃, …, b₁, b₂, b₃, …, dont voici la formation : en supposant r+s=m, le coefficient du terme Dʳ⁡ϕa/1.2…r.Dˢ⁡ψb/1.2…s est composé de tous les produits de m lettres, dont un nombre r des quantités polynômiales en a et un nombre b des quantités polynômiales en b, de manière qua la somme de tous les indices de chaque produit soit égale à n. Quant aux coefficiens numériques de chaque produit, on les obtient en multipliant l’un par l’autre le nombre qui indique celui des permutations qu’on peut faire, entre les quantités polynômiales en a, et le nombre qui indique celui des permutions qu’on peut faire entre les quantités polynômiales en b.

Ainsi, le coefficient de 1/2 D²⁡ϕa.1/6 D³⁡ψb, dans A, est

3a₁²b₁b₂²+3a₁²b₁b₃+6a₁a₂b₁²b₃+2a₁a₃b₁³+a₂²b₁³, c’est-à-dire,

3a₁a₁b₁b₂b₂+3a₁a₁b₁b₁b₃+2.3a₁a₂b₁b₁b₃+2a₁a₃b₁b₁b₁+a₂a₂b₁b₁b₁,

qui contient en effet tous les produits possibles de deux quantités polynômiales en a et de trois en b, de manière que la somme des indices soit égale à 7 ; et qui a des coefficiens numériques qui suivent la loi que nous venons d’indiquer.

On pourrait donc, au moyen de cette loi, qui est d’ailleurs la même pour une fonction quelconque de deux polynômes indépendans, former immédiatement un terme quelconque du développement, par la théorie des combinaisons ; ce qui donnerait une extension considérable à l’analise combinatoire de Hindenburg ; mais le moyen que nous donnerons par la suite sera à la fois plus simple, plus direct et plus, analitique.

On remarquera sans doute que la simplicité de l’énoncé de cette loi, ainsi que de celle du n.ᵒ 10, n’est due qu’au choix que nous avons fait d’indices numériques, pour représenter les quantités polynômiales ; elle n’aurait pu s’énoncer que très-difficilement, avec les lettres dans l’ordre alphabétique, employées par Arbogast et Hindenburg. Ces lettres à indices ont encore l’avantage d’indiquer, de la manière la plus caractéristique, leurs relations avec les dérivées du premier terme du polynôme, puisqu’on a généralement aₙ=Dⁿ⁡a/1.2…n : tant il est vrai que souvent le plus léger changement dans les notations peut avoir l’influence la plus heureuse sur les méthode.

17. Remarque II. La théorie que nous venons d’exposer, contient tout ce qui est nécessaire pour le développement complet des fonctions d’un polynôme, et même, à la rigueur, pour celui d’une fonction quelconque de deux polynômes indépendans ; car il suffirait, pour le développement de ϕ(a+a₁x+a₂x²+…, b+b₁x+b₂x²+…), de remplacer, dans les formules (18), les produits tels que Dʳ⁡ϕa/1.2…r.Dˢ⁡ψb/1.2…s par les dérivées partielles du même ordre de ϕ(a, b). Mais, nous allons encore envisager cette théorie sous un autre point de vue, qui nous facilitera singulièrement l’exposition de celle du retour des fonctions et des séries, à laquelle nous nous proposons de consacrer l’article II.

18. On a, par le n.ᵒ I,

ϕ(a+y)=ϕa+D⁡ϕa.y+1/2 D²⁡ϕa.y²+1/6 D³⁡ϕa.y³+…

Si l’on suppose

y=a₁x+a₂x²+a₃x³+…=x(a₁+a₂x+a₃x²+…),

où les coefficiens a₁, a₂, a₃, … représentent des quantités quelconques et indépendantes, et qu’on substitue cette valeur de y dans l’équation (19), son premier membre se transformera en celui de l’équation (7) ; on aura donc

ϕ(a+y)=ϕ(a+a₁x+a₂x²+a₃x³+…)

=ϕa+D⁡ϕax+1/2 D².ϕax²+1/6 D³.ϕax³+…

Ces deux développements ne diffèrent l’un de l’autre qu’en ce que, dans les dérivées du premier, on suppose D⁡a=1, D²⁡a=0, D³⁡a=0, … et dans, celle du second, D⁡a=a₁, D²⁡a=2a₂, D³⁡a=6a₃, …. La manière de déduire les dérivées suivies d’un point, de celles sans point a été exposée aux n.ᵒˢ 7 et 8.

L’équation (21) fournit donc le moyen de résoudre cette question : y étant une fonction donnée de x, ou un polynôme en x ; développer, selon les puissances de x, une fonction quelconque ϕ(a+y) ? En effet, d’après le n.ᵒ 3, l’équation (20) peut être mise sous la forme

y=xf⁡(α+x),

et l’on a

a₁=f⁡α, a₂=D⁡f⁡α, a₃=1/2 D²⁡f⁡α, a₄=1/6 D³⁡f⁡α, …;

c’est-à-dire, que a₁ doit être considéré comme un premier terme de polynôme.

En substituant ces valeurs dans l’équation (21), on obtient

ϕ{a+xf⁡(α+x)}=ϕ(a+f⁡α.x+D⁡f⁡α.x²+1/2 D²⁡f⁡α.x³+…)

=ϕa+D.ϕa.x+1/2 D².ϕa.x²+1/6 D².ϕa.x³+…,

où, dans le développement du dernier membre, qu’on exécute d’après la règle du n.ᵒ 8, il faut substituer, pour a₁, a₂, a₃, … leurs valeurs (23).

19. Si, dans la question du n.ᵒ précédent, la valeur de y était donnée par l’équation suivante :

y=xψf⁡(α+x)

qui, d’après le n.ᵒ 5 devient

y=x(ψf⁡α+D.ψf⁡α.x+1/2 D².ψf⁡α.x²+1/6 D³.ψf⁡α.x³+…)

il faudrait faire, dans le développement du dernier membre de l’équation (24),

a₁=ψf⁡α, a₂=D.ψf⁡α, a₃=1/2 D².ψf⁡α, a₄=1/6 D³.ψf⁡α, …;

mais, conformément aux principes des n.ᵒˢ 5 et 6, ces dernières dérivées doivent être survies d’un point, et développées d’après le n.ᵒ 7.

20. Si l’on avait à développer, selon les puissances de x ; la fonction ϕ[a+x), z étant donné par l’équation

z=xψ(β+y),

et y par l’équation (22) ; on aurait

ϕ(a+z)=ϕ[a+xψ(β+y)]=ϕ{a+xψ[β+xf⁡(α+x)]}

=ϕ{a+xψ(β+f⁡α.x+D⁡f⁡α.x²+1/2 D²⁡f⁡α.x³+…)}

=ϕa+D.ϕa.x+1/2 D².ϕa.x²+1/6 D³.ϕa+…;

mais ici, dans le développement des dérivées du dernier membre, il faudrait remplacer a₁, a₂, a₃, … par ψβ, D.wβ, 1/2 D².ψβ, …, en observant de mettre, dans le développement de ces dernières dérivées, f⁡α à la place de D⁡β, D⁡f⁡α à la place de 1/2 D²⁡β, 1/2 D²⁡f⁡α à la place de 1/6 D³⁡β; et ainsi de suite.

On pourrait aisément pousser plus loin ces substitutions de fonctions dans les fonctions, ou de séries dans les séries ; et l’on voit que le principe de leur développement par les dérivations est simple et uniforme : il ne reste que la complication des résultats, qui est inhérente à la chose même.

### Article II.

21. Depuis le n.ᵒ 18 de l’article précédent, nous nous sommes occupés de la question suivante : le développement d’une fonction quelconque, selon les puissances d’une fonction donnée de la variable principale, étant supposé connu ; en déduire le développement selon les puissances de la variable principale ? Dans cet article, nous allons résoudre la question inverse, savoir : le développement d’une fonction quelconque, selon les puissances de la variable principale, étant donné, ainsi que la relation entre celle variable et une autre fonction ; en déduire le développement selon, les puissances de cette dernière fonction, Cette question contient le problème général du retour des fonctions et des séries.

22. Proposons-nous de transformer le polynôme

A+A₁x+A₂x²+A₃x³+…,

procédant selon les puissances de la variable principale x, en un polynôme

B+B₁y+B₂y²+B₃y³+…,

procédant selon les puissances de y, dont la valeur est supposée donnée par l’équation (20)

y=x(a₁+a₂x+a₃x²+…).

En comparant le polynôme (30) avec l’équation (21), et le polynôme (31) avec l’équation (19), on obtient

A=ϕa, A₁=D.ϕa, A₂=1/2 D².ϕa, A₃=1/6 D³.ϕa, …;

B=ϕa, B₁=D⁡ϕa, B₂=1/2 D²⁡ϕa, B₃=1/6 D³⁡ϕa, …;

Ici, ce sont les dérivées suivies d’un point qui sont données immédiatement ; et la question se réduit à en déduire celles sans points. On pourrait la résoudre en tirant les valeurs de ces dernières des équations (10), par des éliminations successives ; mais, outre que ce moyen serait trop long, il est peu propre à faire découvrir la loi qui y règne : il est bien plus simple de les former immédiatement de la manière suivante.

On a, d’après les n.ᵒˢ 6 et 7, D.ϕa=D⁡ϕa.a₁, et par conséquent D⁡ϕa=a₁⁻¹D.ϕa; donc, en répétant l’opération indiquée par cette équation, on obtient

{D⁡ϕa=a₁⁻¹D.ϕa, D²⁡ϕa=D⁡(D⁡ϕa)=a₁⁻¹D.(a₁⁻¹D.ϕa), D³⁡ϕa=D⁡(D²⁡ϕa)=a₁⁻¹D.[a₁⁻¹D.(a₁⁻¹D.ϕa)], ……………………;

mais ici il faut observer que la manière dont nous sommes parvenus à ces relations suppose que a₁ est un second terme de polynôme, c’est-à-dire que, d’après le n.ᵒ 1, on a

a₁=D⁡a, a₂=1/2 D²⁡a, a₃=1/6 D³⁡a, …aₙ=1/1.2…n Dⁿ⁡a;

et par conséquent

D⁡a₁=D²⁡a=2a₂, D²⁡a₁=D³⁡a=6a₃, …

Dⁿ⁻¹⁡a₁=Dⁿ⁡a=1.2…n.aₙ;

ainsi, dans le développement des seconds membres des équations (34), il faudra substituer pour les dérivées de a₁ leurs valeurs précédentes.

Mais si, conformément au n.ᵒ 18, on veut considérer a₁ comme premier terme de polynôme, on a

a₂=D⁡a₁, a₃=1/2 D²⁡a₁, a₄=1/6 D³⁡a₁, …aₙ=1/1.2…n Dⁿ⁻¹⁡a₁:

en substituant ces valeurs dans les développemens des équations (34), ce qui revient à y écrire 2D⁡a₁ pour D⁡a₁, 3D²⁡a₁ pour D²⁡a₁, 4D³⁡a₁ pour D³⁡a₁, nDⁿ⁻¹⁡a₁ pour Dⁿ⁻¹⁡a₁, on pourra mettre ces équations sous la forme suivante ;

{D⁡ϕa=a₁⁻¹D.ϕa, D²⁡ϕa=D.(a₁⁻²D.ϕa), D³⁡ϕa=D².(a₁⁻³D.ϕa), …………………………Dⁿ⁡ϕa=Dⁿ⁻¹.(a₁⁻ⁿD.ϕa),

En substituant ces valeurs dans les équations (33), et remplaçant ϕa par A, ou obtient enfin

{B=A, B₁=a₁⁻¹D.A, B₂=1/2 D.(a₁⁻²D.A), B₃=1/6 D².(a₁⁻³D.A), …………………………Bₙ=1/1.2…n Dⁿ⁻¹.(a₁⁻ⁿD.A);

où, d’après l’observation du n.ᵒ 18, et d’après l’observation précédente, a₁ doit être considéré comme un premier terme de polynôme.

23. Si l’on fait attention que l’équation (20) peut être mise sous la forme (22), et que le polynôme (30), d’après le n.ᵒ 3, peut représenter une fonction quelconque ϕ(b+x) ; le problème du n.ᵒ précédent fournit la solution de la question suivante : étant donnée la relation y=xf⁡(α+x), développer la fonction quelconque, ϕ(b+x) suivant les puissances de y.

D’après cela, si l’on substitue, dans le polynôme (31), les valeurs (36), et dans celle-ci pour A et a₁ leurs valeurs ϕb et f⁡α, on aura

ϕ(b+x)=B+B₁y+B₂y²+B₃y³+…

=ϕb+(f⁡α)⁻¹D.ϕby+1/2 D.{(f⁡α)⁻²D.ϕb}y²

+1/6 D².{(f⁡α)⁻³D.ϕb}y³+…

où l’on peut supprimer, si l’on veut, les points qui suivent les signes de dérivation qui affectent ϕb ; car, dans le binôme b+x, on a D.b=1; et par conséquent D.ϕb=D⁡ϕb.

On aurait de même, dans la même hypothèse,

ϕ(b+b₁x+b₂x²+b₃x³+…)=B+B₁y+B₂y²+B₃y³+…

=ϕb+(f⁡α)⁻¹D.ϕb.y+1/2 D.{(f⁡α)⁻²D.ϕb}y²

+1/6 D².{(f⁡α)⁻³D.ϕb}y³+…

mais ici les points, après tous les signes de dérivation, sont indispensables, car on a D⁡b=b₁, D²⁡b=2b₂, D³⁡b=6b₃, …, et par conséquent D.ϕb=D⁡ϕb.b₁, …

On aurait encore, de la même manière, et pour la même valeur de y,

ψϕ(b+x)=ψϕb+(f⁡α)⁻¹D.ψϕb.y+1/2 {(f⁡α)⁻²D.ψϕb}y²

+1/6 {(f⁡α)⁻³D.ψϕb}y³+…,

où les points, après les signes de dérivation sont encore nécessaires : parce qu’on a D.ψϕb=D⁡ψϕb.D⁡ϕb.

24. Si, dans la question du n.ᵒ précédent, la valeur de y était donnée par l’équation (25)

y=xψf⁡(α+x),

il faudrait substituer, dans les équations (36), pour a₁ sa valeur ψf⁡α, conformément au n.ᵒ 19 ; ce qui donnerait, d’après l’observation faite sur l’équation (37)

ϕ(b+x)=ϕb+(ψf⁡α)⁻¹D⁡ϕb.y+1/2 D.{(ψf⁡α)⁻²D⁡ϕb}y²

+1/6 D².{(ψf⁡α)⁻³D⁡ϕb}y³+…

Dans le cas particulier où ψ, représente la puissance −1, l’équation (25) devient

y=x/(f⁡(α+x)),

et alors l’équation (40) se change en

ϕ(b+x)=ϕb+f⁡α.D⁡ϕb.y+1/2 D.{(f⁡α)²D⁡ϕb}y²

+1/6 D².{(f⁡α)³D⁡ϕb}y³+…

On aurait de même, pour la valeur de y (25),

F⁡ϕ(b+x)=F⁡ϕb+(ψf⁡α)⁻¹D.F⁡ϕb.y

+1/2 D.{(ψf⁡α)⁻²D.F⁡ϕb}y²+1/6 D².{(ψf⁡α)⁻³D.F⁡ϕb}y³+…

où la même observation n’a lieu qu’après l’équation (39).

25. Proposons-nous enfin de résoudre la question suivante : étant données les relations

y=xf⁡(α+x), z=xψf⁡(β+y),

développer la fonction ϕ(b+x) selon les puissances de z, sans x ni y.

En comparant les solutions des n.ᵒˢ précédens avec la question du n.ᵒ 20, dont celle-ci est l’inverse, on obtient immédiatement

ϕ(b+x)=ϕb+(ψβ)⁻¹D⁡ϕb.z+1/2 D.{(ψβ)⁻²D⁡ϕb}z²

+1/6 D².{(ψβ)⁻³D⁡ϕb}z³+…;

en observant seulement de mettre, dans le développement des dérivées de ψβ, f⁡α pour D⁡β, Df⁡α pour 1/2 D²⁡β, 1/2 Df⁡α pour 1/6 D³⁡β, et ainsi de suite.

En se conformant à cette observation, on aurait de même

F⁡ϕ(b+x)=F⁡ϕb+(ψβ)⁻¹D.F⁡ϕb.z+1/2 D.{(ψβ)⁻²D.F⁡ϕb}z²

+1/6 D².{(ψβ)⁻³D.F⁡ϕb}z³+…

26. Remarque. La question traitée au n.ᵒ précédent est une espèce de retour double : on pourrait en former de pareilles sur des retours triples, quadruples, etc. : le principe de leurs solutions se déduit aisément de celle du n.ᵒ précédent ; et leur développement par les dérivations s’exécuterait aussi facilement que leur complication naturelle peut le permettre.

27. Depuis le commencement de cet article, nous n’avons fait qu’établir les formules générales du retour des fonctions et des séries ; occupons-nous maintenant de leur développement complet et effectif. Reprenons, à cet effet, les problèmes du n.ᵒ 23, et proposons-nous de développer complètement les coefficiens successifs B, B₁, B₂, B₃, … de l’équation (38).

Comme nous avons vu, aux n.ᵒˢ 18 et 22, que a₁ devait être considéré comme un premier terme de polynôme, dans l’équation (20) ou (22), et que d’ailleurs les quantités a, a₁, a₂, a₃, … peuvent être quelconques ; nous les remplacerons par c, c₁, c₂, c₃, …, afin de conserver la régularité dans les développemens ; ainsi, l’équation (20) ou (22) deviendra

y=xf⁡(α+x)=x(c+c₁x+c₂x²+c₃x³+…)

Au moyen de cette observation, le problème en question se réduit à développer les termes B, B₁, B₂, B₃, … des équations (36), en y substituant c à la place de a₁, et de ϕb à la place de A; ce qui donne

1B₁=c⁻¹D.ϕb, 2B₂=D.(c⁻²D.ϕb),

3B₃=1/2 D².(c⁻³D.ϕb)…

En comparant ces termes avec la formule (17), on voit aisément que leur développement doit s’exécuter par la même règle, en observant cependant qu’ici la fonction ϕa est remplacée par une puissance négative de c, dont l’exposant est égal à l’indice du terme ; et que la fonction ψb est remplacée par D.ϕb. Avec cette attention, on aura, en suivant la règle du n.ᵒ 15, les développemens suivans, analogues à ceux (18)

{B=ϕb, 1B₁=c⁻¹.D⁡ϕb.b₁, 2B₂=c⁻².D⁡ϕb.b₂+c⁻².DD⁡ϕb.b₁²+D⁡(c⁻²).D⁡ϕb.c₁b₁3B₃=c⁻³.D⁡ϕb.b₃+c⁻³.DD⁡ϕb.2b₁b₂+c⁻³.D²⁡D⁡ϕb.b₁³/2.b₁³+D⁡(c⁻³).D⁡ϕb(c₁b₂+c₂b₁)+D⁡(c⁻³).DD⁡ϕb.c₁b₁²+D²⁡(c⁻³/2 D⁡ϕb.c₁²b₁4B₄=c⁻⁴.D⁡ϕb.b₄+c⁻⁴.DD⁡ϕb(2b₁b₃+b₂²)+c⁻⁴.D²D⁡ϕb/2.3b₁²b₂+D⁡(c⁻⁴).D⁡ϕb(c₁b₃+c₂b₂+c₃b₁)+D⁡(c⁻⁴).DD⁡ϕ(2c₁b₂+c₂b₁²)+D²⁡(c⁻²/2.D⁡ϕb(c₁²b₂+2c₁c₂b₁)+c⁻⁴.D³D⁡ϕb/6.b₁⁴+D⁡(c⁻⁴).D²D⁡ϕb/2.c₁b₁³+D²⁡(c⁻⁴)/2.DD⁡ϕb.c₁²b₁²+D³⁡(c⁻⁴)/6.D⁡ϕb.c₁³b₁5B₅=c⁻⁵.D⁡ϕb.b₅+c⁻⁵.DD⁡ϕb(2b₁b₄+2b₂b₃)+c⁻⁵.D²D⁡ϕb/2.(3b₁²b₃+3b₁b₂²)+D⁡(c⁻⁵).D⁡ϕb(c₁b₄+c₂b₃+c₃b₂+c₄b₁)+D⁡(c⁻⁵).DD⁡ϕb.(2c₁b₁b₃+c₁b₂²+2c₂b₁b₂+c₃b₁²)+D²⁡(c⁻⁵)/2.D⁡ϕb.(c₁²b₃+2c₁c₂b₂+2c₁c₃b₁+c₂²b₁)+c⁻⁵.D³D⁡ϕb/2.4b₁³b₂++c⁻⁵.D⁴⁡D⁡ϕb/24.b₁⁵+D⁡(c⁻⁵).D²D⁡ϕb/2 (3c₁b₁²+c₂b₁³)+D⁡(c⁻⁵).D³D⁡ϕb/6.c₁b₁⁴

{+D²⁡(c⁻⁵)/2.DD⁡ϕb(2c₁²b₁b₂+2c₁c₂b₁²)+D²⁡(c⁻⁵)/2.D²D⁡ϕb/2.c₁²b₁³+D³⁡(c⁻⁵)/6.D⁡ϕb(c₁³b₂+3c₁²c₂b₁)+D³⁡(c⁻⁵)/6.DD⁡ϕb.c₁³b₁²+D⁴⁡(c⁻⁵)/24.D⁡ϕb.c₁⁴b₁

et ainsi de suite.

Si l’on effectue les dérivations des puissances négatives de c, qui ne sont qu’indiquées, ainsi que celles de D⁡ϕb, et qu’on ordonne selon les dérivées de ϕb, on obtient

(50){B=ϕb1B₁=D⁡ϕb.c⁻¹.b₁2B₂=D⁡ϕb|c⁻².b₂+2 D²⁡ϕb/2.c⁻².b₁²−2c⁻³.b₁c₂3B₃=D⁡ϕb|c⁻³.b₃+2 D²⁡ϕb/2 −3c⁻⁴(b₁c₂+b₂c₁)+3.4/2.c⁻⁵.b₁c₁²|c⁻³.2b₁b₂+3 D³⁡ϕb/6.c⁻³.b₁³−3c⁻⁴.b₁²c₁4B₄=D⁡ϕb|c⁻⁴.b₄+2 D²⁡ϕb/2 −4c⁻⁵.(b₁c₃+b₂c₂+b₃c₁)+4.5/2.c⁻⁶(2b₁c₁c₂+b₂c₁²)−4.5.6/2.3.c⁻⁷.b₁c₁³|c⁻⁴(2b₁b₃+b₂²)−4c⁻⁵(b₁²c₂+2b₁b₂c₁)+4.5/2.c⁻⁶.b₁²c₁²+3 D³⁡ϕb/6 |c⁻⁴.3b₁²b₂+4 D⁴⁡ϕb/24.c⁻⁴.b₁⁴−4c⁻⁵.b₁³c₁

{5B₅=D⁡ϕb|c⁻⁵.b₅+2 D²⁡ϕb/2 −5c⁻⁶(b₁c₄+b₂c₃+b₃c₂+b₄c₁)+5.6/2 c⁻⁷(2b₁c₁c₃+b₁c₂²+2b₂c₁c₂+b₃c₁²)−5.6.7/2.3.c⁻⁸(3b₁c₁c₃+b₂c₁³)+5.6.7.8/2.3.4.c⁻⁹.b₁c₁⁴|c⁻⁵(2b₁b₄+2b₂b₃)−5c⁻⁶(b₁²c₃+2b₁b₂c₂+2b₁b₃c₁+b₂²c₁)+5.6/2.c⁻⁷(2b₁²c₁c₂+2b₁b₂c₁²)−5.6.7/2.3.c⁻⁸.b₁²c₁³+3 D³⁡ϕb/6 |c⁻⁵(3b₁²b₃+3b₁b₂²)+4.D⁴⁡ϕb/24 −5c⁻⁶(b₁³c₂+3b₁²b₂c₁)+5.6/2 c⁻⁷.b₁³c₁²|c⁻⁵.4b₁³b₂+5.D⁵⁡ϕb/120.c⁻⁵.b₁⁵−5c⁻⁶.b₁⁴c₁

et ainsi de suite.

L’examen de la composition successive des termes fournit encore une règle pratique, pour déduire un terme quelconque de celui qui le précède.

28. Pour déduire le développement de (n+1)Bₙ₋₁ de celui de nBₙ, celui-ci étant ordonné en colonnes, par rapport aux dérivées successives de ϕb, les termes de chaque colonne, par rapport aux puissances de c, et les quantités polynômiales d’après leur ordre de succession ;

1.ᵒ On divisera tous les termes de nBₙ par n, on multipliera chacun par l’exposant de c dans ce terme (abstraction faite du signe), et l’on augmentera cet exposant d’une unité (aussi abstraction faite de son signe) ;

2.ᵒ On ne fera varier, dans chaque ferme de chaque colonne, que les coefficiens composés des quantités polynômiales b₁, b₂, b₃, …, c₁, c₂, c₃, …, d’après la règle du n.ᵒ 8 ; en observant, pour ceux qui contiennent à la fois des b et des c, de ne faire varier d’abord que les c, et ensuite les b, mais dans le dernier terme seulement de chaque coefficient.

3.ᵒ On fera varier de plus, mais dans le dernier terme seulement de chaque colonne, la puissance de c ; et, comme la puissance de c₁ augmente alors d’une unité, on divisera par son exposant ainsi augmenté ;

4.ᵒ Enfin, on fera varier ϕb, dans le tout dernier terme seulement, en mettant (n+1) (Dⁿ⁺¹⁡ϕb)/(1.2…(n+1)) pour n Dⁿ⁡ϕb/1.2…n, et augmentant la puissance de b₁ d’une unité.

Cette règle est analogue à celle du n.ᵒ 15 : dans l’exécution, on n’a pas besoin de faire d’avance la préparation de la première partie ; elle peut se faire à mesure qu’on opère sur chaque terme.

29. Au moyen de la règle précédente, on peut écrire de suite les termes successifs du développement de l’équation (38) tout ordonnés et réduits à leur plus simple expression. Si l’on suppose b₁=1, b₂=0, b₃=0, …, on aura le cas de l’équation (37) ; et il n’en résulte d’autre changement à la règle précédente qu’une simplification dans la seconde partie, parce qu’il n’y a plus que des quantités polynômiales d’une seule espèce ; ainsi, les formules (50) deviendront, pour ce cas (en remplaçant les colonnes par des parenthèses),

B= ϕb1B₁=D⁡ϕb.c⁻¹2B₂=D⁡ϕb(−2c⁻³.c₁)+2 D²⁡ϕb/2.c⁻²3B₃=D⁡ϕb(−3c⁻⁴.c₂+3.4/2.c⁻⁵.c₁²)+2 D²⁡ϕb/2 (−3c⁻⁴.c₁)+3 D³⁡ϕb/6.c⁻⁵, 4B₄=D⁡ϕb(−4c⁻⁵.c₃+4.5/2.c⁻⁶.2c₁c₂−4.5.6/6.c⁻⁷.c₁³)+2 D²⁡ϕb/2 (−4c⁻⁵.c₂+4.5/2.c⁻⁶.c₁²)+3 D³⁡ϕb/6 (−4c⁻⁵.c₁)+4.D⁴⁡ϕb/24.c⁻⁴,

5B₅=D⁡ϕb[−5c⁻⁶.c₄+5.6/2.c⁻⁷(2c₁c₃+c₂²)−5.6.7/6.c⁻⁸.3c₁²c₂+5.6.7.8/24.c⁻⁹.c₁⁴]+2 D²⁡ϕb/2 (−5c⁻⁶.c₂+5.6/2.c⁻⁷.2c₁c₂−5.6.7/2.3.c⁻⁸.c₁³)+3 D³⁡ϕb/6 (−5c⁻⁶.c₂+5.6/2.c⁻⁷.c₁²c²)+4 D⁴⁡ϕb/24 (−5c⁻⁶.c₁)+5.D⁵⁡ϕb/120 c⁻⁵

et ainsi de suite.

Pour le développement de l’équation (39) ; comme on a, d’après le n.ᵒ 5,

ψϕ(b+x)=ψ(ϕb+D.ϕb.x+1/2 D²⁡ϕb.x²+1/6 D³.ϕb.x³+…),

la règle reste la même ; mais, au lieu de ϕb, il faut écrire ψϕb, D⁡ϕb, au lieu de b₁, 1/2 D²⁡ϕb, au lieu de b₂, 1/6 D³⁡ϕb, au lieu de b₃, et ainsi de suite.

Pour le cas de l’équation (40), comme ψf⁡α=ψc, il faudrait, en conservant la même règle, mettre partout ψc à la place de c, D⁡ψc à la place de c₁, 1/2 D²⁡ψc à la place de c₂, 1/6 D³⁡ψc à la place de c₃, …; en observant que ces dérivées doivent être elles-mêmes développées selon la règle du n.ᵒ 8 ; que D.ψc=D⁡ψc.c₁=D⁡ψf⁡α.Df⁡α, …, et que, dans ces derniers développemens, il faut substituer f⁡α pour c, Df⁡α pour c₁, 1/2 D²f⁡α pour c₂, 1/6 D³f⁡α pour c₃, ….

Pour le cas de l’équation (43), il faudrait tenir compte, à la fois, des deux observations précédentes, et écrire F⁡ϕb pour ϕb D⁡ϕb pour b₁ et ψc pour c D.ψc pour c₁, ….

Mais, pour l’équation (42), la règle du n.ᵒ précédent s’emploie sans la moindre restriction, parce que cette équation ne diffère de celle (37) que par le signe des exposans de c ou f⁡α dont cette règle est indépendante.

Pour le développement de l’équation (45), il faudrait remplacer c, c₁, c₂, …, par ψβ, D.ψβ, 1/2 D².ψβ, …, en observant que D⁡β=c, 1/2 D²⁡β=c₁, 1/6 D³⁡β=c₂, …,

Enfin, pour le développement de l’équation (46), il faudrait tenir compte de l’observation précédente, et de plus mettre F⁡ϕb pour ϕb, D⁡ϕb pour b₁, 1/2 D²⁡ϕb pour b₂, ….

Au moyen de ces observations, l’application de la règle du n.ᵒ précédent est générale.

30. Remarque. En effectuant les dérivations de D⁡ϕb, indiquées dans les équations (49), on obtient D²⁡ϕb, D³⁡ϕb/2, D⁴⁡ϕb/6, …: nous avons préféré, dans les équations (50) et (51), d’écrire, à la place de ces résultats, 2 D²⁡ϕb/2, 3 D³⁡ϕb/6, 4 D⁴⁡ϕb/24, … parce que D⁡ϕb, 1/2 D²⁡ϕb, 1/6 D³⁡ϕb, 1/24 D⁴⁡ϕb, … sont les coefficiens du développement de ϕ(b+x), et que, par ce moyen, les coefficiens numériques sont mis en évidence : ainsi, pour le problème du n.ᵒ 22, on a D⁡ϕb=A₁, D²⁡ϕb=2A₂, D³⁡ϕb/2=3A₃, D⁴⁡ϕb/6=4A₄, ….

Nous avons déjà remarqué au n.ᵒ 27 que, d’après les n.ᵒˢ 18 et 22, a₁ devait être considéré comme premier terme de polynôme, et par conséquent comme indépendant de a ; c’est pourquoi, dès le n.ᵒ 23, nous avons remplacé partout cette lettre par b, sous les signes de fonction, afin de ne pas induire en erreur, par une prétendue dépendance qui n’existait plus. Cette observation deviendra encore plus claire par la théorie de l’article suivant.

C’est pour la même raison, et pour conserver la régularité de la loi des développemens, que nous avons remplacé, au n.ᵒ 27, le polynôme a₁+a₂x+a₃x²+… par celui c+c₁x+c₂x²+… Si, au n.ᵒ 18, nous avons préféré la première de ces deux formes, ce n’était que pour mieux faire apercevoir l’identité des développemens de ϕ(a+a₁x+a₂x²+…) et de ϕ{a+x(a₁+a₂x+…)}, et pour rendre plus palpable la dépendance mutuelle des coefficiens des développemens de ϕ(a+y) et de ϕ(a+a₁x+a₂x²+…); dépendance qui nous a tant simplifié l’exposition de la théorie du retour des suites. On aura remarqué sans doute que la loi de cette dépendance est la même que celle du changement de la variable principale, dans la différentiation d’une fonction de deux variables.

### Article III.

31. Les règles des n.ᵒˢ 8, 15 et 28 ne sont que des conclusions d’induction, tirées de l’examen de la formation successive des termes d’un développement ; et, sous ce rapport, elles peuvent laisser-quelque doute sur l’exactitude des résultats qu’elles fournissent. Il est donc nécessaire de démontrer ces règles, afin que le calcul de dérivations soit non seulement un instrument commode et expéditif, mais encore sûr et rigoureux.

Ces mêmes règles n’offrent que le moyen de former successivement les termes du développement, en déduisant chacun de celui qui le précède ; de sorte que, pour avoir, par exemple, le vingtième terme du développement, il faut calculer auparavant les dix-neuf qui sont à sa gauche. Mais souvent on n’a besoin que d’un terme assez éloigné de l’origine du développement pour que le calcul préalable de tous ceux qui le précèdent exige une perte de temps aussi considérable qu’inutile à l’objet qu’on a en vue. Il est donc essentiel d’avoir le moyen de former immédiatement un terme quelconque, indépendamment de tous ceux qui seraient avant lui.

Tels sont les deux objets que nous nous proposons de remplir dans cet article.

32. Si, dans le polynôme a+a₁x+a₂x²+a₃x³+…, on suppose a₂=0, a₃=0, …, l’équation (7) deviendra

ϕ(a+a₁x)=ϕa+D.ϕa.x+1/2 D².ϕa.x²+1/6 D³.ϕa.x³+…;

mais, par la supposition que nous venons de faire, on a [équations (3)], D⁡a=a₁, D²⁡a=0, D³⁡a=0, …; ce qui donne, d’après les n.ᵒˢ 6 et 7, D.ϕa=D⁡ϕa.a₁, D².ϕa=D²⁡ϕa.a₁², D³.ϕa= D³⁡ϕa.a₁³, … En substituant ces valeurs dans l’équation (25), on obtient

ϕ(a+a₁x)=ϕa+D⁡ϕa.a₁x+1/2 D².ϕa.a₁²x²

+1/6 D³⁡ϕa.a₁³x³+…;

résultat identique avec celui qu’on aurait obtenu en mettant a₁x à la place de x dans le théorème de Taylor.

Supposons maintenant que a₁ devienne a₁+a₂x : les puissances de a₁ se changeront en puissances de a₁+a₂x qui, étant elles-mêmes des fonctions de binôme, peuvent être développées comme les équations (52) et (53) ; mais, dans ce cas, ces formules se termineront, parce que D⁡a₁=a₂, D²⁡a₁=0, … donnent, en général, Dⁿ.a₁ⁿ/1.2…n=(D⁡a₁)ⁿ=a₂ⁿ, et Dⁿ⁺¹.a₁ⁿ=0. Substituant donc, avec cette attention, a₁+a₂x pour a₁, dans l’équation (53), on obtient

ϕ{a+x(a₁+a₂x)}=ϕa+D⁡ϕa(a₁+D⁡a₁x)x

+1/2 D².ϕa(a₁²+D.a₁²x+ 1/2D².a₁².x²)x²

+1/6 D³.ϕa(a₁³+D.a₁³x+1/2 D².a₁³.x²+1/6 D³.a₁³.x³)x³+…

En effectuant les dérivations indiquées, d’après les règles ordinaires de la différentiation, et remplaçant D⁡a₁ par a₂, cette équation devient

ϕ{a+x(a₁+a₂x)}=ϕa+D⁡ϕa(a₁+a₂x)x

+1/2 D².ϕa(a₁²+2a₁a₂x+.a₂²x²)x²

+1/6 D³.ϕa(a₁³+3a₁²a₂x+3a₁a₂².x²+a₂³.x³)x³+…

En ordonnant cette équation par rapport aux puissances de x, on obtient

ϕ(a+a₁x+a₂x²)=ϕa+D⁡ϕa.a₁x+(D⁡ϕa.a₂+1/2 D²⁡ϕa.a₁²)x²

+( 1/2 D²⁡ϕa.2a₁a₂+1/6 D³⁡ϕa.a₁³)x³

+( 1/2 D²⁡ϕa.a₂²+1/6 D³⁡ϕa.3a₁²a₂+1/24 D⁴⁡ϕa.a₁⁴)x⁴+…

Si l’on suppose ensuite que a₂ devienne a₂+a₃x, l’équation précédente deviendra, d’après les mêmes principes,

ϕ{a+a₁x+x²(a₂+a₃x)}=ϕa+D⁡ϕa.a₁x

+{D⁡ϕa(a₂+D⁡a₂.x)+1/2 D²⁡ϕa.a₁²}x²

+{ 1/2 D²⁡ϕa.2a₁(a₂+D⁡a₂.x)+1/6 D³⁡ϕa.a₁³}x³

+{ 1/2 D²⁡ϕa(a₂²+D.a₂².x+1/2 D².a₂²x²)

+1/6 D³⁡ϕa.3a₁²(a₂+D⁡a₂.x)+1/24 D⁴⁡ϕa.a₁⁴}x⁴+…

ou, en effectuant les dérivations indiquées, mettant a₃ pour D⁡a₂, et ordonnant par rapport à x,

ϕ(a+a₁x+a₂x²+a₃x³)=ϕa+D⁡ϕa.a₁x

+(D⁡ϕa.a₂+1/2 D²⁡ϕa.a₁²x²)x²

+(D⁡ϕa.a₃+1/2 D²⁡ϕa.2a₁a₂+1/6 D³⁡ϕa.a₁³)x³

+{ 1/2 D²⁡ϕa(2a₁a₃+a₂²)+1/6 D³⁡ϕa.3a₁²a₂+1/24 D⁴⁡ϕa.a₁⁴}x⁴+…

En comparant les coefficiens des seconds membres des équations (53), (54), (55) avec les formules (10), on voit que les deux premiers termes de l’équation (53), les trois premiers de (54) et les quatre premiers de (55) sont déjà complets. En continuant ces substitutions, on obtiendrait chaque fois un terme complet de plus, et l’on arriverait enfin au développement entier de la fonction de polynôme ϕ(a+a₁x+a₂x²+a₃x³+…). Mais, sans aller plus loin, nous pouvons déjà observer, 1.ᵒ qu’on ne fait jamais varier, dans chaque terme, qu’une seule lettre à la fois, ou sa puissance, et que celle lettre est la dernière dans l’ordre des indices ; car, d’après la marche que nous venons de suivre, dans ces développemens successifs, il est évident que les dernières lettres, dans l’ordre des indices, ne proviennent que des variations qu’ont subies les lettres précédentes ; or, si l’on faisait encore varier celles-ci, il en résulterait que les mêmes lettres auraient subi plusieurs variations ; ce qui est contraire à la marche de ces substitutions successives, où l’on ne fait plus attention aux lettres qui ont déjà subi une variation ; et il s’ensuit que, dans chaque terme, on ne doit faire varier que la dernière lettre ou sa puissance ; 2.ᵒ que, dans ces variations successives, chaque lettre est considérée comme un premier terme de polynôme ; c’est-à-dire, qu’on écrit a₂ pour D⁡a₁, a₃ pour D⁡a₂, a₄ pour D⁡a₃, … sans autre coefficient que l’unité.

Voilà donc les deux conditions principales de la première partie de la règle du n.ᵒ 8 justifiées. Mais examinons de plus près la formation de chaque terme du développement, en supposant que toutes les substitutions précédentes, au lieu d’être successives, soient faites à la fois.

33. Le terme Aₙ₊₁, ou le coefficient de xⁿ⁺¹ dans le développement de ϕ(a+a₁x+a₂x²+a₃x³+…) ne peut être composé que des trois parties suivantes : 1.ᵒ du terme correspondant du développement de ϕ(a+a₁x), équation (53), c’est-à-dire, de (Dⁿ⁺¹⁡ϕa)/(1.2…n+1) a₁ⁿ⁺¹, 2.ᵒ des termes provenant de la substitution de a₁+a₂x pour a₁, de a₂+a₃x pour a₂, de a₃+a₄x pour a₃,

et ainsi de suite, dans les dernières lettres (ou leurs puissances) de chaque terme de Aₙ ; 3.ᵒ enfin, de ceux provenant des mêmes substitutions, dans les puissances des dernières lettres des termes de Aₙ₋₁, Aₙ₋₂, …, en remontant. Examinons chacune de ces trois parties :

1.ᵒ La première partie a toujours évidemment lieu ; car il faut qu’elle subsiste quand a₂, a₃, a₄, … deviennent nuls ; nous verrons tout à l’heure comment la règle du n.ᵒ 8 la fournit.

2.ᵒ En faisant la substitution indiquée, dans un terme de Aₙ de la forme ζa₂ᵅa₄, par exemple ; on obtient ζa₂ᵅ(a₄+a₅x), et il en résulte pour Aₙ₊₁, le terme ζa₂ᵅa₅; ce qui revient à faire varier a₄ de D⁡a₄, et à écrire a₅, à la place de cette dérivée. Si le terme avait été de la forme ζa₂ᵅa₄ᵝ, on aurait obtenu ζa₂ᵅ(a₄ᵝ+D.a₄ᵝ.x+1/2 D².a₄ᵝ.x²+…), et il en serait résulté, pour Aₙ₊₁, le terme ζa₂ᵅ.D.a₄ᵝ=βζa₂ᵅa₄ⁱβ−1a₅; cela revient donc encore à différencier a₄ᵝ d’après les règles ordinaires, et à écrire a₅ à la place de D⁡a₄. C’est ce qui constitue, avec l’observation de la fin du n.ᵒ précédent, la première partie de la règle du n.ᵒ 8.

3.ᵒ Il paraîtrait d’abord que, pour trouver les termes de cette troisième partie, on est obligé de recourir aux termes ou coefficiens antérieurs à celui de Aₙ; mais on peut s’en dispenser, au moyen de l’observation suivante. Si Aₙ₊₁ doit contenir un terme provenant d’une puissance de quantité polynômiale, qui a reçu un accroissement, Aₙ contient aussi un terme dû à cette puissance, qui en est la dérivée immédiatement inférieure ; par exemple, si Aₙ₊₁ doit contenir 1/24 D⁴.a₃⁶=6.5.4.3/1.2.3.4 a₃²a₄⁴, Aₙ contiendra 1/6 D³.a₃⁶= 6.5.4/1.2.3 a₃³a₄³: de plus, ces termes, dus aux puissances des quantités polynômiales, sont toujours aisés à reconnaître, en ce que les deux dernières lettres se suivent, dans l’ordre des indices, et réciproquement ; car, on a évidemment (ζ étant un coefficient numérique convenable) ζa₃³a₄²= ζa₃³(D⁡a₃)²=1/2 D².a₃⁵. Il ne reste donc que de savoir déduire d’un semblable terme dans Aₙ son correspondant dans Aₙ₊₁. Soit donc ζa₂ʳa₃ˢ ce terme dans Aₙ ; on a ζa₂ʳa₃ˢ=ζa₂ʳ(D⁡a₂)ˢ (Dˢ.a₂ʳ⁺ˢ)/1.2…s ; or, le terme correspondant dans Aₙ₊₁, sera (Dˢ⁺¹.a₂ʳ⁺ˢ)/(1.2…(s+1))=1/(s+1) D.(Dˢ.a₂ʳ⁺ˢ)/1.2…s=ζ/(s+1) D.a₂ʳ(D⁡a₂)ˢ =ζr/(s+1).a₂ʳ⁻¹a₃ˢ⁺¹ ce qui revient à différencier l’avant-dernière lettre, ou sa puissance, à écrire a₃ pour D⁡a₂, et à diviser le résultat par l’exposant de la puissance de la dernière lettre, qui se trouve augmenté d’une unité. On ne fait donc autre chose qu’exécuter la seconde partie de la règle du n.ᵒ 8. Cette même partie de la règle, appliquée à la fonction ϕa dans le terme Dⁿ⁡ϕa/1.2…n.a₁ⁿ de Aₙ fournit le terme (Dⁿ⁺¹⁡ϕa)/(1.2…(n+1)).a₁ⁿ⁺¹, dont nous avons parlé au commencement de ce n.ᵒ .

La règle du n.ᵒ 8 est donc parfaitement exacte, et fournit le moyen le plus simple pour déduire le développement d’un terme Aₙ₊₁ de celui du terme Aₙ qui le précède immédiatement.

34. Proposons-nous maintenant de développer immédiatement, et indépendamment des termes qui précèdent, un terme quelconque Aₙ=Dⁿ⁡ϕa/1.2…n de l’équation (6) ou (7).

En faisant

ζ=a₁+a₂x+a₃x²+a₄x³+…,

cette équation devient, d’après le n.ᵒ 32, équation (53)

ϕ(a+ζx)=ϕa+dϕa.ζx+1/2 d²ϕa.ζ²x²+1/6 d³ϕa.ζ³x³+…

Mais, ζ étant lui-même un polynôme, ses puissances sont fonctions de polynômes qui, d’après le n.ᵒ 5, deviennent

(58){ζ²=a₁²+D.a₁².x+ 1/2D².a₁².x²+1/6 D³.a₁².x³+…ζ³=a₁³+D.a₁³.x+1/2 D².a₁³.x²+1/6 D³.a₁³.x³+………………………ζⁿ=a₁ⁿ+D.a₁ⁿ.x+1/2 D².a₁ⁿ.x²+1/6 D³.a₁ⁿ.x³+…

où a₁ doit être considéré comme un premier terme de polynôme ; c’est-à-dire, qu’on a D⁡a₁=a₂, 1/2 D²⁡a₁=a₃, 1/6 D²⁡a₁=a₃, 1/6 D³⁡a₁=a₄, … 1/1.2…(n−1) Dⁿ⁻¹⁡a₁=aₙ.

Substituant ces valeurs dans l’équation. (57), on obtient

ϕ(a+a₁x+a₂x²+a₃x³+…)=

A+A₁x+A₂x²+A₃x³+…+Aₙxⁿ+…

=ϕa+D⁡ϕa.a₁x+D⁡ϕa.D⁡a₁+1/2 D²⁡ϕa.a₁²}x²+D⁡ϕa.1/2 D²⁡a₁+1/2 D²⁡ϕa.D.a₁²+1/6 D³⁡ϕa.D.a₁³}…x³+……+D⁡ϕa.Dⁿ⁻¹⁡a₁/1.2…(n−1) +1/2 D²⁡ϕa.Dⁿ⁻²⁡a₁²/1.2…(n−2) +1/6 D³⁡ϕa.Dⁿ⁻³⁡a₁³/1.2…(n−3) +………+Dⁿ⁻¹⁡ϕa/1.2…(n−1) D.a₁ⁿ⁻¹+Dⁿ⁡ϕa/1.2…n.a₁ⁿ}+…+…xⁿ+…+…+…

d’où l’on tire

Aₙ=Dⁿ⁡ϕa/1.2…n=D⁡ϕa.Dⁿ⁻¹⁡a₁/1.2…(n−1)+1/2 D²⁡ϕa.Dⁿ⁻²⁡a₁²/1.2…(n−2)

+1/6 D³⁡ϕa.Dⁿ⁻³⁡a₁³/1.2…(n−3)+…+Dⁿ⁻¹⁡ϕa/1.2…(n−1) D.a₁ⁿ⁻¹+Dⁿ⁡ϕa/1.2…n.a₁ⁿ;

ou, en écrivant cette formule à rebours,

Aₙ=Dⁿ.ϕa/1.2…n=Dⁿ⁡ϕa/1.2…n.a₁ⁿ+Dⁿ⁻¹⁡ϕa/1.2…(n−1).D.a₁ⁿ⁻¹

+Dⁿ⁻²⁡ϕa/1.2…(n−2).1/2 D².a₁ⁿ⁻²+…

+1/2 D²⁡ϕa.Dⁿ⁻²⁡a₁²/1.2…(n−2)+D⁡ϕa.Dⁿ⁻¹⁡a₁/1.2…(n−1).

Les quantités qui restent à développer, dans cette dernière formule, se succèdent dans l’ordre suivant

a₁ⁿ, D.a₁ⁿ⁻¹, 1/2 D².a₁ⁿ⁻², 1/6 D³.a₁ⁿ⁻³, Dⁿ⁻².a₁²/1.2…(n−2), Dⁿ⁻¹.a₁/1.2…(n−1).

Si tous les exposans de a₁, sous les signes de dérivation, étaient les mêmes et égaux à n, on appliquerait immédiatement, au développement de ces quantités, la règle du n.ᵒ 8 ; mais, comme ils vont toujours en diminuant, il est nécessaire, avant tout, de faire subir à chaque terme une préparation qui consiste à diminuer l’exposant de a₁ d’une unité, d’un terme au suivant, et à modifier en conséquence les coefficiens numériques provenant de ces exposans.

Pour trouver la règle de cette préparation, observons que les dérivées (62) se développent elles-mêmes selon la formule (61), et qu’on a, en général,

Dʳ.a₁ⁿ⁻ʳ/1.2…(n−r)=Dʳ⁡a₁ⁿ⁻ʳ/1.2…r.a₂ʳ+Dʳ⁻¹⁡a₁ⁿ⁻ʳ/1.2…(r−1).D.a₂ʳ⁻¹

+Dʳ⁻²⁡a₁ⁿ⁻ʳ/1.2…(r−2).1/2 D².a₂ʳ⁻²+…

+Dʳ⁻ˢ⁡a₁ⁿ⁻ʳ/1.2…(r−s).Dˢ⁡a₂ʳ⁻ˢ/1.2…s+…+D.a₁ⁿ⁻ʳ.Dʳ⁻¹⁡a₂/1.2…(r−1).

Or, pour déduire de ce développement celui de (Dʳ⁺¹⁡a₁ⁿ⁻ʳ⁻¹)/(1.2…(r+1)), il suffit d’en déduire d’abord celui de Dʳ⁡a₁ⁿ⁻ʳ⁻¹/1.2…r et d’appliquer à ce dernier la règle du n.ᵒ 8. À cet effet, on changera, dans tous les termes de la formule (63), n en n−1 ; mais voyons ce qui en résultera pour un terme quelconque. On a, d’après les règles ordinaires de la différentiation,

Dʳ⁻ˢ⁡a₁ⁿ⁻ʳ/1.2…(r−s)=((n−r)(n−r−1)(n−r−2)…(n−2r+s+2)(n−2r+s+1))/1.2…(r−s).a₁ⁿ⁻²ʳ⁺ˢ,

Dʳ⁻ˢ⁡a₁ⁿ⁻ʳ⁻¹/1.2…(r−s)=((n−r−1)(n−r−2)…(n−2r+s+1)(n−2r+s))/1.2.3…(r−s).a₁ⁿ⁻²ʳ⁺ˢ⁻¹.

Ainsi, pour déduire le développement de Dʳ.a₁ⁿ⁻ʳ⁻¹/1.2…r de celui de Dʳ.a₁ⁿ⁻ʳ/1.2…r, il suffit de diviser chaque terme de ce développement par n−r, de le multiplier par l’exposant de a₁ dans ce terme, et de diminuer cet exposant d’une unité ; ce qui fournit la règle pratique suivante.

35. Pour déduire le développement de (Dʳ⁺¹.a₁ⁿ⁻ʳ⁻¹)/(1.2…(r+1)) de celui de Dʳ.a₁ⁿ⁻ʳ/1.2…r, divisez chaque terme de ce dernier développement par n−r, multipliez-le par l’exposant de a₁ dans ce terme (en observant que, dans les termes sans a₁, cet exposant est zéro), et diminuez son exposant d’une unité. Après cette préparation, suivez la règle du n.ᵒ 8.

Pour donner un exemple de cette règle, nous allons l’appliquer au développement de A₆, dans l’équation (6) du (59). Les quantités à développer, dans ce cas, sont

a₁⁶, D.a₁⁵, 1/2 D².a₁⁴, 1/6D³.a₁³, 1/24 D⁴.a₁², 1/120 D⁵.a₁.

La première de ces quantités reste a₁⁶ la dérivée D.a₁⁵ donne 5a₁⁴a₂; pour en déduire celle 1/2 D².a₁⁴, il faut la diviser par 5, multiplier par 4, exposant de a₁, et diminuer cet exposant d’une unité ; ce qui donne 4a₁³a₂; appliquant ensuite la règle du n.ᵒ 8 à ce terme ainsi préparé, on trouve 4a₁³a₃+4.3/2 a₁²a₂². Pour déduire de cette dérivée celle 1/6 D³.a₁³, il faut diviser le tout par 4, multiplier respectivement les deux termes par 3 et 2, exposans de a₁ et diminuer ces exposans d’une unité ; ce qui donne 3a₁²a₃+3a₁a₂²: appliquant la règle du n.ᵒ 8, on obtient 3a₁²a₄+2.3a₁a₂a₃+a₂³. Pour déduire de cette dernière celle 1/24 D⁴.a₁², il faut diviser le tout par 3, multiplier les trois termes respectivement par 2, 1, 0, exposans de a₁, et diminuer ces exposans d’une unité ; ce qui donne 2a₁a₄+2a₂a₃: appliquant la règle du n.ᵒ 8, on obtient 2a₁a₅+2a₂a₄+a₃². Enfin, pour déduire de cette dérivée celle 1/120 D⁵.a₁, il faut diviser tous les termes par 2, les multiplier respectivement par 1, 0, 0, exposans de a₁, et diminuer ces exposans d’une unité ; ce qui donne a₅, dont la dérivée est a₆, d’après la règle du n.ᵒ 8. En rassemblant tous ces termes, et les multipliant par leurs coefficiens respectifs (61), on aura le développement de A₆, écrit en sens inverse.

L’énoncé de ces opérations peut paraître un peu long ; mais leur exécution est très-expéditive. Après s’être exercé à calculer quatre ou cinq termes, on en a tellement l’habitude qu’il n’en coûte plus, pour ainsi dire, que la peine de les écrire.

36. Remarque. La règle précédente donne non seulement le moyen d’écrire immédiatement le coefficient d’une puissance quelconque de x, dans le développement de ϕ(a+a₁x+a₂x²+…), mais encore une partie quelconque de ce coefficient, sans calculer le reste. Ainsi, si l’on demande le coefficient de Dʳ⁡ϕa/1.2…r, dans le développement de Aₙ=Dⁿ.ϕa/1.2…n, l’équation (61) indiquera que ce coefficient est Dⁿ⁻ʳ.a₁ʳ/1.2…(n−r), dont le développement peut s’exécuter immédiatement, d’après la règle précédente et l’équation (61). Cette observation peut avoir les applications les plus utiles, dans la théorie des hasards, et dans celle de la partition des nombres. Nous avons vu au n.ᵒ 10 que le coefficient de Dʳ.ϕa/1.2…n dans Dⁿ.ϕa/1.2…n était composé de tous les produits de r lettres qu’on peut former avec les quantités polynômiales a₁, a₂, a₃, …, de manière que la somme des indices de chaque produit soit égale à n, et que les coefficiens numériques de chaque produit indiquaient le nombre des permutations dont les lettres de ces produits sont susceptibles : nous aurons donc immédiatement tous ces produits, avec leurs coefficiens numériques, en développant la dérivée Dⁿ⁻ʳ.a₁ʳ/1.2…(n−r). De plus, le nombre des termes dont ce développement sera composé indiquera de combien de manières on peut composer le nombre n, avec r nombres, égaux ou inégaux. Ainsi, en supposant n=12, r=8, on aura, pour le coefficient de D⁴.ϕa/1.2…8 dans D¹².ϕa/1.2…12,

1/24 D⁴.a₁⁸=1/24 D⁴.a₁⁸.a₂⁶+1/6 D³.a₁⁸.D.a₂³+1/2 D².a₁⁸.1/2 D².a₂²

+D.a₁⁸.1/6 D³.a₂

=8.7.6.5/1.2.3.4 a₁⁴a₁⁴+8.7.6/1.2.3 a₁⁵.3a₂²a₃+8.7/1.2 a₁⁶(2a₂a₄+a₃²)+8/1 a₁⁷a₅.

Ce coefficient étant composé de cinq termes, fait voir que le nombre 12 peut être formé de cinq manières différentes, par l’addition de huit nombres, savoir : 1+1+1+1+2+2+2+2, 1+1+1+1+1+2+2+3,

1+1+1+1+1+1+2+4, 1+1+1+1+1+1+3+3, 1+1+1+1+1+1+1+5, lesquels sont donnés immédiatement par les indices et exposans des lettres des produits.

37. La règle du n.ᵒ 15 n’est qu’un corollaire de celle du n.ᵒ 8 et de celle du n.ᵒ 12, qui est une suite évidente des équations (12) : en effet, si dans l’équation (15) on suppose a₂=0, a₃=0, …, b₂=0, b₃=0, …, elle deviendra

ϕ(a+a₁x)×ψ(b+b₁x)=A+A₁x+A₂x²+A₃x³+…

or, le premier membre de cette équation devient, d’après l’équation (53),

(ϕa+D⁡ϕa.a₁x+1/2 D²⁡ϕa.a₁²x²+…)(ψb+D⁡ψb.b₁x+ 1/2 D²⁡ψb.b₁²x²+…)

Ce produit étant développé, d’après l’équation (13), donne

ϕ(a+a₁x)×ψ(b+b₁x)=

A+A₁x+A₂x²+A₃x³

=ϕaψb+ϕa.D⁡ψb.b₁+D⁡ϕa.ψb.a₁}x+ϕa.1/2 D²⁡ψb.b₁²+D⁡ϕa.D⁡ψb.a₁b₁+1/2 D²⁡ϕa.ψb.a₁³}x²+ϕa.1/6 D³⁡ψb.b₁³+D⁡ϕa.1/2 D²⁡ψb.a₁b₁²+1/2 D²⁡ϕa.D⁡ψb.a₁²b₁+1/6 D³⁡ϕa.ψb.a₁³}x³

+…+………………………+Aₙxⁿ+…

+…+ϕa.Dⁿ⁡ψb/1.2…n.b₁ⁿ+…+D⁡ϕa.Dⁿ⁻¹⁡ψb/1.2…(n−1).a₁b₁ⁿ⁻¹+…+1/2 D²⁡ϕa.Dⁿ⁻²⁡ψb/1.2…(n−2).a₁²b₁ⁿ⁻²+…+1/6 D³⁡ϕa.Dⁿ⁻³⁡ψb/1.2…(n−3).a₁³b₁ⁿ⁻³+…+………………+Dⁿ⁻¹⁡ϕa/1.2…(n−1).D⁡ψb.a₁ⁿ⁻¹b₁+Dⁿ⁡ϕa/1.2…n.ψb.a₁ⁿ}+…+…+…xⁿ+…+…+…

Ici les colonnes qui forment les coefficient de x ne sont autre chose que les dernières colonnes des équations (18). Or, le produit (65) s’effectuant comme le produit (13), avec la seule différence qu’à la place de a et de b, il faut écrire ϕa et ψb, et D.ϕa=D⁡ϕa₁, D.ψb=D⁡ψb.b₁, …

à la place de D⁡a, D⁡b, …; ; il s’ensuit que, pour déduire la dernière colonne de Aₙ₊₁ de celle de Aₙ [équations (18)], il faut faire varier ψb dans tous les termes de cette colonne, et de plus faire varier ϕa dans le dernier terme Dⁿ⁡ϕa/1.2.3…n.ψb.a₁ⁿ de cette même colonne. Les deux premières parties du n.ᵒ 15 se trouvent donc démontrées.

Les n autres colonnes qui composent Aₙ₊₁ n+l ne peuvent donc provenir que de la variation des quantités polynômiales a₁, a₂, a₃, …, b₁, b₂, b₃, …

c’est-à-dire, de la substitution de a₁+a₂x pour a₁,

de a₂+a₃x pour a₂,

de b₁+b₂x pour b₁,

de b₂+b₃x pour b₂, …,

dans les termes précédens. Il faut donc appliquer ici la règle du n.ᵒ 8, modifiée par la coexistance de deux polynômes indépendans, c’est-à-dire, par la règle du n.ᵒ 12 ; ce qui constitue la première partie de la règle du n.ᵒ 15. Cette règle se trouve donc entièrement justifiée.

38. Passons maintenant au développement immédiat, et indépendant des termes qui précèdent, d’un terme quelconque de l’équation (15), ou du terme général Aₙ=Dⁿ⁡(ϕa.ψb)/1.2…n.

En effectuant complètement le développement indiqué par la dernière des équations (17), d’après les n.ᵒˢ 34 et 35, et l’ordonnant selon la somme des exposans de dérivation, relatifs à ϕa et ψb, on peut le mettre sous la forme suivante :

Dⁿ.(ϕa.ψb)/1.2…n=Dⁿ⁡ϕa/1.2…n.ψb.a₁ⁿ+Dⁿ⁻¹⁡ϕa/1.2…(n−1).ψb.D.a₁ⁿ⁻¹+…

+Dⁿ⁻¹⁡ϕa/1.2…(n−1).D⁡ψb.a₁ⁿ⁻¹b₁+Dⁿ⁻²⁡ϕa/1.2…(n−2).D⁡ψb.D.(a₁ⁿ⁻²b₁)+…

+Dⁿ⁻¹⁡ϕa/1.2…(n−2).1/2 D²⁡ψb.a₁ⁿ⁻²b₁²+Dⁿ⁻³⁡ϕa/1.2…(n−3). 1/2D²⁡ψb.D⁡(a₁ⁿ⁻³b₁²)+…+………………………………………………+…+1/2 D²⁡ϕa.Dⁿ⁻²⁡ψb/1.2…(n−2).a₁²b₁ⁿ⁻²+D⁡ϕa.Dⁿ⁻²⁡ϕa/1.2…(n−2).D.(a₁b₁ⁿ⁻²)+D⁡ϕa.Dⁿ⁻¹⁡ψb/1.2…(n−1).a₁b₁ⁿ⁻¹+ϕa.Dⁿ⁻¹⁡ψb/1.2…(n−1).D.b₁ⁿ⁻¹+ϕa.Dⁿ⁡ψb/1.2…n.b₁ⁿ…+1/2 D²⁡ϕa.ψb.Dⁿ⁻².a₁²/1.2…(n−2)+D⁡ϕa.ψb.Dⁿ⁻¹⁡a₁/1.2…(n−1) …+D⁡ϕa.D⁡ψb.Dⁿ⁻².(a₁b₁)/1.2…(n−2)+ϕa.D⁡ψb.Dⁿ⁻¹⁡b₁/1.2…(n−1) …+ϕa.1/2 D²⁡ψb.Dⁿ⁻²⁡b₁²/1.2…(n−2)

Il n’y a plus, dans cette formule, dont la loi est très-élégante, que des fonctions de quantités polynômiales à développer ; et elles se succèdent par colonnes dans l’ordre suivant :

(68){a₁ⁿ, D.a₁ⁿ⁻¹, 1/2 D².a₁ⁿ⁻², … Dⁿ⁻².a₁²/1.2…(n−2), … Dⁿ⁻¹.a₁/1.2…(n−1), a₁ⁿ⁻¹b₁, D.(a₁ⁿ⁻²b₁), 1/2 D².(a₁ⁿ⁻³b₁), … Dⁿ⁻²⁡(a₁b₁)/1.2…(n−2), … Dⁿ⁻¹.b₁/1.2…(n−1), a₁ⁿ⁻²b₁², D.(a₁ⁿ⁻³b₁²), 1/2D².(a₁ⁿ⁻⁴b₁²), … Dⁿ⁻².b₁/1.2…(n−2) a₁ⁿ⁻²b₁³, D.(a₁ⁿ⁻⁴b₁³), 1/2 D².(a₁ⁿ⁻⁵b₁³), ………………a₁²b₁ⁿ⁻², D.(a₁b₁ⁿ⁻²), 1/2 D².b₁ⁿ⁻², a₁b₁ⁿ⁻¹, D.b₁ⁿ⁻¹b₁ⁿ,

Si tous les exposans de a₁, sous le signe de dérivation, étaient les mêmes que dans les termes correspondons de la première colonne, on pourrait appliquer immédiatement au développement de ces quantités la règle du n.ᵒ 15 ; mais, comme ces exposans vont en diminuant d’une unité, d’une colonne à l’autre, comme au n.ᵒ 34, il est nécessaire de faire subir à chaque terme la même préparation que dans ce n.ᵒ ; c’est-à-dire, qu’il faut soumettre chaque terme à la règle du n.ᵒ 35, et ensuite y appliquer celle du n.ᵒ 15. Par ce moyen, on peut développer immédiatement un terme quelconque Dⁿ.(ϕaψb)/1.2…n de l’équation (15), indépendamment de ceux qui le précèdent.

39. Remarque. On peut faire ici une observation analogue à celle du n.ᵒ 36. Par le procédé du n.ᵒ précédent, on peut aussi calculer immédiatement un terme quelconque de Dⁿ.(ϕaψb)/1.2…n, indépendamment des autres : ainsi le coefficient de Dʳ⁡ϕa/1.2…r.Dˢ⁡ψb/1.2…s sera Dⁿ⁻ʳ⁻ˢ.(a₁ʳb₁ˢ)/1.2…(n−r−s), dont le développement s’exécutera par le n.ᵒ précédent, en remplaçant n par n−r−s, ϕa par a₁ʳ, et ψb par b₁ˢ, et considérant a₁ et b₁ comme des premiers termes de polynômes. Supposant donc n=12, r=7, s=2, on aura, pour le coefficient de D⁷⁡ϕa/1.2…7.D²⁡ψb/1.2, dans D¹².(ϕa.ψb)/1.2…12,

D³.(a₁⁷b₁²)/1.2.3=1/6 D³⁡a₁⁷b₁²a₂³+ 1/2D²⁡a₁⁷b₁².D.a₂²+D⁡a₁⁷b₁². 1/2D².a₂

+1/2 D²⁡a₁⁷b₁²a₂²b₂+D⁡a₁⁷.D⁡b₁².D.(a₂b₂)+a₁⁷.D⁡b₁². 1/2 D²⁡b₂

+D²⁡a₁⁷.1/2 D²⁡b₁².a₂b₂²+a₁⁷. 1/2D²⁡b₁²D.b₂²

=7.6.5/1.2.3 a₁⁴.b₁².a₂³+7.6/1.2 a₁⁵.b₁².2a₂a₃+7a₁⁶.b₁².a₄

+7.6/1.2 a₁⁵.2b₁.a₂b₂+7/1 a₁⁶.2b₁(a₂b₃+a₃b₂)+a₁⁷.2b₁.b₄

+7/1 a₁⁶.a₂b₂²+a₁⁷.2b₂b₃

D’après la remarque du n.ᵒ 17, qui s’applique également ici, le procédé du n.ᵒ précédent donne aussi le moyen de calculer immédiatement un terme quelconque du développement d’une fonction quelconque de deux polynômes indépendans ϕ(a+a₁x+a₂x²+…, b+b₁x+b₂x²+…); il suffit pour cela de remplacer, dans la formule (67) les produits des dérivées de ϕa et ψb par les dérivées partielles correspondantes de ϕ(a, b) ; et le n.ᵒ précédent fait voir avec quelle facilité le calcul des dérivations fournit la solution de ce problème compliqué, et intraitable par les méthodes ordinaires.

40. La règle du n.ᵒ 28 est un corollaire bien simple de celles des n.ᵒˢ 15 et 35 ; en effet, la forme du terme général (48), nBₙ=Dⁿ⁻¹⁡(c⁻ⁿ.D.ϕb)/1.2…(n−1), étant comparée à celle (17) Aₙ₋₁n=Dⁿ⁻¹.(ϕaψb)/1.2…(n−1), fait voir qu’on obtient la première, en remplaçant, dans celle-ci ψb par c⁻ⁿ, et ϕa par D.ϕb. Les règles de développement doivent donc être les mêmes pour l’une et l’autre formes, sauf les différences suivantes : 1.ᵒ l’exposant de c diminuant d’une unité d’un terme à l’autre, il faut faire subir à chaque terme, avant d’en déduire le suivant, la préparation du n.ᵒ 35 ; 2.ᵒ  ϕa étant remplacé par D.ϕb, il, ensuit qu’on a, en général, Dʳ.D.ϕb/1.2…r=(Dʳ⁺¹.ϕb)/1.2…r=(r+1) (Dʳ⁺¹.ϕb)/(1.2…(r+1)), et Dʳ.D.ϕb/1.2…r =(r+1) (Dʳ⁺¹.ϕb)/(1.2…(r+1)) ; ce qui produit les coefficiens numériques égaux aux exposans de dérivation ; 3.ᵒ enfin, nous avons ordonné différemment les termes des équations (50), en transformant les lignes horizontales des équations (18) en colonnes, et réciproquement : il en est résulté que les dernières colonnes des équations (18) sont devenues les derniers termes de chaque colonne des équations (50), en transformant les lignes horizontales des équations (18) en colonnes, et réciproquement : il en est résulté que les dernières colonnes des équations (18) sont devenues les derniers termes de chaque colonne des équations (50) ; ce qui a produit les modifications des 2.ᵐᵉ et 3.ᵐᵉ parties de la règle du n.ᵒ 15.

41. En tenant compte des observations du n.ᵒ précédent, la formule (67) fournit le moyen de développer immédiatement un quelconque des termes (48), indépendamment des précédens : on a, en général,

nBₙ=Dⁿ⁻¹⁡(c⁻ⁿ.D.ϕb)/1.2…(n−1)=

n Dⁿ⁡ϕb/1.2…n.(c⁻ⁿ).b₁ⁿ+(n−1) Dⁿ⁻¹⁡ϕb/1.2…(n−1).(c⁻ⁿ).D.b₁ⁿ⁻¹+…

+(n−1) Dⁿ⁻¹⁡ϕb/1.2…(n−1).D.(c⁻ⁿ).b₁ⁿ⁻¹c₁

+(n−2) Dⁿ⁻²⁡ϕb/1.2…(n−2).D.(c⁻ⁿ).D.(b₁ⁿ⁻²c₁)+…

+(n−2) Dⁿ⁻²⁡ϕb/1.2…(n−2).1/2 D².(c⁻ⁿ).b₁ⁿ⁻²c₁²

+(n−3) Dⁿ⁻³⁡ϕb/1.2…(n−3).1/2 D².(c⁻ⁿ).D.(b₁ⁿ⁻³c₁²)+…

+3.1/6 D³⁡ϕb.Dⁿ⁻³⁡(c⁻ⁿ)/1.2…(n−3).b₁³c₁ⁿ⁻³

+2.1/2 D²⁡ϕb.Dⁿ⁻³⁡(c⁻ⁿ)/1.2…(n−3).D.(b₁²c₁ⁿ⁻³)+…

+2.1/2 D²⁡ϕb.Dⁿ⁻²⁡(c⁻ⁿ)/1.2…(n−2).b₁²c₁ⁿ⁻²

+D⁡ϕb.Dⁿ⁻²⁡(c⁻ⁿ)/1.2…(n−2).D.(b₁c₁ⁿ⁻²)

+D⁡ϕb.Dⁿ⁻¹⁡(c⁻ⁿ)/1.2…(n−1).b₁c₁ⁿ⁻¹

…+3.1/6 D³⁡ϕb.(c⁻ⁿ).Dⁿ⁻³⁡b₁³/1.2…(n−3)+2.1/2 D²⁡ϕb.(c⁻ⁿ).Dⁿ⁻².b₁²/1.2…(n−2)

+D⁡ϕb.(c⁻ⁿ).Dⁿ⁻¹.b₁³/1.2…(n−1)

…+2.1/2 D²⁡ϕb.D⁡(c⁻ⁿ).Dⁿ⁻³⁡(b₁²c₁)/1.2…(n−3)

+D⁡ϕb.D⁡(c⁻ⁿ).Dⁿ⁻²⁡(b₁c₁)/1.2…(n−2)

…+D⁡ϕb.1/2 D².(c⁻ⁿ).Dⁿ⁻³⁡(b₁c₁²)/1.2…(n−3)

où les quantités qui restent à développer sont de la même forme que celles (68), du n.ᵒ 38, et doivent être développées de la même manière.

42. Au moyen du n.ᵒ précédent, nous sommes donc en état de calculer immédiatement un terme quelconque d’une fonction de polynôme, ordonné selon les puissances d’une fonction ou d’un polynôme donné ; ce qui constitue le problème général du retour des fonctions et des séries, étendu aux fonctions de polynômes. De plus, d’après la remarque du n.ᵒ 39, qui est applicable à ce cas, nous pouvons aussi calculer immédiatement une partie quelconque d’un terme, sans calculer le reste de ce terme. Mais, ce qu’il y a de plus remarquable, c’est que cette question difficile est résolue d’une manière si simple qu’on n’a, pour ainsi dire, que la peine d’écrire le résultat.

43. Résumons, en deux mots, l’objet et l’esprit du calcul des dérivations, tel qu’il résulte de ce petit écrit. Le théorème de Taylor donne le développement d’une fonction simple d’un binôme, selon les puissances ascendantes de la variable principale, ou selon les mêmes puissances d’une fonction quelconque donnée de cette variable. Le passage du théorème de Taylor au développement des fonctions de polynômes, ou des fonctions de fonctions, selon les puissances ascendantes de la variable, n’est autre chose que le passage de la différentiation d’une fonction, en regardant la différentielle de la variable principale comme constante, à la différentielle de la même fonction, en ne regardant aucune différentielle comme constante. Quant au passage du développement d’une fonction, selon les puissances ascendantes de la variable à celui selon les puissances ascendantes d’une fonction donnée de cette variable ; (ce qui constitue le retour des fonctions et des séries) ; il n’est autre chose que celui de la différentiation d’une fonction, en changeant de variable principale ou indépendante.

44. Me voici parvenu au terme que je m’étais proposé : celui de déduire la véritable théorie du calcul des dérivations du seul théorème de Taylor, sans l’emploi d’aucun principe nouveau. J’espère que les géomètres verront avec plaisir ce beau corollaire d’un théorème qui a déjà été si fécond. Le cadre étroit dans lequel j’ai resserré l’essence de ce calcul les engagera sans doute à donner quelques momens à la lecture de ce petit écrit ; et j’ose présumer qu’elle les réconciliera avec le calcul des dérivations, dont l’ouvrage d’Arbogast a pu les éloigner. Mon but n’a pas été d’épuiser la matière, mais d’éveiller l’attention des géomètres sur l’utilité, trop méconnue, des dérivations ; et de leur éviter la recherche pénible de nouveaux moyens de développement, en leur présentant ceux qui sont, à la fois, les plus simples et les plus expéditifs qu’on puisse trouver.

Les géomètres auxquels l’Analise combinatoire est familière verront, par nos remarques des n.ᵒˢ 10, 16 et 36, que le calcul des dérivations contient, non seulement les véritables sources des règles de cette analise, et leur extension à des fonctions de plusieurs polynômes indépendans, mais encore les moyens d’exécution les plus commodes et les plus rapides.

