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

## ANALISE TRANSCENDANTE.

### Essai sur un nouveau mode d’exposition des principes du calcul différentiel, du calcul aux différences et de l’interpolation des suites, considérées comme dérivant d’une source commune ;

I. Le calcul différentiel, le calcul aux différences et les diverses méthodes d’interpolation reposent également sur un petit nombre de formules générales, applicables à toutes les fonctions, et qu’on n’a démontrées jusqu’ici que par des considérations souvent compliquées, presque toujours déduites de principes éloignée, ou d’inductions de nature à laisser des doutes sur leur généralité. On remarque surtout, dans l’exposition de ces diverses branches d’analise, une variété de procédé et de raisonnemens qui ne laisse que difficilement apercevoir leur liaison réciproque, et l’identité des principes dont elles, ne sont pourtant, en quelque sorte, que des traductions variées.

Frappés de ces considérations, nous avons pensé faire une chose qui pourrait intéresser les géomètres, en déduisant toutes ces diverses formules de quelques théorèmes nouveaux ou peu connus, que nous démontrerons d’abord, et dont il nous suffira ensuite de traduire les énoncés dans l’algorithme reçu, relatif aux divers cas particuliers, pour en voir éclore, d’une manière tout-à-fait naturelle, les formules connues qui répondent à chacun d’eux.

2. Soit f⁡x une fonction de ce de forme quelconque, de manière que f⁡a, f⁡b, f⁡c, … soient respectivement ce que devient cette fonction, lorsqu’on y fait, tour-à-tour, x=a, x=b, x=c, …. Soient posés successivement

f⁡b−f⁡a/b−a=f₁⁡(a, b), f⁡c−f⁡a/c−a=f₁⁡(a, c), f⁡d−f⁡a/d−a=f₁⁡(a, d), …

f₁⁡(a, c)−f₁⁡(a, b)/c−b=f₂⁡(a, b, c), f₁⁡(a, d)−f₁⁡(a, b)/d−b=f₂⁡(a, b, d), …

f₂⁡(a, b, d)−f₂⁡(a, b, c)/d−c=f₃⁡(a, b, c, d), …

Les fonctions f₁, f₂, f₃, … sont ce que nous appellerons à l’avenir les fonctions interpolaires des différens ordres des quantités a, b, c, d, …

3. La première remarque que nous ferons au sujet de ces sortes de fonctions, c’est qu’elles sont toujours symétriques, de telle sorte qu’on y peut intervertir, comme on voudra, l’ordre des élémens a, b, c, d, … qui concourent à leur formation, sans qu’elles eu éprouvent aucun changement. En effet, on a, en premier lieu,

f₁⁡(a, b)=f⁡b−f⁡a/b−a=f⁡a/a−b+f⁡b/b−a,

fonction évidemment symétrique en a et b, et on pourrait en dire autant de toutes les fonctions interpolaires à deux lettres ou du premier ordre.

On aura donc, semblablement,

f₁⁡(a, c)=f⁡a/a−c+f⁡c/c−a ;

mais,

f₂⁡(a, b, c)=f₁⁡(a, c)−f₁⁡(a, b)/c−b=f₁⁡(a, b)/b−c+f₁⁡(a, c)/c−b ;

mettant donc pour f₁⁡(a, b) et f₁⁡(a, c) les valeurs ci-dessus, il viendra

f₂⁡(a, b, c)=f⁡a/(a−b)(b−c)+f⁡b/(b−a)(b−c)+f⁡a/(a−c)(c−b)+f⁡c/(c−a)(c−b) ;

ou, en réduisant à une seule les deux fractions de même numérateur,

f₂⁡(a, b, c)=f⁡a/(a−b)(a−c)+f⁡b/(b−a)(b−c)+f⁡c/(c−a)(c−b) ;

fonction également symétrique, et on démontrerait la même chose de toutes les fonctions interpolaires à trois lettres, ou du second ordre.

