# 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 $X\operatorname {d} 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 $X\operatorname {d} 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_{0},a_{1},a_{2},$$\ldots a_{10},a_{11},a_{12}$ 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_{0},b_{1},b_{2},\ldots b_{10},b_{11},$$b_{12}\,;$ nous connaîtrons ainsi treize points de la courbe qu’il s’agit de quarrer entre les limites $x=a_{0}$ et $x=a_{12},$ pour lesquelles on a respectivement $y=b_{0},y=b_{12}.$
3. Soient joints les deux points extrêmes $\left(a_{0},b_{0}\right),\left(a_{12},b_{12}\right)$ 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_{12},$ et posant
$12\left({\tfrac {1}{2}}b_{0}+{\tfrac {1}{2}}b_{12}\right)=b'_{12},$
nous aurons
$S_{12}={\tfrac {1}{12}}cb'_{12}.$
4. Soient joints consécutivement les trois points $\left(a_{0},b_{0}\right),$ $\left(a_{6},b_{6}\right),$ $\left(a_{12},b_{12}\right)$ par deux cordes ; ces cordes formeront, avec $c$ et les trois ordonnées $b_{0},b_{6},b_{12},$ deux trapèzes ; en désignant la somme de leurs aires par $S_{6},$ et posant
$6\left({\tfrac {1}{2}}b_{0}+b_{6}+{\tfrac {1}{2}}b_{12}\right)=b'_{6},$
nous aurons
$S_{6}={\tfrac {1}{12}}cb'_{6}.$
5. Soient joints consécutivement les quatre points $\left(a_{0},b_{0}\right),\left(a_{4},b_{4}\right),$$\left(a_{8},b_{8}\right),\left(a_{12},b_{12}\right)$ trois cordes ; ces cordes formeront, avec $c$ et les quatre ordonnées $b_{0},b_{4},b_{8},b_{12},$ trois trapèzes ; en désignant la somme de leurs aires par $S_{4},$ et posant
$4\left({\tfrac {1}{2}}b_{0}+b_{4}+b_{8}+{\tfrac {1}{2}}b_{12}\right)=b'_{4},$
nous aurons
$S_{4}={\tfrac {1}{12}}cb'_{4}.$
6. Soient joints consécutivement les cinq points $\left(a_{0},b_{0}\right),\left(a_{3},b_{3}\right),$$\left(a_{6},b_{6}\right),\left(a_{9},b_{9}\right),\left(a_{12},b_{12}\right)$ par quatre cordes ; ces cordes formeront, avec $c$ et les cinq ordonnées $b_{0},b_{3},b_{6},b_{9},b_{12},$ quatre trapèzes ; en désignant la somme de leurs aires par $S_{3},$ et posant
$3\left({\tfrac {1}{2}}b_{0}+b_{3}+b_{6}+b_{9}+{\tfrac {1}{2}}b_{12}\right)=b'_{3},$
nous aurons
$S_{3}={\tfrac {1}{12}}cb'_{3}.$
7. Soient joints consécutivement les sept points $\left(a_{0},b_{0}\right),\left(a_{2},b_{2}\right),$$\left(a_{4},b_{4}\right),$$\left(a_{6},b_{6}\right),\left(a_{8},b_{8}\right),\left(a_{10},b_{10}\right)\left(a_{12},b_{12}\right)$ par six cordes ; ces cordes formeront, avec $c$ et les sept ordonnées $b_{0},b_{2},b_{4},b_{6},b_{8},b_{10},b_{12},$ six trapèzes ; en désignant la somme de leurs aires par $S_{2},$ et posant
$2\left({\tfrac {1}{2}}b_{0}+b_{2}+b_{4}+b_{6}+b_{8}+b_{10}+{\tfrac {1}{2}}b_{12}\right)=b'_{2},$
nous aurons
$S_{2}={\tfrac {1}{12}}cb'_{2}.$
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_{1}$, et posant
${\tfrac {1}{2}}b_{0}+b_{1}+b_{2}+b_{3}+b_{4}+b_{5}+b_{6}+b_{7}+b_{8}+b_{9}+b_{10}+b_{11}+{\tfrac {1}{2}}b_{12}=b'_{1},$
nous aurons
$S_{1}={\tfrac {1}{12}}cb'_{1}.$
9. Aucune des aires $S_{12},S_{6},S_{4},S_{3},S_{2},S_{1},$ 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 ${\tfrac {1}{12}}c$ pour avoir l’aire demandée n’est aucune des lignes $b'_{12},b'_{6},b'_{4},b'_{3},b'_{2},b'_{1},$ mais une ligne vers laquelle celles-là convergent de plus en plus.
10. Portons les ordonnées $b'_{12},b'_{6},b'_{4},b'_{3},b'_{2},b'_{1},$ sur leurs correspondantes $b_{12},b_{6},b_{4},b_{3},b_{2},b_{1}\,;$ et imaginons une courbe située au-dessus de la première, passant par les six points $\left(a_{1},b'_{1}\right),\left(a_{2},b'_{2}\right),$$\left(a_{3},b'_{3}\right),\left(a_{4},b'_{4}\right),$$\left(a_{6},b'_{6}\right),\left(a_{12},b'_{12}\right)\,;$ cette courbe prolongée rencontrera le prolongement de l’ordonnée $b_{0}$ en quelque point, en désignant par $b'_{0}$ son ordonnée qui répond à celle-là, et conséquemment à l’abscisse $a_{0}$ les ordonnées $b'_{12},b'_{6},b'_{4}$$,b'_{3},b'_{2},b'_{1},$ tendant continuellement vers la ligne par laquelle il faut multiplier ${\tfrac {1}{12}}c$ pour avoir l’aire cherchée, en désignant cette aire par $S_{0},$ nous pourrons prendre sensiblement
$S_{0}={\tfrac {1}{12}}cb'_{0}.$
et tout se réduira à trouver $b'_{0}$ ; 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'_{0}$ 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'_{12},b'_{6},b'_{4},b'_{3},b'_{2},b'_{1},$ qu’elle doit couper perpendiculairement l’ordonnée $b'_{0},$ 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'_{0},$ on doit avoir en même temps $x=0,{\frac {\operatorname {d} y'}{\operatorname {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^{2}+Cx^{4}+Dx^{6}+Ex^{8}+\ldots \,;$
il s’agira donc de déterminer le coefficient $A,$ auquel se réduit $y'$ lorsque $x=0.$
12. En prenant ${\tfrac {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'_{1},b'_{2},b'_{3},b'_{4},b'_{6},b'_{12},$ ce qui donnera
${\begin{alignedat}{1}b'_{1}\ =&A+\quad \ \,B+\quad \ \,C+\quad \ \ D+\quad \ \,E+\quad \ \ \ F,\\b'_{2}\ =&A+\ \ 2^{2}B+\ \ 2^{4}C+\ \ 2^{6}D+\ \ 2^{8}E+\ \ 2^{10}F,\\b'_{3}\ =&A+\ \ 3^{2}B+\ \ 3^{4}C+\ \ 3^{6}D+\ \ 3^{8}E+\ \ 3^{10}F,\\b'_{4}\ =&A+\ \ 4^{2}B+\ \ 4^{4}C+\ \ 4^{6}D+\ \ 4^{8}E+\ \ 4^{10}F,\\b'_{6}\ =&A+\ \ 6^{2}B+\ \ 6^{4}C+\ \ 6^{6}D+\ \ 6^{8}E+\ \ 6^{10}F,\\b'_{12}\ =&A+12^{2}B+12^{4}C+12^{6}D+12^{8}E+12^{10}F\,;\end{alignedat}}$
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'_{1},b'_{2},b'_{3},b'_{4},b'_{6},b'_{12},$ sera, pour $c=12,$ l’intégrale demandée ; nous avons vu d’ailleurs qu’on a
${\begin{aligned}b'_{1}&={\tfrac {1}{2}}b_{0}+\ \ b_{1}+\ \ b_{2}+\ \ b_{3}+b_{4}+b_{5}+b_{6}+b_{7}+b_{8}+b_{9}+b_{10}+b_{11}+{\tfrac {1}{2}}b_{12},\\b'_{2}&=\ \ \ b_{0}+\ \ b_{2}+\ \ b_{4}+\ \ b_{6}+b_{8}+b_{10}+b_{12},\\b'_{3}&={\tfrac {1}{2}}b_{0}+3b_{3}+3b_{6}+3b_{9}+{\tfrac {1}{2}}b_{12},\\b'_{4}&=\ 2b_{0}+4b_{4}+4b_{8}+2b_{12},\\b'_{6}&=\ 3b_{0}+6b_{6}+3b_{12},\\\end{aligned}}$
$b'_{12}=6b_{0}+6b_{12}.\qquad$
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
${\begin{aligned}a'&=A+Ba+Ca^{2}+Da^{3}+Ea^{4}+Fa^{5},\\b'&=A+Bb+Cb^{2}+Db^{3}+Eb^{4}+Fb^{5},\\c'&=A+Bc+Cc^{2}+Dc^{3}+Ec^{4}+Fc^{5},\\d'&=A+Bd+Cd^{2}+Dd^{3}+Ed^{4}+Fd^{5},\\e'&=A+Be+Ce^{2}+De^{3}+Ee^{4}+Fe^{5},\\f'&=A+Bf+Cf^{2}+Df^{3}+Ef^{4}+Ff^{5}.\\\end{aligned}}$
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
${\begin{aligned}A&={\frac {bcdefa'}{(b-a)(c-a)(d-a)(e-a)(f-a)}}\\&+{\frac {cdefab'}{(c-b)(d-b)(e-b)(f-b)(a-b)}}\\&+{\frac {defabc'}{(d-c)(e-c)(f-c)(a-c)(b-c)}}\\&+{\frac {efabcd'}{(e-d)(f-d)(a-d)(b-d)(c-d)}}\\&+{\frac {fabcde'}{(f-e)(a-e)(b-e)(c-e)(d-e)}}\\&+{\frac {abcdef'}{(a-f)(b-f)(c-f)(d-f)(e-f)}}\\\end{aligned}}$
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,\ldots .$ 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={\tfrac {bcda'}{(b-a)(c-a)(d-a)}}+{\tfrac {cdab'}{(c-b)(d-b)(e-b)}}+{\tfrac {dabc'}{(d-c)(a-c)(b-c)}}+{\tfrac {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={\frac {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 $\alpha ,\beta ,\gamma ,\delta ,\varepsilon ,\zeta ,\eta \,;$ nous aurons
${\begin{aligned}a'&={\tfrac {1}{2}}\alpha +\beta +\gamma +\delta +\varepsilon +\zeta +{\tfrac {1}{2}}\eta ,\\b'&=\alpha +2\gamma +2\varepsilon +\eta ,\\c'&={\tfrac {1}{2}}\alpha +3\delta +{\tfrac {1}{2}}\eta ,\\d'&=3\alpha +3\eta \,;\\\end{aligned}}$
qui donnera ; en substituant,
$A={\frac {246(\alpha +\eta )+1296(\beta +\zeta )+162(\gamma +\varepsilon )+1632\delta }{840}}.$
En prenant donc l’intervalle entier qui sépare les ordonnées extrêmes $\alpha$ et $\eta$ pour unité ; on aura finalement pour l’intégrale cherchée
$\int X\operatorname {d} x={\frac {41(\alpha +\eta )+216(\beta +\zeta )+27(\gamma +\varepsilon )+272\delta }{840}}.\qquad$
Si, dans cette dernière formule, on fait $\alpha =\beta =\gamma =\delta =\varepsilon =\zeta =\eta =1,$ elle devient $\int X\operatorname {d} 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 ${\frac {\operatorname {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={\frac {1}{x}}$ nous aurons d’abord
${\begin{array}{llll}\alpha &={\tfrac {6}{6}}=1,00000000,&\alpha +\eta &=1,50000000,\\\beta &={\tfrac {6}{7}}=0.85714286,&\beta +\zeta &=1,40259740,\\\gamma &={\tfrac {6}{8}}=0,75000000,&\gamma +\varepsilon &=1,35000000,\\\delta &={\tfrac {6}{9}}=0,66666667,&2\delta &=1,33333333.\\\end{array}}$
${\begin{array}{ll}\varepsilon &={\tfrac {6}{10}}=0,60000000,\\\zeta &={\tfrac {6}{11}}=0,54545455,\\\eta &={\tfrac {6}{12}}=0,50000000\,;\\\end{array}}$
Ces valeurs étant substituées dans la formule (I), on aura
$\operatorname {Log} .2=0,69314806.$
La valeur rigoureuse est
$\operatorname {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 $={\frac {\varpi }{4}}$ ?
La longueur de l’arc dont la tangente est $t$ est l’intégrale de ${\frac {\operatorname {d} x}{1+x^{2}}},$ prise depuis $x=0$ jusqu’à $x=t$ ; celle de l’arc $={\frac {\varpi }{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={\tfrac {1}{1+x^{2}}},$ il viendra
${\begin{array}{llll}\alpha &={\tfrac {36}{36}}=1,00000000,&\alpha +\eta &=1,50000000,\\\beta &={\tfrac {36}{37}}=0,97297297,&\beta +\zeta &mm=1,56313690,\\\gamma &={\tfrac {36}{40}}=0,90000000,&\gamma +\varepsilon &=1,59230769,\\\delta &={\tfrac {36}{45}}=0,80000000,&2\delta &=1,60000000.\\\varepsilon &={\tfrac {36}{52}}=1,69230769,&&\\\zeta &={\tfrac {36}{62}}=0,59016393,&&\\\eta &={\tfrac {36}{72}}=0,50000000\,;&&\\\end{array}}$
Ces valeurs étant substituées dans la formule (I), on aura
$={\frac {\varpi }{4}}=0,78539271.$
La valeur rigoureuse est
${\frac {\varpi }{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
${\begin{aligned}A&={\tfrac {bcdea'}{(b-a)(c-a)(d-a)(e-a)}}+{\tfrac {cdeab'}{(c-b)(d-b)(e-b)(a-b)}}+{\tfrac {deabc'}{(d-c)(e-c)(a-c)(b-c)}}\\&\qquad +{\tfrac {eabcd'}{(e-d)(a-d)(b-d)(c-d)}}+{\tfrac {abcde'}{(a-e)(b-e)(c-e)(d-e)}}.\\\end{aligned}}$
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={\frac {1728a'-945b'+320c'-54d'+e'}{12600}}\,;$
mais on a, dans le cas actuel,
${\begin{aligned}a'&={\tfrac {1}{2}}\alpha +\beta +\gamma +\delta +\varepsilon +\zeta +\eta +\theta +\iota +\varkappa +\lambda +\mu +{\tfrac {1}{2}}\nu ,\\b'&=\alpha +2\gamma +2\varepsilon +2\eta +2\iota +2\lambda +\nu ,\\c'&={\tfrac {1}{2}}\alpha +3\delta +3\eta +2\varkappa +{\tfrac {1}{2}}\nu ,\\d'&=2\alpha +4\eta +4\iota +2\nu \,;\\e'&=3\alpha +6\eta +3\nu \,;\\\end{aligned}}$
qui donnera, en substituant,
$\int X\operatorname {d} x={\tfrac {49(\alpha +\nu )+288(\beta +\mu )-27(\gamma +\lambda )+448(\delta +\varkappa )-63(\varepsilon +\iota )+288(\zeta +\theta )+134\eta }{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
${\begin{array}{llll}\alpha &={\tfrac {12}{12}}=1,0000000000,&\alpha +\nu &=1,5000000000,\\\beta &={\tfrac {12}{13}}=0,9230769231,&\beta +\mu &=1,4448160535,\\\gamma &={\tfrac {12}{14}}=0,8571428571,&\gamma +\lambda &=1,4025974026,\\\delta &={\tfrac {12}{15}}=0,8000000000,&\delta +\varkappa &=1,3714285714,\\\varepsilon &={\tfrac {12}{16}}=0,7500000000,&\varepsilon +\iota &=1,3500000000,\\\zeta &={\tfrac {12}{17}}=0,7058823529,&\zeta +\theta &=1,3374613003,\\\eta &={\tfrac {12}{18}}=0,6666666667,&2\eta &=1,3333333333\,;\\\theta &={\tfrac {12}{19}}=1,6315789474,&&\\\iota &={\tfrac {12}{20}}=0,6000000000,&&\\\varkappa &={\tfrac {12}{21}}=0,5714285714,&&\\\lambda &={\tfrac {12}{22}}=0,5454545455,&&\\\mu &={\tfrac {12}{23}}=0,5217391304,&&\\\nu &={\tfrac {12}{24}}=0,5000000000\,;&&\\\end{array}}$
ce qui donnera, en substituant dans la formule (II)
$\operatorname {Log} .2=0.6931471816\,;$
la valeur rigoureuse est
$\operatorname {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^{n}$ 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 $\alpha ={\frac {6}{6}},\beta ={\frac {6}{6+h}},$ $\gamma ={\frac {6}{6+2h}},\delta ={\frac {6}{6+3h}},\varepsilon ={\frac {6}{6+4h}},\zeta ={\frac {6}{6+5h}},\eta ={\frac {6}{6+6h}}={\frac {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^{13}$. On aura ainsi $m=10000,n=13,1+h={\frac {10000}{8192}}=1+{\frac {1808}{8192}}$$=1+{\frac {113}{512}}\,;$ on aura donc
$\operatorname {Log} .10000=13\operatorname {Log} .2+\operatorname {Log} (1+{\frac {113}{512}}).$
Pour trouver ce dernier logarithme, on fera $\alpha ={\frac {3072}{3072}},\beta ={\frac {3072}{3185}},$
$\gamma ={\frac {3072}{3298}},\delta ={\frac {3072}{3411}},\varepsilon ={\frac {3072}{3524}},\zeta ={\frac {3072}{3637}},\eta ={\frac {3072}{3350}}$
On trouvera ensuite
${\begin{aligned}41(\alpha +\eta )&=\ \ 74,58720000,\\216(\beta +\zeta )&=390,78145008,\\27(\gamma +\varepsilon )&=\ \ 48,68667756,\\272\ \delta &=244,96745680.\\\end{aligned}}$
La somme $759,02278444$ de ces quatre nombres, divisée par $840,$ et multipliée par $h={\frac {113}{512}}$ donne pour le logarithme de $1+h$ ou $h={\frac {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é $\operatorname {Log} .2=0,6931471816,$ on aura $13\operatorname {Log} .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
$\operatorname {Log} .10=2,3025850963.$
La valeur rigoureuse étant
$\operatorname {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 ${\frac {\operatorname {d} x}{1+x^{2}}},$ prise depuis $x=0$ jusqu’à $x=t.$ En se servant de la formule (II) les quantités qu’il faudra y substituer pour $\alpha ,\beta ,\gamma ,\ldots \lambda ,\mu ,\nu$ seront des fractions ayant pour numérateur commun $144,$ et pour dénominateurs respectifs les nombres $144,$$144+t^{2},144+4t^{2},144+9t^{2},144+16t^{2},$ $144+100t^{2},144+121t^{2},$$144+144t^{2}\,;$ 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 ${\tfrac {\varpi }{4}}$ ?
Les dénominateurs de nos quantités $\alpha ,\beta ,\gamma ,\ldots \lambda ,\mu ,\nu$ sont ici $144,145,148,153,160,169,180,193,208,225,244,265,288,$ ce qui donnera
${\begin{array}{llll}\alpha +\nu &={\frac {144}{144}}+{\frac {144}{288}}&={\frac {3}{2}}&=1,5000000000,\\\beta +\mu &={\frac {144}{145}}+{\frac {144}{265}}&={\frac {11808}{5.29.53}}&=1,5364996746,\\\gamma +\lambda &={\frac {144}{148}}+{\frac {144}{244}}&={\frac {3528}{37.61}}&=1,5631369074,\\\delta +\varkappa &={\frac {144}{153}}+{\frac {144}{225}}&={\frac {672}{17.25}}&=1,5811764706,\\\varepsilon +\iota &={\frac {144}{160}}+{\frac {144}{208}}&={\frac {207}{130}}&=1,5923076923,\\\end{array}}$
${\begin{array}{llll}\zeta +\theta &={\frac {144}{169}}+{\frac {144}{193}}&={\frac {52128}{169.193}}&=1,5981849955,\\2\eta &={\frac {144}{180}}+{\frac {144}{180}}&={\frac {8}{5}}&=1,6000000000,\\\end{array}}$
d’où on conclura finalement
${\frac {\varpi }{4}}=0,7853981727.$
La valeur rigoureuse étant
${\frac {\varpi }{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 $\lambda \,;$ la formule à intégrer sera
${\frac {a^{2}b^{2}\operatorname {d} \lambda }{\left(a^{2}\operatorname {Cos} .^{2}\lambda +b^{2}\operatorname {Sin} .^{2}\lambda \right)^{\tfrac {3}{2}}}},$
depuis $\lambda =0.$ Si l’on veut se contenter de la première formule, on remplacera successivement la lettre $\lambda ,$ dans
${\frac {a^{2}b^{2}}{\left(a^{2}\operatorname {Cos} .^{2}\lambda +b^{2}\operatorname {Sin} .^{2}\lambda \right)^{\tfrac {3}{2}}}},$
d’abord par zéro et ensuite pars ${\tfrac {1}{6}},{\tfrac {2}{6}},{\tfrac {3}{6}},{\tfrac {4}{6}},{\tfrac {5}{6}},{\tfrac {6}{6}}$ de $\lambda .$ On aura ainsi les valeurs de $\alpha ,\beta ,\gamma ,\delta ,\varepsilon ,\zeta ,\eta ,$ d’où on conclura celles de $a',b',c',d',$ et par suite celle de $A$ qui, multipliée par l’arc entier $\lambda ,$ 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 ; $\lambda$ sera un angle droit, et devra conséquemment devenir successivement $0^{\circ },15^{\circ },30^{\circ },45^{\circ },60^{\circ },75^{\circ },90^{\circ }.$ Or, on sait que
${\begin{array}{lll}\operatorname {Sin} .^{2}\ \ 0^{\circ }&=\operatorname {Cos} .^{2}90^{\circ }&=0,\\\operatorname {Cos} .^{2}\ 0^{\circ }&=\operatorname {Sin} .^{2}90^{\circ }&=1,\\\operatorname {Sin} .^{2}15^{\circ }&=\operatorname {Cos} .^{2}75^{\circ }&={\frac {2-{\sqrt {3}}}{4}},\\\operatorname {Cos} .^{2}15^{\circ }&=\operatorname {Sin} .^{2}75^{\circ }&={\frac {2+{\sqrt {3}}}{4}},\\\operatorname {Sin} .^{2}30^{\circ }&=\operatorname {Cos} .^{2}60^{\circ }&={\frac {1}{4}},\\\operatorname {Cos} .^{2}30^{\circ }&=\operatorname {Sin} .^{2}60^{\circ }&={\frac {1}{4}},\\\operatorname {Sin} .^{2}45^{\circ }&=\operatorname {Cos} .^{2}90^{\circ }&={\frac {1}{2}}.\\\end{array}}$
Il viendra donc, en posant pour abréger le quarré $a^{2}-b^{2}$ de l’excentricité $=c^{2},$
${\begin{aligned}\alpha &={\frac {b^{2}}{a}},\\\beta &={\frac {a^{2}b^{2}}{(a^{2}-{\frac {2-{\sqrt {3}}}{4}}c^{2})^{\frac {3}{2}}}},\\\gamma &={\frac {a^{2}b^{2}}{(a^{2}-{\frac {1}{4}}c^{2})^{\frac {3}{2}}}},\\\delta &={\frac {a^{2}b^{2}}{(a^{2}-{\frac {1}{4}}c^{2})^{\frac {3}{2}}}},\\\varepsilon &={\frac {a^{2}b^{2}}{(a^{2}-{\frac {3}{4}}c^{2})^{\frac {3}{2}}}},\\\zeta &={\frac {a^{2}b^{2}}{(a^{2}-{\frac {2+{\sqrt {3}}}{4}}c^{2})^{\frac {3}{2}}}},\\\eta &={\frac {a^{2}}{b}}.\\\end{aligned}}$
Ayant obtenu la valeur de $A,$ on la multipliera par ${\frac {\varpi }{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^{-t^{2}}\operatorname {d} t,$ depuis $t=0$ jusqu’à $t=3$ ?
L’intervalle étant divisé en douze parties égales, on aura
${\begin{array}{rr}{\text{Pour }}t=0{,}00,&{\text{la fonction }}e^{-t^{2}}=1{,}0000000,\\0{,}25&0{,}9394332,\\0{,}50&0{,}7788009,\\0{,}75&0{,}5697828,\\1{,}00&0{,}3678794,\\1{,}25&0{,}2096113,\\1{,}50&0{,}1053992;\\1{,}75&0{,}0467706,\\2{,}00&0{,}0183156,\\2{,}25&0{,}0063297,\\2{,}50&0{,}0013904,\\2{,}75&0{,}0005244,\\3{,}00&0{,}0001234.\\\end{array}}$
En employant la formule (II), on trouve, pour le numérateur de l’intégrale $\ldots \ldots \ldots \ldots \ldots \ldots 620{,}3635233,$
et pour son dénominateur $\ldots 2100\,;$
mais, à cause de l’intervalle $3,$ il faudra diviser par $700$ seulement, ce qui donnera finalement
$\int e^{-t^{2}}\operatorname {d} t=0{,}8862336.$
31. Dans mon Analise des réfractions astronomiques, j’ai donné une table des intégrales de $\int e^{-t^{2}}\operatorname {d} t$ prises jusqu’à l’infini. J’y trouve
$t=0,\ldots \ldots 0{,}88622692,$
$t=2,\ldots \ldots 0{,}00001958\,;$
$0\,$
$\,3,\ldots \ldots 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 ${\frac {e^{x}\operatorname {d} x}{x}},$ depuis $x=1$ jusqu’à une valeur quelconque de $x$ ?
La courbe dont l’équation est $y={\frac {e^{x}}{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\ldots ,$ 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 :
${\begin{array}{lr}\alpha =&\qquad 2{,}718282,\\\beta =&\qquad 3{,}694528,\\\gamma =&\qquad 6{,}695178,\\\end{array}}$
${\begin{array}{lr}\delta =&29{,}682632,\\\varepsilon =&29{,}682632,\\\zeta =&67{,}238132,\\\eta =&156{,}661880,\\\theta =&372{,}619748,\\\iota =&900{,}342659,\\\varkappa =&2202{,}646579,\\\lambda =&5443{,}103792,\\\mu =&13562{,}899285,\\\nu =&34031{,}769385.\end{array}}$
35. On aura d’abord, par la formule (I), pour l’aire comprise entre $1$ et $7,$
${\frac {41(\alpha +\eta )+216(\beta +\zeta )+27(\gamma +\varepsilon )+272\delta }{140}},$
et pour l’aire comprise entre $7$ et $13,$
${\frac {41(\eta +\nu )+216(\theta +\mu )+27(\iota +\lambda )+272\varkappa }{140}}.$
On trouvera ensuite, par la formule (II), pour l’aire totale, comprise entre $1$ et $13,$
${\frac {7}{25}}(\alpha +\nu )+{\frac {288}{175}}(\beta +\zeta +\theta +\mu )-{\frac {27}{175}}(\gamma +\lambda )+{\frac {64}{25}}(\delta +\varkappa )$
$-{\frac {9}{25}}(\varepsilon +\iota )+{\frac {134}{175}}\eta .$
36. On aura ainsi
$1$
$\ 7,\ldots \ldots \quad 189{,}649401,$
$7$
$13,\ldots \ldots 37015{,}696762,$
$1$
$13,\ldots \ldots 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^{me}$ 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 ${\frac {e^{x}\operatorname {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.$^o$ 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
$(\alpha +\nu )-8(\beta +\mu )+27(\gamma +\lambda )-48(\lambda +\varkappa )+43(\varepsilon +\iota )-8(\zeta +\theta )=14\eta .$
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)
${\begin{array}{rll}21621600A&=2985984a'&-1667952b'\\&+\quad 585728c'&-\ \ 104247d'\\&+\qquad 2288e'&-\qquad \quad f'.\\\end{array}}$
ce qui donnera
${\begin{aligned}1801800A&=41833(\alpha +\nu )+248832(\beta +\zeta +\theta +\mu )\\&-29160(\gamma +\lambda )+395264(\delta +\varkappa )\\&-63909(\varepsilon +\iota )+118416\eta .\\\end{aligned}}$
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
${\begin{aligned}30800A&=496(\alpha +\tau )+2673(\beta +\zeta +\theta +\mu +o+\sigma )\\&+243(\gamma +\varepsilon +\iota +\lambda +\zeta +\rho )+3443(\delta +\varpi )\\&\qquad 991(\eta +\nu )+3444\varkappa ,\\\end{aligned}}$
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.e ordonnée
${\begin{aligned}4200A&=49(\alpha +{\text{ȣ}})+288(\beta +\zeta +\theta +\mu +\xi +\sigma +\upsilon +\omega )\\&-27(\gamma +\lambda +o+\psi )+448(\delta +\varkappa +\varpi +\chi )\\&\qquad -63(\varepsilon +\iota +\varrho +\phi )+134(\eta +\tau )+98.\\\end{aligned}}$
L’erreur que cette formule laisse subsister dans le calcul de l’arc ${\frac {\varpi }{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
${\begin{aligned}80850A&=933(\alpha +{\text{ȣ}})+5632(\beta +\zeta +\theta +\mu +\xi +\sigma +\upsilon +\omega )\\&-836(\gamma +\lambda +o+\psi )+9216(\delta +\varkappa +\varpi +\chi )\\&-1760(\varepsilon +\phi )+2792(\eta +\tau )-1762(\iota +\varrho )+1868.\\\end{aligned}}$
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
${\begin{aligned}15765750A&=26542080a'-15567552b'+5963776c'\\&-1216215sd'+45760se'-2106f'+7g'\,;\\\end{aligned}}$
et ensuite
${\begin{aligned}315310500A&=362135(\alpha +{\text{ȣ}})+2211840(\beta +\zeta +\theta +\mu +\xi +\sigma +\upsilon +\omega )\\&-382782(\gamma +\lambda +o+\psi )+3702784(\delta +\varkappa +\varpi +\chi )\\&-788157(\varepsilon +\phi )+1131072(\eta +\tau )-789561(\iota +\varrho )+725674.\\\end{aligned}}$
Par l’usage de cette formule l’erreur du calcul de ${\frac {\varpi }{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 $\alpha ,\beta ,\gamma ,\ldots \sigma ,\tau ,\upsilon ,$
Les vingt suivantes par $\ \alpha ',\beta ',\gamma ',\ldots \sigma ',\tau ',\upsilon ',$
Les vingt autres par $\quad \ \ \alpha '',\beta '',\gamma '',\ldots \sigma '',\tau '',\upsilon '',$
Et enfin la dernière par $\ \ \alpha '''.$ Il viendra ainsi
${\begin{aligned}1512A&=7(\alpha +\alpha ''')+42(\beta +\theta +\mu +\xi +\sigma +\upsilon +\delta '+\varkappa '+\mu '+\sigma '\\&+\beta ''+\delta ''+\theta ''+\varkappa ''+\xi ''+\upsilon '')-6(\varepsilon +o+\gamma '+\eta '+o'+\tau '+\eta ''+\tau '')\\&+69(\delta +\varkappa +\beta '+\delta '+\xi '+\upsilon '+\mu ''+\sigma '')-14(\varepsilon +\iota +\varrho +\iota '+\nu '+\varepsilon ''+\nu ''+\varrho '')\\&+43(\zeta +\zeta '+\varpi '+\varpi '')+21(\eta +\tau +\gamma ''+o'')-5(\lambda +\lambda '')\\&+13(\upsilon +\varepsilon '+\varrho '+\tau '')+70(\varpi +\zeta '')-13(\alpha '+\alpha '')+110\lambda '.\\\end{aligned}}$
Nous n’avons point mis cette formule à l’épreuve ; mais on peut présumer que dans le calcul de ${\frac {\varpi }{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 $\int X\operatorname {d} 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 $\int X\operatorname {d} 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.)