//Created by libAntimony v2.4

// Warnings from automatic translation:
//    Unable to use the formula "K_mCSQN - (CSQN_max + CSQN_buff + delta_Ca_JSR + Ca_JSR_old)" to set the assignment rule for b1:  Loop detected:  b1's definition (K_mCSQN - (CSQN_max + CSQN_buff + delta_Ca_JSR + Ca_JSR_old)) either includes itself directly (i.e. 's5 := 6 + s5') or by proxy (i.e. 's5 := 8*d3' and 'd3 := 9*s5').
//    Unable to use the formula "K_mCSQN * (CSQN_buff + delta_Ca_JSR + Ca_JSR_old)" to set the assignment rule for c1:  Loop detected:  c1's definition (K_mCSQN * (CSQN_buff + delta_Ca_JSR + Ca_JSR_old)) either includes itself directly (i.e. 's5 := 6 + s5') or by proxy (i.e. 's5 := 8*d3' and 'd3 := 9*s5').
//    Unable to use the formula "TRPN_buff + CMDN_buff + delta_Cai + Cai_old" to set the assignment rule for Ca_total:  Loop detected:  Ca_total's definition (TRPN_buff + CMDN_buff + delta_Cai + Cai_old) either includes itself directly (i.e. 's5 := 6 + s5') or by proxy (i.e. 's5 := 8*d3' and 'd3 := 9*s5').

model zeng_1995__environment(time_)

  // Variable initializations:
  time_ = ;
end

model zeng_1995__membrane(V, R, T, F, time_, i_Na, i_Ca_L, i_Ca_T, i_Kr, i_Ks, i_NaCa, i_K1, i_Kp, i_p_Ca, i_Na_b, i_Ca_b, i_NaK, i_ns_Ca)

  // Rate Rules:
  V' = (-(1) / Cm) * (i_Na + i_Ca_L + i_Ca_T + i_Kr + i_Ks + i_K1 + i_Kp + i_NaCa + i_p_Ca + i_Na_b + i_Ca_b + i_NaK + i_ns_Ca + I_st);

  // Variable initializations:
  V = -84.624;
  R = 8314;
  T = 310;
  F = 96500;
  Cm = 1;
  I_st = -100;
  time_ = ;
  i_Na = ;
  i_Ca_L = ;
  i_Ca_T = ;
  i_Kr = ;
  i_Ks = ;
  i_NaCa = ;
  i_K1 = ;
  i_Kp = ;
  i_p_Ca = ;
  i_Na_b = ;
  i_Ca_b = ;
  i_NaK = ;
  i_ns_Ca = ;
end

model zeng_1995__fast_sodium_current_m_gate(m, V, time_)

  // Assignment Rules:
  alpha_m := (0.32 * (V + 47.13)) / (1 - exp(-(0.1) * (V + 47.13)));
  beta_m := 0.08 * exp(- V / 11);

  // Rate Rules:
  m' = alpha_m * (1 - m) - beta_m * m;

  // Variable initializations:
  m = ;
  V = ;
  time_ = ;
end

model zeng_1995__fast_sodium_current_h_gate(h, V, time_)

  // Assignment Rules:
  alpha_h := piecewise( 0.135 * exp((80 + V) / -(6.8)) , V < -(40) , 0 );
  beta_h := piecewise( 3.56 * exp(0.079 * V) + 310000 * exp(0.35 * V) , V < -(40) , 1 / (0.13 * (1 + exp((V + 10.66) / -(11.1)))) );

  // Rate Rules:
  h' = alpha_h * (1 - h) - beta_h * h;

  // Variable initializations:
  h = ;
  V = ;
  time_ = ;
end

model zeng_1995__fast_sodium_current_j_gate(j, V, time_)

  // Assignment Rules:
  alpha_j := piecewise( ((-(127140) * exp(0.2444 * V) - 0.00003474 * exp(-(0.04391) * V)) * (V + 37.78)) / (1 + exp(0.311 * (V + 79.23))) , V < -(40) , 0 );
  beta_j := piecewise( (0.1212 * exp(-(0.01052) * V)) / (1 + exp(-(0.1378) * (V + 40.14))) , V < -(40) , (0.3 * exp(-(0.0000002535) * V)) / (1 + exp(-(0.1) * (V + 32))) );

  // Rate Rules:
  j' = alpha_j * (1 - j) - beta_j * j;

  // Variable initializations:
  j = ;
  V = ;
  time_ = ;
