//Created by libAntimony v2.4
model winslow_model_1999__environment(time_)

  // Variable initializations:
  time_ = ;
end

model winslow_model_1999__membrane(V, R, T, F, C_sc, time_, i_Na, i_Ca, i_Ca_K, i_Kr, i_Ks, i_to1, i_K1, i_Kp, i_NaCa, i_NaK, i_p_Ca, i_Ca_b, i_Na_b, i_Stim)

  // Assignment Rules:
  i_Stim := piecewise( stim_amplitude , (( geq(time_, stim_start)) && (time_ <= stim_end) && ((time_ - stim_start) - floor((time_ - stim_start) / stim_period) * stim_period <= stim_duration )), 0 );

  // Rate Rules:
  V' = (-(1) * 1 * (i_Na + i_Ca + i_Ca_K + i_Kr + i_Ks + i_to1 + i_K1 + i_Kp + i_NaCa + i_NaK + i_p_Ca + i_Na_b + i_Ca_b + i_Stim)) / C_sc;

  // Variable initializations:
  V = -96.1638;
  R = 8.314472;
  T = 310;
  F = 96.4853415;
  C_sc = 0.001;
  time_ = ;
  i_Na = ;
  i_Ca = ;
  i_Ca_K = ;
  i_Kr = ;
  i_Ks = ;
  i_to1 = ;
  i_K1 = ;
  i_Kp = ;
  i_NaCa = ;
  i_NaK = ;
  i_p_Ca = ;
  i_Ca_b = ;
  i_Na_b = ;
  stim_start = 0.1;
  stim_end = 100000000;
  stim_period = 1;
  stim_duration = 0.002;
  stim_amplitude = -21.1268;
end

model winslow_model_1999__fast_sodium_current_m_gate(m, V, time_)

  // Assignment Rules:
  alpha_m := piecewise( 1000 / (0.1 - 0.005 * E0_m) , abs(E0_m) < 0.00001 , (320 * E0_m) / (1 - exp(-(0.1) * E0_m)) );
  beta_m := 80 * exp(- V / 11);
  E0_m := V + 47.13;

  // Rate Rules:
  m' = piecewise( alpha_m * (1 - m) - beta_m * m , geq(V, -(90)) , 0 );

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

model winslow_model_1999__fast_sodium_current_h_gate(h, V, time_)

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

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

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

