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.

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.

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

Saved successfully!

Ooh no, something went wrong!