Module : Méthodes numériques pour l’ingénieur

Chapitre 2 : Intégration Numérique

Partie 1 : Introduction

Itégration Numérique
Préparé par :
Unité Pédagogique de Mathématiques

A.U. 2022/2023
 Objectif du chapitre

Le but de ce chapitre est d’abordrer le calcul numérique de l’intégrale

d’une fonction f (x) sur un domaine fini délimité par des bornes finies a et
b tel que a < b.
Z b
I = f (x)dx.
a

Définition
Graphiquement : Soit f une fonc-
tion continue et positive sur un inter-
valle [a, b].
On appelle intégrale de a à b de la
Rb
fonction f , et on note a f (x)dx l’aire
en unité d’aire de la partie du plan
délimitée par l’axe des abscisses, les
droites d’équations x = a et x = b et
la courbe représentative de f .

2/7
 Motivations
Dans certains cas limités, une telle intǵrale peut être calculée
analytiquement.
Primitive usuelle:
Z π
2 π
cos xdx = [sinx]02 = 1.
0

Changement de variable:
Z π
2
2x cos(x 2 )dx =?
0

On pose: u = x 2 donc du = 2xdx et on obtient:

Z π Z π2
2 4
2
2x cos(x )dx = cos(u)du.
0 0

π2
= sin( ).
4
3/7
 Intégration par partie:
Z π
2
x cos xdx =?
0
On pose
u(x) = x → u 0 (x) = 1,
v 0 (x) = cos x → v (x) = sin(x).
Z π Z π2
2 π
x cos xdx = [x sin x]02 − sin xdx,
0 0

π π
= + [cos x]02 ,
2
π
= − 1.
2

4/7
 Le plus souvent un des cas suivants se présente:
Le calcul analytique de l’intégrale est long et compliqué.
Le résultat de l’intégrale est une fonction compliquée .
L’intégrale n’aR pas d’expression analytique, par expemle la fonction
x 2
Erf (x) = √2π 0 e −t dt.

5/7
 Le plus souvent un des cas suivants se présente:
Le calcul analytique de l’intégrale est long et compliqué.
Le résultat de l’intégrale est une fonction compliquée .
L’intégrale n’aR pas d’expression analytique, par expemle la fonction
x 2
Erf (x) = √2π 0 e −t dt.
Pour surmonter ces problèmes, on préfèrera calculer numériquement la
valeur de l’intégrale I. Cependant il existe une vaste familles de méthodes
dont le but est d’estimer la valeur exacte de l’intégrale.
⇒ les méthodes d’intégrations numériques.

5/7
 On utilisera, dans ce chapitre, quelques méthodes de quadrature de type
Newton-Cotes qui seront définies dans ce diagramme.

6/7
 Diagramme

7/7
 Formule de quadrature
Rappel

Interpolation de Lagrange
Soient f une fonction continue sur [a, b] et (xi , yi = f (xi ))0≤i≤n n + 1 points d’interpolation
tel que a ≤ x0 ≤ x1 ≤ · · · ≤ xn ≤ b, alors f est interpolée par un polynôme d’interpolation
de Lagrange P ∈ Rn [X ] tel que

f (x) = P(x) + Ep (f ),

avec
Pn
P(x) = i=0 f (xi )Li (x) est le polynôme d’interpolation de Lagrange.
Ep (f ) est l’erreur d’interpolation de Lagrange.

1/1
 Formule de quadrature
Rappel

Interpolation de Lagrange
Soient f une fonction continue sur [a, b] et (xi , yi = f (xi ))0≤i≤n n + 1 points d’interpolation
tel que a ≤ x0 ≤ x1 ≤ · · · ≤ xn ≤ b, alors f est interpolée par un polynôme d’interpolation
de Lagrange P ∈ Rn [X ] tel que

f (x) = P(x) + Ep (f ),