On pourrait, par des transformations analogues, étendre progressivement le même principe aux fonctions interpolaires des ordres supérieurs ; mais, afin de ne pas nous borner à une simple induction, à l’égard d’un principe qui est fondamental, dans la théorie qui nous occupe, prouvons que, si la symétrique se soutient jusqu’aux fonctions du (n−1). ⁱᵉᵐᵉ ordre, inclusivement, elle aura lieu également pour celles du n. ⁱᵉᵐᵉ.

Supposons donc que, pour les n lettres a, b, c, …p, q, r, on ait trouvé

fₙ₋₁⁡(a, b, c, …p, q, r)={ f⁡a/(a−b)(a−c)…(a−q)(a−r)+f⁡b/(b−a)(a−c)…(b−q)(b−r) ………………+f⁡q/(q−a)(q−b)…(q−p)(q−r) +f⁡r/(r−a)(r−b)…(r−p)(r−q) };

nous aurons semblablement, pour les n lettres a, b, c, …p, q, s, …

fₙ₋₁⁡(a, b, c, …p, q, s)={ f⁡a/(a−b)(a−c)…(a−q)(a−s)+f⁡b/(b−a)(b−c)…(b−q)(b−s) ………………+f⁡q/(q−a)(q−b)…(q−p)(q−s) +f⁡s/(s−a)(s−b)…(s−p)(s−q) };

fonctions qui sont symétriques l’une et l’autre. Mais, en vertu de la définition des fonctions interpolaires, on a, pour les n+1 lettres a, b, c, …p, q, s,

fₙ⁡(a, b, c, …p, r, s)=fₙ₋₁⁡(a, b, c, …p, q, s)−fₙ₋₁⁡(a, b, c, …p, q, r)/s−r, =fₙ₋₁⁡(a, b, c, …p, q, r)/r−s+fₙ₋₁⁡(a, b, c, …p, q, s)/s−r ;

il viendra donc, en substituant,

fₙ⁡(a, b, c, …q, r, s)

={ f⁡a/(a−b)(a−c)…(a−q)(a−r)+f⁡a/(a−b)(a−c)…(a−q)(a−s) + f⁡b/(b−a)(b−c)…(b−q)(b−r)+f⁡b/(b−a)(b−c)…(b−q)(b−s) +……………………+……………………+ f⁡q/(q−a)(q−b)…(q−p)(q−r)+f⁡q/(q−a)(q−b)…(q−p)(q−s) + f⁡r/(r−a)(r−b)…(r−p)(r−q)+f⁡s/(s−a)(s−b)…(s−p)(s−q) }:

où, en réduisant à une seule les fractions de même numérateur,

fₙ⁡(a, b, c, …p, q, s)={ f⁡a/(a−b)(a−c)…(a−r)(a−s)+f⁡b/(b−a)(b−c)…(b−r)(b−s) +………………+f⁡r/(r−a)(r−b)…(r−q)(r−s) +f⁡r/(s−a)(s−b)…(s−q)(s−r) };

fonction également symétrique. Il demeure donc établi par là que, si la symétrie se soutient jusqu’à l’emploi de la n.ⁱᵉᵐᵉ lettre, inclusivement, elle aura lieu encore après l’introduction de la (n+1). ⁱᵉᵐᵉ, quel que soit n; puis donc que cette symétrique a lieu en effet pour des fonctions interpolaires formées de deux et de trois lettres, il s’ensuit qu’elle doit être regardée comme un fait analitique généralement démontré.

4. D’après le mode de génération des fonctions interpolaires, il existe une relation fort simple entre deux d’entre elles de même ordre, ne différant l’une de l’autre que par un seul des élémens qui les composent. Puisqu’en effet on a

fₙ₊₁⁡(a, b, c, d, …)=fₙ⁡(b, c, d, …)−fₙ⁡(a, c, d, …)/b−a,

on en conclura

fₙ⁡(a, c, d, …)−fₙ⁡(b, c, d, …)=(a−b)fₙ₊₁⁡(a, b, c, d, …).

On peut également obtenir une relation très-remarquable entre trois de ces fonctions interpolaires de même ordre, ne différant que par l’exclusion donnée tour-à-tour à une lettre, sur trois d’entre elles. On a en effet, par la formule (1)

