Chaste: heart/src/odes/ionicmodels/NobleVargheseKohlNoble1998WithSac.cpp Source File

00001 
00013 #include "NobleVargheseKohlNoble1998WithSac.hpp"
00014 #include "HeartConfig.hpp"
00015 
00016 CML_noble_varghese_kohl_noble_1998_basic_with_sac::CML_noble_varghese_kohl_noble_1998_basic_with_sac(
00017         boost::shared_ptr<AbstractIvpOdeSolver> pSolver,
00018         boost::shared_ptr<AbstractStimulusFunction> pIntracellularStimulus)
00019     : AbstractCardiacCell(pSolver, 22, 0, pIntracellularStimulus)
00020 {
00021     mStretch = 1.0;
00022     
00023     mpSystemInfo = OdeSystemInformation<CML_noble_varghese_kohl_noble_1998_basic_with_sac>::Instance();
00024 
00025     Init();
00026 
00027 }
00028 
00029 CML_noble_varghese_kohl_noble_1998_basic_with_sac::~CML_noble_varghese_kohl_noble_1998_basic_with_sac(void)
00030 {
00031 }
00032 
00033 double CML_noble_varghese_kohl_noble_1998_basic_with_sac::GetIIonic()
00034 {
00035     std::vector<double>& rY = rGetStateVariables();
00036     double var_membrane__V = rY[0];
00037     // Units: millivolt; Initial value: -92.849333
00038     double var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1 = rY[1];
00039     // Units: dimensionless; Initial value: 1.03e-5
00040     double var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2 = rY[2];
00041     // Units: dimensionless; Initial value: 2e-7
00042     double var_slow_delayed_rectifier_potassium_current_xs_gate__xs = rY[3];
00043     // Units: dimensionless; Initial value: 0.001302
00044     double var_fast_sodium_current_m_gate__m = rY[4];
00045     // Units: dimensionless; Initial value: 0.0016203
00046     double var_fast_sodium_current_h_gate__h = rY[5];
00047     // Units: dimensionless; Initial value: 0.9944036
00048     double var_L_type_Ca_channel_d_gate__d = rY[6];
00049     // Units: dimensionless; Initial value: 0
00050     double var_L_type_Ca_channel_f_gate__f = rY[7];
00051     // Units: dimensionless; Initial value: 1
00052     double var_L_type_Ca_channel_f2_gate__f2 = rY[8];
00053     // Units: dimensionless; Initial value: 0.9349197
00054     double var_L_type_Ca_channel_f2ds_gate__f2ds = rY[9];
00055     // Units: dimensionless; Initial value: 0.9651958
00056     double var_transient_outward_current_s_gate__s = rY[10];
00057     // Units: dimensionless; Initial value: 0.9948645
00058     double var_transient_outward_current_r_gate__r = rY[11];
00059     // Units: dimensionless; Initial value: 0
00060     double var_intracellular_sodium_concentration__Na_i = rY[14];
00061     // Units: millimolar; Initial value: 7.3321223
00062     double var_intracellular_potassium_concentration__K_i = rY[15];
00063     // Units: millimolar; Initial value: 136.5644281
00064     double var_intracellular_calcium_concentration__Ca_i = rY[16];
00065     // Units: millimolar; Initial value: 1.4e-5
00066     double var_intracellular_calcium_concentration__Ca_ds = rY[17];
00067     // Units: millimolar; Initial value: 1.88e-5
00068 
00069     const double var_membrane__R = 8314.472;
00070     const double var_membrane__T = 310.0;
00071     const double var_membrane__F = 96485.3415;
00072     double var_reversal_potentials__K_i = var_intracellular_potassium_concentration__K_i;
00073     double var_reversal_potentials__R = var_membrane__R;
00074     double var_reversal_potentials__T = var_membrane__T;
00075     double var_reversal_potentials__F = var_membrane__F;
00076     const double var_extracellular_potassium_concentration__K_o = 4.0;
00077     double var_reversal_potentials__K_o = var_extracellular_potassium_concentration__K_o;
00078     double var_reversal_potentials__E_K = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log(var_reversal_potentials__K_o / var_reversal_potentials__K_i);
00079     double var_time_independent_potassium_current__E_K = var_reversal_potentials__E_K;
00080     double var_time_independent_potassium_current__K_o = var_extracellular_potassium_concentration__K_o;
00081     double var_time_independent_potassium_current__R = var_membrane__R;
00082     double var_time_independent_potassium_current__V = var_membrane__V;
00083     double var_time_independent_potassium_current__T = var_membrane__T;
00084     const double var_time_independent_potassium_current__K_mk1 = 10.0;
00085     const double var_time_independent_potassium_current__g_K1 = 0.5;
00086     double var_time_independent_potassium_current__F = var_membrane__F;
00087     double var_time_independent_potassium_current__i_K1 = (((var_time_independent_potassium_current__g_K1 * var_time_independent_potassium_current__K_o) / (var_time_independent_potassium_current__K_o + var_time_independent_potassium_current__K_mk1)) * (var_time_independent_potassium_current__V - var_time_independent_potassium_current__E_K)) / (1.0 + exp((((var_time_independent_potassium_current__V - var_time_independent_potassium_current__E_K) - 10.0) * var_time_independent_potassium_current__F * 1.25) / (var_time_independent_potassium_current__R * var_time_independent_potassium_current__T)));
00088     double var_membrane__i_K1 = var_time_independent_potassium_current__i_K1;
00089     double var_transient_outward_current__s = var_transient_outward_current_s_gate__s;
00090     double var_transient_outward_current__r = var_transient_outward_current_r_gate__r;
00091     const double var_transient_outward_current__g_to = 0.005;
00092     double var_transient_outward_current__V = var_membrane__V;
00093     double var_transient_outward_current__E_K = var_reversal_potentials__E_K;
00094     const double var_transient_outward_current__g_tos = 0.0;
00095     double var_transient_outward_current__i_to = var_transient_outward_current__g_to * (var_transient_outward_current__g_tos + (var_transient_outward_current__s * (1.0 - var_transient_outward_current__g_tos))) * var_transient_outward_current__r * (var_transient_outward_current__V - var_transient_outward_current__E_K);
00096     double var_membrane__i_to = var_transient_outward_current__i_to;
00097     const double var_rapid_delayed_rectifier_potassium_current__g_Kr2 = 0.0013;
00098     const double var_rapid_delayed_rectifier_potassium_current__g_Kr1 = 0.0021;
00099     double var_rapid_delayed_rectifier_potassium_current__xr1 = var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1;
00100     double var_rapid_delayed_rectifier_potassium_current__xr2 = var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2;
00101     double var_rapid_delayed_rectifier_potassium_current__V = var_membrane__V;
00102     double var_rapid_delayed_rectifier_potassium_current__E_K = var_reversal_potentials__E_K;
00103     double var_rapid_delayed_rectifier_potassium_current__i_Kr = ((((var_rapid_delayed_rectifier_potassium_current__g_Kr1 * var_rapid_delayed_rectifier_potassium_current__xr1) + (var_rapid_delayed_rectifier_potassium_current__g_Kr2 * var_rapid_delayed_rectifier_potassium_current__xr2)) * 1.0) / (1.0 + exp((var_rapid_delayed_rectifier_potassium_current__V + 9.0) / 22.4))) * (var_rapid_delayed_rectifier_potassium_current__V - var_rapid_delayed_rectifier_potassium_current__E_K);
00104     double var_membrane__i_Kr = var_rapid_delayed_rectifier_potassium_current__i_Kr;
00105     double var_slow_delayed_rectifier_potassium_current__xs = var_slow_delayed_rectifier_potassium_current_xs_gate__xs;
00106     const double var_extracellular_sodium_concentration__Na_o = 140.0;
00107     double var_reversal_potentials__Na_o = var_extracellular_sodium_concentration__Na_o;
00108     double var_reversal_potentials__Na_i = var_intracellular_sodium_concentration__Na_i;
00109     const double var_reversal_potentials__P_kna = 0.03;
00110     double var_reversal_potentials__E_Ks = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log((var_reversal_potentials__K_o + (var_reversal_potentials__P_kna * var_reversal_potentials__Na_o)) / (var_reversal_potentials__K_i + (var_reversal_potentials__P_kna * var_reversal_potentials__Na_i)));
00111     double var_slow_delayed_rectifier_potassium_current__E_Ks = var_reversal_potentials__E_Ks;
00112     const double var_slow_delayed_rectifier_potassium_current__g_Ks = 0.0026;
00113     double var_slow_delayed_rectifier_potassium_current__V = var_membrane__V;
00114     double var_slow_delayed_rectifier_potassium_current__i_Ks = var_slow_delayed_rectifier_potassium_current__g_Ks * pow(var_slow_delayed_rectifier_potassium_current__xs, 2.0) * (var_slow_delayed_rectifier_potassium_current__V - var_slow_delayed_rectifier_potassium_current__E_Ks);
00115     double var_membrane__i_Ks = var_slow_delayed_rectifier_potassium_current__i_Ks;
00116     double var_L_type_Ca_channel__d = var_L_type_Ca_channel_d_gate__d;
00117     const double var_L_type_Ca_channel__FrICa = 1.0;
00118     double var_L_type_Ca_channel__f = var_L_type_Ca_channel_f_gate__f;
00119     double var_L_type_Ca_channel__K_o = var_extracellular_potassium_concentration__K_o;
00120     double var_L_type_Ca_channel__K_i = var_intracellular_potassium_concentration__K_i;
00121     double var_L_type_Ca_channel__F = var_membrane__F;
00122     const double var_L_type_Ca_channel__P_Ca_L = 0.1;
00123     double var_L_type_Ca_channel__T = var_membrane__T;
00124     const double var_L_type_Ca_channel__P_CaK = 0.002;
00125     double var_L_type_Ca_channel__V = var_membrane__V;
00126     double var_L_type_Ca_channel__f2 = var_L_type_Ca_channel_f2_gate__f2;
00127     double var_L_type_Ca_channel__R = var_membrane__R;
00128     double var_L_type_Ca_channel__i_Ca_L_K_cyt = ((((1.0 - var_L_type_Ca_channel__FrICa) * var_L_type_Ca_channel__P_CaK * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2 * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__K_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__K_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00129     double var_membrane__i_Ca_L_K_cyt = var_L_type_Ca_channel__i_Ca_L_K_cyt;
00130     double var_L_type_Ca_channel__f2ds = var_L_type_Ca_channel_f2ds_gate__f2ds;
00131     double var_L_type_Ca_channel__i_Ca_L_K_ds = (((var_L_type_Ca_channel__FrICa * var_L_type_Ca_channel__P_CaK * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2ds * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__K_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__K_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00132     double var_membrane__i_Ca_L_K_ds = var_L_type_Ca_channel__i_Ca_L_K_ds;
00133     const double var_sodium_potassium_pump__i_NaK_max = 0.7;
00134     double var_sodium_potassium_pump__Na_i = var_intracellular_sodium_concentration__Na_i;
00135     double var_sodium_potassium_pump__K_o = var_extracellular_potassium_concentration__K_o;
00136     const double var_sodium_potassium_pump__K_mNa = 40.0;
00137     const double var_sodium_potassium_pump__K_mK = 1.0;
00138     double var_sodium_potassium_pump__i_NaK = (((var_sodium_potassium_pump__i_NaK_max * var_sodium_potassium_pump__K_o) / (var_sodium_potassium_pump__K_mK + var_sodium_potassium_pump__K_o)) * var_sodium_potassium_pump__Na_i) / (var_sodium_potassium_pump__K_mNa + var_sodium_potassium_pump__Na_i);
00139     double var_membrane__i_NaK = var_sodium_potassium_pump__i_NaK;
00140     const double var_fast_sodium_current__g_Na = 2.5;
00141     double var_fast_sodium_current__h = var_fast_sodium_current_h_gate__h;
00142     double var_fast_sodium_current__V = var_membrane__V;
00143     double var_reversal_potentials__E_mh = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log((var_reversal_potentials__Na_o + (0.12 * var_reversal_potentials__K_o)) / (var_reversal_potentials__Na_i + (0.12 * var_reversal_potentials__K_i)));
00144     double var_fast_sodium_current__E_mh = var_reversal_potentials__E_mh;
00145     double var_fast_sodium_current__m = var_fast_sodium_current_m_gate__m;
00146     double var_fast_sodium_current__i_Na = var_fast_sodium_current__g_Na * pow(var_fast_sodium_current__m, 3.0) * var_fast_sodium_current__h * (var_fast_sodium_current__V - var_fast_sodium_current__E_mh);
00147     double var_membrane__i_Na = var_fast_sodium_current__i_Na;
00148     double var_sodium_background_current__V = var_membrane__V;
00149     double var_reversal_potentials__E_Na = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log(var_reversal_potentials__Na_o / var_reversal_potentials__Na_i);
00150     double var_sodium_background_current__E_Na = var_reversal_potentials__E_Na;
00151     const double var_sodium_background_current__g_bna = 0.0006;
00152     double var_sodium_background_current__i_b_Na = var_sodium_background_current__g_bna * (var_sodium_background_current__V - var_sodium_background_current__E_Na);
00153     double var_membrane__i_b_Na = var_sodium_background_current__i_b_Na;
00154     const double var_persistent_sodium_current__g_pna = 0.004;
00155     double var_persistent_sodium_current__V = var_membrane__V;
00156     double var_persistent_sodium_current__E_Na = var_reversal_potentials__E_Na;
00157     double var_persistent_sodium_current__i_p_Na = ((var_persistent_sodium_current__g_pna * 1.0) / (1.0 + exp((-(var_persistent_sodium_current__V + 52.0)) / 8.0))) * (var_persistent_sodium_current__V - var_persistent_sodium_current__E_Na);
00158     double var_membrane__i_p_Na = var_persistent_sodium_current__i_p_Na;
00159     const double var_L_type_Ca_channel__P_CaNa = 0.01;
00160     double var_L_type_Ca_channel__Na_o = var_extracellular_sodium_concentration__Na_o;
00161     double var_L_type_Ca_channel__Na_i = var_intracellular_sodium_concentration__Na_i;
00162     double var_L_type_Ca_channel__i_Ca_L_Na_cyt = ((((1.0 - var_L_type_Ca_channel__FrICa) * var_L_type_Ca_channel__P_CaNa * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2 * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Na_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Na_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00163     double var_membrane__i_Ca_L_Na_cyt = var_L_type_Ca_channel__i_Ca_L_Na_cyt;
00164     double var_L_type_Ca_channel__i_Ca_L_Na_ds = (((var_L_type_Ca_channel__FrICa * var_L_type_Ca_channel__P_CaNa * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2ds * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Na_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Na_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00165     double var_membrane__i_Ca_L_Na_ds = var_L_type_Ca_channel__i_Ca_L_Na_ds;
00166     double var_sodium_calcium_exchanger__Na_i = var_intracellular_sodium_concentration__Na_i;
00167     const double var_sodium_calcium_exchanger__n_NaCa = 3.0;
00168     const double var_sodium_calcium_exchanger__gamma = 0.5;
00169     double var_sodium_calcium_exchanger__F = var_membrane__F;
00170     double var_sodium_calcium_exchanger__Na_o = var_extracellular_sodium_concentration__Na_o;
00171     const double var_sodium_calcium_exchanger__FRiNaCa = 0.001;
00172     double var_sodium_calcium_exchanger__R = var_membrane__R;
00173     double var_sodium_calcium_exchanger__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00174     double var_sodium_calcium_exchanger__T = var_membrane__T;
00175     double var_sodium_calcium_exchanger__V = var_membrane__V;
00176     const double var_sodium_calcium_exchanger__d_NaCa = 0.0;
00177     const double var_extracellular_calcium_concentration__Ca_o = 2.0;
00178     double var_sodium_calcium_exchanger__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00179     const double var_sodium_calcium_exchanger__k_NaCa = 0.0005;
00180     double var_sodium_calcium_exchanger__i_NaCa_cyt = ((1.0 - var_sodium_calcium_exchanger__FRiNaCa) * var_sodium_calcium_exchanger__k_NaCa * ((exp((var_sodium_calcium_exchanger__gamma * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_o) - (exp(((var_sodium_calcium_exchanger__gamma - 1.0) * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_i))) / ((1.0 + (var_sodium_calcium_exchanger__d_NaCa * ((var_sodium_calcium_exchanger__Ca_i * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa)) + (var_sodium_calcium_exchanger__Ca_o * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa))))) * (1.0 + (var_sodium_calcium_exchanger__Ca_i / 0.0069)));
00181     double var_membrane__i_NaCa_cyt = var_sodium_calcium_exchanger__i_NaCa_cyt;
00182     double var_sodium_calcium_exchanger__Ca_ds = var_intracellular_calcium_concentration__Ca_ds;
00183     double var_sodium_calcium_exchanger__i_NaCa_ds = (var_sodium_calcium_exchanger__FRiNaCa * var_sodium_calcium_exchanger__k_NaCa * ((exp((var_sodium_calcium_exchanger__gamma * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_o) - (exp(((var_sodium_calcium_exchanger__gamma - 1.0) * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_ds))) / ((1.0 + (var_sodium_calcium_exchanger__d_NaCa * ((var_sodium_calcium_exchanger__Ca_ds * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa)) + (var_sodium_calcium_exchanger__Ca_o * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa))))) * (1.0 + (var_sodium_calcium_exchanger__Ca_ds / 0.0069)));
00184     double var_membrane__i_NaCa_ds = var_sodium_calcium_exchanger__i_NaCa_ds;
00185     double var_L_type_Ca_channel__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00186     double var_L_type_Ca_channel__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00187     double var_L_type_Ca_channel__i_Ca_L_Ca_cyt = ((((1.0 - var_L_type_Ca_channel__FrICa) * 4.0 * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2 * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Ca_i * exp((100.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Ca_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00188     double var_membrane__i_Ca_L_Ca_cyt = var_L_type_Ca_channel__i_Ca_L_Ca_cyt;
00189     double var_L_type_Ca_channel__i_Ca_L_Ca_ds = (((var_L_type_Ca_channel__FrICa * 4.0 * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2ds * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Ca_i * exp((100.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Ca_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00190     double var_membrane__i_Ca_L_Ca_ds = var_L_type_Ca_channel__i_Ca_L_Ca_ds;
00191     double var_reversal_potentials__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00192     double var_reversal_potentials__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00193     double var_reversal_potentials__E_Ca = ((0.5 * var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log(var_reversal_potentials__Ca_o / var_reversal_potentials__Ca_i);
00194     double var_calcium_background_current__E_Ca = var_reversal_potentials__E_Ca;
00195     const double var_calcium_background_current__g_bca = 0.00025;
00196     double var_calcium_background_current__V = var_membrane__V;
00197     double var_calcium_background_current__i_b_Ca = var_calcium_background_current__g_bca * (var_calcium_background_current__V - var_calcium_background_current__E_Ca);
00198     double var_membrane__i_b_Ca = var_calcium_background_current__i_b_Ca;
00199 
00201     // new part of the model - addition of a Stretch-activated channel
00203     const double g_sac = 0.035; // uS   
00204     const double E_sac = -10; // mV
00205     double f = (mStretch > 1.0) ? (mStretch-1.0)/0.15 : 0.0; // f = 0 if stretch < 1, scales linearly to f=1 at 15% stretch 
00206     double sac_ionic_current = g_sac * f * (var_membrane__V - E_sac); // if g is uS, this is nA
00207 
00218     double value_in_nA = var_membrane__i_K1+var_membrane__i_to+var_membrane__i_Kr+var_membrane__i_Ks+var_membrane__i_Ca_L_K_cyt+var_membrane__i_Ca_L_K_ds+var_membrane__i_NaK+var_membrane__i_Na+var_membrane__i_b_Na+var_membrane__i_p_Na+var_membrane__i_Ca_L_Na_cyt+var_membrane__i_Ca_L_Na_ds+var_membrane__i_NaCa_cyt+var_membrane__i_NaCa_ds+var_membrane__i_Ca_L_Ca_cyt+var_membrane__i_Ca_L_Ca_ds+var_membrane__i_b_Ca;
00219 
00220 //    std::cout << value_in_nA << " " << sac_ionic_current << "\n";
00221     value_in_nA += sac_ionic_current;
00222 
00223     double value_in_microA = 0.001*value_in_nA;
00224     double estimated_cell_surface_in_cm_square = 9.5e-05 / HeartConfig::Instance()->GetCapacitance();
00225     double value_in_microA_per_cm_square = value_in_microA/estimated_cell_surface_in_cm_square;
00226     return value_in_microA_per_cm_square;
00227 }
00228 
00229 void CML_noble_varghese_kohl_noble_1998_basic_with_sac::EvaluateYDerivatives (
00230         double var_environment__time,
00231         const std::vector<double> &rY,
00232         std::vector<double> &rDY)
00233 {
00234     // Inputs:
00235     // Time units: second
00236     var_environment__time *= 0.001;
00237     double var_membrane__V = rY[0];
00238     // Units: millivolt; Initial value: -92.849333
00239     double var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1 = rY[1];
00240     // Units: dimensionless; Initial value: 1.03e-5
00241     double var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2 = rY[2];
00242     // Units: dimensionless; Initial value: 2e-7
00243     double var_slow_delayed_rectifier_potassium_current_xs_gate__xs = rY[3];
00244     // Units: dimensionless; Initial value: 0.001302
00245     double var_fast_sodium_current_m_gate__m = rY[4];
00246     // Units: dimensionless; Initial value: 0.0016203
00247     double var_fast_sodium_current_h_gate__h = rY[5];
00248     // Units: dimensionless; Initial value: 0.9944036
00249     double var_L_type_Ca_channel_d_gate__d = rY[6];
00250     // Units: dimensionless; Initial value: 0
00251     double var_L_type_Ca_channel_f_gate__f = rY[7];
00252     // Units: dimensionless; Initial value: 1
00253     double var_L_type_Ca_channel_f2_gate__f2 = rY[8];
00254     // Units: dimensionless; Initial value: 0.9349197
00255     double var_L_type_Ca_channel_f2ds_gate__f2ds = rY[9];
00256     // Units: dimensionless; Initial value: 0.9651958
00257     double var_transient_outward_current_s_gate__s = rY[10];
00258     // Units: dimensionless; Initial value: 0.9948645
00259     double var_transient_outward_current_r_gate__r = rY[11];
00260     // Units: dimensionless; Initial value: 0
00261     double var_calcium_release__ActFrac = rY[12];
00262     // Units: dimensionless; Initial value: 0.0042614
00263     double var_calcium_release__ProdFrac = rY[13];
00264     // Units: dimensionless; Initial value: 0.4068154
00265     double var_intracellular_sodium_concentration__Na_i = rY[14];
00266     // Units: millimolar; Initial value: 7.3321223
00267     double var_intracellular_potassium_concentration__K_i = rY[15];
00268     // Units: millimolar; Initial value: 136.5644281
00269     double var_intracellular_calcium_concentration__Ca_i = rY[16];
00270     // Units: millimolar; Initial value: 1.4e-5
00271     double var_intracellular_calcium_concentration__Ca_ds = rY[17];
00272     // Units: millimolar; Initial value: 1.88e-5
00273     double var_intracellular_calcium_concentration__Ca_up = rY[18];
00274     // Units: millimolar; Initial value: 0.4531889
00275     double var_intracellular_calcium_concentration__Ca_rel = rY[19];
00276     // Units: millimolar; Initial value: 0.4481927
00277     double var_intracellular_calcium_concentration__Ca_Calmod = rY[20];
00278     // Units: millimolar; Initial value: 0.0005555
00279     double var_intracellular_calcium_concentration__Ca_Trop = rY[21];
00280     // Units: millimolar; Initial value: 0.0003542
00281 
00282 
00283     // Mathematics
00284     const double var_membrane__R = 8314.472;
00285     const double var_membrane__T = 310.0;
00286     const double var_membrane__F = 96485.3415;
00287     const double var_membrane__Cm = 9.5e-05;
00288     double var_reversal_potentials__K_i = var_intracellular_potassium_concentration__K_i;
00289     double var_reversal_potentials__R = var_membrane__R;
00290     double var_reversal_potentials__T = var_membrane__T;
00291     double var_reversal_potentials__F = var_membrane__F;
00292     const double var_extracellular_potassium_concentration__K_o = 4.0;
00293     double var_reversal_potentials__K_o = var_extracellular_potassium_concentration__K_o;
00294     double var_reversal_potentials__E_K = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log(var_reversal_potentials__K_o / var_reversal_potentials__K_i);
00295     double var_time_independent_potassium_current__E_K = var_reversal_potentials__E_K;
00296     double var_time_independent_potassium_current__K_o = var_extracellular_potassium_concentration__K_o;
00297     double var_time_independent_potassium_current__R = var_membrane__R;
00298     double var_time_independent_potassium_current__V = var_membrane__V;
00299     double var_time_independent_potassium_current__T = var_membrane__T;
00300     const double var_time_independent_potassium_current__K_mk1 = 10.0;
00301     const double var_time_independent_potassium_current__g_K1 = 0.5;
00302     double var_time_independent_potassium_current__F = var_membrane__F;
00303     double var_time_independent_potassium_current__i_K1 = (((var_time_independent_potassium_current__g_K1 * var_time_independent_potassium_current__K_o) / (var_time_independent_potassium_current__K_o + var_time_independent_potassium_current__K_mk1)) * (var_time_independent_potassium_current__V - var_time_independent_potassium_current__E_K)) / (1.0 + exp((((var_time_independent_potassium_current__V - var_time_independent_potassium_current__E_K) - 10.0) * var_time_independent_potassium_current__F * 1.25) / (var_time_independent_potassium_current__R * var_time_independent_potassium_current__T)));
00304     double var_membrane__i_K1 = var_time_independent_potassium_current__i_K1;
00305     double var_transient_outward_current__s = var_transient_outward_current_s_gate__s;
00306     double var_transient_outward_current__r = var_transient_outward_current_r_gate__r;
00307     const double var_transient_outward_current__g_to = 0.005;
00308     double var_transient_outward_current__V = var_membrane__V;
00309     double var_transient_outward_current__E_K = var_reversal_potentials__E_K;
00310     const double var_transient_outward_current__g_tos = 0.0;
00311     double var_transient_outward_current__i_to = var_transient_outward_current__g_to * (var_transient_outward_current__g_tos + (var_transient_outward_current__s * (1.0 - var_transient_outward_current__g_tos))) * var_transient_outward_current__r * (var_transient_outward_current__V - var_transient_outward_current__E_K);
00312     double var_membrane__i_to = var_transient_outward_current__i_to;
00313     const double var_rapid_delayed_rectifier_potassium_current__g_Kr2 = 0.0013;
00314     const double var_rapid_delayed_rectifier_potassium_current__g_Kr1 = 0.0021;
00315     double var_rapid_delayed_rectifier_potassium_current__xr1 = var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1;
00316     double var_rapid_delayed_rectifier_potassium_current__xr2 = var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2;
00317     double var_rapid_delayed_rectifier_potassium_current__V = var_membrane__V;
00318     double var_rapid_delayed_rectifier_potassium_current__E_K = var_reversal_potentials__E_K;
00319     double var_rapid_delayed_rectifier_potassium_current__i_Kr = ((((var_rapid_delayed_rectifier_potassium_current__g_Kr1 * var_rapid_delayed_rectifier_potassium_current__xr1) + (var_rapid_delayed_rectifier_potassium_current__g_Kr2 * var_rapid_delayed_rectifier_potassium_current__xr2)) * 1.0) / (1.0 + exp((var_rapid_delayed_rectifier_potassium_current__V + 9.0) / 22.4))) * (var_rapid_delayed_rectifier_potassium_current__V - var_rapid_delayed_rectifier_potassium_current__E_K);
00320     double var_membrane__i_Kr = var_rapid_delayed_rectifier_potassium_current__i_Kr;
00321     double var_slow_delayed_rectifier_potassium_current__xs = var_slow_delayed_rectifier_potassium_current_xs_gate__xs;
00322     const double var_extracellular_sodium_concentration__Na_o = 140.0;
00323     double var_reversal_potentials__Na_o = var_extracellular_sodium_concentration__Na_o;
00324     double var_reversal_potentials__Na_i = var_intracellular_sodium_concentration__Na_i;
00325     const double var_reversal_potentials__P_kna = 0.03;
00326     double var_reversal_potentials__E_Ks = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log((var_reversal_potentials__K_o + (var_reversal_potentials__P_kna * var_reversal_potentials__Na_o)) / (var_reversal_potentials__K_i + (var_reversal_potentials__P_kna * var_reversal_potentials__Na_i)));
00327     double var_slow_delayed_rectifier_potassium_current__E_Ks = var_reversal_potentials__E_Ks;
00328     const double var_slow_delayed_rectifier_potassium_current__g_Ks = 0.0026;
00329     double var_slow_delayed_rectifier_potassium_current__V = var_membrane__V;
00330     double var_slow_delayed_rectifier_potassium_current__i_Ks = var_slow_delayed_rectifier_potassium_current__g_Ks * pow(var_slow_delayed_rectifier_potassium_current__xs, 2.0) * (var_slow_delayed_rectifier_potassium_current__V - var_slow_delayed_rectifier_potassium_current__E_Ks);
00331     double var_membrane__i_Ks = var_slow_delayed_rectifier_potassium_current__i_Ks;
00332     double var_L_type_Ca_channel__d = var_L_type_Ca_channel_d_gate__d;
00333     const double var_L_type_Ca_channel__FrICa = 1.0;
00334     double var_L_type_Ca_channel__f = var_L_type_Ca_channel_f_gate__f;
00335     double var_L_type_Ca_channel__K_o = var_extracellular_potassium_concentration__K_o;
00336     double var_L_type_Ca_channel__K_i = var_intracellular_potassium_concentration__K_i;
00337     double var_L_type_Ca_channel__F = var_membrane__F;
00338     const double var_L_type_Ca_channel__P_Ca_L = 0.1;
00339     double var_L_type_Ca_channel__T = var_membrane__T;
00340     const double var_L_type_Ca_channel__P_CaK = 0.002;
00341     double var_L_type_Ca_channel__V = var_membrane__V;
00342     double var_L_type_Ca_channel__f2 = var_L_type_Ca_channel_f2_gate__f2;
00343     double var_L_type_Ca_channel__R = var_membrane__R;
00344     double var_L_type_Ca_channel__i_Ca_L_K_cyt = ((((1.0 - var_L_type_Ca_channel__FrICa) * var_L_type_Ca_channel__P_CaK * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2 * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__K_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__K_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00345     double var_membrane__i_Ca_L_K_cyt = var_L_type_Ca_channel__i_Ca_L_K_cyt;
00346     double var_L_type_Ca_channel__f2ds = var_L_type_Ca_channel_f2ds_gate__f2ds;
00347     double var_L_type_Ca_channel__i_Ca_L_K_ds = (((var_L_type_Ca_channel__FrICa * var_L_type_Ca_channel__P_CaK * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2ds * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__K_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__K_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00348     double var_membrane__i_Ca_L_K_ds = var_L_type_Ca_channel__i_Ca_L_K_ds;
00349     const double var_sodium_potassium_pump__i_NaK_max = 0.7;
00350     double var_sodium_potassium_pump__Na_i = var_intracellular_sodium_concentration__Na_i;
00351     double var_sodium_potassium_pump__K_o = var_extracellular_potassium_concentration__K_o;
00352     const double var_sodium_potassium_pump__K_mNa = 40.0;
00353     const double var_sodium_potassium_pump__K_mK = 1.0;
00354     double var_sodium_potassium_pump__i_NaK = (((var_sodium_potassium_pump__i_NaK_max * var_sodium_potassium_pump__K_o) / (var_sodium_potassium_pump__K_mK + var_sodium_potassium_pump__K_o)) * var_sodium_potassium_pump__Na_i) / (var_sodium_potassium_pump__K_mNa + var_sodium_potassium_pump__Na_i);
00355     double var_membrane__i_NaK = var_sodium_potassium_pump__i_NaK;
00356     const double var_fast_sodium_current__g_Na = 2.5;
00357     double var_fast_sodium_current__h = var_fast_sodium_current_h_gate__h;
00358     double var_fast_sodium_current__V = var_membrane__V;
00359     double var_reversal_potentials__E_mh = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log((var_reversal_potentials__Na_o + (0.12 * var_reversal_potentials__K_o)) / (var_reversal_potentials__Na_i + (0.12 * var_reversal_potentials__K_i)));
00360     double var_fast_sodium_current__E_mh = var_reversal_potentials__E_mh;
00361     double var_fast_sodium_current__m = var_fast_sodium_current_m_gate__m;
00362     double var_fast_sodium_current__i_Na = var_fast_sodium_current__g_Na * pow(var_fast_sodium_current__m, 3.0) * var_fast_sodium_current__h * (var_fast_sodium_current__V - var_fast_sodium_current__E_mh);
00363     double var_membrane__i_Na = var_fast_sodium_current__i_Na;
00364     double var_sodium_background_current__V = var_membrane__V;
00365     double var_reversal_potentials__E_Na = ((var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log(var_reversal_potentials__Na_o / var_reversal_potentials__Na_i);
00366     double var_sodium_background_current__E_Na = var_reversal_potentials__E_Na;
00367     const double var_sodium_background_current__g_bna = 0.0006;
00368     double var_sodium_background_current__i_b_Na = var_sodium_background_current__g_bna * (var_sodium_background_current__V - var_sodium_background_current__E_Na);
00369     double var_membrane__i_b_Na = var_sodium_background_current__i_b_Na;
00370     const double var_persistent_sodium_current__g_pna = 0.004;
00371     double var_persistent_sodium_current__V = var_membrane__V;
00372     double var_persistent_sodium_current__E_Na = var_reversal_potentials__E_Na;
00373     double var_persistent_sodium_current__i_p_Na = ((var_persistent_sodium_current__g_pna * 1.0) / (1.0 + exp((-(var_persistent_sodium_current__V + 52.0)) / 8.0))) * (var_persistent_sodium_current__V - var_persistent_sodium_current__E_Na);
00374     double var_membrane__i_p_Na = var_persistent_sodium_current__i_p_Na;
00375     const double var_L_type_Ca_channel__P_CaNa = 0.01;
00376     double var_L_type_Ca_channel__Na_o = var_extracellular_sodium_concentration__Na_o;
00377     double var_L_type_Ca_channel__Na_i = var_intracellular_sodium_concentration__Na_i;
00378     double var_L_type_Ca_channel__i_Ca_L_Na_cyt = ((((1.0 - var_L_type_Ca_channel__FrICa) * var_L_type_Ca_channel__P_CaNa * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2 * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Na_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Na_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00379     double var_membrane__i_Ca_L_Na_cyt = var_L_type_Ca_channel__i_Ca_L_Na_cyt;
00380     double var_L_type_Ca_channel__i_Ca_L_Na_ds = (((var_L_type_Ca_channel__FrICa * var_L_type_Ca_channel__P_CaNa * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2ds * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Na_i * exp((50.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Na_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00381     double var_membrane__i_Ca_L_Na_ds = var_L_type_Ca_channel__i_Ca_L_Na_ds;
00382     double var_sodium_calcium_exchanger__Na_i = var_intracellular_sodium_concentration__Na_i;
00383     const double var_sodium_calcium_exchanger__n_NaCa = 3.0;
00384     const double var_sodium_calcium_exchanger__gamma = 0.5;
00385     double var_sodium_calcium_exchanger__F = var_membrane__F;
00386     double var_sodium_calcium_exchanger__Na_o = var_extracellular_sodium_concentration__Na_o;
00387     const double var_sodium_calcium_exchanger__FRiNaCa = 0.001;
00388     double var_sodium_calcium_exchanger__R = var_membrane__R;
00389     double var_sodium_calcium_exchanger__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00390     double var_sodium_calcium_exchanger__T = var_membrane__T;
00391     double var_sodium_calcium_exchanger__V = var_membrane__V;
00392     const double var_sodium_calcium_exchanger__d_NaCa = 0.0;
00393     const double var_extracellular_calcium_concentration__Ca_o = 2.0;
00394     double var_sodium_calcium_exchanger__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00395     const double var_sodium_calcium_exchanger__k_NaCa = 0.0005;
00396     double var_sodium_calcium_exchanger__i_NaCa_cyt = ((1.0 - var_sodium_calcium_exchanger__FRiNaCa) * var_sodium_calcium_exchanger__k_NaCa * ((exp((var_sodium_calcium_exchanger__gamma * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_o) - (exp(((var_sodium_calcium_exchanger__gamma - 1.0) * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_i))) / ((1.0 + (var_sodium_calcium_exchanger__d_NaCa * ((var_sodium_calcium_exchanger__Ca_i * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa)) + (var_sodium_calcium_exchanger__Ca_o * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa))))) * (1.0 + (var_sodium_calcium_exchanger__Ca_i / 0.0069)));
00397     double var_membrane__i_NaCa_cyt = var_sodium_calcium_exchanger__i_NaCa_cyt;
00398     double var_sodium_calcium_exchanger__Ca_ds = var_intracellular_calcium_concentration__Ca_ds;
00399     double var_sodium_calcium_exchanger__i_NaCa_ds = (var_sodium_calcium_exchanger__FRiNaCa * var_sodium_calcium_exchanger__k_NaCa * ((exp((var_sodium_calcium_exchanger__gamma * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_o) - (exp(((var_sodium_calcium_exchanger__gamma - 1.0) * (var_sodium_calcium_exchanger__n_NaCa - 2.0) * var_sodium_calcium_exchanger__V * var_sodium_calcium_exchanger__F) / (var_sodium_calcium_exchanger__R * var_sodium_calcium_exchanger__T)) * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa) * var_sodium_calcium_exchanger__Ca_ds))) / ((1.0 + (var_sodium_calcium_exchanger__d_NaCa * ((var_sodium_calcium_exchanger__Ca_ds * pow(var_sodium_calcium_exchanger__Na_o, var_sodium_calcium_exchanger__n_NaCa)) + (var_sodium_calcium_exchanger__Ca_o * pow(var_sodium_calcium_exchanger__Na_i, var_sodium_calcium_exchanger__n_NaCa))))) * (1.0 + (var_sodium_calcium_exchanger__Ca_ds / 0.0069)));
00400     double var_membrane__i_NaCa_ds = var_sodium_calcium_exchanger__i_NaCa_ds;
00401     double var_L_type_Ca_channel__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00402     double var_L_type_Ca_channel__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00403     double var_L_type_Ca_channel__i_Ca_L_Ca_cyt = ((((1.0 - var_L_type_Ca_channel__FrICa) * 4.0 * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2 * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Ca_i * exp((100.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Ca_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00404     double var_membrane__i_Ca_L_Ca_cyt = var_L_type_Ca_channel__i_Ca_L_Ca_cyt;
00405     double var_L_type_Ca_channel__i_Ca_L_Ca_ds = (((var_L_type_Ca_channel__FrICa * 4.0 * var_L_type_Ca_channel__P_Ca_L * var_L_type_Ca_channel__d * var_L_type_Ca_channel__f * var_L_type_Ca_channel__f2ds * (var_L_type_Ca_channel__V - 50.0) * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)) / (1.0 - exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T)))) * ((var_L_type_Ca_channel__Ca_i * exp((100.0 * var_L_type_Ca_channel__F) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))) - (var_L_type_Ca_channel__Ca_o * exp(((-(var_L_type_Ca_channel__V - 50.0)) * var_L_type_Ca_channel__F * 2.0) / (var_L_type_Ca_channel__R * var_L_type_Ca_channel__T))));
00406     double var_membrane__i_Ca_L_Ca_ds = var_L_type_Ca_channel__i_Ca_L_Ca_ds;
00407     double var_reversal_potentials__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00408     double var_reversal_potentials__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00409     double var_reversal_potentials__E_Ca = ((0.5 * var_reversal_potentials__R * var_reversal_potentials__T) / var_reversal_potentials__F) * log(var_reversal_potentials__Ca_o / var_reversal_potentials__Ca_i);
00410     double var_calcium_background_current__E_Ca = var_reversal_potentials__E_Ca;
00411     const double var_calcium_background_current__g_bca = 0.00025;
00412     double var_calcium_background_current__V = var_membrane__V;
00413     double var_calcium_background_current__i_b_Ca = var_calcium_background_current__g_bca * (var_calcium_background_current__V - var_calcium_background_current__E_Ca);
00414     double var_membrane__i_b_Ca = var_calcium_background_current__i_b_Ca;
00415     double var_membrane__i_Stim = GetStimulus((1.0/0.001)*var_environment__time);
00416     double var_rapid_delayed_rectifier_potassium_current_xr1_gate__V = var_rapid_delayed_rectifier_potassium_current__V;
00417     double var_rapid_delayed_rectifier_potassium_current_xr1_gate__alpha_xr1 = 50.0 / (1.0 + exp((-(var_rapid_delayed_rectifier_potassium_current_xr1_gate__V - 5.0)) / 9.0));
00418     double var_rapid_delayed_rectifier_potassium_current_xr1_gate__beta_xr1 = 0.05 * exp((-(var_rapid_delayed_rectifier_potassium_current_xr1_gate__V - 20.0)) / 15.0);
00419     double var_rapid_delayed_rectifier_potassium_current_xr2_gate__V = var_rapid_delayed_rectifier_potassium_current__V;
00420     double var_rapid_delayed_rectifier_potassium_current_xr2_gate__alpha_xr2 = 50.0 / (1.0 + exp((-(var_rapid_delayed_rectifier_potassium_current_xr2_gate__V - 5.0)) / 9.0));
00421     double var_rapid_delayed_rectifier_potassium_current_xr2_gate__beta_xr2 = 0.4 * exp(-pow((var_rapid_delayed_rectifier_potassium_current_xr2_gate__V + 30.0) / 30.0, 3.0));
00422     double var_slow_delayed_rectifier_potassium_current_xs_gate__V = var_slow_delayed_rectifier_potassium_current__V;
00423     double var_slow_delayed_rectifier_potassium_current_xs_gate__alpha_xs = 14.0 / (1.0 + exp((-(var_slow_delayed_rectifier_potassium_current_xs_gate__V - 40.0)) / 9.0));
00424     double var_slow_delayed_rectifier_potassium_current_xs_gate__beta_xs = 1.0 * exp((-var_slow_delayed_rectifier_potassium_current_xs_gate__V) / 45.0);
00425     double var_fast_sodium_current_m_gate__V = var_fast_sodium_current__V;
00426     double var_fast_sodium_current_m_gate__E0_m = var_fast_sodium_current_m_gate__V + 41.0;
00427     const double var_fast_sodium_current_m_gate__delta_m = 1e-05;
00428     double var_fast_sodium_current_m_gate__alpha_m = (fabs(var_fast_sodium_current_m_gate__E0_m) < var_fast_sodium_current_m_gate__delta_m) ? 2000.0 : ((200.0 * var_fast_sodium_current_m_gate__E0_m) / (1.0 - exp((-0.1) * var_fast_sodium_current_m_gate__E0_m)));
00429     double var_fast_sodium_current_m_gate__beta_m = 8000.0 * exp((-0.056) * (var_fast_sodium_current_m_gate__V + 66.0));
00430     double var_fast_sodium_current_h_gate__V = var_fast_sodium_current__V;
00431     const double var_fast_sodium_current_h_gate__shift_h = 0.0;
00432     double var_fast_sodium_current_h_gate__alpha_h = 20.0 * exp((-0.125) * ((var_fast_sodium_current_h_gate__V + 75.0) - var_fast_sodium_current_h_gate__shift_h));
00433     double var_fast_sodium_current_h_gate__beta_h = 2000.0 / (1.0 + (320.0 * exp((-0.1) * ((var_fast_sodium_current_h_gate__V + 75.0) - var_fast_sodium_current_h_gate__shift_h))));
00434     double var_L_type_Ca_channel__Ca_ds = var_intracellular_calcium_concentration__Ca_ds;
00435     const double var_L_type_Ca_channel__Km_f2 = 100000.0;
00436     const double var_L_type_Ca_channel__Km_f2ds = 0.001;
00437     const double var_L_type_Ca_channel__R_decay = 20.0;
00438     double var_L_type_Ca_channel_d_gate__V = var_L_type_Ca_channel__V;
00439     double var_L_type_Ca_channel_d_gate__E0_d = (var_L_type_Ca_channel_d_gate__V + 24.0) - 5.0;
00440     double var_L_type_Ca_channel_d_gate__alpha_d = (fabs(var_L_type_Ca_channel_d_gate__E0_d) < 0.0001) ? 120.0 : ((30.0 * var_L_type_Ca_channel_d_gate__E0_d) / (1.0 - exp((-var_L_type_Ca_channel_d_gate__E0_d) / 4.0)));
00441     double var_L_type_Ca_channel_d_gate__beta_d = (fabs(var_L_type_Ca_channel_d_gate__E0_d) < 0.0001) ? 120.0 : ((12.0 * var_L_type_Ca_channel_d_gate__E0_d) / (exp(var_L_type_Ca_channel_d_gate__E0_d / 10.0) - 1.0));
00442     const double var_L_type_Ca_channel_d_gate__speed_d = 3.0;
00443     double var_L_type_Ca_channel_f_gate__V = var_L_type_Ca_channel__V;
00444     double var_L_type_Ca_channel_f_gate__E0_f = var_L_type_Ca_channel_f_gate__V + 34.0;
00445     const double var_L_type_Ca_channel_f_gate__delta_f = 0.0001;
00446     double var_L_type_Ca_channel_f_gate__alpha_f = (fabs(var_L_type_Ca_channel_f_gate__E0_f) < var_L_type_Ca_channel_f_gate__delta_f) ? 25.0 : ((6.25 * var_L_type_Ca_channel_f_gate__E0_f) / (exp(var_L_type_Ca_channel_f_gate__E0_f / 4.0) - 1.0));
00447     double var_L_type_Ca_channel_f_gate__beta_f = 12.0 / (1.0 + exp(((-1.0) * (var_L_type_Ca_channel_f_gate__V + 34.0)) / 4.0));
00448     const double var_L_type_Ca_channel_f_gate__speed_f = 0.3;
00449     double var_L_type_Ca_channel_f2_gate__Km_f2 = var_L_type_Ca_channel__Km_f2;
00450     double var_L_type_Ca_channel_f2_gate__Ca_i = var_L_type_Ca_channel__Ca_i;
00451     double var_L_type_Ca_channel_f2ds_gate__Km_f2ds = var_L_type_Ca_channel__Km_f2ds;
00452     double var_L_type_Ca_channel_f2ds_gate__R_decay = var_L_type_Ca_channel__R_decay;
00453     double var_L_type_Ca_channel_f2ds_gate__Ca_ds = var_L_type_Ca_channel__Ca_ds;
00454     double var_transient_outward_current_s_gate__V = var_transient_outward_current__V;
00455     double var_transient_outward_current_s_gate__alpha_s = 0.033 * exp((-var_transient_outward_current_s_gate__V) / 17.0);
00456     double var_transient_outward_current_s_gate__beta_s = 33.0 / (1.0 + exp((-0.125) * (var_transient_outward_current_s_gate__V + 10.0)));
00457     double var_transient_outward_current_r_gate__V = var_transient_outward_current__V;
00458     double var_sarcoplasmic_reticulum_calcium_pump__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00459     double var_sarcoplasmic_reticulum_calcium_pump__Ca_up = var_intracellular_calcium_concentration__Ca_up;
00460     const double var_sarcoplasmic_reticulum_calcium_pump__alpha_up = 0.4;
00461     const double var_sarcoplasmic_reticulum_calcium_pump__beta_up = 0.03;
00462     const double var_sarcoplasmic_reticulum_calcium_pump__K_srca = 0.5;
00463     const double var_sarcoplasmic_reticulum_calcium_pump__K_xcs = 0.4;
00464     const double var_sarcoplasmic_reticulum_calcium_pump__K_cyca = 0.0003;
00465     double var_sarcoplasmic_reticulum_calcium_pump__K_1 = (var_sarcoplasmic_reticulum_calcium_pump__K_cyca * var_sarcoplasmic_reticulum_calcium_pump__K_xcs) / var_sarcoplasmic_reticulum_calcium_pump__K_srca;
00466     double var_sarcoplasmic_reticulum_calcium_pump__K_2 = var_sarcoplasmic_reticulum_calcium_pump__Ca_i + (var_sarcoplasmic_reticulum_calcium_pump__Ca_up * var_sarcoplasmic_reticulum_calcium_pump__K_1) + (var_sarcoplasmic_reticulum_calcium_pump__K_cyca * var_sarcoplasmic_reticulum_calcium_pump__K_xcs) + var_sarcoplasmic_reticulum_calcium_pump__K_cyca;
00467     double var_sarcoplasmic_reticulum_calcium_pump__i_up = ((var_sarcoplasmic_reticulum_calcium_pump__Ca_i / var_sarcoplasmic_reticulum_calcium_pump__K_2) * var_sarcoplasmic_reticulum_calcium_pump__alpha_up) - (((var_sarcoplasmic_reticulum_calcium_pump__Ca_up * var_sarcoplasmic_reticulum_calcium_pump__K_1) / var_sarcoplasmic_reticulum_calcium_pump__K_2) * var_sarcoplasmic_reticulum_calcium_pump__beta_up);
00468     double var_calcium_translocation__Ca_rel = var_intracellular_calcium_concentration__Ca_rel;
00469     double var_calcium_translocation__Ca_up = var_intracellular_calcium_concentration__Ca_up;
00470     double var_calcium_translocation__i_trans = 50.0 * (var_calcium_translocation__Ca_up - var_calcium_translocation__Ca_rel);
00471     const double var_calcium_release__K_m_rel = 250.0;
00472     const double var_calcium_release__K_leak_rate = 0.05;
00473     double var_calcium_release__Ca_rel = var_intracellular_calcium_concentration__Ca_rel;
00474     double var_calcium_release__i_rel = ((pow(var_calcium_release__ActFrac / (var_calcium_release__ActFrac + 0.25), 2.0) * var_calcium_release__K_m_rel) + var_calcium_release__K_leak_rate) * var_calcium_release__Ca_rel;
00475     double var_calcium_release__V = var_membrane__V;
00476     double var_calcium_release__VoltDep = exp(0.08 * (var_calcium_release__V - 40.0));
00477     const double var_calcium_release__K_m_Ca_cyt = 0.0005;
00478     double var_calcium_release__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00479     double var_calcium_release__CaiReg = var_calcium_release__Ca_i / (var_calcium_release__Ca_i + var_calcium_release__K_m_Ca_cyt);
00480     double var_calcium_release__Ca_ds = var_intracellular_calcium_concentration__Ca_ds;
00481     const double var_calcium_release__K_m_Ca_ds = 0.01;
00482     double var_calcium_release__CadsReg = var_calcium_release__Ca_ds / (var_calcium_release__Ca_ds + var_calcium_release__K_m_Ca_ds);
00483     double var_calcium_release__RegBindSite = var_calcium_release__CaiReg + ((1.0 - var_calcium_release__CaiReg) * var_calcium_release__CadsReg);
00484     double var_calcium_release__ActRate = (0.0 * var_calcium_release__VoltDep) + (500.0 * pow(var_calcium_release__RegBindSite, 2.0));
00485     double var_calcium_release__InactRate = 60.0 + (500.0 * pow(var_calcium_release__RegBindSite, 2.0));
00486     double var_calcium_release__PrecFrac = (1.0 - var_calcium_release__ActFrac) - var_calcium_release__ProdFrac;
00487     double var_calcium_release__SpeedRel = (var_calcium_release__V < (-50.0)) ? 5.0 : 1.0;
00488     const double var_intracellular_calcium_concentration__V_up_ratio = 0.01;
00489     const double var_intracellular_calcium_concentration__V_rel_ratio = 0.1;
00490     const double var_intracellular_calcium_concentration__V_e_ratio = 0.4;
00491     double var_intracellular_calcium_concentration__V_i_ratio = ((1.0 - var_intracellular_calcium_concentration__V_e_ratio) - var_intracellular_calcium_concentration__V_up_ratio) - var_intracellular_calcium_concentration__V_rel_ratio;
00492     const double var_intracellular_calcium_concentration__radius = 0.012;
00493     const double var_intracellular_calcium_concentration__length = 0.074;
00494     double var_intracellular_calcium_concentration__V_Cell = 3.141592654 * pow(var_intracellular_calcium_concentration__radius, 2.0) * var_intracellular_calcium_concentration__length;
00495     double var_intracellular_calcium_concentration__V_i = var_intracellular_calcium_concentration__V_Cell * var_intracellular_calcium_concentration__V_i_ratio;
00496     double var_intracellular_sodium_concentration__V_i = var_intracellular_calcium_concentration__V_i;
00497     double var_intracellular_sodium_concentration__F = var_membrane__F;
00498     double var_intracellular_sodium_concentration__i_Na = var_fast_sodium_current__i_Na;
00499     double var_intracellular_sodium_concentration__i_b_Na = var_sodium_background_current__i_b_Na;
00500     double var_intracellular_sodium_concentration__i_p_Na = var_persistent_sodium_current__i_p_Na;
00501     double var_intracellular_sodium_concentration__i_Ca_L_Na_cyt = var_L_type_Ca_channel__i_Ca_L_Na_cyt;
00502     double var_intracellular_sodium_concentration__i_Ca_L_Na_ds = var_L_type_Ca_channel__i_Ca_L_Na_ds;
00503     double var_intracellular_sodium_concentration__i_NaK = var_sodium_potassium_pump__i_NaK;
00504     double var_intracellular_sodium_concentration__i_NaCa_cyt = var_sodium_calcium_exchanger__i_NaCa_cyt;
00505     double var_intracellular_potassium_concentration__V_i = var_intracellular_calcium_concentration__V_i;
00506     double var_intracellular_potassium_concentration__i_K1 = var_time_independent_potassium_current__i_K1;
00507     double var_intracellular_potassium_concentration__i_Kr = var_rapid_delayed_rectifier_potassium_current__i_Kr;
00508     double var_intracellular_potassium_concentration__i_Ks = var_slow_delayed_rectifier_potassium_current__i_Ks;
00509     double var_intracellular_potassium_concentration__i_Ca_L_K_cyt = var_L_type_Ca_channel__i_Ca_L_K_cyt;
00510     double var_intracellular_potassium_concentration__i_Ca_L_K_ds = var_L_type_Ca_channel__i_Ca_L_K_ds;
00511     double var_intracellular_potassium_concentration__i_to = var_transient_outward_current__i_to;
00512     double var_intracellular_potassium_concentration__i_NaK = var_sodium_potassium_pump__i_NaK;
00513     double var_intracellular_potassium_concentration__F = var_membrane__F;
00514     const double var_intracellular_calcium_concentration__Calmod = 0.02;
00515     const double var_intracellular_calcium_concentration__Trop = 0.05;
00516     const double var_intracellular_calcium_concentration__alpha_Calmod = 100000.0;
00517     const double var_intracellular_calcium_concentration__beta_Calmod = 50.0;
00518     const double var_intracellular_calcium_concentration__alpha_Trop = 100000.0;
00519     const double var_intracellular_calcium_concentration__beta_Trop = 200.0;
00520     const double var_intracellular_calcium_concentration__V_ds_ratio = 0.1;
00521     const double var_intracellular_calcium_concentration__Kdecay = 10.0;
00522     double var_intracellular_calcium_concentration__i_up = var_sarcoplasmic_reticulum_calcium_pump__i_up;
00523     double var_intracellular_calcium_concentration__i_trans = var_calcium_translocation__i_trans;
00524     double var_intracellular_calcium_concentration__i_rel = var_calcium_release__i_rel;
00525     double var_intracellular_calcium_concentration__i_NaCa_cyt = var_sodium_calcium_exchanger__i_NaCa_cyt;
00526     double var_intracellular_calcium_concentration__i_NaCa_ds = var_sodium_calcium_exchanger__i_NaCa_ds;
00527     double var_intracellular_calcium_concentration__i_Ca_L_Ca_cyt = var_L_type_Ca_channel__i_Ca_L_Ca_cyt;
00528     double var_intracellular_calcium_concentration__i_Ca_L_Ca_ds = var_L_type_Ca_channel__i_Ca_L_Ca_ds;
00529     double var_intracellular_calcium_concentration__i_b_Ca = var_calcium_background_current__i_b_Ca;
00530     double var_intracellular_calcium_concentration__F = var_membrane__F;
00531 
00533     // new part of the model
00535     const double g_sac = 0.035; // uS  
00536     const double E_sac = -10; // mV
00537     double f = (mStretch > 0) ? (mStretch-1.0)/0.15 : 0.0; // f = 0 if stretch < 1, scales linearly to f=1 at 15% stretch 
00538     double sac_ionic_current = g_sac * f * (var_membrane__V - E_sac); // if g is uS, this is nA
00539 
00540     double d_dt_membrane__V;
00541     if (mSetVoltageDerivativeToZero)
00542     {
00543         d_dt_membrane__V = 0.0;
00544     }
00545     else
00546     {
00547         d_dt_membrane__V = ((-1.0) / var_membrane__Cm) * (sac_ionic_current + var_membrane__i_Stim + var_membrane__i_K1 + var_membrane__i_to + var_membrane__i_Kr + var_membrane__i_Ks + var_membrane__i_NaK + var_membrane__i_Na + var_membrane__i_b_Na + var_membrane__i_p_Na + var_membrane__i_Ca_L_Na_cyt + var_membrane__i_Ca_L_Na_ds + var_membrane__i_NaCa_cyt + var_membrane__i_NaCa_ds + var_membrane__i_Ca_L_Ca_cyt + var_membrane__i_Ca_L_Ca_ds + var_membrane__i_Ca_L_K_cyt + var_membrane__i_Ca_L_K_ds + var_membrane__i_b_Ca);
00548     }
00549     double d_dt_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1 = (var_rapid_delayed_rectifier_potassium_current_xr1_gate__alpha_xr1 * (1.0 - var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1)) - (var_rapid_delayed_rectifier_potassium_current_xr1_gate__beta_xr1 * var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1);
00550     double d_dt_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2 = (var_rapid_delayed_rectifier_potassium_current_xr2_gate__alpha_xr2 * (1.0 - var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2)) - (var_rapid_delayed_rectifier_potassium_current_xr2_gate__beta_xr2 * var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2);
00551     double d_dt_slow_delayed_rectifier_potassium_current_xs_gate__xs = (var_slow_delayed_rectifier_potassium_current_xs_gate__alpha_xs * (1.0 - var_slow_delayed_rectifier_potassium_current_xs_gate__xs)) - (var_slow_delayed_rectifier_potassium_current_xs_gate__beta_xs * var_slow_delayed_rectifier_potassium_current_xs_gate__xs);
00552     double d_dt_fast_sodium_current_m_gate__m = (var_fast_sodium_current_m_gate__alpha_m * (1.0 - var_fast_sodium_current_m_gate__m)) - (var_fast_sodium_current_m_gate__beta_m * var_fast_sodium_current_m_gate__m);
00553     double d_dt_fast_sodium_current_h_gate__h = (var_fast_sodium_current_h_gate__alpha_h * (1.0 - var_fast_sodium_current_h_gate__h)) - (var_fast_sodium_current_h_gate__beta_h * var_fast_sodium_current_h_gate__h);
00554     double d_dt_L_type_Ca_channel_d_gate__d = var_L_type_Ca_channel_d_gate__speed_d * ((var_L_type_Ca_channel_d_gate__alpha_d * (1.0 - var_L_type_Ca_channel_d_gate__d)) - (var_L_type_Ca_channel_d_gate__beta_d * var_L_type_Ca_channel_d_gate__d));
00555     double d_dt_L_type_Ca_channel_f_gate__f = var_L_type_Ca_channel_f_gate__speed_f * ((var_L_type_Ca_channel_f_gate__alpha_f * (1.0 - var_L_type_Ca_channel_f_gate__f)) - (var_L_type_Ca_channel_f_gate__beta_f * var_L_type_Ca_channel_f_gate__f));
00556     double d_dt_L_type_Ca_channel_f2_gate__f2 = 1.0 - (1.0 * ((var_L_type_Ca_channel_f2_gate__Ca_i / (var_L_type_Ca_channel_f2_gate__Km_f2 + var_L_type_Ca_channel_f2_gate__Ca_i)) + var_L_type_Ca_channel_f2_gate__f2));
00557     double d_dt_L_type_Ca_channel_f2ds_gate__f2ds = var_L_type_Ca_channel_f2ds_gate__R_decay * (1.0 - ((var_L_type_Ca_channel_f2ds_gate__Ca_ds / (var_L_type_Ca_channel_f2ds_gate__Km_f2ds + var_L_type_Ca_channel_f2ds_gate__Ca_ds)) + var_L_type_Ca_channel_f2ds_gate__f2ds));
00558     double d_dt_transient_outward_current_s_gate__s = (var_transient_outward_current_s_gate__alpha_s * (1.0 - var_transient_outward_current_s_gate__s)) - (var_transient_outward_current_s_gate__beta_s * var_transient_outward_current_s_gate__s);
00559     double d_dt_transient_outward_current_r_gate__r = 333.0 * ((1.0 / (1.0 + exp((-(var_transient_outward_current_r_gate__V + 4.0)) / 5.0))) - var_transient_outward_current_r_gate__r);
00560     double d_dt_calcium_release__ActFrac = (var_calcium_release__PrecFrac * var_calcium_release__SpeedRel * var_calcium_release__ActRate) - (var_calcium_release__ActFrac * var_calcium_release__SpeedRel * var_calcium_release__InactRate);
00561     double d_dt_calcium_release__ProdFrac = (var_calcium_release__ActFrac * var_calcium_release__SpeedRel * var_calcium_release__InactRate) - (var_calcium_release__SpeedRel * 1.0 * var_calcium_release__ProdFrac);
00562     double d_dt_intracellular_sodium_concentration__Na_i = ((-1.0) / (1.0 * var_intracellular_sodium_concentration__V_i * var_intracellular_sodium_concentration__F)) * (var_intracellular_sodium_concentration__i_Na + var_intracellular_sodium_concentration__i_p_Na + var_intracellular_sodium_concentration__i_b_Na + (3.0 * var_intracellular_sodium_concentration__i_NaK) + (3.0 * var_intracellular_sodium_concentration__i_NaCa_cyt) + var_intracellular_sodium_concentration__i_Ca_L_Na_cyt + var_intracellular_sodium_concentration__i_Ca_L_Na_ds);
00563     double d_dt_intracellular_potassium_concentration__K_i = ((-1.0) / (1.0 * var_intracellular_potassium_concentration__V_i * var_intracellular_potassium_concentration__F)) * ((var_intracellular_potassium_concentration__i_K1 + var_intracellular_potassium_concentration__i_Kr + var_intracellular_potassium_concentration__i_Ks + var_intracellular_potassium_concentration__i_Ca_L_K_cyt + var_intracellular_potassium_concentration__i_Ca_L_K_ds + var_intracellular_potassium_concentration__i_to) - (2.0 * var_intracellular_potassium_concentration__i_NaK));
00564     double d_dt_intracellular_calcium_concentration__Ca_Trop = (var_intracellular_calcium_concentration__alpha_Trop * var_intracellular_calcium_concentration__Ca_i * (var_intracellular_calcium_concentration__Trop - var_intracellular_calcium_concentration__Ca_Trop)) - (var_intracellular_calcium_concentration__beta_Trop * var_intracellular_calcium_concentration__Ca_Trop);
00565     double d_dt_intracellular_calcium_concentration__Ca_Calmod = (var_intracellular_calcium_concentration__alpha_Calmod * var_intracellular_calcium_concentration__Ca_i * (var_intracellular_calcium_concentration__Calmod - var_intracellular_calcium_concentration__Ca_Calmod)) - (var_intracellular_calcium_concentration__beta_Calmod * var_intracellular_calcium_concentration__Ca_Calmod);
00566     double d_dt_intracellular_calcium_concentration__Ca_i = ((((((-1.0) / (2.0 * 1.0 * var_intracellular_calcium_concentration__V_i * var_intracellular_calcium_concentration__F)) * (((var_intracellular_calcium_concentration__i_Ca_L_Ca_cyt + var_intracellular_calcium_concentration__i_b_Ca) - (2.0 * var_intracellular_calcium_concentration__i_NaCa_cyt)) - (2.0 * var_intracellular_calcium_concentration__i_NaCa_ds))) + (var_intracellular_calcium_concentration__Ca_ds * var_intracellular_calcium_concentration__V_ds_ratio * var_intracellular_calcium_concentration__Kdecay) + ((var_intracellular_calcium_concentration__i_rel * var_intracellular_calcium_concentration__V_rel_ratio) / var_intracellular_calcium_concentration__V_i_ratio)) - d_dt_intracellular_calcium_concentration__Ca_Calmod) - d_dt_intracellular_calcium_concentration__Ca_Trop) - var_intracellular_calcium_concentration__i_up;
00567     double d_dt_intracellular_calcium_concentration__Ca_ds = (((-1.0) * var_intracellular_calcium_concentration__i_Ca_L_Ca_ds) / (2.0 * 1.0 * var_intracellular_calcium_concentration__V_ds_ratio * var_intracellular_calcium_concentration__V_i * var_intracellular_calcium_concentration__F)) - (var_intracellular_calcium_concentration__Ca_ds * var_intracellular_calcium_concentration__Kdecay);
00568     double d_dt_intracellular_calcium_concentration__Ca_up = ((var_intracellular_calcium_concentration__V_i_ratio / var_intracellular_calcium_concentration__V_up_ratio) * var_intracellular_calcium_concentration__i_up) - var_intracellular_calcium_concentration__i_trans;
00569     double d_dt_intracellular_calcium_concentration__Ca_rel = ((var_intracellular_calcium_concentration__V_up_ratio / var_intracellular_calcium_concentration__V_rel_ratio) * var_intracellular_calcium_concentration__i_trans) - var_intracellular_calcium_concentration__i_rel;
00570 
00571     rDY[0] = 0.001*d_dt_membrane__V;
00572     rDY[1] = 0.001*d_dt_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1;
00573     rDY[2] = 0.001*d_dt_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2;
00574     rDY[3] = 0.001*d_dt_slow_delayed_rectifier_potassium_current_xs_gate__xs;
00575     rDY[4] = 0.001*d_dt_fast_sodium_current_m_gate__m;
00576     rDY[5] = 0.001*d_dt_fast_sodium_current_h_gate__h;
00577     rDY[6] = 0.001*d_dt_L_type_Ca_channel_d_gate__d;
00578     rDY[7] = 0.001*d_dt_L_type_Ca_channel_f_gate__f;
00579     rDY[8] = 0.001*d_dt_L_type_Ca_channel_f2_gate__f2;
00580     rDY[9] = 0.001*d_dt_L_type_Ca_channel_f2ds_gate__f2ds;
00581     rDY[10] = 0.001*d_dt_transient_outward_current_s_gate__s;
00582     rDY[11] = 0.001*d_dt_transient_outward_current_r_gate__r;
00583     rDY[12] = 0.001*d_dt_calcium_release__ActFrac;
00584     rDY[13] = 0.001*d_dt_calcium_release__ProdFrac;
00585     rDY[14] = 0.001*d_dt_intracellular_sodium_concentration__Na_i;
00586     rDY[15] = 0.001*d_dt_intracellular_potassium_concentration__K_i;
00587     rDY[16] = 0.001*d_dt_intracellular_calcium_concentration__Ca_i;
00588     rDY[17] = 0.001*d_dt_intracellular_calcium_concentration__Ca_ds;
00589     rDY[18] = 0.001*d_dt_intracellular_calcium_concentration__Ca_up;
00590     rDY[19] = 0.001*d_dt_intracellular_calcium_concentration__Ca_rel;
00591     rDY[20] = 0.001*d_dt_intracellular_calcium_concentration__Ca_Calmod;
00592     rDY[21] = 0.001*d_dt_intracellular_calcium_concentration__Ca_Trop;
00593 }
00594 
00595 
00596 template<>
00597 void OdeSystemInformation<CML_noble_varghese_kohl_noble_1998_basic_with_sac>::Initialise(void)
00598 {
00599     // Time units: second
00600     //
00601     this->mVariableNames.push_back("V");
00602     this->mVariableUnits.push_back("millivolt");
00603     this->mInitialConditions.push_back(-92.849333);
00604 
00605     this->mVariableNames.push_back("xr1");
00606     this->mVariableUnits.push_back("dimensionless");
00607     this->mInitialConditions.push_back(1.03e-5);
00608 
00609     this->mVariableNames.push_back("xr2");
00610     this->mVariableUnits.push_back("dimensionless");
00611     this->mInitialConditions.push_back(2e-7);
00612 
00613     this->mVariableNames.push_back("xs");
00614     this->mVariableUnits.push_back("dimensionless");
00615     this->mInitialConditions.push_back(0.001302);
00616 
00617     this->mVariableNames.push_back("m");
00618     this->mVariableUnits.push_back("dimensionless");
00619     this->mInitialConditions.push_back(0.0016203);
00620 
00621     this->mVariableNames.push_back("h");
00622     this->mVariableUnits.push_back("dimensionless");
00623     this->mInitialConditions.push_back(0.9944036);
00624 
00625     this->mVariableNames.push_back("d");
00626     this->mVariableUnits.push_back("dimensionless");
00627     this->mInitialConditions.push_back(0);
00628 
00629     this->mVariableNames.push_back("f");
00630     this->mVariableUnits.push_back("dimensionless");
00631     this->mInitialConditions.push_back(1);
00632 
00633     this->mVariableNames.push_back("f2");
00634     this->mVariableUnits.push_back("dimensionless");
00635     this->mInitialConditions.push_back(0.9349197);
00636 
00637     this->mVariableNames.push_back("f2ds");
00638     this->mVariableUnits.push_back("dimensionless");
00639     this->mInitialConditions.push_back(0.9651958);
00640 
00641     this->mVariableNames.push_back("s");
00642     this->mVariableUnits.push_back("dimensionless");
00643     this->mInitialConditions.push_back(0.9948645);
00644 
00645     this->mVariableNames.push_back("r");
00646     this->mVariableUnits.push_back("dimensionless");
00647     this->mInitialConditions.push_back(0);
00648 
00649     this->mVariableNames.push_back("ActFrac");
00650     this->mVariableUnits.push_back("dimensionless");
00651     this->mInitialConditions.push_back(0.0042614);
00652 
00653     this->mVariableNames.push_back("ProdFrac");
00654     this->mVariableUnits.push_back("dimensionless");
00655     this->mInitialConditions.push_back(0.4068154);
00656 
00657     this->mVariableNames.push_back("Na_i");
00658     this->mVariableUnits.push_back("millimolar");
00659     this->mInitialConditions.push_back(7.3321223);
00660 
00661     this->mVariableNames.push_back("K_i");
00662     this->mVariableUnits.push_back("millimolar");
00663     this->mInitialConditions.push_back(136.5644281);
00664 
00665     this->mVariableNames.push_back("Ca_i");
00666     this->mVariableUnits.push_back("millimolar");
00667     this->mInitialConditions.push_back(1.4e-5);
00668 
00669     this->mVariableNames.push_back("Ca_ds");
00670     this->mVariableUnits.push_back("millimolar");
00671     this->mInitialConditions.push_back(1.88e-5);
00672 
00673     this->mVariableNames.push_back("Ca_up");
00674     this->mVariableUnits.push_back("millimolar");
00675     this->mInitialConditions.push_back(0.4531889);
00676 
00677     this->mVariableNames.push_back("Ca_rel");
00678     this->mVariableUnits.push_back("millimolar");
00679     this->mInitialConditions.push_back(0.4481927);
00680 
00681     this->mVariableNames.push_back("Ca_Calmod");
00682     this->mVariableUnits.push_back("millimolar");
00683     this->mInitialConditions.push_back(0.0005555);
00684 
00685     this->mVariableNames.push_back("Ca_Trop");
00686     this->mVariableUnits.push_back("millimolar");
00687     this->mInitialConditions.push_back(0.0003542);
00688 
00689     this->mInitialised = true;
00690 }
00691 
00692 // Serialization for Boost>=1.36
00693 #include "SerializationExportWrapperForCpp.hpp"
00694 CHASTE_CLASS_EXPORT(CML_noble_varghese_kohl_noble_1998_basic_with_sac)