//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 shaw_1997__environment(time_) // Variable initializations: time_ = ; end model shaw_1997__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, i_K_ATP) // Rate Rules: V' = -(1.0 / Cm) * (i_Na + i_Ca_L + i_Ca_T + i_Kr + i_Ks + i_K1 + i_Kp + i_NaCa + i_K_ATP + i_p_Ca + i_Na_b + i_Ca_b + i_NaK + i_ns_Ca + I_st); // Variable initializations: V = -84.624; R = 8314.0; T = 310.0; F = 96500.0; Cm = 1.0; I_st = -100.0; 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 = ; i_K_ATP = ; end model shaw_1997__fast_sodium_current_m_gate(m, V, time_) // Assignment Rules: alpha_m := (0.32 * (V + 47.13)) / (1.0 - exp(-0.1 * (V + 47.13))); beta_m := 0.08 * exp(- V / 11.0); // Rate Rules: m' = alpha_m * (1.0 - m) - beta_m * m; // Variable initializations: m = ; V = ; time_ = ; end model shaw_1997__fast_sodium_current_h_gate(h, V, time_) // Assignment Rules: alpha_h := piecewise( 0.135 * exp((80.0 + V) / -6.8) , V < -40.0 , 0.0 ); beta_h := piecewise( 3.56 * exp(0.079 * V) + 310000.0 * exp(0.35 * V) , V < -40.0 , 1.0 / (0.13 * (1.0 + exp((V + 10.66) / -11.1))) ); // Rate Rules: h' = alpha_h * (1.0 - h) - beta_h * h; // Variable initializations: h = ; V = ; time_ = ; end model shaw_1997__fast_sodium_current_j_gate(j, V, time_) // Assignment Rules: alpha_j := piecewise( (-127140.0 * exp(0.2444 * V) - 0.00003474 * exp(-0.04391 * V)) * ((V + 37.78) / (1.0 + exp(0.311 * (V + 79.23)))) , V < -40.0 , 0.0 ); beta_j := piecewise( (0.1212 * exp(-0.01052 * V)) / (1.0 + exp(-0.1378 * (V + 40.14))) , V < -40.0 , (0.3 * exp(-0.0000002535 * V)) / (1.0 + exp(-0.1 * (V + 32.0))) ); // Rate Rules: j' = alpha_j * (1.0 - j) - beta_j * j; // Variable initializations: j = ; V = ; time_ = ; end model shaw_1997__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: shaw_1997__fast_sodium_current_m_gate(m, V, time_); fast_sodium_current_h_gate: shaw_1997__fast_sodium_current_h_gate(h, V, time_); fast_sodium_current_j_gate: shaw_1997__fast_sodium_current_j_gate(j, V, time_); // Assignment Rules: i_Na := g_Na * power(m, 3.0) * h * j * (V - E_Na); E_Na := ((R * T) / F) * ln(Nao / Nai); // Variable initializations: g_Na = 16.0; R = ; F = ; T = ; Nai = ; Nao = ; end model shaw_1997__L_type_Ca_channel_d_gate(d, V, time_) // Assignment Rules: alpha_d := d_infinity / tau_d; beta_d := (1.0 - d_infinity) / tau_d; d_infinity := 1.0 / (1.0 + exp(-((V + 10.0) / 6.24))); tau_d := d_infinity * ((1.0 - exp(-((V + 10.0) / 6.24))) / (0.035 * (V + 10.0))); // Rate Rules: d' = alpha_d * (1.0 - d) - beta_d * d; // Variable initializations: d = ; V = ; time_ = ; end model shaw_1997__L_type_Ca_channel_f_gate(f, V, time_) // Assignment Rules: alpha_f := f_infinity / tau_f; beta_f := (1.0 - f_infinity) / tau_f; f_infinity := 1.0 / (1.0 + exp((V + 35.06) / 8.6)) + 0.6 / (1.0 + exp((50.0 - V) / 20.0)); tau_f := 1.0 / (0.0197 * exp(-power(0.0337 * (V + 10.0), 2.0)) + 0.02); // Rate Rules: f' = alpha_f * (1.0 - f) - beta_f * f; // Variable initializations: f = ; V = ; time_ = ; end model shaw_1997__L_type_Ca_channel_f_Ca_gate(f_Ca, V, time_, Cai) // Assignment Rules: f_Ca := 1.0 / (1.0 + power(Cai / Km_Ca, 2.0)); // Variable initializations: Km_Ca = 0.6; V = ; time_ = ; Cai = ; end model shaw_1997__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: shaw_1997__L_type_Ca_channel_d_gate(d, V, time_); L_type_Ca_channel_f_gate: shaw_1997__L_type_Ca_channel_f_gate(f, V, time_); L_type_Ca_channel_f_Ca_gate: shaw_1997__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.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Cai * Cai * exp((2.0 * V * F) / (R * T)) - gamma_Cao * Cao) / (exp((2.0 * V * F) / (R * T)) - 1.0)); I_CaK := P_K * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Ki * Ki * exp((1.0 * V * F) / (R * T)) - gamma_Ko * Ko) / (exp((1.0 * V * F) / (R * T)) - 1.0)); I_CaNa := P_Na * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Nai * Nai * exp((1.0 * V * F) / (R * T)) - gamma_Nao * Nao) / (exp((1.0 * V * F) / (R * T)) - 1.0)); // 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.0; gamma_Cao = 0.34; R = ; T = ; F = ; Cao = ; Nao = ; Ko = ; Nai = ; Ki = ; end model shaw_1997__T_type_Ca_channel_b_gate(b, time_, V) // Assignment Rules: b_infinity := 1.0 / (1.0 + exp(-((V + 14.0) / 10.8))); tau_b := (3.7 + 6.1) / (1.0 + exp((25.0 + V) / 4.5)); // Rate Rules: b' = (b_infinity - b) / tau_b; // Variable initializations: b = ; time_ = ; V = ; end model shaw_1997__T_type_Ca_channel_g_gate(g, time_, V) // Assignment Rules: g_infinity := 1.0 / (1.0 + exp((V + 60.0) / 5.6)); tau_g := piecewise( 12.0 , V > 0.0 , -0.875 * V + 12.0 ); // Rate Rules: g' = (g_infinity - g) / tau_g; // Variable initializations: g = ; time_ = ; V = ; end model shaw_1997__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: shaw_1997__T_type_Ca_channel_b_gate(b, time_, V); T_type_Ca_channel_g_gate: shaw_1997__T_type_Ca_channel_g_gate(g, time_, V); // Assignment Rules: i_Ca_T := g_Ca_T * power(b, 2.0) * g * (V - E_Ca); E_Ca := ((R * T) / (2.0 * F)) * ln(Cao / Cai); // Variable initializations: g_Ca_T = 0.05; Cai = ; Cao = ; R = ; T = ; F = ; end model shaw_1997__rapid_time_dependent_potassium_current_Xr_gate(Xr, V, time_) // Assignment Rules: Xr_infinity := 1.0 / (1.0 + exp(-((V + 21.5) / 7.5))); tau_Xr := 1.0 / ((0.00138 * (V + 14.2)) / (1.0 - exp(-0.123 * (V + 14.2))) + (0.00061 * (V + 38.9)) / (exp(0.145 * (V + 38.9)) - 1.0)); // Rate Rules: Xr' = (Xr_infinity - Xr) / tau_Xr; // Variable initializations: Xr = ; V = ; time_ = ; end model shaw_1997__rapid_time_dependent_potassium_current_Rr_gate(Rr, V, time_) // Assignment Rules: Rr := 1.0 / (1.0 + exp((V + 9.0) / 22.4)); // Variable initializations: V = ; time_ = ; end model shaw_1997__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: shaw_1997__rapid_time_dependent_potassium_current_Xr_gate(Xr, V, time_); rapid_time_dependent_potassium_current_Rr_gate: shaw_1997__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 shaw_1997__slow_time_dependent_potassium_current_Xs_gate(Xs, tau_Xs, V, time_) // Assignment Rules: tau_Xs := 1.0 / ((0.0000719 * (V + 30.0)) / (1.0 - exp(-0.148 * (V + 30.0))) + (0.000131 * (V + 30.0)) / (exp(0.0687 * (V + 30.0)) - 1.0)); Xs_infinity := 1.0 / (1.0 + exp(-(V + -1.5) / 16.7)); // Rate Rules: Xs' = (Xs_infinity - Xs) / tau_Xs; // Variable initializations: Xs = ; V = ; time_ = ; end model shaw_1997__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: shaw_1997__slow_time_dependent_potassium_current_Xs_gate(Xs, tau_Xs, V, time_); // Assignment Rules: i_Ks := g_Ks * power(Xs, 2.0) * (V - E_Ks); g_Ks := (0.057 + 0.19) / (1.0 + 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.0; // Variable initializations: P_NaK = 0.01833; Ko = ; Ki = ; Nao = ; Nai = ; Cai = ; R = ; T = ; F = ; end model shaw_1997__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.0 + 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.0 + exp(-0.5143 * (V - (E_K1 + 4.753)))); // Variable initializations: V = ; time_ = ; E_K1 = ; end model shaw_1997__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: shaw_1997__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 shaw_1997__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.0 / (1.0 + exp((7.488 - V) / 5.98)); // Variable initializations: g_Kp = 0.0183; time_ = ; V = ; E_K1 = ; end model shaw_1997__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 shaw_1997__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 shaw_1997__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.0 * F)) * ln(Cao / Cai); // Variable initializations: g_Cab = 0.003016; time_ = ; V = ; R = ; T = ; F = ; Cai = ; Cao = ; end model shaw_1997__sodium_potassium_pump(i_NaK, time_, V, R, T, F, Nai, Nao, Ko) // Assignment Rules: i_NaK := I_NaK * f_NaK * (1.0 / (1.0 + power(K_mNai / Nai, 1.5))) * (Ko / (Ko + K_mKo)); f_NaK := 1.0 / ((1.0 + 0.1245 * exp(-0.1 * ((V * F) / (R * T)))) + 0.0365 * sigma * exp(-((V * F) / (R * T)))); sigma := (1.0 / 7.0) * (exp(Nao / 67.3) - 1.0); // Variable initializations: I_NaK = ; K_mNai = 10.0; K_mKo = 1.5; time_ = ; V = ; R = ; T = ; F = ; Nai = ; Nao = ; Ko = ; end model shaw_1997__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.0 / (1.0 + power(K_m_ns_Ca / Cai, 3.0))); i_ns_K := I_ns_K * (1.0 / (1.0 + power(K_m_ns_Ca / Cai, 3.0))); I_ns_Na := P_ns_Ca * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Nai * Nai * exp((1.0 * V * F) / (R * T)) - gamma_Nao * Nao) / (exp((1.0 * V * F) / (R * T)) - 1.0)); I_ns_K := P_ns_Ca * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Ki * Ki * exp((1.0 * V * F) / (R * T)) - gamma_Ko * Ko) / (exp((1.0 * V * F) / (R * T)) - 1.0)); // 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 shaw_1997__ATP_dependent_potassium_current(i_K_ATP, time_, V, Ko, Ki, R, T, F) // Assignment Rules: i_K_ATP := g_K_ATP * (V - E_K); g_K_ATP := G_K_ATP * P_ATP * power(Ko / Ko_normal, n); G_K_ATP := 0.000195 / Nichols_area; P_ATP := 1.0 / (1.0 + power(ATP_i / K_05, H)); E_K := ((R * T) / F) * ln(Ko / Ki); // Variable initializations: Ko_normal = 4.0; K_05 = 0.250; ATP_i = 3.0; n = 0.24; H = 2.0; Nichols_area = 0.005; time_ = ; V = ; Ko = ; Ki = ; R = ; T = ; F = ; end model shaw_1997__Na_Ca_exchanger(i_NaCa, time_, V, R, T, F, Nai, Nao, Cai, Cao) // Assignment Rules: i_NaCa := K_NaCa * (1.0 / (power(K_mNa, 3.0) + power(Nao, 3.0))) * (1.0 / (K_mCa + Cao)) * (1.0 / (1.0 + K_sat * exp((eta - 1.0) * V * (F / (R * T))))) * (exp(eta * V * (F / (R * T))) * power(Nai, 3.0) * Cao - exp((eta - 1.0) * V * (F / (R * T))) * power(Nao, 3.0) * Cai); // Variable initializations: K_NaCa = 2000.0; 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 shaw_1997__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.0) + 4.0 * c1) - b1) / 2.0; Cai := (2.0 / 3.0) * root(power(b, 2.0) - 3.0 * c) * cos(arccos((9.0 * b * c - (2.0 * power(b, 3.0) + 27.0 * d)) / (2.0 * power(power(b, 2.0) - 3.0 * c, 3.0 / 2.0))) / 3.0) - b / 3.0; 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.0; CMDN_max = 50.0; Ca_total = ; Ca_JSR_old = ; Cai_old = ; b1 = ; c1 = ; time_ = ; CSQN_max = ; K_mCSQN = ; CSQN_buff = ; end model shaw_1997__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.0 - exp(-(t / tau_on))) * exp(-(t / tau_off)) , calcium_overload == 0.0 , G_rel_max * (1.0 - exp(-(t / tau_on))) * exp(-(t / tau_off)) ); G_rel_max := piecewise( piecewise( 0.0 , delta_Ca_i2 < delta_Ca_ith , 60.0 ) , calcium_overload == 0.0 , piecewise( 0.0 , CSQN_buff < CSQN_th , 4.0 ) ); K_leak := I_up / Ca_NSR_max; // Variable initializations: CSQN_max = 10.0; K_mCSQN = 0.8; tau_on = 2.0; tau_off = 2.0; tau_tr = 180.0; t = 0.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 shaw_1997__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_K_ATP, 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.0 + i_NaK * 3.0) * (A_cap / (V_myo * F)); Ki' = -(i_CaK + i_Kr + i_Ks + i_K1 + i_Kp + i_K_ATP + i_ns_K + -(i_NaK * 2.0)) * (A_cap / (V_myo * F)); Ko' = (i_CaK + i_Kr + i_Ks + i_K1 + i_Kp + i_K_ATP + i_ns_K + -(i_NaK * 2.0)) * (A_cap / (V_cleft * F)); Ca_NSR' = -((i_leak + i_tr) - i_up); Ca_foot' = - i_CaCa * (A_cap / (2.0 * V_myo * F)) * R_A_V; // Variable initializations: Nai = 10.0; Nao = 140.0; Cao = 1.8; Ki = 145.0; Ko = 5.4; Ca_NSR = 15.0; 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_K_ATP = ; i_Ks = ; i_K1 = ; i_Kp = ; i_ns_K = ; i_tr = ; i_leak = ; i_up = ; end model *shaw_1997____main() // Sub-modules, and any changes to those submodules: environment: shaw_1997__environment(time_); membrane: shaw_1997__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, i_K_ATP); fast_sodium_current: shaw_1997__fast_sodium_current(i_Na, E_Na, time_, V, R, F, T, Nai, Nao, m, h, j); L_type_Ca_channel: shaw_1997__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: shaw_1997__T_type_Ca_channel(i_Ca_T, time_, V, Cai, Cao, R, T, F, b, g); rapid_time_dependent_potassium_current: shaw_1997__rapid_time_dependent_potassium_current(i_Kr, time_, V, R, T, F, Ko, Ki, Xr, Rr); slow_time_dependent_potassium_current: shaw_1997__slow_time_dependent_potassium_current(i_Ks, time_, V, Ko, Ki, Nao, Nai, Cai, R, T, F, Xs); time_independent_potassium_current: shaw_1997__time_independent_potassium_current(i_K1, E_K1, time_, V, Ko, Ki, R, T, F, K1_infinity); plateau_potassium_current: shaw_1997__plateau_potassium_current(i_Kp, time_, V, E_K1); sarcolemmal_calcium_pump: shaw_1997__sarcolemmal_calcium_pump(i_p_Ca, time_, Cai); sodium_background_current: shaw_1997__sodium_background_current(i_Na_b, time_, V, E_Na); calcium_background_current: shaw_1997__calcium_background_current(i_Ca_b, time_, V, R, T, F, Cai, Cao); sodium_potassium_pump: shaw_1997__sodium_potassium_pump(i_NaK, time_, V, R, T, F, Nai, Nao, Ko); non_specific_calcium_activated_current: shaw_1997__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); ATP_dependent_potassium_current: shaw_1997__ATP_dependent_potassium_current(i_K_ATP, time_, V, Ko, Ki, R, T, F); Na_Ca_exchanger: shaw_1997__Na_Ca_exchanger(i_NaCa, time_, V, R, T, F, Nai, Nao, Cai, Cao); calcium_buffers_in_the_myoplasm: shaw_1997__calcium_buffers_in_the_myoplasm(Ca_JSR_new, Cai, time_, CSQN_max, K_mCSQN, CSQN_buff); calcium_fluxes_in_the_SR: shaw_1997__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: shaw_1997__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_K_ATP, i_Ks, i_K1, i_Kp, i_ns_K, i_tr, i_leak, i_up); end