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.