UIT. FSK.

Département de Mathématiques Filière SMP : S3

Module Analyse Numérique

Série N◦ 3 :
Z 2
dx
Exercice 1 On considère l'intégrale suivante, I =
1 x
1. Évaluer numériquement cette intégrale par la méthode des trapèzes pour le pas h = 1/3.
2. Calculer la valeur exacte de I .
(a) Pourquoi la valeur numérique obtenue à la question 1. est-elle supérieure à ln 2?
(b) Est-ce vrai quelque soit h ?
(c) Proposer une autre fonction dont la valeur de l'intégrale évaluée par la méthode des
trapèzes est toujours supérieure à la valeur exacte de l'intégrale.
3. Si on souhaite évaluer I avec la méthode de Simpson, quel pas faut-il choisir pour avoir une
erreur inférieure à 10−4 ?
Z x
Exercice 2 Soit F (x) = te−t dt. Combien faut-il de subdivisions de [0, 1] pour évaluer F (1)
0
à 10−5 près en utilisant
1. la méthode du point milieu
2. la méthode de Simpson
Z 1
(Février 2015) Soit I = e−x dx.
2
Exercice 3
0

1. Combien faut-il de subdivisions de [0, 1] pour évaluer I à 10−5 près en utilisant la méthode
de Simpson
1 1 3
2. Donner la valeur approchée de I en utilisant x0 = 0,x1 = ,x2 = , x3 = et x4 = 1
4 2 4
Z 1
dx
Exercice 4 (janvier 2016)On considère l'intégrale suivante, I =
0 1 + x2
1. Rappeler les formules composites du point milieu, des trapèzes et de Simpson, ainsi que
l'erreur commise dans chacune des méthodes.
2. Combien faut-il de subdivisions de [0, 1] pour évaluer I à 10−3 près en utilisant la méthode
de Simpson
1
3. Donner la valeur approchée de I en utilisant la subdivision x0 = 0,x1 = , et x2 = 1
2
4. quelle est alors l'erreur commise

M.H. EL ALJ 1 Décembre 2016

UIT. FSK. Département de Mathématiques Filière SMP : S3

b
b−a
Z
Exercice 5 Soit, I = f (x)dx avec xk = a + k × pour k = 0, . . . , n
a n
Z xk
1. Calculer P0,k (f )(x)dx . oú P0,k (x) est le polynôme de degré 0 qui passe par le point
xk−1
d'abscisse x̄k .
Z xk
2. Calculer P1,k (f )(x)dx . oú P1,k (x) est le polynôme de degré 1 qui passe par les points
xk−1
d'abscisses xk−1 et xk .
3. Combien faut-il de subdivisions de [0, 1] pour évaluer I à 10−3 près en utilisant pour f (t) =
t cos(t) en utilisant

(a) la méthode des Trapèzes

(b) la méthode de Simpson

Module Analyse Numérique

Corrigé de la série N◦ 3 :
Z 2
dx
Solution de l'e xercice 1 On considère l'intégrale suivante, I =
1 x
1. La formule composite du trapèze est donnée par
n−1
h X
ITc (f ) = [f (a) + f (b)] + h f (xk )
2 k=1

. Dans notre cas on a donc pour le pas h = 1/3.,

    
1 1 4 5 21
ITc (f ) = (f (1) + f (2)) + f +f = = 0.7.
6 3 3 3 30

2. La valeur exacte de I est ln 2 ≈ 0.69

