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

## ANALISE TRANSCENDANTE.

### Formules nouvelles, pour l’intégration approchée de toute fonction différentielle d’une seule variable, entre deux limites données quelconques ;

L’objet que nous nous proposons dans ce mémoire est d’enseigner à déterminer, entre des limites données quelconques, l’intégrale de toute différentielle de la forme Xd⁡x, quelle que puisse être d’ailleurs la forme de la fonction de x désignée par X. La méthode que nous allons faire connaître a cela de particulier qu’elle est, en quelque sorte, étrangère aux principes du calcul intégral et à la notion des infinimens petits ; elle ne suppose que les principes connus de l’algèbre élémentaire ; elle s’étend à toutes les fonctions quelconques, à celles même qui se sont constamment refusées jusqu’ici à tous les moyens d’intégration connus ; elle donne l’intégrale demandée, moyennant un nombre très-limité de termes, avec une précision bien supérieure à tout ce qu’on pourrait se promettre de l’usage des suites infinies.

1. On sait que l’intégration de toute formule Xd⁡x, entre des limites données, x=a et x=a′, par exemple, revient à quarrer l’aire mixtiligne terminée d’une part par la courbe dont l’équation serait y=X, d’une autre par l’axe des x, et enfin par les ordonnées b, b′ de cette courbe répondant respectivement aux abscisses a et a′, C’est même de là que ce problème a été appelé problème des quadratures, et c’est sous ce point de vue que nous l’envisagerons constamment, dans tout ce qui va suivre.

2. Soit fait, pour abréger a′−a=c ; et, pour fixer les idées, imaginons que l’on ait divisé l’intervalle c en douze parties égales ; désignons par a₀, a₁, a₂, …a₁₀, a₁₁, a₁₂ les abscisses qui répondent aux treize points de divisions ; au moyen de l’équation y=X, nous pourrons calculer les ordonnées qui leur correspondent ; représentons-les respectivement par b₀, b₁, b₂, …b₁₀, b₁₁, b₁₂; nous connaîtrons ainsi treize points de la courbe qu’il s’agit de quarrer entre les limites x=a₀ et x=a₁₂, pour lesquelles on a respectivement y=b₀, y=b₁₂.

3. Soient joints les deux points extrêmes (a₀, b₀), (a₁₂, b₁₂) par une corde, cette corde, avec sa projection c et les deux ordonnées extrêmes formera un trapèze ; en désignant son aire par S₁₂, et posant

12( 1/2 b₀+1/2 b₁₂)=b₁₂′,

nous aurons

S₁₂=1/12 cb₁₂′.

4. Soient joints consécutivement les trois points (a₀, b₀), (a₆, b₆), (a₁₂, b₁₂) par deux cordes ; ces cordes formeront, avec c et les trois ordonnées b₀, b₆, b₁₂, deux trapèzes ; en désignant la somme de leurs aires par S₆, et posant

6( 1/2 b₀+b₆+1/2 b₁₂)=b₆′,

nous aurons

S₆=1/12 cb₆′.

5. Soient joints consécutivement les quatre points (a₀, b₀), (a₄, b₄), (a₈, b₈), (a₁₂, b₁₂) trois cordes ; ces cordes formeront, avec c et les quatre ordonnées b₀, b₄, b₈, b₁₂, trois trapèzes ; en désignant la somme de leurs aires par S₄, et posant

4( 1/2 b₀+b₄+b₈+1/2 b₁₂)=b₄′,

nous aurons

S₄=1/12 cb₄′.

6. Soient joints consécutivement les cinq points (a₀, b₀), (a₃, b₃), (a₆, b₆), (a₉, b₉), (a₁₂, b₁₂) par quatre cordes ; ces cordes formeront, avec c et les cinq ordonnées b₀, b₃, b₆, b₉, b₁₂, quatre trapèzes ; en désignant la somme de leurs aires par S₃, et posant

3( 1/2 b₀+b₃+b₆+b₉+1/2 b₁₂)=b₃′,

nous aurons

S₃=1/12 cb₃′.

7. Soient joints consécutivement les sept points (a₀, b₀), (a₂, b₂), (a₄, b₄), (a₆, b₆), (a₈, b₈), (a₁₀, b₁₀)(a₁₂, b₁₂) par six cordes ; ces cordes formeront, avec c et les sept ordonnées b₀, b₂, b₄, b₆, b₈, b₁₀, b₁₂, six trapèzes ; en désignant la somme de leurs aires par S₂, et posant

2( 1/2 b₀+b₂+b₄+b₆+b₈+b₁₀+1/2 b₁₂)=b₂′,

nous aurons

S₂=1/12 cb₂′.

8. Enfin, soient joints consécutivement tous les treize points de la courbe par douze cordes ; ces cordes formeront, avec c et les treize ordonnées, douze trapèzes ; en désignant la somme de leurs aires par S₁, et posant

