Introduction : comparaison entre série et intégrale
modifier
Soit f une fonction infiniment dérivable sur [1, +∞[ et n un entier naturel non nul.
On veut obtenir un développement asymptotique de la somme
∑
i
=
1
n
f
(
i
)
=
f
(
1
)
+
f
(
2
)
+
⋯
+
f
(
n
)
{\displaystyle \sum _{i=1}^{n}f(i)=f(1)+f(2)+\cdots +f(n)}
en la comparant à l'intégrale
∫
1
n
f
(
x
)
d
x
{\displaystyle \int _{1}^{n}f(x)~{\rm {d}}x}
.
La formule d'Euler-Maclaurin donne une expression de la différence
f
(
1
)
+
f
(
2
)
+
⋯
+
f
(
n
)
−
∫
1
n
f
(
x
)
d
x
{\displaystyle f(1)+f(2)+\cdots +f(n)-\int _{1}^{n}f(x)~{\rm {d}}x}
en fonction des valeurs de la fonction et de ses dérivées aux extrémités 1 et n et d'un reste :
∑
i
=
1
n
f
(
i
)
−
∫
1
n
f
(
x
)
d
x
=
f
(
1
)
+
f
(
n
)
2
+
1
6
f
′
(
n
)
−
f
′
(
1
)
2
!
−
1
30
f
‴
(
n
)
−
f
‴
(
1
)
4
!
+
⋯
+
b
2
k
f
(
2
k
−
1
)
(
n
)
−
f
(
2
k
−
1
)
(
1
)
(
2
k
)
!
+
R
k
,
n
.
{\displaystyle \sum _{i=1}^{n}f(i)-\int _{1}^{n}f(x)~{\rm {d}}x={\frac {f(1)+f(n)}{2}}+{\frac {1}{6}}{\frac {f'(n)-f'(1)}{2!}}-{\frac {1}{30}}{\frac {f'''(n)-f'''(1)}{4!}}+\cdots +b_{2k}{\frac {f^{(2k-1)}(n)-f^{(2k-1)}(1)}{(2k)!}}+R_{k,n}.}
Les nombres
b
2
=
1
6
,
b
4
=
−
1
30
,
b
6
=
1
42
,
b
8
=
−
1
30
,
…
b
2
k
{\displaystyle b_{2}={\frac {1}{6}},\quad b_{4}=-{\frac {1}{30}},\quad b_{6}={\frac {1}{42}},\quad b_{8}=-{\frac {1}{30}},\quad \ldots b_{2k}}
qui apparaissent dans la formule sont les nombres de Bernoulli .
La série obtenue n'est en général pas convergente mais on connaît plusieurs expressions du reste R k ,n de la formule qui permettent de majorer l'erreur ainsi faite.
Soient p et q deux entiers relatifs (p < q ), f une fonction continue complexe définie sur [p , q ].
L'énoncé qui suit exprime la somme
f
(
p
)
+
f
(
p
+
1
)
+
⋯
+
f
(
q
)
=
∑
i
=
p
q
f
(
i
)
{\displaystyle {f(p)}+f(p+1)+\cdots +{f(q)}=\sum _{i=p}^{q}f\left(i\right)}
avec l'intégrale
∫
p
q
f
(
x
)
d
x
{\displaystyle \int _{p}^{q}f(x)~{\rm {d}}x}
, les valeurs de f (ainsi que de ses dérivées) aux extrémités f (p ) et f (q ) et d'un reste.
Si f est une fonction complexe continûment dérivable une fois sur le segment [p , q ], la formule d'Euler-Maclaurin s'énonce ainsi :
∑
i
=
p
q
f
(
i
)
=
f
(
p
)
+
f
(
q
)
2
+
∫
p
q
f
(
x
)
d
x
+
R
0
{\displaystyle \sum _{i=p}^{q}f\left(i\right)={\frac {f(p)+f(q)}{2}}+\int _{p}^{q}f(x)~{\rm {d}}x+R_{0}}
avec
R
0
=
∫
p
q
f
′
(
x
)
(
x
−
⌊
x
⌋
−
1
2
)
d
x
{\displaystyle R_{0}=\int _{p}^{q}f'(x)\left(x-\lfloor x\rfloor -{\frac {1}{2}}\right)~{\rm {d}}x}
où
⌊
x
⌋
{\displaystyle \lfloor x\rfloor }
est la partie entière de x , notée aussi E(x ), et
x
−
⌊
x
⌋
{\displaystyle x-\lfloor x\rfloor }
est la partie fractionnaire de x .
Pour une fonction f continûment dérivable 2k fois sur le segment [p , q ] (avec k ≥ 1), la formule d'Euler-Maclaurin s'énonce ainsi :
∑
i
=
p
q
f
(
i
)
=
f
(
p
)
+
f
(
q
)
2
+
∫
p
q
f
(
x
)
d
x
+
∑
j
=
1
k
b
2
j
(
2
j
)
!
(
f
(
2
j
−
1
)
(
q
)
−
f
(
2
j
−
1
)
(
p
)
)
+
R
k
.
{\displaystyle \sum _{i=p}^{q}f\left(i\right)={\frac {f\left(p\right)+f\left(q\right)}{2}}+\int _{p}^{q}f(x)~{\rm {d}}x+\sum _{j=1}^{k}{\frac {b_{2j}}{(2j)!}}\left(f^{(2j-1)}(q)-f^{(2j-1)}(p)\right)+R_{k}.}
Les nombres b 2j désignent les nombres de Bernoulli et le reste Rk s'exprime à l'aide du polynôme de Bernoulli B 2k :
R
k
=
−
1
(
2
k
)
!
∫
p
q
f
(
2
k
)
(
x
)
B
2
k
(
x
−
⌊
x
⌋
)
d
x
{\displaystyle R_{k}=-{1 \over (2k)!}\int _{p}^{q}f^{(2k)}(x)B_{2k}(x-\lfloor x\rfloor )~{\rm {d}}x}
.
La notation B 2k désigne le 2k -ième polynôme de Bernoulli et
B
2
k
(
x
−
⌊
x
⌋
)
{\displaystyle B_{2k}(x-\lfloor x\rfloor )}
en est une version périodisée , de période 1, égale à B 2k (x ) si 0 < x < 1.
B
2
(
x
)
=
x
2
−
x
+
1
6
B
4
(
x
)
=
x
4
−
2
x
3
+
x
2
−
1
30
B
6
(
x
)
=
x
6
−
3
x
5
+
5
2
x
4
−
1
2
x
2
+
1
42
{\displaystyle B_{2}(x)=x^{2}-x+{\frac {1}{6}}\qquad B_{4}(x)=x^{4}-2x^{3}+x^{2}-{\frac {1}{30}}\qquad B_{6}(x)=x^{6}-3x^{5}+{\frac {5}{2}}x^{4}-{\frac {1}{2}}x^{2}+{\frac {1}{42}}}
.
Les nombres de Bernoulli vérifient les égalités
b
2
k
=
B
2
k
(
0
)
=
B
2
k
(
1
)
{\displaystyle b_{2k}=B_{2k}(0)=B_{2k}(1)}
.
D'autres expressions du reste sont données plus loin si la fonction f est 2k + 1 fois dérivable ou 2k + 2 fois dérivable.
Si l'on somme de p à q − 1 les nombres f (i ) , on a :
∑
i
=
p
q
−
1
f
(
i
)
=
∫
p
q
f
(
x
)
d
x
+
f
(
p
)
−
f
(
q
)
2
+
∑
j
=
1
k
b
2
j
(
2
j
)
!
(
f
(
2
j
−
1
)
(
q
)
−
f
(
2
j
−
1
)
(
p
)
)
+
R
k
{\displaystyle \sum _{i=p}^{q-1}f(i)=\int _{p}^{q}f(x)~{\rm {d}}x+{\frac {f(p)-f(q)}{2}}+\sum _{j=1}^{k}{\frac {b_{2j}}{(2j)!}}\left(f^{(2j-1)}(q)-f^{(2j-1)}(p)\right)+R_{k}}
.
Si l'on considère les nombres de Bernoulli d'indice impair : b 1 = –1 / 2 et b 2j +1 = 0 si j > 1, on peut énoncer la formule d'Euler-Maclaurin de la manière suivante[ N 1] :
pour une fonction complexe f qui est r fois continûment dérivable (avec r > 0) :
∑
i
=
p
q
−
1
f
(
i
)
=
∫
p
q
f
(
x
)
d
x
+
∑
j
=
1
r
b
j
j
!
(
f
(
j
−
1
)
(
q
)
−
f
(
j
−
1
)
(
p
)
)
+
R
r
′
{\displaystyle \sum _{i=p}^{q-1}f\left(i\right)=\int _{p}^{q}f(x)~{\rm {d}}x+\sum _{j=1}^{r}{\frac {b_{j}}{j!}}\left(f^{(j-1)}(q)-f^{(j-1)}(p)\right)+R'_{r}}
On a
R
r
′
=
R
[
r
2
]
,
R
0
=
R
1
′
et
R
k
=
R
2
k
′
=
R
2
k
+
1
′
{\displaystyle R'_{r}=R_{[{\frac {r}{2}}]},\qquad R_{0}=R'_{1}\qquad {\text{et}}\qquad R_{k}=R'_{2k}=R'_{2k+1}}
si k > 0.
Avec les notations précédentes, l'expression du reste R'r pour une fonction complexe r fois continûment dérivable (avec r > 0) est la suivante[ N 1] :
R
r
′
=
(
−
1
)
r
−
1
r
!
∫
p
q
f
(
r
)
(
x
)
B
r
(
x
−
⌊
x
⌋
)
d
x
{\displaystyle R'_{r}={\frac {(-1)^{r-1}}{r!}}\int _{p}^{q}f^{(r)}(x)B_{r}(x-\lfloor x\rfloor )~{\rm {d}}x}
.
La notation Br désigne le r -ème polynôme de Bernoulli , et
B
r
(
x
−
⌊
x
⌋
)
{\displaystyle B_{r}(x-\lfloor x\rfloor )}
en est une version périodisée , de période 1, égale à Br (x ) si 0 < x < 1.
B
0
(
x
)
=
1
B
1
(
x
)
=
x
−
1
2
B
2
(
x
)
=
x
2
−
x
+
1
6
B
3
(
x
)
=
x
3
−
3
2
x
2
+
1
2
x
B
4
(
x
)
=
x
4
−
2
x
3
+
x
2
−
1
30
{\displaystyle B_{0}(x)=1\qquad B_{1}(x)=x-{\frac {1}{2}}\qquad B_{2}(x)=x^{2}-x+{\frac {1}{6}}\qquad B_{3}(x)=x^{3}-{\frac {3}{2}}x^{2}+{\frac {1}{2}}x\qquad B_{4}(x)=x^{4}-2x^{3}+x^{2}-{\frac {1}{30}}}
.
(Les polynômes et les nombres de Bernoulli sont reliés par les égalités br = Br (0) = Br (1) si r > 1.).
Une autre formulation équivalente, où l'on somme de p + 1 à q , est donnée par Tenenbaum[ 2] : pour une fonction complexe f qui est r + 1 fois continûment dérivable (avec r ≥ 0) :
∑
p
<
n
⩽
q
f
(
n
)
=
∫
p
q
f
(
x
)
d
x
+
∑
j
=
0
r
(
−
1
)
j
+
1
b
j
+
1
(
j
+
1
)
!
(
f
(
j
)
(
q
)
−
f
(
j
)
(
p
)
)
+
R
r
+
1
′
{\displaystyle \sum _{p<n\leqslant q}f\left(n\right)=\int _{p}^{q}f(x)~{\rm {d}}x+\sum _{j=0}^{r}{\frac {(-1)^{j+1}b_{j+1}}{(j+1)!}}\left(f^{(j)}(q)-f^{(j)}(p)\right)+R'_{r+1}}
.
Le coefficient (–1)j +1 n'intervient dans la formule que pour j = 0. Son rôle est de remplacer le nombre de Bernoulli d'indice 1, b 1 = –1 / 2 , par b' 1 = +1 / 2 .
Si f est 2k fois continûment dérivable sur le segment [p , q ] (avec k ≥ 1), le reste Rk s'exprime de la manière suivante :
R
k
=
−
1
(
2
k
)
!
∫
p
q
f
(
2
k
)
(
x
)
B
2
k
(
x
−
⌊
x
⌋
)
d
x
.
{\displaystyle R_{k}=-{1 \over (2k)!}\int _{p}^{q}f^{(2k)}(x)B_{2k}(x-\lfloor x\rfloor )~{\rm {d}}x.}
Si f est 2k + 1 fois continûment dérivable sur le segment [p , q ] (avec k ≥ 0), le reste Rk s'exprime comme suit[ N 2] :
R
k
=
1
(
2
k
+
1
)
!
∫
p
q
f
(
2
k
+
1
)
(
x
)
B
2
k
+
1
(
x
−
⌊
x
⌋
)
d
x
.
{\displaystyle R_{k}={1 \over (2k+1)!}\int _{p}^{q}f^{(2k+1)}(x){B_{2k+1}(x-\lfloor x\rfloor )}~{\rm {d}}x.}
Si f est une fonction réelle 2k + 2 fois continûment dérivable sur le segment [p , q ] (avec k ≥ 0), le reste peut s'écrire des manières suivantes[ 3] :
R
k
=
−
1
(
2
k
+
2
)
!
∫
p
q
f
(
2
k
+
2
)
(
x
)
(
B
2
k
+
2
(
x
−
⌊
x
⌋
)
−
b
2
k
+
2
)
d
x
=
q
−
p
(
2
k
+
2
)
!
b
2
k
+
2
f
(
2
k
+
2
)
(
ξ
)
,
avec
ξ
∈
]
p
,
q
[
.
{\displaystyle {\begin{aligned}R_{k}&=-{1 \over (2k+2)!}\int _{p}^{q}f^{(2k+2)}(x)(B_{2k+2}(x-\lfloor x\rfloor )-b_{2k+2})~{\rm {d}}x\\&={q-p \over (2k+2)!}b_{2k+2}f^{(2k+2)}(\xi ),\qquad {\text{ avec }}\quad \xi \in \left]p,q\right[.\end{aligned}}}
Si l'on considère les nombres de Bernoulli sans leur signe (on a
b
2
k
=
(
−
1
)
k
−
1
|
b
2
k
|
{\displaystyle b_{2k}=(-1)^{k-1}|b_{2k}|}
), la dernière formule s'écrit[ 4] :
R
k
=
(
−
1
)
k
q
−
p
(
2
k
+
2
)
!
|
b
2
k
+
2
|
f
(
2
k
+
2
)
(
ξ
)
,
avec
ξ
∈
]
p
,
q
[
.
{\displaystyle R_{k}=(-1)^{k}{q-p \over (2k+2)!}|b_{2k+2}|f^{(2k+2)}(\xi ),\qquad {\text{ avec }}\quad \xi \in \left]p,q\right[.}
Remarque : le reste Rk est nul pour tout polynôme de degré au plus 2k + 1.
Si f est une fonction complexe 1 fois continûment dérivable sur le segment [p , q ], le maximum
sup
t
∈
[
0
;
1
]
|
B
1
(
t
)
|
=
1
2
{\displaystyle \sup _{t\in [0\,;\,1]}|B_{1}(t)|={\frac {1}{2}}}
du polynôme de Bernoulli B 1 (t ) permet d'obtenir la majoration :
|
R
0
|
⩽
1
2
∫
p
q
|
f
′
(
t
)
|
d
t
{\displaystyle |R_{0}|\leqslant {\frac {1}{2}}\int _{p}^{q}\left|f'(t)\right|\mathrm {d} t}
.
Si f est une fonction complexe 2k fois continûment dérivable sur le segment [p , q ] (avec k ≥ 1), on peut majorer le reste (ou « terme d'erreur ») de la formule d'Euler-Maclaurin en utilisant la majoration des polynômes de Bernoulli d'indice pair[ N 3] :
sup
t
∈
[
0
;
1
]
|
B
2
k
(
t
)
|
=
|
b
2
k
|
{\displaystyle \sup _{t\in [0\,;\,1]}|B_{2k}(t)|=|b_{2k}|}
:
|
R
k
|
⩽
|
b
2
k
|
(
2
k
)
!
∫
p
q
|
f
(
2
k
)
(
x
)
|
d
x
{\displaystyle |R_{k}|\leqslant {\frac {|b_{2k}|}{(2k)!}}\int _{p}^{q}\left|f^{(2k)}(x)\right|\mathrm {d} x}
.
Par exemple, avec
b
2
=
1
6
{\displaystyle b_{2}={\frac {1}{6}}}
, on a :
|
R
1
|
⩽
1
12
∫
p
q
|
f
″
(
x
)
|
d
x
{\displaystyle |R_{1}|\leqslant {\frac {1}{12}}\int _{p}^{q}\left|f''(x)\right|\mathrm {d} x}
.
L'inégalité peut être réécrite en utilisant la formule due à Euler (pour k ≥ 1) :
|
b
2
k
|
(
2
k
)
!
=
2
(
2
π
)
2
k
∑
i
=
1
∞
1
i
2
k
=
2
(
2
π
)
2
k
ζ
(
2
k
)
{\displaystyle {\frac {|b_{2k}|}{(2k)!}}={\frac {2}{(2\pi )^{2k}}}\sum _{i=1}^{\infty }{\frac {1}{i^{2k}}}={\frac {2}{(2\pi )^{2k}}}\zeta (2k)}
. On déduit que
|
b
2
k
|
(
2
k
)
!
⩽
2
ζ
(
2
)
(
2
π
)
2
k
{\displaystyle {\frac {|b_{2k}|}{(2k)!}}\leqslant {\frac {2\zeta (2)}{(2\pi )^{2k}}}}
(on a l'équivalent :
|
b
2
k
|
(
2
k
)
!
∼
2
(
2
π
)
2
k
{\displaystyle {\frac {|b_{2k}|}{(2k)!}}\sim {\frac {2}{(2\pi )^{2k}}}}
).
|
R
k
|
⩽
2
ζ
(
2
)
(
2
π
)
2
k
∫
p
q
|
f
(
2
k
)
(
t
)
|
d
t
{\displaystyle |R_{k}|\leqslant {\frac {2\zeta (2)}{(2\pi )^{2k}}}\int _{p}^{q}\left|f^{(2k)}(t)\right|\mathrm {d} t}
.
Si f est une fonction complexe 2k + 1 fois continûment dérivable sur le segment [p , q ] (avec k ≥ 1), en utilisant la majoration des polynômes de Bernoulli d'indice impair :
|
B
2
k
+
1
(
t
)
|
⩽
2
(
2
k
+
1
)
!
(
2
π
)
2
k
+
1
{\displaystyle |B_{2k+1}(t)|\leqslant {\frac {2(2k+1)!}{(2\pi )^{2k+1}}}}
(si
t
∈
[
0
;
1
]
{\displaystyle t\in [0\,;\,1]}
) démontrée par Derrick Lehmer , on obtient l'inégalité[ N 4] :
|
R
k
|
⩽
2
(
2
π
)
2
k
+
1
∫
p
q
|
f
(
2
k
+
1
)
(
x
)
|
d
x
{\displaystyle |R_{k}|\leqslant {\frac {2}{(2\pi )^{2k+1}}}\int _{p}^{q}\left|f^{(2k+1)}(x)\right|\mathrm {d} x}
Pour le polynôme de Bernoulli B 3 (t ) , on a le maximum
sup
t
∈
[
0
;
1
]
|
B
3
(
t
)
|
=
3
36
{\displaystyle \sup _{t\in [0\,;\,1]}|B_{3}(t)|={\frac {\sqrt {3}}{36}}}
qui permet d'obtenir la majoration :
|
R
1
|
⩽
3
216
∫
p
q
|
f
(
3
)
(
t
)
|
d
t
{\displaystyle |R_{1}|\leqslant {\frac {\sqrt {3}}{216}}\int _{p}^{q}\left|f^{(3)}(t)\right|\mathrm {d} t}
.
Si f est une fonction réelle 2k + 2 fois continûment dérivable sur le segment [p , q ] (avec k ≥ 0), dont la dérivée d'ordre 2k + 2 est de signe constant[ N 5] , le reste Rk a le même signe que le « premier terme négligé[ 5] » :
R
k
≈
b
2
k
+
2
(
2
k
+
2
)
!
(
f
(
2
k
+
1
)
(
q
)
−
f
(
2
k
+
1
)
(
p
)
)
≈
(
−
1
)
k
(
f
(
2
k
+
1
)
(
q
)
−
f
(
2
k
+
1
)
(
p
)
)
{\displaystyle R_{k}\approx {\frac {b_{2k+2}}{(2k+2)!}}\left(f^{(2k+1)}(q)-f^{(2k+1)}(p)\right)\approx (-1)^{k}\left(f^{(2k+1)}(q)-f^{(2k+1)}(p)\right)}
.
De plus on a les majorations suivantes[ 5] :
|
R
k
|
⩽
2
(
1
−
2
−
(
2
k
+
2
)
)
|
b
2
k
+
2
|
(
2
k
+
2
)
!
|
f
(
2
k
+
1
)
(
q
)
−
f
(
2
k
+
1
)
(
p
)
|
{\displaystyle |R_{k}|\leqslant 2\left(1-2^{-(2k+2)}\right){\frac {|b_{2k+2}|}{(2k+2)!}}\left|f^{(2k+1)}(q)-f^{(2k+1)}(p)\right|}
.
et
|
R
k
+
1
|
⩽
|
b
2
k
+
2
|
(
2
k
+
2
)
!
|
f
(
2
k
+
1
)
(
q
)
−
f
(
2
k
+
1
)
(
p
)
|
{\displaystyle |R_{k+1}|\leqslant {\frac {|b_{2k+2}|}{(2k+2)!}}\left|f^{(2k+1)}(q)-f^{(2k+1)}(p)\right|}
. (le reste suivant R k + 1 n'excède pas (en valeur absolue) le « dernier terme retenu »[ 6] ).
Si f est une fonction réelle 2k + 4 fois continûment dérivable sur le segment [p , q ] (avec k ≥ 0), dont les dérivées d'ordre 2k + 2 et 2k + 4 sont de signe constant et de même signe, alors les restes Rk et R k + 1 sont de signes opposés et le reste Rk est majoré (en valeur absolue) par le premier terme négligé[ 5] :
|
R
k
|
⩽
|
b
2
k
+
2
|
(
2
k
+
2
)
!
|
f
(
2
k
+
1
)
(
q
)
−
f
(
2
k
+
1
)
(
p
)
|
{\displaystyle |R_{k}|\leqslant {\frac {|b_{2k+2}|}{(2k+2)!}}\left|f^{(2k+1)}(q)-f^{(2k+1)}(p)\right|}
.
On démontre la formule
R
k
=
−
1
(
2
k
)
!
∫
p
q
f
(
2
k
)
(
x
)
B
2
k
(
x
−
⌊
x
⌋
)
d
x
.
{\displaystyle R_{k}=-{1 \over (2k)!}\int _{p}^{q}f^{(2k)}(x)B_{2k}(x-\lfloor x\rfloor )~{\rm {d}}x.}
sur l'intervalle [n , n + 1], avec n ∈ ℤ, puis on déduit la formule précédente par sommation sur n ∈ ℤ (p ≤ n ≤ q – 1).
Démonstration
Soit g une fonction continûment dérivable sur [n , n + 1]. En utilisant la propriété des polynômes de Bernoulli :
∀
m
∈
N
B
m
+
1
′
=
(
m
+
1
)
B
m
{\displaystyle \forall m\in \mathbb {N} \quad B_{m+1}'=\left(m+1\right)B_{m}}
, on trouve en faisant une intégration par parties :
∫
n
n
+
1
g
(
t
)
B
m
(
t
−
n
)
d
t
=
[
g
(
t
)
B
m
+
1
(
t
−
n
)
m
+
1
]
n
n
+
1
−
∫
n
n
+
1
g
′
(
t
)
B
m
+
1
(
t
−
n
)
m
+
1
d
t
.
{\displaystyle \int _{n}^{n+1}g(t)B_{m}(t-n)~{\rm {d}}t=\left[g(t){\frac {B_{m+1}(t-n)}{m+1}}\right]_{n}^{n+1}-\int _{n}^{n+1}g'(t){\frac {B_{m+1}(t-n)}{m+1}}~{\rm {d}}t.}
Or, sachant que pour m ≥ 1, on a
B
m
+
1
(
1
)
=
B
m
+
1
(
0
)
=
b
m
+
1
{\displaystyle B_{m+1}(1)=B_{m+1}(0)=b_{m+1}}
, on en déduit :
pour m ≥ 1 :
∫
n
n
+
1
g
(
t
)
B
m
(
t
−
n
)
d
t
=
b
m
+
1
m
+
1
(
g
(
n
+
1
)
−
g
(
n
)
)
−
1
m
+
1
∫
n
n
+
1
g
′
(
t
)
B
m
+
1
(
t
−
n
)
d
t
.
{\displaystyle \int _{n}^{n+1}g(t)B_{m}(t-n)~{\rm {d}}t={\frac {b_{m+1}}{m+1}}\left(g(n+1)-g(n)\right)-{\frac {1}{m+1}}\int _{n}^{n+1}g'(t)B_{m+1}(t-n)~{\rm {d}}t.}
Sachant que (pour
m = 0) on a
B
0
(
t
)
=
1
et
B
1
(
1
)
=
1
2
et
B
1
(
0
)
=
−
1
2
,
{\displaystyle B_{0}(t)=1{\text{ et }}B_{1}(1)={\frac {1}{2}}{\text{ et }}B_{1}(0)=-{\frac {1}{2}},}
on en déduit si
g est dérivable continûment deux fois :
∫
n
n
+
1
g
(
t
)
d
t
=
∫
n
n
+
1
g
(
t
)
B
0
(
t
−
n
)
d
t
=
1
2
(
g
(
n
+
1
)
+
g
(
n
)
)
−
∫
n
n
+
1
g
′
(
t
)
B
1
(
t
−
n
)
d
t
=
g
(
n
+
1
)
+
g
(
n
)
2
−
b
2
2
(
g
′
(
n
+
1
)
−
g
′
(
n
)
)
+
1
2
∫
n
n
+
1
g
″
(
t
)
B
2
(
t
−
n
)
d
t
.
{\displaystyle {\begin{aligned}\int _{n}^{n+1}g(t)~{\rm {d}}t&=\int _{n}^{n+1}g(t)B_{0}(t-n)~{\rm {d}}t\\&={\frac {1}{2}}\left(g(n+1)+g(n)\right)-\int _{n}^{n+1}g'(t)B_{1}(t-n)~{\rm {d}}t\\&={\frac {g(n+1)+g(n)}{2}}-{\frac {b_{2}}{2}}\left(g'(n+1)-g'(n)\right)+{\frac {1}{2}}\int _{n}^{n+1}g''(t){B_{2}(t-n)}~{\rm {d}}t.\end{aligned}}}
Soit f une fonction continûment dérivable 2k fois (k > 0). Par récurrence sur k , on montre :
∫
n
n
+
1
f
(
t
)
d
t
=
f
(
n
)
+
f
(
n
+
1
)
2
+
∑
i
=
2
2
k
(
−
1
)
i
−
1
b
i
i
!
(
f
(
i
−
1
)
(
n
+
1
)
−
f
(
i
−
1
)
(
n
)
)
+
1
(
2
k
)
!
∫
n
n
+
1
f
(
2
k
)
(
t
)
B
2
k
(
t
−
n
)
d
t
.
{\displaystyle \int _{n}^{n+1}f(t)~{\rm {d}}t={\frac {f(n)+f(n+1)}{2}}+\sum _{i=2}^{2k}{\frac {\left(-1\right)^{i-1}b_{i}}{i!}}\left(f^{(i-1)}(n+1)-f^{(i-1)}(n)\right)+{\frac {1}{(2k)!}}\int _{n}^{n+1}f^{(2k)}(t)B_{2k}(t-n)~{\rm {d}}t.}
Enfin, avec la propriété :
∀
j
≥
1
,
b
2
j
+
1
=
0
{\displaystyle \forall j\geq 1,b_{2j+1}=0}
, on en déduit (avec i = 2j ) :
∫
n
n
+
1
f
(
t
)
d
t
=
f
(
n
)
+
f
(
n
+
1
)
2
−
∑
j
=
1
k
b
2
j
(
2
j
)
!
(
f
(
2
j
−
1
)
(
n
+
1
)
−
f
(
2
j
−
1
)
(
n
)
)
+
1
(
2
k
)
!
∫
n
n
+
1
f
(
2
k
)
(
t
)
B
2
k
(
t
−
n
)
d
t
{\displaystyle \int _{n}^{n+1}f(t)~{\rm {d}}t={\frac {f(n)+f(n+1)}{2}}-\sum _{j=1}^{k}{\frac {b_{2j}}{(2j)!}}\left(f^{(2j-1)}(n+1)-f^{(2j-1)}(n)\right)+{\frac {1}{(2k)!}}\int _{n}^{n+1}f^{(2k)}(t)B_{2k}(t-n)~{\rm {d}}t}
donc, avec n = ⌊t ⌋
f
(
n
)
+
f
(
n
+
1
)
2
=
∫
n
n
+
1
f
(
t
)
d
t
+
∑
j
=
1
k
b
2
j
(
2
j
)
!
(
f
(
2
j
−
1
)
(
n
+
1
)
−
f
(
2
j
−
1
)
(
n
)
)
−
1
(
2
k
)
!
∫
n
n
+
1
f
(
2
k
)
(
t
)
B
2
k
(
t
−
⌊
t
⌋
)
d
t
.
{\displaystyle {\frac {f(n)+f(n+1)}{2}}=\int _{n}^{n+1}f(t)~{\rm {d}}t+\sum _{j=1}^{k}{\frac {b_{2j}}{(2j)!}}\left(f^{(2j-1)}(n+1)-f^{(2j-1)}(n)\right)-{\frac {1}{(2k)!}}\int _{n}^{n+1}f^{(2k)}(t)B_{2k}\left(t-\lfloor t\rfloor \right)~{\rm {d}}t.}
D'où par sommation sur sommation sur n ∈ ℤ (p ≤ n ≤ q – 1) :
∑
n
=
p
q
−
1
f
(
n
)
+
f
(
n
+
1
)
2
=
∫
p
q
f
(
t
)
d
t
+
∑
j
=
1
k
b
2
j
(
2
j
)
!
(
f
(
2
j
−
1
)
(
q
)
−
f
(
2
j
−
1
)
(
p
)
)
−
1
(
2
k
)
!
∫
p
q
f
(
2
k
)
(
t
)
B
2
k
(
t
−
⌊
t
⌋
)
d
t
.
{\displaystyle \sum _{n=p}^{q-1}{\frac {f(n)+f(n+1)}{2}}=\int _{p}^{q}f(t)~{\rm {d}}t+\sum _{j=1}^{k}{\frac {b_{2j}}{(2j)!}}\left(f^{(2j-1)}(q)-f^{(2j-1)}(p)\right)-{\frac {1}{(2k)!}}\int _{p}^{q}f^{(2k)}(t)B_{2k}\left(t-\lfloor t\rfloor \right)~{\rm {d}}t.}
Application à l'intégration numérique
modifier
La formule sommatoire peut être utilisée pour approcher des intégrales par un procédé discret, par exemple dans la méthode des trapèzes ou celle de Romberg , ou à l'inverse pour transformer une somme discrète (finie ou non) et lui appliquer les techniques du calcul infinitésimal .
La formule d'Euler-Maclaurin peut aussi être utilisée pour une estimation précise de l'erreur commise dans le calcul numérique d'une intégrale ; en particulier, c'est sur elle que reposent les méthodes d'extrapolation. La méthode de quadrature de Clenshaw-Curtis est essentiellement un changement de variables ramenant une intégrale arbitraire à l'intégration de fonctions périodiques, pour lesquelles la formule sommatoire est très précise (dans ce cas, elle prend la forme d'une transformée en cosinus discrète ).
Dans la méthode des trapèzes , on approxime l'intégrale
∫
n
n
+
1
f
(
x
)
d
x
{\displaystyle \int _{n}^{n+1}f(x)~{\rm {d}}x}
par interpolation linéaire sur chaque intervalle [n , n + 1] :
∫
n
n
+
1
f
(
x
)
d
x
≈
f
(
n
)
+
f
(
n
+
1
)
2
{\displaystyle \int _{n}^{n+1}f(x)~{\rm {d}}x\approx {\frac {f(n)+f(n+1)}{2}}}
.
En sommant sur tous les intervalles de longueur 1, on approxime l'intégrale
∫
p
q
f
(
x
)
d
x
{\displaystyle \int _{p}^{q}f(x)~{\rm {d}}x}
par la somme
f
(
p
)
2
+
f
(
p
+
1
)
+
f
(
p
+
2
)
+
…
+
f
(
q
−
1
)
+
f
(
q
)
2
=
f
(
p
)
+
f
(
q
)
2
+
∑
i
=
p
+
1
q
−
1
f
(
i
)
{\displaystyle {\frac {f(p)}{2}}+f(p+1)+f(p+2)+\ldots +f(q-1)+{\frac {f(q)}{2}}={\frac {f\left(p\right)+f\left(q\right)}{2}}+\sum _{i=p+1}^{q-1}f\left(i\right)}
La formule d'Euler-Maclaurin peut s'écrire :
∫
p
q
f
(
x
)
d
x
=
f
(
p
)
+
f
(
q
)
2
+
∑
i
=
p
+
1
q
−
1
f
(
i
)
−
∑
j
=
1
k
b
2
j
(
2
j
)
!
(
f
(
2
j
−
1
)
(
q
)
−
f
(
2
j
−
1
)
(
p
)
)
−
R
k
{\displaystyle \int _{p}^{q}f(x)~{\rm {d}}x={\frac {f\left(p\right)+f\left(q\right)}{2}}+\sum _{i=p+1}^{q-1}f\left(i\right)-\sum _{j=1}^{k}{\frac {b_{2j}}{(2j)!}}\left(f^{(2j-1)}(q)-f^{(2j-1)}(p)\right)-R_{k}}
Intégration sur un intervalle quelconque
modifier
Un simple changement de variable permet d'obtenir une formule analogue pour une fonction définie sur un segment à bornes non entières. Les restes sont donnés avec le « point moyen »
ξ
∈
]
a
,
b
[
{\displaystyle \;\xi \in ~]a,b[}
pour une fonction dérivable 2k + 2 fois.
En posant
h
=
b
−
a
N
{\displaystyle h={\frac {b-a}{N}}}
, on a[ 7] :
∫
a
b
f
(
x
)
d
x
=
h
f
(
a
)
+
f
(
b
)
2
+
h
∑
m
=
1
N
−
1
f
(
a
+
m
h
)
−
∑
j
=
1
k
b
2
j
(
2
j
)
!
h
2
j
(
f
(
2
j
−
1
)
(
b
)
−
f
(
2
j
−
1
)
(
a
)
)
−
(
b
−
a
)
b
2
k
+
2
(
2
k
+
2
)
!
f
(
2
k
+
2
)
(
ξ
)
h
2
k
+
2
{\displaystyle \int _{a}^{b}f(x)~{\rm {d}}x=h{\frac {f\left(a\right)+f\left(b\right)}{2}}+h\sum _{m=1}^{N-1}f\left(a+mh\right)-\sum _{j=1}^{k}{\frac {b_{2j}}{(2j)!}}h^{2j}\left(f^{(2j-1)}(b)-f^{(2j-1)}(a)\right)-(b-a){\frac {b_{2k+2}}{(2k+2)!}}f^{(2k+2)}(\xi )h^{2k+2}}
.
Le « terme d'erreur » peut également s'écrire :
−
R
k
=
−
(
b
−
a
)
b
2
k
+
2
(
2
k
+
2
)
!
f
(
2
k
+
2
)
(
ξ
)
h
2
k
+
2
=
−
b
2
k
+
2
(
b
−
a
)
2
k
+
3
(
2
k
+
2
)
!
N
2
k
+
2
f
(
2
k
+
2
)
(
ξ
)
{\displaystyle -R_{k}=-(b-a){\frac {b_{2k+2}}{(2k+2)!}}f^{(2k+2)}(\xi )h^{2k+2}=-{\frac {b_{2k+2}(b-a)^{2k+3}}{(2k+2)!N^{2k+2}}}f^{(2k+2)}(\xi )}
.
Si N = 1 , on a une formule où n'interviennent que les extrémités a et b [ 8] :
∫
a
b
f
(
x
)
d
x
=
(
b
−
a
)
f
(
a
)
+
f
(
b
)
2
−
∑
j
=
1
k
(
b
−
a
)
2
j
(
2
j
)
!
b
2
j
(
f
(
2
j
−
1
)
(
b
)
−
f
(
2
j
−
1
)
(
a
)
)
−
(
b
−
a
)
2
k
+
3
b
2
k
+
2
(
2
k
+
2
)
!
f
(
2
k
+
2
)
(
ξ
)
{\displaystyle \int _{a}^{b}f(x)~{\rm {d}}x=(b-a){\frac {f(a)+f(b)}{2}}-\sum _{j=1}^{k}{\frac {(b-a)^{2j}}{(2j)!}}b_{2j}\left(f^{(2j-1)}(b)-f^{(2j-1)}(a)\right)-(b-a)^{2k+3}{\frac {b_{2k+2}}{(2k+2)!}}f^{(2k+2)}(\xi )}
Expressions du reste pour k = 0 et pour k = 1
modifier
Expressions de R 0 et erreur de la méthode des trapèzes
modifier
Les premiers polynômes de Bernoulli sont :
B
0
(
y
)
=
1
,
B
1
(
y
)
=
y
−
1
2
,
B
2
(
y
)
=
y
2
−
y
+
1
6
{\displaystyle B_{0}(y)=1,\quad B_{1}(y)=y-{\frac {1}{2}},\quad B_{2}(y)=y^{2}-y+{\frac {1}{6}}}
.
R
0
=
f
(
p
)
+
f
(
q
)
2
+
∑
n
=
p
+
1
q
−
1
f
(
n
)
−
∫
p
q
f
(
x
)
d
x
{\displaystyle R_{0}={\frac {f\left(p\right)+f\left(q\right)}{2}}+\sum _{n=p+1}^{q-1}f\left(n\right)-\int _{p}^{q}f(x)~{\rm {d}}x}
.
(-R 0 ) est l'erreur faite en approximant l'intégrale
∫
n
n
+
1
f
(
x
)
d
x
≈
f
(
n
)
+
f
(
n
+
1
)
2
{\displaystyle \int _{n}^{n+1}f(x)~{\rm {d}}x\approx {\frac {f(n)+f(n+1)}{2}}}
par la méthode des trapèzes sur chaque intervalle [n , n + 1].
Si f est continûment dérivable une fois (on pose
y
=
x
−
⌊
x
⌋
{\displaystyle y=x-\lfloor x\rfloor }
[ 9] ) :
R
0
=
∫
p
q
f
′
(
x
)
B
1
(
y
)
d
x
=
∫
p
q
f
′
(
x
)
(
y
−
1
2
)
d
x
{\displaystyle R_{0}=\int _{p}^{q}f'(x){B_{1}(y)}~{\rm {d}}x=\int _{p}^{q}f'(x)\left(y-{\tfrac {1}{2}}\right)~{\rm {d}}x}
Si f est une fonction réelle continûment dérivable deux fois[ 10] :
R
0
=
−
1
2
!
∫
p
q
f
″
(
x
)
(
B
2
(
y
)
−
1
6
)
d
x
=
−
1
2
∫
p
q
f
″
(
x
)
(
y
2
−
y
)
d
x
=
q
−
p
12
f
″
(
ξ
)
(
avec
ξ
∈
]
p
,
q
[
)
.
{\displaystyle {\begin{aligned}R_{0}=-{\frac {1}{2!}}\int _{p}^{q}f''(x)\left(B_{2}(y)-{\tfrac {1}{6}}\right)~{\rm {d}}x&=-{\frac {1}{2}}\int _{p}^{q}f''(x)(y^{2}-y)~{\rm {d}}x\\&={\frac {q-p}{12}}f''(\xi )\qquad ({\text{avec }}\;\xi \in \left]p,q\right[).\end{aligned}}}
Les polynômes de Bernoulli qui interviennent sont[ N 6] :
B
2
(
y
)
=
y
2
−
y
+
1
6
,
B
3
(
y
)
=
y
3
−
3
2
y
2
+
1
2
y
=
y
(
y
−
1
)
(
y
−
1
2
)
,
B
4
(
y
)
=
y
4
−
2
y
3
+
y
2
−
1
30
=
(
y
2
−
y
)
2
−
1
30
{\displaystyle B_{2}(y)=y^{2}-y+{\frac {1}{6}},\quad B_{3}(y)=y^{3}-{\frac {3}{2}}y^{2}+{\frac {1}{2}}y=y(y-1)\left(y-{\tfrac {1}{2}}\right),\quad B_{4}(y)=y^{4}-2y^{3}+y^{2}-{\frac {1}{30}}=(y^{2}-y)^{2}-{\frac {1}{30}}}
.
R
1
=
f
(
p
)
+
f
(
q
)
2
+
∑
n
=
p
+
1
q
−
1
f
(
n
)
−
∫
p
q
f
(
x
)
d
x
−
1
12
(
f
′
(
q
)
−
f
′
(
p
)
)
{\displaystyle R_{1}={\frac {f\left(p\right)+f\left(q\right)}{2}}+\sum _{n=p+1}^{q-1}f\left(n\right)-\int _{p}^{q}f(x)~{\rm {d}}x-{\frac {1}{12}}(f'(q)-f'(p))}
.
(–R 1 ) est le terme d'erreur correspondant à la formule de quadrature[ 11] , exacte pour les polynômes de degré inférieur ou égal à trois :
∫
p
q
f
(
x
)
d
x
≈
f
(
p
)
+
f
(
q
)
2
+
∑
n
=
p
+
1
q
−
1
f
(
n
)
+
f
′
(
p
)
−
f
′
(
q
)
12
{\displaystyle \int _{p}^{q}f(x)~{\rm {d}}x\approx {\frac {f\left(p\right)+f\left(q\right)}{2}}+\sum _{n=p+1}^{q-1}f\left(n\right)+{\frac {f'(p)-f'(q)}{12}}}
(dans l'article Calcul numérique d'une intégrale , il s'agit de la formule de Newton-Coates généralisée NC-2-2).
Si f est continûment dérivable deux fois (en posant
y
=
x
−
⌊
x
⌋
{\displaystyle y=x-\lfloor x\rfloor }
[ 9] ) :
R
1
=
−
1
2
!
∫
p
q
f
″
(
x
)
B
2
(
y
)
d
x
=
−
1
2
∫
p
q
f
″
(
x
)
(
y
2
−
y
+
1
6
)
d
x
.
{\displaystyle R_{1}=-{\frac {1}{2!}}\int _{p}^{q}f''(x){B_{2}(y)}~{\rm {d}}x=-{\frac {1}{2}}\int _{p}^{q}f''(x)\left(y^{2}-y+{\tfrac {1}{6}}\right)~{\rm {d}}x.}
Si f est continûment dérivable trois fois :
R
1
=
1
3
!
∫
p
q
f
‴
(
x
)
B
3
(
y
)
d
x
=
1
6
∫
p
q
f
‴
(
x
)
(
y
3
−
3
2
y
2
+
1
2
y
)
d
x
.
{\displaystyle R_{1}={\frac {1}{3!}}\int _{p}^{q}f'''(x){B_{3}(y)}~{\rm {d}}x={\frac {1}{6}}\int _{p}^{q}f'''(x)\left(y^{3}-{\tfrac {3}{2}}y^{2}+{\tfrac {1}{2}}y\right)~{\rm {d}}x.}
Si f est une fonction réelle continûment dérivable quatre fois :
R
1
=
−
1
4
!
∫
p
q
f
(
4
)
(
x
)
(
B
4
(
y
)
+
1
30
)
d
x
=
−
1
24
∫
p
q
f
(
4
)
(
x
)
(
y
2
−
y
)
2
d
x
=
−
q
−
p
720
f
(
4
)
(
ξ
)
(
avec
ξ
∈
]
p
,
q
[
)
.
{\displaystyle R_{1}=-{\frac {1}{4!}}\int _{p}^{q}f^{(4)}(x)\left(B_{4}(y)+{\tfrac {1}{30}}\right)~{\rm {d}}x=-{\frac {1}{24}}\int _{p}^{q}f^{(4)}(x)(y^{2}-y)^{2}~{\rm {d}}x=-{\frac {q-p}{720}}f^{(4)}(\xi )\qquad ({\text{avec }}\;\xi \in \left]p,q\right[).}
Le problème de Bâle demandait de déterminer la somme
1
+
1
4
+
1
9
+
1
16
+
1
25
+
⋯
=
∑
n
=
1
∞
1
n
2
.
{\displaystyle 1+{\frac {1}{4}}+{\frac {1}{9}}+{\frac {1}{16}}+{\frac {1}{25}}+\cdots =\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}.}
Euler calcula cette somme à 20 décimales en utilisant seulement quelques termes de la formule d'Euler-Maclaurin. Ce calcul le convainquit probablement qu'elle valait π2 / 6 , résultat qu'il publia en 1735 (mais avec des arguments incorrects ; il lui fallut six ans de plus pour trouver une démonstration rigoureuse )[ 12] .
Si f est un polynôme de degré d et si l'on applique la formule sommatoire avec p = 0, q = n et k choisi tel que d ≤ 2k +1 , le reste Rk disparaît.
Par exemple, si f (x ) = x 3 , on peut prendre k = 1 pour obtenir
∑
i
=
0
n
i
3
=
∫
0
n
x
3
d
x
+
0
3
+
n
3
2
+
1
6
3
n
2
−
0
2
!
=
n
4
4
+
n
3
2
+
n
2
4
=
n
4
+
2
n
3
+
n
2
4
=
(
n
(
n
+
1
)
2
)
2
{\displaystyle \sum _{i=0}^{n}i^{3}=\int _{0}^{n}x^{3}\mathrm {d} x+{\frac {0^{3}+n^{3}}{2}}+{\frac {1}{6}}{\frac {3n^{2}-0}{2!}}={\frac {n^{4}}{4}}+{\frac {n^{3}}{2}}+{\frac {n^{2}}{4}}={\frac {n^{4}+2n^{3}+n^{2}}{4}}=\left({\frac {n(n+1)}{2}}\right)^{2}}
.
Développements asymptotiques de fonctions définies par une série
modifier
Pour déterminer des développements asymptotiques de sommes et de séries, la forme la plus utile de la formule sommatoire est sans doute (pour a et b entiers) :
∑
n
=
a
b
f
(
n
)
∼
∫
a
b
f
(
x
)
d
x
+
f
(
a
)
+
f
(
b
)
2
+
∑
j
=
1
∞
B
2
j
(
2
j
)
!
(
f
(
2
j
−
1
)
(
b
)
−
f
(
2
j
−
1
)
(
a
)
)
{\displaystyle \sum _{n=a}^{b}f(n)\sim \int _{a}^{b}f(x)~{\rm {d}}x+{\frac {f(a)+f(b)}{2}}+\sum _{j=1}^{\infty }\,{\frac {B_{2j}}{(2j)!}}\left(f^{(2j-1)}(b)-f^{(2j-1)}(a)\right)}
Ce développement reste souvent valide même lorsque l'on prend les limites quand a → –∞ ou b → +∞ , ou les deux. Dans de nombreux cas, l'intégrale de droite peut être calculée de manière exacte avec des fonctions élémentaires, alors que ce n'est pas le cas de la somme.
L'écriture précédente doit être interprétée comme une série formelle , car, le plus souvent, cette série est divergente ; la formule ne peut en général pas être exploitée directement sous cette forme. Toutefois, Euler avait déjà remarqué qu'on obtenait une précision numérique remarquable en tronquant la formule au plus petit terme de la série, ce qui fut précisé et expliqué par les travaux d'Émile Borel [ 13] .
Par exemple :
∑
k
=
0
∞
1
(
z
+
k
)
2
⏟
=
ψ
1
(
z
)
∼
∫
0
∞
1
(
z
+
x
)
2
d
x
⏟
=
1
/
z
+
1
2
z
2
+
∑
j
=
1
∞
B
2
j
z
2
j
+
1
.
{\displaystyle \underbrace {\sum _{k=0}^{\infty }{\frac {1}{(z+k)^{2}}}} _{=\psi _{1}(z)}\sim \underbrace {\int _{0}^{\infty }{\frac {1}{(z+x)^{2}}}~{\rm {d}}x} _{=1/z}+{\frac {1}{2z^{2}}}+\sum _{j=1}^{\infty }{\frac {B_{2j}}{z^{2j+1}}}.\,}
Ici, le membre de gauche est égal à
ψ
1
(
z
)
{\displaystyle \psi _{1}(z)}
, c'est-à-dire à la fonction polygamma d'ordre 1 (appelée aussi fonction trigamma) définie à partir de la fonction Gamma :
ψ
1
(
z
)
=
d
2
d
z
2
ln
Γ
(
z
)
=
∑
k
=
0
∞
1
(
z
+
k
)
2
{\displaystyle \psi _{1}(z)={\frac {\mathrm {d} ^{2}}{\mathrm {d} z^{2}}}\ln \Gamma (z)=\sum _{k=0}^{\infty }{\frac {1}{(z+k)^{2}}}}
.
La formule d'Euler Maclaurin amène à un développement asymptotique de
ψ
1
(
z
)
{\displaystyle \psi _{1}(z)}
, lequel permet une estimation précise de l'erreur de la formule de Stirling pour la fonction Gamma :
ln
Γ
(
z
)
=
z
ln
z
−
z
−
1
2
ln
z
+
1
2
ln
(
2
π
)
+
∑
n
=
1
N
B
2
n
2
n
(
2
n
−
1
)
z
2
n
−
1
+
R
N
(
z
)
{\displaystyle \ln \Gamma (z)=z\ln z-z-{\frac {1}{2}}\ln {z}+{\frac {1}{2}}\ln(2\pi )+\sum _{n=1}^{N}{\frac {B_{2n}}{2n(2n-1)z^{2n-1}}}+R_{N}(z)}
.
Γ
(
z
)
=
2
π
z
(
z
e
)
z
(
1
+
O
(
1
z
)
)
.
{\displaystyle \Gamma (z)={\sqrt {\frac {2\pi }{z}}}\,{\left({\frac {z}{e}}\right)}^{z}\left(1+O\left({\frac {1}{z}}\right)\right).}
↑ a et b Cohen 2007 , p. 23. La formule avec son reste est donnée dans le cas, plus général, où les extrémités de l'intervalle sont des nombres réels a et a +N qui diffèrent d'un nombre entier naturel N . Cette formule est un cas particulier d'une formule donnée page 22 où les nombres de Bernoulli bj sont remplacés par les valeurs des polynômes de Bernoulli
B
r
(
x
−
⌊
x
⌋
)
{\displaystyle B_{r}(x-\lfloor x\rfloor )}
aux extrémités a et b .
↑ Dieudonné 1980 , p. 303. Selon Dieudonné, il suffit que f admette une dérivée (2k + 1)-ème continue par morceaux.
↑ Dieudonné 1980 , p. 301. Dieudonné note Bk les coefficients
|
b
2
k
|
{\displaystyle |b_{2k}|}
.
↑ (en) Derrick Lehmer , « On the maxima and minima of Bernoulli polynomials », The American Mathematical Monthly , vol. 47, no 8, 1940 , p. 533–538 (DOI 10.2307/2303833 , JSTOR 2303833 ) .
↑ Selon Bourbaki, FVR , p. VI.26, il suffit que f soit 2k + 1 fois dérivable et que la dérivée d'ordre 2k + 1 soit monotone.
↑ Deheuvels 1980 , p. 189. Noter que dans la notation de Deheuvels, les nombres de Bernoulli sont considérés sans leur signe (
B
2
k
=
|
b
2
k
|
{\displaystyle B_{2k}=|b_{2k}|}
).
↑ « Formule d'Euler-MacLaurin [sic] », sur bibmath.net/dico .
↑ Tenenbaum 2008 , p. 23.
↑ Rombaldi 2005 , p. 527 ; Hardy 1949 , p. 325.
↑ Deheuvels 1980 , p. 185-187 et p. 195.
↑ a b et c Cohen 2007 , p. 25.
↑ Tenenbaum 2008 , p. 26.
↑ Rombaldi 2005 , p. 330.
↑ Deheuvels 1980 , p. 188.
↑ a et b Hardy 1949 , p. 318.
↑ Hardy 1949 , p. 319.
↑ Rombaldi 2005 , p. 328.
↑ (en) David J. Pengelley , « Dances between continuous and discrete: Euler's summation formula » , dans Robert Bradley et Ed Sandifer, Proceedings, Euler 2K+2 Conference (Rumford, Maine, 2002) , Euler Society, 2003 (lire en ligne ) .
↑ Voir cette conférence de F. Pham
: document utilisé comme source pour la rédaction de cet article.
N. Bourbaki , Éléments de mathématique , Fonctions d'une variable réelle , chap. VI
(en) Henri Cohen , Number Theory , vol. II : Analytic and Modern Tools , Springer-Verlag , 2007
Paul Deheuvels , L'Intégrale , PUF , coll. « mathématiques », 1980
Jean-Pierre Demailly , Analyse numérique et équations différentielles , Presses universitaires de Grenoble
Jean Dieudonné , Calcul infinitésimal , Paris, Hermann , 1980
(en) G. H. Hardy , Divergent Series , OUP , 1949
Jean-Étienne Rombaldi, Interpolation et approximation , Vuibert , 2005
Gérald Tenenbaum , Introduction à la théorie analytique et probabiliste des nombres , Belin , 2008