function dydt = ode_fhn(t,y)

v=y(1);  
w=y(2);  

if t>1 & t<1.01
    I=0.05;
else
    I=0;
end


    e=0.001;     % Parametro
    f = v*(v-.1)*(1-v)-w;
    g= v - 0.5*w;
vp=f+I;     % Ecuacion para v prima
wp=g;       % Ecuacion para w prima

dydt = [vp/e; wp];