1/2 b₀+b₁+b₂+b₃+b₄+b₅+b₆+b₇+b₈+b₉+b₁₀+b₁₁+1/2 b₁₂=b₁′,

nous aurons

S₁=1/12 cb₁′.

9. Aucune des aires S₁₂, S₆, S₄, S₃, S₂, S₁,

n’est l’aire demandée mais il résulte évidemment de notre procédé que ces aires convergent de plus en plus vers celle-là. Donc aussi la ligne par laquelle il faut multiplier 1/12 c pour avoir l’aire demandée n’est aucune des lignes b₁₂′, b₆′, b₄′, b₃′, b₂′, b₁′,

mais une ligne vers laquelle celles-là convergent de plus en plus.

10. Portons les ordonnées b₁₂′, b₆′, b₄′, b₃′, b₂′, b₁′,

sur leurs correspondantes b₁₂, b₆, b₄, b₃, b₂, b₁;

et imaginons une courbe située au-dessus de la première, passant par les six points (a₁, b₁′), (a₂, b₂′), (a₃, b₃′), (a₄, b₄′), (a₆, b₆′), (a₁₂, b₁₂′); cette courbe prolongée rencontrera le prolongement de l’ordonnée b₀ en quelque point, en désignant par b₀′ son ordonnée qui répond à celle-là, et conséquemment à l’abscisse a₀ les ordonnées b₁₂′, b₆′, b₄′ , b₃′, b₂′, b₁′,

tendant continuellement vers la ligne par laquelle il faut multiplier 1/12 c pour avoir l’aire cherchée, en désignant cette aire par S₀, nous pourrons prendre sensiblement

S₀=1/12 cb₀′.

et tout se réduira à trouver b₀′ ; problème qui rentre dans les méthodes connues d’interpolation. Par la nature même de ces méthodes, et de l’espèce d’arbitraire auquel elles sont inévitablement assujetties, la valeur que nous trouverons pour b₀′ ne sera point proprement la véritable ; mais sa différence avec elle sera comparable à celle qui existe entre le rayon et le sinus-verse d’un très-petit angle, tel que serait, par exemple, celui d’une minute ou même d’une seconde. Effectivement nous verrons bientôt que, dans tous les cas ordinaires d’intégration, cette différence n’est sensible qu’à la dixième ou à la douzième décimale. D’ailleurs on peut la diminuer à volonté, en augmentant le nombre des parties égales de c qu’on pourra porter à 18, 24, 30, 36, 48 ou 60 au lieu de 12.

11. Il est facile de voir, par la nature de la courbe dont les ordonnées sont b₁₂′, b₆′, b₄′, b₃′, b₂′, b₁′,

qu’elle doit couper perpendiculairement l’ordonnée b₀′, c’est-à-dire, en d’autres termes, qu’en prenant y′ pour le symbole général des ordonnées de cette courbe, et faisant répondre l’origine à l’ordonnée b₀′, on doit avoir en même temps x=0, d⁡y′/d⁡x=0; ce qui exîge que l’expression de y′ ne renferme point la première puissance de x. Pour plus d’uniformité, nous en exclurons également toutes les autres puissances impaires, et nous poserons simplement

y′=A+Bx²+Cx⁴+Dx⁶+Ex⁸+…;

il s’agira donc de déterminer le coefficient A, auquel se réduit y′ lorsque x=0.

12. En prenant 1/12 c pour unité, il faudra donc qu’aux valeurs 1, 2, 3, 4, 6, 12 de x, répondent pour y′ les valeurs b₁′, b₂′, b₃′, b₄′, b₆′, b₁₂′, ce qui donnera

b₁′ =A+ B+ C+ D+ E+ F, b₂′ =A+ 2²B+ 2⁴C+ 2⁶D+ 2⁸E+ 2¹⁰F, b₃′ =A+ 3²B+ 3⁴C+ 3⁶D+ 3⁸E+ 3¹⁰F, b₄′ =A+ 4²B+ 4⁴C+ 4⁶D+ 4⁸E+ 4¹⁰F, b₆′ =A+ 6²B+ 6⁴C+ 6⁶D+ 6⁸E+ 6¹⁰F, b₁₂′ =A+12²B+12⁴C+12⁶D+12⁸E+12¹⁰F;

et, en éliminant, entre ces six équations, les cinq coefficiens B, C, D, E, F, la valeur de A que l’on tirera de l’équation finale, en fonction de b₁′, b₂′, b₃′, b₄′, b₆′, b₁₂′, sera, pour c=12, l’intégrale demandée ; nous avons vu d’ailleurs qu’on a