end

model zeng_1995__fast_sodium_current(i_Na, E_Na, time_, V, R, F, T, Nai, Nao, m, h, j)

  // Sub-modules, and any changes to those submodules:
  fast_sodium_current_m_gate: zeng_1995__fast_sodium_current_m_gate(m, V, time_);
  fast_sodium_current_h_gate: zeng_1995__fast_sodium_current_h_gate(h, V, time_);
  fast_sodium_current_j_gate: zeng_1995__fast_sodium_current_j_gate(j, V, time_);

  // Assignment Rules:
  i_Na := g_Na * power(m, 3) * h * j * (V - E_Na);
  E_Na := ((R * T) / F) * ln(Nao / Nai);

  // Variable initializations:
  g_Na = 16;
  R = ;
  F = ;
  T = ;
  Nai = ;
  Nao = ;
end

model zeng_1995__L_type_Ca_channel_d_gate(d, V, time_)

  // Assignment Rules:
  alpha_d := d_infinity / tau_d;
  beta_d := (1 - d_infinity) / tau_d;
  d_infinity := 1 / (1 + exp(-(V + 10) / 6.24));
  tau_d := (d_infinity * (1 - exp(-(V + 10) / 6.24))) / (0.035 * (V + 10));

  // Rate Rules:
  d' = alpha_d * (1 - d) - beta_d * d;

  // Variable initializations:
  d = ;
  V = ;
  time_ = ;
end

model zeng_1995__L_type_Ca_channel_f_gate(f, V, time_)

  // Assignment Rules:
  alpha_f := f_infinity / tau_f;
  beta_f := (1 - f_infinity) / tau_f;
  f_infinity := 1 / (1 + exp((V + 35.06) / 8.6)) + 0.6 / (1 + exp((50 - V) / 20));
  tau_f := 1 / (0.0197 * exp(-power(0.0337 * (V + 10), 2)) + 0.02);

  // Rate Rules:
  f' = alpha_f * (1 - f) - beta_f * f;

  // Variable initializations:
  f = ;
  V = ;
  time_ = ;
end

model zeng_1995__L_type_Ca_channel_f_Ca_gate(f_Ca, V, time_, Cai)

  // Assignment Rules:
  f_Ca := 1 / (1 + power(Cai / Km_Ca, 2));

  // Variable initializations:
  Km_Ca = 0.6;
  V = ;
  time_ = ;
  Cai = ;
end

model zeng_1995__L_type_Ca_channel(i_Ca_L, i_CaCa, i_CaK, i_CaNa, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, time_, V, Cai, R, T, F, Cao, Nao, Ko, Nai, Ki, d, f, f_Ca)

  // Sub-modules, and any changes to those submodules:
  L_type_Ca_channel_d_gate: zeng_1995__L_type_Ca_channel_d_gate(d, V, time_);
  L_type_Ca_channel_f_gate: zeng_1995__L_type_Ca_channel_f_gate(f, V, time_);
  L_type_Ca_channel_f_Ca_gate: zeng_1995__L_type_Ca_channel_f_Ca_gate(f_Ca, V, time_, Cai);

  // Assignment Rules:
  i_Ca_L := i_CaCa + i_CaK + i_CaNa;
  i_CaCa := d * f * f_Ca * I_CaCa;
  i_CaK := d * f * f_Ca * I_CaK;
  i_CaNa := d * f * f_Ca * I_CaNa;
  I_CaCa := (((P_Ca * power(2, 2) * V * power(F, 2)) / (R * T)) * (gamma_Cai * Cai * exp((2 * V * F) / (R * T)) - gamma_Cao * Cao)) / (exp((2 * V * F) / (R * T)) - 1);
  I_CaK := (((P_K * power(1, 2) * V * power(F, 2)) / (R * T)) * (gamma_Ki * Ki * exp((1 * V * F) / (R * T)) - gamma_Ko * Ko)) / (exp((1 * V * F) / (R * T)) - 1);
  I_CaNa := (((P_Na * power(1, 2) * V * power(F, 2)) / (R * T)) * (gamma_Nai * Nai * exp((1 * V * F) / (R * T)) - gamma_Nao * Nao)) / (exp((1 * V * F) / (R * T)) - 1);

  // Variable initializations:
  gamma_Nai = 0.75;
  gamma_Nao = 0.75;
  gamma_Ki = 0.75;
  gamma_Ko = 0.75;
  P_Ca = 0.00054;
  P_Na = 0.000000675;
  P_K = 0.000000193;
  gamma_Cai = 1;
  gamma_Cao = 0.34;
  R = ;
  T = ;
  F = ;
  Cao = ;
  Nao = ;
  Ko = ;
  Nai = ;
  Ki = ;
