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

## ANALISE TRANSCENDANTE.

1. Les produits dont les facteurs procèdent suivant une progression arithmétique, et que j’ai nommés facultés numériques, n’ont pas été inutiles au progrès de l’analise. Ils ont servi à exprimer, par un seul terme, et à trouver, d’une manière fort simple, les valeurs numériques de toutes les fonctions transcendantes qui dépendent du cercle, aussi bien que quelques classes, très-nombreuses, d’intégrales définies, Il s’en faut de beaucoup que cette mine soit épuisée. Le langage de l’analise transcendante a été borné, jusqu’ici, aux seules idées de fonctions exponentielles et de fonctions circulaires ; et il est naturel de considérer cette extrême pénurie, comme une des causes principales de l’impossibilité où nous nous trouvions de résoudre le plus grand nombre des problèmes qui se présentaient à nous. Les facultés numériques viennent, fort à propos, pour enrichir ce langage, et pour étendre ainsi le domaine de la science.

2. J’ai prouvé, dans un premier mémoire, que toute faculté était réductible à la forme très-simple 1^{y|1} ou y! ; mais, comme les facultés de cette dernière forme ne dispensent pas de la considération des autres ; afin de faire correspondre une différence de dénomination à une différence de symboles, j’appellerai, à l’avenir, factorielles les fonctions de la forme générale a^{y|r}, et je réserverai exclusivement le nom de facultés numériques, ou simplement de facultés, pour désigner les fonctions de la forme 1^{y|1} ou y!, auxquelles se réduisent les premières, dans le cas particulier où l’on a a=1 et r=1.

3. La factorielle a^{m|r} ou a(a+r)(a+2r)(a+3r)…(a+mr−r)

peut toujours être développée en une série de la forme

aᵐ+Aaᵐ⁻¹r+Baᵐ⁻²r²+…+Mrᵐ.

J’ai fait voir ailleurs que, dans le cas d’un exposant infiniment petit, les coefficiens A, B, C, … devenaient ces nombres même dont l’usage, dans le calcul sommatoire, a été remarqué par leur illustre inventeur Jacques Bernoulli. Mettant dm à la place de m, et désignant par −B₁dm, −B₂dm, −B₃dm, … les valeurs que reçoivent les coefficiens A, B, C, …, dans le cas d’un exposant infiniment petit, on aura

1/2=B₁, 1/3=B₁−2B₂, 1/4=B₁−3B₂+3B₃, 1/5=B₁−4B₂+6B₃−4B₄,

et en général

1/(n+1)=B₁−n/1 B₂+n/1.n−1/2 B₃−n/1.n−1/2.n−2/3 B₄+…± n/1 Bₙ.

En faisant le calcul de ces nombres, on verra que tous ceux d’un indice impair, tels que B₁, B₃, B₅, … sont égaux à zéro, à l’exception du premier B₁ qui est 1/2 ; et que tous ceux d’un indice pair, savoir B₂, B₄, B₆, … sont alternativement positifs et négatifs. Leurs valeurs sont

B₂=+1/12, B₄=−1/120, B₆=+1/252, B₈=−1/240, …

4. Les nombres de Bernoulli nous mènent naturellement aux deux fonctions que j’ai désignées par Λt et Γt. La première Λt, par laquelle nous exprimons la série

B₁t+B₂t²+B₄t⁴+B₆t⁶+…,

sert à trouver la première dériver de la factorielle a^{y|r}, dans laquelle nous regardons l’exposant y comme la variable de la fonction. En faisant, pour abréger, a+ry=t, on a

D

(a^{y|r})=a^{y|r}(Log.t−Λ r/t ).

La seconde Γt, par laquelle nous représentons la série

B₂t+1/3 B₄t³+1/5 B₆t⁵+1/7 B₈t⁷+…,

est liée avec la première, par l’équation linéaire très-simple

B₁t+t²D

Γt=Λt.

Elle est essentielle pour trouver le logarithme naturel de la factorielle a^{y|r} On a en effet,

