// Created by libAntimony v2.8.0
function GAMMAF(VV, theta, sigma)
  1/(1 + exp(-(VV - theta)/sigma));
end


model *model_0000001()

  // Compartments and Species:
  compartment compartment_0000001;

  // Assignment Rules:
  INa := gNa*Minfs^pms*hhs*(V - VNa);
  Minfs := GAMMAF(V, thetam, sigmam);
  INaP := gNaP*Pinfs*(V - VNa);
  Pinfs := GAMMAF(V, thetap, sigmap);
  IKdr := gKdr*nns^pns*(V - VK);
  IA := gA*Ainfs^3*bbs*(V - VK);
  Ainfs := GAMMAF(V, thetaa, sigmaa);
  Iz := gZ*zzs*(V - VK);

  // Rate Rules:
  hhs' = (phi*(GAMMAF(V, thetah, sigmah) - hhs))/(1 + 7.5*GAMMAF(V, t_tauh, -6));
  V' = (-gL*(V - VL) - INa - INaP - IKdr - IA - Iz + Iapp)/Cm;
  nns' = (phi*(GAMMAF(V, thetan, sigman) - nns))/(1 + 5*GAMMAF(V, t_taun, -15));
  bbs' = (GAMMAF(V, thetab, sigmab) - bbs)/tauBs;
  zzs' = (GAMMAF(V, thetaz, sigmaz) - zzs)/tauZs;

  // Compartment initializations:
  compartment_0000001 = 1;

  // Variable initializations:
  Cm = 1;
  pms = 3;
  pns = 4;
  VNa = 55;
  t_tauh = -40.5;
  t_taun = -27;
  thetaa = -50;
  sigmaa = 20;
  thetab = -80;
  sigmab = -6;
  tauBs = 15;
  sigmam = 9.5;
  sigmah = -7;
  sigman = 10;
  sigmaz = 5;
  gNa = 35;
  gKdr = 6;
  gL = 0.05;
  Iapp = 0.661914;
  gA = 1.4;
  gNaP = 0.3;
  gZ = 1;
  thetaz = -39;
  tauZs = 75;
  phi = 10;
  thetah = -45;
  thetam = -30;
  thetan = -35;
  thetap = -47;
  sigmap = 3;
  VK = -90;
  VL = -70;
  hhs = 0.98786;
  V = -71.81327;
  nns = 0.02457;
  bbs = 0.203517;
  zzs = 0.00141;

  // Other declarations:
  var INa, Minfs, hhs, V, INaP, Pinfs, IKdr, nns, IA, Ainfs, bbs, Iz, zzs;
  const compartment_0000001, Cm, pms, pns, VNa, t_tauh, t_taun, thetaa, sigmaa;
  const thetab, sigmab, tauBs, sigmam, sigmah, sigman, sigmaz, gNa, gKdr;
  const gL, Iapp, gA, gNaP, gZ, thetaz, tauZs, phi, thetah, thetam, thetan;
  const thetap, sigmap, VK, VL;

  // Unit definitions:
  unit time_unit = 1e-3 second;
  unit mV = 1e-3 volt;

  // Display Names:
  time_unit is "ms";
  V is "VVs";
end