b₁′=1/2 b₀+ b₁+ b₂+ b₃+b₄+b₅+b₆+b₇+b₈+b₉+b₁₀+b₁₁+1/2 b₁₂, b₂′= b₀+ b₂+ b₄+ b₆+b₈+b₁₀+b₁₂, b₃′=1/2 b₀+3b₃+3b₆+3b₉+1/2 b₁₂, b₄′= 2b₀+4b₄+4b₈+2b₁₂, b₆′= 3b₀+6b₆+3b₁₂,

b₁₂′=6b₀+6b₁₂.

13. Les équations qu’il s’agit de résoudre étant au nombre de six, nous allons les présenter sous la forme plus générale que voici

a′=A+Ba+Ca²+Da³+Ea⁴+Fa⁵, b′=A+Bb+Cb²+Db³+Eb⁴+Fb⁵, c′=A+Bc+Cc²+Dc³+Ec⁴+Fc⁵, d′=A+Bd+Cd²+Dd³+Ed⁴+Fd⁵, e′=A+Be+Ce²+De³+Ee⁴+Fe⁵, f′=A+Bf+Cf²+Df³+Ef⁴+Ff⁵.

Les quantités a′, b′, c′, d′, e′, f′, de même que a, b, c, d, e, f, sont regardées comme données, et il s’agit uniquement d’obtenir la valeur de A ; de sorte que les cinq autres quantités B, C, D, E, F sont tout à fait indifférentes au problème qui nous occupe. Or, on trouve

A=bcdefa′/(b−a)(c−a)(d−a)(e−a)(f−a) +cdefab′/(c−b)(d−b)(e−b)(f−b)(a−b) +defabc′/(d−c)(e−c)(f−c)(a−c)(b−c) +efabcd′/(e−d)(f−d)(a−d)(b−d)(c−d) +fabcde′/(f−e)(a−e)(b−e)(c−e)(d−e) +abcdef′/(a−f)(b−f)(c−f)(d−f)(e−f)

La loi que suit cette expression générale est évidente, et on peut aisément l’étendre au cas où l’on aurait un plus grand nombre d’équations.

14. J’appellerai diviseur général le nombre des parties égales dans lesquelles on aura divisé l’intervalle c qui sépare les ordonnées extrêmes qui terminent l’espace mixtiligne qu’il s’agit de quarrer ; nombre qui a constamment été supposé 12 dans ce qui précède. Le choix de ce nombre n’est point indifférent ; et à grandeur à peu près égale, on doit donner la préférence à celui qui a le plus grand nombre de petits diviseurs, tel que 6, 12, 18, 24, 30, 36, 48, 60, …. Nous allons voir, au surplus, que, dans les applications pratiques, il doit être à peu près superflu d’aller au-delà de 24; attendu qu’en se bornant à ce nombre, on peut, dans les cas ordinaires, obtenir les intégrales avec douze chiffres décimaux exacts, au moins.

15. Le diviseur général étant choisi, le nombre et la nature des parties aliquotes à employer sont encore arbitraires. Il convient de ne jamais donner l’exclusion aux aliquotes 1, 2, 3; et le plus exact sera de les employer toutes ; mais il en résultera nécessairement plus de peine pour le calculateur ; d’ailleurs en n’allant pas même au-delà de 6, on peut obtenir des résultats qui, pour la précision, excèdent déjà les besoins ordinaires de l’analise.

16. Première formule. Prenons d’abord pour diviseur général le nombre 6, en employant tous les aliquotes 1, 2, 3, 6 ; nous aurons simplement ici

A=bcda′/(b−a)(c−a)(d−a)+cdab′/(c−b)(d−b)(e−b)+dabc′/(d−c)(a−c)(b−c)+abcd′/(a−d)(b−d)(c−d)

Or, on a, dans le cas actuel, a=1, b=4, c=9, d=36, ce qui donne, en substituant,

A=(1296a′−567b′+112c′−d′)/840 ;

et, en prenant pour unité l’intervalle entre deux ordonnées consécutives, ce sera là l’intégrale demandée.

17. Pour plus de simplicité, appelons les sept ordonnées du cas actuel α, β, γ, δ, ε, ζ, η;

nous aurons

a′=1/2 α+β+γ+δ+ε+ζ+1/2 η, b′=α+2γ+2ε+η, c′=1/2 α+3δ+1/2 η, d′=3α+3η;

qui donnera ; en substituant,

A=(246(α+η)+1296(β+ζ)+162(γ+ε)+1632δ)/840.

En prenant donc l’intervalle entier qui sépare les ordonnées extrêmes α et η pour unité ; on aura finalement pour l’intégrale cherchée

∫Xd⁡x=(41(α+η)+216(β+ζ)+27(γ+ε)+272δ)/840.

Si, dans cette dernière formule, on fait α=β=γ=δ=ε=ζ=η=1, elle devient ∫Xd⁡x=1 ainsi que cela doit être.

18. Exemple premier. On demande le logarithme naturel de deux ?