Log.(a^{y|r})=−y(1−Log.a)+( t/r−1/2 )Log.t/a−Γ r/a+Γ r/t.

5. Le logarithme naturel de la factorielle a^{y|r} que, pour abréger, nous représenterons simplement par Y, est remarquable par la forme de ses dérivées successives. On a d’abord

D

Log.Y=Log.t−Λ r/t :

sur quoi on peut remarquer que c’est l’expression de la somme de fractions

r/a+r/(a+r)+r/(a+2r)+…+r/t−r,

augmentée de Log.a−Λ r/a. Si ensuite, pour abréger, on désigne simplement par Σ rⁿ/tⁿ. la somme infinie de fractions

rⁿ/tⁿ+rⁿ/((t+r)ⁿ)+rⁿ/((t+2r)ⁿ)+rⁿ/((t+3r)ⁿ)+…,

on aura ; en faisant toujours a^{y|r}=Y,

D²

Log.X=+Σ r²/t²,

D³

Log.Y=−2Σ r³/t³,

D⁴

Log.Y=+6Σ r⁴/t⁴,

D⁵

Log.Y=−24Σ r⁵/t⁵,

D⁶

Log.Y=+120Σ r⁶/t⁶,

6. Ces sommes de fractions se trouvent facilement, par les formules connues. On a, en effet,

Σ r²/t² =r/t+r²/2t²+2B₂ r³/t³+4B₄ r⁵/t⁵+6B₆ r⁷/t⁷+…, Σ r³/t³ =r²/2t²+r³/2t³+3B₂ r⁴/t⁴+10B₄ r⁶/t⁶+21B₆ r⁸/t⁸+…, Σ r⁴/t⁴ =r³/3t³+r⁴/2t⁴+4B₂ r⁵/t⁵+20B₄ r⁷/t⁷+56B₆ r⁹/t⁹+…, …………………………

Toutes ces séries peuvent être regardées comme convergentes à volonté, attendu qu’on n’aura qu’à calculer à part quelques-uns des premiers termes de Σ rⁿ/tⁿ, et employer ensuite la formule, pour trouver la somme des autres. Lorsque, dans ces formules, on suppose à r une valeur imaginaire, on est conduit à une suite de théorèmes du plus grand intérêt dans l’analise, et sur lesquels nous reviendrons en son lieu.

7. En désignant par la série

Ay+By²+Cy³+Dy⁴+…,

ordonnée suivant les puissances de l’exposant y, le logarithme naturel de la factorielle a^{y|r}, les coefficienss A, B, C, D, … seront ce que deviennent les dérivées de ce logarithme, dans le cas de y=0, ou de t=a, respectivement divisées par 1, 2, 6, 24, 120, …. On aura ainsi

A=Log.a−Λ r/a, B=+1/2 Σ r²/a², C=−1/3 Σ r³/a³,

D=+1/4 Σ r⁴/a⁴, E=−1/5 Σ r⁵/a⁵, ………………

De Log.a^{y|r} on passe facilement à cette factorielle elle-même ; et si l’on représente par

1+A′y+B′y²+C′y³+D′y⁴+…

la série qui l’exprime, on aura

A′=A, 2B′=AA′+2B, 3C′=AB′+2BA′+3C, 4D′=AC′+2BB′+3CA′+4D, ……………

8. On peut remarquer, au sujet des nombres de Bernonlli, que les séries que nous allons désigner, et dont nous ferons un usage fréquent, dans le calcul sommatoire, sont toutes parfaitement sommables. En faisant, pour abréger, eʸ=x, on a

