You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
72 4. RESOLUCIÓN DE ECUACIONES NO LINEALES<br />
Observemos que si f ′ (r) = 0 entonces f ′ (xn) = 0 para xn cercano a r (esto lo justificaremos con<br />
más precisión después).<br />
Usando el desarrollo de Taylor de orden 2 centrado en la raíz r se tiene,<br />
0 = f(r) = f(xn) − (xn − r)f ′ (xn) + 1<br />
2 (xn − r) 2 f ′′ (ξ)<br />
donde ξ es un valor intermedio entre xn y r. Entonces<br />
Reemplazando en la igualdad 4.1 queda<br />
enf ′ (xn) − f(xn) = 1<br />
2 f ′′ (ξ)e 2 n<br />
en+1 = 1 f<br />
2<br />
′′ (ξ)<br />
f ′ (xn) e2 n<br />
Con todo esto podemos demostrar el siguiente teorema.<br />
Teorema 4.5. (de convergencia) Si r es un cero simple de f (i.e. f ′ (r) = 0) y sea I =<br />
[r − α, r + α] un intervalo tal que |f ′ (x)| ≥ δ > 0 y |f ′′ (x)| ≤ M en I. Entonces,<br />
Existe ε > 0 tal que Iε = [r − ε, r + ε] ⊂ I y se tiene que |en| → 0 y<br />
siempre que x0 ∈ Iε.<br />
(4.2)<br />
|en+1| ≤ 1 M<br />
2 δ |en| 2 , (4.3)<br />
Demostración. Como las cotas para f ′ y f ′′ siguen siendo ciertas para cualquier subintervalo de<br />
I, podemos elegir ε > 0 tal que<br />
1<br />
2<br />
M<br />
ε = λ < 1.<br />
δ<br />
Entonces, si x0 ∈ Iε tenemos que |e0| = |x0 − r| < ε y usando (4.2) obtenemos<br />
En particular, x1 ∈ Iε. Análogamente,<br />
|e1| = |x1 − r| ≤ λ|e0|.<br />
|e2| = |x2 − r| ≤ λ|e1| ≤ λ 2 |e0|