Le logarithme d’un nombre quelconque n est l’intégrale de d⁡x/x depuis x=1 jusqu’à x=n. En divisant donc en six parties égales l’intervalle compris entre un et deux, et remarquant qu’ici X=1/x nous aurons d’abord

α=6/6=1, 00000000, α+η=1, 50000000, β=6/7=0.85714286, β+ζ=1, 40259740, γ=6/8=0, 75000000, γ+ε=1, 35000000, δ=6/9=0, 66666667, 2δ=1, 33333333.

ε=6/10=0, 60000000, ζ=6/11=0, 54545455, η=6/12=0, 50000000;

Ces valeurs étant substituées dans la formule (I), on aura

Log⁡.2=0, 69314806.

La valeur rigoureuse est

Log⁡.2=0, 69314718.

La différence est donc +0, 00000088, moindre qu’un millionième.

19. Exemple II. On demande la longueur du demi-quadrans =ϖ/4 ?

La longueur de l’arc dont la tangente est t est l’intégrale de d⁡x/(1+x²), prise depuis x=0 jusqu’à x=t ; celle de l’arc =ϖ/4 sera donc cette même intégrale, prise entre zéro et un ; divisant donc cet intervalle en six parties égales, et remarquant qu’ici X=1/(1+x²), il viendra

α=36/36=1, 00000000, α+η=1, 50000000, β=36/37=0, 97297297, β+ζmm=1, 56313690, γ=36/40=0, 90000000, γ+ε=1, 59230769, δ=36/45=0, 80000000, 2δ=1, 60000000.ε=36/52=1, 69230769, ζ=36/62=0, 59016393, η=36/72=0, 50000000;

Ces valeurs étant substituées dans la formule (I), on aura

=ϖ/4=0, 78539271.

La valeur rigoureuse est

ϖ/4=0, 78539816.

Ainsi l’erreur est −0, 00000545, moindre que un cent millième.

20. Cette première formule (17) est la plus simple et la plus aisée de toutes ; c’est celle qui exige le moins de calculs ; mais c’est aussi c’elle qui donne les résultats les moins approchés. Celles qui vont suivre seront beaucoup plus exactes.

21. Deuxième formule. Prenons 12 pour diviseur général, mais n’admettons d’abord que les parties aliquotes 1, 2, 3, 4, 6 ; cela donnera

A=bcdea′/(b−a)(c−a)(d−a)(e−a)+cdeab′/(c−b)(d−b)(e−b)(a−b)+deabc′/(d−c)(e−c)(a−c)(b−c) +eabcd′/(e−d)(a−d)(b−d)(c−d)+abcde′/(a−e)(b−e)(c−e)(d−e).

Or, nous avons ici a=1, b=4, c==9, d=16, e==36, nous aurons donc, en prenant pour unité l’intervalle entre les deux ordonnées extrêmes,

A=(1728a′−945b′+320c′−54d′+e′)/12600 ;

mais on a, dans le cas actuel,

a′=1/2 α+β+γ+δ+ε+ζ+η+θ+ι+ϰ+λ+μ+1/2 ν, b′=α+2γ+2ε+2η+2ι+2λ+ν, c′=1/2 α+3δ+3η+2ϰ+1/2 ν, d′=2α+4η+4ι+2ν;e′=3α+6η+3ν;

qui donnera, en substituant,

∫Xd⁡x=(49(α+ν)+288(β+μ)−27(γ+λ)+448(δ+ϰ)−63(ε+ι)+288(ζ+θ)+134η)/2100.

22. Exemple I. On demande de nouveau le logarithme de deux ?

En divisant en douze parties égales l’intervalle entre un et deux, nous aurons

α=12/12=1, 0000000000, α+ν=1, 5000000000, β=12/13=0, 9230769231, β+μ=1, 4448160535, γ=12/14=0, 8571428571, γ+λ=1, 4025974026, δ=12/15=0, 8000000000, δ+ϰ=1, 3714285714, ε=12/16=0, 7500000000, ε+ι=1, 3500000000, ζ=12/17=0, 7058823529, ζ+θ=1, 3374613003, η=12/18=0, 6666666667, 2η=1, 3333333333;θ=12/19=1, 6315789474, ι=12/20=0, 6000000000, ϰ=12/21=0, 5714285714, λ=12/22=0, 5454545455, μ=12/23=0, 5217391304, ν=12/24=0, 5000000000;

ce qui donnera, en substituant dans la formule (II)

Log⁡.2=0.6931471816;

la valeur rigoureuse est

Log⁡.2=0, 6931471806;

l’erreur est donc

0, 0000000010,

c’est-à-dire, que cette valeur est exacte dans les huit premiers chiffres décimaux,

