FST A.

U:2024-2025
TD IGL3+IDS3.

Série 2: Correction d’exercice 3 et 4

Exercice 3. Soit f une fonction continue et de dérivée première continue sur l’intervalle [0, 1]. Pour
l’intégration de f sur [0, 1], on considère la formule de quadrature :
Z 1
(1) f (t)dt = A1 · f (0) + A2 · f ′ (0) + A3 · f (α) + E(f ),
0

où f′ est la dérivée première de f, α ∈]0, 1[, A1 , A2 , A3 ∈ R et E(f ) est le terme d’erreur.

1) Calculer A1 , A2 , A3 et α pour que la formule d’intégration (1) ait le plus grand degré de précision
possible.
2) Si n est le degré de précision trouvé pour la formule (1), calculer E(f ) pour f (x) = xn+1 .
3) Par un changement de variable convenable (x = (b − a)t + a), construire à partir de la formule
Z b
d’intégration (1), une deuxième formule pour calculer f (x)dx.
a Z 10
x
4) En utilisant cette deuxieme formule, calculer une valeur approchée pour exp(− )dx. Calculer
0 5
la valeur exacte de cette intégrale puis donner l’erreur absolue commise.
5) On utilise maintenant la methode de Simpson composée pour calculer cette intégrale, pour cela on
divise l’intervalle [0, 10] en n sous intervalles égaux de longueur h. Quelle est la valeur minimale de n
qui donne une erreur d’integration inferieure ou égale à celle obtenue en 4).
Correction.

1. Z 1
f (t)dt = A1 · f (0) + A2 · f ′ (0) + A3 · f (α) + E(f ).
0
Z 1
f (t) = 1, E(f ) = 0 ⇐⇒ A1 · f (0) + A2 · f ′ (0) + A3 · f (α) =
1dt ⇐⇒ A1 + A3 = 1 (1)
Z 1 0
′ ′ 1
f (t) = t, f (t) = 1 E(f ) = 0 ⇐⇒ A1 ·f (0)+A2 ·f (0)+A3 ·f (α) = tdt ⇐⇒ A2 +A3 α = (2)
0Z 2
1
1
f (t) = t2 , f ′ (t) = 2t E(f ) = 0 ⇐⇒ A1 ·f (0)+A2 ·f ′ (0)+A3 ·f (α) = t2 dt ⇐⇒ A3 α2 = (3)
Z0 1 3
1
3 ′ 2 ′
f (t) = t , f (t) = 3t E(f ) = 0 ⇐⇒ A1 ·f (0)+A2 ·f (0)+A3 ·f (α) = t3 dt ⇐⇒ A3 α3 = (4)
0 4
1 α3 1 3 1 16 16
(3) =⇒ A3 = , (4) =⇒ = =⇒ α = , donc A3 = · =
3α2 3α2 4 4 3 9 27
11 1 1 16 3 1
(1) =⇒ A1 = 1 − A3 = , (2) =⇒ A2 = − A3 α = − = .
27 2 2 27 4 18
La formule de quadrature s’écrit :
Z 1  
11 1 ′ 16 3
(1) f (t)dt = · f (0) + · f (0) + ·f + E(f ).
0 27 18 27 4
Cette formule est exacte pour les polynômes de dégré ⩽ 3.
2. Déterminons le degré de précision
Z 1
16 34
 
1 11 1 16 3 3
f (t) = t4 , I = f (t)dt = , J(f ) = · f (0) + · f ′ (0) + ·f = · 4 =
0 5 27 18 27 4 27 4 16
1 3 1
E(f ) = I − J(f ) = − = ̸= 0, La méthode est de degré de précision 3
5 16 80

1
 Z b
3. Pour f (x)dx, on effectue le changement de variable x = (b − a)t + a =⇒ dx = (b − a)dt, x =
a
a =⇒ t = 0, x = b =⇒ t = 1
Z b Z 1 Z 1
f (x)dx = f ((b − a)t + a)(b − a)dt = g(t)dt,
a 0 0

g(t) = (b − a)f ((b − a)t + a), g ′ (t) = (b − a)2 f ′ ((b − a)t + a)

Z b  
11 1 ′ 16 3
f (x)dx = · g(0) + · g (0) + ·g + E(g)
a 27 18 27 4
Z b  
11 1 2 ′ 16 3
f (x)dx = (b − a) · f (a) + (b − a) · f (a) + (b − a) · f (b − a) + a + R(f ) (2)
a 27 18 27 4
Z 10  x  x 1  x
4. exp − dx, [a, b] = [0, 10], f (x) = exp − , f ′ (x) = − exp −
0 5 5 5 5
La formule (2) donne
Z 10  x  
110 100 ′ 160 15
exp − dx ≃ f (0) + f (0) + f
0 5 27 18 27 2
Z 10  x  
110 20 160 3
exp − dx ≃ − + exp − ≃ 4.2852
0 5 27 18 27 2
Valeur exacte :
Z 10  x h  x i10
exp − dx = −5 exp − = −5e−2 + 5 ≃ 4.3233
0 5 5 0

