INTEGRATION NUMERIQUE

Le but de ce chapitre est de donner des méthodes permettant de calculer des valeurs approchées
d’intégrales en effet, le nombre de fonctions dont on sait calculer une primitive est en fait très
faible. Par ailleurs la connaissance d’une primitive 𝐹 ne suffit pas lorsque l’on ne sait pas
2 𝑑𝑡
calculer les valeurs de 𝐹. Par exemple, on sait que ∫1 𝑡 = ln 2 ; mais connait-on
𝑏
explicitement ln 2 ? Il est donc intéressant concernant ∫𝑎 𝑓(𝑡) 𝑑𝑡 :
 Sur le plan pratique, de pouvoir obtenir une approximation lorsque les primitives de 𝑓
ne sont pas calculables.
 Sur le plan théorique, de connaître des méthodes permettant d’obtenir des encadrements
d’amplitude aussi petite que souhaitée.
Lorsque la fonction 𝑓 est de classe 𝐶 𝑛 sur l’intervalle réel 𝐼 = [𝑎, 𝑏], on note :

𝑀𝑖 = max{|𝑓 (𝑖) (𝑥)|; 𝑥 ∈ [𝑎, 𝑏]} 𝑝𝑜𝑢𝑟 𝑖 = 0,1, … , 𝑛

𝑏−𝑎
On subdivise l’intervalle [𝑎, 𝑏] en 𝑛 intervalles 𝑛 ∈ 𝑁 ∗ de même longueur ℎ = et on note,
𝑛
pour tout 𝑖 ∈ {0,1, … , 𝑛} 𝑥𝑖 = 𝑎 + 𝑖ℎ.

1. Méthode des rectangles

1.1 principe

𝑓(𝑥) On remplace 𝑓 par la fonction escalier qui

prend, sur chaque segment de la
subdivision, la même valeur à l’extrémité
gauche de ce segment que 𝑓. Cela revient
donc à interpoler la fonction 𝑓 sur le
segment [𝑥𝑖 , 𝑥𝑖+1 ] par le polynôme de
Lagrange de degré 0 qui vaut 𝑓(𝑥𝑖 ).
𝑎 𝑥1 𝑥2 𝑏 𝑥

𝑎 𝑥1 𝑥2 𝑏

Proposition 1

La valeur approchée de l’intégrale de f sur I par la méthode des rectangles à gauche est
alors donnée par :

1
INTEGRATION [Link]
 𝑛−1
𝑏−𝑎
𝑅𝑛 = ∑ 𝑓(𝑥𝑖 )
𝑛
𝑖=0

La valeur approchée de l’intégrale de f sur I par la méthode des rectangles à droite est
alors donnée par :
𝑛
𝑏−𝑎
𝑅𝑛 = ∑ 𝑓(𝑥𝑖 )
𝑛
𝑖=1

Démonstration
L’aire du rectangle de base [𝑥𝑖 , 𝑥𝑖+1 ] est 𝑓(𝑥𝑖 ). (𝑥𝑖+1 − 𝑥𝑖 ), donc :
𝑛−1 𝑛−1
𝑏−𝑎
𝑅𝑛 = ∑ 𝑓(𝑥𝑖 )(𝑥𝑖+1 − 𝑥𝑖 ) = ∑ 𝑓(𝑥𝑖 )
𝑛
𝑖=0 𝑖=0

𝑏−𝑎
Car 𝑥𝑖+1 − 𝑥𝑖 = ℎ = 𝑛

1.2 Evaluation de l’erreur

Proposition 2

Si 𝑓 est de classe 𝐶 1 sur [𝑎, 𝑏], alors on a, pour tous entier naturel 𝑛 non nul :

𝑏
𝑀1
|𝑅𝑛 − ∫ 𝑓(𝑡)𝑑𝑡| ≤ (𝑏 − 𝑎)2
𝑎 2𝑛

𝑏
On en déduit que 𝑅𝑛 converge vers ∫𝑎 𝑓(𝑡)𝑑𝑡

2. Méthode des trapèzes

2.1 principe

On remplace la courbe représentative de 𝑓,