B₂.y²/1.2+B₄.y⁴/1.2.3.4+B₆.y⁶/1.2.3.4.5.6+…=Log.e¹ᐟ²ʸ−e⁻¹ᐟ²ʸ/y, B₂. y/1 +B₄. y³/1.2.3 +B₆.y⁵/1.2.3.4.5+…=−1/y+1/2.(x+1)/x−1, B₂ +B₄. y²/1.2 +B₆.y⁴/1.2.3.4+…=+1/y²−x/(x−1)², B₄. y/1 +B₆. y³/1.2.3 +B₈.y⁵/1.2.3.4.5+…=−2/y³+(x(x+1))/(x−1)³, B₄ +B₆. y²/1.2 +B₈.y⁴/1.2.3.4+…=+6/y⁴−(x(x²+4x+1))/(x−1)⁴, B₆. y/1 +B₈. y³/1.2.3 +B₁₀.y⁵/1.2.3.4.5+…=−24/y⁵+(x(x³+11x²+11x+1))/(x−1)⁵, B₆ +B₈. y²/1.2 +B₁₀.y⁴/1.2.3.4+…=+120/y⁶−(x(x⁴+26x³+66x²+26x+1))/(x−1)⁶, B₈. y/1 +B₁₀.y³/1.2.3 +B₁₂.y⁵/1.2.3.4.5+…=−720/y⁷+(x(x⁵+57x⁴+302x³+302x²+57x+1))/(x−1)⁷,

et ainsi des autres. La loi que suivent les coefficiens des polynômes fonctions de x qui entrent dans les seconds membres se présente assez naturellement ; on a, par exemple, pour le dernier

57=5.1+2.26, 302=4.26+3.66, 302=3.66+4.26, 57=2.26+5.1;

et ainsi des autres. La seconde de ces séries a été donnée par Euler (Inst. calculi differentialis, part. II, chap. VI, §. 163.). De celle-ci j’ai déduit la première, par intégration, et les autres, par des différentiations successives.

9. Si, dans la première de ces séries, on suppose y=2ϕ√(−1),

on est conduit à celle qui suit :

Log.ϕ/Sin.ϕ=4B₂.ϕ²/1.2−16B₄.ϕ⁴/1.2.3.4+64B₆.ϕ⁶/1.2.3.4.5.6+…

Cette série, peu connue, est peut-être la plus convergente de toutes celles qui font connaître le logarithme du sinus d’un angle proposé. La supposition de y imaginaire, appliquée aux autres séries, conduit aussi à des théorèmes fort intéressans.

10. Ayant trouvé le logarithme naturel de la factorielle a^{y|r} égal à

−y(1−Log.a)+( t/r−1/2 )Log.t/a−Γ r/a+Γ r/t,

la lettre t désignant toujours a+ry, il importe d’examiner ce que devient cette expression, dans le cas d’un exposant imaginaire. Soit donc y=p+iq, la lettre i désignant la racine quarrée de moins un ; on aura

Log.(a^{p+iq|r})=(p+iq)(Log.a−1)+(p−1/2+a/r+iq)Log.(a+pr+iqr)/a

−Γ r/a+Γ r/(a+pr+iqr).

Ici, si l’on fait

k²=(a+pr)²+(qr)²,

Tang.ϕ=qr/(a+pr),

on aura

a+pr+iqr=k(Cos.ϕ+iSin.ϕ),

Log.(a+pr+iqr)=Log.k+iϕ,

le logarithme de a^{p+iq|r} prendra donc la forme d’un binôme M+iN, dans lequel on aura

M=p(−1+Log.a)+(p−1/2+a/r )Log.k/a−qϕ−Γ r/a

+B₂.r/k Cos.ϕ+B₄.r³/3k³ Cos⁡.3ϕ+B₆.r⁵/5k⁵ Cos⁡.5ϕ+…,

N=q(−1+Log.a)+(p−1/2+a/r )ϕ+qLog.k/a

−B₂.r/k Sin.ϕ−B₄.r³/3k³ Sin⁡.3ϕ−B₆.r⁵/5k⁵ Sin⁡.5ϕ−…;

les deux séries sont convergentes à volonté.

11. Ayant ainsi

Log.(a^{p+iq|r})=M+iN,

il est visible qu’on aura

Log.(a^{p−iq|r})=M−iN,

ce qui donnera

Log.(a^{p+iq|r})(a^{p−iq|r})=2M,

Log.(a^{p+iq|r})/a^{p−iq|r}=2iN;

d’où on conclura

(a^{p+iq|r})/a^{p−iq|ra}=e²ⁱᴺ=Cos⁡.2N+iSin⁡.2N.

