La solution de la série d’exercice N04

Exercice 1 :

Peut-on trouver des coefficients réels𝐴0 , 𝐴1 , 𝐴2 𝑒𝑡𝐴3 tels que

1
1 2
∫ 𝜑(𝑥) 𝑑𝑥 = 𝐴0 𝜑(0) + 𝐴1 𝜑 ( ) + 𝐴2 𝑓( ) + 𝐴3 𝜑(1)
0 3 3
Pour tout polynôme 𝜑de degré inférieur ou égal à 3 ? Si la réponse est oui, calculer ces
coefficients.

Oui on peut trouver les coefficients réels𝐴0 , 𝐴1 , 𝐴2 𝑒𝑡𝐴3 tels que

1
1 2
∫ 𝜑(𝑥) 𝑑𝑥 = 𝐴0 𝜑(0) + 𝐴1 𝜑 ( ) + 𝐴2 𝑓( ) + 𝐴3 𝜑(1)
0 3 3
Car On cherche une formule de quadrature à 4 points, elle doit donc être exacte sur l’ensemble
des polynômes de degré au plus 3.

Donc

Pour 𝑘 = 0; 𝜑(𝑥) = 𝑥 0 = 1, 𝐼(𝜑) = 𝐴0 + 𝐴1 + 𝐴2 + 𝐴3

1 2
Pour 𝑘 = 1; 𝜑(𝑥) = 𝑥1 = 𝑥, 𝐼(𝜑) = 𝐴1 + 3 𝐴2 + 𝐴3
3

1 4
Pour 𝑘 = 2; 𝜑(𝑥) = 𝑥 2 = , 𝐼(𝜑) = 𝐴1 + 9 𝐴2 + 𝐴3
9

1 8
Pour 𝑘 = 3; 𝜑(𝑥) = 𝑥 3 , 𝐼(𝜑) = 𝐴1 + 27 𝐴2 + 𝐴3
27

Pour calculer𝐴0 , 𝐴1 , 𝐴2 𝑒𝑡𝐴3 , il faut donc résoudre le système suivant

𝐴0 + 𝐴1 + 𝐴2 + 𝐴3 = 𝟏 … . (𝟏)
1 2 1
𝐴1 + 𝐴2 + 𝐴3 = … . . (2)
3 3 2
1 4 1
𝐴1 + 𝐴2 + 𝐴3 = … (3)
9 9 3
1 8 1
{ 27 𝐴1 + 𝐴2 + 𝐴3 = … (4)
27 4
2 2 1 3
(2) − (3) → 𝐴1 + 𝐴2 = → 𝐴1 + 𝐴2 = … . (5)
9 9 6 4
8 10 1 27
(2) − (4) → 𝐴1 + 𝐴2 = → 4 𝐴1 + 5𝐴2 = … . (6)
27 27 4 8
 21 3
(6) − 4(5) → 𝐴2 = → 𝐴2 = 8 On remplace la valeur de 𝐴2 dans l’équation (5)
8

3 3 3
Donc 𝐴1 + 8 = 4 → 𝐴1 = 8

D’après l’équation (2) on trouve

13 23 1 1
+ + 𝐴3 = → 𝐴3 =
38 38 2 8
1
AU fian d’prés (1) on trouve que 𝐴0 = 8

Finalement, la formule de quadrature est

1
1 3 1 3 2 1
∫ 𝜑(𝑥) 𝑑𝑥 = 𝜑(0) + 𝜑 ( ) + 𝑓( ) + 𝜑(1)
0 8 8 3 8 3 8

2. Appliquer la formule de quadrature ainsi définie à la fonction f :

1
𝑥 ↦ 𝑓 (𝑥) = 1/(2 − 𝑥) ; calculer la valeur exacte de l’intégrale∫0 𝑓(𝑥) 𝑑𝑥 puis donner un
majorant de l’erreur.

D’après la premier question on trouve que

1
1 3 1 3 2 1 11 33 33 1
𝐼(𝑓) = ∫ 𝑓(𝑥) 𝑑𝑥 ≃ 𝑓(0) + 𝑓 ( ) + 𝑓 ( ) + 𝑓(1) = + + + =
0 8 8 3 8 3 8 82 85 84 8
1
3 9 9 111
∫ 𝑓(𝑥) 𝑑𝑥 ≃ + = = 0.69375
0 16 8 20 160
1
- calculer la valeur exacte de l’intégrale∫0 𝑓(𝑥) 𝑑𝑥

1 1 1
I=∫0 𝑓(𝑥) 𝑑𝑥 = ∫0 𝑑𝑥 = −ln (2 − 𝑥)|10 = ln (2)
2−𝑥

111
L’erreur 𝐸(𝑓) = |𝐼 − 𝐼(𝑓)| = |ln(2) − 160| = 0.00060

𝑏
Exercice2 : On cherche à approcher l’intégrale 𝐼(𝑓) = ∫𝑎 𝑓(𝑥) 𝑑𝑥 𝑎 < 𝑏 par la formule des
trapèzes (cas n = 1 des formules de Newton-Cotes fermé). On approxime ainsi cette intégrale par
(1) (1) (1)
la formule de quadrature suivante(𝑏 − 𝑎) (𝑤0 𝑓(𝑎) + 𝑤1 𝑓(𝑏)) avec 𝑤𝑖 =
1 𝑏
∫ 𝑙 (𝑥) 𝑑𝑥
𝑏−𝑎 𝑎 𝑖

(1) (1)
et 𝑙𝑖 (𝑥) polynôme de Lagrange nécessaire pour l’interpolation. Calculer les poids 𝑤0 et 𝑤1
- Pour cette formule de quadrature, les nœuds sont les points a et b donc pour calculer les poids
(1) (1)
𝑤0 et 𝑤0 , il faut connaitre les les polynôme de Lagrange 𝑙0 (𝑥)et 𝑙1 (𝑥) nécessaire pour le
calcul du polynôme d’interpolation passant par les points : (x0, y0) = (a, f(a)) et (x1, y1)= (b, f(b)).
Pour rappel,
1
(𝑥 − 𝑥𝑗 )
𝑙𝑖 (𝑥) = ∏
𝑗=0
(𝑥𝑖 − 𝑥𝑗 )
𝑗≠𝑖

Dans notre cas, on a alors

(𝑥−𝑥1 ) (𝑥−𝑏) (𝑥−𝑥0 ) (𝑥−𝑎)

𝑙0 (𝑥) = (𝑥 = (𝑎−𝑏) ; 𝑙1 (𝑥) = (𝑥 = (𝑏−𝑎)
0 −𝑥1 ) 1 −𝑥0 )

D’où

𝑏 𝑏 (𝑥 𝑏
(1) 1 1 − 𝑏) 1 1 (𝑥 − 𝑏)2 1
𝑤0 = ∫ 𝑙0 (𝑥) 𝑑𝑥 = ∫ 𝑑𝑥 = | =
𝑏−𝑎 𝑎 𝑏 − 𝑎 𝑎 (𝑎 − 𝑏) (𝑏 − 𝑎) (𝑎 − 𝑏) 2 𝑎
2

𝑏 𝑏 (𝑥 𝑏
(1) 1 1 − 𝑎) 1 1 (𝑥 − 𝑎)2 1
𝑤1 = ∫ 𝑙01 (𝑥) 𝑑𝑥 = ∫ 𝑑𝑥 = | =
𝑏−𝑎 𝑎 𝑏 − 𝑎 𝑎 (𝑏 − 𝑎) (𝑏 − 𝑎) (𝑏 − 𝑎) 2 𝑎
2

La formule des trapèzes est donc

𝑏 (𝑏−𝑎)
∫𝑎 𝑓(𝑥) 𝑑𝑥 ≃ (𝑓(𝑎) + 𝑓(𝑏)) .
2

Exercice 3 :
𝑏
On cherche à approcher l’intégrale 𝐼(𝑓) = ∫𝑎 𝑓(𝑥) 𝑑𝑥 𝑎 < 𝑏 par la formule de Simpson (cas n =
𝑏−𝑎 𝑎+𝑏
2 des formules de Newton-Cotes fermé) 𝐼 = (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏))Montrer que la
6 2
formule de Simpson est de degré de précision 3.

