Método de Simpson
Método de Simpson
Método de Simpson
Partiendo del método del trapecio, vamos a resolver los ejercicios de aproximación al valor de
integral definida.
x dx
2
n=4
En su forma general:
x
1
x dx
2
f (a) + 4 f (a + x) + 2 f (a + 2x) + 4 f (a + 3x) + ... + 4 f (b− x) + f (b)
0
3
b−a
x =
n
1− 0 1
x = =
4 4
1 1
4 f (0) + 4 f (0 + 1 ) + 2 f (0 + 2( 1 )) + 4 f (0 + 3( 1 )) + ... + 4 f (1 − 1 ) + f (1)
0
2
x dx
3 4 4 4 4
1 1
4 f (0) + 4 f ( 1 ) + 2 f ( 2 ) + 4 f ( 3 ) + f (1)
2
x dx
0
3 4 4 4
1
1
x dx 12 0 + 4( 116) + 2( 416) + 4( 916) + 1
2
0
1
1
x dx 12 416 + 816 + 3616 + 1616
2
0
1
1
x dx 12 6416
2
0
1
1
x dx 12 4
2
0
1
4 1
x dx 12 = 3
2
0
1
x dx 0.333333
2
Pero para evitar todo este desarrollo, utilizando la sumatoria de f(x)*patrón en los corchetes, lo
podemos dejar como:
1
1
x dx 12 4
2
0
1
4 1
x dx 12 3
2
0
1
x dx 0.333333
2
0
Ejemplo 2:
4
1
x dx
1
n=6
x
4
1
xdx f (a) + 4 f (a + x) + 2 f (a + 2x) + 4 f (a + 3x) + ... + 4 f (b− x) + f (b)
1
3
b−a
x =
n
4 −1
x =
6
3
x =
6
1
x =
2
x a=1 3/2 4/2 5/2 6/2 7/2 b=4
f(x) 1 2/3 1/2 2/5 1/3 2/7 1/4
Patrón 1 4 2 4 2 4 1
f(x)*patrón 1 8/3 1 8/5 2/3 8/7 1/4 Ʃ(f(x)*patrón)=8.326190
4 1
1
1 xdx 32 1 + 4( 2 3 ) + 2( 1 2) + 4( 2 5 ) + 2( 13 ) + 4( 2 7 ) + 1 4
4
1 1
xdx 6 8.326190
1
4
1 8.326190
xdx
1
6
4
1
xdx 1.387698
1