You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
3. MÉTODOS MULTIPASO LINEALES 189<br />
Vamos a ver que para los métodos multipaso no es suficiente pedir que τ tienda a cero para que<br />
el método sea convergente. En este caso serán dos las condiciones (consistencia y estabilidad)<br />
que aseguren la convergencia de xn+k a la solución exacta.<br />
Definición 8.18. Se dice que un método multipaso es consistente si el error de truncamiento<br />
local, τ, tiende a cero.<br />
En forma análoga a la Definición 8.7 de orden para los métodos de un paso, se tiene que un<br />
método multipaso es de orden p si τ = O(h p ).<br />
En lo que sigue, estudiaremos un algoritmo que nos permite establecer el orden de un método<br />
multipaso. Para esto, necesitamos asumir que x, la solución de la ecuación diferencial x ′ = f(t, x)<br />
admite p + 1 derivadas continuas. En tal caso, de la expresión (8.24) se puede escribir hτ en<br />
términos de las derivadas de x como sigue<br />
hτ = C0x(t) + C1hx ′ (t) + C2h 2 x ′′ (t) + · · · + Cph p x (p) (t) + O(h p+1 ). (8.25)<br />
En efecto, consideremos el desarrollo de Taylor para x y x ′ ,<br />
<br />
x(t + jh) = x(t) + x ′ (t)jh + x′′ (t)<br />
2 (jh) 2 + · · ·<br />
x ′ (t + jh) = x ′ (t) + x ′′ (t)jh + x′′′ (t)<br />
2 (jh) 2 + · · ·<br />
Ahora, reemplazamos en la expresión (8.24) y agrupamos las potencias de h. De esta forma, se<br />
obtiene (8.25) con<br />
C0 = α0 + · · · + αk,<br />
C1 = α1 + 2α2 + 3α3 + · · · + kαk − β0 − β1 − β2 − · · · − βk,<br />
y en general para cualquier q ≥ 1,<br />
Cq = 1<br />
q! (α1 + 2 q α2 + 3 q α3 + · · · + k q 1<br />
αk) −<br />
(q − 1)! (β1 + 2 q−1 β2 + · · · + k q−1 βk). (8.26)<br />
De la expresión (8.26) se obtiene inmediatamente<br />
Proposición 8.19. Para un método de multipaso se tiene<br />
i) El método es consitente si y sólo si C0 = C1 = 0.<br />
i) El método es de orden p si y sólo si los valores Cq, q ≥ 1 definidos por (8.26) verifican<br />
C0 = C1 = · · · = Cp = 0 y Cp+1 = 0.