23. Remarque. Avec le seul logarithme de deux on peut facilement trouver tous les autres. Soit m un nombre absolument quelconque dont il faille chercher le logarithme. Soit 2ⁿ la puissance de deux qui lui est immédiatement inférieure ; et soit 1+h le quotient qu’on obtient en divisant le premier de ces deux nombres par le second. Comme h sera certainement un nombre moindre que l’unité, la formule (I) suffira pour déterminer le logarithme de 1+h; divisant donc à en six parties égales on aura α=6/6, β=6/(6+h), γ=6/(6+2h), δ=6/(6+3h), ε=6/(6+4h), ζ=6/(6+5h), η=6/(6+6h)=1/(1+h). Il faudra multiplier par h l’intégrale obtenue par la formule ; on aura alors le logarithme de 1+h, auquel ajoutant n fois celui de deux, on aura celui de m avec une erreur qui ne tombera pas au-dessus de la huitième ou même de la neuvième décimale.

24. Exemple II. On demande le logarithme naturel de 10000 ?

La puissance de deux immédiatement inférieure à 10000 est 8192=2¹³. On aura ainsi m=10000, n=13, 1+h=10000/8192=1+1808/8192 =1+113/512 ; on aura donc

Log⁡.10000=13Log⁡.2+Log⁡(1+113/512).

Pour trouver ce dernier logarithme, on fera α=3072/3072, β=3072/3185,

γ=3072/3298, δ=3072/3411, ε=3072/3524, ζ=3072/3637, η=3072/3350

On trouvera ensuite

41(α+η)= 74, 58720000, 216(β+ζ)=390, 78145008, 27(γ+ε)= 48, 68667756, 272 δ=244, 96745680.

La somme 759, 02278444 de ces quatre nombres, divisée par 840, et multipliée par h=113/512 donne pour le logarithme de 1+h ou h=625/512

0, 1994270243. La valeur rigoureuse est 0, 1994270243; la différence est donc seulement de quatre unités décimales du dixième ordre.

25. D’un autre côté, ayant trouvé Log⁡.2=0, 6931471816, on aura 13Log⁡.2=9, 0109133608. Ajoutant celui qu’on vient de trouver, il en résultera le logarithme de 10000 égal à 9, 2103403851. En prenant le quart de ce logarithme, on aura, pour le module de nos tables vulgaires

Log⁡.10=2, 3025850963.

La valeur rigoureuse étant

Log⁡.10=2, 3025850930;

on voit que l’erreur est au-dessous de quatre unités décimales du neuvième ordre.

26. Exemple III On demande la longueur de tout arc dont on connaît la tangente ?

L’arc dont la tangente est t est l’intégrale de d⁡x/(1+x²), prise depuis x=0 jusqu’à x=t. En se servant de la formule (II) les quantités qu’il faudra y substituer pour α, β, γ, …λ, μ, ν seront des fractions ayant pour numérateur commun 144, et pour dénominateurs respectifs les nombres 144, 144+t², 144+4t², 144+9t², 144+16t², 144+100t², 144+121t², 144+144t²; et l’arc cherché sera la valeur qui en résultera pour A, multipliée par t. L’exemple suivant nous fera juger du degré d’exactitude de ce procédé.

27. Exemple IV. On demande, suivant la formule précédente, la longueur de l’arc ϖ/4 ?

Les dénominateurs de nos quantités α, β, γ, …λ, μ, ν sont ici 144, 145, 148, 153, 160, 169, 180, 193, 208, 225, 244, 265, 288, ce qui donnera

α+ν=144/144+144/288 =3/2 =1, 5000000000, β+μ=144/145+144/265 =11808/5.29.53 =1, 5364996746, γ+λ=144/148+144/244 =3528/37.61 =1, 5631369074, δ+ϰ=144/153+144/225 =672/17.25 =1, 5811764706, ε+ι=144/160+144/208 =207/130 =1, 5923076923,

ζ+θ=144/169+144/193 =52128/169.193 =1, 5981849955, 2η=144/180+144/180 =8/5 =1, 6000000000,

d’où on conclura finalement

ϖ/4=0, 7853981727.

La valeur rigoureuse étant

ϖ/4=0, 7853981634.

il s’ensuit que l’erreur tombe au-dessous d’une unité décimale du huitième ordre.

28. Exemple V. Rectification générale de l’ellipse.

Ce célèbre problème qui exerce, depuis plus d’un siècle, le génie de nos plus grands analistcs, rentre de lui-même dans nos formules générales, dont il ne présente qu’un cas très-particulier. Soient a, b les deux demi-axes, et proposons-nous de rectifier l’arc compris depuis l’extrémité de a jusqu’au point dont la normale fait avec a un angle λ; la formule à intégrer sera

a²b²d⁡λ/((a²Cos.²λ+b²Sin.²λ)³ᐟ²),

depuis λ=0. Si l’on veut se contenter de la première formule, on remplacera successivement la lettre λ, dans

a²b²/((a²Cos.²λ+b²Sin.²λ)³ᐟ²),

