Université de la Manouba

Institut Supérieur des Arts Multimédias

Correction série N°3

A.U 2023-2024 2 IM.

Exercice 1 :
Soit f une fonction C ∞ (R, R).
1. On considère l’approximation
Z 1
2 1 1
f (x)dx ≈ (2f (− ) − f (0) + 2f ( ))
−1 3 2 2

Quel est le degré d’exactitude de cette formule de quadrature ?

b−a
2. On se donne les points {xi }i=n
i=0 de subdivision de l’intervalle [a, b] : xi = a + ih avec h = n .
À l’aide d’un changement de variable affine, en déduire une formule de quadrature pour
l’intégrale Z xi+1
f (x)dx
xi
En tirer une formule de quadrature composite pour l’intégrale
Z b
f (x)dx
a

3. Écrire l’algorithme pour approcher

Z b
f (x)dx
a

Réponse :
R1 2 1
1. — pour k = 0, f (X) = X k = X 0 = 1, −1 1dx = [x]1−1 = 1 − 1(−1) = 2 et (2f (− ) −
3 2
1 2
f (0) + 2f ( )) = (2 ∗ 1 − 1 + 2 ∗ 1) = 2 ⇒ exacte pour k = 0.
2 3
R1 2 2 1
— pour k = 1, f (X) = X k = X, −1 xdx = [ x2 ]1−1 = 12 − 21 = 0 et (2f (− ) − f (0) +
3 2
1 2 −1 1
2f ( )) = (2 ∗ − 0 + 2 ∗ ) = 0 ⇒ exacte pour k = 1.
2 3 2 2
R1 3 2 1
— pour k = 2, f (X) = X = X 2 , −1 x2 dx = [ x3 ]1−1 = 13 + 13 = 32 et (2f (− ) − f (0) +
k
3 2
1 2 1 2 1 2
2f ( )) = (2 ∗ − 0 + 2 ∗ ) = 3 ⇒ exacte pour k = 2.
2 3 4 4
R1 4 2 1
— pour k = 3, f (X) = X k = X 3 , −1 x3 dx = [ x4 ]1−1 = 14 − 41 = 0 et (2f (− ) − f (0) +
3 2
1 2 −1 3 1
2f ( )) = (2 ∗ − 0 + 2 ∗ ) = 0 ⇒ exacte pour k = 3.
2 3 8 8
R1 5 2 1
— pour k = 4, f (X) = X = X , −1 x4 dx = [ x5 ]1−1 = 15 + 15 = 52 et (2f (− ) − f (0) +
k 4
3 2
1 2 1 4 1 1 2 1
2f ( )) = (2 ∗ − 0 + 2 ∗ ) = 6 , 5 ̸= 6 ⇒ n’est pas exacte pour k = 4.
2 3 16 16
donc le degré d’exactitude est 3.

1
 2.
xi+1 1
xi+1 − xi xi+1 − xi
Z Z
xi+1 + xi
f (x)dx = f( t+ )dt
xi 2 −1 2 2

xi = a + i b−a b−a
n , xi+1 = a + (i + 1) n , xi+1 − xi = b−a
n =h
xi+1 −xi
2 = h2
xi+1 +xi
2 = 2a+(2i+1)h
2 = a + (i + 12 )h
Z xi+1 Z 1
h h 1
⇒ f (x)dx = f ( t + a + (i + )h)dt
xi 2 −1 2 2
 
h2 h 1 1 h 1 h1 1
≃ 2f ( (− ) + a + (i + )h) − f ( (0) + a + (i + )h) + 2f ( + a + (i + )h)
23 2 2 2 2 2 22 2
 
h 1 1 3
≃ 2f (a + (i + )h) − f (a + (i + )h) + 2f (a + (i + )h)
3 4 2 4
Z b Z x1 Z x2 Z xn n−1
X Z xi+1
f (x)dx = f (x)dx + f (x)dx + · · · + f (x)dx = f (x)dx
a x0 x1 xn−1 i=0 xi

n−1
X h 1 1 3
≃ (2f (a + (i + )h) − f (a + (i + )h) + 2f (a + (i + )h))
3 4 2 4
i=0
n−1 n−1 n−1
" #
h X 1 X 1 X 3
≃ 2 f (a + (i + )h) − f (a + (i + )h) + 2 f (a + (i + )h)
3 4 2 4
i=0 i=0 i=0

3. Algorithme :
Algorithm 1: Algorithme Calcul intégrale
Result: In
a, b et n;
S1 = 0 ; S2 = 0 ; S3 = 0 ; h = b−an ;
for i = 0 : n − 1 do
S1 = S1 + f (a + (i + 14 )h);
S2 = S2 + f (a + (i + 12 )h);
S3 = S3 + f (a + (i + 34 )h);
end
In = h3 (2S1 − S2 + 2S3 );