(a) On observe eectivement que I T calculer en 1. est supérieur à ln(2) ≈ 0.69. On peut
se convaincre à l'aide d'un dessin que les trapèzes sont au-dessus de la courbe y = 1/x
, l'aire sous les trapèzes sera donc supérieure à l'aire sous la courbe.
(b) Cela reste vrai quelque soit le pas h choisi car la fonction x → 1/x est convexe sur
]0; +∞[ ce qui signie qu'une corde dénie par deux points de la courbe y = 1/x sera
toujours au-dessus de la courbe et par le raisonnement précédant l'aire sous les trapèzes
sera supérieure à l'aire exacte.
(c) Il sut de trouver une fonction dont la dérivée seconde soit positive,
par exemple x → x2 .
h4
3. une condition susante pour avoir |E S | ≤ 10−4 est d'avoir max fx∈[1,2] (x) ≤ 10−4 .
(4)
24 × 120
24
Or f (4) (x) = ce qui implique maxx∈[1,2] |f (x)|(4) (x) = 24.
x5
M.H. EL ALJ 2 Décembre 2016
UIT. FSK. Département de Mathématiques Filière SMP : S3

Donc on devrait avoir

24 × 120 × 10−4
h4 ≤ = 0.012 ⇒ h ≤ 0.3309
24
Donc h = 41 donnera la précision voulue. C'est à dire que
        
1 5 6 7
ln 2 ≈ f (1) + 2 f +f +f + f (2)
4 × 6    4 4  4 
1 9 11 13 15
+ 4 +f +f +f
4 ×6 8 8 8  8 
1 4 4 4 4 4 1 8 8 8 8
= +2 + + + + + + +
24 4 5 6 7 8 6 9 11 13 15
= 0.693154
Z x
Solution de l'e xercice 2 Soit F (x) = te−t dt.On cherche le nombre de subdivisions de [0, 1]
0
pour évaluer F (1) à 10−5 près:
1. la méthode du point milieu
une condition susante pour avoir
h2
|EPc M | ≤ 10−5 est d'avoir (b − a) maxt∈[1,2] |f (2) (t)| ≤ 10−5 .
24
Or (b − a) = 1, dans ce cas et f (2) (t) = (t − 2)e−t
ce qui implique maxt∈[0,1] |f (2) (t)| = 2.
Donc on devrait avoir
h2 ≤ 12 × 10−5 ⇒ h ≤ 0.0109544
Donc h = 1
92
donnera la précision voulue. On aura donc besoin de 92 points de subdivision.
2. la méthode de Simpson
une condition susante pour avoir
h4
|ESc | ≤ 10−5 est d'avoir maxt∈[1,2] |f (4) (t)| ≤ 10−5 . Or f (4) (t) = (t − 4)e−t
2880
ce qui implique maxt∈[0,1] |f (2) (t)| = 4. Donc on devrait avoir
2880 × 10−5
h4 ≤ ⇒ h ≤ 0.29
4
Donc h = 14 donnera la précision voulue. On aura donc besoin de seulement 4 points de
subdivision dans cette méthode.
Solution de l'e xercice 3 (Sessionn de rattrapage Février 2015)
Z 1
soit I = e −x2
dx où f (x) = e−x
2
∀x ∈ [0, 1]
0

1. Nombre de subdivisions pour évaluer I à 10−4 près par la méthode de Simpson

Dérivées successives de f (x) = e−x2 :
f (x) = −2xe−x )
0 2

f ” (t) = (4x2 − 2)e−x

2

f (3) (x) = (−8x3 + 12x)e−x

2

f (4) (x) = (16x4 − 48x2 + 12)e−x

2

M.H. EL ALJ 3 Décembre 2016

UIT. FSK. Département de Mathématiques Filière SMP : S3

et on a (16x4 − 48x2 + 12)e−x ≤ |(16x4 − 48x2 + 12)| e−x

2 2

