Intégrales de Wallis

John Wallis, mathématicien anglais, est né en 1616 et est mort en 1703. Wallis est donc antérieur à Newton.
1) Définition.
On pose
Z π/2
∀n ∈ N, Wn = sinn t dt.
0
h πi
Wn existe pour tout entier naturel n car la fonction t 7→ sinn t est continue sur 0, .
2
2) Autres expressions de Wn .
π
Le changement de variables u = − t fournit
2
Z π/2
∀n ∈ N, Wn = cosn t dt.
0
i πh h π i
Soit ε un réel de 0, . La fonction u 7→ Arcsin u = t est de classe C1 sur 0, Arcsin −ε et on peut poser
2 Z π/2−ε Z Arcsin(π/2−ε) 2
un
t = Arcsin u ou encore u = sin t pour obtenir sinn t dt = √ du. Quand ε tend vers 0
Z π/2−ε 0 0 1 − u2
Z Arcsin(π/2−ε)
un
par valeurs supérieures, sinn t dt tend vers Wn et il en est de même de √ du de sorte que
Z1 0 0 1 − u2
un un
l’intégrale √ du converge. Comme la fonction u 7→ √ est positive sur [0, 1[, on en déduit que cette
0 1 − u2 1 − u2
fonction est intégrable sur [0, 1[. Quand ε tend vers 0, on obtient alors
Z1
un
∀n ∈ N, Wn = √ du.
0 1 − u2
On peut aussi poser u = sin t dans l’intégrale définissant W2n+1 pour obtenir
Z π/2 Z π/2 Z1
2n+1 2

n
n
W2n+1 = cos t dt = 1 − sin t cos t dt = 1 − u2 du.
0 0 0
Z1
∀n ∈ N, W2n+1 = (1 − u2 )n du.
0

3) Signe et sens de variation de la suite (Wn )n∈N .

i πh
Pour tout entier naturel n et tout réel t de 0, ,
2 h πi
Soit n ∈ N. La fonction t 7→ sinn t, est continue, positive et non nulle (mais pas strictement positive) sur 0, . Donc,
2
Wn > 0.
hLa fonction t 7→ sinn t − sinn+1 t = sint (1 − sin t), est continue, positive et non nulle (mais pas strictement positive) sur
πi
0, . Donc, Wn − Wn+1 > 0 et finalement 0 < Wn+1 < Wn .
2
La suite (Wn )n∈N est strictement positive et strictement décroissante.
4) Limite.
1ère idée. On montre « à la main » que lim Wn = 0.
i πh n→+∞

Soit ε > 0. Soit a un réel de 0, . Pour tout naturel n, on a

2
Za Z π/2 π 
n
0 6 Wn = sin t dt + sinn t dt 6 a sinn a + −a .
0 a 2
π π ε ε
[ de sorte que 0 < −a < . Pour tout entier naturel n, on a alors 0 6 Wn 6 a sinn a+ .
On choisit et on fixe a dans ]0,
2i 2 2 2
πh
Maintenant, puisque a est dans 0, , sin a est dans ]0, 1[ et donc lim a sinn a = 0.
2 n→+∞

1
 ε ε ε
Il existe ainsi un entier naturel n0 tel que, pour n > n0 , a sinn a 6
et donc Wn 6 + = ε.
2 2 2
On a montré que ∀ε > 0, ∃n0 ∈ N/ (∀n ∈ N), (n > n0 ⇒ 0 6 Wn 6 ε) et donc

lim Wn = 0.
n→+∞

2ème idée.
h π Oni utilise le théorème de convergence dominée pour atteindre le même h but.
πi
Pour t ∈ 0, et n ∈ N, on pose fn (t) = sinn t (avec la convention usuelle ∀t ∈ 0, , f0 (t) = 1).
2 h πi 2
• Chaque fonction fn est intégrable sur le segment 0, car continue sur ce segment.
h π2i
• La suite de fonction fn converge simplement sur 0, vers la fonction f définie par :
 2
 π
h πi  0 si t <
2
∀t ∈ 0, , f(t) = .
2  1 si t = π

2 h πi
De plus, f est continue par morceaux sur 0, .
h πi 2 h πi
• ∀n ∈ N, ∀t ∈ 0, , |fn (t)| 6 1 = ϕ(t) où ϕ est une fonction continue et intégrable sur 0, .
2 2
Z π/2
D’après le théorème de convergence dominée, lim Wn = f(t) dt = 0.
n→+∞ 0

3ème idée. L’équivalent de Wn obtenu en 11) fournit en particulier lim Wn = 0.

n→+∞

5) Premières valeurs.
Z π/2 Z π/2
π
W0 = dt = et W1 = sin t dt = 1.
0 2 0
π
W0 = et W1 = 1.
2
6) Relation de récurrence.
h πi
Soit n un entier naturel. Les deux fonctions t 7→ − cos t et t 7→ sinn+1 t sont de classe C1 sur 0, . On peut donc
2
effectuer une intégration par parties qui fournit

Z π/2 Z π/2
n+1 n+1
π/2
cos2 t sinn t dt