fₙ⁡(a, c, d, …)−fₙ⁡(b, c, d, …)=(a−b)fₙ₊₁⁡(a, b, c, d, …), fₙ⁡(a, b, d, …)−fₙ⁡(a, c, d, …)=(b−c)fₙ₊₁⁡(a, b, c, d, …).;

d’où, en retranchant du produit de la première par b−c le produit de la seconde par a−b et réduisant,

(a−b)fₙ⁡(a, b, d, …)+(b−c)fₙ⁡(b, c, d, …)+(c−a)fₙ⁡(a, c, d, …)=0;

formule que l’on retiendra facilement dans sa mémoire, en remarquant 1.ᵒ que la somme des trois coefficiens binômes est nulle ; 2.ᵒ que la lettre qui manque dans chacun de ces coefficiens est aussi celle qui manque dans la fonction qu’il multiplie. Remarquons, au surplus, que, bien que nous ayons supposé qu’il s’agissait des trois premières lettres, le théorème n’en demeure pas moins établi généralement pour trois lettres quelconques ; puisqu’à cause de la symétrie des fonctions interpolaires, on peut toujours amener ces trois lettres à être les trois premières.

5. Les fonctions interpolaires des différens ordres sont, comme l’on voit, complètement déterminées, tant que les élémens a, b, c, d, … dont elles se composent sent tous différens les uns des autres ; mais, si tous ou partie d’entre eux sont égaux, il arrive que toutes ou partie des fonctions f₁, f₂, f₃, qu’il faut successivement calculer, pour parvenir à la valeur de fₙ⁡(a, b, c, d, …) se présentent sous la forme 0/0. Or, on sait que, lorsqu’une fonction quelconque, qui n’est susceptible que d’une valeur unique, se présente sous cette forme, en vertu de l’égalité de deux des élémens p et q dont elle se compose, trois cas seulement peuvent se présenter. 1.ᵒ il peut arriver que la différence p−q décroissant indéfiniment, la fonction décroisse aussi indéfiniment, de manière à pouvoir devenir moindre que toute grandeur donnée ; ce qu’on exprime en disant qu’elle devient nulle quand p=q; 2.ᵒ il peut arriver, au contraire, que la différence p−q décroissant indéfiniment, la fonction croisse indéfiniment, de manière à pouvoir surpasser toute grandeur donnée ; ce qu’on exprime en disant qu’elle devient infinie, quand p=q; 3.ᵒ enfin il peut se faire que p−q décroissant indéfiniment, la fonction ne décroisse ni ne croisse indéfiniment, mais seulement de manière à tendre sans cesse vers une grandeur constante, vers une limite finie, dont elle pourra différer de moins de toute grandeur donnée, en prenant p−q d’une petitesse suffisante ; et alors cette limite fixe sera dite la valeur de la fonction qui répond p=q.

Or, nous allons voir qu’à l’égard des fonctions interpolaires, ce troisième cas est, généralement parlant, le seul qui puisse avoir lieu, ou, en d’autres termes, qu’en général ces sortes de fonctions ne sauraient jamais devenir nulles, par l’effet de l’égalité de tous ou partie des élémens dont elles se composent.

Pour le prouver, imaginons qu’on partage l’intervalle b−a en un nombre arbitraire m de parties égales ; posons, pour abréger, b−a=k, et soient

a₁=a+k/m, a₂=a₁+k/m, a₃=a₂+k/m, …b=aₘ₋₁+k/m ;

nous aurons (1)

fₙ₋₁⁡(a₁, c, d, …)−fₙ₋₁⁡(a, c, d, …)=k/m fₙ⁡(a₁, a, c, d, …), fₙ₋₁⁡(a₂, c, d, …)−fₙ₋₁⁡(a₁, c, d, …)=k/m fₙ⁡(a₁, a₂, c, d, …), fₙ₋₁⁡(a₃, c, d, …)−fₙ₋₁⁡(a₂, c, d, …)=k/m fₙ⁡(a₂, a₃, c, d, …), ………………………………………………fₙ₋₁⁡(b, c, d, …)−fₙ₋₁⁡(aₘ₋₁, c, d, …)=k/m fₙ⁡(aₘ₋₁, b, c, d, …);

