// Created by libAntimony v2.8.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";
end