Exercice 2 :
On considère l’intégrale Z 2
1
I= dx
1 x
1. Calculer la valeur exacte de I.
2. Évaluer numériquement cette intégrale par la méthode des trapèzes avec m = 3 sous-
intervalles.
3. Pourquoi la valeur numérique obtenue à la question précédente est-elle supérieure à ln(2) ?
Est-ce vrai quelque soit m ? Justifier la réponse. (On pourra s’aider par un dessin.)
4. Quel nombre de sous-intervalles m faut-il choisir pour avoir une erreur inférieure à 10−4 ?
On rappelle que l’erreur de quadrature associée s’écrit, si f ∈ C 2 ([a, b]),

(b − a)3 ′′
|Em (f )| = |f (α)|, α ∈ [a, b]
12m2

2
Réponse :
R2 1 2
1. I = 1 x dx = [ln(x)]1 = ln(2) − ln(1) = ln(2) ≃ 0.693
h Pm−1 i
2. I ≃ b−a
2m f (a) + f (b) + 2 i=1 f (x i , a = 1, b = 2 et m = 3
)
h i
⇒ I ≃ 16 f (1) + f (2) + 2 2i=1 f (xi ) , x0 = 1, x1 = 43 , x2 = 53 , x3 = 2
P

⇒ I ≃ 61 11 + 12 + 2( 34 + 35 ) = 61 1 + 12 + 27 7
   
10 = 10 = 0.7
3. La fonction f est convexe sur [1, 2] donc surface du trapèze est plus grande que celle entre
la courbe est l’axe des x donc l’intégrale.
3
De plus l’erreur Em (f ) = − (b−a)
12m2
f ′′ (α), α ∈ [1, 2] et f ′′ (x) = 2
x3
donc Em (f ) < 0 d’où la
valeur numérique de l’intégrale est supérieure à ln(2).

f (x)

x10 x1 x2 x23

4.
(2 − 1)3 ′′
|Em (f )| = |f (α)| ⩽ 10−4 , α ∈ [1, 2]
12m2
1 2 1 −4 2 104 100
⇒ 2
| 3
| ⩽ 2
⩽ 10 ⇒ m ⩾ ⇒ m ⩾ √ ≃ 40.825
12m α 6m 6 6
d’où m = 41

Exercice 3 :
Soit f une fonction C ∞ (R, R). On se donne les points {xi }i=2n
i=0 de subdivision de l’intervalle [a, b] :
xi = a + ih avec h = b−a2n . Le but de l’exercice est de trouver une formule de quadrature à 2n + 1
points basée sur la formule de Simpson pour approcher
Z b
f (x)dx, (1)
a

On propose dans un premier temps (question 1 à 5) de construire la formule de quadrature à 3

points de Simpson : Z 1
g(x)dx ≈ αg(−1) + βg(0) + αg(1), (2)
−1
où les réels α et β sont à déterminer.
1. Sous quelle condition (portant sur α et β) la formule de quadrature (2) est exacte pour une
fonction g constante ?
2. Sous quelle condition (portant sur α et β) la formule de quadrature (2) est exacte pour une
fonction g polynomiale de degré au plus 2 ?
3. En déduire le choix de α et β rendant la formule de quadrature (2) exacte pour une fonction
g polynomiale de degré au plus 2.

3
 4. La formule de quadrature est-elle exacte pour tout polynôme de degré 3 ? La formule de
quadrature est-elle exacte pour tout polynôme de degré 4 ?
5. À l’aide d’un changement de variable affine, en déduire une formule de quadrature exacte
sur l’espace des polynôme de degré au plus 3 pour l’intégrale suivante :
Z x2i+2
f (x)dx.
x2i

6. En déduire une formule de quadrature à 2n points, notée F , pour le calcul approché de (1).
Cette formule de quadrature est-elle stable ?
7. Écrire l’algorithme du calcul de F .
8. Soit x un élément de [x2i , x2i+2 ]. Écrire une formule de Taylor f (x) = P (x) + R(x) à l’ordre
3 pour f en x, avec P ∈ P3 .
Majorer R sur [x2i , x2i+2 ] en fonction de h.
9. En déduire une estimation d’erreur entre (1) et F .

Réponse :
R1
1
1. pour g(x) = c on a −1 cdx = c[x]−1 = 2c, αg(−1) + βg(0) + αg(1) = c(2α + β) ⇒ 2α + β = 2
R1 2
2. — pour g(x) = X on a −1 xdx = [ x2 ]1−1 = 0, αg(−1) + βg(0) + αg(1) = −α + β ∗ 0 + α = 0.
R1 3
— pour g(x) = X 2 on a −1 x2 dx = [ x3 ]1−1 = 23 , αg(−1) + βg(0) + αg(1) = α + β ∗ 0 + α =
2α ⇒ 2α = 32 ⇒ α = 13 .
α = 31
 
