Exercice
Soit la fonction :
1
𝑓(𝑥) =
1 + 𝑥2
1. Déterminer le polynôme d’interpolation de Lagrange 𝑃𝑛 (𝑥) pour les points d’appui
d’abscisses −2, −1, 0, 1, 2.
2. Discuter l’erreur d’interpolation.
2 𝑑𝑥
3. Déterminer ∫−2 1+𝑥 2
2
4. Déterminer ∫−2 𝑃𝑛 (𝑥)𝑑𝑥
5. Déterminer par la méthode des rectangles, des trapèzes puis par celle de Simpson
2
∫−2 𝑓(𝑥)𝑑𝑥 avec les points d’appui −2, −1, 0, 1, 2. Comparer les résultats.
Corrigé
1. Déterminer le polynôme d’interpolation de Lagrange pour les points d’appui
d’abscisses −2, −1, 0, 1, 2.
Ici le nombre de points d’appui étant 5, alors le polynôme de Lagrange est un polynôme de
degré 𝑛 = 4.
Donc
𝑥0 𝑥1 𝑥2 𝑥3 𝑥4
𝑥𝑖 −2 −1 0 1 2
𝑓(𝑥𝑖 ) 1 1 1 1 1
5 2 2 5
𝑓(𝑥0 ) 𝑓(𝑥1 ) 𝑓(𝑥2 ) 𝑓(𝑥3 ) 𝑓(𝑥4 )
Le polynôme de Lagrange de degré 4 est de la forme :
𝒏 𝒏
𝒙 − 𝒙𝒋
𝑷𝒏 (𝒙) = ∑ 𝒇(𝒙𝒊 ) 𝑳𝒊 (𝒙) 𝒂𝒗𝒆𝒄 𝑳𝒊 (𝒙) = ∏
𝒙𝒊 − 𝒙𝒋
𝒊=𝟎 𝒋=𝟎
𝒊≠𝒋
4 4
𝑥 − 𝑥𝑗
𝑃4 (𝑥) = ∑ 𝑓(𝑥𝑖 ) 𝐿𝑖 (𝑥) 𝑎𝑣𝑒𝑐 𝐿𝑖 (𝑥) = ∏
𝑥𝑖 − 𝑥𝑗
𝑖=0 𝑗=0
𝑖≠𝑗
𝑛=4
𝑖 = 0,1,2,3,4 𝑒𝑡 𝑗 = 0,1,2,3,4
𝑖 = 0 𝑒𝑡 𝑗 = 1,2,3,4
(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )(𝑥 − 𝑥3 )(𝑥 − 𝑥4 )
𝐿0 (𝑥) =
(𝑥0 − 𝑥1 )(𝑥0 − 𝑥2 )(𝑥0 − 𝑥3 )(𝑥0 − 𝑥4 )
(𝑥 + 1)(𝑥 − 0)(𝑥 − 1)(𝑥 − 2) 𝑥(𝑥 + 1)(𝑥 − 1)(𝑥 − 2)
𝐿0 (𝑥) = =
(−2 + 1)(−2 − 0)(−2 − 1)(−2 − 2) (−1)(−2)(−3)(−4)
𝟏
𝑳𝟎 (𝒙) = 𝒙(𝒙 + 𝟏)(𝒙 − 𝟏)(𝒙 − 𝟐)
𝟐𝟒
𝑖 = 1 𝑒𝑡 𝑗 = 0,2,3,4
(𝑥 − 𝑥0 )(𝑥 − 𝑥2 )(𝑥 − 𝑥3 )(𝑥 − 𝑥4 )
𝐿1 (𝑥) =
(𝑥1 − 𝑥0 )(𝑥1 − 𝑥2 )(𝑥1 − 𝑥3 )(𝑥1 − 𝑥4 )
(𝑥 + 2)(𝑥 − 0)(𝑥 − 1)(𝑥 − 2) 𝑥(𝑥 + 2)(𝑥 − 1)(𝑥 − 2)
𝐿1 (𝑥) = =
(−1 + 2)(−1 − 0)(−1 − 1)(−1 − 2) (1)(−1)(−2)(−3)
𝟏
𝑳𝟏 (𝒙) = − 𝒙(𝒙 + 𝟐)(𝒙 − 𝟏)(𝒙 − 𝟐)
𝟔
𝑖 = 2 𝑒𝑡 𝑗 = 0,1,3,4
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥3 )(𝑥 − 𝑥4 )
𝐿2 (𝑥) =
(𝑥2 − 𝑥0 )(𝑥2 − 𝑥1 )(𝑥2 − 𝑥3 )(𝑥2 − 𝑥4 )
(𝑥 + 2)(𝑥 + 1)(𝑥 − 1)(𝑥 − 2) (𝑥 + 2)(𝑥 + 1)(𝑥 − 1)(𝑥 − 2)
𝐿2 (𝑥) = =
(0 + 2)(0 + 1)(0 − 1)(0 − 2) (2)(1)(−1)(−2)
𝟏
𝑳𝟐 (𝒙) = (𝒙 + 𝟐)(𝒙 + 𝟏)(𝒙 − 𝟏)(𝒙 − 𝟐)
𝟒
𝑖 = 3 𝑒𝑡 𝑗 = 0,1,2,4
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )(𝑥 − 𝑥4 )
𝐿3 (𝑥) =
(𝑥3 − 𝑥0 )(𝑥3 − 𝑥1 )(𝑥3 − 𝑥2 )(𝑥3 − 𝑥4 )
(𝑥 + 2)(𝑥 + 1)(𝑥 − 0)(𝑥 − 2) 𝑥(𝑥 + 2)(𝑥 + 1)(𝑥 − 2)
𝐿3 (𝑥) = =
(1 + 2)(1 + 1)(1 − 0)(1 − 2) (3)(2)(1)(−1)
𝟏
𝑳𝟑 (𝒙) = − 𝒙(𝒙 + 𝟐)(𝒙 + 𝟏)(𝒙 − 𝟐)
𝟔
𝑖 = 4 𝑒𝑡 𝑗 = 0,1,2,3
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )(𝑥 − 𝑥3 )
𝐿4 (𝑥) =
(𝑥4 − 𝑥0 )(𝑥4 − 𝑥1 )(𝑥4 − 𝑥2 )(𝑥4 − 𝑥3 )
(𝑥 + 2)(𝑥 + 1)(𝑥 − 0)(𝑥 − 1) 𝑥(𝑥 + 2)(𝑥 + 1)(𝑥 − 1)
𝐿4 (𝑥) = =
(2 + 2)(2 + 1)(2 − 0)(2 − 1) (4)(3)(2)(1)
𝟏
𝑳𝟒 (𝒙) = 𝒙(𝒙 + 𝟐)(𝒙 + 𝟏)(𝒙 − 𝟏)
𝟐𝟒
Finalement
𝒏
𝑷𝒏 (𝒙) = ∑ 𝒇(𝒙𝒊 ) 𝑳𝒊 (𝒙)
𝒊=𝟎
𝑃4 (𝑥) = 𝑓(𝑥0 )𝐿0 (𝑥) + 𝑓(𝑥1 )𝐿1 (𝑥) + 𝑓(𝑥2 )𝐿2 (𝑥) + 𝑓(𝑥3 )𝐿3 (𝑥) + 𝑓(𝑥4 )𝐿4 (𝑥)
1 𝟏 1 𝟏
𝑃4 (𝑥) = . 𝒙(𝒙 + 𝟏)(𝒙 − 𝟏)(𝒙 − 𝟐) − . 𝒙(𝒙 + 𝟐)(𝒙 − 𝟏)(𝒙 − 𝟐)
5 𝟐𝟒 2 𝟔
𝟏 1𝟏
+ 1. (𝒙 + 𝟐)(𝒙 + 𝟏)(𝒙 − 𝟏)(𝒙 − 𝟐) − 𝒙(𝒙 + 𝟐)(𝒙 + 𝟏)(𝒙 − 𝟐)
𝟒 2𝟔
1 𝟏
+ . 𝒙(𝒙 + 𝟐)(𝒙 + 𝟏)(𝒙 − 𝟏)
5 𝟐𝟒
1 1
𝑃4 (𝑥) = 𝑥(𝑥 + 1)(𝑥 − 1)(𝑥 − 2) − 𝑥(𝑥 + 2)(𝑥 − 1)(𝑥 − 2)
120 12
1 1
+ (𝑥 + 2)(𝑥 + 1)(𝑥 − 1)(𝑥 − 2) − 𝑥(𝑥 + 2)(𝑥 + 1)(𝑥 − 2)
4 12
1
+ 𝑥(𝑥 + 2)(𝑥 + 1)(𝑥 − 1)
120
1 1 1
𝑃4 (𝑥) = 𝑥(𝑥 2 − 1)(𝑥 − 2) − 𝑥(𝑥 2 − 4)(𝑥 − 1) + (𝑥 2 − 1)(𝑥 2 − 4)
120 12 4
1 1
− 𝑥(𝑥 2 − 4)(𝑥 + 1) + 𝑥(𝑥 + 2)(𝑥 2 − 1)
12 120
1 1 1
𝑃4 (𝑥) = 𝑥(𝑥 3 − 2𝑥 2 − 𝑥 + 2) − 𝑥(𝑥 3 − 𝑥 2 − 4𝑥 + 4) + (𝑥 4 − 5𝑥 2 + 4)
120 12 4
1 1
− 𝑥(𝑥 3 + 𝑥 2 − 4𝑥 − 4) + 𝑥(𝑥 3 + 2𝑥 2 − 𝑥 − 2)
12 120
1 1 1
𝑃4 (𝑥) = (𝑥 4 − 2𝑥 3 − 𝑥 2 + 2𝑥) − (𝑥 4 − 𝑥 3 − 4𝑥 2 + 4𝑥) + (𝑥 4 − 5𝑥 2 + 4)
120 12 4
1 4 1
− (𝑥 + 𝑥 3 − 4𝑥 2 − 4𝑥) + (𝑥 4 + 2𝑥 3 − 𝑥 2 − 2𝑥)
12 120
1
𝑃4 (𝑥) = [(1 − 10 + 30 − 10 + 1)𝑥 4 + (−2 + 10 − 10 + 2)𝑥 3
120
+ (−1 + 40 − 150 + 40 − 1)𝑥 2 + 120]
1
𝑃4 (𝑥) = [12𝑥 4 − 72𝑥 2 + 120]
120
𝟏 𝟒 𝟑 𝟐
𝑷𝟒 (𝒙) = 𝒙 − 𝒙 +𝟏
𝟏𝟎 𝟓
f(x) avec P(x)
1,2
1
0,8
0,6
0,4
0,2
0
-2 -1,5 -1 -0,5 0 0,5 1 1,5 2
f(x) P(x)
2. Discuter l’erreur d’interpolation.
L’erreur théorique sur une interpolation est donnée en point d’abscisse 𝑥 par :
1
𝐸(𝑥) = 𝑓(𝑥) − 𝑃𝑛 (𝑥) = 𝛾𝑛+1 (𝑥) − 𝑓 (𝑛+1) (𝜉𝑥 ) 𝑜ù 𝜉 ∈ 𝐼 = (min 𝑥𝑖 , max 𝑥𝑖 )
(𝑛 + 1)!
Elle vérifie
𝑛
1
|𝐸(𝑥)| ≤ |𝛾𝑛+1 (𝑥)| 𝑀 ⏟ |𝑓 (𝑛+1) (𝑡)|
𝑜ù 𝛾𝑛+1 (𝑥) = ∏ 𝑥 − 𝑥𝑘 𝑒𝑡 𝑀𝑛+1 = 𝑚𝑎𝑥
(𝑛 + 1)! 𝑛+1 𝑡∈𝐼 𝑘=0
1
Comme ici on a 5 points d’appui, cette erreur est majorée par : |𝐸(𝑥)| ≤ |𝛾5 (𝑥)| (5)! 𝑀5
On a clairement :
4
𝛾5 (𝑥) = ∏ 𝑥 − 𝑥𝑘 = (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )(𝑥 − 𝑥3 )(𝑥 − 𝑥4 )
𝑘=0
= (𝑥 + 2)(𝑥 + 1)(𝑥 − 0)(𝑥 − 1)(𝑥 − 2)
𝜸𝟓 (𝒙) = 𝒙(𝒙 + 𝟐)(𝒙 + 𝟏)(𝒙 − 𝟏)(𝒙 − 𝟐) = 𝒙(𝒙𝟐 − 𝟒)(𝒙𝟐 − 𝟏)
⏟ |𝑓 (5) (𝑡)|
Il reste à calculer 𝑀5 = 𝑚𝑎𝑥
𝑡∈𝐼
On a :
1
𝑓(𝑥) =
1 + 𝑥2
Dérivée première de 𝑓(𝑥) :
−𝟐𝒙
𝒇′(𝒙) =
(𝟏 + 𝒙𝟐 )𝟐
Dérivée seconde de 𝑓(𝑥) :
−2(1 + 𝑥 2 )2 + 2𝑥. 2(1 + 𝑥 2 ). 2𝑥 −2(1 + 𝑥 2 ) + 8𝑥 2
𝑓 ′′(𝑥) = =
(1 + 𝑥 2 )4 (1 + 𝑥 2 )3
′′ (𝒙)
−𝟐 + 𝟔𝒙𝟐
𝒇 =
(𝟏 + 𝒙𝟐 )𝟑
Dérivée troisième de 𝑓(𝑥) :
12𝑥(1 + 𝑥 2 )3 − (−2 + 6𝑥 2 )x3x2𝑥(1 + 𝑥 2 )2
𝑓 (3) (𝑥) =
(1 + 𝑥 2 )6
(3) 12𝑥(1 + 𝑥 2 ) − 6𝑥(−2 + 6𝑥 2 )
𝑓 (𝑥) =
(1 + 𝑥 2 )4
(3) 12𝑥 + 12𝑥 3 + 12𝑥 − 36𝑥 2 (𝒙 − 𝒙𝟑 )
𝑓 (𝑥) = = 𝟐𝟒
(1 + 𝑥 2 )4 (𝟏 + 𝒙𝟐 )𝟒
Dérivée quatrième de 𝑓(𝑥) :
(1 − 3𝑥 2 )(1 + 𝑥 2 )4 − (𝑥 − 𝑥 3 )x4x2𝑥(1 + 𝑥 2 )3
𝑓 (4) (𝑥) = 24
(1 + 𝑥 2 )8
(4)
(1 − 3𝑥 2 )(1 + 𝑥 2 ) − 8𝑥(𝑥 − 𝑥 3 ) 1 + 𝑥 2 − 3𝑥 2 − 3𝑥 4 − 8𝑥 2 + 8𝑥 4
𝑓 (𝑥) = 24 = 24
(1 + 𝑥 2 )5 (1 + 𝑥 2 )5
(4) 𝟏 − 𝟏𝟎𝒙𝟐 + 𝟓𝒙𝟒
𝑓 (𝑥) = 𝟐𝟒
(𝟏 + 𝒙𝟐 )𝟓
Dérivée cinquième de 𝑓(𝑥) :
(5)
(−20𝑥 + 20𝑥 3 )(1 + 𝑥 2 )5 − (1 − 10𝑥 2 + 5𝑥 4 )x5x2𝑥(1 + 𝑥 2 )4
𝑓 (𝑥) = 24
(1 + 𝑥 2 )10
(−20𝑥 + 20𝑥 3 )(1 + 𝑥 2 ) − 10𝑥(1 − 10𝑥 2 + 5𝑥 4 )
𝑓 (5) (𝑥) = 24
(1 + 𝑥 2 )6
(−20𝑥 − 20𝑥 3 + 20𝑥 3 + 20𝑥 5 ) + (−10𝑥 + 100𝑥 3 − 50𝑥 5 )
𝑓 (5) (𝑥) = 24
(1 + 𝑥 2 )6
−30𝑥 + 100𝑥 3 − 30𝑥 5 𝟐𝟒𝟎𝒙(𝟏 − 𝟏𝟎𝒙𝟐 + 𝟑𝒙𝟒 )
𝑓 (5) (𝑥) = 24 = −
(1 + 𝑥 2 )6 (𝟏 + 𝒙𝟐 )𝟔
L’étude de 𝑓 (5) (𝑥) montre que 𝑀5 = 100.
Dérivée cinquième de f(x)
100
50
0
-2 -1 0 1 2
-50
-100
Finalement
1 100
|𝐸(𝑥)| ≤ |𝛾5 (𝑥)| 𝑀5 = |𝑥(𝑥2 − 4)(𝑥2 − 1)|
5! 5!
𝟓
|𝑬(𝒙)| ≤ |𝒙(𝒙𝟐 − 𝟒)(𝒙𝟐 − 𝟏)|
𝟔
-----------------------------------------------------------------------------------------------------------------
Calcul de
𝑑𝑥
𝐼=∫
1 + 𝑥2
On fait un changement de variable et on pose
sin 𝑢
𝑥 = tan 𝑢 =
cos 𝑢
sin 𝑢 cos 𝑢 cos 𝑢 − sin 𝑢 (− sin 𝑢) cos 2 𝑢 + sin2 𝑢
𝑑𝑥 = 𝑑(tan 𝑢) = 𝑑 ( )= 𝑑𝑢 = 𝑑𝑢
cos 𝑢 cos2 𝑢 cos 2 𝑢
1
𝑑𝑥 = 𝑑𝑢
cos 2 𝑢
Alors
1 1
𝐼=∫ 2 cos 2 𝑢 𝑑𝑢
sin 𝑢
1 + (cos 𝑢 )
1 1 cos 2 𝑢 1
𝐼=∫ 2 2 2
𝑑𝑢 = ∫ 𝑑𝑢 = ∫ 𝑑𝑢 = 𝑢 + 𝐶
cos 𝑢 + sin 𝑢 cos 𝑢 cos 𝑢 + sin 𝑢 cos 2 𝑢
⏟ 2 2
cos 2 𝑢 1
𝑥 = tan 𝑢 ⇒ 𝑢 = arctan 𝑥
On a alors :
𝑑𝑥
𝐼=∫ = arctan 𝑥 + 𝐶
1 + 𝑥2
La valeur exacte de l’intégrale est :
2
𝑑𝑥
∫ 2
= 𝑎𝑟𝑐𝑡𝑛 𝑥|2−2 = 1,10714872 − (−1,10714872) = 𝟐, 𝟐𝟏𝟒𝟐𝟗𝟕𝟒𝟒
−2 1 + 𝑥
1 4 3 2
𝑃4 (𝑥) = 𝑥 − 𝑥 +1
10 5
∫ 𝑓(𝑥)𝑑𝑥 ≈ ∫ 𝑃𝑛 (𝑥)𝑑𝑥
2 2
1 4 3 2 1 𝑥5 3 𝑥3 25 23
∫ ( 𝑥 − 𝑥 + 1) 𝑑𝑥 = − + 𝑥⌋ = 2 ( − + 2) = 2x1,04 = 𝟐. 𝟎𝟖
−2 10 5 10 5 5 3 −2
50 5
Méthode des rectangles
𝑥0 𝑥1 𝑥2 𝑥3 𝑥4
𝑥 −2 −1 0 1 2
𝑓(𝑥) 1 1 1 1 1
5 2 2 5
Méthodes des rectangles à gauche
𝑛−1
𝑏−𝑎
𝑅𝑛𝑔 = ∑ 𝑓(𝑥𝑖 )
𝑛
𝑖=0
On a :
𝑏 − 𝑎 2 − (−2)
𝑎 = −2, 𝑏=2 𝑒𝑡 𝑛 = 4 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = =1
𝑛 4
3
𝑅4𝑔 = 1. ∑ 𝑓(𝑥𝑖 )
𝑖=0
𝑅4𝑔 = 𝑓(𝑥0 ) + 𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 )
1 1 1
𝑅4𝑔 = + + 1 + = 𝟐, 𝟐
5 2 2
Méthodes des rectangles à droite
𝑛
𝑏−𝑎
𝑅𝑛𝑑 = ∑ 𝑓(𝑥𝑖 )
𝑛
𝑖=1
𝑅4𝑑 = 1. ∑ 𝑓(𝑥𝑖 )
𝑖=1
𝑅4𝑑 = 𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 ) + 𝑓(𝑥4 )
1 1 1
𝑅4𝑑 = + 1 + + = 𝟐, 𝟐
2 2 5
Méthode des trapèzes
𝑥0 𝑥1 𝑥2 𝑥3 𝑥4
𝑥 −2 −1 0 1 2
𝑓(𝑥) 1 1 1 1 1
5 2 2 5
On a :
𝑏 − 𝑎 2 − (−2)
𝑎 = −2, 𝑏=2 𝑒𝑡 𝑛 = 4 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = =1
𝑛 4
𝑛−1
𝑏 − 𝑎 𝑓(𝑎) + 𝑓(𝑏)
𝑇𝑛 = ( + ∑ 𝑓(𝑥𝑖 ))
𝑛 2
𝑖=1
𝑛−1
ℎ
𝑇𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ))
2
𝑖=1
ℎ
𝑇4 = [𝑓(𝑎) + 𝑓(𝑏) + 2(𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 ))]
2
1 1 1 1 1
𝑇4 = [ + + 2 ( + 1 + )]
2 5 5 2 2
1
𝑇4 = + 2 = 𝟐, 𝟎𝟐
5
Méthode de Simpson
Intervalle 2
Intervalle 1
𝑥0 = 𝑎 𝜉0 𝑥1 𝜉2 𝑥2 = 𝑏
𝑥 −2 −1 0 1 2
𝑓(𝑥) 1 1 1 1 1
5 2 2 5
𝑓(𝑎) 𝑓(𝜉0 ) 𝑓(𝑥1 ) 𝑓(𝜉1 ) 𝑓(𝑏)
On a :
𝑏 − 𝑎 𝑏 − 𝑎 2 − (−2)
𝑎 = −2, 𝑏=2 𝑒𝑡 𝑛 = 2 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = = =2
𝑛 2 2
𝑛−1 𝑛−1
ℎ
𝑆𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6
𝑖=1 𝑖=0
2
𝑆2 = (𝑓(𝑎) + 𝑓(𝑏) + 2𝑓(𝑥1 ) + 4(𝑓(𝜉0 ) + 𝑓(𝜉1 )))
6
1 1 1 1 1 1 2 2
𝑆2 = . ( + + 2x1 + 4 ( + )) = . ( + 6) = 2 + = 𝟐, 𝟏𝟑𝟑𝟑𝟑
3 5 5 2 2 3 5 15
Intervalle 2 Intervalle 4
Intervalle 1 Intervalle 3
𝑥0 𝜉0 𝑥1 𝜉1 𝑥2 𝜉2 𝑥3 𝜉3 𝑥4
𝑥 −2 3 −1 1 0 1 1 3 2
− −
2 2 2 2
𝑓(𝑥) 1 4 1 4 1 4 1 4 1
5 13 2 5 5 2 13 5
𝑓(𝑥0 ) 𝑓(𝜉0 ) 𝑓(𝑥1 ) 𝑓(𝜉1 ) 𝑓(𝑥2 ) 𝑓(𝜉2 ) 𝑓(𝑥3 ) 𝑓(𝜉3 ) 𝑓(𝑥4 )
On a :
𝑏 − 𝑎 𝑏 − 𝑎 2 − (−2)
𝑎 = −2, 𝑏=2 𝑒𝑡 𝑛 = 4 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = = =1
𝑛 4 4
𝑛−1 𝑛−1
ℎ
𝑆𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6
𝑖=1 𝑖=0
1
𝑆4 = (𝑓(𝑎) + 𝑓(𝑏) + 2(𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 )) + 4(𝑓(𝜉0 ) + 𝑓(𝜉1 ) + 𝑓(𝜉2 ) + 𝑓(𝜉3 )))
6
1 1 1 1 1 4 4 4 4
𝑆4 = ( + + 2 ( + 1 + ) + 4 ( + + + ))
6 5 5 2 2 13 5 5 13
1 2 4 4 4 4
𝑆4 = ( + 4 + 4 ( + + + )) = 2.210256
6 5 13 5 5 13
𝑺𝟒 = 𝟐. 𝟐𝟏𝟎𝟐𝟓𝟔
Résultats
Valeur exacte
2
𝑑𝑥
∫ 2
= 𝑎𝑟𝑐𝑡𝑛 𝑥|2−2 = 𝟐, 𝟐𝟏𝟒𝟐𝟗𝟕𝟒𝟒
−2 1 + 𝑥
Méthode des rectangles
𝑅4𝑔 = 𝑅4𝑔 = 𝟐, 𝟐
Méthode des trapèzes
𝑇4 = 𝟐, 𝟎𝟐
Méthode de Simpson
𝑆2 = 𝟐, 𝟏𝟑𝟑𝟑𝟑
𝑺𝟒 = 𝟐. 𝟐𝟏𝟎𝟐𝟓𝟔