You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
140 7. INTEGRACIÓN NUMÉRICA<br />
Luego, 1<br />
−1<br />
<br />
2<br />
p(x) dx = 3 3a0 + a2<br />
9<br />
8<br />
7<br />
6<br />
5<br />
4<br />
3<br />
2<br />
1<br />
(a+h, f(a+h))<br />
(a+2h, f(a+2h))<br />
(a+3h, f(a+3h))<br />
0<br />
−1 −0.5 0 0.5 1 1.5 2<br />
Figura 7.4. Regla de Simpson abierta<br />
= 2<br />
3<br />
<br />
1 2f(− 2 ) − f(0) + 2f(1 2 ) . Al pasar a un intervalo [a, b] por<br />
medio del Lema 7.3 y escribiendo la fórmula en términos del paso h = b−a<br />
4 obtenemos,<br />
b<br />
a<br />
f(x) dx ∼ 4h<br />
<br />
2f(a + h) − f(a + 2h) + 2f(a + 3h) .<br />
3<br />
Ejemplo 7.5. Consideremos nuevamente la función f(x) = x 3 − 4x + 4 en el intervalo [−1, 2].<br />
Queremos hallar una aproximación de la integral de f en dicho intervalo por medio de la regla<br />
de Simpson abierta.<br />
Como h = b−a<br />
4<br />
2<br />
−1<br />
= 3<br />
4<br />
, entonces a + h = − 1<br />
4<br />
x 3 − 4x + 4 dx ∼ 1 2f(− 1<br />
4<br />
que vuelve a ser un cálculo exacto.<br />
, a + 2h = 1<br />
2<br />
5<br />
, a + 3h = 4 , así tenemos<br />
) − f(1 ) + 2f(5<br />
2 4 ) = 2 319 17<br />
−<br />
64 8<br />
624 39<br />
+ 261 = =<br />
64 64 4 ,<br />
Es claro que la fórmula de Simpson es exacta para polinomios de P2 y que la fórmula de<br />
Trapecios lo es para polinomios de P1. En general, una regla como la dada en (7.1), construida<br />
interpolando en n + 1 puntos, es exacta para polinomios de Pn, sin embargo, veremos que para<br />
ciertas elecciones de puntos la regla resulta exacta también para polinomios de Pk para algún<br />
k > n. Por otra parte, veremos luego que el grado de polinomios para el cual una regla es exacta<br />
se relaciona con el error cometido al usar dicha regla. Esto motiva la siguiente definición.