2α + β = 2
3. ⇒
α = 13 β = 34
R1 4
4. — pour g(x) = X 3 on a −1 x3 dx = [ x4 ]1−1 = 0, αg(−1) + βg(0) + αg(1) = −α + β ∗ 0 + α =
0 ⇒ la formule est exacte
R 1 pour k = 3.
5
— pour g(x) = X on a −1 x dx = [ x5 ]1−1 = 25 , αg(−1) + βg(0) + αg(1) = α + β ∗ 0 + α =
4 4

2α = 23 ̸= 25 ⇒ la formule n’est pas exacte pour k = 4.

5.
x2i+2 1
x2i+2 − x2i x2i+2 − x2i
Z Z
x2i+2 + x2i
f (x)dx = f( t+ )dt
x2i 2 −1 2 2

x2i = a + 2i b−a b−a

2n , x2i+2 = a + (2i + 2) 2n ,
x2i+2 −x2i b−a
2 = 2n = h
x2i+2 +x2i 2a+(4i+2)h
2 = 2 = a + (2i + 1)h
Z x2i+2 Z 1
⇒ f (x)dx = h f (ht + a + (2i + 1)h)dt
x2i −1
 
1 4 1
≃ h f (h ∗ (−1) + a + (2i + 1)h) + f (h ∗ 0 + a + (2i + 1)h) + f (h ∗ 1 + a + (2i + 1)h)
3 3 3
h
≃ [f (a + 2ih) + 4f (a + (2i + 1)h) + f (a + (2i + 2)h)]
3
6.
Z b Z x2 Z x4 Z x2n n−1
X Z x2i+2
f (x)dx = f (x)dx + f (x)dx + · · · + f (x)dx = f (x)dx
a x0 x2 x2n−2 i=0 x2i

n−1
X h
≃ [f (a + 2ih) + 4f (a + (2i + 1)h) + f (a + (2i + 2)h)]
3
i=0

4
 n−1 n−1 n−1
!
h X X X
≃ f (a + 2ih) + 4 f (a + (2i + 1)h) + f (a + (2i + 2)h)
3
i=0 i=0 i=0
n−1 n−1 n
!
h X X X
≃ f (a + 2ih) + 4 f (a + (2i + 1)h) + f (a + 2ih)
3
i=0 i=0 i=1
n−1 n−1
!
h X X
≃ f (a) + f (b) + 2 f (a + 2ih) + 4 f (a + (2i + 1)h) , (F )
3
i=1 i=0

les poids de la formule de quadrature (F ) sont positifs donc la formule est stable.
7. Algorithme :
Algorithm 2: Algorithme Calcul intégrale
Result: In
a, b et n;
h = b−a
2n ;
S1 = 0 ; S2 = f (a + h);
for i = 1 : n − 1 do
S1 = S1 + f (a + 2ih);
S2 = S2 + f (a + (2i + 1)h);
end
In = h3 (f (a) + f (b) + 2S1 + 4S2 );

8. Soit x un élément de [x2i , x2i+2 ]. Une formule de Taylor à l’ordre 3 pour f en x s’écrit :
f (x) = P (x) + R(x)

(x − x2i )2 ′′ (x − x2i )3 ′′′ (x − x2i )4 (4)

f (x) = f (x2i ) + (x − x2i )f ′ (x2i ) + f (x2i ) + f (x2i ) + f (α),
| 2!
{z 3! } | 4! {z }
P (x) R(x)

α ∈ [x2i , x2i+2 ]
x2i+2 −x2i
On peut majorer R sur [x2i , x2i+2 ] en fonction de h = 2 :

(x − x2i )4 (4) (2h)4 2h4

|R(X)| = f (α) ⩽ maxα∈[x2i ,x2i+2 ] f (4) (α) = maxα∈[x2i ,x2i+2 ] f (4) (α)
4! 4! 3

9.
Z b n−1
X Z x2i+2
f (x)dx − F = (P (x) + R(x))dx − F
a i=0 x2i

n−1
X Z x2i+2 n−1
X Z x2i+2 n−1
X Z x2i+2
⩽ (P (x) − F )dx + R(x)dx ⩽ |R(x)| dx
i=0 x2i i=0 x2i i=0 x2i

n−1
X 2h4 2h4 4nh5
⩽ 2h maxα∈[x2i ,x2i+2 ] f (4) (α) ⩽ 2hn maxα∈[a,b] f (4) (α) = maxα∈[a,b] f (4) (α)
3 3 3
i=0