Calculo Aproximado de Integrales Definidas
Calculo Aproximado de Integrales Definidas
Calculo Aproximado de Integrales Definidas
Sea y = f(x) una función continua en un intervalo [a, b] de la recta real. Estudiemos
distintas aproximaciones de I f x dx .
b
Sea Pn = {a = x 0 , x 1 , ..., x n = b} una partición del intervalo [a, b]. Sea h la amplitud de
b-a
cada subintervalo de la partición, h x i x i 1 , i 1, ..., n , es decir, h .
n
Imagen 1
Las expresiones (1) y (2) se llaman fórmulas de los rectángulos. El error cometido al
aplicarlas disminuye a medida que aumenta n.
Imagen 2
La suma de las áreas de estos trapecios constituye una aproximación de I mejor que la
dada por las fórmulas (1) y (2):
y 0 y1 y y2 y yn h
I f x dx h y 0 y n 2 y1 ... y n 1
b
h 1 h ... n 1
a 2 2 2 2
h h n 1
f (x 0 ) f (x n ) 2 f (x1 ) ... f (x n 1 ) f (x 0 ) f (x n ) 2 f (x k )
2 2 k 1
h n 1
2
f (a) f (b) 2
k 1
f (a kh)
(3)
La estimación por el método de los trapecios es exacta para polinomios de primer grado
(tienen nula la segunda derivada). El factor clave como indicador del error es h 2 . Si
reducimos h a la mitad, el error se reduce en un factor 4.
Imagen 3
La parábola que pasa por A 0 x 0 , y 0 , A 1 x 1 , y1 y A 2 x 2 , y 2 tiene de ecuación
y a0 a1x a 2 x 2 , pudiendo considerar x 0 0, x 1 h y x 2 2h .
Relacionemos a 0 , a1 y a 2 con y 0 , y1 e y 2 para sustituir en la expresión anterior:
A i 1 parábola x i 1 0 yi 1 a 0
A i parábola x i h yi a 0 a1h a 2 h 2
A i 1 parábola x i 1 2h yi 1 a 0 a1 2h a 2 2h
2
6a 0 6a1h 8a 2 h a 0 4 a 0 a1h a 2 h a 0 a1 2h a 2 4h
h h 2 2
2
3 3 y
i1 yi yi1
h
4y1 y0 y 2
3
Análogamente se obtendrían:
f x dx 3 y 4 y 3 y 4 , …, f x dx y n 2 4 y n 1 y n
x4 h xn h
x2
2 x n -2 3
h
I f x dx y 0 4y1 y 2 y 2 4y3 y 4 ... y n 2 4y n 1 y n
b
a 3
h
y0 4 y1 y3 ... y n 1 2 y 2 y 4 ... y n 2 y n
3
h (n 2) / 2 (n 2) / 2
3
f (a) 4 f a (2k 1)h 2 f a 2kh f (b) (4)
k 0 k 1
El método de Simpson es exacto para todo polinomio de grado tres (tiene nula la
derivada cuarta).
Cabe esperar que este método sea, en general, más preciso que el de los trapecios, pues,
al ir disminuyendo h, el factor h 4 tiende a cero más rápidamente que h 2 . Al reducir h en
un factor 2, el error se reduce en un factor 16.