11.08.2013 Views

Apunte

Apunte

Apunte

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

80 4. RESOLUCIÓN DE ECUACIONES NO LINEALES<br />

es decir,<br />

0 = f(xn) + (xn+1 − xn) f(xn) − f(xn−1)<br />

xn − xn−1<br />

xn − xn−1<br />

xn+1 = xn − f(xn)<br />

f(xn) − f(xn−1) .<br />

Observemos que esta fórmula es análoga a la de Newton reemplazando f ′ por un cociente<br />

incremental.<br />

La ventaja es que no hay que calcular la derivada de f (esto es de gran ayuda en un caso en que<br />

f ′ sea difícil de calcular).<br />

La desventaja es, según veremos, que la convergencia de la sucesión es más lenta que la que<br />

produce el método de Newton.<br />

Observemos que la iteración del método de la secante también puede escribirse como<br />

xn+1 = f(xn)xn−1 − f(xn−1)xn<br />

.<br />

f(xn) − f(xn−1)<br />

Analicemos el orden de convergencia de este método, según la definición 4.8. Tenemos<br />

Luego,<br />

Es decir,<br />

f(r) = 0 y en = r − xn.<br />

en+1 = r − xn+1 = r − f(xn)xn−1 − f(xn−1)xn<br />

f(xn) − f(xn−1)<br />

= en−1f(xn) − enf(xn−1)<br />

f(xn) − f(xn−1)<br />

= en−1(f(xn) − f(r)) − en(f(r) − f(xn−1))<br />

f(xn) − f(xn−1)<br />

= −enen−1 f(xn)−f(r)<br />

en+1 = enen−1<br />

xn−r<br />

f(xn−1)−f(r)<br />

xn−1−r<br />

f(r)−f(xn−1)<br />

+ enen−1<br />

f(xn) − f(xn−1)<br />

f(r)−f(xn)<br />

−<br />

xn−r<br />

f(xn) − f(xn−1)<br />

xn−1−r<br />

.<br />

. (4.7)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!