d’abord par zéro et ensuite pars 1/6, 2/6, 3/6, 4/6, 5/6, 6/6 de λ. On aura ainsi les valeurs de α, β, γ, δ, ε, ζ, η, d’où on conclura celles de a′, b′, c′, d′, et par suite celle de A qui, multipliée par l’arc entier λ, exprimé en parties du rayon, fera connaître la longueur de l’arc cherché.

29. Supposons, par exemple, qu’il soit question d’assigner la longueur du quart de l’ellipse ; λ sera un angle droit, et devra conséquemment devenir successivement 0^{∘}, 15^{∘}, 30^{∘}, 45^{∘}, 60^{∘}, 75^{∘}, 90^{∘}. Or, on sait que

Sin.² 0^{∘}=Cos.²90^{∘}=0, Cos.² 0^{∘}=Sin.²90^{∘}=1, Sin.²15^{∘}=Cos.²75^{∘}=2−√3/4, Cos.²15^{∘}=Sin.²75^{∘}=(2+√3)/4, Sin.²30^{∘}=Cos.²60^{∘}=1/4, Cos.²30^{∘}=Sin.²60^{∘}=1/4, Sin.²45^{∘}=Cos.²90^{∘}=1/2.

Il viendra donc, en posant pour abréger le quarré a²−b² de l’excentricité =c²,

α=b²/a, β=a²b²/(a²−2−√3/4 c²)³ᐟ², γ=a²b²/(a²−1/4 c²)³ᐟ², δ=a²b²/(a²−1/4 c²)³ᐟ², ε=a²b²/(a²−3/4 c²)³ᐟ², ζ=a²b²/((a²−(2+√3)/4 c²)³ᐟ²), η=a²/b.

Ayant obtenu la valeur de A, on la multipliera par ϖ/2=1, 5708, et l’on aura ainsi une valeur du quart d’ellipse qui, dans les cas ordinaires, ne sera pas fautive d’un cent-millième.

30. Exemple VI. On demande l’intégrale de e⁻ᵗ²d⁡t, depuis t=0 jusqu’à t=3 ?

L’intervalle étant divisé en douze parties égales, on aura

Pour t=0, 00, la fonction e⁻ᵗ²=1, 0000000, 0, 250, 9394332, 0, 500, 7788009, 0, 750, 5697828, 1, 000, 3678794, 1, 250, 2096113, 1, 500, 1053992;1, 750, 0467706, 2, 000, 0183156, 2, 250, 0063297, 2, 500, 0013904, 2, 750, 0005244, 3, 000, 0001234.

En employant la formule (II), on trouve, pour le numérateur de l’intégrale ………………620, 3635233,

et pour son dénominateur …2100;

mais, à cause de l’intervalle 3, il faudra diviser par 700 seulement, ce qui donnera finalement

∫e⁻ᵗ²d⁡t=0, 8862336.

31. Dans mon Analise des réfractions astronomiques, j’ai donné une table des intégrales de ∫e⁻ᵗ²d⁡t prises jusqu’à l’infini. J’y trouve

t=0, ……0, 88622692,

t=2, ……0, 00001958;

0

3, ……0, 88620734;

la différence avec la précédente n’excède guère un quarante-millième d’unité.

32. Exemple VII. On demande l’intégrale de eˣd⁡x/x, depuis x=1 jusqu’à une valeur quelconque de x ?

La courbe dont l’équation est y=eˣ/x, n’a aucun de ses points situé dans les angles des coordonnées de signes contraires ; mais elle a, dans chacun des angles des coordonnées de mêmes signes, une partie qui présente deux branches infinies. Les axes sont les asymptotes de la partie située dans l’angle des coordonnées négatives ; quant à l’autre partie, elle n’a qu’une seule asymptote qui est l’axe des y; elle a une ordonnée minimum qui répond au point pour lequel on a x=1, et conséquemment y=e; et, à partir de ce point jusqu’à l’origine, l’accroissement de l’ordonnée est très-rapide, et s’élève à un ordre d’infini qu’il n’est pas même possible de déterminer ; de sorte que cette branche ne peut approcher indéfiniment d’aucune courbe connue, à moins peut-être que ce ne soit la branche curviligne de la Logistique ordinaire.

33. Cherchons, par la première formule, l’aire de la courbe, d’abord entre 1 et 7 puis entre 7 et 13; et nous chercherons ensuite, par la seconde formule, l’aire totale entre 1 et 13, laquelle doit être rigoureusement égale à la somme des deux premières. La différence que nous trouverons entre les deux résultats nous mettra, à même d’apprécier l’erreur que notre méthode laisse subsister, dans le cas particulier de ce problème.

34. On trouve, dans les Tables logarithmiques de Schulze (Berlin, 1778), une table des puissances de e=2, 71828…, depuis la première jusqu’à la vingt-quatrième. Divisant donc les treize premiers par leurs exposans respectifs, nous aurons nos treize ordonnées ainsi qu’il suit :