Dans le cas, très-fréquent, de p=0, lequel donne

k²=a²+q²r², Tang.ϕ=qr/a ;

on aura

M=( a/r−1/2 )Log.k/a−qϕ−Γ r/a+B₂.r/k Cos.ϕ+B₄.r³/3k³ Cos⁡.3ϕ+…,

N=q(−1+Log.a)+( a/r−1/2 )ϕ+qLog.k/a−B₂.r/k Sin.ϕ−B₄.r³/3k³ Sin⁡.3ϕ−…;

et par suite

Log.a^{iq|r}=M+iN,

Log.a^{−iq|r}=M−iN,

d’où

Log.a^{iq|r}×a^{−iq|r}=2M, Log.a^{iq|r}/a^{−iq|r}=2iN.

12. Quelles que soient la base a et la différence r, la factorielle a^{y|r}, dans les deux cas de r positif et de r négatif, qui doivent toujours être soigneusement distingués, sera réductible à la faculté 1^{y|1}, pour laquelle nous avons proposé la notation très-simple y!, devenue nécessaire, par l’usage très-fréquent de ce genre particulier de fonctions, dans la plupart des opérations de haute analise. Nous avons observé, dans un premier mémoire, que le passage des factorielles aux facultés s’exécutait au moyen des deux formules

a^{y|r}=(( a/r−1+y)!)/( a/r−1)! rʸ, a^{y|−r}=(a/r)!/( a/r−y)! rʸ.

La supposition de a=1, r=1 donne aux formules précédemment calculées une très-grande simplicité. On a alors

Log.y!=−y+(y+1/2 )Log.(1+y)−Γ1+Γ 1/(1+y),

formule qu’on peut rendre convergente à volonté, en y introduisant un nombre entier arbitraire h, qu’il suffira de prendre de 4 à 6.

Il viendra ainsi

Log.y!=−h−y−Log.(1+y)^{h|1}+(h+1/2+y)Log.(1+h+y)−Γ1+Γ 1/(1+h+y).

Sur quoi on doit observer qu’il s’agit toujours ici de logarithmes naturels.

13. Si, dans cette supposition de a=r=1, l’exposant y prenait la forme du binôme imaginaire p+iq, en posant alors

k²=(1+p)²+q², Tang.ϕ=q/(1+p) ;

M=−p+(p+1/2 )Log.k−qϕ−Γ1+B₂.Cos.ϕ/k+B₄ Cos⁡.3ϕ/3k³+…,

N=−q+(p+1/2 )ϕ+qLog.k−B₂.Sin.ϕ/k−B₄.Sin⁡.3ϕ/3k³−…,

il viendrait

Log.(p+iq)!=M+iN, Log.(p−iq)!=M−iN.

Posant, de plus, p=0, ce qui donne

k²=1+q², Tang.ϕ=q

M=1/2 Log.k−qϕ−Γ1+B₂.Cos.ϕ/k+B₄ Cos⁡.3ϕ/3k³+…,

N=−q+1/2 ϕ+qLog.k−B₂.Sin.ϕ/k−B₄.Sin⁡.3ϕ/3k³+…,

on aura

Log.(+iq)!=M+iN, Log.(−iq)!=M−iN,

ce qui donnera encore

Log.(+iq)!(−iq)!=2M, Log.((+iq)!)/(−iq)!=−2iN.

14. Le logarithme de y! pouvant toujours être développé en une série de la forme

Ay+By²+Cy³+Dy⁴+…,

on aura dans le cas actuel de a=r=1,

− A=Λ1=B₁+B₂+B₄+B₆+…, +2B=1+1/4+1/9+1/16+…=1/1+1/2+2B₂+4B₄+6B₆+…, −3C=1+1/8+1/27+1/64+…=1/2+1/2+3B₂+10B₄+21B₆+…, +4D=1+1/16+1/81+1/256+…=1/3+1/2+4B₂+20B₄+56B₆+…, −5E=1+1/32+1/243+1/1024+…=1/4+1/2+5B₂+35B₄+126B₆+…, +6F=1+1/64+1/729+1/4096+…=1/5+1/2+6B₂+56B₄+252B₆+…, …………………………………

