Méthode d’intégration de Simpson : introduction

La méthode de Simpson, comme d’autres techniques d’approximation d’intégrales, fait appel à des concepts mathématiques avancés. Dans cet article, nous allons explorer cette méthode qui se distingue par son utilisation de courbes paraboliques pour estimer des intégrales.

Comprendre la méthode

Considérons une fonction ( f ) définie sur l’intervalle ([a;b]). Pour appliquer la méthode de Simpson, cet intervalle est divisé en un nombre pair d’intervalles. Sur chaque intervalle, nous choisissons trois points :

• ( (x{2k}; f(x{2k})) )
• ( (x{2k+1}; f(x{2k+1})) ), où ( m = \frac{x{2k} + x{2(k+1)}}{2} )
• ( (x{2(k+1)}; f(x{2(k+1)})) )

L’idée est de construire une parabole passant par ces trois points, et l’aire sous cette parabole servira d’approximation de l’intégrale sur l’intervalle considéré.

Formule de la méthode de Simpson

Avant de plonger dans la formule, il est essentiel d’évoquer le polynôme d’interpolation de Lagrange. Pour ( n + 1 ) points ((x_k; f(x_k))), la formule de Lagrange est donnée par :

[
L(x) = \sum_{i=0}^n f(xi) \prod{\substack{k=0\k \neq i}}^n \frac{x – x_k}{x_k – x_i}
]

En appliquant cela à trois points, nous obtenons un polynôme qui nous permet d’estimer l’intégrale.

Approximations et calculs

En posant ( t = x – x_k ), le polynôme devient :

[
P_k(t) = \frac{(t-h)(t-2h)}{2h^2}f(xk) – \frac{t(t-2h)}{h^2}f(x{k+1}) + \frac{t(t-h)}{2h^2}f(x_{k+2})
]

L’intégrale sur l’intervalle peut alors être calculée, conduisant à l’approximation suivante :

[
\inta^b f(x)~\text{d}x \approx \frac{b-a}{6n}\sum{k=0}^{n-1}\left( f(xk) + 4f(x{k+1}) + f(x_{k+2}) \right)
]

Implémentation en Python

Pour ceux qui souhaitent mettre en pratique cette méthode, voici un exemple d’implémentation en Python :

python
from math import log, exp

def f(x):
return 5 log(x) exp(-x/2)

def simpson(a, b, n):
S = 0
h = (b – a) / (2 n)
for k in range(n):
x0 = a + 2 k h
x1 = x0 + h
x2 = x0 + 2 h
S += f(x0) + 4 * f(x1) + f(x2)

return h * S / 3

Pour utiliser cette fonction, il suffit d’appeler simpson(1, 10, 10) pour obtenir une approximation de l’intégrale sur l’intervalle de 1 à 10.

Conclusion

La méthode de Simpson est une technique puissante pour l’approximation d’intégrales, particulièrement utile dans les domaines nécessitant des calculs rapides et efficaces. Pour ceux qui souhaitent approfondir leurs compétences, il est recommandé de comparer différentes méthodes d’intégration et d’anticiper les coûts associés à leur utilisation. En intégrant des outils numériques, il est possible d’éviter les frais liés à des calculs manuels fastidieux.

Pour réserver votre place dans des cours ou des ateliers sur les mathématiques appliquées, vous pouvez explorer des options en ligne sur des plateformes telles que GetYourGuide ou Booking.