//Created by libAntimony v2.4
model shorten_2000__environment(time_, Ca_e, K_e, K_i, V_tau, k_tau)

  // Variable initializations:
  time_ = ;
  Ca_e = 20.0;
  K_e = 5.6;
  K_i = 140.0;
  V_tau = -60.0;
  k_tau = 22.0;
end

model shorten_2000__membrane(V, R, T, F, time_, i_Ca_L, i_Ca_T, i_K_DR, i_K_Ca, i_K_IR)

  // Rate Rules:
  V' = -(i_Ca_L + i_Ca_T + i_K_DR + i_K_Ca + i_K_IR) / Cm;

  // Variable initializations:
  V = ;
  R = 8.314;
  T = 310.0;
  F = 96845.0;
  Cm = 7.0;
  time_ = ;
  i_Ca_L = ;
  i_Ca_T = ;
  i_K_DR = ;
  i_K_Ca = ;
  i_K_IR = ;
end

model shorten_2000__L_type_calcium_current_m_gate(m_L, V, V_tau, k_tau, time_)

  // Assignment Rules:
  m_L_infinity := 1.0 / (1.0 + exp((V_m_L - V) / k_m_L));
  tau_m_L := tau_m_L_max / (exp((V - V_tau) / k_tau) + 2.0 * exp((2.0 * (V_tau - V)) / k_tau));

  // Rate Rules:
  m_L' = (m_L_infinity - m_L) / tau_m_L;

  // Variable initializations:
  m_L = ;
  tau_m_L_max = 27.0;
  V_m_L = -12.0;
  k_m_L = 12.0;
  V = ;
  V_tau = ;
  k_tau = ;
  time_ = ;
end

model shorten_2000__L_type_calcium_current(i_Ca_L, V_Ca, time_, V, V_tau, k_tau, R, T, F, Ca_e, Ca_i, m_L)

  // Sub-modules, and any changes to those submodules:
  L_type_calcium_current_m_gate: shorten_2000__L_type_calcium_current_m_gate(m_L, V, V_tau, k_tau, time_);

  // Assignment Rules:
  i_Ca_L := g_Ca_L * power(m_L, 2.0) * V_Ca;
  V_Ca := V * (((Ca_i - Ca_e) * exp(-(2.0) * ((F * V) / (R * T)))) / (1.0 - exp(-(2.0) * ((F * V) / (R * T)))));

  // Variable initializations:
  g_Ca_L = 9.0;
  R = ;
  T = ;
  F = ;
  Ca_e = ;
  Ca_i = ;
end

model shorten_2000__T_type_calcium_current_m_gate(m_T, V, V_tau, k_tau, time_)

  // Assignment Rules:
  m_T_infinity := 1.0 / (1.0 + exp((V_m_T - V) / k_m_T));
  tau_m_T := tau_m_T_max / (exp((V - V_tau) / k_tau) + 2.0 * exp((2.0 * (V_tau - V)) / k_tau));

  // Rate Rules:
  m_T' = (m_T_infinity - m_T) / tau_m_T;

  // Variable initializations:
  m_T = ;
  tau_m_T_max = 10.0;
  V_m_T = -30.0;
  k_m_T = 10.5;
  V = ;
  V_tau = ;
  k_tau = ;
  time_ = ;
end

model shorten_2000__T_type_calcium_current_h_gate(h_T, V, time_)

  // Assignment Rules:
  h_T_infinity := 1.0 / (1.0 + exp((V - V_h_T) / k_h_T));

  // Rate Rules:
  h_T' = (h_T_infinity - h_T) / tau_h_T;

  // Variable initializations:
  h_T = ;
  tau_h_T = 15.0;
  V_h_T = -57.0;
  k_h_T = 5.0;
  V = ;
  time_ = ;
end

model shorten_2000__T_type_calcium_current(i_Ca_T, time_, V, V_tau, k_tau, V_Ca, m_T, h_T)

  // Sub-modules, and any changes to those submodules:
  T_type_calcium_current_m_gate: shorten_2000__T_type_calcium_current_m_gate(m_T, V, V_tau, k_tau, time_);
  T_type_calcium_current_h_gate: shorten_2000__T_type_calcium_current_h_gate(h_T, V, time_);

  // Assignment Rules:
  i_Ca_T := g_Ca_T * power(m_T, 2.0) * h_T * V_Ca;

  // Variable initializations:
  g_Ca_T = 10.0;
  V_Ca = ;