end

model zeng_1995__T_type_Ca_channel_b_gate(b, time_, V)

  // Assignment Rules:
  b_infinity := 1 / (1 + exp(-(V + 14) / 10.8));
  tau_b := (3.7 + 6.1) / (1 + exp((25 + V) / 4.5));

  // Rate Rules:
  b' = (b_infinity - b) / tau_b;

  // Variable initializations:
  b = ;
  time_ = ;
  V = ;
end

model zeng_1995__T_type_Ca_channel_g_gate(g, time_, V)

  // Assignment Rules:
  g_infinity := 1 / (1 + exp((V + 60) / 5.6));
  tau_g := piecewise( 12 , V > 0 , -(0.875) * V + 12 );

  // Rate Rules:
  g' = (g_infinity - g) / tau_g;

  // Variable initializations:
  g = ;
  time_ = ;
  V = ;
end

model zeng_1995__T_type_Ca_channel(i_Ca_T, time_, V, Cai, Cao, R, T, F, b, g)

  // Sub-modules, and any changes to those submodules:
  T_type_Ca_channel_b_gate: zeng_1995__T_type_Ca_channel_b_gate(b, time_, V);
  T_type_Ca_channel_g_gate: zeng_1995__T_type_Ca_channel_g_gate(g, time_, V);

  // Assignment Rules:
  i_Ca_T := g_Ca_T * power(b, 2) * g * (V - E_Ca);
  E_Ca := ((R * T) / (2 * F)) * ln(Cao / Cai);

  // Variable initializations:
  g_Ca_T = 0.05;
  Cai = ;
  Cao = ;
  R = ;
  T = ;
  F = ;
end

model zeng_1995__rapid_time_dependent_potassium_current_Xr_gate(Xr, V, time_)

  // Assignment Rules:
  Xr_infinity := 1 / (1 + exp(-(V + 21.5) / 7.5));
  tau_Xr := 1 / ((0.00138 * (V + 14.2)) / (1 - exp(-(0.123) * (V + 14.2))) + (0.00061 * (V + 38.9)) / (exp(0.145 * (V + 38.9)) - 1));

  // Rate Rules:
  Xr' = (Xr_infinity - Xr) / tau_Xr;

  // Variable initializations:
  Xr = ;
  V = ;
  time_ = ;
end

model zeng_1995__rapid_time_dependent_potassium_current_Rr_gate(Rr, V, time_)

  // Assignment Rules:
  Rr := 1 / (1 + exp((V + 9) / 22.4));

  // Variable initializations:
  V = ;
  time_ = ;
end

model zeng_1995__rapid_time_dependent_potassium_current(i_Kr, time_, V, R, T, F, Ko, Ki, Xr, Rr)

  // Sub-modules, and any changes to those submodules:
  rapid_time_dependent_potassium_current_Xr_gate: zeng_1995__rapid_time_dependent_potassium_current_Xr_gate(Xr, V, time_);
  rapid_time_dependent_potassium_current_Rr_gate: zeng_1995__rapid_time_dependent_potassium_current_Rr_gate(Rr, V, time_);

  // Assignment Rules:
  i_Kr := g_Kr * Xr * Rr * (V - E_Kr);
  g_Kr := 0.02614 * root(Ko / 5.4);
  E_Kr := ((R * T) / F) * ln(Ko / Ki);

  // Variable initializations:
  R = ;
  T = ;
  F = ;
  Ko = ;
  Ki = ;
