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 65 entries in the algebraic variable array.
There are a total of 36 entries in each of the rate and state variable arrays.
There are a total of 100 entries in the constant variable array.
*/
/*
* VOI is time in component environment (ms).
* CONSTANTS[0] is lambda in component environment (dimensionless).
* CONSTANTS[1] is dlambdadt in component environment (per_ms).
* STATES[0] is V in component membrane (mV).
* CONSTANTS[2] is R in component membrane (gas_constant_units).
* CONSTANTS[3] is T in component membrane (kelvin).
* CONSTANTS[4] is F in component membrane (faradays_constant_units).
* ALGEBRAIC[17] is i_Na in component fast_sodium_current (uA_per_mm2).
* ALGEBRAIC[22] is i_Ca_L_Ca in component L_type_Ca_channel (uA_per_mm2).
* ALGEBRAIC[24] is i_Ca_L_K in component L_type_Ca_channel (uA_per_mm2).
* ALGEBRAIC[28] is i_K in component time_dependent_potassium_current (uA_per_mm2).
* ALGEBRAIC[34] is i_K1 in component time_independent_potassium_current (uA_per_mm2).
* ALGEBRAIC[38] is i_NaCa in component Na_Ca_exchanger (uA_per_mm2).
* ALGEBRAIC[37] is i_Kp in component plateau_potassium_current (uA_per_mm2).
* ALGEBRAIC[39] is i_p_Ca in component sarcolemmal_calcium_pump (uA_per_mm2).
* ALGEBRAIC[41] is i_Na_b in component sodium_background_current (uA_per_mm2).
* ALGEBRAIC[43] is i_Ca_b in component calcium_background_current (uA_per_mm2).
* ALGEBRAIC[46] is i_NaK in component sodium_potassium_pump (uA_per_mm2).
* ALGEBRAIC[63] is i_ns_Ca in component non_specific_calcium_activated_current (uA_per_mm2).
* CONSTANTS[5] is Cm in component membrane (uF_per_mm2).
* ALGEBRAIC[6] is I_stim in component membrane (uA_per_mm2).
* CONSTANTS[6] is stim_start in component membrane (ms).
* CONSTANTS[7] is stim_end in component membrane (ms).
* CONSTANTS[8] is stim_period in component membrane (ms).
* CONSTANTS[9] is stim_duration in component membrane (ms).
* CONSTANTS[10] is stim_amplitude in component membrane (uA_per_mm2).
* ALGEBRAIC[14] is E_Na in component fast_sodium_current (mV).
* CONSTANTS[11] is g_Na in component fast_sodium_current (mS_per_mm2).
* STATES[1] is Nai in component ionic_concentrations (mM).
* CONSTANTS[12] is Nao in component ionic_concentrations (mM).
* 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[0] is alpha_m in component fast_sodium_current_m_gate (per_ms).
* ALGEBRAIC[8] is beta_m in component fast_sodium_current_m_gate (per_ms).
* ALGEBRAIC[1] is alpha_h in component fast_sodium_current_h_gate (per_ms).
* ALGEBRAIC[9] is beta_h in component fast_sodium_current_h_gate (per_ms).
* ALGEBRAIC[2] is alpha_j in component fast_sodium_current_j_gate (per_ms).
* ALGEBRAIC[10] is beta_j in component fast_sodium_current_j_gate (per_ms).
* CONSTANTS[13] is P_Ca in component L_type_Ca_channel (mm_per_ms).
* CONSTANTS[14] is P_K in component L_type_Ca_channel (mm_per_ms).
* ALGEBRAIC[23] is p_k in component L_type_Ca_channel (mm_per_ms).
* CONSTANTS[15] is i_Ca_L_Ca_half in component L_type_Ca_channel (uA_per_mm2).
* ALGEBRAIC[20] is i_Ca_L_Ca_max in component L_type_Ca_channel (uA_per_mm2).
* STATES[5] is O in component L_type_Ca_channel (dimensionless).
* STATES[6] is O_Ca in component L_type_Ca_channel (dimensionless).
* ALGEBRAIC[3] is alpha in component L_type_Ca_channel (per_ms).
* ALGEBRAIC[11] is beta in component L_type_Ca_channel (per_ms).
* ALGEBRAIC[21] is gamma in component L_type_Ca_channel (per_ms).
* ALGEBRAIC[16] is alpha_a in component L_type_Ca_channel (per_ms).
* ALGEBRAIC[19] is beta_b in component L_type_Ca_channel (per_ms).
* CONSTANTS[16] is a in component L_type_Ca_channel (dimensionless).
* CONSTANTS[17] is b in component L_type_Ca_channel (dimensionless).
* CONSTANTS[18] is g in component L_type_Ca_channel (per_ms).
* CONSTANTS[19] is f in component L_type_Ca_channel (per_ms).
* CONSTANTS[20] is g_ in component L_type_Ca_channel (per_ms).
* CONSTANTS[21] is f_ in component L_type_Ca_channel (per_ms).
* CONSTANTS[22] is omega in component L_type_Ca_channel (per_ms).
* STATES[7] is C0 in component L_type_Ca_channel (dimensionless).
* STATES[8] is C1 in component L_type_Ca_channel (dimensionless).
* STATES[9] is C2 in component L_type_Ca_channel (dimensionless).
* STATES[10] is C3 in component L_type_Ca_channel (dimensionless).
* STATES[11] is C4 in component L_type_Ca_channel (dimensionless).
* STATES[12] is C_Ca0 in component L_type_Ca_channel (dimensionless).
* STATES[13] is C_Ca1 in component L_type_Ca_channel (dimensionless).
* STATES[14] is C_Ca2 in component L_type_Ca_channel (dimensionless).
* STATES[15] is C_Ca3 in component L_type_Ca_channel (dimensionless).
* STATES[16] is C_Ca4 in component L_type_Ca_channel (dimensionless).
* STATES[17] is Ca_SS in component calcium_subsystem (mM).
* CONSTANTS[23] is Cao in component ionic_concentrations (mM).
* STATES[18] is Ko in component ionic_concentrations (mM).
* STATES[19] is Ki in component ionic_concentrations (mM).
* STATES[20] is y in component L_type_Ca_channel_y_gate (dimensionless).
* ALGEBRAIC[4] is y_infinity in component L_type_Ca_channel_y_gate (dimensionless).
* ALGEBRAIC[12] is tau_y in component L_type_Ca_channel_y_gate (ms).
* ALGEBRAIC[25] is g_K in component time_dependent_potassium_current (mS_per_mm2).
* CONSTANTS[24] is g_K_max in component time_dependent_potassium_current (mS_per_mm2).
* ALGEBRAIC[26] is E_K in component time_dependent_potassium_current (mV).
* CONSTANTS[25] is P_NaK in component time_dependent_potassium_current (dimensionless).
* STATES[21] is X in component time_dependent_potassium_current_X_gate (dimensionless).
* ALGEBRAIC[27] is Xi in component time_dependent_potassium_current_Xi_gate (dimensionless).
* ALGEBRAIC[5] is alpha_X in component time_dependent_potassium_current_X_gate (per_ms).
* ALGEBRAIC[13] is beta_X in component time_dependent_potassium_current_X_gate (per_ms).
* ALGEBRAIC[30] is E_K1 in component time_independent_potassium_current (mV).
* ALGEBRAIC[29] is g_K1 in component time_independent_potassium_current (mS_per_mm2).
* CONSTANTS[26] is g_K1_max in component time_independent_potassium_current (mS_per_mm2).
* ALGEBRAIC[33] is K1_infinity in component time_independent_potassium_current_K1_gate (dimensionless).
* ALGEBRAIC[31] is alpha_K1 in component time_independent_potassium_current_K1_gate (per_ms).
* ALGEBRAIC[32] is beta_K1 in component time_independent_potassium_current_K1_gate (per_ms).
* ALGEBRAIC[35] is E_Kp in component plateau_potassium_current (mV).
* CONSTANTS[27] is g_Kp in component plateau_potassium_current (mS_per_mm2).
* ALGEBRAIC[36] is Kp in component plateau_potassium_current (dimensionless).
* CONSTANTS[28] is k_NaCa in component Na_Ca_exchanger (uA_per_mm2).
* CONSTANTS[29] is K_mNa in component Na_Ca_exchanger (mM).
* CONSTANTS[30] is K_mCa in component Na_Ca_exchanger (mM).
* CONSTANTS[31] is k_sat in component Na_Ca_exchanger (dimensionless).
* CONSTANTS[32] is eta in component Na_Ca_exchanger (dimensionless).
* STATES[22] is Cai in component calcium_subsystem (mM).
* CONSTANTS[33] is K_mpCa in component sarcolemmal_calcium_pump (mM).
* CONSTANTS[34] is I_pCa in component sarcolemmal_calcium_pump (uA_per_mm2).
* CONSTANTS[35] is g_Nab in component sodium_background_current (mS_per_mm2).
* ALGEBRAIC[40] is E_NaN in component sodium_background_current (mV).
* CONSTANTS[36] is g_Cab in component calcium_background_current (mS_per_mm2).
* ALGEBRAIC[42] is E_CaN in component calcium_background_current (mV).
* CONSTANTS[37] is I_NaK in component sodium_potassium_pump (uA_per_mm2).
* ALGEBRAIC[44] is f_NaK in component sodium_potassium_pump (dimensionless).
* CONSTANTS[38] is K_mNai in component sodium_potassium_pump (mM).
* CONSTANTS[39] is K_mKo in component sodium_potassium_pump (mM).
* CONSTANTS[90] is sigma in component sodium_potassium_pump (dimensionless).
* ALGEBRAIC[55] is i_ns_Na in component non_specific_calcium_activated_current (uA_per_mm2).
* ALGEBRAIC[61] is i_ns_K in component non_specific_calcium_activated_current (uA_per_mm2).
* ALGEBRAIC[52] is I_ns_Na in component non_specific_calcium_activated_current (uA_per_mm2).
* ALGEBRAIC[59] is I_ns_K in component non_specific_calcium_activated_current (uA_per_mm2).
* CONSTANTS[40] is K_m_ns_Ca in component non_specific_calcium_activated_current (mM).
* CONSTANTS[41] is P_ns_Ca in component non_specific_calcium_activated_current (mm_per_ms).
* ALGEBRAIC[48] is EnsCa in component non_specific_calcium_activated_current (mV).
* ALGEBRAIC[50] is VnsCa in component non_specific_calcium_activated_current (mV).
* CONSTANTS[42] is Am in component calcium_subsystem (per_mm).
* CONSTANTS[43] is V_myo in component calcium_subsystem (dimensionless).
* ALGEBRAIC[45] is RyR_open in component calcium_subsystem (dimensionless).
* STATES[23] is P_O1 in component calcium_subsystem (dimensionless).
* STATES[24] is P_O2 in component calcium_subsystem (dimensionless).
* STATES[25] is P_C1 in component calcium_subsystem (dimensionless).
* STATES[26] is P_C2 in component calcium_subsystem (dimensionless).
* CONSTANTS[44] is v1 in component calcium_subsystem (per_ms).
* CONSTANTS[45] is v2 in component calcium_subsystem (per_ms).
* CONSTANTS[46] is v3 in component calcium_subsystem (mM_per_ms).
* CONSTANTS[47] is nCa in component calcium_subsystem (dimensionless).
* CONSTANTS[48] is mCa in component calcium_subsystem (dimensionless).
* CONSTANTS[49] is k_a_plus in component calcium_subsystem (per_mM4_per_ms).
* CONSTANTS[50] is k_a_minus in component calcium_subsystem (per_ms).
* CONSTANTS[51] is k_b_plus in component calcium_subsystem (per_mM3_per_ms).
* CONSTANTS[52] is k_b_minus in component calcium_subsystem (per_ms).
* CONSTANTS[53] is k_c_plus in component calcium_subsystem (per_ms).
* CONSTANTS[54] is k_c_minus in component calcium_subsystem (per_ms).
* CONSTANTS[55] is k_htrpn_plus in component calcium_subsystem (per_mM_per_ms).
* CONSTANTS[56] is k_htrpn_minus in component calcium_subsystem (per_ms).
* CONSTANTS[57] is k_ltrpn_plus in component calcium_subsystem (per_mM_per_ms).
* CONSTANTS[58] is k_ltrpn_minus in component calcium_subsystem (per_ms).
* CONSTANTS[59] is tau_tr in component calcium_subsystem (ms).
* STATES[27] is Ca_JSR in component calcium_subsystem (mM).
* STATES[28] is Ca_NSR in component calcium_subsystem (mM).
* CONSTANTS[96] is V_JSR in component calcium_subsystem (dimensionless).
* CONSTANTS[94] is V_NSR in component calcium_subsystem (dimensionless).
* CONSTANTS[91] is V_SS in component calcium_subsystem (dimensionless).
* CONSTANTS[60] is K_mup in component calcium_subsystem (mM).
* CONSTANTS[61] is K_mCMDN in component calcium_subsystem (mM).
* CONSTANTS[62] is K_mCSQN in component calcium_subsystem (mM).
* CONSTANTS[63] is tau_xfer in component calcium_subsystem (ms).
* CONSTANTS[64] is HTRPN_tot in component calcium_subsystem (mM).
* CONSTANTS[65] is LTRPN_tot in component calcium_subsystem (mM).
* STATES[29] is HTRPNCa in component calcium_subsystem (mM).
* STATES[30] is LTRPNCa in component calcium_subsystem (mM).
* CONSTANTS[66] is CSQN_tot in component calcium_subsystem (mM).
* CONSTANTS[67] is CMDN_tot in component calcium_subsystem (mM).
* ALGEBRAIC[64] is Bi in component calcium_subsystem (dimensionless).
* ALGEBRAIC[56] is B_SS in component calcium_subsystem (dimensionless).
* ALGEBRAIC[57] is B_JSR in component calcium_subsystem (dimensionless).
* ALGEBRAIC[47] is J_rel in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[49] is J_leak in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[51] is J_up in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[53] is J_tr in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[54] is J_xfer in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[62] is J_trpn in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[58] is J_htrpn in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[60] is J_ltrpn in component calcium_subsystem (mM_per_ms).
* CONSTANTS[68] is Tref in component Myofilaments (N_per_mm2).
* CONSTANTS[69] is beta0 in component Myofilaments (dimensionless).
* CONSTANTS[70] is a in component Myofilaments (dimensionless).
* STATES[31] is Q1 in component Myofilaments (dimensionless).
* STATES[32] is Q2 in component Myofilaments (dimensionless).
* STATES[33] is Q3 in component Myofilaments (dimensionless).
* CONSTANTS[71] is A1 in component Myofilaments (dimensionless).
* CONSTANTS[72] is A2 in component Myofilaments (dimensionless).
* CONSTANTS[73] is A3 in component Myofilaments (dimensionless).
* CONSTANTS[74] is alpha1 in component Myofilaments (dimensionless).
* CONSTANTS[75] is alpha2 in component Myofilaments (dimensionless).
* CONSTANTS[76] is alpha3 in component Myofilaments (dimensionless).
* CONSTANTS[77] is Ca50ref in component Myofilaments (mM).
* CONSTANTS[78] is zp in component Myofilaments (dimensionless).
* CONSTANTS[79] is beta1 in component Myofilaments (dimensionless).
* CONSTANTS[80] is alpha0 in component Myofilaments (per_ms).
* CONSTANTS[81] is alphar1 in component Myofilaments (per_ms).
* CONSTANTS[82] is alphar2 in component Myofilaments (per_ms).
* CONSTANTS[83] is nRel in component Myofilaments (dimensionless).
* CONSTANTS[84] is Kz in component Myofilaments (dimensionless).
* CONSTANTS[85] is nHill in component Myofilaments (dimensionless).
* CONSTANTS[86] is kon in component Myofilaments (per_mM_per_ms).
* CONSTANTS[87] is koff in component Myofilaments (per_ms).
* CONSTANTS[88] is gamma_trpn in component Myofilaments (dimensionless).
* CONSTANTS[89] is TRPN_tot in component Myofilaments (mM).
* ALGEBRAIC[7] is T0 in component Myofilaments (N_per_mm2).
* CONSTANTS[92] is T0max in component Myofilaments (N_per_mm2).
* STATES[34] is z in component Myofilaments (dimensionless).
* CONSTANTS[99] is z_max in component Myofilaments (dimensionless).
* ALGEBRAIC[15] is Q in component Myofilaments (dimensionless).
* STATES[35] is Cab in component Myofilaments (mM).
* CONSTANTS[97] is Ca50 in component Myofilaments (mM).
* CONSTANTS[98] is CaTRPN50 in component Myofilaments (mM).
* CONSTANTS[93] is K_2 in component Myofilaments (dimensionless).
* CONSTANTS[95] is K_1 in component Myofilaments (dimensionless).
* ALGEBRAIC[18] is Tension in component Myofilaments (N_per_mm2).
* RATES[0] is d/dt V in component membrane (mV).
* 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[7] is d/dt C0 in component L_type_Ca_channel (dimensionless).
* RATES[8] is d/dt C1 in component L_type_Ca_channel (dimensionless).
* RATES[9] is d/dt C2 in component L_type_Ca_channel (dimensionless).
* RATES[10] is d/dt C3 in component L_type_Ca_channel (dimensionless).
* RATES[11] is d/dt C4 in component L_type_Ca_channel (dimensionless).
* RATES[5] is d/dt O in component L_type_Ca_channel (dimensionless).
* RATES[12] is d/dt C_Ca0 in component L_type_Ca_channel (dimensionless).
* RATES[13] is d/dt C_Ca1 in component L_type_Ca_channel (dimensionless).
* RATES[14] is d/dt C_Ca2 in component L_type_Ca_channel (dimensionless).
* RATES[15] is d/dt C_Ca3 in component L_type_Ca_channel (dimensionless).
* RATES[16] is d/dt C_Ca4 in component L_type_Ca_channel (dimensionless).
* RATES[6] is d/dt O_Ca in component L_type_Ca_channel (dimensionless).
* RATES[20] is d/dt y in component L_type_Ca_channel_y_gate (dimensionless).
* RATES[21] is d/dt X in component time_dependent_potassium_current_X_gate (dimensionless).
* RATES[25] is d/dt P_C1 in component calcium_subsystem (dimensionless).
* RATES[23] is d/dt P_O1 in component calcium_subsystem (dimensionless).
* RATES[24] is d/dt P_O2 in component calcium_subsystem (dimensionless).
* RATES[26] is d/dt P_C2 in component calcium_subsystem (dimensionless).
* RATES[29] is d/dt HTRPNCa in component calcium_subsystem (mM).
* RATES[30] is d/dt LTRPNCa in component calcium_subsystem (mM).
* RATES[22] is d/dt Cai in component calcium_subsystem (mM).
* RATES[17] is d/dt Ca_SS in component calcium_subsystem (mM).
* RATES[27] is d/dt Ca_JSR in component calcium_subsystem (mM).
* RATES[28] is d/dt Ca_NSR in component calcium_subsystem (mM).
* RATES[1] is d/dt Nai in component ionic_concentrations (mM).
* RATES[19] is d/dt Ki in component ionic_concentrations (mM).
* RATES[18] is d/dt Ko in component ionic_concentrations (mM).
* RATES[35] is d/dt Cab in component Myofilaments (mM).
* RATES[34] is d/dt z in component Myofilaments (dimensionless).
* RATES[31] is d/dt Q1 in component Myofilaments (dimensionless).
* RATES[32] is d/dt Q2 in component Myofilaments (dimensionless).
* RATES[33] is d/dt Q3 in component Myofilaments (dimensionless).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1;
CONSTANTS[1] = 0;
STATES[0] = -84.1638;
CONSTANTS[2] = 8.3145e3;
CONSTANTS[3] = 310;
CONSTANTS[4] = 9.6845e4;
CONSTANTS[5] = 0.01;
CONSTANTS[6] = 100;
CONSTANTS[7] = 9000;
CONSTANTS[8] = 750;
CONSTANTS[9] = 1;
CONSTANTS[10] = -100;
CONSTANTS[11] = 0.128;
STATES[1] = 10.2042;
CONSTANTS[12] = 140;
STATES[2] = 0.0328302;
STATES[3] = 0.988354;
STATES[4] = 0.99254;
CONSTANTS[13] = 33.75e-6;
CONSTANTS[14] = 1e-9;
CONSTANTS[15] = -4.58e-3;
STATES[5] = 9.84546e-21;
STATES[6] = 0;
CONSTANTS[16] = 2;
CONSTANTS[17] = 2;
CONSTANTS[18] = 2;
CONSTANTS[19] = 0.3;
CONSTANTS[20] = 0;
CONSTANTS[21] = 0;
CONSTANTS[22] = 0.01;
STATES[7] = 0.997208;
STATES[8] = 6.38897e-5;
STATES[9] = 1.535e-9;
STATES[10] = 1.63909e-14;
STATES[11] = 6.56337e-20;
STATES[12] = 2.72826e-3;
STATES[13] = 6.99215e-7;
STATES[14] = 6.71989e-11;
STATES[15] = 2.87031e-15;
STATES[16] = 4.59752e-20;
STATES[17] = 1.36058e-4;
CONSTANTS[23] = 1.8;
STATES[18] = 5.4;
STATES[19] = 143.727;
STATES[20] = 0.998983;
CONSTANTS[24] = 0.001128;
CONSTANTS[25] = 0.01833;
STATES[21] = 0.000928836;
CONSTANTS[26] = 7.5e-3;
CONSTANTS[27] = 8.28e-5;
CONSTANTS[28] = 50;
CONSTANTS[29] = 87.5;
CONSTANTS[30] = 1.38;
CONSTANTS[31] = 0.1;
CONSTANTS[32] = 0.35;
STATES[22] = 9.94893e-11;
CONSTANTS[33] = 0.5e-3;
CONSTANTS[34] = 1.15e-2;
CONSTANTS[35] = 1.41e-5;
CONSTANTS[36] = 6.032e-5;
CONSTANTS[37] = 0.013;
CONSTANTS[38] = 10;
CONSTANTS[39] = 1.5;
CONSTANTS[40] = 1.2e-3;
CONSTANTS[41] = 1.75e-9;
CONSTANTS[42] = 546.69;
CONSTANTS[43] = 0.92;
STATES[23] = 1.19168e-3;
STATES[24] = 6.30613e-9;
STATES[25] = 0.762527;
STATES[26] = 0.236283;
CONSTANTS[44] = 1.8;
CONSTANTS[45] = 0.58e-4;
CONSTANTS[46] = 1.8e-3;
CONSTANTS[47] = 4;
CONSTANTS[48] = 3;
CONSTANTS[49] = 1.215e10;
CONSTANTS[50] = 0.1425;
CONSTANTS[51] = 4.05e7;
CONSTANTS[52] = 1.93;
CONSTANTS[53] = 0.018;
CONSTANTS[54] = 0.0008;
CONSTANTS[55] = 20;
CONSTANTS[56] = 0.066e-3;
CONSTANTS[57] = 40;
CONSTANTS[58] = 0.04;
CONSTANTS[59] = 34.48;
STATES[27] = 1.17504;
STATES[28] = 1.243891;
CONSTANTS[60] = 0.5e-3;
CONSTANTS[61] = 2.38e-3;
CONSTANTS[62] = 0.8;
CONSTANTS[63] = 3.125;
CONSTANTS[64] = 0.14;
CONSTANTS[65] = 0.07;
STATES[29] = 0.13598;
STATES[30] = 0.00635;
CONSTANTS[66] = 15;
CONSTANTS[67] = 0.05;
CONSTANTS[68] = 56.2;
CONSTANTS[69] = 4.9;
CONSTANTS[70] = 0.35;
STATES[31] = 0;
STATES[32] = 0;
STATES[33] = 0;
CONSTANTS[71] = -29;
CONSTANTS[72] = 138;
CONSTANTS[73] = 129;
CONSTANTS[74] = 0.03;
CONSTANTS[75] = 0.13;
CONSTANTS[76] = 0.625;
CONSTANTS[77] = 1.05e-3;
CONSTANTS[78] = 0.85;
CONSTANTS[79] = -4;
CONSTANTS[80] = 8e-3;
CONSTANTS[81] = 2e-3;
CONSTANTS[82] = 1.75e-3;
CONSTANTS[83] = 3;
CONSTANTS[84] = 0.15;
CONSTANTS[85] = 3;
CONSTANTS[86] = 100;
CONSTANTS[87] = 0.2;
CONSTANTS[88] = 2;
CONSTANTS[89] = 0.07;
STATES[34] = 0;
STATES[35] = 0;
CONSTANTS[90] = (1.00000/7.00000)*(exp(CONSTANTS[12]/67.3000) - 1.00000);
CONSTANTS[91] = 5.82800e-05*CONSTANTS[43];
CONSTANTS[92] = CONSTANTS[68]*(1.00000+ CONSTANTS[69]*(CONSTANTS[0] - 1.00000));
CONSTANTS[93] = (( 1.00000*CONSTANTS[82]*pow(CONSTANTS[78], CONSTANTS[83]))/(pow(CONSTANTS[78], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83])))*(1.00000 - ( CONSTANTS[83]*pow(CONSTANTS[84], CONSTANTS[83]))/(pow(CONSTANTS[78], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83])));
CONSTANTS[94] = 0.0810000*CONSTANTS[43];
CONSTANTS[95] = ( 1.00000*CONSTANTS[82]*pow(CONSTANTS[78], CONSTANTS[83] - 1.00000)*CONSTANTS[83]*pow(CONSTANTS[84], CONSTANTS[83]))/pow(pow(CONSTANTS[78], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83]), 2.00000);
CONSTANTS[96] = 0.00464000*CONSTANTS[43];
CONSTANTS[97] = CONSTANTS[77]*(1.00000+ CONSTANTS[79]*(CONSTANTS[0] - 1.00000));
CONSTANTS[98] = ( CONSTANTS[97]*CONSTANTS[89])/(CONSTANTS[97]+ (CONSTANTS[87]/CONSTANTS[86])*(1.00000 - ( (1.00000+ CONSTANTS[69]*(CONSTANTS[0] - 1.00000))*0.500000)/CONSTANTS[88]));
CONSTANTS[99] = ( CONSTANTS[80]*1.00000)/pow(CONSTANTS[98]/CONSTANTS[89], CONSTANTS[85]) - ( CONSTANTS[93]*1.00000)/(CONSTANTS[81]+ CONSTANTS[95]*1.00000+CONSTANTS[80]/pow(CONSTANTS[98]/CONSTANTS[89], CONSTANTS[85]));
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[5] = CONSTANTS[19]*STATES[11] - CONSTANTS[18]*STATES[5];
RATES[6] = CONSTANTS[21]*STATES[16] - CONSTANTS[20]*STATES[6];
RATES[25] = - CONSTANTS[49]*pow(STATES[17], CONSTANTS[47])*STATES[25]+ CONSTANTS[50]*STATES[23];
RATES[23] = ( CONSTANTS[49]*pow(STATES[17], CONSTANTS[47])*STATES[25] - ( CONSTANTS[50]*STATES[23]+ CONSTANTS[51]*pow(STATES[17], CONSTANTS[48])*STATES[23]+ CONSTANTS[53]*STATES[23]))+ CONSTANTS[52]*STATES[24]+ CONSTANTS[54]*STATES[26];
RATES[24] = CONSTANTS[51]*pow(STATES[17], CONSTANTS[48])*STATES[23] - CONSTANTS[52]*STATES[24];
RATES[26] = CONSTANTS[53]*STATES[23] - CONSTANTS[54]*STATES[26];
RATES[34] = CONSTANTS[80]*pow(STATES[35]/CONSTANTS[98], CONSTANTS[85])*(1.00000 - STATES[34])+ - STATES[34]*CONSTANTS[81]+( - CONSTANTS[82]*pow(STATES[34], CONSTANTS[83]))/(pow(STATES[34], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83]));
RATES[31] = CONSTANTS[71]*CONSTANTS[1] - 1.00000*CONSTANTS[74]*STATES[31];
RATES[32] = CONSTANTS[72]*CONSTANTS[1] - 1.00000*CONSTANTS[75]*STATES[32];
RATES[33] = CONSTANTS[73]*CONSTANTS[1] - 1.00000*CONSTANTS[76]*STATES[33];
ALGEBRAIC[0] = ( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300)));
ALGEBRAIC[8] = 0.0800000*exp(- STATES[0]/11.0000);
RATES[2] = ALGEBRAIC[0]*(1.00000 - STATES[2]) - ALGEBRAIC[8]*STATES[2];
ALGEBRAIC[1] = (STATES[0]<- 40.0000 ? 0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[9] = (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[1]*(1.00000 - STATES[3]) - ALGEBRAIC[9]*STATES[3];
ALGEBRAIC[2] = (STATES[0]<- 40.0000 ? ( ( - 127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( - 0.0439100*STATES[0]))*(STATES[0]+37.7800))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))) : 0.00000);
ALGEBRAIC[10] = (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[2]*(1.00000 - STATES[4]) - ALGEBRAIC[10]*STATES[4];
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+55.0000)/7.50000))+0.100000/(1.00000+exp((- STATES[0]+21.0000)/6.00000));
ALGEBRAIC[12] = 20.0000+600.000/(1.00000+exp((STATES[0]+30.0000)/9.50000));
RATES[20] = (ALGEBRAIC[4] - STATES[20])/ALGEBRAIC[12];
ALGEBRAIC[5] = ( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)));
ALGEBRAIC[13] = ( 0.000131000*(STATES[0]+30.0000))/(- 1.00000+exp( 0.0687000*(STATES[0]+30.0000)));
RATES[21] = ALGEBRAIC[5]*(1.00000 - STATES[21]) - ALGEBRAIC[13]*STATES[21];
ALGEBRAIC[7] = ( CONSTANTS[92]*STATES[34])/CONSTANTS[99];
ALGEBRAIC[15] = STATES[31]+STATES[32]+STATES[33];
ALGEBRAIC[18] = (ALGEBRAIC[15]<0.00000 ? ( ALGEBRAIC[7]*( CONSTANTS[70]*ALGEBRAIC[15]+1.00000))/(1.00000 - ALGEBRAIC[15]) : ( ALGEBRAIC[7]*(1.00000 - (CONSTANTS[70]+2.00000)*ALGEBRAIC[15]))/(1.00000+ALGEBRAIC[15]));
RATES[35] = CONSTANTS[86]*STATES[22]*(CONSTANTS[89] - STATES[35]) - CONSTANTS[87]*(1.00000 - ALGEBRAIC[18]/( CONSTANTS[88]*CONSTANTS[68]))*STATES[35];
ALGEBRAIC[3] = 0.400000*exp((STATES[0]+12.0000)/10.0000);
ALGEBRAIC[11] = 0.0500000*exp((STATES[0]+12.0000)/- 13.0000);
ALGEBRAIC[21] = ( 0.187500*STATES[17])/1.00000;
RATES[7] = ( ALGEBRAIC[11]*STATES[8]+ CONSTANTS[22]*STATES[12]) - ( 4.00000*ALGEBRAIC[3]+ALGEBRAIC[21])*STATES[7];
RATES[8] = ( 4.00000*ALGEBRAIC[3]*STATES[7]+ 2.00000*ALGEBRAIC[11]*STATES[9]+ (CONSTANTS[22]/CONSTANTS[17])*STATES[13]) - (ALGEBRAIC[11]+ 3.00000*ALGEBRAIC[3]+ ALGEBRAIC[21]*CONSTANTS[16])*STATES[8];
RATES[9] = ( 3.00000*ALGEBRAIC[3]*STATES[8]+ 3.00000*ALGEBRAIC[11]*STATES[10]+ (CONSTANTS[22]/pow(CONSTANTS[17], 2.00000))*STATES[14]) - ( ALGEBRAIC[11]*2.00000+ 2.00000*ALGEBRAIC[3]+ ALGEBRAIC[21]*pow(CONSTANTS[16], 2.00000))*STATES[9];
RATES[10] = ( 2.00000*ALGEBRAIC[3]*STATES[9]+ 4.00000*ALGEBRAIC[11]*STATES[11]+ (CONSTANTS[22]/pow(CONSTANTS[17], 3.00000))*STATES[15]) - ( ALGEBRAIC[11]*3.00000+ALGEBRAIC[3]+ ALGEBRAIC[21]*pow(CONSTANTS[16], 3.00000))*STATES[10];
RATES[11] = ( ALGEBRAIC[3]*STATES[10]+ CONSTANTS[18]*STATES[5]+ (CONSTANTS[22]/pow(CONSTANTS[17], 4.00000))*STATES[16]) - ( ALGEBRAIC[11]*4.00000+CONSTANTS[19]+ ALGEBRAIC[21]*pow(CONSTANTS[16], 4.00000))*STATES[11];
ALGEBRAIC[16] = ALGEBRAIC[3]*CONSTANTS[16];
ALGEBRAIC[19] = ALGEBRAIC[11]/CONSTANTS[17];
RATES[12] = ( ALGEBRAIC[19]*STATES[13]+ ALGEBRAIC[21]*STATES[7]) - ( 4.00000*ALGEBRAIC[16]+CONSTANTS[22])*STATES[12];
RATES[13] = ( 4.00000*ALGEBRAIC[16]*STATES[12]+ 2.00000*ALGEBRAIC[19]*STATES[14]+ ALGEBRAIC[21]*CONSTANTS[16]*STATES[8]) - (ALGEBRAIC[19]+ 3.00000*ALGEBRAIC[16]+CONSTANTS[22]/CONSTANTS[17])*STATES[13];
RATES[14] = ( 3.00000*ALGEBRAIC[16]*STATES[13]+ 3.00000*ALGEBRAIC[19]*STATES[15]+ ALGEBRAIC[21]*pow(CONSTANTS[16], 2.00000)*STATES[9]) - ( ALGEBRAIC[19]*2.00000+ 2.00000*ALGEBRAIC[16]+CONSTANTS[22]/pow(CONSTANTS[17], 2.00000))*STATES[14];
RATES[15] = ( 2.00000*ALGEBRAIC[16]*STATES[14]+ 4.00000*ALGEBRAIC[19]*STATES[16]+ ALGEBRAIC[21]*pow(CONSTANTS[16], 3.00000)*STATES[10]) - ( ALGEBRAIC[19]*3.00000+ALGEBRAIC[16]+CONSTANTS[22]/pow(CONSTANTS[17], 3.00000))*STATES[15];
RATES[16] = ( ALGEBRAIC[16]*STATES[15]+ CONSTANTS[20]*STATES[6]+ ALGEBRAIC[21]*pow(CONSTANTS[16], 4.00000)*STATES[11]) - ( ALGEBRAIC[19]*4.00000+CONSTANTS[21]+CONSTANTS[22]/pow(CONSTANTS[17], 4.00000))*STATES[16];
ALGEBRAIC[49] = CONSTANTS[45]*(STATES[28] - STATES[22]);
ALGEBRAIC[51] = ( CONSTANTS[46]*pow(STATES[22], 2.00000))/(pow(CONSTANTS[60], 2.00000)+pow(STATES[22], 2.00000));
ALGEBRAIC[53] = (STATES[28] - STATES[27])/CONSTANTS[59];
RATES[28] = ( (ALGEBRAIC[51] - ALGEBRAIC[49])*CONSTANTS[43])/CONSTANTS[94] - ( ALGEBRAIC[53]*CONSTANTS[96])/CONSTANTS[94];
ALGEBRAIC[58] = CONSTANTS[55]*STATES[22]*(CONSTANTS[64] - STATES[29]) - CONSTANTS[56]*STATES[29];
RATES[29] = ALGEBRAIC[58];
ALGEBRAIC[20] = ( (( CONSTANTS[13]*4.00000*STATES[0]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.00100000*exp(( 2.00000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 0.341000*CONSTANTS[23]))/(exp(( 2.00000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[22] = ALGEBRAIC[20]*STATES[20]*(STATES[5]+STATES[6]);
ALGEBRAIC[56] = 1.00000/(1.00000+( CONSTANTS[67]*CONSTANTS[61])/pow(CONSTANTS[61]+STATES[17], 2.00000));
ALGEBRAIC[45] = STATES[23]+STATES[24];
ALGEBRAIC[47] = CONSTANTS[44]*ALGEBRAIC[45]*(STATES[27] - STATES[17]);
ALGEBRAIC[54] = (STATES[17] - STATES[22])/CONSTANTS[63];
RATES[17] = ALGEBRAIC[56]*((( ALGEBRAIC[47]*CONSTANTS[96])/CONSTANTS[91] - ( ALGEBRAIC[54]*CONSTANTS[43])/CONSTANTS[91]) - ( ALGEBRAIC[22]*CONSTANTS[42])/( 2.00000*CONSTANTS[91]*CONSTANTS[4]));
ALGEBRAIC[57] = 1.00000/(1.00000+( CONSTANTS[66]*CONSTANTS[62])/pow(CONSTANTS[62]+STATES[27], 2.00000));
RATES[27] = ALGEBRAIC[57]*(ALGEBRAIC[53] - ALGEBRAIC[47]);
ALGEBRAIC[14] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log(CONSTANTS[12]/STATES[1]);
ALGEBRAIC[17] = CONSTANTS[11]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*(STATES[0] - ALGEBRAIC[14]);
ALGEBRAIC[38] = (( (( (( CONSTANTS[28]*1.00000)/(pow(CONSTANTS[29], 3.00000)+pow(CONSTANTS[12], 3.00000)))*1.00000)/(CONSTANTS[30]+CONSTANTS[23]))*1.00000)/(1.00000+ CONSTANTS[31]*exp(( (CONSTANTS[32] - 1.00000)*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))))*( exp(( CONSTANTS[32]*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))*pow(STATES[1], 3.00000)*CONSTANTS[23] - exp(( (CONSTANTS[32] - 1.00000)*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))*pow(CONSTANTS[12], 3.00000)*STATES[22]);
ALGEBRAIC[40] = ALGEBRAIC[14];
ALGEBRAIC[41] = CONSTANTS[35]*(STATES[0] - ALGEBRAIC[40]);
ALGEBRAIC[44] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))+ 0.0365000*CONSTANTS[90]*exp(( - STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])));
ALGEBRAIC[46] = ( (( CONSTANTS[37]*ALGEBRAIC[44]*1.00000)/(1.00000+pow(CONSTANTS[38]/STATES[1], 1.50000)))*STATES[18])/(STATES[18]+CONSTANTS[39]);
ALGEBRAIC[48] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log((STATES[18]+CONSTANTS[12])/(STATES[19]+STATES[1]));
ALGEBRAIC[50] = STATES[0] - ALGEBRAIC[48];
ALGEBRAIC[52] = ( (( CONSTANTS[41]*pow(1.00000, 2.00000)*ALGEBRAIC[50]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.750000*STATES[1]*exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 0.750000*CONSTANTS[12]))/(exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[55] = ( ALGEBRAIC[52]*1.00000)/(1.00000+pow(CONSTANTS[40]/STATES[22], 3.00000));
RATES[1] = ( - (ALGEBRAIC[17]+ALGEBRAIC[41]+ALGEBRAIC[55]+ ALGEBRAIC[38]*3.00000+ ALGEBRAIC[46]*3.00000)*CONSTANTS[42])/( CONSTANTS[43]*CONSTANTS[4]);
ALGEBRAIC[60] = CONSTANTS[57]*STATES[22]*(CONSTANTS[65] - STATES[30]) - CONSTANTS[58]*STATES[30];
RATES[30] = ALGEBRAIC[60];
ALGEBRAIC[23] = CONSTANTS[14]/(1.00000+ALGEBRAIC[20]/CONSTANTS[15]);
ALGEBRAIC[24] = ( (( ALGEBRAIC[23]*STATES[20]*(STATES[5]+STATES[6])*STATES[0]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( STATES[19]*exp(( STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - STATES[18]))/(exp(( STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[25] = CONSTANTS[24]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[26] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log((STATES[18]+ CONSTANTS[25]*CONSTANTS[12])/(STATES[19]+ CONSTANTS[25]*STATES[1]));
ALGEBRAIC[27] = 1.00000/(1.00000+exp((STATES[0] - 56.2600)/32.1000));
ALGEBRAIC[28] = ALGEBRAIC[25]*ALGEBRAIC[27]*pow(STATES[21], 2.00000)*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[30] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log(STATES[18]/STATES[19]);
ALGEBRAIC[29] = CONSTANTS[26]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[31] = 1.02000/(1.00000+exp( 0.238500*((STATES[0] - ALGEBRAIC[30]) - 59.2150)));
ALGEBRAIC[32] = ( 0.491240*(exp( 0.0803200*((STATES[0]+5.47600) - ALGEBRAIC[30]))+exp( 0.0617500*(STATES[0] - (ALGEBRAIC[30]+594.310)))))/(1.00000+exp( - 0.514300*((STATES[0] - ALGEBRAIC[30])+4.75300)));
ALGEBRAIC[33] = ALGEBRAIC[31]/(ALGEBRAIC[31]+ALGEBRAIC[32]);
ALGEBRAIC[34] = ALGEBRAIC[29]*ALGEBRAIC[33]*(STATES[0] - ALGEBRAIC[30]);
ALGEBRAIC[35] = ALGEBRAIC[30];
ALGEBRAIC[36] = 1.00000/(1.00000+exp((7.48800 - STATES[0])/5.98000));
ALGEBRAIC[37] = CONSTANTS[27]*ALGEBRAIC[36]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[59] = ( (( CONSTANTS[41]*pow(1.00000, 2.00000)*ALGEBRAIC[50]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.750000*STATES[19]*exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 0.750000*STATES[18]))/(exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[61] = ( ALGEBRAIC[59]*1.00000)/(1.00000+pow(CONSTANTS[40]/STATES[22], 3.00000));
RATES[19] = ( - (ALGEBRAIC[24]+ALGEBRAIC[28]+ALGEBRAIC[34]+ALGEBRAIC[37]+ALGEBRAIC[61]+ - ALGEBRAIC[46]*2.00000)*CONSTANTS[42])/( CONSTANTS[43]*CONSTANTS[4]);
RATES[18] = ( (ALGEBRAIC[24]+ALGEBRAIC[28]+ALGEBRAIC[34]+ALGEBRAIC[37]+ALGEBRAIC[61]+ - ALGEBRAIC[46]*2.00000)*CONSTANTS[42])/( CONSTANTS[43]*CONSTANTS[4]);
ALGEBRAIC[39] = ( CONSTANTS[34]*STATES[22])/(CONSTANTS[33]+STATES[22]);
ALGEBRAIC[42] = (( CONSTANTS[2]*CONSTANTS[3])/( 2.00000*CONSTANTS[4]))*log(CONSTANTS[23]/STATES[22]);
ALGEBRAIC[43] = CONSTANTS[36]*(STATES[0] - ALGEBRAIC[42]);
ALGEBRAIC[63] = ALGEBRAIC[55]+ALGEBRAIC[61];
ALGEBRAIC[6] = (VOI>=CONSTANTS[6]&&VOI<=CONSTANTS[7]&&(VOI - CONSTANTS[6]) - floor((VOI - CONSTANTS[6])/CONSTANTS[8])*CONSTANTS[8]<=CONSTANTS[9] ? CONSTANTS[10] : 0.00000);
RATES[0] = (ALGEBRAIC[17]+ALGEBRAIC[22]+ALGEBRAIC[24]+ALGEBRAIC[28]+ALGEBRAIC[38]+ALGEBRAIC[34]+ALGEBRAIC[37]+ALGEBRAIC[39]+ALGEBRAIC[41]+ALGEBRAIC[43]+ALGEBRAIC[46]+ALGEBRAIC[63]+ALGEBRAIC[6])/CONSTANTS[5];
ALGEBRAIC[64] = 1.00000/(1.00000+( CONSTANTS[67]*CONSTANTS[61])/pow(CONSTANTS[61]+STATES[22], 2.00000));
ALGEBRAIC[62] = ALGEBRAIC[58]+ALGEBRAIC[60];
RATES[22] = ALGEBRAIC[64]*((ALGEBRAIC[49]+ALGEBRAIC[54]) - (ALGEBRAIC[51]+ALGEBRAIC[62]+( ((ALGEBRAIC[43] - 2.00000*ALGEBRAIC[38])+ALGEBRAIC[39])*CONSTANTS[42])/( 2.00000*CONSTANTS[43]*CONSTANTS[4])));
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = ( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300)));
ALGEBRAIC[8] = 0.0800000*exp(- STATES[0]/11.0000);
ALGEBRAIC[1] = (STATES[0]<- 40.0000 ? 0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[9] = (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[2] = (STATES[0]<- 40.0000 ? ( ( - 127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( - 0.0439100*STATES[0]))*(STATES[0]+37.7800))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))) : 0.00000);
ALGEBRAIC[10] = (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[4] = 1.00000/(1.00000+exp((STATES[0]+55.0000)/7.50000))+0.100000/(1.00000+exp((- STATES[0]+21.0000)/6.00000));
ALGEBRAIC[12] = 20.0000+600.000/(1.00000+exp((STATES[0]+30.0000)/9.50000));
ALGEBRAIC[5] = ( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)));
ALGEBRAIC[13] = ( 0.000131000*(STATES[0]+30.0000))/(- 1.00000+exp( 0.0687000*(STATES[0]+30.0000)));
ALGEBRAIC[7] = ( CONSTANTS[92]*STATES[34])/CONSTANTS[99];
ALGEBRAIC[15] = STATES[31]+STATES[32]+STATES[33];
ALGEBRAIC[18] = (ALGEBRAIC[15]<0.00000 ? ( ALGEBRAIC[7]*( CONSTANTS[70]*ALGEBRAIC[15]+1.00000))/(1.00000 - ALGEBRAIC[15]) : ( ALGEBRAIC[7]*(1.00000 - (CONSTANTS[70]+2.00000)*ALGEBRAIC[15]))/(1.00000+ALGEBRAIC[15]));
ALGEBRAIC[3] = 0.400000*exp((STATES[0]+12.0000)/10.0000);
ALGEBRAIC[11] = 0.0500000*exp((STATES[0]+12.0000)/- 13.0000);
ALGEBRAIC[21] = ( 0.187500*STATES[17])/1.00000;
ALGEBRAIC[16] = ALGEBRAIC[3]*CONSTANTS[16];
ALGEBRAIC[19] = ALGEBRAIC[11]/CONSTANTS[17];
ALGEBRAIC[49] = CONSTANTS[45]*(STATES[28] - STATES[22]);
ALGEBRAIC[51] = ( CONSTANTS[46]*pow(STATES[22], 2.00000))/(pow(CONSTANTS[60], 2.00000)+pow(STATES[22], 2.00000));
ALGEBRAIC[53] = (STATES[28] - STATES[27])/CONSTANTS[59];
ALGEBRAIC[58] = CONSTANTS[55]*STATES[22]*(CONSTANTS[64] - STATES[29]) - CONSTANTS[56]*STATES[29];
ALGEBRAIC[20] = ( (( CONSTANTS[13]*4.00000*STATES[0]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.00100000*exp(( 2.00000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 0.341000*CONSTANTS[23]))/(exp(( 2.00000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[22] = ALGEBRAIC[20]*STATES[20]*(STATES[5]+STATES[6]);
ALGEBRAIC[56] = 1.00000/(1.00000+( CONSTANTS[67]*CONSTANTS[61])/pow(CONSTANTS[61]+STATES[17], 2.00000));
ALGEBRAIC[45] = STATES[23]+STATES[24];
ALGEBRAIC[47] = CONSTANTS[44]*ALGEBRAIC[45]*(STATES[27] - STATES[17]);
ALGEBRAIC[54] = (STATES[17] - STATES[22])/CONSTANTS[63];
ALGEBRAIC[57] = 1.00000/(1.00000+( CONSTANTS[66]*CONSTANTS[62])/pow(CONSTANTS[62]+STATES[27], 2.00000));
ALGEBRAIC[14] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log(CONSTANTS[12]/STATES[1]);
ALGEBRAIC[17] = CONSTANTS[11]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*(STATES[0] - ALGEBRAIC[14]);
ALGEBRAIC[38] = (( (( (( CONSTANTS[28]*1.00000)/(pow(CONSTANTS[29], 3.00000)+pow(CONSTANTS[12], 3.00000)))*1.00000)/(CONSTANTS[30]+CONSTANTS[23]))*1.00000)/(1.00000+ CONSTANTS[31]*exp(( (CONSTANTS[32] - 1.00000)*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))))*( exp(( CONSTANTS[32]*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))*pow(STATES[1], 3.00000)*CONSTANTS[23] - exp(( (CONSTANTS[32] - 1.00000)*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))*pow(CONSTANTS[12], 3.00000)*STATES[22]);
ALGEBRAIC[40] = ALGEBRAIC[14];
ALGEBRAIC[41] = CONSTANTS[35]*(STATES[0] - ALGEBRAIC[40]);
ALGEBRAIC[44] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))+ 0.0365000*CONSTANTS[90]*exp(( - STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])));
ALGEBRAIC[46] = ( (( CONSTANTS[37]*ALGEBRAIC[44]*1.00000)/(1.00000+pow(CONSTANTS[38]/STATES[1], 1.50000)))*STATES[18])/(STATES[18]+CONSTANTS[39]);
ALGEBRAIC[48] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log((STATES[18]+CONSTANTS[12])/(STATES[19]+STATES[1]));
ALGEBRAIC[50] = STATES[0] - ALGEBRAIC[48];
ALGEBRAIC[52] = ( (( CONSTANTS[41]*pow(1.00000, 2.00000)*ALGEBRAIC[50]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.750000*STATES[1]*exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 0.750000*CONSTANTS[12]))/(exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[55] = ( ALGEBRAIC[52]*1.00000)/(1.00000+pow(CONSTANTS[40]/STATES[22], 3.00000));
ALGEBRAIC[60] = CONSTANTS[57]*STATES[22]*(CONSTANTS[65] - STATES[30]) - CONSTANTS[58]*STATES[30];
ALGEBRAIC[23] = CONSTANTS[14]/(1.00000+ALGEBRAIC[20]/CONSTANTS[15]);
ALGEBRAIC[24] = ( (( ALGEBRAIC[23]*STATES[20]*(STATES[5]+STATES[6])*STATES[0]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( STATES[19]*exp(( STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - STATES[18]))/(exp(( STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[25] = CONSTANTS[24]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[26] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log((STATES[18]+ CONSTANTS[25]*CONSTANTS[12])/(STATES[19]+ CONSTANTS[25]*STATES[1]));
ALGEBRAIC[27] = 1.00000/(1.00000+exp((STATES[0] - 56.2600)/32.1000));
ALGEBRAIC[28] = ALGEBRAIC[25]*ALGEBRAIC[27]*pow(STATES[21], 2.00000)*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[30] = (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log(STATES[18]/STATES[19]);
ALGEBRAIC[29] = CONSTANTS[26]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[31] = 1.02000/(1.00000+exp( 0.238500*((STATES[0] - ALGEBRAIC[30]) - 59.2150)));
ALGEBRAIC[32] = ( 0.491240*(exp( 0.0803200*((STATES[0]+5.47600) - ALGEBRAIC[30]))+exp( 0.0617500*(STATES[0] - (ALGEBRAIC[30]+594.310)))))/(1.00000+exp( - 0.514300*((STATES[0] - ALGEBRAIC[30])+4.75300)));
ALGEBRAIC[33] = ALGEBRAIC[31]/(ALGEBRAIC[31]+ALGEBRAIC[32]);
ALGEBRAIC[34] = ALGEBRAIC[29]*ALGEBRAIC[33]*(STATES[0] - ALGEBRAIC[30]);
ALGEBRAIC[35] = ALGEBRAIC[30];
ALGEBRAIC[36] = 1.00000/(1.00000+exp((7.48800 - STATES[0])/5.98000));
ALGEBRAIC[37] = CONSTANTS[27]*ALGEBRAIC[36]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[59] = ( (( CONSTANTS[41]*pow(1.00000, 2.00000)*ALGEBRAIC[50]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.750000*STATES[19]*exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 0.750000*STATES[18]))/(exp(( ALGEBRAIC[50]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[61] = ( ALGEBRAIC[59]*1.00000)/(1.00000+pow(CONSTANTS[40]/STATES[22], 3.00000));
ALGEBRAIC[39] = ( CONSTANTS[34]*STATES[22])/(CONSTANTS[33]+STATES[22]);
ALGEBRAIC[42] = (( CONSTANTS[2]*CONSTANTS[3])/( 2.00000*CONSTANTS[4]))*log(CONSTANTS[23]/STATES[22]);
ALGEBRAIC[43] = CONSTANTS[36]*(STATES[0] - ALGEBRAIC[42]);
ALGEBRAIC[63] = ALGEBRAIC[55]+ALGEBRAIC[61];
ALGEBRAIC[6] = (VOI>=CONSTANTS[6]&&VOI<=CONSTANTS[7]&&(VOI - CONSTANTS[6]) - floor((VOI - CONSTANTS[6])/CONSTANTS[8])*CONSTANTS[8]<=CONSTANTS[9] ? CONSTANTS[10] : 0.00000);
ALGEBRAIC[64] = 1.00000/(1.00000+( CONSTANTS[67]*CONSTANTS[61])/pow(CONSTANTS[61]+STATES[22], 2.00000));
ALGEBRAIC[62] = ALGEBRAIC[58]+ALGEBRAIC[60];
}