11.08.2013 Views

Apunte

Apunte

Apunte

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

El primer método a analizar es<br />

xtn+3<br />

3. MÉTODOS MULTIPASO LINEALES 193<br />

− xtn+1 = h[7<br />

3 fn+2 − 2<br />

3 fn+1 + 1<br />

3 fn],<br />

Para estudiar la consistencia y la condición de la raíz escribimos los polinomios p y q.<br />

p(z) = z 3 − z y q(z) = 7<br />

3 z2 − 2 1<br />

z +<br />

3 3 .<br />

La consistencia se verifica puesto que p(1) = 0 y p ′ (1) = 3z 2 −1|z=1 = 2 = q(1). Para la condición<br />

de la raíz, vemos que las raíces de p son {−1, 0, 1}, todas distintas y de módulo menor o igual a<br />

1. Luego, el método es convergente.<br />

Para el método de Milne tenemos<br />

p(z) = z 2 − 1 y q(z) = 1<br />

3 z2 + 4 1<br />

z +<br />

3 3 .<br />

= q(1). Para la<br />

condición de la raíz, notar que las raíces de p son {−1, 1} sin multiplicidad y de módulo menor<br />

o igual a 1. Luego, el método es convergente.<br />

La consistencia se verifica puesto que p(1) = 0 y p ′ (1) = 2z|z=1 = 2 = 6<br />

3<br />

Para finalizar analizaremos el orden de convergencia de un método.<br />

Ejemplo 8.24. Se quiere determinar el orden de convergencia de<br />

(a) el método hallado en el Ejemplo 8.17.<br />

(b) el método de Milne dado por la ecuación (8.23).<br />

Las constantes Cq con q ≥ 1 dependen de αj y βj, j = 0, . . . , k.<br />

Para el ítem (a), k = 3, (α0, α1, α2, α3) = (0, −1, 0, 1) y (β0, β1, β2, β3) = ( 1<br />

3<br />

, − 2<br />

3<br />

Como α0 = α2 = 0 y β3 = 0 omitimos escribirlos. Calculamos los coeficientes Cq:<br />

C0 = α1 + α3 = 0,<br />

C1 = α1 + 3α3 − (β0 + β1 + β2) = −1 + 3 − ( 6<br />

3 ) = 0,<br />

C2 = 1<br />

2 (α1 + 32α3) − (β1 + 2β2) = 1<br />

12<br />

2 (−1 + 9) − 3 = 0.<br />

C3 = 1<br />

6 (α1 + 33α3) − 1<br />

2 (β1 + 22β2) = 1<br />

1 2 28<br />

6 (−1 + 27) − 2 (− 3 + 3 ) = 0,<br />

C4 = 1<br />

24 (α1 + 34α3) − 1<br />

6 (β1 + 23β2) = 1<br />

1 2 56<br />

24 (−1 + 81) − 6 (− 3 + 3<br />

Luego, el método tiene orden de convergencia p = 3.<br />

) = 10<br />

3<br />

7 , 3 , 0).<br />

− 3 = 0.

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

Saved successfully!

Ooh no, something went wrong!