avec
Pn
P(x) = i=0 f (xi )Li (x) est le polynôme d’interpolation de Lagrange.
Ep (f ) est l’erreur d’interpolation de Lagrange.

L’intégrale I (f ) peut s’écrire comme suit:

Z b Z b Z b
I (f ) = f (x)dx = P(x)dx + Ep (f )dx,
a a a
Z b n
X
= f (xi )Li (x)dx + Eq (f ),
a i=0
Rb
où Eq (f ) = a Ep (f )dx est l’erreur d’intégration.
Rb
Si on pose Wi = a Li (x)dx, on trouve
Z b Xn
I (f ) = f (x)dx = Wi f (xi ) + Eq (f ).
a i=0

1/1
 Formule de quadrature
Rappel

Interpolation de Lagrange
Soient f une fonction continue sur [a, b] et (xi , yi = f (xi ))0≤i≤n n + 1 points d’interpolation
tel que a ≤ x0 ≤ x1 ≤ · · · ≤ xn ≤ b, alors f est interpolée par un polynôme d’interpolation
de Lagrange P ∈ Rn [X ] tel que

f (x) = P(x) + Ep (f ),

avec
Pn
P(x) = i=0 f (xi )Li (x) est le polynôme d’interpolation de Lagrange.
Ep (f ) est l’erreur d’interpolation de Lagrange.

L’intégrale I (f ) peut s’écrire comme suit:

Z b Z b Z b
I (f ) = f (x)dx = P(x)dx + Ep (f )dx,
a a a
Z b n
X
= f (xi )Li (x)dx + Eq (f ),
a i=0
Rb
où Eq (f ) = a Ep (f )dx est l’erreur d’intégration.
Rb
Si on pose Wi = a Li (x)dx, on trouve
Z b Xn
I (f ) = f (x)dx = Wi f (xi ) + Eq (f ).
a i=0
Notre objectif est d’approcher I (f ) par une somme finie.
1/1
 Formule de quadrature

On appelle formule de quadrature ou formule d’intégration numérique toute

formule permettant de calculer une approximation de I (f ).

Définition
On dit qu’une formule de quadrature de type interpolation toute formule
s’écrie sous la forme suivante :
Z b n
X
I (f ) = f (x)dx ' Iq (f ) = Wi f (xi ),
a i=0

où
Rb
Wi = a
Li (x)dx : poids de la formule de quadrature.
xi : nœuds ou points d’intégration.

2/1
 Formule de quadrature

On appelle formule de quadrature ou formule d’intégration numérique toute

formule permettant de calculer une approximation de I (f ).

Définition
On dit qu’une formule de quadrature de type interpolation toute formule
s’écrie sous la forme suivante :
Z b n
X
I (f ) = f (x)dx ' Iq (f ) = Wi f (xi ),
a i=0

où
Rb
Wi = a
Li (x)dx : poids de la formule de quadrature.
xi : nœuds ou points d’intégration.
Exemple : Soit f une fonction continue.
Les formules suivantes sont des formules de quadrature:
R1
1
−1
f (x)dx ' Iq (f ) = 13 f (−1) + 34 f (0) + 13 f (1).
R1
2
0
f (x)dx ' Iq (f ) = 0.5f (0) + 0.4238f (1).

2/1
 Formule de quadrature
Définition

Définitions
L’erreur de quadrature est donnée par

Eq (f ) = |I (f ) − Iq (f )| .

La formule de quadrature est exacte lorsque l’erreur de quadrature Eq (f ) vaut zéro.

3/1
 Formule de quadrature
Définition

Définitions
L’erreur de quadrature est donnée par

Eq (f ) = |I (f ) − Iq (f )| .

La formule de quadrature est exacte lorsque l’erreur de quadrature Eq (f ) vaut zéro.

Exercice : Parmi les formules de quadrature précédentes, indiquez laquelle est exacte, avec
1
f (x) = .
1.2626 + x