Les valeurs numériques de toutes ces sommes de puissances sont connues et calculées ; quant à celle de Λ1, elle est

0, 57721 56649…

On sait de plus que les sommes à indice pair sont réductibles aux puissances paires de ϖ ; d’où l’on obtient

B=+2B₂.ϖ²/1.2, D=−2³B₄.ϖ⁴/1.2.3.4, F=+2⁵B₆.ϖ⁶/1.2.3.4.5.6, H=−2⁷B₈.ϖ⁸/1.2.3.4.5.6.7.8, ………………

Quant aux autres coefficiens A, C, E, … l’analise ne nous offre pas les mêmes moyens de les obtenir ; il faudrait, pour y parvenir, interpoler la série des nombres de Bernoulli, d’après une loi probablement fort simple, mais que nous ne connaissons pas encore.

15. Ayant

Log.(+y)!=Ay+By²+Cy³+Dy⁴+…,

on doit avoir

Log.(−y)!=−Ay+By²−Cy³+Dy⁴−…,

de là résulte

Log.(+y)!×(−y)!=2By²+2Dy⁴+2Fy⁶+…;

en changeant y en iy, on aura de même

Log.(+iy)!×(−iy)!=−2By²+2Dy⁴−2Fy⁶+…;

16. On a vu (8) que

B₂.y²/1.2+B₄.y⁴/1.2.3.4+B₆.y⁶/1.2.3.4.5.6+…=Log.e¹ᐟ²ʸ−e⁻¹ᐟ²ʸ/y,

En remplaçant y par iy, on aura, après les réductions connues,

−B₂.y²/1.2+B₄.y⁴/1.2.3.4−B₆.y⁶/1.2.3.4.5.6+…=Log.( 2/y Sin.y/2 ).

Les réductions appliquées aux valeurs des logarithmes de (+y)!(−y)! et de (+iy)!(−iy)! qu’on vient de trouver, conduisent aux deux théorèmes très-importans qui suivent :

THÉORÈME I. (+y)!(−y)!=ϖy/Sin.ϖy.

THÉORÈME II. (+iy)!(−iy)!=2ϖy/e^{ϖy}−e^{−ϖy}.

Les principes généraux étant posés, proposons-nous la solution générale des deux problèmes qui suivent :

17. PROBLÈME I. Évaluer numériquement le produit

{1−x²/a² }{1−x²/((a+r)²) }{1−x²/((a+2r)²) }{1−x²/((a+3r)²) }…

continué à l’infini ?

Solution. Ce produit se décompose dans ceux-ci

a−x/a.(a−x+r)/(a+r).(a−x+2r)/(a+2r).(a−x+3r)/(a+3x) …,

(a+x)/a.(a+x+r)/(a+r).(a+x+2r)/(a+2r).(a+x+3r)/(a+3x) …;

dans mon Analise des réfractions, je les ai réduit respectivement à

(a−x)^{x/r|r}/(infin.r)ˣᐟʳ, ((a+x)^{−x/r|r})/(infin.r)⁻ˣᐟʳ,

ce qui rend le produit cherché égal au simple produit des deux factorielles

(a−x)^{x/r|r}, (a+x)^{−x/r|r};

il ne reste donc plus qu’à réduire ces factorielles aux facultés, ce qui se fait à l’aide des formules ci-dessus (12). En posant, pour abréger, a/r−1=f, on trouvera

(a−x)^{x/r|r}=f!/(f−x/r )! rˣᐟʳ, (a+x)^{−x/r|r}=f!/((f+x/r )!) r⁻ˣᐟʳ;

d’où on conclura

(a−x)^{x/r|r}⋅(a+x)^{−x/r|r}=f!f!/((f−x/r )!(f+x/r )!) ;

la solution du problème proposé sera donc réduite à la détermination des trois facultés f!, (f−x/r )!, (f+x/r )! dont on trouvera les valeurs numériques toutes calculées, dans la table donnée à la page 6 de ce volume.

