Université de Caen 1er semestre 2024–2025
UFR Sciences Portail MIM
L3 Maths Calcul formel
TP 6 : Intégration numérique
Pour calculer de manière approchée l’intégrale de f sur [a, b], on approxime l’aire sous la courbe par celle
de n rectangles de largeur (b − a)/n et de hauteur donnée par la valeur de f à gauche de chaque intervalle. On
note In0 (f, a, b) la valeur obtenue par cette méthode des rectangles, qui est également la ne somme de Riemann à
gauche avec pas constant pour cette intégrale. On peut procéder de même en remplaçant les rectangles par des
trapèzes, cf les figures ci-dessous. On note In1 (f, a, b) la valeur approchée obtenue par la méthode des trapèzes.
Figure 1 – Rectangles Figure 2 – Trapèzes
Dans la méthode des rectangles (resp. des trapèzes) on approche f par une fonction constante par morceaux
(resp. affine par morceaux). De manière analogue, on peut approcher f sur chaque sous-intervalle [xi , xi+1 ]
par son polynôme d’interpolation de Lagrange aux trois points xi , (xi + xi+1 )/2 et xi+1 , qui est de degré 2.
C’est la méthode de Simpson 1 . En calculant l’intégrale de ces polynômes sous chaque sous-intervalle, on obtient
la formule suivante pour les valeurs approchées In2 (f, a, b) correspondantes, où les xi sont les n + 1 points
régulièrement espacés dans [a, b] :
n−1 n−1
!
2 b−a X X x +x
i i+1
In (f, a, b) = f (x0 ) + f (xn ) + 2 f (xi ) + 4 f 2 .
6n i=1 i=0
Exercice 1.
a. Écrire l’expression mathématique de In0 (f, a, b), sous la même forme que l’expression de In2 (f, a, b) ci-
dessus.
b. Écrire une procédure Rectangles(f,a,b,n) qui calcule la valeur de In0 (f, a, b). On pourra utiliser une
procédure auxiliaire qui fournit la liste des n + 1 points a = x0 < x1 < . . . < xn = b régulièrement
espacés dans [a, b].
c. Tester la procédure Rectangles avec la la fonction f (x) = sin(x) sur l’intervalle 0, π2 , pour n =
1, 3, 10, 30, 50. Comparer avec la valeur exacte I de l’intégrale.
d. Écrire l’expression mathématique de In1 (f, a, b).
e. Écrire une procédure Trapèzes(f,a,b,n) qui calcule la valeur de In1 (f, a, b).
f. Écrire une procédure Simpson(f,a,b,n) qui calcule la valeur de In2 (f, a, b).
g. Répéter les tests avec les procédures Trapèzes et Simpson.
Quelle méthode semble avoir la meilleure convergence ?
h. Remplacer l’intervalle 0, π2 par [0, π]. Que constate-t-on ?
Expliquer en faisant un schéma dans le cas n = 4.
1. Thomas Simpson, mathématicien anglais du 18e siècle. L’histoire ne dit pas s’il aimait les donuts.
Exercice 2. On souhaite étudier l’erreur relative E p (n) = |I − Inp |/|I| commise par chacune des trois méthodes
(p = 0, 1, 2), en fonction du nombre d’intervalles n.
a. On dispose d’une classe Intégrale avec une méthode approx(m) qui calcule les valeurs approchées sur
m sous-intervalles. Lors de l’instantiation de la classe on fournit la fonction à intégrer, l’intervalle, ainsi
que la valeur exacte de l’intégrale et la méthode à utiliser pour le calcul approché.
Ajouter à cette classe une méthode erreur(m) qui calcule l’erreur relative E p (m) commise par le calul
approché sur m sous-intervalles.
b. On étudie la fonction g(x) = 4/(1 + x2 ) sur [0, 1].
(i) Tracer sur un même graphique les erreurs E p (n) en faisant varier n dans l’ensemble {2k ; k = 0, . . . , 6}.
L’erreur semble décroître très vite quand n croît, mais il n’est pas évident de voir si la décroissance
est polynômiale ou exponentielle.
(ii) Modifier le script pour utiliser une échelle log-log : [Link](’log’), [Link](’log’). On fait
maintenant varier n dans l’ensemble {2k ; k = 0, . . . , 15}. Qu’observe-t-on ? Qu’est ce que cela signifie
pour le comportement asymptotique de E p (n) ?
(iii) Afficher les quotients ln(E p (n))/ ln(n), k
π pour n dans l’ensemble {2 ; k = 5, . . . , 15}.p Exécuter également
les calculs avec f (x) = sin(x) sur 0, 2 . Conjecturer un équivalent simple pour E (n), à une constante
multiplicative près, pour chacune des trois méthodes.
√
c. On étudie maintenant la fonction x sur l’intervalle [0, 1].
(i) Afficher le graphe des erreurs relatives pour les méthodes des trapèzes et de Simpson, toujours en
échelle log-log. On ajoutera sur le graphe les courbes n−2 , n−4 en pointillés (linestyle=’dashed’)
pour comparer. Que constate-t-on ?
(ii) Recommencer sur les intervalles [, 1], en prenant successivement = 0.0001, 0.001, 0.01, 0.1. À quoi
le phénomène observé peut-il être dû ?
Exercice 3. Les fonctions de Bessel de première espèce sont définies, pour n ∈ N, x ∈ R, par la formule
1 π
Z
Jn (x) = cos(nt − x sin(t))dt.
π 0
a. Tracer sur une même figure les graphes de J0 , J1 , J2 , J3 sur l’intervalle [0, 30].
b. Conjecturer les valeurs C et α dans le comportement asymptotique Jn (x) ∼ Cxα pour J1 , J2 et J3 .