sur chaque segment de la subdivision, par
le segment qui joint (𝑥𝑖 , 𝑓(𝑥𝑖 ))
à (𝑥𝑖+1 , 𝑓(𝑥𝑖+1 )). Cela revient donc à
interpoler la fonction 𝑓 sur le
segment [𝑥𝑖 , 𝑥𝑖+1 ] par le polynôme de
Lagrange de degré 1 aux points 𝑥𝑖 et 𝑥𝑖+1
𝑎 𝑥1 𝑥2 𝑏

Proposition 3

2
INTEGRATION [Link]
 La valeur approchée de l’intégrale de 𝑓 sur 𝐼 par la méthode des trapèzes est alors
donnée par :
𝑛−1
𝑏 − 𝑎 𝑓(𝑎) + 𝑓(𝑏)
𝑇𝑛 = ( + ∑ 𝑓(𝑥𝑖 ))
𝑛 2
𝑖=1
Démonstration
L’aire du trapèze de hauteur [𝑥𝑖 , 𝑥𝑖+1 ] est :

(𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 )) (𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 ))

(𝑥𝑖+1 − 𝑥𝑖 ) =ℎ
2 2
𝑛−1
(𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 ))
𝑇𝑛 = ∑ ℎ
2
𝑖=0
(𝑓(𝑥0 ) + 𝑓(𝑥1 )) (𝑓(𝑥1 ) + 𝑓(𝑥2 ))
=ℎ +ℎ +⋯
2 2
(𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 ))
+ℎ
2
On en déduit que :
𝑛−1 𝑛−1
(𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 )) 𝑓(𝑎) + 𝑓(𝑏)
𝑇𝑛 = ∑ ℎ = ℎ( + ∑ 𝑓(𝑥𝑖 ))
2 2
𝑖=0 𝑖=1

2.2 Evaluation de l’erreur

Proposition 4
Si 𝑓 est de classe 𝐶 2 sur [𝑎, 𝑏], alors on a, pour tout entier naturel 𝑛 non nul :

𝑏
𝑀2
|𝑇𝑛 − ∫ 𝑓(𝑡)𝑑𝑡| ≤ (𝑏 − 𝑎)2
𝑎 12𝑛2

𝑏
On en déduit que 𝑇𝑛 converge vers ∫𝑎 𝑓(𝑡)𝑑𝑡

3. Méthode Simpson
3.1 principe

On remplace 𝑓, sur chaque

segment [𝑥𝑖 , 𝑥𝑖+1 ] de la subdivision, par la
fonction polynômiale de degré inférieur ou
égal à 2 qui prend les mêmes valeurs que 𝑓
aux extrémités et au milieu 𝜉𝑖 de ce segment.
Cette méthode consiste à remplacer 𝑓 sur
le segment [𝑥𝑖 , 𝑥𝑖+1 ] par son polynôme
𝑥0 = 𝑎 𝜉0 𝑥1 𝜉1 𝑥2 𝜉2 𝑥3 = 𝑏 d’interpolation 𝑃𝑖 de Lagrange de degré 2
𝑎 ayant les mêmes valeurs que 𝑓 aux bornes
𝑏
de l’intervalle et en son milieu.
3
INTEGRATION [Link]
 Proposition 5

La valeur approchée de l’intégrale de 𝑓 sur 𝐼 par la méthode de Simpson est alors donnée
par :
𝑛−1
𝑏−𝑎
𝑆𝑛 = ∑(𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 ) + 4𝑓(𝜉𝑖 ))
6𝑛
𝑖=0
Ou
𝑛−1 𝑛−1
𝑏−𝑎
𝑆𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6𝑛
𝑖=1 𝑖=0
On a
𝑛−1 𝑥𝑖+1
𝑆𝑛 = ∑ ∫ 𝑃𝑖 (𝑡) 𝑑𝑡
𝑖=0 𝑥𝑖
On est donc ramené au calcul de l’intégrale de Riemann d’un polynôme de degré
inférieur ou égal à 2.

