TD Maple n7 : Méthodes d'intégration numérique

Lycée Louis le Grand

PCSI 1

Lundis 11 et 18 février 2008

1 Calculs de primitives et d'intégrales

1. Caluler les intégrales et les primitives suivantes :

R∞ 2 R∞ x9
R1
e−x dx x4 e−x dx 1 √1
R R
0 0 1+x20
dx 0 cos4 x+sin4 x
dx chx ch2x
dx
2. Soit Z x
1
In (x) = dt
0 (1 + t2 )n
A l'aide d'un intégration par partie (et d'une petite astuce...), déterminer une relation entre In
et In−1 . On pourra charger le package student avec l'instruction with(student).

2 Intégration numérique

Dans cette partie, on étudie plusieurs méthodes numériques permettant de calculer de manière
approchée l'intégrale d'une fonction f [a, b]. Pour cela, on découpe l'intervalle
continue sur un intervalle
b−a
[a, b] en N petits intervalles de longueurs identiques h = . On construit ainsi une subdivision
N
régulière de l'intervalle t0 = a < t1 < .. < tN = b où ti = a + ih.

2.1 Méthode des rectangles

La méthode des rectangles à gauche consiste à interpoler la fonction f sur chaque intervalle [ti , ti+1 ]
par sa valeur en ti , c'est-à-dire qu'on considère que sur le petit intervalle [ti , ti+1 ], la fonction f est
constante égale à f (ti ).

1
 Sur l'intervalle [ti , ti+1 ], l'intégrale est donc approchée par l'air du rectangle soit (ti+1 − ti )f (ti ) =
b−a
hf (ti ) = f (ti ). Donc l'intégrale sur l'intervalle [a, b] est approchée par :
N
b N −1
b−a X
Z
f (t)dt ≈ f (ti )
a N
i=0
Rb
1. Ecrire une procédure rectgauche(f,a,b,N), qui calcule la valeur approchée de
a f (t)dt avec
N intervalles de discrétisation par la méthode des rectangles à gauches.

2. De même, écrire une procédure

Rb rectdroite(f,a,b,N), qui calcule la valeur approchée de

a f (t)dt avec N pas par la méthode des rectangles à droite.

3. Dans quels cas est-on sûr que rectgauche(f,a,b,N) surestime l'intégrale ?

4. Tester vos procédures sur la fonction f (t) = exp(t) sur l'intervalle [0, 1] et pour N = 10.
5. On note Z b
f (t)dt − rectgauche(f,a,b,N)
a
l'erreur commise lorsqu'on utilise la méthode des rectangles à gauche avec N pas. Vérier (avec
1

les mains) que l'erreur commise pour la fonction f (t) = exp(t) est un O = O(h). Pour
N
cela, on pourra tracer le logarithme de l'erreur en fonction de lnh pour N = 10, 20, 30, 40, 50, 60
et calculer la pente de la droite obtenue. On dit alors que la méthode est d'ordre 1.

2.2 Méthode du point milieu

Cette méthode consiste elle aussi à interpoler la fonction f sur l'intervalle[ti , ti+1 ] par une seule
t + ti+1 
i
valeur, la valeur prise au milieu de l'intervalle. Sur [ti , ti+1 ], f (t) ≈ cste = f .
2

2
 Rb
1. Ecrire une procédure pointmilieu(f,a,b,N), qui calcule la valeur approchée de a f (t)dt avec
N pas par la méthode du point milieu.

2. Vérier que cette méthode est aussi d'ordre 2 c'est-à-dire que l'erreur commise en utilisant cette
1
méthode est un O( ) = O(h2 ).
N2
3. Que peut-on dire de l'odre de la méthode si, sur l'intervalle [ti , ti+1 ], on avait pris f (t) ≈
1
f (λti + (1 − λ)ti+1 ) avec λ ∈ [0, 1], λ 6= 2?

2.3 Méthode des trapèzes

Pour obetenir de meilleures approximations le l'intégrale, on peut augmenter la précision de l'inter-

polation de f sur chaque intervalle [ti , ti+1 ], ie augmenter le nombre de points utilisés pour approcher
f sur cet intervalle. La méthode des trapèzes consiste à utiliser les points ti et ti+1 en approchant la
fonction f par la droite ane passant par les points (ti , f (ti )) et (ti+1 , f (ti+1 )).

1. Calculer l'aire du trapèze délimité par la fonction ane sur [ti , ti+1 ]. En déduire une nouvelle
formule pour le calcul approché de l'intégrale sur l'intervalle [a, b].
2. Ecrire la procédure tapeze(f,a,b,N) et tester cette procédure.

3. Vérier que la méthode des trapèzes est d'ordre 2.

2.4 Méthode de Simpson

La méthode de Simpson est une méthode d'ordre 4. Elle repose sur une interpolation à trois points
de la fonction f sur l'intervalle [ti , ti+1 ].

3
 ti + ti+1
1. On note ti+ 1 = le milieu du segment [ti , ti+1 ]. Quel est le polynôme qui s'annule en
2 2
ti+ 1 et en ti+1 et qui vaut 1 en ti ?
2

2. En déduire l'équation du polynôme de degré 2 qui coïncide avec f aux points ti , ti+1 et ti+ 1 .
2

Comme il est très facile d'intégrer des polynômes, on va approcher l'intégrale de f sur l'intervalle
[ti , ti+1 ] par l'air sous la parabole.

3. A l'aide de Maple, exprimer l'aire sous la parabole sur l'intervalle [ti , ti+1 ] en fonction de ti ,
ti+1 , f (ti ), f (ti+1 ) et f (ti+ 1 ).
2

4. En déduire une nouvelle formule pour le calcul approché de l'intégrale de f sur l'intervalle [a, b].
5. Ecrire la procédure simpson(f,a,b,N) puis la tester.

6. Vérier que cette méthode donne une approximation en O(h4 ) de l'intégrale de f. On peut
ainsi construire des méthodes aussi précises que l'on veut en augmentant le degré du polynôme
interpolateur.