- La formule de Simpson est une formule de quadrature de type interpolation à 3 points.

D’après l’un des théorèmes du cours, elle est donc exacte pour tous les polynômes de degré au
plus 2. Elle est donc de degré de précision au moins 2, ce que confirme les calculs suivants :
𝑏
Degré 0 : 𝑓(𝑥) = 1 → 𝐼(𝑓) = ∫𝑎 1 𝑑𝑥 = 𝑏 – 𝑎 ,

𝑏−𝑎 𝑎+𝑏 𝑏−𝑎

𝐼= 6
(𝑓(𝑎) + 4𝑓 ( 2 ) + 𝑓(𝑏)) = 6
(1 + 4 + 1) = 𝑏 – a

La formule de Simpson a au moins un degré de précision supérieur à 0.

 𝑏 𝑏2 −𝑎 2
Degré 1 : 𝑓(𝑥) = 𝑥 → 𝐼(𝑓) = ∫𝑎 𝑥 𝑑𝑥 = 2
,

𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 𝑎+𝑏 3(𝑏 − 𝑎)

𝐼= (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏)) = (𝑎 + 4 + 𝑏) = (𝑎 + 𝑏)
6 2 6 2 6
𝑏 2 − 𝑎2
=
2

La formule de Simpson a au moins un degré de précision supérieur à 1.

𝑏 𝑏3 −𝑎3
Degré 2 : 𝑓(𝑥) = 𝑥 2 → 𝐼(𝑓) = ∫𝑎 𝑥 2 𝑑𝑥 = 3
,

𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 2 (𝑎 + 𝑏)2

𝐼= (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏)) = (𝑎 + 4 + 𝑏2)
6 2 6 4
𝑏−𝑎 𝑏 3 − 𝑎3
= (2𝑎2 + 2𝑎𝑏 + 2𝑏 2 ) =
6 3

La formule de Simpson a au moins un degré de précision supérieur à 2. Maintenant pour le degré

supérieur, on a
𝑏 𝑏4 −𝑎4
Degré 3 : 𝑓(𝑥) = 𝑥 3 → 𝐼(𝑓) = ∫𝑎 𝑥 3 𝑑𝑥 = 4
,

𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 3 (𝑎 + 𝑏)3

𝐼= (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏)) = (𝑎 + 4 + 𝑏3)
6 2 6 8
𝑏−𝑎 3 1 3 3 2 3 1 𝑏−𝑎 3
= (𝑎 + 𝑎 + 𝑎 𝑏 + 𝑎𝑏 2 + 𝑏 3 + 𝑏 3 ) = (𝑎 + 𝑎𝑏 2 + 𝑎𝑏 2 + 𝑏 3 )
6 2 2 2 2 4
𝑏 4 − 𝑎4
=
4

La formule de Simpson a au moins un degré de précision supérieur à 3.

𝑏 𝑏5 −𝑎5
Degré 4 : 𝑓(𝑥) = 𝑥 4 → 𝐼(𝑓) = ∫𝑎 𝑥 4 𝑑𝑥 = 5
,

𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 (𝑎 + 𝑏)4

𝐼= (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏)) = (𝑎4 + 4 + 𝑏4)
6 2 6 16
𝑏−𝑎 1 3 1
= (𝑎4 + 𝑎4 + 𝑎3 𝑏 + 𝑎2 𝑏 2 + 𝑎𝑏 3 + 𝑏 4 + 𝑏 4 )
6 4 2 4
𝑏−𝑎
= (5𝑎4 + 4𝑎3 𝑏 + 6𝑎2 𝑏 2 + 4𝑎𝑏 3 + 5 𝑏 4 )
24
5 5 𝑎𝑏 4 1 2 3 1 3 2 𝑎4 𝑏 5 5
= 𝑏 − + 𝑎 𝑏 − 𝑎 𝑏 + − 𝑎
24 24 12 12 24 24

Le résultat obtenu par la formule de quadrature est différent de la valeur de l’intégrale pour
𝑓(𝑥) = 𝑥 4 . On en conclut que la formule de Simpson est de degré de précision 3.
 1
2
Exercice 4 : On veut approcher l’intégrale𝐼 = ∫02 𝑒 −𝑥 𝑑𝑥 .

1- Donner une valeur approchée de I en utilisant la méthode : trapèze -Simpson

Trapèze
1 1 2
ℎ 𝑏−𝑎 −0
𝐼 (𝑓) = 2 (𝑓(𝑎) + 𝑓(𝑏)) = (𝑓(𝑎) + 𝑓(𝑏)) = 2
(𝑒 −0 + 𝑒 −(2) ) = 0.4447
2 2

Simpson

1 1
(𝑏 − 𝑎 ) 𝑎+𝑏 −0 0+
𝐼 (𝑓) = (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏)) = 2 (𝑓(0) + 4𝑓 ( 2) + 𝑓 (1))
6 2 6 2 2

= 0.461371

2- Utiliser les méthodes, des trapèzes et de Simpson composites avec 3 sous-intervalles pour
approcher I.

Trapèze composite
2
𝑓(𝑎) + 𝑓(𝑏)
𝐼 (𝑓) = ℎ [ + ∑ 𝑓 (𝑥𝑖 )]
2
𝑖=1

1
𝑏−𝑎 −0 1 1 1 1
Avec ℎ = = 2 3 = 6 ; 𝑥0 = 𝑎 = 0; 𝑥1 = 𝑥0 + ℎ = 6 ; 𝑥2 = 𝑥1 + ℎ = 3 ; 𝑥3 = 𝑏 = 2.
𝑛

D’où

1
1 𝑓(0) + 𝑓 (2) 1 1
𝐼 (𝑓) = [ + 𝑓 ( ) + 𝑓 ( )] = 459474
6 2 6 3

Simpson composite
k=3 intervalles→ 𝑛 = 2𝑘 = 6 sous-intervalles
1
𝑏−𝑎 2−0 1 1
ℎ= = = ; 𝑥0 = 𝑎 = 0; 𝑥1 = 𝑥0 + ℎ = ; 𝑥 = 𝑥1 + ℎ
𝑛 6 12 12 2
1 1 1 5
= ; 𝑥3 = 𝑥2 + ℎ = ; 𝑥4 = 𝑥3 + ℎ = ; 𝑥5 = 𝑥4 + ℎ = ;𝑥 = 𝑏
6 4 3 12 6
1
=
2
 ℎ
𝐼 (𝑓) ≈ [𝑦 + 𝑦𝑛 + 2 ∑ 𝑦𝑖 + 4 ∑ 𝑦𝑖 ]
3 0
𝑖 𝑝𝑎𝑖𝑟 𝑖 𝑖𝑚𝑝𝑎𝑖𝑟
ℎ
= [𝑓(𝑥0 ) + 𝑓(𝑥6 ) + 2(𝑓(𝑥2 ) + 𝑓(𝑥4 )) + 4(𝑓(𝑥1 ) + 𝑓(𝑥3 ) + 𝑓(𝑥5 ))]
3
= 0.461282

3- Donner dans les deux cas un majorant de l’erreur.