Lemme 6
Soit 𝑃 ∈ 𝑅2 [𝑋] 𝑒𝑡 𝑐 < 𝑑, alors :
𝑑
𝑑−𝑐 𝑐+𝑑
∫ 𝑃(𝑡)𝑑𝑡 = [𝑃(𝑐) + 𝑃(𝑑) + 4𝑃 ( )]
𝑐 6 2

Démonstration
On fait un changement de variable 𝑡 = (𝑑 − 𝑐)𝑢 + 𝑐 pour se ramener à une intégrale
entre 0 et 1
𝑢 = 0 ⇒ 𝑡 = 𝑐; 𝑢 = 1 ⇒ 𝑡 = 𝑑 𝑒𝑡 𝑑𝑡 = (𝑑 − 𝑐)𝑑𝑢
𝑑 1
∫ 𝑃(𝑡)𝑑𝑡 = (𝑑 − 𝑐) ∫ 𝑃((𝑑 − 𝑐)𝑢 + 𝑐)𝑑𝑢
𝑐 0

Les polynômes de Lagrange aux points 0,1 et ½ sont :

𝑥0 = 𝑎 𝑥1 = 𝜉 𝑥2 = 𝑏
0 1 1
2
𝒏 𝒏
𝒙 − 𝒙𝒋
𝑷𝒏 (𝒙) = ∑ 𝒇(𝒙𝒊 ) 𝑳𝒊 (𝒙) 𝒂𝒗𝒆𝒄 𝑳𝒊 (𝒙) = ∏
𝒙𝒊 − 𝒙𝒋
𝒊=𝟎 𝒋=𝟎
𝒊≠𝒋
(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )
𝑖=0 𝐿0 (𝑥) =
(𝑥0 − 𝑥1 )(𝑥0 − 𝑥2 )
1
(𝑥 − 2) (𝑥 − 1)
𝐿0 = = (2𝑥 − 1)(𝑥 − 1) = 2𝑥 2 − 3𝑥 + 1
1
(0 − 2) (0 − 1)

4
INTEGRATION [Link]
 𝐿0 (𝑥) = 2𝑥 2 − 3𝑥 + 1

(𝑥 − 𝑥0 )(𝑥 − 𝑥2 )
𝑖=1 𝐿1 =
(𝑥1 − 𝑥0 )(𝑥1 − 𝑥2 )
(𝑥 − 0)(𝑥 − 1)
𝐿1 = = −4𝑥(𝑥 − 1) = −4𝑥 2 + 4𝑥
1 1
(2 − 0) (2 − 1)

𝐿1 (𝑥) = 4𝑥 − 4𝑥 2

(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )
𝑖=2 𝐿2 =
(𝑥2 − 𝑥0 )(𝑥2 − 𝑥1 )

1
(𝑥 − 0) (𝑥 − )
𝐿2 = 2 = 𝑥(2𝑥 − 1) = 2𝑥 2 − 𝑥
1
(1 − 0) (1 − )
2

𝐿2 (𝑥) = 2𝑥 2 − 𝑥
1 1 1
2 3 2 3 1
∫ 𝐿0 (𝑥) = ∫ (2𝑥 2 − 3𝑥 + 1)𝑑𝑥 = [ 𝑥 3 − 𝑥 2 + 𝑥] = ( − + 1) =
0 0 3 2 0 3 2 6

1 1
4 3 4 2 1 4 2 4
∫ 𝐿1 (𝑥) = ∫ (−4𝑥 2 + 4𝑥)𝑑𝑥 = [− 𝑥 + 𝑥 ] = (− + 2) = =
0 0 3 2 0 3 3 6

1 1 1
2 𝑥2 2 1 1
∫ 𝐿2 (𝑥) = ∫ (2𝑥 2 − 𝑥)𝑑𝑥 = [ 𝑥 3 − ] = ( − ) =
0 0 3 2 0 3 2 6

Leurs intégrales entre 0 et 1 sont respectivement : 1/6, 4/6 et 1/6.

