Ecuaciones de Fitzhugh - Nagumo

Considerando las ecuaciones de Fitzhugh-Nagumo:
dv/dt= -v*(v-a)*(v-1)-w
dw/dt=b*(v-c*w)
siendo a, b y c parámetros que cumplen:
0<a<1, b>1 y c>0.
La variable v representa el potencial de acción y la variable w modela los canales de sodio y potasio.
Ilustrar (utilizando Matlab) la existencia de condiciones iniciales (v,0), cuyas correspondientes órbitas se aproximan a equilibrios distintos. Se supone que a=1/4, b=3/2 y c=10.
Lo que me sucede, es que no entiendo como meter las condiciones iniciales que me pide en matlab. Sé representar las órbitas para los puntos de equilibrio pero sin condiciones iniciales y por mas que he probado, no he conseguido nada.

function ode45caballero8 global a b c fprintf('CONDICIONES: 0<a<1, b>1 , c>0\n') fprintf('\n EJEMPLO: a=0.25  b=1.5  c=10  I=[0 1]  y0=[-4 0]\n\n') a=input('Ingrese el valor de a: '); b=input('Ingrese el valor de b: '); c=input('Ingrese el valor de c: '); I=input('Ingrese el intervalo I='); y0=input('Ingrese ondición inicial y0='); [t,y]=ode45(@potencial_modelado1,[I(1) I(2)],[y0(1) y0(2)]); figure(gcf) plot(t,y(:,1), t,y(:,2)); end function dydt=potencial_modelado1(t,y) global a b c dydt=zeros(2,1); dydt(1)=-y(1)*(y(1)-a)*(y(1)-1)-y(2);  %dv/dt= -v*(v-a)*(v-1)-w dydt(2)=b*(y(1)-c*y(2));               %dw/dt=b*(v-c*w) end

EJECUCION EN MATLAB

>> ode45caballero8 CONDICIONES: 0<a<1, b>1 , c>0 EJEMPLO: a=0.25    b=1.5    c=10  I=[0 1]    y0=[-4 0] Ingrese el valor de a: 0.25 Ingrese el valor de b: 1.5 Ingrese el valor de c: 3 Ingrese el intervalo I=[0 1] Ingrese condicion inicial y0=[1 2]

José Jeremías Caballero
Servicios de Programación en MatLab
https://programacioncaballero.blogspot.com/
Lima - Perú

Para mayor informacion, contactarse a jjcc94@hotmail.com