α=2, 718282, β=3, 694528, γ=6, 695178,

δ=29, 682632, ε=29, 682632, ζ=67, 238132, η=156, 661880, θ=372, 619748, ι=900, 342659, ϰ=2202, 646579, λ=5443, 103792, μ=13562, 899285, ν=34031, 769385.

35. On aura d’abord, par la formule (I), pour l’aire comprise entre 1 et 7,

(41(α+η)+216(β+ζ)+27(γ+ε)+272δ)/140,

et pour l’aire comprise entre 7 et 13,

(41(η+ν)+216(θ+μ)+27(ι+λ)+272ϰ)/140.

On trouvera ensuite, par la formule (II), pour l’aire totale, comprise entre 1 et 13,

7/25 (α+ν)+288/175 (β+ζ+θ+μ)−27/175 (γ+λ)+64/25 (δ+ϰ)

−9/25 (ε+ι)+134/175 η.

36. On aura ainsi

1

7, ……189, 649401,

7

13, ……37015, 696762,

1

13, ……37198, 443648.

Cette dernière est un peu moindre que la somme des deux autres ; et elle doit naturellement être réputée plus exacte ; la différence est 8, 902515 ; c’est environ la 4200ᵐᵉ partie de l’intégrale entière. Cette différence est un peu plus sensible que celles de tous les problèmes précédens ; mais il faut considérer aussi à quelle intégrale on avait à faire. Celle de eˣd⁡x/x tient une des premières places parmi ces intégrales éminemment réfractaires, qui se sont constamment jusqu’ici montrées rebelles à tous les moyens d’intégration connus, sans même en exclure l’emploi des séries infinies. Hinc ergo natura hujus functionis transcendentis parum cognoscitur, dit Euler (Calc. intég., vol. 1, n.ᵒ 228).

37. Égalant entre elles la somme des deux formules qui nous ont donné les aires partielles et celle qui nous a donné l’aire totale, on est conduit à cette nouvelle égalité très-remarquable

(α+ν)−8(β+μ)+27(γ+λ)−48(λ+ϰ)+43(ε+ι)−8(ζ+θ)=14η.

C’est l’équation de condition ; pour que le point du milieu, ou bien tout autre de nos treize points, se trouve sur la courbe déterminée par les douze autres. Elle est rigoureusement satisfaite, dans le cas où nos treize ordonnées sont égales entre elles : elle est rigoureusement remplie encore dans une infinité d’autres courbes dont il serait trop long de faire ici l’énumération. Elle est remplie, quoiqu’avec une différence presque insensible, lorsque la portion de courbe qui est comprise entre les limites de l’intégrale, est sans asymptote, sans imaginaires, sans point d’inflexion ni de rebroussement ; lorsqu’enfin elle ne s’écarte pas trop de quelque courbe rentrante ; telle que sont les ellipses de différens degrés. Ces sortes d’équation de condition, nouvelles dans l’analise, sont essentielles dans la théorie de l’interpolation ; elles pourront être le sujet d’un mémoire particulier.

38. Troisième formule. En conservant le diviseur général 12, ajoutons aux cinq premiers aliquotes 1, 2, 3, 4, 6, le nombre 12 lui-même. Nous aurons alors a=1, b=4, c=9, d=16, e=36, f=144; et ensuite (13)

21621600A=2985984a′−1667952b′+585728c′− 104247d′+2288e′−f′.

ce qui donnera

1801800A=41833(α+ν)+248832(β+ζ+θ+μ)−29160(γ+λ)+395264(δ+ϰ)−63909(ε+ι)+118416η.

Le calcul d’après cette formule est beaucoup plus compliqué, mais aussi elle réduit à peu près au quart l’erreur de l’autre.

39. Quatrième formule. En prenant pour diviseur général 18, et pour ses parties aliquotes 1,2,3,6,9, d’où a=1, b=4, c=9, d=36, e=81, on obtient d’abord (21)

277200A=24057a′−10935b′+2310c′−33d′+e′,

et ensuite

30800A=496(α+τ)+2673(β+ζ+θ+μ+o+σ)+243(γ+ε+ι+λ+ζ+ρ)+3443(δ+ϖ)991(η+ν)+3444ϰ,

40. Cinquième formule. En prenant pour diviseur général 24, et pour ses aliquotes 1, 2, 3, 4, 6, d’où a=1, b=4, c=9, d=16, e=36, on obtient d’abord

25200A=1728a′−945b′=320c′−54d′+e′.

Cette formule paraît être, au dénominateur près, qui est double, identique avec notre seconde formule ; elle ne l’est pourtant pas ; parce que le nombre des ordonnées étant double, les lettres a′, b′, c′, d′, e′, en acquièrent des valeurs entièrement différentes. On trouve, en effet, en désignant par ȣ la 25.ᵉ ordonnée

