Unha mellor aproximación que a conseguida co método dos rectángulos ou co dos trapecios lógrase co método de Simpson, que consiste en aproximar a función mediante un polinomio de segundo grao que pase polos extremos e polo punto medio da función en cada subintervalo.
Para integrar a función f(x) no intervalo
facendo n subintervalos de amplitude
, para cada subintervalo
, chamando ao punto medio
, temos que a parábola está definida por: (utilizamos o polinomio de interpolación de Lagrange)
![[Graphics:Images/index_gr_5.gif]](imaxe/index_gr_5.gif)
integrando a parábola entre
e
, e substituíndo
por h obtemos
![[Graphics:Images/index_gr_10.gif]](imaxe/index_gr_10.gif)
Algoritmo
![[Graphics:Images/index_gr_11.gif]](imaxe/index_gr_11.gif)
![[Graphics:Images/index_gr_12.gif]](imaxe/index_gr_12.gif)
A integral aproximada é: 2.0000067844418011042
![[Graphics:Images/index_gr_16.gif]](imaxe/index_gr_16.gif)