Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
148 7. INTEGRACIÓN NUMÉRICA<br />
Entonces, si f ∈ C k+1 [a, b] y h = máxj(xj+1 − xj) se tiene que<br />
|R(f)| = |I(f) − Q(f)| ≤<br />
Demostración. Por el Teorema 7.10 sabemos que<br />
(1 + M)(b − a)<br />
(k + 1)!<br />
h k+1 f k+1 ∞ .<br />
|Rj(f)| ≤ (1 + M)(xj+1 − xj) k+2<br />
f<br />
(k + 1)!<br />
k+1 ∞ .<br />
En consecuencia, acotando (xj+1 − xj) k+1 por h k+1 se tiene<br />
n−1 <br />
|R(f)| ≤ |Rj(f)| ≤<br />
j=0<br />
≤<br />
n−1<br />
(1 + M) <br />
(k + 1)!<br />
j=0<br />
n−1<br />
(1 + M) <br />
(k + 1)!<br />
j=0<br />
(xj+1 − xj) k+2 f k+1 ∞<br />
(xj+1 − xj)h k+1 f k+1 ∞ ,<br />
pero n−1<br />
j=0 (xj+1 − xj) = b − a, lo que concluye la demostración.<br />
El teorema anterior puede aplicarse a las reglas compuestas usuales (tales como Trapecios y<br />
Simpson) para obtener una estimación del error en términos de h = máxj(xj+1 − xj). Sin<br />
embargo, un análisis mas detallado permite dar una expresión más precisa del error en esos<br />
casos particulares. Para esto vamos a necesitar el siguiente lema.<br />
Lema 7.18. Sea g ∈ C[a, b] y sean {a0, . . . , ak} constantes con el mismo signo y {t0, . . . , tk} ∈<br />
[a, b], entonces se tiene<br />
para algún η ∈ [a, b].<br />
k<br />
ajg(tj) = g(η)<br />
j=0<br />
Demostración. Sea m = mín g(x) y M = máx g(x) en [a, b]. Podemos suponer que aj ≥ 0 para<br />
todo j = 1, . . . , k, luego<br />
k<br />
j=0<br />
para cada j : m ≤ g(tj) ≤ M<br />
(aj ≥ 0) m aj ≤ ajg(tj) ≤ M aj<br />
(sumando)<br />
k<br />
m aj ≤<br />
k<br />
g(tj)aj ≤<br />
k<br />
M<br />
Ahora, definimos la función G : [a, b] → IR,<br />
j=0<br />
G(x) = g(x)<br />
j=0<br />
k<br />
aj,<br />
j=0<br />
aj<br />
j=0<br />
aj