end

model shorten_2000__voltage_sensitive_K_current_n_gate(n, V, time_)

  // Assignment Rules:
  n_infinity := 1.0 / (1.0 + exp((V_n - V) / k_n));

  // Rate Rules:
  n' = (n_infinity - n) / tau_n;

  // Variable initializations:
  n = ;
  tau_n = 20.0;
  V_n = -20.0;
  k_n = 4.5;
  V = ;
  time_ = ;
end

model shorten_2000__voltage_sensitive_K_current(i_K_DR, V_K, time_, V, R, T, F, K_e, K_i, n)

  // Sub-modules, and any changes to those submodules:
  voltage_sensitive_K_current_n_gate: shorten_2000__voltage_sensitive_K_current_n_gate(n, V, time_);

  // Assignment Rules:
  i_K_DR := g_K_DR * n * (V - V_K);
  V_K := V * (((K_i - K_e) * exp(-(1.0) * ((F * V) / (R * T)))) / (1.0 - exp(-(1.0) * ((F * V) / (R * T)))));

  // Variable initializations:
  g_K_DR = 0.1;
  R = ;
  T = ;
  F = ;
  K_e = ;
  K_i = ;
end

model shorten_2000__Ca_activated_K_current(i_K_Ca, time_, V, V_K, Ca_i)

  // Assignment Rules:
  i_K_Ca := g_K_Ca * (power(Ca_i, 4.0) / (power(Ca_i, 4.0) + power(Kc, 4.0))) * V_K;

  // Variable initializations:
  g_K_Ca = 0.09;
  Kc = 0.5;
  time_ = ;
  V = ;
  V_K = ;
  Ca_i = ;
end

model shorten_2000__inward_rectifier_current(i_K_IR, time_, V)

  // Assignment Rules:
  i_K_IR := piecewise( I_max * (1.0 - exp((V_K_IR - V) / ks)) , V > V_K_IR , g_K_IR * (V - V_K_IR) );

  // Variable initializations:
  I_max = 2.7;
  g_K_IR = 0.3;
  V_K_IR = -71.5;
  ks = 6.5;
  time_ = ;
  V = ;
end

model shorten_2000__cytosolic_calcium(Ca_i, i_Ca_L, i_Ca_T, time_)

  // Assignment Rules:
  J_in := - alpha * (i_Ca_L + i_Ca_T);
  J_eff := (v_p * power(Ca_i, 2.0)) / (power(Ca_i, 2.0) + power(K_p, 2.0));

  // Rate Rules:
  Ca_i' = f_cyt * beta * (J_in - J_eff);

  // Variable initializations:
  Ca_i = ;
  K_p = 0.08;
  f_cyt = 0.005;
  v_p = 0.025;
  alpha = 0.0074;
  beta = 0.4;
  i_Ca_L = ;
  i_Ca_T = ;
  time_ = ;
end

model *shorten_2000____main()

  // Sub-modules, and any changes to those submodules:
  environment: shorten_2000__environment(time_, Ca_e, K_e, K_i, V_tau, k_tau);
  membrane: shorten_2000__membrane(V, R, T, F, time_, i_Ca_L, i_Ca_T, i_K_DR, i_K_Ca, i_K_IR);
  L_type_calcium_current: shorten_2000__L_type_calcium_current(i_Ca_L, V_Ca, time_, V, V_tau, k_tau, R, T, F, Ca_e, Ca_i, m_L);
  T_type_calcium_current: shorten_2000__T_type_calcium_current(i_Ca_T, time_, V, V_tau, k_tau, V_Ca, m_T, h_T);
  voltage_sensitive_K_current: shorten_2000__voltage_sensitive_K_current(i_K_DR, V_K, time_, V, R, T, F, K_e, K_i, n);
  Ca_activated_K_current: shorten_2000__Ca_activated_K_current(i_K_Ca, time_, V, V_K, Ca_i);
  inward_rectifier_current: shorten_2000__inward_rectifier_current(i_K_IR, time_, V);
  cytosolic_calcium: shorten_2000__cytosolic_calcium(Ca_i, i_Ca_L, i_Ca_T, time_);
end