end

model zeng_1995__slow_time_dependent_potassium_current_Xs_gate(Xs, tau_Xs, V, time_)

  // Assignment Rules:
  tau_Xs := 1 / ((0.0000719 * (V + 30)) / (1 - exp(-(0.148) * (V + 30))) + (0.000131 * (V + 30)) / (exp(0.0687 * (V + 30)) - 1));
  Xs_infinity := 1 / (1 + exp(-(V - 1.5) / 16.7));

  // Rate Rules:
  Xs' = (Xs_infinity - Xs) / tau_Xs;

  // Variable initializations:
  Xs = ;
  V = ;
  time_ = ;
end

model zeng_1995__slow_time_dependent_potassium_current(i_Ks, time_, V, Ko, Ki, Nao, Nai, Cai, R, T, F, Xs)

  // Sub-modules, and any changes to those submodules:
  slow_time_dependent_potassium_current_Xs_gate: zeng_1995__slow_time_dependent_potassium_current_Xs_gate(Xs, tau_Xs, V, time_);

  // Assignment Rules:
  i_Ks := g_Ks * power(Xs, 2) * (V - E_Ks);
  g_Ks := (0.057 + 0.19) / (1 + exp((-(7.2) + P_Ca) / 0.6));
  E_Ks := ((R * T) / F) * ln((Ko + P_NaK * Nao) / (Ki + P_NaK * Nai));
  P_Ca := -log(Cai) + 3;

  // Variable initializations:
  P_NaK = 0.01833;
  Ko = ;
  Ki = ;
  Nao = ;
  Nai = ;
  Cai = ;
  R = ;
  T = ;
  F = ;
end

model zeng_1995__time_independent_potassium_current_K1_gate(K1_infinity, V, time_, E_K1)

  // Assignment Rules:
  K1_infinity := alpha_K1 / (alpha_K1 + beta_K1);
  alpha_K1 := 1.02 / (1 + exp(0.2385 * ((V - E_K1) - 59.215)));
  beta_K1 := (0.49124 * exp(0.08032 * ((V + 5.476) - E_K1)) + exp(0.06175 * (V - (E_K1 + 594.31)))) / (1 + exp(-(0.5143) * (V - (E_K1 + 4.753))));

  // Variable initializations:
  V = ;
  time_ = ;
  E_K1 = ;
end

model zeng_1995__time_independent_potassium_current(i_K1, E_K1, time_, V, Ko, Ki, R, T, F, K1_infinity)

  // Sub-modules, and any changes to those submodules:
  time_independent_potassium_current_K1_gate: zeng_1995__time_independent_potassium_current_K1_gate(K1_infinity, V, time_, E_K1);

  // Assignment Rules:
  i_K1 := g_K1 * K1_infinity * (V - E_K1);
  E_K1 := ((R * T) / F) * ln(Ko / Ki);
  g_K1 := 0.75 * root(Ko / 5.4);

  // Variable initializations:
  Ko = ;
  Ki = ;
  R = ;
  T = ;
  F = ;
end

model zeng_1995__plateau_potassium_current(i_Kp, time_, V, E_K1)

  // Assignment Rules:
  i_Kp := g_Kp * Kp * (V - E_Kp);
  E_Kp := E_K1;
  Kp := 1 / (1 + exp((7.488 - V) / 5.98));

  // Variable initializations:
  g_Kp = 0.00552;
  time_ = ;
  V = ;
  E_K1 = ;
end

model zeng_1995__sarcolemmal_calcium_pump(i_p_Ca, time_, Cai)

  // Assignment Rules:
  i_p_Ca := (I_pCa * Cai) / (K_mpCa + Cai);

  // Variable initializations:
  K_mpCa = 0.5;
  I_pCa = 1.15;
  time_ = ;
  Cai = ;
end

model zeng_1995__sodium_background_current(i_Na_b, time_, V, E_Na)

  // Assignment Rules:
  i_Na_b := g_Nab * (V - E_NaN);
  E_NaN := E_Na;

  // Variable initializations:
  g_Nab = 0.00141;
  time_ = ;
  V = ;
  E_Na = ;
end

