Université Sultan Moulay Slimane, ENSA Khouribga, Année
Universitaire 2022-23, API 2, S2, TD3 : Analyse
Numérique-différentiation et intégration Numériques
Exercice 1 :
1. A partir d’un polynôme de degré 2 passant par les points (x0 , f (x0 )), (x1 , f (x1 ))
et (x2 , f (x2 )), obtenir les formules aux différences avant, centrée et arrière d’ordre
2 pour le calcul de f 0 (x). Déduire également les termes d’erreur.
2. On définit
−y(x + 2h) + 4y(x + h) − 3y(x)
S(h) =
2h
a) Montrer que
y 0 (x) − S(h) = c1 h2 + c2 h3 + c3 h4 + ...
b) Calculer y 0 (0.398) en utilisant l’approximation S(h) et le tableau suivant:
x 0.398 0.399 0.400 0.401 0.402
f(x) 0.408591 0.409671 0.410752 0.411834 0.412915
Donner une estimation de l’erreur.
3. En vous servant des développements de Taylor appropriés, donner l’ordre de
précision de l’approximation
f (x + 3h) − 3f (x + 2h) + 3f (x + h) − f (x)
f (3) (x) ≈
h3
Exercice 2 :
On suppose que la fonction
2
f (x) =
1 + x2
est connue pour le trois points de collocation x1 = −1, x2 = 0, et x3 = 1.
Interpoler la fonction avec une spline cubique naturelle et approximer l’intégrale
Z 1
2
2
=π
−1 1 + x
par l’intégrale de la spline. Le résultat est-il plus précis que le résultat de la méthode
de Simpson 1/3 simple?
Exercice 3 :
La longueur L d’une ellipse avec ses deux axes a et b est donnée par la formule L =
2 2
4aE(m) où m = a a−b
2 et
Z π
2 1
E(m) = (1 − m sin2 Φ) 2 .
0
La fonction E(m) est un intégral elliptique, certains de ses valeurs sont données par
le tableau suivant:
m 0 0.1 0.2 0.3 0.4 0.5
E(m) 1.57080 1.53076 1.48904 1.44536 1.39939 1.35064
1. Calculer L en utilisant une interpolation quadratique du tableau ci-dessus pour
a = 5 et b = 4.
2. Calculer L en utilisant le méthode de Romberg à E(m), de sorte qu’une valeur
de Romberg est obtenue avec une erreur inférieure à 5 × 10−5 .
Page 2
