25/04/2026 12:42 TP3_INT_Num
Méthodes Numériques
TD-TP3 : Intégration numérique
Objectif :
–––––––––
Le but de ce TP est l’implémentation et la comparaison des méthodes d'intégration numérique de base :
b
Rectangles, Trapèzes, Simpson, Gauss-Legendre, pour calculer une approximation de ∫a f (x)dx quant on
ne peut pas trouver la primitive de f , ou quand f n'est connue que par des points de mesure.
Soient a et b deux réels (a < b ) et f : [a, b] → R . On considère une subdivision uniforme de l'intervalle
b−a
[a, b] en n sous-intervalles [xi , xi+1 ] pour i = 0, ⋯ , n − 1 , on note h =
n
le pas de la subdivision,
b
,
x i = a + ih i = 0, ⋯ , n . On note I (f ) = ∫
a
f (x)dx .
Méthode des rectangles :
On remplace l'aire sous la courbe par l'aire des rectangles.
Formule version rectange à gauche :
n−1
I (f ) ≈ h ∑ f (x i ), O(h)
i=0
Formule version rectange à droite :
n
I (f ) ≈ h ∑ f (x i ), O(h)
i=1
Formule version point milieu :
n−1
1
2
I (f ) ≈ h ∑ f (a + (i + )h), O(h )
2
i=0
localhost:8888/nbconvert/html/Activites_Enseignement/CoursENSAM/Cours_MNI1_25_26/TPs/TP3_INT_Num.ipynb?download=false 1/5
25/04/2026 12:42 TP3_INT_Num
Méthode des trapèzes
On remplace la courble par des segments.
Formule :
n−1
h
I (f ) ≈ [f (a) + f (b) + 2 ∑ f (x i )] := Th (f )
2
i=1
L'erreur :
Si f est de classe C 2 alors il existe η ∈ [a, b] tel que:
′′
(b − a)f (η)
2
E(f ) = I (f ) − Th (f ) = − h
12
L'erreur est O(h2 ).
Méthode de Simpson
On remplace la courble par des paraboles. n doit être pair.
Formule :
n n
−1 −1
h⎡ ⎤
2 2
I (f ) ≈ ⎢f (a) + f (b) + 4 ∑ f (x 2i+1 ) + ∑ f (x 2i )⎥ := Sh (f )
3 ⎣ ⎦
i=0 i=1
L'erreur :
Si f est de classe C 4 alors il existe η ∈ [a, b] tel que:
(4)
(b − a)f (η)
4
E(f ) = I (f ) − Sh (f ) = − h
180
L'erreur est O(h4 ).
localhost:8888/nbconvert/html/Activites_Enseignement/CoursENSAM/Cours_MNI1_25_26/TPs/TP3_INT_Num.ipynb?download=false 2/5
25/04/2026 12:42 TP3_INT_Num
Exercice 1
On considère l'intégrale :
2
x−1
I = ∫ e dx
1
1) Calculer la valeur exacte de I .
2) Considérons une subdivision uniforme de l'intervalle [1, 2] en 4 sous-intervalles [xi , xi+1 ], où
1
,
x i = 1 + ih i = 0, 1, ⋯ , 4 h = , 4
, donner une valeur approchée de I , par la méthode des trapèzes
puis par celle de Simpson.
3) Déduire une approximation de e par les deux méthodes.
4) Donner une majoration d'erreur d'approximation de I .
5) Combien faut il de sous intervalles pour assurer une erreur inférieure à 10−3 par les deux méthodes?
Exercice 2
Soient les donn\'ees fournies dans le tableau suivant :
1 1 3
x 0 1
4 2 4
f (x) 1 0.941176 0.8 0.64 0.5
1) Déterminer par la méthode des trapèzes puis par celle de Simpson une valeur approchée de
1
I = ∫
0
f (x) dx sur la base du tableau.
1
2) Les points d'appui donnés dans le tableau ci-dessus correspondent à la fonction x ↦ , donner une
1+x2
valeur approchée de π par les deux méthodes.
−4
3) Trouver le nombre minimal n de subdivisions nécessaires de l'intervalle [0, 1], pour évaluer à 10 près,
la valeur de π par la méthode des trapèzes.
Exercice 3 (TP) : Calcul de π
Objectif :
Calculer une valeur approcheé de π par les trois méthodes et comparer à la valeur de référence [Link]
1) Coder les quatres fonctions :
rectangles(f,a,b,n)
trapezes(f,a,b,n)
simpson(f,a,b,n)
localhost:8888/nbconvert/html/Activites_Enseignement/CoursENSAM/Cours_MNI1_25_26/TPs/TP3_INT_Num.ipynb?download=false 3/5
25/04/2026 12:42 TP3_INT_Num
In [ ]: import numpy as np
def f(x):
return 4/(1+x**2)
def rectangles(f,a,b,n):
return I
def trapezes(f,a,b,n):
return
def Simpson(f,a,b,n):
return
2) Tester les méthodes pour n = 4, 8, 16, 20, 80
In [ ]:
3) Tracer sur un même graphe l'erreur obtenue par les trois méthodes en fonction de h en échelle log-log.
Conclure.
In [ ]:
4) Vérifie numériquement l'ordre en utilisant la formule :
e h1
log( )
e h2
p ≈
e h1
log( )
e h1
In [ ]:
−−−−−
5) On peut aussi calculer une valeur approchée de π en intégrant la fonction x ↦ √1 − x
2
sur [0, 1] :
1
−−−−−
2
π = 4∫ √1 − x dx.
0
Reprendre les questions (1)-(4) de l'exercice puis comparer les résultats obtenus.
localhost:8888/nbconvert/html/Activites_Enseignement/CoursENSAM/Cours_MNI1_25_26/TPs/TP3_INT_Num.ipynb?download=false 4/5
25/04/2026 12:42 TP3_INT_Num
Exercice 4 (Facultatif) : Gauss-Legendre
1
Pour approcher l'intégrale ∫−1 f (x) dx, on considère la formule de quadrature J (f ) suivante :
1
J (f ) = A0 f (t1 ) + A1 f ( )
–
√3
1) Déterminer les poids d'intégration A0 et A1 , ainsi que le point d'intégration t1 de sorte que la formule de
quadrature J (f ) soit de précision la plus élevée possible.
2) Donner le degré d'exactitude de la formule.
3) Déduire la formule sur un intervalle [a, b] quelconque.
4) Utiliser la formule pour calculer une valeur approchée de π avec 101 point, comparer avec la méthode de
Simpson composite.
localhost:8888/nbconvert/html/Activites_Enseignement/CoursENSAM/Cours_MNI1_25_26/TPs/TP3_INT_Num.ipynb?download=false 5/5