L′ erreur : |R(f )| = |4.3233 − 4.2852| ≃ 0.0381

(b − a)5 M4
5. L’erreur de Simpson : |Rn (f )| ⩽ , M4 = sup |f 4 (x)|
2880n4 x∈[0,10]

x 1 x 1 x 1 x 1
f (x) = e− 5 , f ′ (x) = − e− 5 , f ′′ (x) = 2 e− 5 , · · · , f (4) (x) = 4 e− 5 , M4 =
5 5 5 54
Pour que l’erreur soit inférieur à 0.0381 il suffit que

(b − a)5 M4 4 105
⩽ 0.0381 =⇒ n ⩾ ,
2880n4 0.0381 ∗ 2880 ∗ 54
n ⩾ 1.458 On prend n = 2
Exercice 4. Soit la formule d’intégration numérique sur [0, 1] suivante :
Z 1
f (t)dt = αf (0) + βf (1) + E(f ) (1).
0

1. Déterminer α et β pour que cette formule soit exacte : E(f ) = 0, sur l’espace des polynômes le
plus grand possible.
2. On suppose que f est de classe C 2 sur [0, 1]. Montrer qu’il existe η ∈ [0, 1] tel que

f ′′ (η)
E(f ) = −
12
3. Par un changement de variable affine, construire à partir de la formule d’intégration (1), une
Z b
deuxième formule pour calculer f (x)dx.
a

2
Correction.
1  1
Z
Pour f (x) = 1, E(f ) = 0 ⇐⇒ αf (0) + βf (1) = 1dx = x 0 = 1 ⇐⇒ α + β = 1
Z 1 0
 x2 1 1 1
• f (x) = x, E(f ) = 0 ⇐⇒ αf (0) + βf (1) = xdx = 0
= ⇐⇒ β =
0 2 2 2
Alors on trouve α = 21 , β = 21 . La formule de quadrature est de la forme :
Z 1
1 1
f (t)dt = f (0) + f (1) + E(f ) (1).
0 2 2

Cette formule est d’ordre au moins 1 : exacte pour les polynôme de degré ⩽ 1.
Déterminons le degréZ 1 de précision :
1 1 1 1
• Pour f (x) = x2 ; x2 dx = ; et f (0) + f (1) = =⇒ E(f ) ̸= 0.
0 3 2 2 2
La formule de quadrature est de dégré de précision 1.
2) Soit P le polynôme qui interpole f aux points x0 = 0, x1 = 1, deg(P ) = 1, puisque la formule de
quadrature est exacte pour les polynômes de degré ⩽ 1, alors
Z 1
1 1 1 1
P (t)dt = P (0) + P (1) = f (0) + f (1)
0 2 2 2 2
Z 1 Z 1 Z 1
Alors E(f ) = f (t)dt − P (t)dt = (f (t) − P (t))dt.
0 0 0
f (2) (ξt )) f (2) (ξt ))
Or f (t) − P (t) = (t − 0)(t − 1) = t(t − 1), or t(t − 1) ne change pas de signe sur [0, 1]
2! 2
parsuite, d’après la formule de la moyenne ∃ une constante η ∈ [0, 1] tel que :
1
f (2) (ξt )) f ′′ (η) 1
f ′′ (η)
Z Z
t(t − 1)dt = t(t − 1)dt = − ,
0 2 2 0 12
Z b
3) Une formule pour f (x)dx. On effectue un changement de variable affine :
a

ϕ : [0, 1] −→ [a, b]
t −→ x = αt + β
 
t = 0 ⇐⇒ x = a ⇐⇒ a = β (1) β=a
=⇒
t = 1 ⇐⇒ x = b ⇐⇒ b = α + β (2) α=b−a
On effectue le changement de variable x = (b − a)t + a =⇒ dx = (b − a)dt :
Z b Z 1 Z 1
f (x)dx = f ((b − a)t + a) (b − a)dt = (b − a) g(t)dt,
a 0 0

où g(t) = f ((b − a)t + a) , on utilise la formule de quadrature (1), on obtient :

Z b  
1 1
f (x)dx = (b − a) g(0) + g(1) + R(g)
a 2 2
b
(b − a)
Z
f (x)dx = (f (a) + f (b)) + R2 (f ).
a 2