18. PROBLÈME II. Évaluer numériquement le produit

{1+x²/a² }{1+x²/((a+r)²) }{1+x²/((a+2r)²) }{1+x²/((a+3r)²) }…

continué à l’infini ?

Solution. En continuant de faire a/r−1=f, il suffira de remplacer x par ix, dans la formule qu’on vient de trouver. Le produit demandé deviendra égal à

f!f!/((f−ix/r )!(f+ix/r )!).

Les logarithmes des deux facultés du dénominateur sont réductibles aux formes M−iN et M+iN, ce qui rend le logarithme de leur produit égal à 2M. Il serait fort à désirer que quelque calculateur courageux voulût calculer les binômes M+iN=Log.(iy)!, pour toutes les valeurs de y depuis 0 jusqu’à 1, de même que nous devons à M. le professeur Bessel une table des logarithmes de y!, dans le cas d’une base réelle. En attendant, la série (13), qui a l’avantage d’être convergente à volonté, nous fournit un moyen très-expéditif de trouver la valeur numérique de 2M, logarithme du produit, (f−ix/r )!(f+ix/r )!. Il faudra, pour en faire usage, déterminer l’angle ϕ et le coefficient k de manière qu’on ait

Tang.ϕ=x/a, k²=(a²+x²)/r²,

ce qui donne

k=a/rCos.ϕ ;

et on aura

M=−a/r+1+( a/r−1/2 )Log.k−x/r ϕ−Γ1+B₂ Cos.ϕ/k+B₄ Cos⁡.3ϕ/3k³+…,

19. Appliquant la solution de ces deux problèmes au cas particulier de a=r=ϖ, qui donne f=0, et qui rend (5) la variable y=x/ϖ, d’où ϖy=x, on sera conduit aux théorèmes très-connus, démontrés par EULER (Introd. in. anali. infin., 1.ʳᵉ partie, n.ᵒˢ 156 et suivans.) ; savoir :

Sin.ϕ=ϕ(1−ϕ²/ϖ² )(1−ϕ²/4ϖ² )(1−ϕ²/9ϖ² )…,

e^{ϕ}−e^{−ϕ}/2=ϕ(1−ϕ²/ϖ² )(1−ϕ²/4ϖ² )(1−ϕ²/9ϖ² )…

20. Si dans (17) on fait a=ϖ/2, r=ϖ, on aura, d’un côté, le produit

(1−4x²/ϖ² )(1−4x²/9ϖ² )(1−4x²/25ϖ² )(1−4x²/49ϖ² )…

continué à l’infini, lequel, par conséquent, sera égal à ce que devient la fraction

f!f!/((f−x/r )!(f+x/r )!),

par cette supposition qui donne

f=−1/2=−1+1/2, f+x/r =−1/2+x/r=−( 1/2−x/r ), f−x/r =−1/2−x/r=−1+( 1/2−x/r ).

La faculté f! devient ainsi 2(1/2)!=√ϖ. La faculté (f−x/r )! deviendra

( 1/2−x/ϖ )!/ 1/2−x/ϖ ,

et le produit des deux facultés ( 1/2+x/ϖ ) et !( 1/2−x/ϖ )! sera (16)

ϖ/2−x/Sin.( ϖ/2−x)= ϖ/2−x/Cos.x.

Après avoir employé toutes ces réductions, on sera conduit au théorème très-connu,

Cos.φ=(1−4φ²/ϖ² )(1−4φ²/9ϖ² )(1−4φ²/25ϖ² )…

21. Et si, dans cette dernière formule, on change ϕ en iφ, elle deviendra

e^{ϕ}−e^{−ϕ}/2=(1+4φ²/ϖ² )(1−4φ²/9ϖ² )(1−4φ²/25ϖ² )(1−4φ²/49ϖ² )…;

formule connue depuis l’analise d’Euler.

22. Les formules (17) et (19) nous conduisent aux deux théorèmes qui suivent

h!(1−h)!=h(1−h)/Sin.hϖ.ϖ;