2 2 2
𝑓(𝑥) = 𝑒 −𝑥 → 𝑓 ′ (𝑥) = −2𝑥𝑒 −𝑥 → 𝑓 (2) (𝑥) = (4𝑥 2 − 2)𝑒 −𝑥 → 𝑓 (3) (𝑥)
3 2 2
= −8𝑥 (𝑥 2 − ) 𝑒 −𝑥 → 𝑓 (4) (𝑥) = (16𝑥 4 − 4𝑥 2 + 12)𝑒 −𝑥
2
Et les majorants

𝑀2 = 𝑠𝑢𝑝[0;1/2] |𝑓 (2) (𝑥)| = |𝑓 (2) (0)| = 2; 𝑀4 = 𝑠𝑢𝑝[0;1/2] |𝑓 (4) (𝑥)| = |𝑓 (4) (0)| = 12

Ceci nous permet l’erreur dans la formule des trapèzes de majorer l’erreur commise avec les deux
méthodes composites

Erreur dans la formule des trapèzes

3
1
(𝑏 − 𝑎) (
3 − 0)
|𝐸(𝐼)| ≤ 𝑀2 =2 2 ≈ 0.0023
12𝑛2 12 ∗ 32
Erreur dans la formule des Simpson
5
1
(𝑏 − 𝑎)5 (2 − 0)
|𝐸(𝐼)| ≤ 𝑀4 = 12 ≈ 1.6 ∗ 10−6
180𝑛4 180 ∗ 64
4- En combien de sous-intervalles faut-il découper [0, 1/2] dans la méthode des trapèzes
composites pour avoir un majorant du même ordre de grandeur que celui obtenu avec la
méthode de Simpson pour trois sous-intervalles.

Soit nt le nombre nécessaire de sous-intervalles nécessaire pour obtenir une précision de

1.6*10−6 avec la méthode des trapèzes.

(𝑏 − 𝑎)3
|𝐸(𝐼)| ≤ 𝑀2 = 1.6 ∗ 10−6 → 𝑀2 (𝑏 − 𝑎)3 ≤ 12𝑛2 ∗ 1.6 ∗ 10−6
12𝑛2

3
1
𝑀2 (𝑏 − 𝑎) 3
𝑀2 (𝑏 − ∗ 𝑎)3 10−6 2
√ 2( − 0) ∗ 106
→ 𝑛2 ≥ −6
→𝑛≥√ = = 114.108
12 ∗ 1.6 ∗ 10 12 ∗ 1.6 ∗ 10−6 12 ∗ 1.6
Donc n=115

Exercice5 :
𝜋
1- Déterminer par la méthode des trapèzes puis par celle de Simpson∫0 𝑓(𝑥) 𝑑𝑥sur la base
2

du tableau suivant 𝜋 𝜋 𝜋
x 0 3𝜋
8 4 8 2

F(x) 0 0.382683 0.707107 0.923880 1

Trapèze composite
3
𝑓(𝑎) + 𝑓(𝑏)
𝐼 (𝑓) = ℎ [ + ∑ 𝑓 (𝑥𝑖 )]
2
𝑖=1
𝑏 − 𝑎 𝑓(𝑎) + 𝑓(𝑏)
= [ + 𝑓 (𝑥1 ) + 𝑓 (𝑥2 ) + 𝑓 (𝑥3 )]
𝑛 2
𝜋
𝜋 𝑓(0)+𝑓( 2 ) 𝜋 𝜋 3
= 8[ + 𝑓 ( 8 ) + 𝑓 (4 ) + 𝑓 (8)] = 0.987116
2

Simpson composite
𝑛 = 4 sous-intervalles
𝜋
𝑏−𝑎 2−0 𝜋
ℎ= = = .
𝑛 4 8
ℎ
𝐼 (𝑓) ≈ [𝑦 + 𝑦𝑛 + 2 ∑ 𝑦𝑖 + 4 ∑ 𝑦𝑖 ]
3 0
𝑖 𝑝𝑎𝑖𝑟 𝑖 𝑖𝑚𝑝𝑎𝑖𝑟
ℎ
= [𝑓(𝑥0 ) + 𝑓(𝑥4 ) + 2𝑓(𝑥2 ) + 4(𝑓(𝑥1 ) + 𝑓(𝑥3 ))]
3
𝜋
8 𝜋 𝜋 𝜋 3𝜋
= [𝑓(0) + 𝑓 (2 ) + 2𝑓 ( 4 ) + 4 (𝑓 ( 8 ) + 𝑓 ( 8 ))] = 1.000135
3

2- Ces points d’appui sont ceux donnant sin(x) comparé alors les résultats obtenus avec la
valeur exacte

Les points d’appui donnés dans cet exercice correspondent à la fonction sin x. Et 𝐼 =
𝜋
∫02 𝑓 (𝑥) 𝑑𝑥 = 1.On constate donc que l’approximation de I donnée par la méthode de Simpson
est meilleure que celle par les trapèzes, puisque |S − I| = 0.000135 et |T − I| = 0.012884
Exercice 6 : Considère la fonction f définie sur R par 𝑓 (𝑥) = 𝑒𝑥𝑝(𝑥 2 ). On souhaite calculer
1
une valeur approchée de𝐼 = ∫0 𝑓(𝑥) 𝑑𝑥par la méthode de Simpson.
On note A(n; f ) la valeur approchée fournie par la méthode de Simpson combinée avec n sous
intervalles.
(a) Déterminer une valeur de n aussi petite que possible assurant|𝐼 − 𝐴(𝑛, 𝑓 )| ≤ 10−3 .