3/1
 Formule de quadrature
Définition

Définitions
L’erreur de quadrature est donnée par

Eq (f ) = |I (f ) − Iq (f )| .

La formule de quadrature est exacte lorsque l’erreur de quadrature Eq (f ) vaut zéro.

Exercice : Parmi les formules de quadrature précédentes, indiquez laquelle est exacte, avec
1
f (x) = .
1.2626 + x
Solution :

3/1
 Formule de quadrature
Définition

Définitions
L’erreur de quadrature est donnée par

Eq (f ) = |I (f ) − Iq (f )| .

La formule de quadrature est exacte lorsque l’erreur de quadrature Eq (f ) vaut zéro.

Exercice : Parmi les formules de quadrature précédentes, indiquez laquelle est exacte, avec
1
f (x) = .
1.2626 + x
Solution Z:
1 h i1 2.2626
1 On a f (x)dx = log(1.2626 + x) = log( ) = 2.1536.
−1 −1 0.2626
1
Or Iq (f ) = 3f (−1) + 43 f (0) + 31 f (1) = 13 0.2626
1
+ 3 1.2626 + 13 2.2626
4 1 1
= 2.4727.
Alors
Eq (f ) = I (f ) − Iq (f ) = 2.1536 − 2.4727 = 0.3191

D’où, cette formule n’est pas exacte.

3/1
 Formule de quadrature
Définition

Définitions
L’erreur de quadrature est donnée par

Eq (f ) = |I (f ) − Iq (f )| .

La formule de quadrature est exacte lorsque l’erreur de quadrature Eq (f ) vaut zéro.

Exercice : Parmi les formules de quadrature précédentes, indiquez laquelle est exacte, avec
1
f (x) = .
1.2626 + x
Solution Z:
1 h i1 2.2626
1 On a f (x)dx = log(1.2626 + x) = log( ) = 2.1536.
−1 −1 0.2626
1
Or Iq (f ) = 3f (−1) + 43 f (0) + 31 f (1) = 13 0.2626
1
+ 3 1.2626 + 13 2.2626
4 1 1
= 2.4727.
Alors
Eq (f ) = I (f ) − Iq (f ) = 2.1536 − 2.4727 = 0.3191

D’où, cette formule n’est pas exacte.

Z 1 h i1 2.2626
2 On a f (x)dx = log(1.2626 + x) = log( ) = 0.5833.
0 0 1.2626
1 1
Or Iq (f ) = 0.5f (0) + 0.4238f (1) = 0.5 1.2626 + 0.4238 2.2626 = 0.5833
Alors
Eq (f ) = I (f ) − Iq (f ) = 0.5833 − 0.5833 = 0

D’où, cette formule est exacte.

3/1
 Formule de quadrature
Degré de précision

4/1
 Formule de quadrature
Degré de précision

Définition
Une formule de quadrature est dite de degré de précision (degré
d’exactitude) n, s’il elle est exacte pour tout polynôme Pk (x) = x k ,
k = 0, 1, · · · , n et non exacte pour Pn+1 (x) = x n+1 .

4/1
 Formule de quadrature
Degré de précision

Définition
Une formule de quadrature est dite de degré de précision (degré
d’exactitude) n, s’il elle est exacte pour tout polynôme Pk (x) = x k ,
k = 0, 1, · · · , n et non exacte pour Pn+1 (x) = x n+1 .

Proposition
Une formule de quadrature a un degré de précision au moins égale n si et
seulement si c’est une formule de quadrature interpolatoire à n + 1 points.

4/1
 Formule de quadrature
Degré de précision

Définition
Une formule de quadrature est dite de degré de précision (degré
d’exactitude) n, s’il elle est exacte pour tout polynôme Pk (x) = x k ,
k = 0, 1, · · · , n et non exacte pour Pn+1 (x) = x n+1 .