4200A=49(α+ȣ)+288(β+ζ+θ+μ+ξ+σ+υ+ω)−27(γ+λ+o+ψ)+448(δ+ϰ+ϖ+χ)−63(ε+ι+ϱ+ϕ)+134(η+τ)+98.

L’erreur que cette formule laisse subsister dans le calcul de l’arc ϖ/4, laquelle est soustractive, n’est sensible qu’à la onzième décimale.

41. Sixième formule. Le diviseur général étant encore 24 ; si l’on emploie les parties aliquotes 1, 2, 3, 4, 6, 8 ; ce qui donne a=1, b=4, c=9, d=16, e=36, f=64; on aura d’abord (13)

970200A=67584a′−38808b′+14336c′−2772d′+88e′−3f′;

et ensuite

80850A=933(α+ȣ)+5632(β+ζ+θ+μ+ξ+σ+υ+ω)−836(γ+λ+o+ψ)+9216(δ+ϰ+ϖ+χ)−1760(ε+ϕ)+2792(η+τ)−1762(ι+ϱ)+1868.

L’usage de cette formule réduit au quart l’erreur que la précédente avait laissé subsister.

42. Septième formule. Le diviseur général étant toujours 24, prenons ses aliquotes 1, 2, 3, 4, 6, 8, 12 ; ce qui donnera a=1, b=4, c=9, d=16, e=36, f=64, g=144; on aura d’abord

15765750A=26542080a′−15567552b′+5963776c′−1216215sd′+45760se′−2106f′+7g′;

et ensuite

315310500A=362135(α+ȣ)+2211840(β+ζ+θ+μ+ξ+σ+υ+ω)−382782(γ+λ+o+ψ)+3702784(δ+ϰ+ϖ+χ)−788157(ε+ϕ)+1131072(η+τ)−789561(ι+ϱ)+725674.

Par l’usage de cette formule l’erreur du calcul de ϖ/4, ne devient sensible que sur la douzième décimale.

43. Huitième formule. Prenons pour diviseur général 60, et pour ses aliquotes 1, 2, 3, 4, 5; d’où a=1, b=4, c=9, d=16, e=25; cela donnera d’abord

126A=210a′−120b′+45c′−10d′+e′.

Comme nous avons ici 61 ordonnées, nous représenterons

Les vingt premières par α, β, γ, …σ, τ, υ,

Les vingt suivantes par α′, β′, γ′, …σ′, τ′, υ′,

Les vingt autres par α″, β″, γ″, …σ″, τ″, υ″,

Et enfin la dernière par α‴. Il viendra ainsi

1512A=7(α+α‴)+42(β+θ+μ+ξ+σ+υ+δ′+ϰ′+μ′+σ′+β″+δ″+θ″+ϰ″+ξ″+υ″)−6(ε+o+γ′+η′+o′+τ′+η″+τ″)+69(δ+ϰ+β′+δ′+ξ′+υ′+μ″+σ″)−14(ε+ι+ϱ+ι′+ν′+ε″+ν″+ϱ″)+43(ζ+ζ′+ϖ′+ϖ″)+21(η+τ+γ″+o″)−5(λ+λ″)+13(υ+ε′+ϱ′+τ″)+70(ϖ+ζ″)−13(α′+α″)+110λ′.

Nous n’avons point mis cette formule à l’épreuve ; mais on peut présumer que dans le calcul de ϖ/4, l’erreur qu’entraînerait son usage tomberait au-delà de la vingtième décimale.

44. Nous ne craignons pas d’avancer qu’à l’aide de ces diverses formules, toute intégrale quelconque, de la forme ∫Xd⁡x, peut-être évaluée numériquement, entre les limites données, moyennant un nombre fini et très-limité de termes, indépendamment du calcul intégral, et de toute notion d’infiniment petit, par les seuls moyens que fournit l’algèbre élémentaire, et avec toute la précision que l’on veut donner à son calcul ; pourvu seulement que la fonction X ne devienne ni infinie ni imaginaire, dans l’étendue de l’intégration. Toute fonction intégrale d’une seule variable, telle que ∫Xd⁡x doit donc être comprise désormais dans la classe des quantités entièrement connues.

* ↑ Je dois l’idée, très-ingénieuse, qui sert de fondement à cette nouvelle méthode d’intégration à M. d’Obenheim, ancien sous-directeur des fortifications professeur de mathématiques à l’école d’artillerie de Strasbourg. Il l’a exposée dans un ouvrage qu’il vient de publier sous le titre de Balistique ou Indication de quelques expériences propres à compléter la théorie du mouvement des projectiles de l’artillerie ; mais je crois pouvoir en revendiquer les développemens et applications qui vont suivre, lesquels sont entièrement mon ouvrage. (Note de M. Kramp.)