D’après le cours, on a

(𝑏 − 𝑎)5
|𝐼 − 𝐴(𝑛, 𝑓 )| ≤ 𝑀4
180𝑛4
2
Calculons les dérivées de 𝑓 (𝑥) = 𝑒 𝑥 .
2 2 2
𝑜𝑛𝑎 𝑓′ (𝑥) = 2𝑥𝑒 𝑥 → 𝑓′′ (𝑥) = 2(2𝑥 2 + 1)𝑒 𝑥 → 𝑓 (3) (𝑥) = 2(4𝑥 3 + 6𝑥)𝑒 𝑥
2
→ 𝑓 (4) (𝑥) = 2(8𝑥 4 + 24𝑥 2 + 6)𝑒 𝑥 → 𝑀4 = 𝑚𝑎𝑥[0 ;1] |𝑓 (4) (𝑥)| = |𝑓 (4) (1)| = 76𝑒

Pour avoir la propriété demandée il suffit donc d’avoir

1 −3 4
76𝑒 3
76𝑒 ≤ 10 → 𝑛 ≥ 10 → 𝑛 = 6 →= 2𝑘 → 𝑘 = 3
180𝑛4 180
On vérifie immédiatement que n = 6 est la plus petite valeur de n satisfaisant la condition.

(b) D’après le cours on a