model zeng_1995__calcium_background_current(i_Ca_b, time_, V, R, T, F, Cai, Cao)

  // Assignment Rules:
  i_Ca_b := g_Cab * (V - E_CaN);
  E_CaN := ((R * T) / (2 * F)) * ln(Cao / Cai);

  // Variable initializations:
  g_Cab = 0.003016;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Cai = ;
  Cao = ;
end

model zeng_1995__sodium_potassium_pump(i_NaK, time_, V, R, T, F, Nai, Nao, Ko)

  // Assignment Rules:
  i_NaK := (((I_NaK * f_NaK * 1) / (1 + root(power(K_mNai / Nai, 3)))) * Ko) / (Ko + K_mKo);
  f_NaK := 1 / (1 + 0.1245 * exp((-(0.1) * V * F) / (R * T)) + 0.0365 * sigma * exp((- V * F) / (R * T)));
  sigma := (1 / 7) * (exp(Nao / 67.3) - 1);

  // Variable initializations:
  I_NaK = ;
  K_mNai = 10;
  K_mKo = 1.5;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Nai = ;
  Nao = ;
  Ko = ;
end

model zeng_1995__non_specific_calcium_activated_current(i_ns_Ca, i_ns_Na, i_ns_K, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, R, T, F, Cao, Nao, Ko, Nai, Ki, time_, Cai, V)

  // Assignment Rules:
  i_ns_Ca := i_ns_Na + i_ns_K;
  i_ns_Na := (I_ns_Na * 1) / (1 + power(K_m_ns_Ca / Cai, 3));
  i_ns_K := (I_ns_K * 1) / (1 + power(K_m_ns_Ca / Cai, 3));
  I_ns_Na := (((P_ns_Ca * power(1, 2) * V * power(F, 2)) / (R * T)) * (gamma_Nai * Nai * exp((1 * V * F) / (R * T)) - gamma_Nao * Nao)) / (exp((1 * V * F) / (R * T)) - 1);
  I_ns_K := (((P_ns_Ca * power(1, 2) * V * power(F, 2)) / (R * T)) * (gamma_Ki * Ki * exp((1 * V * F) / (R * T)) - gamma_Ko * Ko)) / (exp((1 * V * F) / (R * T)) - 1);

  // Variable initializations:
  P_ns_Ca = 1.75e-7;
  gamma_Nai = ;
  gamma_Nao = ;
  gamma_Ki = ;
  gamma_Ko = ;
  R = ;
  T = ;
  F = ;
  Cao = ;
  Nao = ;
  Ko = ;
  Nai = ;
  Ki = ;
  K_m_ns_Ca = 1.2;
  time_ = ;
  Cai = ;
  V = ;
end

model zeng_1995__Na_Ca_exchanger(i_NaCa, time_, V, R, T, F, Nai, Nao, Cai, Cao)

  // Assignment Rules:
  i_NaCa := ((((((K_NaCa * 1) / (power(K_mNa, 3) + power(Nao, 3))) * 1) / (K_mCa + Cao)) * 1) / (1 + K_sat * exp(((eta - 1) * V * F) / (R * T)))) * (exp((eta * V * F) / (R * T)) * power(Nai, 3) * Cao - exp(((eta - 1) * V * F) / (R * T)) * power(Nao, 3) * Cai);

  // Variable initializations:
  K_NaCa = 2000;
  K_mNa = 87.5;
  K_mCa = 1.38;
  K_sat = 0.1;
  eta = 0.35;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Nai = ;
  Nao = ;
  Cai = ;
  Cao = ;
end

