You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2. ESTIMACIÓN DEL ERROR 141<br />
Definición 7.6. Decimos que una fórmula de cuadratura Q(f) =<br />
exactitud k, si<br />
b<br />
a<br />
n<br />
Ajf(xj) tiene grado de<br />
j=0<br />
p(x) dx = Q(p) para todo polinomio p ∈ Pk y no para Pk+1.<br />
Observación 7.7. Toda fórmula de cuadratura<br />
b<br />
a<br />
f(x) dx ∼ Q(f) =<br />
Es decir, Q(αf + g) = αQ(f) + Q(g) para todo α ∈ IR, f y g funciones.<br />
n<br />
Ajf(xj) es lineal.<br />
En virtud de este resultado, podemos reducir el estudio del grado de exactitud al comportamiento<br />
de la fórmula sobre una base del espacio de polinomios. Esto queda expresado en la siguiente<br />
observación.<br />
Observación 7.8. Una fórmula de cuadratura Q tiene grado de exactitud k si y solo si es exacta<br />
para la base de Pk, B = {1, x, x2 , . . . , xk } y no lo es para el polinomio xk+1 . Esto es, la igualdad<br />
b<br />
x<br />
a<br />
m n<br />
dx = Ajx<br />
j=0<br />
m j<br />
debe verificarse para todo m = 0, . . . , k y no para m = k + 1.<br />
Además, gracias al Lema 7.3, una fórmula de cuadratura tiene grado de exactitud k independientemente<br />
del intervalo [a, b] para el cual está calculada.<br />
Ejemplo 7.9. Se quiere calcular el grado de exactitud de la fórmula de Simpson cerrada.<br />
Es claro que las fórmulas de Simpson, tanto abiertas como cerradas, son exactas para polinomios<br />
de P2. Luego, por Observación 7.8, basta ver qué sucede con x 3 , x 4 , . . ., y esto puede hacerse,<br />
sin perder generalidad, en el intervalo [−1, 1]. Tenemos<br />
1<br />
−1 1<br />
−1<br />
x 3 dx = 0 = 2<br />
6 [(−1)3 + 4(0) 3 + (1) 3 ],<br />
x 4 dx = 2 2 1<br />
= =<br />
5 3 3 [(−1)4 + 4(0) 4 + (1) 4 ].<br />
Luego, la fórmula de Simpson cerrada tiene grado de exactitud k = 3. Lo mismo sucederá para<br />
la fórmula abierta.<br />
En lo que sigue notaremos por I(f) =<br />
2. Estimación del error<br />
b<br />
a<br />
j=0<br />
f(x) dx. Una vez obtenida una regla de cuadratura<br />
Q(f) surge la pregunta de qué se puede decir sobre el error cometido R(f) = I(f) − Q(f).