( 1/2−h)!( 1/2+h)!=(( 1/2−h)( 1/2+h))/Cos.hϖ.ϖ.

Ces deux formules, qui sont identiques entre elles, procureront à ceux qui voudront s’occuper de la construction d’une table des facultés, pour les fractions décimales comprises entre 0 et 1, l’avantage précieux de diminuer leur travail de moitié, en ne les obligeant à le pousser que jusqu’à h=1/2=0, 5; ce qui donnera, en outre, une grande convergence à la série qu’on est obligé d’employer pour le calcul de cette table. En changeant h en ih, on aura pareillement

(ih)!(1−ih)!=2h(1−ih)/e^{hϖ}−e^{−hϖ}.ϖ;

( 1/2−ih)!( 1/2+ih)!=(2( 1/2−ih)( 1/2+ih))/(e^{hϖ}+e^{−hϖ}).ϖ=(2( 1/4+h²))/(e^{hϖ}+e^{−hϖ}).ϖ.

23. Le théorème suivant mérite d’être remarqué ; il concerne le produit x!y! de deux facultés dans lesquelles la somme x+y des exposans est un nombre entier quelconque, pair ou impair.

Dans le cas d’une somme paire, soient x=a+h, y=a−h, et par conséquent x+y=2a ; la lettre a pourra alors désigner un nombre entier quelconque, et h une fraction moindre que l’unité. Cela posé, on a

(a+h)!=(+h)!(1+h)^{a|1},

(a−h)!=(−h)!(1−h)^{a|1},

ou

(a+h)!=(+h)!(1+h)(2+h)(3+h)…(a+h),

(a−h)!=(−h)!(1−h)(2−h)(3−h)…(a−h),

on aura donc (16)

x!y!=hϖ/Sin.hϖ (1−h²)(4−h²)(9−h²)…(a²−h²).

Dans le cas d’une somme impaire, soient posés

x=a+1/2+h, y=a+1/2−h,

ce qui donnera

x+y=2a+1,

a pouvant désigner un nombre entier quelconque. On aura alors

x!=(a+1/2+h)!=( 1/2+h)!( 3/2+h)^{a|1},

y!=(a+1/2−h)!=( 1/2+h)!( 3/2−h)^{a|1},

ou

x!=( 1/2+h)( 3/2+h)( 5/2+h)( 7/2+h)…( (2a+1)/2+h),

y!=( 1/2−h)( 3/2−h)( 5/2−h)( 7/2−h)…( (2a+1)/2−h);

donc (22)

x!y!=ϖ/Cos.hϖ ( 1/4−h²)( 9/4−h²)( 25/4−h²)( 49/4−h²)…{ ((2a+1)²)/4−h²}.

Si, dans le cas de x+y pair, h se change en ih, ce qui donnera toujours x+y=2a, a étant un nombre entier quelconque, il viendra

x!y!=2hϖ/e^{hϖ}−e^{−hϖ} (1+h²)(4+h²)(9+h²)…(a²+h²).

Si le même changement arrive, dans le cas de x+y impair, en sorte qu’on ait toujours x+y=2a+1, a étant un nombre entier quelconque, il viendra

x!y!=2ϖ/(e^{hϖ}+e^{−hϖ}) ( 1/4+h²)( 9/4+h²)( 25/4+h²)( 49/4+h²)…{ ((2a+1)²)/4+h²}.

24. L’analise que nous venons de développer n’est nullement bornée au cas proposé ; et si l’on demandait soit la valeur du produit

{1+xⁿ/aⁿ }{1+xⁿ/((a+r)ⁿ) }{1+xⁿ/((a+2r)ⁿ) }{1+xⁿ/((a+3r)ⁿ) }…,

soit celle du produit

{1−xⁿ/aⁿ }{1−xⁿ/((a+r)ⁿ) }{1−xⁿ/((a+2r)ⁿ) }{1−xⁿ/((a+3r)ⁿ) }…,

continué à l’infini, on la trouverait encore, en suivant rigoureusement les mêmes principes.