𝑘−1 𝑘−1
ℎ𝑛
𝐴(𝑛, 𝑓 ) = {𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓2𝑘 + 4 ∑ 𝑓2𝑘+1
3
𝑖=1 𝑖=0

𝑏−𝑎 𝑏−𝑎
Où a=0 , b=1et ℎ𝑛 = =
𝑛 2𝑘

𝑘−1 𝑘−1
ℎ𝑛
𝐴(𝑛, 𝑓̃ ) = {𝑓̃(𝑎) + 𝑓̃(𝑏) + 2 ∑ 𝑓̃2𝑘 + 4 ∑ 𝑓̃2𝑘+1
3
𝑖=1 𝑖=0

D’où l’on déduit

|𝐴(𝑛, 𝑓 ) − 𝐴𝑛 (𝑓̃ )| =
𝑘−1 𝑘−1
ℎ𝑛
= [(𝑓 − 𝑓̃)(𝑎) + (𝑓 − 𝑓̃)(𝑏) + 2 ∑(𝑓2𝑘+1 − 𝑓̃2𝑘 ) + 4 ∑(𝑓2𝑘+1 − 𝑓̃2𝑘+1 )]
3
𝑖=1 0=1

puis, en utilisant que la valeur absolue d’une somme est majorée par la somme des valeurs
absolues
 𝑘−1 𝑘−1
ℎ𝑛
|𝐴(𝑛, 𝑓 ) − 𝐴𝑛 (𝑓̃ )| ≤ [𝜀 + 𝜀 + 2 ∑ 𝜀 + 4 ∑ 𝜀 ]
3
𝑖=1 𝑖=0
𝑏−𝑎 𝑏−𝑎
≤ [𝜀 + 𝜀 + 2(𝑘 − 1 + 1 − 1)𝜀 + 4(𝑘 − 1 + 0 + 1)𝜀] = 6𝑘𝜀
3 ∗ 2𝑘 6𝑘
≤ (𝑏 − 𝑎)𝜀 = 𝜀
(c) On a

|𝐼 − 𝐴(𝑣, 𝑓̃ )| ≤ |𝐼 − 𝐴(𝑣, 𝑓 )| + |𝐴(𝑣, 𝑓 ) − 𝐴(𝑣, 𝑓̃ )| ≤ 10−3 + 𝜀

(d) Le résultat fourni par la calculatrice est 𝐴(𝑣, 𝑓̃ )𝑎𝑣𝑒𝑐 𝜀 = 10−12. Le résultat obtenu vérifie

Donc |𝐼 − 𝐴(𝑣, 𝑓̃ )| ≤ 10−3 + 10−12 La perte de précision de 10-12 due au calcul est
négligeable devant l’erreur de 10-3 due à la méthode.

Exercice 7 :
9
Utiliser la méthode de Simpson 3/8 avec 6 intervalles pour évaluer : ∫1 √𝑥𝑑𝑥. Comparer le
résultat avec la valeur exacte

D’après le cours on a
𝑏
3ℎ
∫ 𝑓(𝑥) 𝑑𝑥 ≈ (𝑓 + 3𝑓1 + 3𝑓2 + 2𝑓3 + 3𝑓4 + 3𝑓5 + 𝑓6 )
8 0
𝑎

𝑏−𝑎 4 7 11 19 23
ℎ= = ; 𝑥0 = 1; 𝑥𝑖 = 𝑥0 + 𝑛ℎ; 𝑥1 = ; 𝑥2 = ; 𝑥3 = 5; 𝑥4 = ; 𝑥5 = ;𝑥 = 9
𝑛 3 3 3 3 3 6

𝑏
1 7 11 19 23
∫ 𝑓(𝑥) 𝑑𝑥 ≈ (√1 + 3√ + 3√ + 2√5 + 3√ + 3√ + √9) = 17.3278663
𝑎 2 3 3 3 3
1

3 9 3
9 2 2
La valeur exacte de 𝐼 = ∫1 √𝑥𝑑𝑥 = 3 𝑥 2 | = 3 (92 − 1) = 17.3333333
1

𝑏
L’erreur absolue |𝐸| = |𝐼 − ∫𝑎 𝑓(𝑥)𝑑𝑥| = |17.3333333 − 17.3278663| = 0.54 ∗ 10−2

2
Exercice 8:On considère l’intégrale𝐼 = ∫1 ln (𝑥) 𝑑𝑥

1- La méthode des trapèzes composite à n+1 points (n sous-intervalles) pour calculer l’intégrale
d’une fonction f sur l’intervalle [a,b] s’écrit
 𝑛−1
𝑓(𝑥0 ) + 𝑓(𝑥𝑛 ) 𝑏−𝑎
𝐼 (𝑓) = ℎ [ + ∑ 𝑓 (𝑥𝑖 )] , 𝑎𝑣𝑒𝑐 ℎ =
2 𝑛
𝑖=1

2 1 5 6 7
Ici on a 𝑓 (𝑥) = ∫1 ln (𝑥) 𝑑𝑥 ; 𝑛 = 4 𝑑’𝑜ù ℎ = 4 ; 𝑥0 = 1; 𝑥1 = , 𝑥2 = ; 𝑥3 = ; 𝑥4 = 2.
4 4 4
et on obtient

𝑓(1) + 𝑓(2) 5 6 7
𝐼 (𝑓) = ℎ [ + 𝑓 ( ) + 𝑓 ( ) + 𝑓 ( )] = 0.3836995094
2 4 4 4

1-une primitive de ln (𝑥).est 𝐹(𝑥) = xln(𝑥) − 𝑥.La exacte est alors 𝐼 = [xln(𝑥) − 𝑥]12 =
2ln(2) − 1 = 0.386294361.

La valeur numérique obtenue est inférieure à celle exacte quelque soit le pas h choisi car la
fonction f est concave, ce qui signifie qu’une corde définie par deux points de la courbe y
=ln(x) sera toujours en-dessous de la courbe, donc l’aire sous les trapèzes sera inférieure à l’aire
exacte. Pour bien visualiser la construction considérons n= 2 :

2. L’erreur est majorée par

(𝑏 − 𝑎)3
|𝐸(𝐼)| ≤ 𝑀2
12𝑛2
1 1
𝑓 (𝑥) = ln(𝑥) → 𝑓 ′ (𝑥) = 𝑥 → 𝑓 ′′ (𝑥) = − 𝑥 2 → 𝑀2 = max|𝑓 ′′ (𝑥)| = |𝑓 ′′ (1)| = 1.
[1;2]

(𝑏 − 𝑎)3 −2 2
102 100
|𝐸(𝐼)| ≤ 𝑀2 2
≤ 10 → 𝑛 ≥ →𝑛=√ = 2.86
12𝑛 12 12

À partir de 3 sous-intervalles, l’erreur de quadrature est inférieure à 10−2 .

Exercice 9 :

D’après le cours l’erreur dans la formule des Simpson s’écrit :

(𝑏 − 𝑎)5
|𝐸(𝐼)| ≤ 𝑀4 ≤ 0.5 ∗ 10−8
180𝑛4
 𝑓 (𝑡) = 𝑡𝑒 −𝑡 → 𝑓 ′ (𝑥) = (1 − 𝑡)𝑒 −𝑡 → 𝑓 ′′ (𝑥) = (𝑡 − 2)𝑒 −𝑡 → 𝑓 (3) (𝑥) = (3 − 𝑡)𝑒 −𝑡
→ 𝑓 (4) (𝑥) = (𝑡 − 4)𝑒 −𝑡 → 𝑀4 = max|𝑓 (4) (𝑥)| = |𝑓 (4) (0)| = 4
[0;1]

D’où

(𝑏 − 𝑎)5 4 (𝑏 − 𝑎)5 4 4
|𝐸(𝐼)| ≤ 𝑀4 ≤ 0.5 ∗ 10 −8
→ 𝑛 ≥ √ ∗ 108 = √ ∗ 108 = 45.91
180𝑛4 180 ∗ 0.5 90

À partir de 46 sous-intervalles, l’erreur de quadrature est inférieure à 0.5 ∗ 10−8 .

Exercice 10 :
1 4 𝑤 2
Soit ∫−1 𝑔(𝑥) 𝑑𝑥 ≈ 3 𝑔 (− 2 ) + 3 𝑔(𝑤) … (1) ;où0 < 𝑤 ≤ 1 est à déterminer

1. Déterminer w pour que la formule de quadrature (1) soit exacte pour toute fonction g
polynomiale de degré𝑛 > 1 et donner la plus grande valeur den.

-Une formule de quadrature est dite de degré de précision n si elle est exacte pour

𝑥 𝑘 , 𝑘 = 0, . . . , 𝑛 Et non exacte pour𝑥 𝑛+1 .

1
Degré 0 : 𝑝(𝑥) = 𝑥0 = 1 → 𝐼(𝑝) = ∫−1 1 𝑑𝑥 = 2 ,

4 𝑤 2 4 2
𝐼 = 𝑝 (− ) + 𝑝(𝑤) = + = 2
3 2 3 3 3
La formule de quadrature a au moins un degré de précision supérieur à 0.
1
Degré 1 : 𝑝(𝑥) = 𝑥 → 𝐼(𝑝) = ∫−1 𝑥 𝑑𝑥 = 0 ,

4 𝑤 2 4 𝑤 2
𝐼 = 𝑝 (− ) + 𝑝(𝑤) = (− ) + 𝑤 = 0
3 2 3 3 2 3
La formule de quadrature a au moins un degré de précision supérieur à 1.
1 2
Degré 2 : 𝑝(𝑥) = 𝑥 2 → 𝐼(𝑝) = ∫−1 𝑥 2 𝑑𝑥 = 3 ,

4 𝑤 2 4 𝑤2 2
𝐼 = 𝑝 (− ) + 𝑝(𝑤) = ( ) + 𝑤2 = 𝑤2
3 2 3 3 4 3

2
si 𝑤 ≠ √3 la formule de quadrature est exacte pour les polynômes de degré au plus 1car →
2
𝐼(𝑝) ≠ 𝐼.La formule de quadrature a au moins un degré de précision supérieur à 2 si 𝑤 = √
3
 2
* si 𝑤 = √3

1
Degré 3 : 𝑝(𝑥) = 𝑥 3 → 𝐼(𝑝) = ∫−1 𝑥 3 𝑑𝑥 = 0 ,

3
2
4 𝑤 2 4 𝑤3 2 3 𝑤3 √3
𝐼 = 𝑝 (− ) + 𝑝(𝑤) = (− ) + 𝑤 = = ≠ 𝐼(𝑃) = 0
3 2 3 3 8 3 2 2

La formule de quadrature a un degré de précision 2.

Donc on peut résumer les résultats dans le tableau suivant :

𝑝(𝑥) = 1 4 𝑤 Degré d’exactitude

𝑘 ∫ 𝑝(𝑥)𝑑𝑥 𝑝 (− )
𝑥 −1 3 2
k 2
+ 𝑝(𝑤)
3
0 1 2 2 Au moins 0

1 x 0 0 Au moins 1

2 𝑥2 2/3 𝑤2 2 2
1 si 𝑤 ≠ √3 Au moins 2 si 𝑤 = √3

3 𝑥3 0 𝑤 3 /2 2

2 2
Si 𝑤 ≠ √3la formule de quadrature est exacte pour les polynômes de degré au plus 1, si 𝑤 = √3
la formule est exacte pour les polynômes de degré au plus 2.
(𝑥 −𝑥 )
2. Par le changement de variable 𝑦 = 𝑥𝑖 + (𝑥 + 1) 𝑖+12 𝑖 on déduit la formule de quadrature
𝑥𝑖+1
𝑥𝑖+1 − 𝑥𝑖 1 𝑥𝑖+1 − 𝑥𝑖
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓 (𝑥𝑖 + (𝑥 + 1) ) 𝑑𝑥
𝑥𝑖 2 −1 2

𝑥𝑖+1 − 𝑥𝑖 2 𝑥𝑖+1 − 𝑥𝑖 1 𝑥𝑖+1 − 𝑥𝑖

≈ [𝑓 (𝑥𝑖 + (1 + √ ) ) + 𝑓 (𝑥𝑖 + (1 − √ ) )]
3 3 2 6 2
𝑏−𝑎
3. Si 𝐻 = 𝑥𝑖+1 − 𝑥𝑖 = 𝑛 an (i.e. si on considère une subdivision de l’intervalle [a;b]
équirépartie) alors on trouve la formule de quadrature composite (i.e. sur n sous-intervalles et à
2n points)
 𝑥𝑖+1 𝑛−1
𝐻 2 1
∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ [𝑓 (𝑥𝑖 + 𝐻 (1 + √ )) + 2𝑓 (𝑥𝑖 + 𝐻 (1 − √ ))]
𝑥𝑖 3 3 6
𝑖=0

𝐻 2 1
∑𝑛−1
𝑖=0 [𝑓 (𝑎 + 𝐻 (𝑖 + 1 + √3)) + 2𝑓 (𝑎 + 𝐻 (𝑖 + 1 − √6))]
3

Cette formule de quadrature est stable puisque tous les coefficients sont positifs

Exercice 14 :
Soit f une fonction C1(R,R).
1 2 1 1
1. On considère l’approximation ∫−1 𝑓(𝑥) 𝑑𝑥 ≈ 3 (2𝑓 (− 2) − 𝑓(0) + 2𝑓(2))
Quel est le degré de précision de cette formule de quadrature ?
-Une formule de quadrature est dite de degré de précision n si elle est exacte pour

𝑥 𝑘 , 𝑘 = 0, . . . , 𝑛 Et non exacte pour𝑥 𝑛+1 .

1
Degré 0 : 𝑝(𝑥) = 𝑥0 = 1 → 𝐼(𝑝) = ∫−1 1 𝑑𝑥 = 2 ,

2 1 1 2
𝐼 = (2𝑝 (− ) − 𝑝(0) + 2𝑝( )) = (2 − 1 + 2) = 2
3 2 2 3

La formule de quadrature a au moins un degré de précision supérieur à 0.

1
Degré 1 : 𝑝(𝑥) = 𝑥 → 𝐼(𝑝) = ∫−1 𝑥 𝑑𝑥 = 0 ,

2 1 1 2 1 1
𝐼 = (2𝑝 (− ) − 𝑝(0) + 2𝑝( )) = (−2 + 0 + 2 ) = 0
3 2 2 3 2 2

La formule de quadrature a au moins un degré de précision supérieur à 1.

1 2
Degré 2 : 𝑝(𝑥) = 𝑥 2 → 𝐼(𝑝) = ∫−1 𝑥 2 𝑑𝑥 = 3 ,

2 1 1 2 1 1 2
𝐼 = (2𝑝 (− ) − 𝑝(0) + 2𝑝( )) = (2 + 0 + 2 ) =
3 2 2 3 4 4 3

La formule de quadrature a au moins un degré de précision supérieur à 2.

1
Degré 3 : 𝑝(𝑥) = 𝑥 3 → 𝐼(𝑝) = ∫−1 𝑥 3 𝑑𝑥 = 0 ,

2 1 1 2 1 1
𝐼 = (2𝑝 (− ) − 𝑝(0) + 2𝑝( )) = (−2 + 0 + 2 ) = 0
3 2 2 3 8 8
La formule de quadrature a au moins un degré de précision supérieur à 3.
 1 2
Degré 4 : 𝑝(𝑥) = 𝑥 4 → 𝐼(𝑝) = ∫−1 𝑥 4 𝑑𝑥 = 5 ,

2 1 1 2 1 1 1 2
𝐼= (2𝑝 (− ) − 𝑝(0) + 2𝑝( )) = (2 + 0 + 2 ) = ≠ 𝐼(𝑝) =
3 2 2 3 16 16 6 5
La formule de quadrature et donc exacte de degré 3.

Donc on peut résumer les dans le tableau suivant :

1
k 𝑝(𝑥) = 2 1 1 Degré d’exactitude
𝑘 ∫ 𝑝(𝑥)𝑑𝑥 (2𝑝 (− ) − 𝑝(0) + 2𝑓( ))
𝑥 −1 3 2 2

0 1 2 2 au moins 0

1 x 0 0 au moins 1

2 𝑥2 2/3 2/3 au moins 2

3 𝑥3 0 0 au moins 3

𝑥4 2/5 1/6 3

2. Soit 𝑥 = 𝑚𝑡 + 𝑞, alors
𝑥𝑖+1 1
𝑥 = 𝑚𝑡 + 𝑞
∫ 𝑓(𝑥)𝑑𝑥 = 𝑚 ∫ 𝑓(𝑚𝑡 + 𝑞)𝑑𝑡 𝑎𝑣𝑒𝑐 {𝑥 𝑖 = 𝑚𝑡 + 𝑞
𝑥𝑖 −1 𝑖+1

𝑥 −𝑥 𝑏−𝑎
d’où le changement de variable 𝑥 = 𝑥𝑖 + (𝑡 + 1) 𝑖+12 𝑖 = 𝑛 . On déduit la formule de
quadrature (exacte sur l’espace des polynômes de degré au plus 3)
𝑥𝑖+1
𝑥𝑖+1 − 𝑥𝑖 1 𝑥𝑖+1 − 𝑥𝑖
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓 (𝑥𝑖 + (𝑡 + 1) ) 𝑑𝑡
𝑥𝑖 2 −1 2
𝑥𝑖+1 + 𝑥𝑖 𝑥𝑖+1 + 𝑥𝑖
𝑥𝑖+1 − 𝑥𝑖 𝑥𝑖+1 + 𝑥𝑖 +𝑥𝑖+1
≈ [2𝑓 (𝑥𝑖 + 2 )−𝑓( ) + 2𝑓 ( 2 )]
3 2 2 2

𝑏−𝑎
Soit ℎ = 𝑥𝑖+1 − 𝑥𝑖 = La formule précédente se réécrit
𝑛
𝑥𝑖+1
ℎ ℎ ℎ 3ℎ
∫ 𝑓(𝑥)𝑑𝑥 ≈ [2𝑓 (𝑥𝑖 + ) − 𝑓 (𝑥𝑖 + ) + 2𝑓 (𝑥𝑖 + )]
𝑥𝑖 3 4 3 4

et la formule de quadrature composite déduite de cette approximation est

 𝑏 𝑛−1 𝑥𝑖+1 𝑛−1
ℎ ℎ ℎ 3ℎ
∫ 𝑓(𝑥)𝑑𝑥 = ∑ ∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ [2𝑓 (𝑥𝑖 + ) − 𝑓 (𝑥𝑖 + ) + 2𝑓 (𝑥𝑖 + )]
𝑎 3 4 3 4
𝑖=0 𝑥𝑖 𝑖=0