Soit 𝑄(𝑢) = 𝑃((𝑑 − 𝑐)𝑢 + 𝑐), 𝑄 est un polynôme de degré inférieur ou égal à 2, on a
donc :
1
𝑄 = 𝑄(0)𝐿0 + 𝑄 ( ) 𝐿1 + 𝑄(1)𝐿2
2
Et
1
1 1
∫ 𝑄(𝑢)𝑑𝑢 = [𝑄(0) + 𝑄(1) + 4𝑄 ( )]
0 6 2
Comme
1 𝑐+𝑑
𝑄(0) = 𝑃(𝑐), 𝑄(1) = 𝑃(𝑑) 𝑒𝑡 𝑄 ( ) = 𝑃 ( )
2 2
On obtient le Lemme.
On applique le Lemme à 𝑆𝑛 et on obtient :
𝑛−1 𝑥𝑖+1 𝑛−1
𝑥𝑖 + 𝑥𝑖+1
𝑆𝑛 = ∑ ∫ 𝑃𝑖 (𝑡) = ∑ (𝑃𝑖 (𝑥𝑖 ) + 𝑃𝑖 (𝑥𝑖+1 ) + 4𝑃𝑖 (𝜉𝑖 ))
2
𝑖=0 𝑥𝑖 𝑖=0

5
INTEGRATION [Link]
 𝑛−1
𝑥𝑖 + 𝑥𝑖+1
𝑆𝑛 = ∑ (𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 ) + 4𝑓(𝜉𝑖 )) par définition de 𝑃𝑖
2
𝑖=0
𝑛−1 𝑛−1 𝑛−1
ℎ
𝑆𝑛 = (∑ 𝑓(𝑥𝑖 ) + ∑ 𝑓(𝑥𝑖+1 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6
𝑖=0 𝑖=0 𝑖=0

𝑛−1 𝑛−1
𝑏−𝑎
𝑆𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6𝑛
𝑖=1 𝑖=0

3.2 Evaluation de l’erreur

Proposition 4
Si 𝑓 est de classe 𝐶 3 sur [𝑎, 𝑏], alors on a, pour tout entier naturel 𝑛 non nul, on a :

𝑏
𝑀3
|𝑆𝑛 − ∫ 𝑓(𝑡)𝑑𝑡| ≤ (𝑏 − 𝑎)4
𝑎 192𝑛3

𝑏
On en déduit que 𝑇𝑛 converge vers ∫𝑎 𝑓(𝑡)𝑑𝑡

6
INTEGRATION [Link]
 APPLICATION
Exercice 1
𝜋
Déterminer par la méthode des trapèzes puis par celle de Simpson ∫02 𝑓(𝑥)𝑑𝑥 sur la base du
tableau suivant :

𝑥 0 𝜋 𝜋 3𝜋 𝜋
8 4 8 2
𝑓(𝑥) 0 0,382683 0,707107 0,923880 1

Ces points d’appui sont ceux donnant 𝑠𝑖𝑛 𝑥, comparer alors les résultats obtenus avec la
valeur exacte.

Corrigé exercice 1

𝑥0 𝑥1 𝑥2 𝑥3 𝑥4
𝑥 0 𝜋 𝜋 3𝜋 𝜋
8 4 8 2
𝑓(𝑥) 0 0,382683 0,707107 0,923880 1
On a :
𝜋
𝜋 𝑏−𝑎 2−0 𝜋
𝑎 = 0, 𝑏= 𝑒𝑡 𝑛 = 4 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = =
2 𝑛 4 8

 Méthode des trapèzes

𝑛−1
𝑏 − 𝑎 𝑓(𝑎) + 𝑓(𝑏)
𝑇𝑛 = ( + ∑ 𝑓(𝑥𝑖 ))
𝑛 2
𝑖=1

𝑛−1
ℎ
𝑇𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ))
2
𝑖=1

ℎ
𝑇4 = [𝑓(𝑎) + 𝑓(𝑏) + 2(𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 ))]
2

𝑥0 = 𝑎 𝑥1 𝑥2 𝑥3 𝑥4 = 𝑏
𝑥 0 𝜋 𝜋 3𝜋 𝜋
8 4 8 2
𝑓(𝑥) 0 0,382683 0,707107 0,923880 1