≤ |(16x2 (x2 − 1)| + |(12(1 − x2 )| + |20x2 | ≤ 48

1
(1 − 0)4
e−x dx − ISC (f ) ≤
Z
2 48
4
max |f (4) (x)| ≤ 4
≤ 10−4
0 120 × 24n x∈[0,1] 120 × 24n
104 √
ce qui entraine que n ≥
4
d'où n ≥ 4 166.66 c'est à dire n ≥ 4.
60
2. Valeur de I en utilisant la subdivision x0 = 1,x1 = 0.25,x2 = 0.5,x3 = 0.75 et x4 = 1.
( 3 4
)
ISC (f ) =
1−0 X X xk−1 + xk
f (0) + f (1) + 2 f (xk ) + 4 f( )
6n k=1 k=1
2

c'est à dire
(
ISC (f ) =
1  
f (x0 ) + f (x4 ) + 2 f (x1 ) + f (x2 ) + f (x3 )
24
 )
x0 + x1 x1 + x2 x2 + x3 x3 + x4
+ 4 f( ) + f( ) + f( ) + f( )
2 2 2 2
(
1  
= f (0) + f (1) + 2 f (0.25) + f (0.5) + f (0.75)
24
)
 
+ 4 f (.125) + f (0.375) + f (0.625) + f (0.875)
(
1  
= 1 + 2.71828183 + 2 1.06449446 + 1.28402542 + 1.75505466
24
)
 
+ 4 1.01574771 + 1.15099294 + 1.4779042 + 2.15033792 = 1.46272341

Solution de l'e xercice 4 (contrôle nal de janvier 2015)

1. Formules composites et estimations d'erreurs

a Formules composites du point milieu, des trapèzes et de Simpson

n n
C (f ) = b − a X f xk−1 + xk
 
b−aX
IPM = f (x̄k ) (1)
n k=1 2 n k=1
n 
( n−1
)
ITC (f ) =

b − a X f (x k−1 ) + f (x k ) b − a X
= f (a) + f (b) + 2 f (xk ) (2)
n k=1 2 2n k=1
( n−1 n
)
ISC (f ) =
b − a X X x k−1 + x k
f (a) + f (b) + 2 f (xk ) + 4 f( ) (3)
6n k=1 k=1
2

M.H. EL ALJ 4 Décembre 2016

UIT. FSK. Département de Mathématiques Filière SMP : S3

b Estimations d'erreurs du point milieu, des trapèzes et de Simpson

b
(b − a)3
C (f ) ≤
Z
f (x) dx − IPM max |f 00 (x)|
2 x∈[a,b]
(4)
a 24n
Z b
(b − a)3
f (x) dx − ITC (f ) ≤ max |f 00 (x)|
2 x∈[a,b]
(5)
a 12n
Z b
(b − a)4
f (x) dx − ISC (f ) ≤ max |f (4) (x)|
4 x∈[a,b]
(6)
a 120 × 24 × n
2. pour évaluer I à 10−3 par la formule de simpson composite, il sut d'avoir
h4
|ESc | ≤ 10−3 c'est à dire que maxx∈[1,2] |f (4) (x)| ≤ 10−3 .
2880
−24 144x + 48 96
Or f (4)
(x) = 2 3
+ 2 4
− d'oú maxx∈[0,1] |f (4) (x)| ≤= 312. Donc on
(1 + x ) (1 + x ) (1 + x2 )5
devrait avoir
h4 2880 × 10−5 1
× 312 ≤ 10−3 =⇒ h ≤ ( ) 4 ≈ 0.31
2880 312
Donc il sut de prendre h = 14
Pour obtenir la précision voulue. On aura donc besoin de seulement 4 points de subdivision
dans cette méthode.
3. La valeur approchée de I est donnée par
( n−1 n
)
ISC (f ) =
b−a X X xk−1 + xk
f (a) + f (b + 2 f (xk ) + 4 f( )
6n k=1 k=1
2
( 2 )
b−a X
= [f (xk−1 ) + 4 × f (x̄k ) + f (xk )]
6n k=1
       
h 1 1 3
= f (0) + 4 × f +2×f +4×f + f (1)
6 4 2 4
 
1 16 5 16 1
= 1+4× +2× +4× + ≈ 0.86039
12 17 4 25 2
π
4. La valeur exacte de I est
4
l'erreur commise est
b
f (x) dx − ISC (f ) =
Z
π
− 0.86039 < 0.075
a 4
Solution de l'e xercice 5 (contrôle nal de janvier 2015)
1. Formules composites du point milieu
on se place dans l'intervalle [xk−1 , xk ] pour k donné, avec 1 ≤ k ≤ n et on considère le
x + xk
polynôme constant P0,k (f )(x) = f ( k−1 ) Alors:
2
Z xk 
xk−1 + xk
k
IPM (f ) = P0,k (f )(x)dx = [xk − xk−1 ] f
xk−1 2
Z b n x n
C
Z  
k
xk−1 + xk
P0 (f )(x)dx = P0,k (f )(x)dx =
X X
IPM (f ) = [xk − xk−1 ] f
a k=1 xk−1 k=1
2
n   n
b−aX xk−1 + xk b−aX
= f = f (x̄k )
n k=1 2 n k=1
M.H. EL ALJ 5 Décembre 2016
UIT. FSK. Département de Mathématiques Filière SMP : S3

b n
(b − a)3
 
b−aX
Z
xk + xk−1
et on a f (x) dx − f ≤ 2
max |f 00 (x)|
a n k=1 2 24n x∈[a,b]

2. - Formule composite des trapèzes

on se place dans l'intervalle [xk−1 , xk ] pour k donné, avec 1 ≤ k ≤ n et on considère le
polynôme d'interpolation de degré 1 passant par les points xk−1 et xk donné par:
   
f (xk ) − f (xk−1 ) f (xk ) − f (xk−1 ) xk f (xk−1 ) − xk−1 f (xk )
P1,k (f )(x) = (x−xk )+f (xk ) = x+
xk − xk−1 xk − xk−1 xk − xk−1
Alors par intégration:
xk  2
xk − x2k−1
 
f (xk ) − f (xk−1 )
Z
k
IT (f ) = P0,k (f )(x)dx = + (xk f (xk−1 ) − xk−1 f (xk ))
xk−1 xk − xk−1 2
 
f (xk−1 ) + f (xk )
= (xk − xk−1 )
2
Z b n Z xk n 
C

xk − xk−1
P0 (f )(x)dx = P0,k (f )(x)dx =
X X
IT (f ) = {f (xk−1 ) + f (xk )}
a k=1 x k−1 k=1
2
n   ( n−1
)
b − a X f (xk−1 ) + f (xk ) b−a X
= = f (a) + f (b) + 2 f (xk )
n k=1 2 2n k=1

b n
b−aX (b − a)3
Z
et on a f (x) dx − (f (xk ) + f (xk−1 )) ≤ 2
max |f 00 (x)|
a 2n k=1 12n x∈[a,b]

3. nombres de subdivisions pour une évaluation à 10−4 près

Dérivées successives de f (t) = t cos(t) :

0
f (t) = cos(t) − t sin(t)
f ” (t) = −2 sin(t) − t cos(t)
f (3) (t) = −3 cos(t) + t sin(t)
f (4) (t) = 4 sin(t) + t cos(t)

(a) par
Z la méthode du point milieu
1
(1 − 0)3
C (f ) ≤
t cos(t) dt − IPM 2
max | − 2 sin(t) − t cos(t)| ≤ 10−4
0 24n t∈[0,1]
avec max | − 2 sin(t) − t cos(t)| ≤ 3,
t∈[0,1]
3
si et seulement si 2
≤ 10−4
24n
104 √
ce qui entraine que n2 ≥ d'où n ≥ 1250 c'est à dire n ≥ 36.
8
(b) par
Z la éthode de Simpson
1
(1 − 0)4
t cos(t) dt − ISC (f ) ≤ max |4 sin(t) − t cos(t)| ≤ 10−4
4 t∈[0,1]
0 120 × 24n
avec max |4 sin(t) − t cos(t)| ≤ 5,
t∈[0,1]
5
si et seulement si ≤ 10−4
120 × 24n4
5 × 104 √
ce qui entraine que n ≥
4
d'où n ≥ 4 17.36 c'est à dire n ≥ 3.
2880
M.H. EL ALJ 6 Décembre 2016