function p=memoscparam(ver)

if ver==1
  % Corinto, Ascoli, BCModel, p. 2725, G=3.3 mS, GN1=-0.4 mS, GN2=-1.2 mS, Goff=0.06 mS, Gon=1.9 mS
  p.G=3.3e-3; p.GN1=-0.4e-3; p.GN2=-1.2e-3;    
  p.Goff=0.06e-3; p.Gon=1.9e-3;
  p.alpha=0.74;   p.beta=0.0333;  p.gamma=0.12;
  p.i0=8.9189; p.eta1=1; p.eta2=-1; %p.i0=1e-2;     
  p.Region=4; % correct for initial conditions: [0.006,0.02,-0.3,0,0]' 
  %if the initial condition wstart is changed then p.Region must be set properly
end