model zeng_1995__calcium_buffers_in_the_myoplasm(Ca_JSR_new, Cai, time_, CSQN_max, K_mCSQN, CSQN_buff)

  // Assignment Rules:
  Ca_JSR_new := (root(power(b1, 2) + 4 * c1) - b1) / 2;
  Cai := (2 / 3) * root(power(b, 2) - 3 * c) * cos(arccos((9 * b * c - (2 * power(b, 3) + 27 * d)) / (2 * power(power(b, 2) - 3 * c, 3 / 2))) / 3) - b / 3;
  TRPN_buff := (TRPN_max * Cai) / (Cai + K_mTRPN);
  CMDN_buff := (CMDN_max * Cai) / (Cai + K_mCMDN);
  delta_Ca_JSR := Ca_JSR_old - Ca_JSR_new;
  delta_Cai := Cai_old - Cai;
  b := (CMDN_max + TRPN_max + K_mTRPN + K_mCMDN) - Ca_total;
  c := (K_mCMDN * K_mTRPN + TRPN_max * K_mCMDN + CMDN_max * K_mTRPN) - Ca_total * (K_mTRPN + K_mCMDN);
  d := - K_mTRPN * K_mCMDN * Ca_total;

  // Variable initializations:
  K_mTRPN = 0.5;
  K_mCMDN = 2.38;
  TRPN_max = 70;
  CMDN_max = 50;
  Ca_total = ;
  Ca_JSR_old = ;
  Cai_old = ;
  b1 = ;
  c1 = ;
  time_ = ;
  CSQN_max = ;
  K_mCSQN = ;
  CSQN_buff = ;
end

model zeng_1995__calcium_fluxes_in_the_SR(i_rel, i_up, i_leak, i_tr, CSQN_buff, CSQN_max, K_mCSQN, time_, Cai, Ca_JSR_new, Ca_NSR)

  // Assignment Rules:
  i_rel := G_rel * (Ca_JSR_new - Cai);
  i_up := (I_up * Cai) / (Cai + K_mup);
  i_leak := K_leak * Ca_NSR;
  i_tr := (Ca_NSR - Ca_JSR_new) / tau_tr;
  CSQN_buff := (CSQN_max * Ca_JSR_new) / (Ca_JSR_new + K_mCSQN);
  G_rel := piecewise( ((G_rel_max * (delta_Ca_i2 - delta_Ca_ith)) / ((K_mrel + delta_Ca_i2) - delta_Ca_ith)) * (1 - exp(- t / tau_on)) * exp(- t / tau_off) , calcium_overload == 0 , G_rel_max * (1 - exp(- t / tau_on)) * exp(- t / tau_off) );
  G_rel_max := piecewise( 0 , ( (calcium_overload == 0)) && (delta_Ca_i2 < delta_Ca_ith ), 60 , ( (calcium_overload == 0)) && (geq(delta_Ca_i2, delta_Ca_ith) ), 0 , CSQN_buff < CSQN_th , 4 );
  K_leak := I_up / Ca_NSR_max;

  // Variable initializations:
  CSQN_max = 10;
  K_mCSQN = 0.8;
  tau_on = 2;
  tau_off = 2;
  tau_tr = 180;
  t = 0;
  K_mrel = 0.8;
  delta_Ca_i2 = ;
  delta_Ca_ith = 0.18;
  CSQN_th = 0.7;
  K_mup = 0.92;
  I_up = 0.005;
  Ca_NSR_max = 0.15;
  calcium_overload = ;
  time_ = ;
  Cai = ;
  Ca_JSR_new = ;
  Ca_NSR = ;
end

model zeng_1995__ionic_concentrations(Nai, Nao, Cao, Ki, Ko, Ca_NSR, time_, F, i_Na, i_CaNa, i_CaCa, i_Na_b, i_ns_Na, i_NaCa, i_NaK, i_CaK, i_Kr, i_Ks, i_K1, i_Kp, i_ns_K, i_tr, i_leak, i_up)

  // Rate Rules:
  Nai' = (-(i_Na + i_CaNa + i_Na_b + i_ns_Na + i_NaCa * 3 + i_NaK * 3) * A_cap) / (V_myo * F);
  Ki' = (-(i_CaK + i_Kr + i_Ks + i_K1 + i_Kp + i_ns_K + - i_NaK * 2) * A_cap) / (V_myo * F);
  Ko' = ((i_CaK + i_Kr + i_Ks + i_K1 + i_Kp + i_ns_K + - i_NaK * 2) * A_cap) / (V_cleft * F);
  Ca_NSR' = -((i_leak + i_tr) - i_up);
  Ca_foot' = ((- i_CaCa * A_cap) / (2 * V_myo * F)) * R_A_V;

  // Variable initializations:
  Nai = 10;
  Nao = 140;
  Cao = 1.8;
  Ki = 145;
  Ko = 5.4;
  Ca_NSR = 15;
  Ca_foot = ;
  A_cap = 0.000153;
  R_A_V = ;
  V_myo = ;
  V_cleft = ;
  time_ = ;
  F = ;
  i_Na = ;
  i_CaNa = ;
  i_CaCa = ;
  i_Na_b = ;
  i_ns_Na = ;
  i_NaCa = ;
  i_NaK = ;
  i_CaK = ;
  i_Kr = ;
  i_Ks = ;
  i_K1 = ;
  i_Kp = ;
  i_ns_K = ;
  i_tr = ;
  i_leak = ;
  i_up = ;