Wn+2 = sin t sin t dt = − cos t sin t 0
+ (n + 1)
0 0
Z π/2
= (n + 1) (1 − sin2 t) sinn t dt = (n + 1) (Wn − Wn+2 ) .
0

Donc,
n+1
∀n ∈ N, (n + 2)Wn+2 = (n + 1)Wn ou encore Wn+2 = Wn .
n+2
7) Calcul de Wn .
Soit p un entier naturel non nul.

 
2p
(2p − 1)(2p − 3) . . . 1 (2p)(2p − 1)(2p − 2)(2p − 3) . . . 1 π (2p)! π p π
W2p = W0 = 2
= 2p 2 = 2p ,
(2p)(2p − 2) . . . 2 ((2p)(2p − 2) . . . 2) 2 2 p! 2 2 2
ce qui reste vrai pour p = 0.
 
2p
(2p − 1)(2p − 3) . . . 1 π (2p)! π p π
∀p ∈ N∗ , W2p = × et ∀p ∈ N∗ , W2p = 2p 2 = 2p .
(2p)(2p − 2) . . . 2 2 2 p! 2 2 2
De même, si p un entier naturel non nul,

(2p)(2p − 2) . . . 2 ((2p)(2p − 2) . . . 2)2 22p p!2

W2p+1 = W1 = = ,
(2p + 1)(2p − 1) . . . 1 (2p + 1)(2p)(2p − 1) . . . 1 (2p + 1)!

2
ce qui reste vrai pour p = 0.

(2p)(2p − 2) . . . 2 22p p!2

∀p ∈ N∗ , W2p+1 = et ∀p ∈ N, W2p+1 = .
(2p + 1)(2p − 1) . . . 1 (2p + 1)!
8) Wn+1 est équivalent à Wn .
D ?après 3), la suite (Wn )n∈N est strictement positive et strictement décroissante. Donc, pour tout naturel n, on a
Wn+2 < Wn+1 < Wn . Après division par le réel strictement positif Wn , on obtient d’après 6)
n+1 Wn+2 Wn+1 Wn
= < < = 1.
n+2 Wn Wn Wn
Wn+1
Quand n tend vers l’infini, le théorème des gendarmes montre que lim = 1 ou encore
n→+∞ Wn
Wn+1 ∼ Wn .
n→+∞

9) Formule de Wallis.
W2p+1 W2p+1 (2 × 4 × . . . × (2p))2 2
D’après 8), lim = 1. D’autre part, d’après 7), = 2
. On obtient ainsi
p→+∞ W2p W2p (3 × 5 × . . . × (2p − 1)) (2p + 1)π
une première version de la formule de Wallis
1 √ 1 × 3 × . . . × (2p − 1)
√ = lim p .
π p→+∞ 2 × 4 × . . . × (2p)
ou encore, en élevant au carré et après simplification, on obtient (avec une formulation médiocre car les produits appa-
raissant au numérateur et au dénominateur sont divergents) :
4 3 × 3 × 5 × 5 × 7 × 7 × ...
= .
π 2 × 4 × 4 × 6 × 6 × ...
10) La suite ((n + 1)Wn Wn+1 )n∈N est constante.
D’après 6), pour tout entier naturel n, on a (n + 2)Wn+2 = (n + 1)Wn et donc (n + 2)Wn+1 Wn+2 = (n + 1)Wn Wn+1 .
π
Ainsi, la suite ((n + 1)Wn Wn+1 )n∈N est constante et donc, pour tout entier naturel n, (n + 1)Wn Wn+1 = W0 W1 = .
2
π
∀n ∈ N, (n + 1)Wn Wn+1 = .
2
11) Equivalent simple de Wn quand n tend vers +∞.
D’après 8), Wn+1 ∼ Wn et donc, d’après 10),
n→+∞

