texto: metodos numericos para ecuaciones diferenciales ordinarias
texto: metodos numericos para ecuaciones diferenciales ordinarias
texto: metodos numericos para ecuaciones diferenciales ordinarias
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
c = h<br />
2 f(x,y)<br />
En consecuencia<br />
donde<br />
R1(x + h<br />
2<br />
T (2) <br />
(x,y) = f x + h h<br />
,y +<br />
2 2 f(x,y)<br />
<br />
h h2<br />
,y + f(x,y)) =<br />
2 8<br />
− R1<br />
<br />
x + h h<br />
,y +<br />
2 2 f(x,y)<br />
<br />
∂2f h2<br />
∂x2(ǫ,µ) +<br />
4 f(x,y) ∂2f h2<br />
(ǫ,µ) +<br />
∂x∂y<br />
8 (f(x,y))2 ∂2f (ǫ,µ)<br />
∂y2 El método que resulta de sustituir T (2) (x,y) por f x + h h ,y + 2 2f(x,y) en el método de<br />
Taylor de orden dos es un método específico de Runge-Kutta, conocido con el nombre de<br />
Método del Punto medio<br />
3.1. Método del punto medio<br />
yn+1 = yn + hf(xn + h<br />
2 ,yn + h<br />
2 f(xn,yn)) (3.3)<br />
n = 0, 1, 2,... (3.4)<br />
Solo tres parámetros se encuentran en af(x+b,y+c) y los tres se requieren en la igualdad<br />
con T (2) ,necesitamos una forma mas compleja <strong>para</strong> cumplir las condiciones que requiere<br />
cualquiera de los métodos de Taylor de orden superior.<br />
La forma más apropiada de cuatro parámetros con que se aproxima<br />
es<br />
T (3) (x,y) = f(x,y) + h<br />
2 f ′ (x,y) + h2<br />
6 f ′′ (x,y)<br />
af(x,y) + bf(x + c,y + df(x,y)) (3.5)<br />
y ni siquiera con esto se tiene la suficiente flexibilidad <strong>para</strong> igualar el término<br />
h 2<br />
6<br />
resultante de la expansión de h2<br />
6 f ′′ (x,y)<br />
∂f<br />
∂y (x,y)<br />
2<br />
f(x,y)<br />
Por lo tanto, lo mejor que podemos lograr utilizando son métodos con el error local<br />
de truncamiento O(h 2 ), sin embargo el hecho que tenga cuatro parámetros, da cierta<br />
flexibilidad en su elección <strong>para</strong> derivar varios métodos O(h 2 ). Uno de los más importantes<br />
es el método modificado de Euler, que corresponde a seleccionar a = b = 1<br />
2<br />
y se representa de la siguiente forma.<br />
22<br />
y c = d = h