// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real unit_compartment = 1 L;
  real V(time) mV;
  real V_neg(time) mV;
  real E(time) mV;
  real I = 0;
  real i_Na(time);
  real i_K(time);
  real i_L(time);
  real m(time);
  real h(time);
  real n(time);
  real E_R = -75 mV;
  real Cm = 1;
  real g_Na = 120;
  real g_K = 36;
  real g_L = .3;
  real E_Na = -190 mV;
  real E_K = -63 mV;
  real E_L = -85.613 mV;
  real V_Na(time) mV;
  real V_K(time) mV;
  real V_L(time) mV;
  real alpha_m(time);
  real beta_m(time);
  real alpha_h(time);
  real beta_h(time);
  real alpha_n(time);
  real beta_n(time);

  // equations
  when (time=time.min) V = 0;
  V:time = (I-(i_Na+i_K+i_L))/Cm;
  V_neg = (-1)*V;
  E = V+E_R;
  i_Na = g_Na*m^3*h*(V-V_Na);
  i_K = g_K*n^4*(V-V_K);
  i_L = g_L*(V-V_L);
  when (time=time.min) m = .052932;
  m:time = alpha_m*(1-m)-beta_m*m;
  when (time=time.min) h = .59612;
  h:time = alpha_h*(1-h)-beta_h*h;
  when (time=time.min) n = .31768;
  n:time = alpha_n*(1-n)-beta_n*n;
  V_Na = E_Na-E_R;
  V_K = E_K-E_R;
  V_L = E_L-E_R;
  alpha_m = .1*(V+25)/(exp((V+25)/10)-1);
  beta_m = 4*exp(V/18);
  alpha_h = .07*exp(V/20);
  beta_h = 1/(exp((V+30)/10)+1);
  alpha_n = .01*(V+10)/(exp((V+10)/10)-1);
  beta_n = .125*exp(V/80);

  // variable properties
  unit_compartment.sbmlRole="compartment";
  V.sbmlRole="parameter";
  V_neg.sbmlRole="parameter";
  E.sbmlRole="parameter";
  I.sbmlRole="parameter";
  i_Na.sbmlRole="parameter";
  i_K.sbmlRole="parameter";
  i_L.sbmlRole="parameter";
  m.sbmlRole="parameter";
  h.sbmlRole="parameter";
  n.sbmlRole="parameter";
  E_R.sbmlRole="parameter";
  Cm.sbmlRole="parameter";
  g_Na.sbmlRole="parameter";
  g_K.sbmlRole="parameter";
  g_L.sbmlRole="parameter";
  E_Na.sbmlRole="parameter";
  E_K.sbmlRole="parameter";
  E_L.sbmlRole="parameter";
  V_Na.sbmlRole="parameter";
  V_K.sbmlRole="parameter";
  V_L.sbmlRole="parameter";
  alpha_m.sbmlRole="parameter";
  beta_m.sbmlRole="parameter";
  alpha_h.sbmlRole="parameter";
  alpha_h.sbmlName="auxiliary alpha_h";
  beta_h.sbmlRole="parameter";
  alpha_n.sbmlRole="parameter";
  beta_n.sbmlRole="parameter";
}