El trabajo para el próximo martes es resolver por Euler y programar RK de cuarto orden para el ejercicio 14 pág. 88 del capitúlo 4, de Paviz Moin sobre ODEs.
Grafiquen las concentraciones (Euler vs. RK)
Cualquier cosa, por el blogspot
Este modelo es muy simpático, deben hacer un hstep adecuado y verán que las soluciones se comportan "raras".
Si pueden busquen en la red "stiff Differential Equations" para tener una idea de este tipo de ecuaciones muy comunes en la vida práctica.
Suscribirse a:
Enviar comentarios (Atom)
2 comentarios:
Hola profesor aun no entiendo cual es la EDO que tenemos que usar en el método de Runge-Kutta de 4ºorden, tengo entendido que el método consiste en yi+1=yi+hstep*F(x,y)
donde
F(x,y)=(k1+2k2+2k3+k4)/6
donde
k1=f(xi,yi)
k2=f(x1+[htep/2],yi+[0.5*htep*k1])
k3=f(x1+[htep/2],yi+[0.5*htep*k2])
k2=f(x1+[htep],yi+[htep*k3])
despuúes se realiza un proceso iterativo hasta llegar al último intervalo.
Esperó su respuesta, gracias.
ajá así es
RK es un método explícito evaluado en puntos muy particulares, los que has escrito.
la f(xi,yi) que escribiste es el lado derecho de la ODE a resolver.
En el caso del ejercicio, este es un sistema de 4 ecuaciones, es decir, 4 lado derecho.
??? más o menos claro? lo que hicimos en clase el jueves fue para el modelo predador-presa... 2 ecuaciones, hubo 2 lado derechos a los que aplicamos Euler.
Ahora debes hacer lo propio, pero para RK y para 4 ecuaciones ;)
Publicar un comentario