Intégration Numérique
Oualid Oukrizi
Informatique
CPGE REDA SLAOUI AGADIR
2025/2026
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 1 / 16
Introduction
Approcher une intégrale définie :
Z b
I= f (x ) dx
a
Calcul utile quand la solution analytique est impossible.
Lien avec la discrétisation de l’aire sous une courbe.
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 2 / 16
Conditions de validité du calcul intégral numérique
Hypothèses nécessaires
Pour que le programme de calcul intégral fonctionne correctement, il faut
que :
La fonction f (x ) soit continue sur l’intervalle [a, b].
f (x ) soit bornée sur [a, b] (pas de valeurs infinies).
L’intervalle d’intégration [a, b] soit bien défini avec a < b.
Le nombre de subdivisions n choisi pour la méthode numérique
(rectangles, trapèzes, Simpson, . . .) soit suffisamment grand pour
garantir une bonne précision.
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 3 / 16
Méthode des rectangles
Il n’est pas toujours possible de calculer une primitive pour une fonction
Z b
f : [a, b] → R. On va déterminer des valeurs approchées de f (t)dt.
a
Les trois méthodes d’approximation que l’on va étudier sont toutes basées
sur le même principe :
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 4 / 16
On divise l’intervalle [a, b] en n sous-intervalles en posant
xk = a + k × b−a n pour 0 ≤ k ≤ n. Alors x0 = a et xn = b et chaque
sous-intervalle [xk , xk+1 ] est de longueur constante xk+1 − xk = b−a
n .
Figure – Illustration de la subdivision de l’intervalle
Sur chaque sous-intervalle [xk , xk+1 ], on approche l’aire sous la courbe par
l’aire d’une figure géomé trique simple.
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 5 / 16
La méthode des rectangles (à gauche) consiste à approcher l’aire sous
la courbe par l’aire de rectangles. La hauteur de chaque rectangle est la
valeur à gauche de f sur le sous-intervalle (voir les figures ci-dessous).
Pour un intervalle élémentaire, cela revient à approcher xxkk+1 f (t) dt par
R
(xk+1 − xk )f (xk ).
Figure – Illustration de la méthode des rectangles à gauche
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 6 / 16
Implémentation Python : Rectangles
Rectangles à droite :
Z xk+1
f (x )dx ≈ h · f (xk+1 )
xk
1 def rect_gauche (f , a , b , n ) :
2 h = (b - a) / n
3 somme = 0
4
5 for k in range ( n ) :
6 somme = somme + f ( a + k * h )
7
8 return somme * h
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 7 / 16
La méthode des rectangles (à droit) consiste à approcher l’aire sous la
courbe par l’aire de rectangles. La hauteur de chaque rectangle est la
valeur à gauche de f sur Rle sous-intervalle. Pour un intervalle élémentaire,
cela revient à approcher xxkk+1 f (t) dt par (xk+1 − xk )f (xk+1 ).
Figure – Illustration deIntégration
O. Oukrizi (CPGE REDA SLAOUI) la méthode
Numériquedes rectangles à droite
2025/2026 8 / 16
Implémentation Python : Rectangles
1 def rect_droit (f , a , b , n ) :
2 h = (b - a) / n
3 somme = 0
4
5 for k in range ( n ) :
6 somme = somme + f ( a + ( k +1) * h )
7
8 return somme * h
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 9 / 16
R xi+1
La méthode du point milieu : on approche xi f (t) dt par
xi +xi+1
(xi+1 − xi ) · f . 2
Géométriquement, cela signifie qu’on approche l’intégrale de f par l’aire
des rectangles hachurés en bleu :
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 10 / 16
Implémentation Python : Rectangles
Z xk+1
xk + xk+1
f (x )dx ≈ h · f
xk 2
1 def rect _poin t_mili eu (f , a , b , n ) :
2 h = (b - a) / n
3 somme = 0
4
5 for k in range ( n ) :
6 somme = somme + f ( a + k *( h /2) )
7
8 return somme * h
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 11 / 16
Méthode des trapèzes
On approche l’aire sous la courbe d’un intervalle élémentaire [xi , xi+1 ], par
l’aire d’un trapèze :
Z xi+1
(xi+1 − xi )
f (t) dt ≈ · [f (xi ) + f (xi+1 )]
xi 2
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 12 / 16
Implémentation Python : Trapèzes
1 def trapezes (f , a , b , n ) :
2 h = (b - a) / n
3 return ( h /2) * ( f ( a ) + 2* sum ( f ( a + k * h ) for k in
range (1 , n ) ) + f ( b ) )
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 13 / 16
Fonction quad - Intégration numérique
Syntaxe et utilisation
Définition
La fonction quad est une fonction d’intégration numérique qui utilise
l’algorithme quadrature adaptative pour calculer une approximation
numérique d’une intégrale.
Syntaxe Python (SciPy)
1 from scipy import integrate
2
3 result , error = integrate . quad (f , a , b , args =() ,
full_output =0 , ...)
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 14 / 16
Paramètres
f : fonction à intégrer
a, b : bornes d’intégration
args : arguments supplémentaires pour f
full_output : afficher info détaillées
Valeurs de retour
result : valeur approximative
error : estimation de l’erreur
info : informations détaillées (si full_output=1)
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 15 / 16
Exemples d’utilisation
Applications pratiques
Exemple 1 : Intégrale simple
1 import numpy as np
2 from scipy import integrate
3
4 # D é finition de la fonction
5 def f ( x ) :
6 return x **2 + 3* x + 2
7
8 # Calcul de l ’ int é grale de 0 à 1
9 result , error = integrate . quad (f , 0 , 1)
0 print ( f " R é sultat : { result }")
1 print ( f " Erreur estim é e : { error }")
Sortie : Résultat : 3.8333, Erreur estimée : 4.25e-14
O. Oukrizi (CPGE REDA SLAOUI) Intégration Numérique 2025/2026 16 / 16