%
% the reduced mec model from Koppel and co (sec 3.2.1)
%
%
% C v' = -Gp*pinf(v)(v-ENa) - GL*(v-EL) - Gh*(.65rf+.35rs)(v-Eh)
%
%  rf' = (rfinf(v)-rf)/taurf(v)
%
%  rs' = (rsinf(v)-rs)/taurs(v)
%

v = -120:40;

pinf = 1./(1+exp(-(v+38)/6.5));

rfinf = 1./(1+exp((v+79.2)/9.78));

taurf = 1 + 0.51./(exp((v-1.7)/10)+exp(-(v+340)/52));

rsinf = 1./(1+exp((v+2.83)/15.9)).^58;

taurs = 1 + 5.6./(exp((v-1.7)/14)+exp(-(v+260)/43));