en ajoutant, réduisant et se rappelant que (1)

fₙ₋₁⁡(b, c, d, …)−fₙ₋₁⁡(a, c, d, …)=kfₙ⁡(a, b, c, d, …),

il viendra, en divisant par k,

fₙ⁡(a, b, c, d, …)=1/m {fₙ⁡(a, a₁, c, d, …)+fₙ⁡(a₁, a₂, c, d, …)+fₙ⁡(a₂, a₃, c, d, …)+……………+fₙ⁡(aₘ₋₁, b, c, d, …)};(3)

c’est-à-dire que fₙ⁡(a, b, c, d, …) est la moyenne arithmétique entre toutes les autres fonctions interpolaires du même ordre déduites de celle-là, en y mettant successivement a et a₁, a₁ et a₂, a₂ et a₃, …, aₘ₋₁ et b, en place de a et b.

Il suit de là qu’excepté le cas où toutes les fonctions dont la somme divisée par m compose le second membre de l’équation (3) seraient égales entre elles, auquel cas chacune d’elles serait égale à celle qui forme le premier membre, il y en aura toujours de plus grandes et de plus petites que celle-là. Or, on peut toujours prendre le nombre m assez grand pour rendre et conséquemment les différences consécutives a₁−a, a₂−a₁, a₃−a₂, …b−aₘ₋₁ d’une petitesse indéfinie ; donc, si une fonction interpolaire telle que fₙ⁡(a, b, c, d, …) pouvait devenir nulle ou infinie, lorsqu’on y fait a=b, on pourrait toujours prendre pour m un assez grand nombre pour rendre toutes les fonctions qui composent le second membre de l’équation (3) moindres ou plus grandes qu’une grandeur donnée quelconque, et, en particulier, moindres ou plus grandes que celle qui forme son premier membre, ce qui rendrait cette équation absurde.

On voit donc que la fonction fₙ⁡(a, b, c, d, …) continue à avoir une valeur déterminée, qui n’est ni nulle ni infinie, lorsqu’on suppose les deux premiers élémens a et b égaux entre eux ; et il en sera de même encore, dans le cas de l’égalité entre deux autres quelconques de ses élémens ; puisqu’en vertu de la symétrie de la fonction, on pourra toujours amener ces deux élémens à être les premiers. Il en serait encore évidemment de même dans le cas de l’égalité de plus de deux élémens entre eux, ou de l’égalité entre plusieurs groupes d’élémens ; et, bien qu’alors le procédé que nous avons indiqué pour parvenir à ces sortes de fonctions n’en assigne plus la valeur, on pourra néanmoins les représenter et les employer dans les calculs, comme si cette valeur était connue. Observons d’ailleurs que rien n’empêche que les fonctions interpolaires, comme toutes les autres fonctions, ne deviennent nulles ou infinies, pour certaines valeurs particulières des élémens qui les composent ; mais cela n’arrivera jamais, tant que ces élémens conserveront leur indétermination.

6. De la formule (1) on peut conclure

f⁡x=f⁡a+(x−a)f₁⁡(a, x), f₁⁡(a, x)=f₁⁡(a, b)+(x−b)f₂⁡(a, b, x), f₂⁡(a, b, x)=f₂⁡(a, b, c)+(x−c)f₃⁡(a, b, c, x), ……………………………………fₙ₋₁⁡(a, b, c, …p, x)=fₙ₋₁⁡(a, b, c, …p, q)+(x−q)fₙ⁡(a, b, c, …p, q, x)

ce qui donne, par des substitutions successives

