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.

1. MÉTODOS DE UN PASO 173<br />

xi+1 = xi + hf(ti, xi) + 1<br />

2 h2 [ft(ti, xi) + fx(ti, xi)f(ti, xi)]. (8.7)<br />

La fórmula de iteración (8.7) se conoce con el nombre de Método de Taylor de orden 2 y<br />

corresponde a la expresión (8.4) si consideramos Φ(t, x, h) = f(t, x) + h<br />

2 [ft(t, x) + fx(t, x)f(t, x)].<br />

El método de Taylor de orden 3 se consigue derivando x ′′ . En lo que sigue escribimos simplemente<br />

x en lugar de x(t). La fórmula para x ′′′ queda en términos de f y sus derivadas,<br />

x ′′′ (t) = ftt(t, x) + 2ftx(t, x)f(t, x) + fxx(t, x)f 2 (t, x)<br />

+ft(t, x)fx(t, x) + f 2 t (t, x)f(t, x).<br />

Continuando de esta manera se pueden calcular las derivadas de x hasta orden k en términos<br />

de f y sus derivadas. Esto nos da una manera de calcular métodos de aproximación de un paso<br />

(solo usamos ti, xi y h para calcular la siguiente aproximación xi+1).<br />

Para formalizar este procedimiento introducimos dos notaciones. Llamamos D j f a la derivada<br />

total de f respecto de t, es decir<br />

D j f(t, x) = d(j)<br />

f(t, x(t)).<br />

dt<br />

Por ejemplo, si j = 0, D 0 f(t, x) = f(t, x) y para j = 1, D 1 f(t, x) = ft(t, x) + fx(t, x)f(t, x).<br />

Ahora, ponemos<br />

Tk(t, x, h) =<br />

Así, de (8.6) obtenemos la aproximación<br />

k hj−1 j=1<br />

x(t + h) ∼ x(t) + hTk(t, x, h)<br />

de donde se tiene la siguiente fórmula de recurrencia,<br />

xi+1 = xi + hTk(ti, xi, h).<br />

j! Dj f(t, x). (8.8)<br />

Estos métodos se conocen como métodos de Taylor de orden k. Su mayor problema es que<br />

requieren encontrar y evaluar derivadas sucesivas de f(t, x).<br />

Ejemplo 8.4. Aplicar el método de Euler y Taylor de orden 2 para el problema<br />

x ′ (t) = t+2<br />

t+1 x,<br />

x(0) = 1.

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

Saved successfully!

Ooh no, something went wrong!