Generated Code
The following is c code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
There are a total of 68 entries in the algebraic variable array.
There are a total of 25 entries in each of the rate and state variable arrays.
There are a total of 42 entries in the constant variable array.
*/
/*
* VOI is time in component environment (millisecond).
* STATES[0] is V in component membrane (millivolt).
* CONSTANTS[0] is R in component membrane (joule_per_mole_kelvin).
* CONSTANTS[1] is T in component membrane (kelvin).
* CONSTANTS[2] is F in component membrane (coulomb_per_millimole).
* CONSTANTS[3] is Cm in component membrane (picoF).
* CONSTANTS[4] is i_stim in component membrane (picoA_per_picoF).
* ALGEBRAIC[14] is i_Na in component fast_sodium_current (picoA_per_picoF).
* ALGEBRAIC[34] is i_K1 in component time_independent_potassium_current (picoA_per_picoF).
* ALGEBRAIC[41] is i_to in component transient_outward_K_current (picoA_per_picoF).
* ALGEBRAIC[43] is i_Kur_d in component ultrarapid_delayed_rectifier_K_current (picoA_per_picoF).
* ALGEBRAIC[44] is i_Kr in component rapid_delayed_rectifier_K_current (picoA_per_picoF).
* ALGEBRAIC[45] is i_Ks in component slow_delayed_rectifier_K_current (picoA_per_picoF).
* ALGEBRAIC[46] is i_Ca in component sarcolemmal_Ca_current (picoA_per_picoF).
* ALGEBRAIC[48] is i_Cl_Ca in component Ca_activated_Cl_current (picoA_per_picoF).
* ALGEBRAIC[55] is i_p_Ca in component Ca_pump_current (picoA_per_picoF).
* ALGEBRAIC[50] is i_NaK in component sodium_potassium_pump (picoA_per_picoF).
* ALGEBRAIC[51] is i_NaCa in component Na_Ca_exchanger_current (picoA_per_picoF).
* ALGEBRAIC[52] is i_B_Na in component background_currents (picoA_per_picoF).
* ALGEBRAIC[54] is i_B_Ca in component background_currents (picoA_per_picoF).
* ALGEBRAIC[0] is E_Na in component fast_sodium_current (millivolt).
* CONSTANTS[5] is g_Na in component fast_sodium_current (nanoS_per_picoF).
* STATES[1] is Na_i in component intracellular_ion_concentrations (millimolar).
* CONSTANTS[6] is Na_o in component standard_ionic_concentrations (millimolar).
* STATES[2] is m in component fast_sodium_current_m_gate (dimensionless).
* STATES[3] is h in component fast_sodium_current_h_gate (dimensionless).
* STATES[4] is j in component fast_sodium_current_j_gate (dimensionless).
* ALGEBRAIC[1] is alpha_m in component fast_sodium_current_m_gate (per_millisecond).
* ALGEBRAIC[15] is beta_m in component fast_sodium_current_m_gate (per_millisecond).
* ALGEBRAIC[2] is alpha_h in component fast_sodium_current_h_gate (per_millisecond).
* ALGEBRAIC[16] is beta_h in component fast_sodium_current_h_gate (per_millisecond).
* ALGEBRAIC[3] is alpha_j in component fast_sodium_current_j_gate (per_millisecond).
* ALGEBRAIC[17] is beta_j in component fast_sodium_current_j_gate (per_millisecond).
* ALGEBRAIC[27] is E_K in component time_independent_potassium_current (millivolt).
* CONSTANTS[7] is g_K1 in component time_independent_potassium_current (nanoS_per_picoF).
* CONSTANTS[8] is K_o in component standard_ionic_concentrations (millimolar).
* STATES[5] is K_i in component intracellular_ion_concentrations (millimolar).
* CONSTANTS[9] is g_to in component transient_outward_K_current (nanoS_per_picoF).
* STATES[6] is oa in component transient_outward_K_current_oa_gate (dimensionless).
* STATES[7] is oi in component transient_outward_K_current_oi_gate (dimensionless).
* ALGEBRAIC[4] is alpha_oa in component transient_outward_K_current_oa_gate (per_millisecond).
* ALGEBRAIC[18] is beta_oa in component transient_outward_K_current_oa_gate (per_millisecond).
* ALGEBRAIC[28] is tau_oa in component transient_outward_K_current_oa_gate (millisecond).
* ALGEBRAIC[35] is oa_infinity in component transient_outward_K_current_oa_gate (dimensionless).
* ALGEBRAIC[5] is alpha_oi in component transient_outward_K_current_oi_gate (per_millisecond).
* ALGEBRAIC[19] is beta_oi in component transient_outward_K_current_oi_gate (per_millisecond).
* ALGEBRAIC[29] is tau_oi in component transient_outward_K_current_oi_gate (millisecond).
* ALGEBRAIC[36] is oi_infinity in component transient_outward_K_current_oi_gate (dimensionless).
* ALGEBRAIC[42] is g_Kur_d in component ultrarapid_delayed_rectifier_K_current (nanoS_per_picoF).
* STATES[8] is ua in component ultrarapid_delayed_rectifier_K_current_ua_gate (dimensionless).
* STATES[9] is ui in component ultrarapid_delayed_rectifier_K_current_ui_gate (dimensionless).
* ALGEBRAIC[6] is alpha_ua in component ultrarapid_delayed_rectifier_K_current_ua_gate (per_millisecond).
* ALGEBRAIC[20] is beta_ua in component ultrarapid_delayed_rectifier_K_current_ua_gate (per_millisecond).
* ALGEBRAIC[30] is tau_ua in component ultrarapid_delayed_rectifier_K_current_ua_gate (millisecond).
* ALGEBRAIC[37] is ua_infinity in component ultrarapid_delayed_rectifier_K_current_ua_gate (dimensionless).
* ALGEBRAIC[7] is alpha_ui in component ultrarapid_delayed_rectifier_K_current_ui_gate (per_millisecond).
* ALGEBRAIC[21] is beta_ui in component ultrarapid_delayed_rectifier_K_current_ui_gate (per_millisecond).
* ALGEBRAIC[31] is tau_ui in component ultrarapid_delayed_rectifier_K_current_ui_gate (millisecond).
* ALGEBRAIC[38] is ui_infinity in component ultrarapid_delayed_rectifier_K_current_ui_gate (dimensionless).
* CONSTANTS[10] is g_Kr in component rapid_delayed_rectifier_K_current (nanoS_per_picoF).
* STATES[10] is xr in component rapid_delayed_rectifier_K_current_xr_gate (dimensionless).
* ALGEBRAIC[8] is alpha_xr in component rapid_delayed_rectifier_K_current_xr_gate (per_millisecond).
* ALGEBRAIC[22] is beta_xr in component rapid_delayed_rectifier_K_current_xr_gate (per_millisecond).
* ALGEBRAIC[32] is tau_xr in component rapid_delayed_rectifier_K_current_xr_gate (millisecond).
* ALGEBRAIC[39] is xr_infinity in component rapid_delayed_rectifier_K_current_xr_gate (dimensionless).
* CONSTANTS[11] is g_Ks in component slow_delayed_rectifier_K_current (nanoS_per_picoF).
* STATES[11] is xs in component slow_delayed_rectifier_K_current_xs_gate (dimensionless).
* ALGEBRAIC[9] is alpha_xs in component slow_delayed_rectifier_K_current_xs_gate (per_millisecond).
* ALGEBRAIC[23] is beta_xs in component slow_delayed_rectifier_K_current_xs_gate (per_millisecond).
* ALGEBRAIC[33] is tau_xs in component slow_delayed_rectifier_K_current_xs_gate (millisecond).
* ALGEBRAIC[40] is xs_infinity in component slow_delayed_rectifier_K_current_xs_gate (dimensionless).
* CONSTANTS[12] is g_Ca in component sarcolemmal_Ca_current (nanoS_per_picoF).
* STATES[12] is Ca_i in component intracellular_ion_concentrations (millimolar).
* STATES[13] is d in component sarcolemmal_Ca_current_d_gate (dimensionless).
* STATES[14] is f in component sarcolemmal_Ca_current_f_gate (dimensionless).
* STATES[15] is f_Ca in component sarcolemmal_Ca_current_f_Ca_gate (dimensionless).
* ALGEBRAIC[10] is d_infinity in component sarcolemmal_Ca_current_d_gate (dimensionless).
* ALGEBRAIC[24] is tau_d in component sarcolemmal_Ca_current_d_gate (millisecond).
* ALGEBRAIC[11] is f_infinity in component sarcolemmal_Ca_current_f_gate (dimensionless).
* ALGEBRAIC[25] is tau_f in component sarcolemmal_Ca_current_f_gate (millisecond).
* ALGEBRAIC[12] is f_Ca_infinity in component sarcolemmal_Ca_current_f_Ca_gate (dimensionless).
* CONSTANTS[13] is tau_f_Ca in component sarcolemmal_Ca_current_f_Ca_gate (millisecond).
* CONSTANTS[14] is g_Cl_Ca in component Ca_activated_Cl_current (nanoS_per_picoF).
* ALGEBRAIC[47] is E_Cl in component Ca_activated_Cl_current (millivolt).
* ALGEBRAIC[57] is Fn in component Ca_release_current_from_JSR (dimensionless).
* STATES[16] is Cl_i in component intracellular_ion_concentrations (millimolar).
* CONSTANTS[15] is Cl_o in component standard_ionic_concentrations (millimolar).
* CONSTANTS[16] is q_Ca in component Ca_activated_Cl_current_q_Ca_gate (dimensionless).
* ALGEBRAIC[59] is q_Ca_infinity in component Ca_activated_Cl_current_q_Ca_gate (dimensionless).
* CONSTANTS[17] is Km_Na_i in component sodium_potassium_pump (millimolar).
* CONSTANTS[18] is Km_K_o in component sodium_potassium_pump (millimolar).
* CONSTANTS[19] is i_NaK_max in component sodium_potassium_pump (picoA_per_picoF).
* ALGEBRAIC[49] is f_NaK in component sodium_potassium_pump (dimensionless).
* CONSTANTS[40] is sigma in component sodium_potassium_pump (dimensionless).
* CONSTANTS[20] is I_NaCa_max in component Na_Ca_exchanger_current (picoA_per_picoF).
* CONSTANTS[21] is K_mNa in component Na_Ca_exchanger_current (millimolar).
* CONSTANTS[22] is K_mCa in component Na_Ca_exchanger_current (millimolar).
* CONSTANTS[23] is K_sat in component Na_Ca_exchanger_current (dimensionless).
* CONSTANTS[24] is Ca_o in component standard_ionic_concentrations (millimolar).
* CONSTANTS[25] is g_B_Na in component background_currents (nanoS_per_picoF).
* CONSTANTS[26] is g_B_Ca in component background_currents (nanoS_per_picoF).
* ALGEBRAIC[53] is E_Ca in component background_currents (millivolt).
* CONSTANTS[27] is i_p_Ca_max in component Ca_pump_current (picoA_per_picoF).
* ALGEBRAIC[56] is i_rel in component Ca_release_current_from_JSR (picoA_per_picoF).
* CONSTANTS[28] is K_rel in component Ca_release_current_from_JSR (per_millisecond).
* CONSTANTS[29] is V_rel in component Ca_release_current_from_JSR (micrometre_3).
* STATES[17] is Ca_rel in component intracellular_ion_concentrations (millimolar).
* STATES[18] is u in component Ca_release_current_from_JSR_u_gate (dimensionless).
* STATES[19] is v in component Ca_release_current_from_JSR_v_gate (dimensionless).
* STATES[20] is w in component Ca_release_current_from_JSR_w_gate (dimensionless).
* CONSTANTS[41] is tau_u in component Ca_release_current_from_JSR_u_gate (millisecond).
* ALGEBRAIC[60] is u_infinity in component Ca_release_current_from_JSR_u_gate (dimensionless).
* ALGEBRAIC[61] is tau_v in component Ca_release_current_from_JSR_v_gate (millisecond).
* ALGEBRAIC[63] is v_infinity in component Ca_release_current_from_JSR_v_gate (dimensionless).
* ALGEBRAIC[13] is tau_w in component Ca_release_current_from_JSR_w_gate (millisecond).
* ALGEBRAIC[26] is w_infinity in component Ca_release_current_from_JSR_w_gate (dimensionless).
* ALGEBRAIC[58] is i_tr in component transfer_current_from_NSR_to_JSR (picoA_per_picoF).
* CONSTANTS[30] is tau_tr in component transfer_current_from_NSR_to_JSR (millisecond).
* STATES[21] is Ca_up in component intracellular_ion_concentrations (millimolar).
* CONSTANTS[31] is I_up_max in component Ca_uptake_current_by_the_NSR (picoA_per_picoF).
* ALGEBRAIC[62] is i_up in component Ca_uptake_current_by_the_NSR (picoA_per_picoF).
* CONSTANTS[32] is K_up in component Ca_uptake_current_by_the_NSR (millimolar).
* ALGEBRAIC[65] is i_up_leak in component Ca_leak_current_by_the_NSR (picoA_per_picoF).
* CONSTANTS[33] is Ca_up_max in component Ca_leak_current_by_the_NSR (millimolar).
* CONSTANTS[34] is CMDN_max in component Ca_buffers (millimolar).
* CONSTANTS[35] is TRPN_max in component Ca_buffers (millimolar).
* CONSTANTS[36] is CSQN_max in component Ca_buffers (millimolar).
* ALGEBRAIC[66] is J_Ca_CMDN in component Ca_buffers (millimolar_per_millisecond).
* ALGEBRAIC[67] is J_Ca_TRPN in component Ca_buffers (millimolar_per_millisecond).
* ALGEBRAIC[64] is J_Ca_CSQN in component Ca_buffers (millimolar_per_millisecond).
* STATES[22] is Ca_CMDN in component Ca_buffers (millimolar).
* STATES[23] is Ca_TRPN in component Ca_buffers (millimolar).
* STATES[24] is Ca_CSQN in component Ca_buffers (millimolar).
* CONSTANTS[37] is V_i in component intracellular_ion_concentrations (micrometre_3).
* CONSTANTS[38] is V_rel in component intracellular_ion_concentrations (micrometre_3).
* CONSTANTS[39] is V_up in component intracellular_ion_concentrations (micrometre_3).
* RATES[0] is d/dt V in component membrane (millivolt).
* RATES[2] is d/dt m in component fast_sodium_current_m_gate (dimensionless).
* RATES[3] is d/dt h in component fast_sodium_current_h_gate (dimensionless).
* RATES[4] is d/dt j in component fast_sodium_current_j_gate (dimensionless).
* RATES[6] is d/dt oa in component transient_outward_K_current_oa_gate (dimensionless).
* RATES[7] is d/dt oi in component transient_outward_K_current_oi_gate (dimensionless).
* RATES[8] is d/dt ua in component ultrarapid_delayed_rectifier_K_current_ua_gate (dimensionless).
* RATES[9] is d/dt ui in component ultrarapid_delayed_rectifier_K_current_ui_gate (dimensionless).
* RATES[10] is d/dt xr in component rapid_delayed_rectifier_K_current_xr_gate (dimensionless).
* RATES[11] is d/dt xs in component slow_delayed_rectifier_K_current_xs_gate (dimensionless).
* RATES[13] is d/dt d in component sarcolemmal_Ca_current_d_gate (dimensionless).
* RATES[14] is d/dt f in component sarcolemmal_Ca_current_f_gate (dimensionless).
* RATES[15] is d/dt f_Ca in component sarcolemmal_Ca_current_f_Ca_gate (dimensionless).
* RATES[18] is d/dt u in component Ca_release_current_from_JSR_u_gate (dimensionless).
* RATES[19] is d/dt v in component Ca_release_current_from_JSR_v_gate (dimensionless).
* RATES[20] is d/dt w in component Ca_release_current_from_JSR_w_gate (dimensionless).
* RATES[22] is d/dt Ca_CMDN in component Ca_buffers (millimolar).
* RATES[23] is d/dt Ca_TRPN in component Ca_buffers (millimolar).
* RATES[24] is d/dt Ca_CSQN in component Ca_buffers (millimolar).
* RATES[1] is d/dt Na_i in component intracellular_ion_concentrations (millimolar).
* RATES[5] is d/dt K_i in component intracellular_ion_concentrations (millimolar).
* RATES[16] is d/dt Cl_i in component intracellular_ion_concentrations (millimolar).
* RATES[12] is d/dt Ca_i in component intracellular_ion_concentrations (millimolar).
* RATES[21] is d/dt Ca_up in component intracellular_ion_concentrations (millimolar).
* RATES[17] is d/dt Ca_rel in component intracellular_ion_concentrations (millimolar).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -83.53;
CONSTANTS[0] = 8.3143;
CONSTANTS[1] = 310.0;
CONSTANTS[2] = 96.4867;
CONSTANTS[3] = 100.0;
CONSTANTS[4] = -2900.0;
CONSTANTS[5] = 7.8;
STATES[1] = 11.75;
CONSTANTS[6] = 140.0;
STATES[2] = 0.001972;
STATES[3] = 0.9791;
STATES[4] = 0.9869;
CONSTANTS[7] = 0.15;
CONSTANTS[8] = 5.4;
STATES[5] = 138.4;
CONSTANTS[9] = 0.19824;
STATES[6] = 0.07164;
STATES[7] = 0.9980;
STATES[8] = 0.05869;
STATES[9] = 0.9987;
CONSTANTS[10] = 0.06984;
STATES[10] = 0.0000007433;
CONSTANTS[11] = 0.0561;
STATES[11] = 0.01791;
CONSTANTS[12] = 0.24;
STATES[12] = 0.0001024;
STATES[13] = 0.000004757;
STATES[14] = 0.9999;
STATES[15] = 0.7484;
CONSTANTS[13] = 2.0;
CONSTANTS[14] = 0.3;
STATES[16] = 29.26;
CONSTANTS[15] = 132.0;
CONSTANTS[16] = 0.0;
CONSTANTS[17] = 10.0;
CONSTANTS[18] = 1.5;
CONSTANTS[19] = 0.6;
CONSTANTS[20] = 1600.0;
CONSTANTS[21] = 87.5;
CONSTANTS[22] = 1.38;
CONSTANTS[23] = 0.1;
CONSTANTS[24] = 1.8;
CONSTANTS[25] = 0.000674;
CONSTANTS[26] = 0.00113;
CONSTANTS[27] = 0.275;
CONSTANTS[28] = 30.0;
CONSTANTS[29] = 96.48;
STATES[17] = 1.502;
STATES[18] = 0.0;
STATES[19] = 1.0;
STATES[20] = 0.9993;
CONSTANTS[30] = 180.0;
STATES[21] = 1.502;
CONSTANTS[31] = 0.005;
CONSTANTS[32] = 0.00092;
CONSTANTS[33] = 15.0;
CONSTANTS[34] = 0.045;
CONSTANTS[35] = 0.35;
CONSTANTS[36] = 10.0;
STATES[22] = 0.001856;
STATES[23] = 0.007022;
STATES[24] = 6.432;
CONSTANTS[37] = 13668.0;
CONSTANTS[38] = 96.48;
CONSTANTS[39] = 1109.52;
CONSTANTS[40] = (1.00000/7.00000)*(exp(CONSTANTS[6]/67.3000) - 1.00000);
CONSTANTS[41] = 8.00000;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[12] = 0.290000+ 0.800000*pow(1.00000+exp((STATES[12] - 0.000120000)/6.00000e-05), -1.00000);
RATES[15] = (ALGEBRAIC[12] - STATES[15])/CONSTANTS[13];
ALGEBRAIC[1] = 0.320000*((STATES[0]+47.1300)/(1.00000 - exp( -0.100000*(STATES[0]+47.1300))));
ALGEBRAIC[15] = 0.0800000*exp(STATES[0]/-11.0000);
RATES[2] = ALGEBRAIC[1]*(1.00000 - STATES[2]) - ALGEBRAIC[15]*STATES[2];
ALGEBRAIC[2] = (STATES[0]<-40.0000 ? 0.135000*exp((STATES[0]+80.0000)/-6.80000) : 0.00000);
ALGEBRAIC[16] = (STATES[0]<-40.0000 ? 3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0]) : 1.00000/( 0.130000*(1.00000+exp((STATES[0]+10.6600)/-11.1000))));
RATES[3] = ALGEBRAIC[2]*(1.00000 - STATES[3]) - ALGEBRAIC[16]*STATES[3];
ALGEBRAIC[3] = (STATES[0]<-40.0000 ? (( -127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( -0.0439100*STATES[0]))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))))*(STATES[0]+37.7800) : 0.00000);
ALGEBRAIC[17] = (STATES[0]<-40.0000 ? ( 0.121200*exp( -0.0105200*STATES[0]))/(1.00000+exp( -0.137800*(STATES[0]+40.1400))) : ( 0.300000*exp( -2.53500e-07*STATES[0]))/(1.00000+exp( -0.100000*(STATES[0]+32.0000))));
RATES[4] = ALGEBRAIC[3]*(1.00000 - STATES[4]) - ALGEBRAIC[17]*STATES[4];
ALGEBRAIC[10] = pow(1.00000+exp((STATES[0]+10.0000)/-6.00000), -1.00000);
ALGEBRAIC[24] = (1.00000 - exp((STATES[0]+10.0000)/-6.24000))/( 0.0350000*(STATES[0]+10.0000)*(1.00000+exp((STATES[0]+10.0000)/-6.24000)));
RATES[13] = (ALGEBRAIC[10] - STATES[13])/ALGEBRAIC[24];
ALGEBRAIC[11] = pow(1.00000+exp((STATES[0]+24.6000)/6.20000), -1.00000);
ALGEBRAIC[25] = 400.000*pow(1.00000+ 4.50000*exp( -0.000700000*pow(STATES[0] - 9.00000, 2.00000)), -1.00000);
RATES[14] = (ALGEBRAIC[11] - STATES[14])/ALGEBRAIC[25];
ALGEBRAIC[13] = (6.00000 - 6.00000*exp((STATES[0] - 7.90000)/-5.00000))/( (1.00000+ 0.300000*exp((STATES[0] - 7.90000)/-5.00000))*(STATES[0] - 7.90000));
ALGEBRAIC[26] = 1.00000 - pow(1.00000+exp((STATES[0] - 40.0000)/-17.0000), -1.00000);
RATES[20] = (ALGEBRAIC[26] - STATES[20])/ALGEBRAIC[13];
ALGEBRAIC[4] = 0.650000*pow(exp((STATES[0]+18.0000)/-8.50000)+exp((STATES[0] - 16.0000)/-59.0000), -1.00000);
ALGEBRAIC[18] = 1.20000*pow(2.20000+exp((STATES[0]+75.0000)/18.0000), -1.00000);
ALGEBRAIC[28] = pow(ALGEBRAIC[4]+ALGEBRAIC[18], -1.00000);
ALGEBRAIC[35] = pow(1.00000+exp((STATES[0]+0.500000)/-10.5000), - (1.00000/3.00000));
RATES[6] = (ALGEBRAIC[35] - STATES[6])/ALGEBRAIC[28];
ALGEBRAIC[5] = pow(6.20000+exp((STATES[0]+105.200)/9.85000), -1.00000);
ALGEBRAIC[19] = pow(7.54000+exp((STATES[0] - 8.89000)/-12.8700), -1.00000);
ALGEBRAIC[29] = pow(ALGEBRAIC[5]+ALGEBRAIC[19], -1.00000);
ALGEBRAIC[36] = pow(1.00000+exp((STATES[0]+43.3770)/6.45000), -1.00000);
RATES[7] = (ALGEBRAIC[36] - STATES[7])/ALGEBRAIC[29];
ALGEBRAIC[6] = 1.47000*pow(exp((STATES[0]+33.2000)/-30.6300)+exp((STATES[0] - 27.6000)/-30.6500), -1.00000);
ALGEBRAIC[20] = 0.420000*pow(exp((STATES[0]+26.6400)/2.49000)+exp((STATES[0]+44.4100)/20.3600), -1.00000);
ALGEBRAIC[30] = pow(ALGEBRAIC[6]+ALGEBRAIC[20], -1.00000);
ALGEBRAIC[37] = pow(1.00000+exp((STATES[0]+2.81000)/-9.49000), - (1.00000/3.00000));
RATES[8] = (ALGEBRAIC[37] - STATES[8])/ALGEBRAIC[30];
ALGEBRAIC[7] = pow(21.0000+exp((STATES[0] - 185.000)/-28.0000), -1.00000);
ALGEBRAIC[21] = exp((STATES[0] - 158.000)/16.0000);
ALGEBRAIC[31] = pow(ALGEBRAIC[7]+ALGEBRAIC[21], -1.00000);
ALGEBRAIC[38] = pow(1.00000+exp((STATES[0] - 99.4500)/27.4800), -1.00000);
RATES[9] = (ALGEBRAIC[38] - STATES[9])/ALGEBRAIC[31];
ALGEBRAIC[8] = 0.0400000*((STATES[0] - 248.000)/(1.00000 - exp((STATES[0] - 248.000)/-28.0000)));
ALGEBRAIC[22] = 0.0280000*((STATES[0]+163.000)/(exp((STATES[0]+163.000)/21.0000) - 1.00000));
ALGEBRAIC[32] = pow(ALGEBRAIC[8]+ALGEBRAIC[22], -1.00000);
ALGEBRAIC[39] = pow(1.00000+exp((STATES[0]+7.65400)/-5.37700), -1.00000);
RATES[10] = (ALGEBRAIC[39] - STATES[10])/ALGEBRAIC[32];
ALGEBRAIC[9] = 1.00000e-05*((STATES[0]+28.5000)/(1.00000 - exp((STATES[0]+28.5000)/-115.000)));
ALGEBRAIC[23] = 0.000230000*((STATES[0]+28.5000)/(exp((STATES[0]+28.5000)/3.30000) - 1.00000));
ALGEBRAIC[33] = pow(ALGEBRAIC[9]+ALGEBRAIC[23], -1.00000);
ALGEBRAIC[40] = pow(1.00000+exp((STATES[0] - 13.0000)/-12.0000), -0.500000);
RATES[11] = (ALGEBRAIC[40] - STATES[11])/ALGEBRAIC[33];
ALGEBRAIC[47] = (( CONSTANTS[0]*CONSTANTS[1])/( -1.00000*CONSTANTS[2]))*log(CONSTANTS[15]/STATES[16]);
ALGEBRAIC[48] = CONSTANTS[14]*CONSTANTS[16]*(STATES[0] - ALGEBRAIC[47]);
RATES[16] = ALGEBRAIC[48]/( CONSTANTS[37]*CONSTANTS[2]);
ALGEBRAIC[27] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[8]/STATES[5]);
ALGEBRAIC[34] = ( CONSTANTS[7]*(STATES[0] - ALGEBRAIC[27]))/(1.00000+exp( 0.0700000*(STATES[0]+80.0000)));
ALGEBRAIC[41] = CONSTANTS[9]*pow(STATES[6], 3.00000)*STATES[7]*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[42] = 0.00855000+0.0779000/(1.00000+exp((STATES[0]+11.0000)/-16.0000));
ALGEBRAIC[43] = ALGEBRAIC[42]*pow(STATES[8], 3.00000)*STATES[9]*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[44] = CONSTANTS[10]*STATES[10]*(0.0700000+0.580000/(1.00000+exp((STATES[0]+15.0000)/22.4000)))*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[45] = CONSTANTS[11]*pow(STATES[11], 2.00000)*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[49] = pow(1.00000+ 0.124500*exp( -0.100000*(( CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1])))+ 0.0365000*CONSTANTS[40]*exp(- (( CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]))), -1.00000);
ALGEBRAIC[50] = CONSTANTS[19]*ALGEBRAIC[49]*(1.00000/(1.00000+pow(CONSTANTS[17]/STATES[1], 1.50000)))*(CONSTANTS[8]/(CONSTANTS[8]+CONSTANTS[18]));
RATES[5] = ( 2.00000*ALGEBRAIC[50] - (ALGEBRAIC[34]+ALGEBRAIC[41]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[45]))/( CONSTANTS[37]*CONSTANTS[2]);
ALGEBRAIC[0] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[6]/STATES[1]);
ALGEBRAIC[14] = CONSTANTS[5]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*(STATES[0] - ALGEBRAIC[0]);
ALGEBRAIC[51] = ( CONSTANTS[20]*( exp(( 0.350000*CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[1], 3.00000)*CONSTANTS[24] - exp(( -0.650000*CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]))*pow(CONSTANTS[6], 3.00000)*STATES[12]))/( (pow(CONSTANTS[21], 3.00000)+pow(CONSTANTS[6], 3.00000))*(CONSTANTS[22]+CONSTANTS[24])*(1.00000+ CONSTANTS[23]*exp(( -0.650000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))));
ALGEBRAIC[52] = CONSTANTS[25]*(STATES[0] - ALGEBRAIC[0]);
RATES[1] = ( -3.00000*ALGEBRAIC[50] - ( 3.00000*ALGEBRAIC[51]+ALGEBRAIC[52]+ALGEBRAIC[14]))/( CONSTANTS[37]*CONSTANTS[2]);
ALGEBRAIC[46] = CONSTANTS[12]*STATES[13]*STATES[14]*STATES[15]*(STATES[0] - 65.0000);
ALGEBRAIC[55] = CONSTANTS[27]*(STATES[12]/(0.000500000+STATES[12]));
ALGEBRAIC[53] = (( CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]))*log(CONSTANTS[24]/STATES[12]);
ALGEBRAIC[54] = CONSTANTS[26]*(STATES[0] - ALGEBRAIC[53]);
RATES[0] = - (ALGEBRAIC[14]+ALGEBRAIC[34]+ALGEBRAIC[41]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[46]+ALGEBRAIC[48]+ALGEBRAIC[55]+ALGEBRAIC[51]+ALGEBRAIC[50]+ALGEBRAIC[52]+ALGEBRAIC[54]+CONSTANTS[4])/CONSTANTS[3];
RATES[24] = 0.480000*STATES[17]*(1.00000 - STATES[24]/CONSTANTS[36]) - 0.400000*(STATES[24]/CONSTANTS[36]);
ALGEBRAIC[56] = CONSTANTS[28]*pow(STATES[18], 2.00000)*STATES[19]*STATES[20]*(STATES[17] - STATES[12]);
ALGEBRAIC[57] = 1.00000e-12*CONSTANTS[29]*ALGEBRAIC[56] - 5.00000e-13*( (1.00000/( 2.00000*CONSTANTS[2]))*ALGEBRAIC[46] - (1.00000/( 5.00000*CONSTANTS[2]))*ALGEBRAIC[51]);
ALGEBRAIC[60] = pow(1.00000+exp((ALGEBRAIC[57] - 3.41750e-13)/-1.36700e-15), -1.00000);
RATES[18] = (ALGEBRAIC[60] - STATES[18])/CONSTANTS[41];
ALGEBRAIC[61] = 1.91000+ 2.09000*pow(1.00000+exp((ALGEBRAIC[57] - 3.41750e-13)/-1.36700e-15), -1.00000);
ALGEBRAIC[63] = 1.00000 - pow(1.00000+exp((ALGEBRAIC[57] - 6.83500e-14)/-1.36700e-15), -1.00000);
RATES[19] = (ALGEBRAIC[63] - STATES[19])/ALGEBRAIC[61];
RATES[22] = 200.000*STATES[12]*(1.00000 - STATES[22]/CONSTANTS[34]) - 0.476000*(STATES[22]/CONSTANTS[34]);
ALGEBRAIC[58] = (STATES[21] - STATES[17])/CONSTANTS[30];
ALGEBRAIC[62] = CONSTANTS[31]/(1.00000+CONSTANTS[32]/STATES[12]);
ALGEBRAIC[65] = CONSTANTS[31]*(STATES[21]/CONSTANTS[33]);
RATES[21] = ALGEBRAIC[62] - (ALGEBRAIC[65]+ ALGEBRAIC[58]*(CONSTANTS[38]/CONSTANTS[39]));
ALGEBRAIC[64] = RATES[24];
RATES[17] = ALGEBRAIC[58] - (ALGEBRAIC[56]+ 31.0000*ALGEBRAIC[64]);
RATES[23] = 78.4000*STATES[12]*(1.00000 - STATES[23]/CONSTANTS[35]) - 0.392000*(STATES[23]/CONSTANTS[35]);
ALGEBRAIC[66] = RATES[22];
ALGEBRAIC[67] = RATES[23];
RATES[12] = (( 2.00000*ALGEBRAIC[51] - (ALGEBRAIC[55]+ALGEBRAIC[46]+ALGEBRAIC[54]))/( 2.00000*CONSTANTS[37]*CONSTANTS[2])+( CONSTANTS[39]*(ALGEBRAIC[65] - ALGEBRAIC[62])+ ALGEBRAIC[56]*CONSTANTS[38])/CONSTANTS[37]) - ( CONSTANTS[35]*ALGEBRAIC[67]+ CONSTANTS[34]*ALGEBRAIC[66]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[12] = 0.290000+ 0.800000*pow(1.00000+exp((STATES[12] - 0.000120000)/6.00000e-05), -1.00000);
ALGEBRAIC[1] = 0.320000*((STATES[0]+47.1300)/(1.00000 - exp( -0.100000*(STATES[0]+47.1300))));
ALGEBRAIC[15] = 0.0800000*exp(STATES[0]/-11.0000);
ALGEBRAIC[2] = (STATES[0]<-40.0000 ? 0.135000*exp((STATES[0]+80.0000)/-6.80000) : 0.00000);
ALGEBRAIC[16] = (STATES[0]<-40.0000 ? 3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0]) : 1.00000/( 0.130000*(1.00000+exp((STATES[0]+10.6600)/-11.1000))));
ALGEBRAIC[3] = (STATES[0]<-40.0000 ? (( -127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( -0.0439100*STATES[0]))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))))*(STATES[0]+37.7800) : 0.00000);
ALGEBRAIC[17] = (STATES[0]<-40.0000 ? ( 0.121200*exp( -0.0105200*STATES[0]))/(1.00000+exp( -0.137800*(STATES[0]+40.1400))) : ( 0.300000*exp( -2.53500e-07*STATES[0]))/(1.00000+exp( -0.100000*(STATES[0]+32.0000))));
ALGEBRAIC[10] = pow(1.00000+exp((STATES[0]+10.0000)/-6.00000), -1.00000);
ALGEBRAIC[24] = (1.00000 - exp((STATES[0]+10.0000)/-6.24000))/( 0.0350000*(STATES[0]+10.0000)*(1.00000+exp((STATES[0]+10.0000)/-6.24000)));
ALGEBRAIC[11] = pow(1.00000+exp((STATES[0]+24.6000)/6.20000), -1.00000);
ALGEBRAIC[25] = 400.000*pow(1.00000+ 4.50000*exp( -0.000700000*pow(STATES[0] - 9.00000, 2.00000)), -1.00000);
ALGEBRAIC[13] = (6.00000 - 6.00000*exp((STATES[0] - 7.90000)/-5.00000))/( (1.00000+ 0.300000*exp((STATES[0] - 7.90000)/-5.00000))*(STATES[0] - 7.90000));
ALGEBRAIC[26] = 1.00000 - pow(1.00000+exp((STATES[0] - 40.0000)/-17.0000), -1.00000);
ALGEBRAIC[4] = 0.650000*pow(exp((STATES[0]+18.0000)/-8.50000)+exp((STATES[0] - 16.0000)/-59.0000), -1.00000);
ALGEBRAIC[18] = 1.20000*pow(2.20000+exp((STATES[0]+75.0000)/18.0000), -1.00000);
ALGEBRAIC[28] = pow(ALGEBRAIC[4]+ALGEBRAIC[18], -1.00000);
ALGEBRAIC[35] = pow(1.00000+exp((STATES[0]+0.500000)/-10.5000), - (1.00000/3.00000));
ALGEBRAIC[5] = pow(6.20000+exp((STATES[0]+105.200)/9.85000), -1.00000);
ALGEBRAIC[19] = pow(7.54000+exp((STATES[0] - 8.89000)/-12.8700), -1.00000);
ALGEBRAIC[29] = pow(ALGEBRAIC[5]+ALGEBRAIC[19], -1.00000);
ALGEBRAIC[36] = pow(1.00000+exp((STATES[0]+43.3770)/6.45000), -1.00000);
ALGEBRAIC[6] = 1.47000*pow(exp((STATES[0]+33.2000)/-30.6300)+exp((STATES[0] - 27.6000)/-30.6500), -1.00000);
ALGEBRAIC[20] = 0.420000*pow(exp((STATES[0]+26.6400)/2.49000)+exp((STATES[0]+44.4100)/20.3600), -1.00000);
ALGEBRAIC[30] = pow(ALGEBRAIC[6]+ALGEBRAIC[20], -1.00000);
ALGEBRAIC[37] = pow(1.00000+exp((STATES[0]+2.81000)/-9.49000), - (1.00000/3.00000));
ALGEBRAIC[7] = pow(21.0000+exp((STATES[0] - 185.000)/-28.0000), -1.00000);
ALGEBRAIC[21] = exp((STATES[0] - 158.000)/16.0000);
ALGEBRAIC[31] = pow(ALGEBRAIC[7]+ALGEBRAIC[21], -1.00000);
ALGEBRAIC[38] = pow(1.00000+exp((STATES[0] - 99.4500)/27.4800), -1.00000);
ALGEBRAIC[8] = 0.0400000*((STATES[0] - 248.000)/(1.00000 - exp((STATES[0] - 248.000)/-28.0000)));
ALGEBRAIC[22] = 0.0280000*((STATES[0]+163.000)/(exp((STATES[0]+163.000)/21.0000) - 1.00000));
ALGEBRAIC[32] = pow(ALGEBRAIC[8]+ALGEBRAIC[22], -1.00000);
ALGEBRAIC[39] = pow(1.00000+exp((STATES[0]+7.65400)/-5.37700), -1.00000);
ALGEBRAIC[9] = 1.00000e-05*((STATES[0]+28.5000)/(1.00000 - exp((STATES[0]+28.5000)/-115.000)));
ALGEBRAIC[23] = 0.000230000*((STATES[0]+28.5000)/(exp((STATES[0]+28.5000)/3.30000) - 1.00000));
ALGEBRAIC[33] = pow(ALGEBRAIC[9]+ALGEBRAIC[23], -1.00000);
ALGEBRAIC[40] = pow(1.00000+exp((STATES[0] - 13.0000)/-12.0000), -0.500000);
ALGEBRAIC[47] = (( CONSTANTS[0]*CONSTANTS[1])/( -1.00000*CONSTANTS[2]))*log(CONSTANTS[15]/STATES[16]);
ALGEBRAIC[48] = CONSTANTS[14]*CONSTANTS[16]*(STATES[0] - ALGEBRAIC[47]);
ALGEBRAIC[27] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[8]/STATES[5]);
ALGEBRAIC[34] = ( CONSTANTS[7]*(STATES[0] - ALGEBRAIC[27]))/(1.00000+exp( 0.0700000*(STATES[0]+80.0000)));
ALGEBRAIC[41] = CONSTANTS[9]*pow(STATES[6], 3.00000)*STATES[7]*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[42] = 0.00855000+0.0779000/(1.00000+exp((STATES[0]+11.0000)/-16.0000));
ALGEBRAIC[43] = ALGEBRAIC[42]*pow(STATES[8], 3.00000)*STATES[9]*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[44] = CONSTANTS[10]*STATES[10]*(0.0700000+0.580000/(1.00000+exp((STATES[0]+15.0000)/22.4000)))*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[45] = CONSTANTS[11]*pow(STATES[11], 2.00000)*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[49] = pow(1.00000+ 0.124500*exp( -0.100000*(( CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1])))+ 0.0365000*CONSTANTS[40]*exp(- (( CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]))), -1.00000);
ALGEBRAIC[50] = CONSTANTS[19]*ALGEBRAIC[49]*(1.00000/(1.00000+pow(CONSTANTS[17]/STATES[1], 1.50000)))*(CONSTANTS[8]/(CONSTANTS[8]+CONSTANTS[18]));
ALGEBRAIC[0] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[6]/STATES[1]);
ALGEBRAIC[14] = CONSTANTS[5]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*(STATES[0] - ALGEBRAIC[0]);
ALGEBRAIC[51] = ( CONSTANTS[20]*( exp(( 0.350000*CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[1], 3.00000)*CONSTANTS[24] - exp(( -0.650000*CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]))*pow(CONSTANTS[6], 3.00000)*STATES[12]))/( (pow(CONSTANTS[21], 3.00000)+pow(CONSTANTS[6], 3.00000))*(CONSTANTS[22]+CONSTANTS[24])*(1.00000+ CONSTANTS[23]*exp(( -0.650000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))));
ALGEBRAIC[52] = CONSTANTS[25]*(STATES[0] - ALGEBRAIC[0]);
ALGEBRAIC[46] = CONSTANTS[12]*STATES[13]*STATES[14]*STATES[15]*(STATES[0] - 65.0000);
ALGEBRAIC[55] = CONSTANTS[27]*(STATES[12]/(0.000500000+STATES[12]));
ALGEBRAIC[53] = (( CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]))*log(CONSTANTS[24]/STATES[12]);
ALGEBRAIC[54] = CONSTANTS[26]*(STATES[0] - ALGEBRAIC[53]);
ALGEBRAIC[56] = CONSTANTS[28]*pow(STATES[18], 2.00000)*STATES[19]*STATES[20]*(STATES[17] - STATES[12]);
ALGEBRAIC[57] = 1.00000e-12*CONSTANTS[29]*ALGEBRAIC[56] - 5.00000e-13*( (1.00000/( 2.00000*CONSTANTS[2]))*ALGEBRAIC[46] - (1.00000/( 5.00000*CONSTANTS[2]))*ALGEBRAIC[51]);
ALGEBRAIC[60] = pow(1.00000+exp((ALGEBRAIC[57] - 3.41750e-13)/-1.36700e-15), -1.00000);
ALGEBRAIC[61] = 1.91000+ 2.09000*pow(1.00000+exp((ALGEBRAIC[57] - 3.41750e-13)/-1.36700e-15), -1.00000);
ALGEBRAIC[63] = 1.00000 - pow(1.00000+exp((ALGEBRAIC[57] - 6.83500e-14)/-1.36700e-15), -1.00000);
ALGEBRAIC[58] = (STATES[21] - STATES[17])/CONSTANTS[30];
ALGEBRAIC[62] = CONSTANTS[31]/(1.00000+CONSTANTS[32]/STATES[12]);
ALGEBRAIC[65] = CONSTANTS[31]*(STATES[21]/CONSTANTS[33]);
ALGEBRAIC[64] = RATES[24];
ALGEBRAIC[66] = RATES[22];
ALGEBRAIC[67] = RATES[23];
ALGEBRAIC[59] = 1.00000 - pow(1.00000+pow(ALGEBRAIC[57]/1.10000e-10, 3.00000), -1.00000);
}