f⁡x={f⁡a+(x−a)f₁⁡(a, b)+(x−a)(x−b)f₂⁡(a, b, c)+(x−a)(x−b)(x−c)f₃⁡(a, b, c, d)+…………………………………+(x−a)(x−b)(x−c)…(x−p)fₙ₋₁⁡(a, b, c, …p, q)+(x−a)(x−b)(x−c)…(x−p)(x−q)fₙ⁡(a, b, c, …p, q, r).

Cette valeur de f⁡x contient une fonction interpolaire de x, à son dernier terme ; mais, si les élémens a, b, c, …p, q, ont des valeurs comprises entre des limites peu étendues, et que x soit lui-même compris entre eux, la série sera assez convergente pour qu’on puisse se permettre d’en négliger le dernier terme, et on aura alors sensiblement

f⁡x={f⁡a+a(x−a)f₁⁡(a, b)+a(x−a)(x−b)f₂⁡(a, b, c)+a(x−a)(x−b)(x−b)f₃⁡(a, b, c, d)+a(x−a)(x−b)(x−c)f₃⁡(a, b, c, d)+…………………………+a(x−a)(x−b)(x−c)…(x−p)fₙ₋₁⁡(a, b, c, …p, q)}.

On reconnaît ici la formule ordinaire d’interpolation ; ce qui justifie la dénomination d’interpolaires que nous avons donnée aux fonctions dont elle se compose.

7. Posons y=f⁡x, et désignons

Δy

y

x

x+Δx

Δ²y

Δy

x

x+Δx

Δ³y

Δ²y

x

x+Δx

Δⁿy

Δⁿ⁻¹y

x

x+Δx

nous aurons, par la définition des fonctions interpolaire,

y =f⁡x,

Δy =f⁡(x+Δx)−f⁡x=Δx.f₁⁡(x, x+Δx),

Δ²y=Δx{f₁⁡(x+Δx, x+2Δx)−f₁⁡(x, x+Δx)}.

L’analogie conduit à soupçonner qu’en général

Δⁿy=n!Δxⁿ{fₙ⁡(x, x+Δx, x+2Δx, …x+nΔx).

Pour changer ce soupçon en certitude, il suffit de prouver qu’il devra en être ainsi ; si l’on a

Δⁿ⁻¹y=(n−1)!Δxⁿ⁻¹fₙ₋₁⁡(x, x+Δx, x+2Δx, …x+(n−1)Δx);

or, en adoptant cette hypothèse, on aura

Δⁿy=(n−1)!Δxⁿ⁻¹{fₙ₋₁⁡(x, x+Δx, x+2Δx, …x+nΔx)

−fₙ₋₁⁡(x, x+Δx, …x+(n−1)Δx)};

c’est-à-dire,

Δⁿy=n!Δxⁿ{fₙ⁡(x, x+Δx, x+2Δx, …x+nΔx);

comme nous l’avions annoncé.

8. Reprenons la formule (4)

f⁡x=f⁡a+(x−a)f₁⁡(a, b)+(x−a)(x−b)f₂⁡(a, b, c)

+(x−a)(x−b)(x−c)f₃⁡(a, b, c, d)+…

supposons les élémens a, b, c, … équidifférens, et soient posés

b−a=c−b=d−c=…=k;

soient faits, de plus

yₚ=f⁡a, yₚ₊₁=f⁡b, yₚ₊₂=f⁡c, …, yₚ₊ₙ=f⁡x;

nous aurons

Δyₚ/k=f⁡b−f⁡a/b−a=f₁⁡(a, b), Δ²yₚ/k²=2 f₁⁡(b, c)−f₁⁡(a, b)/c−a

=2f₂⁡(a, b, c), Δ³yₚ/k³=…

(Δyₚ₊₁)/k=f⁡c−f⁡b/c−b=f₁⁡(b, c), (Δ²yₚ₊₁)/k²=2 f₁⁡(c, d)−f₁⁡(b, c)/d−b

=2f₂⁡(b, c, d), …

(Δyₚ₊₃)/k=f⁡d−f⁡c/d−c=f₁⁡(c, d), …………

d’où, en substituant dans la formule (4),

yₚ₊ₙ=yₚ+x−a/1.Δyₚ/k+x−a/1.x−b/2.Δ²yₚ/k²

+x−a/1.x−b/2.x−c/3.Δ³yₚ/k³+…

Si nous posons x−a=nk, nous aurons

x−b=(x−a)−(b−a)=nk−k=(n−1)k, x−c=(x−a)−(c−a)=nk−2k=(n−2)k, …………………,

de sorte qu’alors cette formule deviendra

yₚ₊ₙ=yₚ+n/1.Δyₚ+n/1.n−1/2.Δ²yₚ+n/1.n−1/2.n−2/3.Δ³yₚ+…

et, par suite, en supposant p nul,

yₙ=y+n/1.Δy+n/1.n−1/2.Δ²y+n/1.n−1/2.n−2/3.Δ³y+…

ce sont les formules connues d’interpolation, pour le cas de l’équidifférence des valeurs de la variable indépendante.

9. D’après ce qui a été démontré ci-dessus (5), on a

fₙ⁡(a, b, c, …r, s)={ f⁡a/(a−b)(a−c)…(a−r)(a−s)+f⁡b/(b−a)(b−c)…(b−r)(b−s) +f⁡c/(c−a)(c−b)…(c−r)(c−s) +……………………+f⁡s/(s−a)(s−b)…(s−q)(s−r) }

supposant alors, comme ci-dessus,

b−a=c−b=d−c=…=s−r=k,

et posant en outre

yₚ=f⁡a, yₚ₊₁=f⁡b, yₚ₊₂=f⁡c, …yₚ₊ₙ₋₁=f⁡r, yₚ₊ₙ=f⁡s,

on trouvera, en substituant et renversant l’ordre des termes du second membre,

fₙ⁡(a, a+k, a+2k, a+3k, …a+nk)

=(yₚ₊ₙ)/n!kⁿ−(yₚ₊ₙ₋₁)/(n−1)!1!kⁿ+(yₚ₊ₙ₋₂)/(n−2)!2!kⁿ−(yₚ₊ₙ₋₃)/(n−3)!3!kⁿ+

…± (yₚ₊₁)/1!(n−1)!kⁿ ∓ yₚ/n!kⁿ ;

mais si, dans la formule (5), on change y en yₚ, x en a et qu’on y fasse, en outre, Δa=k, elle donnera

fₙ⁡(a, a+k, a+2k, a+3k, …a+nk)=Δⁿyₚ/n!kⁿ ;

donc, en égalant ces deux valeurs et multipliant par n!,

Δⁿyₚ=yₚ₊ₙ−n/1 yₚ₊ₙ₋₁+n/1.n−1/2 yₚ₊ₙ₋₂

−n/1.n−1/2.n−2/3 yₚ₊ₙ₋₃+…;

d’où, en supposant p nul

Δⁿy=yₙ−n/1 yₙ₋₁+n/1.n−1/2 yₙ₋₂

−n/1.n−1/2.n−2/3 yₙ₋₃+…;

formules inverses des formules (7) et (8).

10. Revenons au cas ou plusieurs des élémens de la fonction fₙ sont égaux entre eux. Nous avons observé (5) que, bien qu’alors cette fonction se présentât sous la forme 0/0, elle n’en avait pas moins une valeur déterminée qui, en général, n’était ni nulle ni infinie, et que, bien qu’alors notre procédé ne fut plus propre à en assigner la valeur, les mêmes notations n’en étaient pas moins propres à la représenter ; en continuant donc à les employer, nous aurons d’abord (1)

f₁⁡(b, b)−f₁⁡(a, b)=(b−a)f₂⁡(a, b, b),

f₁⁡(a, b)−f₁⁡(a, a)=(b−a)f₂⁡(a, a, b),

d’où, en ajoutant et réduisant,

f₁⁡(b, b)−f₁⁡(a, a)=(b−a){f₂⁡(a, b, b)+f₂⁡(a, a, b)}.

Nous aurons aussi

f₂⁡(b, b, b)−f₂⁡(a, b, b)=(b−a)f₃⁡(a, b, b, b), f₂⁡(a, b, b)−f₂⁡(a, a, b)=(b−a)f₃⁡(a, a, b, b), f₂⁡(a, a, b)−f₂⁡(a, a, a)=(b−a)f₃⁡(a, a, a, b),

d’où en ajoutant encore et réduisant

f₂⁡(b, b, b)−f₂⁡(a, a, a)=(b−a){f₃⁡(a, b, b, b)+f₃⁡(a, a, b, b)+f₃⁡(a, a, a, b)}.

En continuant de la même manière, on aura en général

fₙ₋₁⁡(b, b, b, …b)−fₙ₋₁⁡(a, a, a, …a)=(b−a){fₙ⁡(a, b, b, …b, b, b)+fₙ⁡(a, a, b, …b, b, b)+fₙ⁡(a, a, a, …b, b, b)+……………+fₙ⁡(a, a, a, …a, b, b)+fₙ⁡(a, a, a, …a, a, b)};

si, dans cette formule, on fait b=a+k, d’où b−a=k, et qu’on change ensuite a en x, elle deviendra

fₙ₋₁⁡(x+k, x+k, x+k, …x+k, x+k)−fₙ₋₁⁡(x, x, x, …x, x)

=k{fₙ⁡(x, x+k, x+k, …x+k, x+k, x+k)+fₙ⁡(x, x, x+k, …x+k, x+k, x+k)+fₙ⁡(x, x, x, …x+k, x+k, x+k)+…………………+fₙ⁡(x, x, x, …x, x+k, x+k)+fₙ⁡(x, x, x, …x, x+k)}

formule qui va nous servir tout-à-l’heure.

11. En conservant toujours les mêmes notations, on a, par la formule fondamentale,

f⁡(x+k)=f⁡x+kf₁⁡(x, x+k), f₁⁡(x, x+k)=f₁⁡(x, x)+kf₂⁡(x, x, x+k), f₂⁡(x, x, x+k)=f₂⁡(x, x, x)+kf₃⁡(x, x, x, x+k), ………………………………………fₙ₋₁⁡(x, x, …x, x+k)=fₙ₋₁⁡(x, x, x, …x, x)+kfₙ⁡(x, x, x, …x, x, x+k);

prenant la somme des produits respectifs de ces équations par 1, k, k², k³, …kⁿ⁻¹, et réduisant, on trouvera

f⁡(x+k)={f⁡x+kf₁⁡(x, x)+k²f₂⁡(x, x, x)+k³f₃⁡(x, x, x, x)+……………………+kⁿ⁻¹fₙ₋₁⁡(x, x, x, …x, x)+kⁿfₙ⁡(x, x, x, …x, x, x+k)};

développement de f⁡(x+k) en série où on pourra toujours prendre k assez petit, sans être nul, pour que la série soit convergente, puisque les fonctions qui multiplient les diverses puissances de k ont toujours des valeurs finies. En outre, si l’on parvient à assigner deux limites entre lesquelles se trouve compris le dernier terme du développement, le seul qui renferme à sous le signe de fonction, on connaîtra aussi par là les limites de l’erreur commise en bornant la série aux seuls termes qui précèdent celui-là.

Examinons présentement, d’une manière plus particulière, la nature des multiplicateurs des diverses puissances de k, dans la formule (12). Soit, en général, désigné par f⁽ᵖ⁾⁡x, la limite vers laquelle tend sans cesse le rapport

(f⁽ᵖ⁻¹⁾⁡(x+k)−f⁽ᵖ⁻¹⁾⁡x)/k,

lorsque k tend vers zéro ; où, en d’autres termes, ce que devient ce rapport, lorsque k devient infiniment petit ou nul. Alors f′⁡x, f″⁡x, f‴⁡x, …, seront ce qu’on appelle les dérivées successives de la fonction f⁡x. On aura d’abord

f₁⁡(x, x)=f′⁡x;

car on a

f₁⁡(x, x+k)=(f⁡(x+k)−f⁡x)/k ;

et les deux quantités f₁⁡(x, x) et f′⁡x sont respectivement ce que deviennent les deux membres de cette équation, lorsqu’on suppose k=0.

On aura en conséquence,

(f′⁡(x+k)−f′⁡x)/k=(f₁⁡(x+k, x+k)−f₁⁡(x, x))/k ;

mais, en vertu de la formule (11)

f₁⁡(x+k, x+k)−f₁⁡(x, x)=k{f₂⁡(x, x+k, x+k)+f₂⁡(x, x, x+k)}

donc

(f′⁡(x+k)−f′⁡x)/k=f₂⁡(x, x+k, x+k)+f₂⁡(x, x, x+k);

donc aussi f″⁡x=2f₂⁡(x, x, x); car ce sont là les limites respectives vers lesquelles tendent les deux membres de cette équation, à mesure que k tend vers zéro.

On aura, en conséquence,

(f″⁡(x+k)−f″⁡x)/k=2.(f₂⁡(x+k, x+k, x+k)−f₂⁡(x, x, x))/k ;

mais, en vertu de la formule (11)

f₂⁡(x+k, x+k, x+k)−f₂⁡(x, x, x)=k{f₃⁡(x, x+k, x+k, x+k)+f₃⁡(x, x, x+k, x+k)+f₃⁡(x, x, x, x+k)};

donc

(f″⁡(x+k)−f″⁡x)/k=

2.{f₃⁡(x, x+k, x+k, x+k)+f₃⁡(x, x, x+k, x+k)+f₃⁡(x, x, x, x+k)};

donc aussi f‴⁡x=1.2.3.f₃⁡(x, x, x, x); car ce sont là les limites respectives vers lesquelles tendent les deux membres de cette équation à mesure que k tend vers zéro.

Comme rien n’empêche de poursuivre ce raisonnement aussi loin qu’on voudra, il s’ensuit qu’on a

f₁⁡(x, x)=f′⁡x/1, f₂⁡(x, x, x)=f″⁡x/1.2, f₃⁡(x, x, x, x)=f‴⁡x/1.2.3, …

et, en général,

fₙ⁡(x, x, x, …x, x)=f⁽ⁿ⁾⁡x/n! ;

valeurs qui, substituées dans la formule (12), donneront

f⁡(x+k)=f⁡x+k/1 f′⁡x+k²/1.2 f″⁡x+k³/1.2.3 f‴⁡x+…+kⁿ⁻¹/(n−1)! f⁽ⁿ⁻¹⁾⁡x

+ kⁿ/f ₙ(x, x, x, …x, x, x+k).

On reconnaît ici la série de Taylor, qui se trouve ainsi démontrée.

12. On a cette suite d’équations (1)

f₁⁡(x, x)=f₁⁡(x, x+k)−kf₂⁡(x, x, x+k), f₂⁡(x, x, x+k)=f₂⁡(x, x+k, x+2k)−2kf₃⁡(x, x, x+k, x+2k), f₃⁡(x, x, x+k, x+2k)=f₃⁡(x, x+k, x+2k, x+3k)−3kf₄⁡(x, x, x+k, x+2k, x+3k),

En prenant la somme des produits respectifs de ces équations par +1, −k, +2k², −3k³, +…, et réduisant, il viendra

f₁⁡(x, x)=f₁⁡(x, x+k)−kf₂⁡(x, x+k, x+2k)

+1.2k²f₃⁡(x, x+k, x+2k, x+3k)−…

±(n−1)!kⁿ⁻¹fₙ⁡(x, x+k, x+2k, …x+nk)

±n!kⁿfₙ₊₁⁡(x, x+k, x+2k, …x+nk);

mais on a par la formule (5), en posant y=f⁡x, d’où f₁⁡(x, x)=f′⁡x=d⁡y/d⁡x,

f₁⁡(x, x+k)=Δy/k, f₂⁡(x, x+k, x+2k)=Δ²y/1.2k², f₃⁡(x, x+k, x+2k, x+3k)=Δ³y/1.2.3k³,

il viendra donc, en substituant et multipliant par k

k d⁡y/d⁡x=Δy/1−Δ²y/2+Δ³y/3−…

± Δⁿy/n ∓n!kⁿ⁺¹fₙ₊₁⁡(x, x, x+k, x+2k, …x+nk).

On pourrait étendre indéfiniment ces recherches, mais ce qui précède suffit pour atteindre le but que nous nous étions proposé.

* ↑ Ce mémoire n’ayant point été rédigé par l’auteur, mais seulement d’après des notes très-sommaires qu’il avait fournies, le lecteur voudra bien ne point lui attribuer les négligences de rédaction ou même les erreurs qui pourraient s’y être glissées. J. D. G.
