// Created by libAntimony v2.13.0
model *Plant1981_BurstingNerveCells()

  // Compartments and Species:
  compartment COMpartment;
  species $V_membrane in COMpartment, $h1 in COMpartment, $x1 in COMpartment;
  species $n1 in COMpartment, $c in COMpartment;

  // Assignment Rules:
  i_Na := g_I*m_infinity^3*h1*(V_I - V_membrane);
  i_Ca := g_T*x1*(V_I - V_membrane);
  i_K := g_K*n1^4*(V_K - V_membrane);
  i_K_Ca := (g_K_Ca*c/(K_p + c))*(V_K - V_membrane);
  i_L := g_L*(V_L - V_membrane);
  h_infinity := alpha_h/(alpha_h + beta_h);
  tau_h := 12.5/(alpha_h + beta_h);
  x_infinity := 1/(exp(0.15*(-V_membrane - 50)) + 1);
  n_infinity := alpha_n/(alpha_n + beta_n);
  tau_n := 12.5/(alpha_n + beta_n);
  m_infinity := alpha_m/(alpha_m + beta_m);
  Vs := a*V_membrane + b;
  alpha_m := 0.1*(50 - Vs)/(exp((50 - Vs)/10) - 1);
  beta_m := 4*exp((25 - Vs)/18);
  alpha_h := 0.07*exp((25 - Vs)/20);
  beta_h := 1/(exp((55 - Vs)/10) + 1);
  alpha_n := 0.01*(55 - Vs)/(exp((55 - Vs)/10) - 1);
  beta_n := 0.125*exp((45 - Vs)/80);

  // Rate Rules:
  V_membrane' = i_Na + i_Ca + i_K + i_K_Ca + i_L;
  h1' = (h_infinity - h1)/tau_h;
  x1' = (x_infinity - x1)/tau_x;
  n1' = (n_infinity - n1)/tau_n;
  c' = f*(K_c*x1*(V_Ca - V_membrane) - c);

  // Species initializations:
  V_membrane = -55;
  h1 = 0.9;
  x1 = 0.27;
  n1 = 0.03;
  c = 0.4;

  // Compartment initializations:
  COMpartment = 1;

  // Variable initializations:
  tau_x = 235;
  f = 0.0003;
  K_c = 0.0085;
  V_Ca = 140;
  g_I = 4;
  V_I = 30;
  V_K = -75;
  V_L = -40;
  g_K = 0.3;
  g_T = 0.01;
  g_K_Ca = 0.03;
  g_L = 0.003;
  K_p = 0.5;
  a = 1.209;
  b = 78.714;
  g_Ca = 0.004;

  // Other declarations:
  var i_Na, i_Ca, i_K, i_K_Ca, i_L, h_infinity, tau_h, x_infinity, n_infinity;
  var tau_n, m_infinity, Vs, alpha_m, beta_m, alpha_h, beta_h, g_Ca, alpha_n;
  var beta_n;
  const COMpartment, tau_x, f, K_c, V_Ca, g_I, V_I, V_K, V_L, g_K, g_T, g_K_Ca;
  const g_L, K_p, a, b;

  // Display Names:
  V_membrane is "V";

  // SBO terms:
  V_membrane.sboTerm = 259
  h1.sboTerm = 247
  x1.sboTerm = 247
  n1.sboTerm = 247
  c.sboTerm = 247

  // CV terms:
  V_membrane hypernym "http://identifiers.org/pato/PATO:0001462"
  h1 hypernym "http://identifiers.org/chebi/CHEBI:29101"
  x1 hypernym "http://identifiers.org/chebi/CHEBI:29108"
  n1 hypernym "http://identifiers.org/chebi/CHEBI:29103"
  c hypernym "http://identifiers.org/chebi/CHEBI:29108"
end

Plant1981_BurstingNerveCells is "Plant1981_BurstingNerveCells"

Plant1981_BurstingNerveCells model_entity_is "http://identifiers.org/biomodels.db/MODEL6762427183"
Plant1981_BurstingNerveCells model_entity_is "http://identifiers.org/biomodels.db/BIOMD0000000304"
Plant1981_BurstingNerveCells description "http://identifiers.org/pubmed/7252375"
Plant1981_BurstingNerveCells taxon "http://identifiers.org/taxonomy/6499",
                                   "http://identifiers.org/bto/BTO:0000022",
                                   "http://identifiers.org/bto/BTO:0000938"
Plant1981_BurstingNerveCells hypernym "http://identifiers.org/go/GO:0006816",
                                      "http://identifiers.org/go/GO:0006814",
                                      "http://identifiers.org/go/GO:0006813",
                                      "http://identifiers.org/go/GO:0019228"
Plant1981_BurstingNerveCells property "http://identifiers.org/mamo/MAMO_0000046"