𝜋
𝑇4 = [0 + 1 + 2(0,382683 + 0,707107 + 0,923880)]
16
𝜋 𝜋
𝑇4 = [0 + 1 + 2x2,01367] = x5,02734 = 0,987116
16 8

𝑻𝟒 = 𝟎, 𝟗𝟖𝟕𝟏𝟏𝟔

7
INTEGRATION [Link]
 𝜋
𝜋
2 𝜋
∫ 𝑠𝑖𝑛 𝑥 𝑑𝑥 = −𝑐𝑜𝑠 𝑥|02 = − (𝑐𝑜𝑠 − 𝑐𝑜𝑠 0) = −(0 − 1) = 1
0 2

 Méthode de Simpson

Intervalle 2
Intervalle 1
𝑥0 = 𝑎 𝜉0 𝑥1 𝜉2 𝑥2 = 𝑏
𝑥 0 𝜋 𝜋 3𝜋 𝜋
8 4 8 2
𝑓(𝑥) 0 0,382683 0,707107 0,923880 1
𝑓(𝑎) 𝑓(𝜉0 ) 𝑓(𝑥1 ) 𝑓(𝜉1 ) 𝑓(𝑏)
On a :
𝜋
𝜋 𝑏−𝑎 𝑏−𝑎 2−0 𝜋
𝑎 = 0, 𝑏= 𝑒𝑡 𝑛 = 2 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = = =
2 𝑛 2 2 4
𝑛−1 𝑛−1
ℎ
𝑆𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6
𝑖=1 𝑖=0

1 𝜋
𝑆2 = . (𝑓(𝑎) + 𝑓(𝑏) + 2𝑓(𝑥1 ) + 4(𝑓(𝜉0 ) + 𝑓(𝜉1 )))
6 4

Intervalle 2
Intervalle 1
𝑥0 = 𝑎 𝜉0 𝑥1 𝜉1 𝑥2 = 𝑏
𝑥 0 𝜋 𝜋 3𝜋 𝜋
8 4 8 2
𝑓(𝑥) 0 0,382683 0,707107 0,923880 1
𝑓(𝑎) 𝑓(𝜉0 ) 𝑓(𝑥1 ) 𝑓(𝜉1 ) 𝑓(𝑏)

𝜋
𝑆2 = (0 + 1 + 2x0,707107 + 4(0,382683 + 0,923880))
24
𝜋
𝑆2 = (0 + 1 + 2x0,707107 + 4x1,306563)
24
𝜋
𝑆2 = x7,640466
24

𝑺𝟐 = 𝟏, 𝟎𝟎𝟎𝟏𝟑𝟓

|𝐸 − 𝑺𝟐 | |1 − 𝟏, 𝟎𝟎𝟎𝟏𝟑𝟓|
𝐸𝑟𝑟𝑒𝑢𝑟 𝑒𝑛% = x100 = x100 = 0,135%
𝐸 1

8
INTEGRATION [Link]
 Exercice 2

On lance une fusée verticalement du sol et l’on mesure pendant les premières 80 secondes
l’accélération 𝛾 :

𝑡 (𝑒𝑛 𝑠) 0 10 20 30 40 50 60 70 80
𝛾 (𝑒𝑛 𝑚/𝑠 2 ) 30 31,63 33,44 35,47 37,75 40,33 43,29 46,70 50,67

Calculer la vitesse 𝑉 de la fusée à l’instant 𝑡 = 80 𝑠, par la méthode des trapèzes puis par
Simpson.

Corrigé exercice 2

On sait que l’accélération 𝛾 est la dérivée de la vitesse 𝑉, donc :

𝑡 80
𝑉(𝑡) = 𝑉(0) + ∫ 𝛾(𝓈) 𝑑𝓈 ⇒ 𝑉(80) = 0 + ∫ 𝛾(𝓈) 𝑑𝓈
0 ⏟0
𝐼

1. Calculons I par la méthode des trapèzes.