model winslow_model_1999__fast_sodium_current_j_gate(j, V, time_)

  // Assignment Rules:
  alpha_j := piecewise( (1000 * -(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( (121.2 * exp(-(0.01052) * V)) / (1 + exp(-(0.1378) * (V + 40.14))) , V < -(40) , (300 * exp(-(0.0000002535) * V)) / (1 + exp(-(0.1) * (V + 32))) );

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

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

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

  // Sub-modules, and any changes to those submodules:
  fast_sodium_current_m_gate: winslow_model_1999__fast_sodium_current_m_gate(m, V, time_);
  fast_sodium_current_h_gate: winslow_model_1999__fast_sodium_current_h_gate(h, V, time_);
  fast_sodium_current_j_gate: winslow_model_1999__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 = 12.8;
  Nao = ;
  Nai = ;
  R = ;
  F = ;
  T = ;
end

model winslow_model_1999__rapid_activating_delayed_rectifiyer_K_current_X_kr_gate(X_kr, V, time_)

  // Assignment Rules:
  K12 := exp(-(5.495) + 0.1691 * V);
  K21 := exp(-(7.677) - 0.0128 * V);
  X_kr_inf := K12 / (K12 + K21);
  tau_X_kr := 0.001 / (K12 + K21) + tau_factor * 0.027;

  // Rate Rules:
  X_kr' = (X_kr_inf - X_kr) / tau_X_kr;

  // Variable initializations:
  X_kr = 0.51;
  tau_factor = 1;
  V = ;
  time_ = ;
end

model winslow_model_1999__rapid_activating_delayed_rectifiyer_K_current(i_Kr, E_K, time_, V, Ko, Ki, R, T, F, X_kr)

  // Sub-modules, and any changes to those submodules:
  rapid_activating_delayed_rectifiyer_K_current_X_kr_gate: winslow_model_1999__rapid_activating_delayed_rectifiyer_K_current_X_kr_gate(X_kr, V, time_);

  // Assignment Rules:
  i_Kr := g_Kr * f_Ko * R_V * X_kr * (V - E_K);
  E_K := ((R * T) / F) * ln(Ko / Ki);
  f_Ko := root(Ko / 4);
  R_V := 1 / (1 + 1.4945 * exp(0.0446 * V));

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

model winslow_model_1999__slow_activating_delayed_rectifiyer_K_current_X_ks_gate(X_ks, V, time_)

  // Assignment Rules:
  tau_X_ks := 0.001 / ((0.0000719 * (V - 10)) / (1 - exp(-(0.148) * (V - 10))) + (0.000131 * (V - 10)) / (exp(0.0687 * (V - 10)) - 1));
  X_ks_infinity := 1 / (1 + exp(-(V - 24.7) / 13.6));

  // Rate Rules:
  X_ks' = (X_ks_infinity - X_ks) / tau_X_ks;

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

model winslow_model_1999__slow_activating_delayed_rectifiyer_K_current(i_Ks, time_, V, Ko, Ki, Nao, Nai, R, T, F, X_ks)

  // Sub-modules, and any changes to those submodules:
  slow_activating_delayed_rectifiyer_K_current_X_ks_gate: winslow_model_1999__slow_activating_delayed_rectifiyer_K_current_X_ks_gate(X_ks, V, time_);

  // Assignment Rules:
  i_Ks := g_Ks * power(X_ks, 2) * (V - E_Ks);
  E_Ks := ((R * T) / F) * ln((Ko + 0.01833 * Nao) / (Ki + 0.01833 * Nai));

  // Variable initializations:
  g_Ks = 0.0027134;
  Ko = ;
  Ki = ;
  Nao = ;
  Nai = ;
  R = ;
  T = ;
  F = ;
end

model winslow_model_1999__transient_outward_potassium_current_X_to1_gate(X_to1, V, time_)

  // Assignment Rules:
  alpha_X_to1 := 45.16 * exp(0.03577 * V);
  beta_X_to1 := 98.9 * exp(-(0.06237) * V);

  // Rate Rules:
  X_to1' = alpha_X_to1 * (1 - X_to1) - beta_X_to1 * X_to1;

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

model winslow_model_1999__transient_outward_potassium_current_Y_to1_gate(Y_to1, V, time_)

  // Assignment Rules:
  alpha_Y_to1 := (5.415 * exp(-(V + 33.5) / 5)) / (1 + 0.051335 * exp(-(V + 33.5) / 5));
  beta_Y_to1 := (5.415 * exp((V + 33.5) / 5)) / (1 + 0.051335 * exp((V + 33.5) / 5));

  // Rate Rules:
  Y_to1' = alpha_Y_to1 * (1 - Y_to1) - beta_Y_to1 * Y_to1;

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

model winslow_model_1999__transient_outward_potassium_current(i_to1, time_, V, E_K, X_to1, Y_to1)

  // Sub-modules, and any changes to those submodules:
  transient_outward_potassium_current_X_to1_gate: winslow_model_1999__transient_outward_potassium_current_X_to1_gate(X_to1, V, time_);
  transient_outward_potassium_current_Y_to1_gate: winslow_model_1999__transient_outward_potassium_current_Y_to1_gate(Y_to1, V, time_);

  // Assignment Rules:
  i_to1 := g_to1 * X_to1 * Y_to1 * (V - E_K);

  // Variable initializations:
  g_to1 = 0.23815;
  E_K = ;
end

model winslow_model_1999__time_independent_potassium_current_K1_gate(K1_infinity_V, V, time_, E_K, R, T, F)

  // Assignment Rules:
  K1_infinity_V := 1 / (2 + exp(((1.5 * F) / (R * T)) * (V - E_K)));

  // Variable initializations:
  V = ;
  time_ = ;
  E_K = ;
  R = ;
  T = ;
  F = ;
end

model winslow_model_1999__time_independent_potassium_current(i_K1, time_, V, R, T, F, E_K, Ko, K1_infinity_V)

  // Sub-modules, and any changes to those submodules:
  time_independent_potassium_current_K1_gate: winslow_model_1999__time_independent_potassium_current_K1_gate(K1_infinity_V, V, time_, E_K, R, T, F);

  // Assignment Rules:
  i_K1 := ((g_K1 * K1_infinity_V * Ko) / (Ko + K_mK1)) * (V - E_K);

  // Variable initializations:
  g_K1 = 2.8;
  K_mK1 = 13;
  Ko = ;
end

model winslow_model_1999__plateau_potassium_current_Kp_gate(Kp_V, V, time_)

  // Assignment Rules:
  Kp_V := 1 / (1 + exp((7.488 - V) / 5.98));

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

model winslow_model_1999__plateau_potassium_current(i_Kp, time_, V, E_K, Kp_V)

  // Sub-modules, and any changes to those submodules:
  plateau_potassium_current_Kp_gate: winslow_model_1999__plateau_potassium_current_Kp_gate(Kp_V, V, time_);

  // Assignment Rules:
  i_Kp := g_Kp * Kp_V * (V - E_K);

  // Variable initializations:
  g_Kp = 0.002216;
  E_K = ;
end

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

  // Assignment Rules:
  i_NaCa := ((K_NaCa * 5000) / ((power(K_mNa, 3) + power(Nao, 3)) * (K_mCa + Cao) * (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_mCa = 1.38;
  K_mNa = 87.5;
  K_NaCa = 0.3;
  K_sat = 0.2;
  eta = 0.35;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Nai = ;
  Nao = ;
  Cai = ;
  Cao = ;
end

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

  // Assignment Rules:
  i_NaK := (((I_NaK * f_NaK) / (1 + power(K_mNai / Nai, 1.5))) * 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 = 0.693;
  K_mNai = 10;
  K_mKo = 1.5;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Nai = ;
  Nao = ;
  Ko = ;
end

model winslow_model_1999__sarcolemmal_calcium_pump(i_p_Ca, time_, Cai)

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

  // Variable initializations:
  K_mpCa = 0.00005;
  I_pCa = 0.05;
  time_ = ;
  Cai = ;
end

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

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

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

model winslow_model_1999__sodium_background_current(i_Na_b, time_, V, E_Na)

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

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

model winslow_model_1999__L_type_Ca_current_y_gate(y, V, time_)

  // Assignment Rules:
  y_infinity := 0.8 / (1 + exp((V + 12.5) / 5)) + 0.2;
  tau_y := (20 + 600 / (1 + exp((V + 20) / 9.5))) / 1000;

  // Rate Rules:
  y' = (y_infinity - y) / tau_y;

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

model winslow_model_1999__L_type_Ca_current(i_Ca, i_Ca_K, time_, V, C_sc, Ca_ss, Cao, Ko, Ki, R, T, F, y)

  // Sub-modules, and any changes to those submodules:
  L_type_Ca_current_y_gate: winslow_model_1999__L_type_Ca_current_y_gate(y, V, time_);

  // Assignment Rules:
  i_Ca := i_Ca_max * y * (O + O_Ca);
  i_Ca_K := ((((p_prime_k / (1 * 1)) * y * (O + O_Ca) * V * power(F, 2)) / (R * T)) * (Ki * exp((V * F) / (R * T)) - Ko)) / (exp((V * F) / (R * T)) - 1);
  p_prime_k := P_K / (1 + i_Ca_max / i_Ca_half);
  i_Ca_max := ((((P_Ca / (1 * 1)) * 4 * V * power(F, 2) * 1000) / (R * T)) * (0.001 * exp((2 * V * F) / (R * T)) - 0.341 * Cao)) / (exp((2 * V * F) / (R * T)) - 1);
  alpha := 400 * exp((V + 2) / 10);
  beta := 50 * exp(-(V + 2) / 13);
  gamma := (103.75 * Ca_ss) / 1;
  alpha_a := alpha * a;
  beta_b := beta / b;

  // Rate Rules:
  O' = f * C4 - g * O;
  O_Ca' = fprime * C_Ca4 - gprime * O_Ca;
  C0' = (beta * C1 + omega * C_Ca0) - (4 * alpha + gamma) * C0;
  C1' = (4 * alpha * C0 + 2 * beta * C2 + (omega / b) * C_Ca1) - (beta + 3 * alpha + gamma * a) * C1;
  C2' = (3 * alpha * C1 + 3 * beta * C3 + (omega / power(b, 2)) * C_Ca2) - (beta * 2 + 2 * alpha + gamma * power(a, 2)) * C2;
  C3' = (2 * alpha * C2 + 4 * beta * C4 + (omega / power(b, 3)) * C_Ca3) - (beta * 3 + alpha + gamma * power(a, 3)) * C3;
  C4' = (alpha * C3 + g * O + (omega / power(b, 4)) * C_Ca4) - (beta * 4 + f + gamma * power(a, 4)) * C4;
  C_Ca0' = (beta_b * C_Ca1 + gamma * C0) - (4 * alpha_a + omega) * C_Ca0;
  C_Ca1' = (4 * alpha_a * C_Ca0 + 2 * beta_b * C_Ca2 + gamma * a * C1) - (beta_b + 3 * alpha_a + omega / b) * C_Ca1;
  C_Ca2' = (3 * alpha_a * C_Ca1 + 3 * beta_b * C_Ca3 + gamma * power(a, 2) * C2) - (beta_b * 2 + 2 * alpha_a + omega / power(b, 2)) * C_Ca2;
  C_Ca3' = (2 * alpha_a * C_Ca2 + 4 * beta_b * C_Ca4 + gamma * power(a, 3) * C3) - (beta_b * 3 + alpha_a + omega / power(b, 3)) * C_Ca3;
  C_Ca4' = (alpha_a * C_Ca3 + gprime * O_Ca + gamma * power(a, 4) * C4) - (beta_b * 4 + fprime + omega / power(b, 4)) * C_Ca4;

  // Variable initializations:
  P_Ca = 3.125e-4;
  P_K = 5.79e-7;
  i_Ca_half = -0.265;
  O = 9.84546e-21;
  O_Ca = 0;
  a = 2;
  b = 2;
  g = 2000;
  f = 300;
  gprime = 7000;
  fprime = 7;
  omega = 10;
  C0 = 0.997208;
  C1 = 6.38897e-5;
  C2 = 1.535e-9;
  C3 = 1.63909e-14;
  C4 = 6.56337e-20;
  C_Ca0 = 0.00272826;
  C_Ca1 = 6.99215e-7;
  C_Ca2 = 6.71989e-11;
  C_Ca3 = 2.87031e-15;
  C_Ca4 = 4.59752e-20;
  C_sc = ;
  Ca_ss = ;
  Cao = ;
  Ko = ;
  Ki = ;
  R = ;
  T = ;
  F = ;
end

model winslow_model_1999__RyR_channel(J_rel, time_, Ca_ss, Ca_JSR)

  // Assignment Rules:
  J_rel := v1 * (P_O1 + P_O2) * (Ca_JSR - Ca_ss);

  // Rate Rules:
  P_O1' = (k_a_plus * power(Ca_ss, n) * P_C1 - (k_a_minus * P_O1 + k_b_plus * power(Ca_ss, m) * P_O1 + k_c_plus * P_O1)) + k_b_minus * P_O2 + k_c_minus * P_C2;
  P_O2' = k_b_plus * power(Ca_ss, m) * P_O1 - k_b_minus * P_O2;
  P_C1' = - k_a_plus * power(Ca_ss, n) * P_C1 + k_a_minus * P_O1;
  P_C2' = k_c_plus * P_O1 - k_c_minus * P_C2;

  // Variable initializations:
  v1 = 1800;
  k_a_plus = 1.215e13;
  k_a_minus = 576;
  k_b_plus = 4.05e9;
  k_b_minus = 1930;
  k_c_plus = 100;
  k_c_minus = 0.8;
  P_O1 = 0;
  P_O2 = 0;
  P_C1 = 0.47;
  P_C2 = 0.53;
  n = 4;
  m = 3;
  time_ = ;
  Ca_ss = ;
  Ca_JSR = ;
end

model winslow_model_1999__SERCA2a_pump(J_up, time_, Cai, Ca_NSR)

  // Assignment Rules:
  J_up := (K_SR * (Vmaxf * fb - Vmaxr * rb)) / (1 + fb + rb);
  fb := power(Cai / K_fb, N_fb);
  rb := power(Ca_NSR / K_rb, N_rb);

  // Variable initializations:
  K_fb = 0.000168;
  K_rb = 3.29;
  Vmaxf = 0.0813;
  Vmaxr = 0.318;
  K_SR = 1;
  N_fb = 1.2;
  N_rb = 1;
  time_ = ;
  Cai = ;
  Ca_NSR = ;
end

model winslow_model_1999__intracellular_Ca_fluxes(J_tr, J_xfer, J_trpn, time_, Ca_ss, Cai, Ca_NSR, Ca_JSR)

  // Assignment Rules:
  J_tr := (Ca_NSR - Ca_JSR) / tau_tr;
  J_xfer := (Ca_ss - Cai) / tau_xfer;
  J_trpn := HTRPN_tot * J_HTRPNCa + LTRPN_tot * J_LTRPNCa;
  J_HTRPNCa := d(HTRPNCa)/ d(unknown- element);
  J_LTRPNCa := d(LTRPNCa)/ d(unknown- element);

  // Rate Rules:
  HTRPNCa' = k_htrpn_plus * Cai * (1 - HTRPNCa) - k_htrpn_minus * HTRPNCa;
  LTRPNCa' = k_ltrpn_plus * Cai * (1 - LTRPNCa) - k_ltrpn_minus * LTRPNCa;

  // Variable initializations:
  tau_tr = 0.0005747;
  tau_xfer = 0.0267;
  HTRPNCa = 0.98;
  LTRPNCa = 0.078;
  HTRPN_tot = 0.14;
  LTRPN_tot = 0.07;
  k_htrpn_plus = 20000;
  k_htrpn_minus = 0.066;
  k_ltrpn_plus = 40000;
  k_ltrpn_minus = 40;
  time_ = ;
  Ca_ss = ;
  Cai = ;
  Ca_NSR = ;
  Ca_JSR = ;
  d = ;
  unknown = ;
  element = ;
end

model winslow_model_1999__intracellular_ion_concentrations(Nai, Cai, Ki, Ca_ss, Ca_JSR, Ca_NSR, time_, F, C_sc, i_Na, i_Ca, i_Na_b, i_NaCa, i_NaK, i_Ca_K, i_Kr, i_Ks, i_K1, i_Kp, i_to1, i_p_Ca, i_Ca_b, J_up, J_rel, J_xfer, J_trpn, J_tr)

  // Assignment Rules:
  beta_i := 1 / (1 + (CMDN_tot * K_mCMDN) / power(K_mCMDN + Cai, 2) + (EGTA_tot * K_mEGTA) / power(K_mEGTA + Cai, 2));
  beta_SS := 1 / (1 + (CMDN_tot * K_mCMDN) / power(K_mCMDN + Ca_ss, 2) + (EGTA_tot * K_mEGTA) / power(K_mEGTA + Ca_ss, 2));
  beta_JSR := 1 / (1 + (CSQN_tot * K_mCSQN) / power(K_mCSQN + Ca_JSR, 2));

  // Rate Rules:
  Nai' = (-(0) * (i_Na + i_Na_b + i_NaCa * 3 + i_NaK * 3) * A_cap * 1) / (V_myo * F);
  Cai' = beta_i * ((J_xfer - (J_up + J_trpn)) + ((2 * i_NaCa - (i_p_Ca + i_Ca_b)) * A_cap * 1) / (2 * V_myo * F));
  Ki' = (-(0) * (i_Ca_K + i_Kr + i_Ks + i_K1 + i_Kp + i_to1 + i_NaK * -(2)) * A_cap * 1) / (V_myo * F);
  Ca_ss' = beta_SS * (((J_rel * V_JSR) / V_SS - (J_xfer * V_myo) / V_SS) - (i_Ca * A_cap * 1) / (2 * V_SS * F));
  Ca_JSR' = beta_JSR * (J_tr - J_rel);
  Ca_NSR' = (J_up * V_myo) / V_NSR - (J_tr * V_JSR) / V_NSR;

  // Variable initializations:
  Nai = 10;
  Cai = 0.00008;
  Ki = 157.8;
  Ca_ss = 0.00011;
  Ca_JSR = 0.257;
  Ca_NSR = 0.257;
  A_cap = 0.0001534;
  V_myo = 0.00002584;
  V_JSR = 0.00000016;
  V_NSR = 0.0000021;
  V_SS = 0.0000000012;
  K_mCMDN = 0.00238;
  K_mEGTA = 0.00015;
  K_mCSQN = 0.8;
  CMDN_tot = 0.05;
  EGTA_tot = 0;
  CSQN_tot = 15;
  time_ = ;
  F = ;
  C_sc = ;
  i_Na = ;
  i_Ca = ;
  i_Na_b = ;
  i_NaCa = ;
  i_NaK = ;
  i_Ca_K = ;
  i_Kr = ;
  i_Ks = ;
  i_K1 = ;
  i_Kp = ;
  i_to1 = ;
  i_p_Ca = ;
  i_Ca_b = ;
  J_up = ;
  J_rel = ;
  J_xfer = ;
  J_trpn = ;
  J_tr = ;
end

model winslow_model_1999__standard_ionic_concentrations(Nao, Cao, Ko)

  // Variable initializations:
  Nao = 138;
  Cao = 2;
  Ko = 4;
end

model *winslow_model_1999____main()

  // Sub-modules, and any changes to those submodules:
  environment: winslow_model_1999__environment(time_);
  membrane: winslow_model_1999__membrane(V, R, T, F, C_sc, time_, i_Na, i_Ca, i_Ca_K, i_Kr, i_Ks, i_to1, i_K1, i_Kp, i_NaCa, i_NaK, i_p_Ca, i_Ca_b, i_Na_b, i_Stim);
  fast_sodium_current: winslow_model_1999__fast_sodium_current(i_Na, E_Na, time_, V, Nao, Nai, R, F, T, m, h, j);
  rapid_activating_delayed_rectifiyer_K_current: winslow_model_1999__rapid_activating_delayed_rectifiyer_K_current(i_Kr, E_K, time_, V, Ko, Ki, R, T, F, X_kr);
  slow_activating_delayed_rectifiyer_K_current: winslow_model_1999__slow_activating_delayed_rectifiyer_K_current(i_Ks, time_, V, Ko, Ki, Nao, Nai, R, T, F, X_ks);
  transient_outward_potassium_current: winslow_model_1999__transient_outward_potassium_current(i_to1, time_, V, E_K, X_to1, Y_to1);
  time_independent_potassium_current: winslow_model_1999__time_independent_potassium_current(i_K1, time_, V, R, T, F, E_K, Ko, K1_infinity_V);
  plateau_potassium_current: winslow_model_1999__plateau_potassium_current(i_Kp, time_, V, E_K, Kp_V);
  Na_Ca_exchanger: winslow_model_1999__Na_Ca_exchanger(i_NaCa, time_, V, R, T, F, Nai, Nao, Cai, Cao);
  sodium_potassium_pump: winslow_model_1999__sodium_potassium_pump(i_NaK, time_, V, R, T, F, Nai, Nao, Ko);
  sarcolemmal_calcium_pump: winslow_model_1999__sarcolemmal_calcium_pump(i_p_Ca, time_, Cai);
  calcium_background_current: winslow_model_1999__calcium_background_current(i_Ca_b, time_, V, R, T, F, Cai, Cao);
  sodium_background_current: winslow_model_1999__sodium_background_current(i_Na_b, time_, V, E_Na);
  L_type_Ca_current: winslow_model_1999__L_type_Ca_current(i_Ca, i_Ca_K, time_, V, C_sc, Ca_ss, Cao, Ko, Ki, R, T, F, y);
  RyR_channel: winslow_model_1999__RyR_channel(J_rel, time_, Ca_ss, Ca_JSR);
  SERCA2a_pump: winslow_model_1999__SERCA2a_pump(J_up, time_, Cai, Ca_NSR);
  intracellular_Ca_fluxes: winslow_model_1999__intracellular_Ca_fluxes(J_tr, J_xfer, J_trpn, time_, Ca_ss, Cai, Ca_NSR, Ca_JSR);
  intracellular_ion_concentrations: winslow_model_1999__intracellular_ion_concentrations(Nai, Cai, Ki, Ca_ss, Ca_JSR, Ca_NSR, time_, F, C_sc, i_Na, i_Ca, i_Na_b, i_NaCa, i_NaK, i_Ca_K, i_Kr, i_Ks, i_K1, i_Kp, i_to1, i_p_Ca, i_Ca_b, J_up, J_rel, J_xfer, J_trpn, J_tr);
  standard_ionic_concentrations: winslow_model_1999__standard_ionic_concentrations(Nao, Cao, Ko);
end