Proposition
Une formule de quadrature a un degré de précision au moins égale n si et
seulement si c’est une formule de quadrature interpolatoire à n + 1 points.

Exemple : Soit f une fonction continue.

La formule de quadrature
Z 1
1 4 1
f (x)dx ' Iq (f ) = f (−1) + f (0) + f (1)
−1 3 3 3

est de type interpolation à 3 points, donc elle est exacte pour tout polynôme de
degré ≤ 2.

4/1
 Formule de quadrature
Exercice

Soient f une fonction continue sur [−1, 1] et I (f ) l’intégrale suivante

Z 1
I (f ) = f (x)dx
−1

On considère la formule de quadrature Iq (f ) suivante approchant la valeur

de I (f ):
Iq (f ) = αf (−1) + βf (1),
avec α et β deux réels.
1. Trouver α et β.
2. Montrer que cette formule est de degré de précision 1.

5/1
 Formule de quadrature
solution
1. D’après la proposition, on remarque que cette formule de quadrature est de type interpolation à 2
points, donc elle est exacte pour tout polynôme de degré ≤ 1.
Trouvons α et β : Soit Pk (x) = x k avec k = 0, 1.
Pour k = 0, on a f (x) = P0 (x) = 1, d’une part
Z 1 Z 1
f (x)dx = 1 dx = 2.
−1 −1

D’autre part, Iq (f ) = αf (−1) + βf (1) = α × 1 + β × 1 = α + β.

Puisque la formule de quadrature est exacte pour un polynôme de degré 0, alors on a
I (f ) = Iq (f ), par suite
α+β =2 (1)
Pour k = 1, on a f (x) = P1 (x) = x, d’une part
Z 1 Z 1 h x 2 i1
f (x)dx = x dx = = 0.
−1 −1 2 −1

D’autre part, Iq (f ) = αf (−1) + βf (1) = α × (−1) + β × 1 = −α + β.