𝑥0 = 𝑎 𝑥1 𝑥2 𝑥3 𝑥4 𝑥5 𝑥6 𝑥7 𝑥8 = 𝑏
𝑡 (𝑒𝑛 𝑠) 0 10 20 30 40 50 60 70 80
𝛾 (𝑒𝑛 𝑚/𝑠 2 ) 30 31,63 33,44 35,47 37,75 40,33 43,29 46,70 50,67

D’après le tableau des valeurs on remarque que :

𝑏−𝑎
𝑎 = 0, 𝑏 = 80 𝑛=8 ℎ= = 10
𝑛
𝑛−1
𝑏 − 𝑎 𝑓(𝑎) + 𝑓(𝑏)
𝑇𝑛 = ( + ∑ 𝑓(𝑥𝑖 ))
𝑛 2
𝑖=1

7
ℎ
𝑇8 = (𝛾(𝑎) + 𝛾(𝑏) + 2 ∑ 𝛾(𝑥𝑖 ))
2
𝑖=1
𝑥0 = 𝑎 𝑥1 𝑥2 𝑥3 𝑥4 𝑥5 𝑥6 𝑥7 𝑥8 = 𝑏
𝑡 (𝑒𝑛 𝑠) 0 10 20 30 40 50 60 70 80
𝛾 (𝑒𝑛 𝑚/𝑠 2 ) 30 31,63 33,44 35,47 37,75 40,33 43,29 46,70 50,67

10
𝑇8 = (30 + 50,67
2
+ 2(31,63 + 33,44 + 35,47 + 37,75 + 40,33 + 43,29 + 46,70))

9
INTEGRATION [Link]
 10
𝑇8 = (30 + 50,67 + 2𝑥268,61) = 5x617,89 = 3089,45 𝑚/𝑠
2

𝑽(𝟖𝟎) ≈ 𝑻𝟖 = 𝟑𝟎𝟖𝟗, 𝟒𝟓 𝒎/𝒔

2. Calculons I par la méthode de Simpson.

On a :
𝑏 − 𝑎 80 − 0
𝑎 = 0, 𝑏 = 80 𝑒𝑡 𝑛 = 4 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = = 20
𝑛 4
𝑛−1 𝑛−1
ℎ
𝑆𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6
𝑖=1 𝑖=0

Intervalle 2 Intervalle 4
Intervalle 1 Intervalle 3
𝑡 (𝑒𝑛 𝑠) 0 10 20 30 40 50 60 70 80
2
𝛾 (𝑒𝑛 𝑚/𝑠 ) 30 31,63 33,44 35,47 37,75 40,33 43,29 46,70 50,67
𝛾(𝑎) 𝛾(𝜉0 ) 𝛾(𝑥1 ) 𝛾(𝜉1 ) 𝛾(𝑥2 ) 𝛾(𝜉2 ) 𝛾(𝑥3 ) 𝛾(𝜉3 ) 𝛾(𝑏)
3 3
ℎ
𝑆𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ) + 4 ∑ 𝑓(𝜉𝑖 ))
6
𝑖=1 𝑖=0

1
𝑆𝑛 = x20 (𝑓(𝑎) + 𝑓(𝑏) + 2(𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 ))
6
+ 4(𝑓(𝜉0 ) + 𝑓(𝜉1 ) + 𝑓(𝜉2 ) + 𝑓(𝜉3 )))

1
𝑆𝑛 = x20 (𝛾(0) + 𝛾(80) + 2(𝛾(20) + 𝛾(40) + 𝛾(60))
6
+ 4(𝛾(10) + 𝛾(30) + 𝛾(50) + 𝛾(70)))

1
𝑆𝑛 = x20(30 + 50,67 + 2(33,44 + 37,75 + 43,29)
6
+ 4(31,63 + 35,47 + 40,33 + 46,70))

1
𝑆𝑛 = x20(30 + 50,67 + 2x114,48 + 4x154,13)
6
1
𝑆𝑛 = x20x926,15
6

𝑽(𝟖𝟎) ≈ 𝑺𝒏 = 𝟑𝟎𝟖𝟕, 𝟏𝟔𝟔𝟔𝟕 𝒎/𝒔

10
INTEGRATION [Link]
Exercice 3

Calculer à l’aide de la méthode des trapèzes l’intégrale :