π 2
= (n + 1)Wn Wn+1 ∼ nWn .
2 n→+∞
r
p
2
π
Puisque Wn > 0, on en déduit que Wn = Wn ∼ .
n→+∞ 2n
r
π
Wn ∼ .
n→+∞ 2n
12) Série entière associée à Wn .
r
π X
Puisque Wn ∼ , la série entière Wn xn a un rayon de convergence égal à 1. Pour x ∈] − 1, 1[, on pose
n→+∞ 2n
+∞
X
f(x) = Wn x n .
n=0
h πi
Pour tout réel t de 0, , tout réel x tel que |x| < 1 et tout entier naturel n, on a
2
|(cosn t) xn | 6 |x|n .

Puisque la série géométrique de terme général |x|n converge pour x donné tel que |x| < 1, la série de fonctions
h π i de terme
n n
général t 7→ cos tx est normalement convergente et donc uniformément convergente sur le segment 0, . D’après le
2
théorème d’intégration terme à terme sur un segment, pour x ∈] − 1, 1[ fixé, on obtient

3
 +∞
X +∞
X Z π/2 ! Z π/2 +∞
X
! Z π/2
n n n n 1
f(x) = Wn x = cos t dt x = (x cos t) dt = dt.
0 0 0 1 − x cos t
n=0 n=0 n=0

t 2du
Soit x ∈] − 1, 1[ fixé. Calculons l’intégrale précédente en posant u = tan et donc dt = . On obtient
2 1 + u2
Z π/2 Z1 Z1
1 1 2du 2
f(x) = dt = 2 1 + u2
du = du
0 1 − x cos t 0 1−u 0 (1 − x) + u2 (1 + x)
1−x
1 + u2
 1
Z1 r r
2 1 2 1+x Arctan r u  2 1+x
= !2 du =  =√ Arctan .
1+x 0
r 1+x 1−x  1 − x  1 − x 2 1−x
1−x
u2 + 1+x 0
1+x

Donc,
+∞
X
r
n 2 1+x
∀x ∈] − 1, 1[, Wn x = √ Arctan .
1− x2 1−x
n=0
r
π
Maintenant, pour x = 1, puisque Wn ∼, la série de terme général Wn diverge ou encore f n’est pas défini en 1.
2n
n→+∞
D’autre part, pour x = −1, la suite Wn est positive et tend vers 0 en décroissant. Donc la série de terme général (−1)n Wn
converge en vertu du critère spécial aux séries alternées.
Montrons alors que la somme f est continue en −1.
Soit x ∈ [−1, 0]. La suite (−1)n Wn xn = Wn (−x)n = Wn |x|n est positive et tend vers 0 en décroissant (produit de
deux suites positives décroissantes). La série de terme général Wn xn est donc une série alternée. D’après une majoration
classique du reste à l’ordre n d’une série alternée, on a
 +∞ 
 X  
k

Wk x  6 Wn+1 xn+1  = Wn+1 |x|n+1 6 Wn+1 ,


 
k=n+1

et donc
 +∞  
 X 
k
sup  Wk x  , x ∈ [−1, 0] 6 Wn+1 → 0.

  n→+∞
k=n+1

La série entière de somme f converge uniformément vers f sur [−1, 0]. On en déduit que f est continue sur [−1, 0] et en
particulier que
r
2 1+x
f(−1) = lim f(x) = lim = √ Arctan .
x→−1 x→−1 1 − x2 1−x
x>−1 x>−1

Maintenant, quand x tend vers −1,

r r
2 1+x 2 1+x 2
√ Arctan ∼√ = ∼ 1.
1 − x2 1−x 1 − x2 1−x 1−x

On a montré que
+∞ r ! +∞
X 2 1+x X
n
∀x ∈] − 1, 1[, Wn x = √ Arctan et (−1)n Wn = 1.
1 − x2 1−x
n=0 n=0

13) Volume de la boule unité en dimension n.

Soit E un espace euclidien de dimension n strictement positive et Bn (R) la boule de centre O et de rayon strictement
positif R. Le volume de Bn (R) est
R R
Vn (R) = ··· dx1 . . . dxn .
x2 2
1 +...+xn 6R
2

4
On effectue déjà le changement
R de Rvariables x1 = Ry1n, ... , xn = Ryn . Le jacobien de ce changement de variables linéaire
est Rn et donc Vn (R) = Rn ··· dy1 ...dyn = R Vn (1).
y2 2
1 +...+yn 61

(∀R > 0), (∀n ∈ N∗ ), Vn (R) = Rn Vn (1).

Il reste à calculer Vn (1). Soit n un naturel supérieur ou égal à 2.

 
Z Z Z xn =1 Z Z Z xn =1 q
Vn (1) = ··· dx1 . . . dxn = ··· dx1 . . . dxn−1  dxn = Vn−1 ( 1 − x2n ) dxn
 

xn =−1 xn =−1
x2 2
1 +...+xn 61 x2 2 2
1 +...+xn−1 61−xn
Z1 p
= ( 1 − x2 )n−1 Vn−1 (1) dx = In−1 Vn−1 (1),
−1
Z 1 p n Z 1 p n
où In = 1 − x2 dx = 2 1 − x2 dx ou encore, en posant x = cos t :
−1 0

Z0 p n Z π/2
In = 2 1 − cos2 t (− sin t) dt = 2 sinn+1 t dt = 2Wn+1 .
π/2 0

Z1
En tenant compte de V1 (1) = dx = 2, on obtient
−1

V1 (1) = 2 et ∀n > 2, Vn (1) = 2Wn Vn−1 (1).

n
Y
Par suite, pour n > 2, Vn (1) = Vn (1)(2W2 )(2W3 )...(2Wn ) = 2n Wk (puisque W1 = 1). On retrouve en particulier
k=1

4 3
∀R > 0, V1 (R) = 2R, V2 (R) = πR2 et V3 (R) = πR .
3
π
Plus généralement, en tenant compte de l’égalité, valable pour tout entier naturel n, Wn Wn+1 = , on a pour p
2(n + 1)
entier naturel non nul donné
2p p p
2p
Y
2p
Y
2p
Y π πp
V2p (1) = 2 Wk = 2 W2k−1 W2k = 2 = .
2(2k) p!
k=1 k=1 k=1

πp πp R2p
∀p ∈ N∗ , V2p (1) = et ∀p ∈ N∗ , ∀R > 0, V2p (R) = .
p! p!

La formule donnant V2p+1 (R) est moins jolie et n’est pas donnée ici.