%
%
%       Biofisica 2005 - FCEyN - UBA
%       FHN
%       Jacobo Sitt
%


function fhn


    ci=[0 0]; % Condiciones iniciales


[t,y]=ode23('ode_fhn',[0 5],ci,odeset('maxstep',.01)); % Integrador

v=y(:,1); % variable v
w=y(:,2); % variable w


vv=-1:.01:1.5;
[w1,w2]=nulc_fhn(0,vv);
plot (vv,w1,vv,w2,v,w) % Grafico del diagrama de fases
xlabel('v')
ylabel('w')


figure
plot(t,v,t,w) % Grafico de la evolucion temporal

xlabel('t')