𝜋
𝐼 = ∫ sin 𝑥 2 𝑑𝑥
0
avec le nombre de points d’appui 𝑛 = 5 puis 𝑛 = 10.

Corrigé exercice 3

Soit :
𝜋
𝐼 = ∫ sin 𝑥 2 𝑑𝑥
0

 Nombre de points d’appui 𝒏 = 𝟓

On a :
𝑏−𝑎 𝜋−0 𝜋
𝑎 = 0, 𝑏=𝜋 𝑒𝑡 𝑛 = 5 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = =
𝑛 5 5

𝑥 0 𝜋 2𝜋 3𝜋 4𝜋 𝜋
5 5 5 5
sin(𝑥)2 0 0,38460898 0,99996522 -0,39995242 0,03335532 -0,43030122

𝑛−1
ℎ
𝑇𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ))
2
𝑖=1

4
ℎ
𝑇5 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ))
2
𝑖=1
ℎ
𝑇5 = [𝑓(𝑎) + 𝑓(𝑏) + 2(𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 ) + 𝑓(𝑥4 ))]
2

ℎ 𝜋 2𝜋 3𝜋 4𝜋
𝑇5 = [𝑓(0) + 𝑓(𝜋) + 2 (𝑓 ( ) + 𝑓 ( ) + 𝑓 ( ) + 𝑓 ( ))]
2 5 5 5 5

1 𝜋
𝑇5 = x [sin(0)2 + sin(𝜋)2
2 5
𝜋 2 2𝜋 2 3𝜋 2 4𝜋 2
+ 2 (sin ( ) + sin ( ) + sin ( ) + sin ( ) )]
5 5 5 5
1 𝜋
𝑇5 = x [0 + (−0,43030122)
2 5
+ 2(0,38460898 + 0,99996522 + (−0,39995242) + 0,03335532)]
𝜋
𝑇5 = x1,60565298 = 0,50443076
10

𝑻𝟓 = 𝟎, 𝟓𝟎𝟒𝟒𝟑𝟎𝟕𝟔
11
INTEGRATION [Link]
  Nombre de points d’appui 𝒏 = 𝟏𝟎

On a :
𝑏−𝑎 𝜋−0 𝜋
𝑎 = 0, 𝑏=𝜋 𝑒𝑡 𝑛 = 10 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑙𝑒𝑠 ⇒ ℎ= = =
𝑛 10 10
𝑥 0 𝜋 2𝜋 3𝜋 4𝜋 5𝜋 6𝜋 7𝜋 8𝜋 9𝜋 𝜋
10 10 10 10 10 10 10 10 10
sin(𝑥)2 0 0,0985 0,3846 0,7759 0,9999 0,6242 -0,3999 -0,9923 0,0333 0,9901 -0,4303

𝑛−1
ℎ
𝑇𝑛 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ))
2
𝑖=1

9
ℎ
𝑇10 = (𝑓(𝑎) + 𝑓(𝑏) + 2 ∑ 𝑓(𝑥𝑖 ))
2
𝑖=1
ℎ
𝑇10 = [𝑓(𝑎) + 𝑓(𝑏)
2
+ 2(𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 ) + 𝑓(𝑥4 ) + 𝑓(𝑥5 ) + 𝑓(𝑥6 ) + 𝑓(𝑥7 )
+ 𝑓(𝑥8 ) + 𝑓(𝑥9 ))]

1 𝜋
𝑇10 = x [0 + (−0,4303)
2 10
+ 2(0,0985 + 0,3846 + 0,7759 + 0,9999 + 0,6242 + (−0,3999)
+ (−0,9923) + 0,0333 + 0,9901 + (−0,4303))]

𝜋
𝑇10 = [(−0,4303) + 2x5,0291] = 0,7224
20

𝑻𝟏𝟎 = 𝟎, 𝟕𝟐𝟐𝟒

Alors que la valeur ‘exacte’ est approximativement 0,772651. Avec ce pas plus petit
l’approximation numérique est meilleure

12
INTEGRATION [Link]
Exercice 4

Corrigé exercice 4

13
INTEGRATION [Link]