%
% I_NaP+I_K+I_KM model of bursting
%


init v=-2.5  n=0.01 nm=0.01

par El=-4.0  ENa=3.0  EK=-4.5  EKm=-4.5
par tauV=0.152 NaV12=-20.0  Nak=15.0
par KV12=-25.0  Kk=5.0
par a=0.005  b=0.8  c=2.0 d=0.5
par I=5.0
par alpha=0.05  vm=-20  Km=5.0

minf(v)=1/(1+exp((NaV12-20*v)/Nak))
ninf(v)=1/(1+exp((KV12-20*v)/Kk))
nminf(v)=1/(1+exp((vm-20*v)/Km))

% nminf1(v)=1-nminf(v)

v'=10*(a*I+b*(El-v)+c*minf(v)*(ENa-v)+n*(EK-v)+d*nm*(EKm-v))
n'=(ninf(v)-n)/tauV
nm'=alpha*(nminf(v)-nm)

@ total=200,method=gear, bound=500, xp=v,yp=n, xlo=-4,xhi=0,yhi=1,ylo=-0.3
done