function y=HH(t,x)

v=x(1); % potencial de membrana
n=x(2); %Subunidad n
m=x(3); %Subunidad m
h=x(4); %Subunidad h

% Parametros
gNa=120;
gK=36;
gL=0.3;
VNa=115;
VK=-12;
VL=10.6;
Cm=1;
An=0.01*((10-v)/(exp((10-v)*0.1)-1));
Bn=0.125*exp(-v/80);
Am=0.1*((25-v)/(exp((25-v)*0.1)-1));
Bm=4*exp(-v/18);
Ah=0.07*exp(-v/20);
Bh=1/(1+exp((30-v)*0.1));

% inyeccion de corriente en un intervalo
if t>2 & t<4
    I=-5;
else
    I=0;
end

%ecuaiones de HH
y(1)=-( gNa*m^3*h*(v-VNa)+ gK*n^4*(v-VK)+gL*(v-VL)+ I ) / Cm;  %Potencial de membrana%
y(2)=An*(1-n)-Bn*n;  %Subunidad n%
y(3)=Am*(1-m)-Bm*m;  %Subunidad m%
y(4)=Ah*(1-h)-Bh*h;  %Subunidad h%

y = y';