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.