Chaste: /tmp/release_1/heart/src/odes/NobleVargheseKohlNoble1998.hpp Source File

00001 #ifndef _CML_noble_varghese_kohl_noble_1998_basic_
00002 #define _CML_noble_varghese_kohl_noble_1998_basic_
00003 
00004 // Model: noble_varghese_kohl_noble_1998_basic
00005 // Processed by pycml - CellML Tools in Python
00006 //     (translate: 4024, pycml: 4024)
00007 // on Wed Jul 16 14:55:26 2008
00008 
00009 #include <cmath>
00010 #include <cassert>
00011 #include "AbstractCardiacCell.hpp"
00012 #include "Exception.hpp"
00013 #include "AbstractStimulusFunction.hpp"
00014 #include "OdeSystemInformation.hpp"
00015 
00016 class CML_noble_varghese_kohl_noble_1998_basic : public AbstractCardiacCell
00017 {
00018 public:
00019     CML_noble_varghese_kohl_noble_1998_basic(AbstractIvpOdeSolver *pSolver,
00020                                              AbstractStimulusFunction *pIntracellularStimulus)
00021         : AbstractCardiacCell(pSolver, 22, 0, pIntracellularStimulus)
00022     {
00023         mpSystemInfo = OdeSystemInformation<CML_noble_varghese_kohl_noble_1998_basic>::Instance();
00024 
00025         Init();
00026 
00027     }
00028 
00029     ~CML_noble_varghese_kohl_noble_1998_basic(void)
00030     {
00031     }
00032 
00033     double 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         
00210         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;
00211         double value_in_microA = 0.001*value_in_nA;
00212         double value_in_microA_per_cm_squared = value_in_microA/0.000095;   
00213         return value_in_microA_per_cm_squared;
00214     }
00215 
00216     void EvaluateYDerivatives (
00217             double var_environment__time,
00218             const std::vector<double> &rY,
00219             std::vector<double> &rDY)
00220     {
00221         // Inputs:
00222         // Time units: second
00223         var_environment__time *= 0.001;
00224         double var_membrane__V = rY[0];
00225         // Units: millivolt; Initial value: -92.849333
00226         double var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1 = rY[1];
00227         // Units: dimensionless; Initial value: 1.03e-5
00228         double var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2 = rY[2];
00229         // Units: dimensionless; Initial value: 2e-7
00230         double var_slow_delayed_rectifier_potassium_current_xs_gate__xs = rY[3];
00231         // Units: dimensionless; Initial value: 0.001302
00232         double var_fast_sodium_current_m_gate__m = rY[4];
00233         // Units: dimensionless; Initial value: 0.0016203
00234         double var_fast_sodium_current_h_gate__h = rY[5];
00235         // Units: dimensionless; Initial value: 0.9944036
00236         double var_L_type_Ca_channel_d_gate__d = rY[6];
00237         // Units: dimensionless; Initial value: 0
00238         double var_L_type_Ca_channel_f_gate__f = rY[7];
00239         // Units: dimensionless; Initial value: 1
00240         double var_L_type_Ca_channel_f2_gate__f2 = rY[8];
00241         // Units: dimensionless; Initial value: 0.9349197
00242         double var_L_type_Ca_channel_f2ds_gate__f2ds = rY[9];
00243         // Units: dimensionless; Initial value: 0.9651958
00244         double var_transient_outward_current_s_gate__s = rY[10];
00245         // Units: dimensionless; Initial value: 0.9948645
00246         double var_transient_outward_current_r_gate__r = rY[11];
00247         // Units: dimensionless; Initial value: 0
00248         double var_calcium_release__ActFrac = rY[12];
00249         // Units: dimensionless; Initial value: 0.0042614
00250         double var_calcium_release__ProdFrac = rY[13];
00251         // Units: dimensionless; Initial value: 0.4068154
00252         double var_intracellular_sodium_concentration__Na_i = rY[14];
00253         // Units: millimolar; Initial value: 7.3321223
00254         double var_intracellular_potassium_concentration__K_i = rY[15];
00255         // Units: millimolar; Initial value: 136.5644281
00256         double var_intracellular_calcium_concentration__Ca_i = rY[16];
00257         // Units: millimolar; Initial value: 1.4e-5
00258         double var_intracellular_calcium_concentration__Ca_ds = rY[17];
00259         // Units: millimolar; Initial value: 1.88e-5
00260         double var_intracellular_calcium_concentration__Ca_up = rY[18];
00261         // Units: millimolar; Initial value: 0.4531889
00262         double var_intracellular_calcium_concentration__Ca_rel = rY[19];
00263         // Units: millimolar; Initial value: 0.4481927
00264         double var_intracellular_calcium_concentration__Ca_Calmod = rY[20];
00265         // Units: millimolar; Initial value: 0.0005555
00266         double var_intracellular_calcium_concentration__Ca_Trop = rY[21];
00267         // Units: millimolar; Initial value: 0.0003542
00268         
00269         
00270         // Mathematics
00271         const double var_membrane__R = 8314.472;
00272         const double var_membrane__T = 310.0;
00273         const double var_membrane__F = 96485.3415;
00274         const double var_membrane__Cm = 9.5e-05;
00275         double var_reversal_potentials__K_i = var_intracellular_potassium_concentration__K_i;
00276         double var_reversal_potentials__R = var_membrane__R;
00277         double var_reversal_potentials__T = var_membrane__T;
00278         double var_reversal_potentials__F = var_membrane__F;
00279         const double var_extracellular_potassium_concentration__K_o = 4.0;
00280         double var_reversal_potentials__K_o = var_extracellular_potassium_concentration__K_o;
00281         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);
00282         double var_time_independent_potassium_current__E_K = var_reversal_potentials__E_K;
00283         double var_time_independent_potassium_current__K_o = var_extracellular_potassium_concentration__K_o;
00284         double var_time_independent_potassium_current__R = var_membrane__R;
00285         double var_time_independent_potassium_current__V = var_membrane__V;
00286         double var_time_independent_potassium_current__T = var_membrane__T;
00287         const double var_time_independent_potassium_current__K_mk1 = 10.0;
00288         const double var_time_independent_potassium_current__g_K1 = 0.5;
00289         double var_time_independent_potassium_current__F = var_membrane__F;
00290         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)));
00291         double var_membrane__i_K1 = var_time_independent_potassium_current__i_K1;
00292         double var_transient_outward_current__s = var_transient_outward_current_s_gate__s;
00293         double var_transient_outward_current__r = var_transient_outward_current_r_gate__r;
00294         const double var_transient_outward_current__g_to = 0.005;
00295         double var_transient_outward_current__V = var_membrane__V;
00296         double var_transient_outward_current__E_K = var_reversal_potentials__E_K;
00297         const double var_transient_outward_current__g_tos = 0.0;
00298         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);
00299         double var_membrane__i_to = var_transient_outward_current__i_to;
00300         const double var_rapid_delayed_rectifier_potassium_current__g_Kr2 = 0.0013;
00301         const double var_rapid_delayed_rectifier_potassium_current__g_Kr1 = 0.0021;
00302         double var_rapid_delayed_rectifier_potassium_current__xr1 = var_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1;
00303         double var_rapid_delayed_rectifier_potassium_current__xr2 = var_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2;
00304         double var_rapid_delayed_rectifier_potassium_current__V = var_membrane__V;
00305         double var_rapid_delayed_rectifier_potassium_current__E_K = var_reversal_potentials__E_K;
00306         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);
00307         double var_membrane__i_Kr = var_rapid_delayed_rectifier_potassium_current__i_Kr;
00308         double var_slow_delayed_rectifier_potassium_current__xs = var_slow_delayed_rectifier_potassium_current_xs_gate__xs;
00309         const double var_extracellular_sodium_concentration__Na_o = 140.0;
00310         double var_reversal_potentials__Na_o = var_extracellular_sodium_concentration__Na_o;
00311         double var_reversal_potentials__Na_i = var_intracellular_sodium_concentration__Na_i;
00312         const double var_reversal_potentials__P_kna = 0.03;
00313         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)));
00314         double var_slow_delayed_rectifier_potassium_current__E_Ks = var_reversal_potentials__E_Ks;
00315         const double var_slow_delayed_rectifier_potassium_current__g_Ks = 0.0026;
00316         double var_slow_delayed_rectifier_potassium_current__V = var_membrane__V;
00317         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);
00318         double var_membrane__i_Ks = var_slow_delayed_rectifier_potassium_current__i_Ks;
00319         double var_L_type_Ca_channel__d = var_L_type_Ca_channel_d_gate__d;
00320         const double var_L_type_Ca_channel__FrICa = 1.0;
00321         double var_L_type_Ca_channel__f = var_L_type_Ca_channel_f_gate__f;
00322         double var_L_type_Ca_channel__K_o = var_extracellular_potassium_concentration__K_o;
00323         double var_L_type_Ca_channel__K_i = var_intracellular_potassium_concentration__K_i;
00324         double var_L_type_Ca_channel__F = var_membrane__F;
00325         const double var_L_type_Ca_channel__P_Ca_L = 0.1;
00326         double var_L_type_Ca_channel__T = var_membrane__T;
00327         const double var_L_type_Ca_channel__P_CaK = 0.002;
00328         double var_L_type_Ca_channel__V = var_membrane__V;
00329         double var_L_type_Ca_channel__f2 = var_L_type_Ca_channel_f2_gate__f2;
00330         double var_L_type_Ca_channel__R = var_membrane__R;
00331         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))));
00332         double var_membrane__i_Ca_L_K_cyt = var_L_type_Ca_channel__i_Ca_L_K_cyt;
00333         double var_L_type_Ca_channel__f2ds = var_L_type_Ca_channel_f2ds_gate__f2ds;
00334         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))));
00335         double var_membrane__i_Ca_L_K_ds = var_L_type_Ca_channel__i_Ca_L_K_ds;
00336         const double var_sodium_potassium_pump__i_NaK_max = 0.7;
00337         double var_sodium_potassium_pump__Na_i = var_intracellular_sodium_concentration__Na_i;
00338         double var_sodium_potassium_pump__K_o = var_extracellular_potassium_concentration__K_o;
00339         const double var_sodium_potassium_pump__K_mNa = 40.0;
00340         const double var_sodium_potassium_pump__K_mK = 1.0;
00341         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);
00342         double var_membrane__i_NaK = var_sodium_potassium_pump__i_NaK;
00343         const double var_fast_sodium_current__g_Na = 2.5;
00344         double var_fast_sodium_current__h = var_fast_sodium_current_h_gate__h;
00345         double var_fast_sodium_current__V = var_membrane__V;
00346         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)));
00347         double var_fast_sodium_current__E_mh = var_reversal_potentials__E_mh;
00348         double var_fast_sodium_current__m = var_fast_sodium_current_m_gate__m;
00349         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);
00350         double var_membrane__i_Na = var_fast_sodium_current__i_Na;
00351         double var_sodium_background_current__V = var_membrane__V;
00352         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);
00353         double var_sodium_background_current__E_Na = var_reversal_potentials__E_Na;
00354         const double var_sodium_background_current__g_bna = 0.0006;
00355         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);
00356         double var_membrane__i_b_Na = var_sodium_background_current__i_b_Na;
00357         const double var_persistent_sodium_current__g_pna = 0.004;
00358         double var_persistent_sodium_current__V = var_membrane__V;
00359         double var_persistent_sodium_current__E_Na = var_reversal_potentials__E_Na;
00360         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);
00361         double var_membrane__i_p_Na = var_persistent_sodium_current__i_p_Na;
00362         const double var_L_type_Ca_channel__P_CaNa = 0.01;
00363         double var_L_type_Ca_channel__Na_o = var_extracellular_sodium_concentration__Na_o;
00364         double var_L_type_Ca_channel__Na_i = var_intracellular_sodium_concentration__Na_i;
00365         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))));
00366         double var_membrane__i_Ca_L_Na_cyt = var_L_type_Ca_channel__i_Ca_L_Na_cyt;
00367         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))));
00368         double var_membrane__i_Ca_L_Na_ds = var_L_type_Ca_channel__i_Ca_L_Na_ds;
00369         double var_sodium_calcium_exchanger__Na_i = var_intracellular_sodium_concentration__Na_i;
00370         const double var_sodium_calcium_exchanger__n_NaCa = 3.0;
00371         const double var_sodium_calcium_exchanger__gamma = 0.5;
00372         double var_sodium_calcium_exchanger__F = var_membrane__F;
00373         double var_sodium_calcium_exchanger__Na_o = var_extracellular_sodium_concentration__Na_o;
00374         const double var_sodium_calcium_exchanger__FRiNaCa = 0.001;
00375         double var_sodium_calcium_exchanger__R = var_membrane__R;
00376         double var_sodium_calcium_exchanger__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00377         double var_sodium_calcium_exchanger__T = var_membrane__T;
00378         double var_sodium_calcium_exchanger__V = var_membrane__V;
00379         const double var_sodium_calcium_exchanger__d_NaCa = 0.0;
00380         const double var_extracellular_calcium_concentration__Ca_o = 2.0;
00381         double var_sodium_calcium_exchanger__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00382         const double var_sodium_calcium_exchanger__k_NaCa = 0.0005;
00383         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)));
00384         double var_membrane__i_NaCa_cyt = var_sodium_calcium_exchanger__i_NaCa_cyt;
00385         double var_sodium_calcium_exchanger__Ca_ds = var_intracellular_calcium_concentration__Ca_ds;
00386         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)));
00387         double var_membrane__i_NaCa_ds = var_sodium_calcium_exchanger__i_NaCa_ds;
00388         double var_L_type_Ca_channel__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00389         double var_L_type_Ca_channel__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00390         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))));
00391         double var_membrane__i_Ca_L_Ca_cyt = var_L_type_Ca_channel__i_Ca_L_Ca_cyt;
00392         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))));
00393         double var_membrane__i_Ca_L_Ca_ds = var_L_type_Ca_channel__i_Ca_L_Ca_ds;
00394         double var_reversal_potentials__Ca_o = var_extracellular_calcium_concentration__Ca_o;
00395         double var_reversal_potentials__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00396         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);
00397         double var_calcium_background_current__E_Ca = var_reversal_potentials__E_Ca;
00398         const double var_calcium_background_current__g_bca = 0.00025;
00399         double var_calcium_background_current__V = var_membrane__V;
00400         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);
00401         double var_membrane__i_b_Ca = var_calcium_background_current__i_b_Ca;
00402         double var_membrane__i_Stim = GetStimulus((1.0/0.001)*var_environment__time);
00403         double var_rapid_delayed_rectifier_potassium_current_xr1_gate__V = var_rapid_delayed_rectifier_potassium_current__V;
00404         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));
00405         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);
00406         double var_rapid_delayed_rectifier_potassium_current_xr2_gate__V = var_rapid_delayed_rectifier_potassium_current__V;
00407         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));
00408         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));
00409         double var_slow_delayed_rectifier_potassium_current_xs_gate__V = var_slow_delayed_rectifier_potassium_current__V;
00410         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));
00411         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);
00412         double var_fast_sodium_current_m_gate__V = var_fast_sodium_current__V;
00413         double var_fast_sodium_current_m_gate__E0_m = var_fast_sodium_current_m_gate__V + 41.0;
00414         const double var_fast_sodium_current_m_gate__delta_m = 1e-05;
00415         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)));
00416         double var_fast_sodium_current_m_gate__beta_m = 8000.0 * exp((-0.056) * (var_fast_sodium_current_m_gate__V + 66.0));
00417         double var_fast_sodium_current_h_gate__V = var_fast_sodium_current__V;
00418         const double var_fast_sodium_current_h_gate__shift_h = 0.0;
00419         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));
00420         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))));
00421         double var_L_type_Ca_channel__Ca_ds = var_intracellular_calcium_concentration__Ca_ds;
00422         const double var_L_type_Ca_channel__Km_f2 = 100000.0;
00423         const double var_L_type_Ca_channel__Km_f2ds = 0.001;
00424         const double var_L_type_Ca_channel__R_decay = 20.0;
00425         double var_L_type_Ca_channel_d_gate__V = var_L_type_Ca_channel__V;
00426         double var_L_type_Ca_channel_d_gate__E0_d = (var_L_type_Ca_channel_d_gate__V + 24.0) - 5.0;
00427         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)));
00428         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));
00429         const double var_L_type_Ca_channel_d_gate__speed_d = 3.0;
00430         double var_L_type_Ca_channel_f_gate__V = var_L_type_Ca_channel__V;
00431         double var_L_type_Ca_channel_f_gate__E0_f = var_L_type_Ca_channel_f_gate__V + 34.0;
00432         const double var_L_type_Ca_channel_f_gate__delta_f = 0.0001;
00433         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));
00434         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));
00435         const double var_L_type_Ca_channel_f_gate__speed_f = 0.3;
00436         double var_L_type_Ca_channel_f2_gate__Km_f2 = var_L_type_Ca_channel__Km_f2;
00437         double var_L_type_Ca_channel_f2_gate__Ca_i = var_L_type_Ca_channel__Ca_i;
00438         double var_L_type_Ca_channel_f2ds_gate__Km_f2ds = var_L_type_Ca_channel__Km_f2ds;
00439         double var_L_type_Ca_channel_f2ds_gate__R_decay = var_L_type_Ca_channel__R_decay;
00440         double var_L_type_Ca_channel_f2ds_gate__Ca_ds = var_L_type_Ca_channel__Ca_ds;
00441         double var_transient_outward_current_s_gate__V = var_transient_outward_current__V;
00442         double var_transient_outward_current_s_gate__alpha_s = 0.033 * exp((-var_transient_outward_current_s_gate__V) / 17.0);
00443         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)));
00444         double var_transient_outward_current_r_gate__V = var_transient_outward_current__V;
00445         double var_sarcoplasmic_reticulum_calcium_pump__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00446         double var_sarcoplasmic_reticulum_calcium_pump__Ca_up = var_intracellular_calcium_concentration__Ca_up;
00447         const double var_sarcoplasmic_reticulum_calcium_pump__alpha_up = 0.4;
00448         const double var_sarcoplasmic_reticulum_calcium_pump__beta_up = 0.03;
00449         const double var_sarcoplasmic_reticulum_calcium_pump__K_srca = 0.5;
00450         const double var_sarcoplasmic_reticulum_calcium_pump__K_xcs = 0.4;
00451         const double var_sarcoplasmic_reticulum_calcium_pump__K_cyca = 0.0003;
00452         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;
00453         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;
00454         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);
00455         double var_calcium_translocation__Ca_rel = var_intracellular_calcium_concentration__Ca_rel;
00456         double var_calcium_translocation__Ca_up = var_intracellular_calcium_concentration__Ca_up;
00457         double var_calcium_translocation__i_trans = 50.0 * (var_calcium_translocation__Ca_up - var_calcium_translocation__Ca_rel);
00458         const double var_calcium_release__K_m_rel = 250.0;
00459         const double var_calcium_release__K_leak_rate = 0.05;
00460         double var_calcium_release__Ca_rel = var_intracellular_calcium_concentration__Ca_rel;
00461         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;
00462         double var_calcium_release__V = var_membrane__V;
00463         double var_calcium_release__VoltDep = exp(0.08 * (var_calcium_release__V - 40.0));
00464         const double var_calcium_release__K_m_Ca_cyt = 0.0005;
00465         double var_calcium_release__Ca_i = var_intracellular_calcium_concentration__Ca_i;
00466         double var_calcium_release__CaiReg = var_calcium_release__Ca_i / (var_calcium_release__Ca_i + var_calcium_release__K_m_Ca_cyt);
00467         double var_calcium_release__Ca_ds = var_intracellular_calcium_concentration__Ca_ds;
00468         const double var_calcium_release__K_m_Ca_ds = 0.01;
00469         double var_calcium_release__CadsReg = var_calcium_release__Ca_ds / (var_calcium_release__Ca_ds + var_calcium_release__K_m_Ca_ds);
00470         double var_calcium_release__RegBindSite = var_calcium_release__CaiReg + ((1.0 - var_calcium_release__CaiReg) * var_calcium_release__CadsReg);
00471         double var_calcium_release__ActRate = (0.0 * var_calcium_release__VoltDep) + (500.0 * pow(var_calcium_release__RegBindSite, 2.0));
00472         double var_calcium_release__InactRate = 60.0 + (500.0 * pow(var_calcium_release__RegBindSite, 2.0));
00473         double var_calcium_release__PrecFrac = (1.0 - var_calcium_release__ActFrac) - var_calcium_release__ProdFrac;
00474         double var_calcium_release__SpeedRel = (var_calcium_release__V < (-50.0)) ? 5.0 : 1.0;
00475         const double var_intracellular_calcium_concentration__V_up_ratio = 0.01;
00476         const double var_intracellular_calcium_concentration__V_rel_ratio = 0.1;
00477         const double var_intracellular_calcium_concentration__V_e_ratio = 0.4;
00478         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;
00479         const double var_intracellular_calcium_concentration__radius = 0.012;
00480         const double var_intracellular_calcium_concentration__length = 0.074;
00481         double var_intracellular_calcium_concentration__V_Cell = 3.141592654 * pow(var_intracellular_calcium_concentration__radius, 2.0) * var_intracellular_calcium_concentration__length;
00482         double var_intracellular_calcium_concentration__V_i = var_intracellular_calcium_concentration__V_Cell * var_intracellular_calcium_concentration__V_i_ratio;
00483         double var_intracellular_sodium_concentration__V_i = var_intracellular_calcium_concentration__V_i;
00484         double var_intracellular_sodium_concentration__F = var_membrane__F;
00485         double var_intracellular_sodium_concentration__i_Na = var_fast_sodium_current__i_Na;
00486         double var_intracellular_sodium_concentration__i_b_Na = var_sodium_background_current__i_b_Na;
00487         double var_intracellular_sodium_concentration__i_p_Na = var_persistent_sodium_current__i_p_Na;
00488         double var_intracellular_sodium_concentration__i_Ca_L_Na_cyt = var_L_type_Ca_channel__i_Ca_L_Na_cyt;
00489         double var_intracellular_sodium_concentration__i_Ca_L_Na_ds = var_L_type_Ca_channel__i_Ca_L_Na_ds;
00490         double var_intracellular_sodium_concentration__i_NaK = var_sodium_potassium_pump__i_NaK;
00491         double var_intracellular_sodium_concentration__i_NaCa_cyt = var_sodium_calcium_exchanger__i_NaCa_cyt;
00492         double var_intracellular_potassium_concentration__V_i = var_intracellular_calcium_concentration__V_i;
00493         double var_intracellular_potassium_concentration__i_K1 = var_time_independent_potassium_current__i_K1;
00494         double var_intracellular_potassium_concentration__i_Kr = var_rapid_delayed_rectifier_potassium_current__i_Kr;
00495         double var_intracellular_potassium_concentration__i_Ks = var_slow_delayed_rectifier_potassium_current__i_Ks;
00496         double var_intracellular_potassium_concentration__i_Ca_L_K_cyt = var_L_type_Ca_channel__i_Ca_L_K_cyt;
00497         double var_intracellular_potassium_concentration__i_Ca_L_K_ds = var_L_type_Ca_channel__i_Ca_L_K_ds;
00498         double var_intracellular_potassium_concentration__i_to = var_transient_outward_current__i_to;
00499         double var_intracellular_potassium_concentration__i_NaK = var_sodium_potassium_pump__i_NaK;
00500         double var_intracellular_potassium_concentration__F = var_membrane__F;
00501         const double var_intracellular_calcium_concentration__Calmod = 0.02;
00502         const double var_intracellular_calcium_concentration__Trop = 0.05;
00503         const double var_intracellular_calcium_concentration__alpha_Calmod = 100000.0;
00504         const double var_intracellular_calcium_concentration__beta_Calmod = 50.0;
00505         const double var_intracellular_calcium_concentration__alpha_Trop = 100000.0;
00506         const double var_intracellular_calcium_concentration__beta_Trop = 200.0;
00507         const double var_intracellular_calcium_concentration__V_ds_ratio = 0.1;
00508         const double var_intracellular_calcium_concentration__Kdecay = 10.0;
00509         double var_intracellular_calcium_concentration__i_up = var_sarcoplasmic_reticulum_calcium_pump__i_up;
00510         double var_intracellular_calcium_concentration__i_trans = var_calcium_translocation__i_trans;
00511         double var_intracellular_calcium_concentration__i_rel = var_calcium_release__i_rel;
00512         double var_intracellular_calcium_concentration__i_NaCa_cyt = var_sodium_calcium_exchanger__i_NaCa_cyt;
00513         double var_intracellular_calcium_concentration__i_NaCa_ds = var_sodium_calcium_exchanger__i_NaCa_ds;
00514         double var_intracellular_calcium_concentration__i_Ca_L_Ca_cyt = var_L_type_Ca_channel__i_Ca_L_Ca_cyt;
00515         double var_intracellular_calcium_concentration__i_Ca_L_Ca_ds = var_L_type_Ca_channel__i_Ca_L_Ca_ds;
00516         double var_intracellular_calcium_concentration__i_b_Ca = var_calcium_background_current__i_b_Ca;
00517         double var_intracellular_calcium_concentration__F = var_membrane__F;
00518         double d_dt_membrane__V;
00519         if (mSetVoltageDerivativeToZero)
00520         {
00521             d_dt_membrane__V = 0.0;
00522         }
00523         else
00524         {
00525             d_dt_membrane__V = ((-1.0) / var_membrane__Cm) * (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);
00526         }
00527         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);
00528         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);
00529         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);
00530         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);
00531         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);
00532         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));
00533         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));
00534         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));
00535         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));
00536         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);
00537         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);
00538         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);
00539         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);
00540         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);
00541         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));
00542         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);
00543         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);
00544         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;
00545         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);
00546         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;
00547         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;
00548         
00549         rDY[0] = 0.001*d_dt_membrane__V;
00550         rDY[1] = 0.001*d_dt_rapid_delayed_rectifier_potassium_current_xr1_gate__xr1;
00551         rDY[2] = 0.001*d_dt_rapid_delayed_rectifier_potassium_current_xr2_gate__xr2;
00552         rDY[3] = 0.001*d_dt_slow_delayed_rectifier_potassium_current_xs_gate__xs;
00553         rDY[4] = 0.001*d_dt_fast_sodium_current_m_gate__m;
00554         rDY[5] = 0.001*d_dt_fast_sodium_current_h_gate__h;
00555         rDY[6] = 0.001*d_dt_L_type_Ca_channel_d_gate__d;
00556         rDY[7] = 0.001*d_dt_L_type_Ca_channel_f_gate__f;
00557         rDY[8] = 0.001*d_dt_L_type_Ca_channel_f2_gate__f2;
00558         rDY[9] = 0.001*d_dt_L_type_Ca_channel_f2ds_gate__f2ds;
00559         rDY[10] = 0.001*d_dt_transient_outward_current_s_gate__s;
00560         rDY[11] = 0.001*d_dt_transient_outward_current_r_gate__r;
00561         rDY[12] = 0.001*d_dt_calcium_release__ActFrac;
00562         rDY[13] = 0.001*d_dt_calcium_release__ProdFrac;
00563         rDY[14] = 0.001*d_dt_intracellular_sodium_concentration__Na_i;
00564         rDY[15] = 0.001*d_dt_intracellular_potassium_concentration__K_i;
00565         rDY[16] = 0.001*d_dt_intracellular_calcium_concentration__Ca_i;
00566         rDY[17] = 0.001*d_dt_intracellular_calcium_concentration__Ca_ds;
00567         rDY[18] = 0.001*d_dt_intracellular_calcium_concentration__Ca_up;
00568         rDY[19] = 0.001*d_dt_intracellular_calcium_concentration__Ca_rel;
00569         rDY[20] = 0.001*d_dt_intracellular_calcium_concentration__Ca_Calmod;
00570         rDY[21] = 0.001*d_dt_intracellular_calcium_concentration__Ca_Trop;
00571     }
00572 
00573 };
00574 
00575 
00576 template<>
00577 void OdeSystemInformation<CML_noble_varghese_kohl_noble_1998_basic>::Initialise(void)
00578 {
00579     // Time units: second
00580     // 
00581     this->mVariableNames.push_back("V");
00582     this->mVariableUnits.push_back("millivolt");
00583     this->mInitialConditions.push_back(-92.849333);
00584 
00585     this->mVariableNames.push_back("xr1");
00586     this->mVariableUnits.push_back("dimensionless");
00587     this->mInitialConditions.push_back(1.03e-5);
00588 
00589     this->mVariableNames.push_back("xr2");
00590     this->mVariableUnits.push_back("dimensionless");
00591     this->mInitialConditions.push_back(2e-7);
00592 
00593     this->mVariableNames.push_back("xs");
00594     this->mVariableUnits.push_back("dimensionless");
00595     this->mInitialConditions.push_back(0.001302);
00596 
00597     this->mVariableNames.push_back("m");
00598     this->mVariableUnits.push_back("dimensionless");
00599     this->mInitialConditions.push_back(0.0016203);
00600 
00601     this->mVariableNames.push_back("h");
00602     this->mVariableUnits.push_back("dimensionless");
00603     this->mInitialConditions.push_back(0.9944036);
00604 
00605     this->mVariableNames.push_back("d");
00606     this->mVariableUnits.push_back("dimensionless");
00607     this->mInitialConditions.push_back(0);
00608 
00609     this->mVariableNames.push_back("f");
00610     this->mVariableUnits.push_back("dimensionless");
00611     this->mInitialConditions.push_back(1);
00612 
00613     this->mVariableNames.push_back("f2");
00614     this->mVariableUnits.push_back("dimensionless");
00615     this->mInitialConditions.push_back(0.9349197);
00616 
00617     this->mVariableNames.push_back("f2ds");
00618     this->mVariableUnits.push_back("dimensionless");
00619     this->mInitialConditions.push_back(0.9651958);
00620 
00621     this->mVariableNames.push_back("s");
00622     this->mVariableUnits.push_back("dimensionless");
00623     this->mInitialConditions.push_back(0.9948645);
00624 
00625     this->mVariableNames.push_back("r");
00626     this->mVariableUnits.push_back("dimensionless");
00627     this->mInitialConditions.push_back(0);
00628 
00629     this->mVariableNames.push_back("ActFrac");
00630     this->mVariableUnits.push_back("dimensionless");
00631     this->mInitialConditions.push_back(0.0042614);
00632 
00633     this->mVariableNames.push_back("ProdFrac");
00634     this->mVariableUnits.push_back("dimensionless");
00635     this->mInitialConditions.push_back(0.4068154);
00636 
00637     this->mVariableNames.push_back("Na_i");
00638     this->mVariableUnits.push_back("millimolar");
00639     this->mInitialConditions.push_back(7.3321223);
00640 
00641     this->mVariableNames.push_back("K_i");
00642     this->mVariableUnits.push_back("millimolar");
00643     this->mInitialConditions.push_back(136.5644281);
00644 
00645     this->mVariableNames.push_back("Ca_i");
00646     this->mVariableUnits.push_back("millimolar");
00647     this->mInitialConditions.push_back(1.4e-5);
00648 
00649     this->mVariableNames.push_back("Ca_ds");
00650     this->mVariableUnits.push_back("millimolar");
00651     this->mInitialConditions.push_back(1.88e-5);
00652 
00653     this->mVariableNames.push_back("Ca_up");
00654     this->mVariableUnits.push_back("millimolar");
00655     this->mInitialConditions.push_back(0.4531889);
00656 
00657     this->mVariableNames.push_back("Ca_rel");
00658     this->mVariableUnits.push_back("millimolar");
00659     this->mInitialConditions.push_back(0.4481927);
00660 
00661     this->mVariableNames.push_back("Ca_Calmod");
00662     this->mVariableUnits.push_back("millimolar");
00663     this->mInitialConditions.push_back(0.0005555);
00664 
00665     this->mVariableNames.push_back("Ca_Trop");
00666     this->mVariableUnits.push_back("millimolar");
00667     this->mInitialConditions.push_back(0.0003542);
00668     
00669     this->mInitialised = true;
00670 }
00671 
00672 
00673 #endif