Tweet
Hola a todos
Estoy tratando de dibujar en matlab el campo de direcciones (direction field) de esta ecuación diferencial:
Para ello he seguido estas instrucciones (pág 3): http://www.slideshare.net/josholucho...curvassolucion
Las órdenes que ejecuto en la consola del programa son estas:
Que son las mismas que indican en ese tutorial, pero que valen para la ecuación diferencial que ponen de ejemplo pero no para la mía xD
Como podéis comprobar, cambié la variable t por x porque daba un error con la función (al final era que en la multiplicación y división hay que añadir un '.' tal que así '.*', './')
Lo he probado varias veces, he revisado los comandos, he cambiado el rango... la verdad es que no se donde puede estar el problema
A ver si me podéis echar una manita. ¡Gracias!
P.D.: lo estoy dibujando en matlab porque en mi facultad todos los profesores están que no cagan con él y es lo que nos enseñan en métodos numéricos, pero estoy abierto a otras posibilidades. Por ejemplo hasta ahora usaba gnuplot para dibujar funciones y hacer ajustes lineales, cuadráticos y tal... y tan contento!
Estoy tratando de dibujar en matlab el campo de direcciones (direction field) de esta ecuación diferencial:
Para ello he seguido estas instrucciones (pág 3): http://www.slideshare.net/josholucho...curvassolucion
Las órdenes que ejecuto en la consola del programa son estas:
Código:
f = inline ('2.*exp(-x) - ((x+1)./x).*y','y','x')
[y,x]=meshgrid (0:0.5:7,-3:0.5:3)
[n,m]=size(y)
dy=ones(n,m)
z=f(y,x)
dx=z
hold on,quiver(y,x,dy,dx)
Como podéis comprobar, cambié la variable t por x porque daba un error con la función (al final era que en la multiplicación y división hay que añadir un '.' tal que así '.*', './')
Lo he probado varias veces, he revisado los comandos, he cambiado el rango... la verdad es que no se donde puede estar el problema
A ver si me podéis echar una manita. ¡Gracias!
P.D.: lo estoy dibujando en matlab porque en mi facultad todos los profesores están que no cagan con él y es lo que nos enseñan en métodos numéricos, pero estoy abierto a otras posibilidades. Por ejemplo hasta ahora usaba gnuplot para dibujar funciones y hacer ajustes lineales, cuadráticos y tal... y tan contento!