25. Jusqu’ici nous avons supposé que les facteurs de nos factorielles constituaient toujours une progression arithmétique du premier ordre ; c’est-à-dire, une progression ayant ses premières différences constantes ; et ces sortes de fonctions peuvent être appelées Factorielles du premier ordre. On peut aussi imaginer une suite de facteurs constituant une progression arithmétique du second ordre ; c’est-à-dire, une progression ayant ses secondes différences constantes, telle que

a, a+b, a+2b+c, a+3b+3c, a+4b+6c, a+5b+10c, ………

Le terme qui répondrait à l’indice n serait alors,

a+n−1/1 b+n−1/1.n−3/2 c.

Un semblable produit pourrait être appelé Factorielles du second ordre ; et, pour peu qu’on suive le développement de la plupart de nos séries, on verra que ces factorielles, de même que celles des ordres supérieurs, c’est-à-dire, celles dans lesquelles ce sont les différences d’un ordre plus élevé que le second, qui sont constantes, doivent se rencontrer très-fréquemment. Heureusement toutes ces factorielles sont réductibles à celles du premier ordre, moyennant une décomposition analitique fort simple. On peut toujours, en effet, pour le second ordre, déterminer les deux premiers termes A, A′, et les deux premières différences r, r′, de manière que le terme général

a+n−1/1 b+n−1/1.n−3/2 c.

devienne équivalent au produit

[A+(n−1)r][A′+(n−1)r′],

indépendamment de l’indice n. Il faudra, pour cela, résoudre les trois équations

AA′=a, Ar′+A′r=b−1/2 c, rr′=1/2 c;

on aura alors

a=AA′, a+b=(A+r)(A′+r′), a+2b+c=(A+2r)(A′+2r′), …………………

en sorte que la factorielle proposée du second ordre deviendra le simple produit

A^{n|r}×A^{′n|r′}

* ↑ Voyez la page première de ce volume.
* ↑ Voyez Élémens d’arithmétique universelle, page 360, n.ᵒˢ 557 et suivans.
* ↑ La lettre D est employée ici comme signe de dérivation ; en sorte qu’en général D ϕ(x)=d.ϕ(x)/dx.
* ↑ Depuis l’impression de la table de M. Bessel, nous nous sommes aperçu que M. Legendre dans ses Exercices de calcul intégral (Paris, 1811), avait publié une table du même genre, et nous venons d’apprendre qu’une pareille table venait aussi d’être calculée par M. Gauss. Voilà donc trois géomètres du premier ordre qui, faute de moyens rapides de communication, ont consommé un temps précieux en de pénibles calculs, pour parvenir aux mêmes résultats. M. Legendre, dans sa table, désigne par a−1 ce que M. Kramp représente par y, et par, ce que M. Kramp désigne par y. Cette table, calculée par une méthode analogue à celle qui a été indiquée dans la note de la page 12 de ce volume, ne contient les logarithmes de la faculté y! qu’à sept décimales seulement, et encore la septième décimale n’y est pas toujours exacte ; on n’y trouve pas non plus les différences des logarithmes qu’elle renferme ; mais ces logarithmes y sont calculés pour les valeurs de y de millième en millième, ce qui rend à la fois les interpolations plus faciles et moins fréquemment nécessaires. Les détails dans lesquels entre M. Legendre, sur le calcul de cette table, et sur la nature, les propriétés et les usages des nombres qu’elle renferme, sont d’ailleurs du plus grand intérêt. J. D. G.
* ↑ Les inconnues A, A′, r, r′ de ces trois équations étant au nombre de quatre, on pourra disposer de l’une d’elles pour rendre les valeurs des autres les plus simples possibles, J. D. G.
* ↑ Ce que M. Kramp appelle ici Factorielles de différens ordres est ce que Vandermonde avait déjà appelé Puissances de différens ordres, avec cette différence seulement que les factorielles de l’ordre n sont des puissances de l’ordre n+1. Ainsi, suivant le langage de M. Kramp, les puissances du premier ordre, c’est-à-dire, les simples puissances que l’on considère dans les élémens, sont des factorielles de l’ordre zéro. J. D. G.