Puisque la formule de quadrature est exacte pour un polynôme de degré 1, alors on a
I (f ) = Iq (f ), par suite
−α + β = 0. (2)
Les deux équations (1) et (2) donnent le système linéaire suivant :
(
2=α+β
⇔α=β=1
0 = −α + β
6/1
 Formule de quadrature
solution

Donc la formule de quadrature est :

Z 1
f (x) dx ' f (−1) + f (1).
−1

2. Montrons que le degré de précision de cette formule est 1 i.e. on montre que
la formule n’est pas exacte pour un polynôme de degré 2.
Posons P2 (x) = x 2 :
on a Iq (P2 ) = P2 (−1) + P2 (1) = 2.
Z 1 h x 3 i1 2
P2 (x) dx = = 6= Iq (P2 ).
−1 3 −1 3

Alors le degré de précision est 1.

7/1
 Méthodes simples du rectangle
Rectangle à gauche

La méthode simple du rectangle à gauche consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f (a)

Soit f une fonction continue sur [a, b], alors la Z b

formule simple du rectangle à gauche est: s
I (f ) = f (x)dx w IRg (f )
s
a
IRg (f ) = (b − a)f (a)

1/1
 Méthodes simples du rectangle
Rectangle à gauche

La méthode simple du rectangle à gauche consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f (a)

Soit f une fonction continue sur [a, b], alors la Z b

formule simple du rectangle à gauche est: s
I (f ) = f (x)dx w IRg (f )
s
a
IRg (f ) = (b − a)f (a)

Proposition
Si f ∈ C 1 ([a, b]), alors l’erreur du rectangle à gauche est
majorée par:

s s (b − a)2 0
ERg (f ) = I (f ) − IRg (f ) ≤ max f (x)
2 x∈[a,b]

1/1
 Méthodes simples du rectangle
Rectangle à gauche

La méthode simple du rectangle à gauche consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f (a)

Soit f une fonction continue sur [a, b], alors la Z b

formule simple du rectangle à gauche est: s
I (f ) = f (x)dx w IRg (f )
s
a
IRg (f ) = (b − a)f (a)

Proposition
Si f ∈ C 1 ([a, b]), alors l’erreur du rectangle à gauche est
majorée par:

s s (b − a)2 0
ERg (f ) = I (f ) − IRg (f ) ≤ max f (x)
2 x∈[a,b]

Degré de précision: Le degré de précision de la méthode du rectangle à gauche est 0.

1/1
 Méthodes simples du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Utiliser la méthode simple du rectangle à gauche pour calculer I (f ).
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :

2/1
 Méthodes simples du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Utiliser la méthode simple du rectangle à gauche pour calculer I (f ).
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1.
Z 1
1
I (f ) = dx
0 1+x
= [ln(1 + x)]10
= ln(2) − ln(1)
= ln(2) ' 0.693

2/1
 Méthodes simples du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Utiliser la méthode simple du rectangle à gauche pour calculer I (f ).
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1.
Z 1
1
I (f ) = dx
0 1+x
= [ln(1 + x)]10
= ln(2) − ln(1)
= ln(2) ' 0.693

s
2. IRg (f ) = (1 − 0)f (0) = 1

2/1
 Méthodes simples du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Utiliser la méthode simple du rectangle à gauche pour calculer I (f ).
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1.
Z 1
1
I (f ) = dx
0 1+x
= [ln(1 + x)]10
= ln(2) − ln(1)
= ln(2) ' 0.693

s
2. IRg (f ) = (1 − 0)f (0) = 1
s
3. ERg (f ) = |I (f ) − IRg (f )| = |ln(2) − 1| '
0.307
2/1
 Méthodes simples du rectangle
Rectangle à droite

La méthode simple du rectangle à droite consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f (b).

Soit f une fonction continue sur [a, b], alors la

formule simple du rectangle à droite est: Z b
s
s
I (f ) = f (x)dx w IRd (f )
IRd (f ) = (b − a)f (b). a

3/1
 Méthodes simples du rectangle
Rectangle à droite

La méthode simple du rectangle à droite consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f (b).

Soit f une fonction continue sur [a, b], alors la

formule simple du rectangle à droite est: Z b
s
s
I (f ) = f (x)dx w IRd (f )
IRd (f ) = (b − a)f (b). a

Proposition
Si f ∈ C 1 ([a, b]), alors l’erreur du rectangle à droite est
majorée par

s s (b − a)2 0
ERd (f ) = I (f ) − IRd (f ) ≤ max f (x)
2 x∈[a,b]

3/1
 Méthodes simples du rectangle
Rectangle à droite

La méthode simple du rectangle à droite consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f (b).

Soit f une fonction continue sur [a, b], alors la

formule simple du rectangle à droite est: Z b
s
s
I (f ) = f (x)dx w IRd (f )
IRd (f ) = (b − a)f (b). a

Proposition
Si f ∈ C 1 ([a, b]), alors l’erreur du rectangle à droite est
majorée par

s s (b − a)2 0
ERd (f ) = I (f ) − IRd (f ) ≤ max f (x)
2 x∈[a,b]

Degré de précision: Le degré de précision de la méthode du rectangle à droite est 0.

3/1
 Méthodes simples du rectangle à droite
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Utiliser la méthode simple du rectangle à droite pour calculer I (f ).

2 Quelle est l’erreur d’intégration pour cette méthode.

Solution :

4/1
 Méthodes simples du rectangle à droite
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Utiliser la méthode simple du rectangle à droite pour calculer I (f ).

2 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
s 1
1. IRd (f ) = (1 − 0)f (1) = 2

4/1
 Méthodes simples du rectangle à droite
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Utiliser la méthode simple du rectangle à droite pour calculer I (f ).

2 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
s 1
1. IRd (f ) = (1 − 0)f (1) = 2
2.
s s
ERd (f ) = |I (f ) − IRd (f )|
1
= ln(2) −
2
' 0.193

4/1
 Méthodes simples du rectangle
Rectangle au milieu

La méthode simple du rectangle au milieu consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f ( a+b
2
)

Soit f une fonction continue sur [a, b], alors la formule

simple du rectangle au milieu est: Z b
s
a+b
I (f ) = f (x)dx w IRm (f )
s a
IRm (f ) = (b − a)f ( )
2

5/1
 Méthodes simples du rectangle
Rectangle au milieu

La méthode simple du rectangle au milieu consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f ( a+b
2
)

Soit f une fonction continue sur [a, b], alors la formule

simple du rectangle au milieu est: Z b
s
a+b
I (f ) = f (x)dx w IRm (f )
s a
IRm (f ) = (b − a)f ( )
2

Proposition
Si f ∈ C 2 ([a, b]), alors l’erreur du rectangle au milieu est
majorée par

s s (b − a)3 00
ERm (f ) = |I (f ) − IRm (f )| ≤ max |f (x)|
24 x∈[a,b]

5/1
 Méthodes simples du rectangle
Rectangle au milieu

La méthode simple du rectangle au milieu consiste à approximer l’intégrale I (f )

par l’aire d’un rectangle dont les deux côtés sont (b − a) et f ( a+b
2
)

Soit f une fonction continue sur [a, b], alors la formule

simple du rectangle au milieu est: Z b
s
a+b
I (f ) = f (x)dx w IRm (f )
s a
IRm (f ) = (b − a)f ( )
2

Proposition
Si f ∈ C 2 ([a, b]), alors l’erreur du rectangle au milieu est
majorée par

s s (b − a)3 00
ERm (f ) = |I (f ) − IRm (f )| ≤ max |f (x)|
24 x∈[a,b]

Degré de précision: Le degré de précision de la méthode du rectangle au milieu est 1.

5/1
 Méthodes simples du rectangle au milieu

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.

6/1
 Méthodes simples du rectangle au milieu

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.
La méthode du rectangle au milieu est :

s 0+1
IRm (f ) = (1 − 0)f ( )
2
1
= f( )
2
2
=
3
' 0.666

6/1
 Méthodes simples du rectangle au milieu

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.
La méthode du rectangle au milieu est :

s 0+1
IRm (f ) = (1 − 0)f ( )
2
1
= f( )
2
2
=
3
' 0.666
s
L’erreur est : ERm (f ) =
s 2
|I (f ) − IRm (f )| = | ln(2) − | ' 0.026
3

6/1
 Méthode simple du trapèze

La méthode simple de trapèse consiste à approximer l’intégrale I (f ) par l’aire

d’un trapèze dont ces deux bases sont f (a) et f (b) et son hauteur est (b − a)

Soit f une fonction continue sur [a, b], alors la for-

Z b
mule simple de trapèze est: I (f ) = f (x)dx w ITs (f )
a
s b−a
IT (f ) = (f (a) + f (b))
2

7/1
 Méthode simple du trapèze

La méthode simple de trapèse consiste à approximer l’intégrale I (f ) par l’aire

d’un trapèze dont ces deux bases sont f (a) et f (b) et son hauteur est (b − a)

Soit f une fonction continue sur [a, b], alors la for-

Z b
mule simple de trapèze est: I (f ) = f (x)dx w ITs (f )
a
s b−a
IT (f ) = (f (a) + f (b))
2

Proposition
Si f ∈ C 2 ([a, b]), alors l’erreur de trapèze est majorée par

s s (b − a)3 00
ET (f ) = |I (f ) − IT (f )| ≤ max |f (x)|
12 x∈[a,b]

7/1
 Méthode simple du trapèze

La méthode simple de trapèse consiste à approximer l’intégrale I (f ) par l’aire

d’un trapèze dont ces deux bases sont f (a) et f (b) et son hauteur est (b − a)

Soit f une fonction continue sur [a, b], alors la for-

Z b
mule simple de trapèze est: I (f ) = f (x)dx w ITs (f )
a
s b−a
IT (f ) = (f (a) + f (b))
2

Proposition
Si f ∈ C 2 ([a, b]), alors l’erreur de trapèze est majorée par

s s (b − a)3 00
ET (f ) = |I (f ) − IT (f )| ≤ max |f (x)|
12 x∈[a,b]

Degré de précision: Le degré de précision de la méthode de trapèze est 1.

7/1
 Méthode simple du trapèze

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.

8/1
 Méthode simple du trapèze

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.
La méthode du trapèze est :

1−0
ITs (f ) = (f (0) + f (1))
2
1 1
= (1 + )
2 2
3
=
4
= 0.75

8/1
 Méthode simple du trapèze

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.
La méthode du trapèze est :

1−0
ITs (f ) = (f (0) + f (1))
2
1 1
= (1 + )
2 2
3
=
4
= 0.75

L’erreur est: ETs (f ) = |I (f ) − ITs (f )| =

| ln(2) − 0.75| ' 0.056

8/1
 Méthode simple de Simpson

La méthode simple de Simpson est obtenue en interpolant f par un polynôme

P ∈ R2 [X ] aux points d’abscisse a, a+b
2
et b, puis en intégrant P sur [a, b].

Soit f une fonction continue sur [a, b], alors la formule

Z b
simple de Simpson est: I (f ) w P(x)dx = ISs (f )
a
b−a
 
s a+b
IS (f ) = f (a) + 4f ( ) + f (b)
6 2

9/1
 Méthode simple de Simpson

La méthode simple de Simpson est obtenue en interpolant f par un polynôme

P ∈ R2 [X ] aux points d’abscisse a, a+b
2
et b, puis en intégrant P sur [a, b].

Soit f une fonction continue sur [a, b], alors la formule

Z b
simple de Simpson est: I (f ) w P(x)dx = ISs (f )
a
b−a
 
s a+b
IS (f ) = f (a) + 4f ( ) + f (b)
6 2

Proposition
Si f ∈ C 4 ([a, b]), alors l’erreur de Simpson est majorée par

s s (b − a)5 (4)
ES (f ) = |I (f ) − IS (f )| ≤ max |f (x)|
2880 x∈[a,b]

9/1
 Méthode simple de Simpson

La méthode simple de Simpson est obtenue en interpolant f par un polynôme

P ∈ R2 [X ] aux points d’abscisse a, a+b
2
et b, puis en intégrant P sur [a, b].

Soit f une fonction continue sur [a, b], alors la formule

Z b
simple de Simpson est: I (f ) w P(x)dx = ISs (f )
a
b−a
 
s a+b
IS (f ) = f (a) + 4f ( ) + f (b)
6 2

Proposition
Si f ∈ C 4 ([a, b]), alors l’erreur de Simpson est majorée par

s s (b − a)5 (4)
ES (f ) = |I (f ) − IS (f )| ≤ max |f (x)|
2880 x∈[a,b]

Degré de précision: Le degré de précision de la méthode de Simpson est 3.

9/1
 Méthode simple de Simpson

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.

10 / 1
 Méthode simple de Simpson

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.
La méthode de Simpson est :
 
1−0 0+1
ISs (f ) = f (0) + 4f ( ) + f (1)
6 2
 
1 2 1
= 1 + 4. +
6 3 2
25
= ' 0.694
36

10 / 1
 Méthode simple de Simpson

Exemple :
Z 1
1
Soit l’intégrale I (f ) = dx
0 1 + x
La valeur exacte de I (f ) vaut ln(2) = 0.693.
La méthode de Simpson est :
 
1−0 0+1
ISs (f ) = f (0) + 4f ( ) + f (1)
6 2
 
1 2 1
= 1 + 4. +
6 3 2
25
= ' 0.694
36

L’erreur est:
ESs (f ) = |I (f ) − ISs (f )| = | ln(2) − 0.694| '
0.001

10 / 1
 Récapitulation

Degré du Nombre de Degré

Nom de la Méthode
polynôme points d’exactitude
Rectangle à gauche 0 1 0
Rectangle à droite 0 1 0
Rectangle au milieu 0 1 1
Trapèze 1 2 1
Simpson 2 3 3

11 / 1
 Exercice

2
Soit la fonction f (t) = e −t . Donner une valeur approchée de
Z 1
I (f ) = f (t)dt par la méthode du
−1
rectangle du milieu.
trapèze.
Simpson.

12 / 1
 Solution
Rectangle du milieu.

s −1 + 1
IRm (f ) = (1 − (−1))f ( ) = 2f (0)
2
=2

13 / 1
 Solution
Rectangle du milieu.

s −1 + 1
IRm (f ) = (1 − (−1))f ( ) = 2f (0)
2
=2

Trapèze.

s 1 − (−1)
IT (f ) = (f (−1) + f (1)) = f (−1) + f (1)
2
−1
= 2e ' 0.735

13 / 1
 Solution
Rectangle du milieu.

s −1 + 1
IRm (f ) = (1 − (−1))f ( ) = 2f (0)
2
=2

Trapèze.

s 1 − (−1)
IT (f ) = (f (−1) + f (1)) = f (−1) + f (1)
2
−1
= 2e ' 0.735

Simpson.

s 1 − (−1) −1 + 1
IS (f ) = (f (−1) + 4f ( ) + f (1))
6 2
1
= (f (−1) + 4f (0) + f (1))
3
1 −1 −1 1 −1
= (e + 4 + e ) = (4 + 2e ) ' 1.578
3 3

13 / 1
 Solution
Rectangle du milieu.

s −1 + 1
IRm (f ) = (1 − (−1))f ( ) = 2f (0)
2
=2

Trapèze.

s 1 − (−1)
IT (f ) = (f (−1) + f (1)) = f (−1) + f (1)
2
−1
= 2e ' 0.735

Simpson.

s 1 − (−1) −1 + 1
IS (f ) = (f (−1) + 4f ( ) + f (1))
6 2
1
= (f (−1) + 4f (0) + f (1))
3
1 −1 −1 1 −1
= (e + 4 + e ) = (4 + 2e ) ' 1.578
3 3

13 / 1
 Solution
Rectangle du milieu.

s −1 + 1
IRm (f ) = (1 − (−1))f ( ) = 2f (0)
2
=2

Trapèze.

s 1 − (−1)
IT (f ) = (f (−1) + f (1)) = f (−1) + f (1)
2
−1
= 2e ' 0.735

Simpson.

s 1 − (−1) −1 + 1
IS (f ) = (f (−1) + 4f ( ) + f (1))
6 2
1
= (f (−1) + 4f (0) + f (1))
3
1 −1 −1 1 −1
= (e + 4 + e ) = (4 + 2e ) ' 1.578
3 3

Est-il possible d’améliorer ces résultats en utilisant d’autres méthodes?

13 / 1
 Solution
Rectangle du milieu.

s −1 + 1
IRm (f ) = (1 − (−1))f ( ) = 2f (0)
2
=2

Trapèze.

s 1 − (−1)
IT (f ) = (f (−1) + f (1)) = f (−1) + f (1)
2
−1
= 2e ' 0.735

Simpson.

s 1 − (−1) −1 + 1
IS (f ) = (f (−1) + 4f ( ) + f (1))
6 2
1
= (f (−1) + 4f (0) + f (1))
3
1 −1 −1 1 −1
= (e + 4 + e ) = (4 + 2e ) ' 1.578
3 3

Est-il possible d’améliorer ces résultats en utilisant d’autres méthodes? Méthodes

composites d’intégration

13 / 1