end

model *zeng_1995____main()

  // Sub-modules, and any changes to those submodules:
  environment: zeng_1995__environment(time_);
  membrane: zeng_1995__membrane(V, R, T, F, time_, i_Na, i_Ca_L, i_Ca_T, i_Kr, i_Ks, i_NaCa, i_K1, i_Kp, i_p_Ca, i_Na_b, i_Ca_b, i_NaK, i_ns_Ca);
  fast_sodium_current: zeng_1995__fast_sodium_current(i_Na, E_Na, time_, V, R, F, T, Nai, Nao, m, h, j);
  L_type_Ca_channel: zeng_1995__L_type_Ca_channel(i_Ca_L, i_CaCa, i_CaK, i_CaNa, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, time_, V, Cai, R, T, F, Cao, Nao, Ko, Nai, Ki, d, f, f_Ca);
  T_type_Ca_channel: zeng_1995__T_type_Ca_channel(i_Ca_T, time_, V, Cai, Cao, R, T, F, b, g);
  rapid_time_dependent_potassium_current: zeng_1995__rapid_time_dependent_potassium_current(i_Kr, time_, V, R, T, F, Ko, Ki, Xr, Rr);
  slow_time_dependent_potassium_current: zeng_1995__slow_time_dependent_potassium_current(i_Ks, time_, V, Ko, Ki, Nao, Nai, Cai, R, T, F, Xs);
  time_independent_potassium_current: zeng_1995__time_independent_potassium_current(i_K1, E_K1, time_, V, Ko, Ki, R, T, F, K1_infinity);
  plateau_potassium_current: zeng_1995__plateau_potassium_current(i_Kp, time_, V, E_K1);
  sarcolemmal_calcium_pump: zeng_1995__sarcolemmal_calcium_pump(i_p_Ca, time_, Cai);
  sodium_background_current: zeng_1995__sodium_background_current(i_Na_b, time_, V, E_Na);
  calcium_background_current: zeng_1995__calcium_background_current(i_Ca_b, time_, V, R, T, F, Cai, Cao);
  sodium_potassium_pump: zeng_1995__sodium_potassium_pump(i_NaK, time_, V, R, T, F, Nai, Nao, Ko);
  non_specific_calcium_activated_current: zeng_1995__non_specific_calcium_activated_current(i_ns_Ca, i_ns_Na, i_ns_K, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, R, T, F, Cao, Nao, Ko, Nai, Ki, time_, Cai, V);
  Na_Ca_exchanger: zeng_1995__Na_Ca_exchanger(i_NaCa, time_, V, R, T, F, Nai, Nao, Cai, Cao);
  calcium_buffers_in_the_myoplasm: zeng_1995__calcium_buffers_in_the_myoplasm(Ca_JSR_new, Cai, time_, CSQN_max, K_mCSQN, CSQN_buff);
  calcium_fluxes_in_the_SR: zeng_1995__calcium_fluxes_in_the_SR(i_rel, i_up, i_leak, i_tr, CSQN_buff, CSQN_max, K_mCSQN, time_, Cai, Ca_JSR_new, Ca_NSR);
  ionic_concentrations: zeng_1995__ionic_concentrations(Nai, Nao, Cao, Ki, Ko, Ca_NSR, time_, F, i_Na, i_CaNa, i_CaCa, i_Na_b, i_ns_Na, i_NaCa, i_NaK, i_CaK, i_Kr, i_Ks, i_K1, i_Kp, i_ns_K, i_tr, i_leak, i_up);
end