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 106 entries in the algebraic variable array.
There are a total of 22 entries in each of the rate and state variable arrays.
There are a total of 93 entries in the constant variable array.
*/
/*
* VOI is time in component environment (ms).
* STATES[0] is V in component membrane (mV).
* CONSTANTS[0] is R in component membrane (mJ_per_mole_K).
* CONSTANTS[1] is T in component membrane (kelvin).
* CONSTANTS[2] is F in component membrane (C_per_mole).
* CONSTANTS[3] is Cm in component membrane (uF).
* ALGEBRAIC[0] is I_Stim in component membrane (uA).
* CONSTANTS[4] is stim_period in component membrane (ms).
* CONSTANTS[5] is stim_duration in component membrane (ms).
* CONSTANTS[6] is stim_amplitude in component membrane (uA).
* ALGEBRAIC[25] is i_Na in component sodium_current (nanoA).
* ALGEBRAIC[27] is i_t in component Ca_independent_transient_outward_K_current (nanoA).
* ALGEBRAIC[28] is i_ss in component steady_state_outward_K_current (nanoA).
* ALGEBRAIC[32] is i_f in component hyperpolarisation_activated_current (nanoA).
* ALGEBRAIC[29] is i_K1 in component inward_rectifier (nanoA).
* ALGEBRAIC[33] is i_B_Na in component background_currents (nanoA).
* ALGEBRAIC[35] is i_B_K in component background_currents (nanoA).
* ALGEBRAIC[21] is i_NaK in component sodium_potassium_pump (uA).
* ALGEBRAIC[105] is I_LCC in component convert_hinch (uA).
* ALGEBRAIC[38] is I_NaCa in component convert_hinch (uA).
* ALGEBRAIC[41] is I_pCa in component convert_hinch (uA).
* ALGEBRAIC[44] is I_CaB in component convert_hinch (uA).
* ALGEBRAIC[11] is FVRT in component membrane (dimensionless).
* ALGEBRAIC[20] is FVRT_Ca in component membrane (dimensionless).
* CONSTANTS[7] is V_myo in component cell_geometry (um3).
* CONSTANTS[8] is V_SR in component cell_geometry (um3).
* CONSTANTS[9] is V_myo_uL in component cell_geometry (uL).
* CONSTANTS[10] is V_SR_uL in component cell_geometry (uL).
* ALGEBRAIC[104] is I_LCC in component LCC_current (mM_per_ms).
* ALGEBRAIC[37] is I_NaCa in component Na_Ca_Exchanger (mM_per_ms).
* ALGEBRAIC[40] is I_pCa in component Sarcolemmal_Ca_pump (mM_per_ms).
* ALGEBRAIC[43] is I_CaB in component Background_Ca_current (mM_per_ms).
* ALGEBRAIC[100] is I_RyR in component RyR_current (mM_per_ms).
* ALGEBRAIC[101] is I_RyR in component convert_hinch (mM_per_ms).
* CONSTANTS[11] is i_NaK_max in component sodium_potassium_pump (uA).
* CONSTANTS[12] is K_m_K in component sodium_potassium_pump (mM).
* CONSTANTS[13] is K_m_Na in component sodium_potassium_pump (mM).
* CONSTANTS[82] is sigma in component sodium_potassium_pump (dimensionless).
* CONSTANTS[14] is K_o in component standard_ionic_concentrations (millimolar).
* CONSTANTS[15] is Na_o in component standard_ionic_concentrations (millimolar).
* STATES[1] is Na_i in component intracellular_ion_concentrations (mM).
* STATES[2] is Ca_i in component intracellular_ion_concentrations (mM).
* STATES[3] is Ca_SR in component intracellular_ion_concentrations (mM).
* STATES[4] is K_i in component intracellular_ion_concentrations (mM).
* STATES[5] is TRPN in component intracellular_ion_concentrations (mM).
* ALGEBRAIC[30] is i_f_Na in component hyperpolarisation_activated_current (nanoA).
* ALGEBRAIC[31] is i_f_K in component hyperpolarisation_activated_current (nanoA).
* ALGEBRAIC[39] is I_SERCA in component SERCA (mM_per_ms).
* ALGEBRAIC[45] is I_SR in component SR_Ca_leak_current (mM_per_ms).
* ALGEBRAIC[23] is Ca_b in component intracellular_ion_concentrations (mM).
* CONSTANTS[16] is Ca_TRPN_Max in component troponin (mM).
* ALGEBRAIC[61] is J_TRPN in component troponin (mM_per_ms).
* ALGEBRAIC[46] is beta_CMDN in component calmodulin_Ca_buffer (dimensionless).
* ALGEBRAIC[24] is E_Na in component sodium_current (millivolt).
* CONSTANTS[17] is g_Na in component sodium_current (microS).
* CONSTANTS[83] is g_Na_endo in component sodium_current (microS).
* STATES[6] is m in component sodium_current_m_gate (dimensionless).
* STATES[7] is h in component sodium_current_h_gate (dimensionless).
* STATES[8] is j in component sodium_current_j_gate (dimensionless).
* ALGEBRAIC[26] is E_K in component Ca_independent_transient_outward_K_current (millivolt).
* CONSTANTS[18] is g_t in component Ca_independent_transient_outward_K_current (microS).
* CONSTANTS[84] is g_t_endo in component Ca_independent_transient_outward_K_current (microS).
* CONSTANTS[19] is a_endo in component Ca_independent_transient_outward_K_current (dimensionless).
* CONSTANTS[20] is b_endo in component Ca_independent_transient_outward_K_current (dimensionless).
* STATES[9] is r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
* STATES[10] is s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
* STATES[11] is s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
* CONSTANTS[21] is g_ss in component steady_state_outward_K_current (microS).
* STATES[12] is r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
* STATES[13] is s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
* CONSTANTS[22] is g_K1 in component inward_rectifier (microS).
* CONSTANTS[23] is g_f in component hyperpolarisation_activated_current (microS).
* CONSTANTS[24] is f_Na in component hyperpolarisation_activated_current (dimensionless).
* CONSTANTS[85] is f_K in component hyperpolarisation_activated_current (dimensionless).
* STATES[14] is y in component hyperpolarisation_activated_current_y_gate (dimensionless).
* ALGEBRAIC[36] is i_B in component background_currents (nanoA).
* ALGEBRAIC[34] is i_B_Ca in component background_currents (nanoA).
* CONSTANTS[25] is g_B_Na in component background_currents (microS).
* CONSTANTS[26] is g_B_Ca in component background_currents (microS).
* CONSTANTS[27] is g_B_K in component background_currents (microS).
* CONSTANTS[28] is E_Ca in component background_currents (millivolt).
* CONSTANTS[29] is Ca_o in component standard_ionic_concentrations (millimolar).
* CONSTANTS[30] is g_D in component CaRU (um3_per_ms).
* CONSTANTS[31] is J_R in component CaRU (um3_per_ms).
* CONSTANTS[32] is J_L in component CaRU (um3_per_ms).
* CONSTANTS[33] is N in component CaRU (dimensionless).
* CONSTANTS[34] is K_mNa in component Na_Ca_Exchanger (mM).
* CONSTANTS[35] is K_mCa in component Na_Ca_Exchanger (mM).
* CONSTANTS[36] is eta in component Na_Ca_Exchanger (dimensionless).
* CONSTANTS[37] is k_sat in component Na_Ca_Exchanger (dimensionless).
* CONSTANTS[38] is g_NCX in component Na_Ca_Exchanger (mM_per_ms).
* CONSTANTS[39] is g_SERCA in component SERCA (mM_per_ms).
* CONSTANTS[40] is K_SERCA in component SERCA (mM).
* CONSTANTS[41] is g_pCa in component Sarcolemmal_Ca_pump (mM_per_ms).
* CONSTANTS[42] is K_mpCa in component Sarcolemmal_Ca_pump (mM).
* ALGEBRAIC[42] is E_Ca in component Background_Ca_current (mV).
* CONSTANTS[43] is g_CaB in component Background_Ca_current (mM_per_mV_ms).
* CONSTANTS[44] is g_SRl in component SR_Ca_leak_current (per_ms).
* CONSTANTS[45] is k_CMDN in component calmodulin_Ca_buffer (mM).
* CONSTANTS[46] is B_CMDN in component calmodulin_Ca_buffer (mM).
* STATES[15] is z in component tropomyosin (dimensionless).
* ALGEBRAIC[54] is z_max in component tropomyosin (dimensionless).
* CONSTANTS[47] is k_on in component troponin (per_mM_per_ms).
* CONSTANTS[48] is k_Ref_off in component troponin (per_ms).
* CONSTANTS[49] is gamma_trpn in component troponin (dimensionless).
* CONSTANTS[50] is alpha_0 in component tropomyosin (per_ms).
* CONSTANTS[51] is alpha_r1 in component tropomyosin (per_ms).
* CONSTANTS[52] is alpha_r2 in component tropomyosin (per_ms).
* CONSTANTS[53] is n_Rel in component tropomyosin (dimensionless).
* CONSTANTS[54] is K_z in component tropomyosin (dimensionless).
* CONSTANTS[55] is n_Hill in component tropomyosin (dimensionless).
* CONSTANTS[56] is Ca_50ref in component tropomyosin (mM).
* CONSTANTS[57] is z_p in component tropomyosin (dimensionless).
* CONSTANTS[58] is beta_1 in component tropomyosin (dimensionless).
* ALGEBRAIC[51] is Ca_50 in component tropomyosin (mM).
* ALGEBRAIC[52] is Ca_TRPN_50 in component tropomyosin (mM).
* CONSTANTS[86] is K_2 in component tropomyosin (per_ms).
* CONSTANTS[90] is K_1 in component tropomyosin (per_ms).
* ALGEBRAIC[53] is alpha_Tm in component tropomyosin (per_ms).
* ALGEBRAIC[47] is beta_Tm in component tropomyosin (per_ms).
* CONSTANTS[59] is beta_0 in component filament_overlap (dimensionless).
* ALGEBRAIC[49] is lambda in component Myofilaments (dimensionless).
* ALGEBRAIC[60] is k_off in component troponin (per_ms).
* ALGEBRAIC[59] is Tension in component Cross_Bridges (N_per_mm2).
* CONSTANTS[60] is T_ref in component length_independent_tension (N_per_mm2).
* ALGEBRAIC[48] is ExtensionRatio in component Myofilaments (dimensionless).
* CONSTANTS[87] is dExtensionRatiodt in component Myofilaments (per_ms).
* ALGEBRAIC[50] is lambda_prev in component Myofilaments (dimensionless).
* ALGEBRAIC[55] is overlap in component filament_overlap (dimensionless).
* ALGEBRAIC[56] is T_Base in component length_independent_tension (N_per_mm2).
* ALGEBRAIC[57] is T_0 in component isometric_tension (N_per_mm2).
* ALGEBRAIC[58] is Q in component Cross_Bridges (dimensionless).
* CONSTANTS[61] is a in component Cross_Bridges (dimensionless).
* STATES[16] is Q_1 in component Cross_Bridges (dimensionless).
* STATES[17] is Q_2 in component Cross_Bridges (dimensionless).
* STATES[18] is Q_3 in component Cross_Bridges (dimensionless).
* CONSTANTS[62] is A_1 in component Cross_Bridges (dimensionless).
* CONSTANTS[63] is A_2 in component Cross_Bridges (dimensionless).
* CONSTANTS[64] is A_3 in component Cross_Bridges (dimensionless).
* CONSTANTS[65] is alpha_1 in component Cross_Bridges (per_ms).
* CONSTANTS[66] is alpha_2 in component Cross_Bridges (per_ms).
* CONSTANTS[67] is alpha_3 in component Cross_Bridges (per_ms).
* ALGEBRAIC[1] is m_infinity in component sodium_current_m_gate (dimensionless).
* ALGEBRAIC[12] is tau_m in component sodium_current_m_gate (second).
* ALGEBRAIC[2] is h_infinity in component sodium_current_h_gate (dimensionless).
* ALGEBRAIC[13] is tau_h in component sodium_current_h_gate (second).
* ALGEBRAIC[3] is j_infinity in component sodium_current_j_gate (dimensionless).
* ALGEBRAIC[14] is tau_j in component sodium_current_j_gate (second).
* ALGEBRAIC[15] is tau_r in component Ca_independent_transient_outward_K_current_r_gate (second).
* ALGEBRAIC[4] is r_infinity in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
* ALGEBRAIC[16] is tau_s_endo in component Ca_independent_transient_outward_K_current_s_gate (second).
* ALGEBRAIC[5] is s_infinity in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
* ALGEBRAIC[17] is tau_s_slow_endo in component Ca_independent_transient_outward_K_current_s_slow_gate (second).
* ALGEBRAIC[6] is s_slow_infinity in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
* ALGEBRAIC[18] is tau_r_ss in component steady_state_outward_K_current_r_ss_gate (second).
* ALGEBRAIC[7] is r_ss_infinity in component steady_state_outward_K_current_r_ss_gate (dimensionless).
* CONSTANTS[88] is tau_s_ss in component steady_state_outward_K_current_s_ss_gate (second).
* ALGEBRAIC[8] is s_ss_infinity in component steady_state_outward_K_current_s_ss_gate (dimensionless).
* ALGEBRAIC[19] is tau_y in component hyperpolarisation_activated_current_y_gate (second).
* ALGEBRAIC[9] is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless).
* ALGEBRAIC[81] is C_oc in component DS_Calcium_Concentrations (mM).
* ALGEBRAIC[79] is C_co in component DS_Calcium_Concentrations (mM).
* CONSTANTS[68] is V_L in component CaRU_Transitions (mV).
* CONSTANTS[69] is del_VL in component CaRU_Transitions (mV).
* CONSTANTS[70] is phi_L in component CaRU_Transitions (dimensionless).
* CONSTANTS[71] is t_L in component CaRU_Transitions (ms).
* CONSTANTS[72] is tau_L in component CaRU_Transitions (ms).
* CONSTANTS[89] is t_R in component CaRU_Transitions (ms).
* CONSTANTS[73] is tau_R in component CaRU_Transitions (ms).
* CONSTANTS[74] is phi_R in component CaRU_Transitions (dimensionless).
* CONSTANTS[75] is theta_R in component CaRU_Transitions (dimensionless).
* CONSTANTS[76] is K_RyR in component CaRU_Transitions (mM).
* CONSTANTS[77] is K_L in component CaRU_Transitions (mM).
* CONSTANTS[78] is a in component CaRU_Transitions (dimensionless).
* CONSTANTS[79] is b in component CaRU_Transitions (dimensionless).
* CONSTANTS[80] is c in component CaRU_Transitions (dimensionless).
* CONSTANTS[81] is d in component CaRU_Transitions (dimensionless).
* ALGEBRAIC[62] is expVL in component CaRU_Transitions (dimensionless).
* ALGEBRAIC[63] is alpha_p in component CaRU_Transitions (per_ms).
* CONSTANTS[91] is alpha_m in component CaRU_Transitions (per_ms).
* ALGEBRAIC[82] is beta_poc in component CaRU_Transitions (per_ms).
* ALGEBRAIC[64] is beta_pcc in component CaRU_Transitions (per_ms).
* CONSTANTS[92] is beta_m in component CaRU_Transitions (per_ms).
* ALGEBRAIC[80] is epsilon_pco in component CaRU_Transitions (per_ms).
* ALGEBRAIC[65] is epsilon_pcc in component CaRU_Transitions (per_ms).
* ALGEBRAIC[66] is epsilon_m in component CaRU_Transitions (per_ms).
* ALGEBRAIC[83] is mu_poc in component CaRU_Transitions (per_ms).
* ALGEBRAIC[67] is mu_pcc in component CaRU_Transitions (per_ms).
* ALGEBRAIC[84] is mu_moc in component CaRU_Transitions (per_ms).
* ALGEBRAIC[68] is mu_mcc in component CaRU_Transitions (per_ms).
* ALGEBRAIC[85] is denom in component CaRU_states (per_ms3).
* ALGEBRAIC[86] is y_oc in component CaRU_states (dimensionless).
* ALGEBRAIC[87] is y_co in component CaRU_states (dimensionless).
* ALGEBRAIC[88] is y_oo in component CaRU_states (dimensionless).
* ALGEBRAIC[90] is y_cc in component CaRU_states (dimensionless).
* ALGEBRAIC[69] is y_ci in component CaRU_states (dimensionless).
* ALGEBRAIC[70] is y_oi in component CaRU_states (dimensionless).
* ALGEBRAIC[71] is y_ic in component CaRU_states (dimensionless).
* ALGEBRAIC[72] is y_io in component CaRU_states (dimensionless).
* ALGEBRAIC[92] is y_ii in component CaRU_states (dimensionless).
* ALGEBRAIC[93] is r_1 in component CaRU_reduced_states (per_ms).
* ALGEBRAIC[95] is r_2 in component CaRU_reduced_states (per_ms).
* ALGEBRAIC[73] is r_3 in component CaRU_reduced_states (per_ms).
* ALGEBRAIC[74] is r_4 in component CaRU_reduced_states (per_ms).
* ALGEBRAIC[97] is r_5 in component CaRU_reduced_states (per_ms).
* ALGEBRAIC[75] is r_6 in component CaRU_reduced_states (per_ms).
* ALGEBRAIC[76] is r_7 in component CaRU_reduced_states (per_ms).
* ALGEBRAIC[77] is r_8 in component CaRU_reduced_states (per_ms).
* STATES[19] is z_1 in component CaRU_reduced_states (dimensionless).
* STATES[20] is z_2 in component CaRU_reduced_states (dimensionless).
* STATES[21] is z_3 in component CaRU_reduced_states (dimensionless).
* ALGEBRAIC[78] is z_4 in component CaRU_reduced_states (dimensionless).
* ALGEBRAIC[10] is C_cc in component DS_Calcium_Concentrations (mM).
* ALGEBRAIC[22] is C_oo in component DS_Calcium_Concentrations (mM).
* ALGEBRAIC[96] is J_Loo in component LCC_and_RyR_fluxes (um3_mM_per_ms).
* ALGEBRAIC[94] is J_Loc in component LCC_and_RyR_fluxes (um3_mM_per_ms).
* ALGEBRAIC[89] is J_Rco in component LCC_and_RyR_fluxes (um3_mM_per_ms).
* ALGEBRAIC[91] is J_Roo in component LCC_and_RyR_fluxes (um3_mM_per_ms).
* ALGEBRAIC[98] is J_R1 in component RyR_current (um3_mM_per_ms).
* ALGEBRAIC[99] is J_R3 in component RyR_current (um3_mM_per_ms).
* ALGEBRAIC[102] is J_L1 in component LCC_current (um3_mM_per_ms).
* ALGEBRAIC[103] is J_L2 in component LCC_current (um3_mM_per_ms).
* RATES[0] is d/dt V in component membrane (mV).
* RATES[1] is d/dt Na_i in component intracellular_ion_concentrations (mM).
* RATES[4] is d/dt K_i in component intracellular_ion_concentrations (mM).
* RATES[5] is d/dt TRPN in component intracellular_ion_concentrations (mM).
* RATES[2] is d/dt Ca_i in component intracellular_ion_concentrations (mM).
* RATES[3] is d/dt Ca_SR in component intracellular_ion_concentrations (mM).
* RATES[15] is d/dt z in component tropomyosin (dimensionless).
* RATES[16] is d/dt Q_1 in component Cross_Bridges (dimensionless).
* RATES[17] is d/dt Q_2 in component Cross_Bridges (dimensionless).
* RATES[18] is d/dt Q_3 in component Cross_Bridges (dimensionless).
* RATES[6] is d/dt m in component sodium_current_m_gate (dimensionless).
* RATES[7] is d/dt h in component sodium_current_h_gate (dimensionless).
* RATES[8] is d/dt j in component sodium_current_j_gate (dimensionless).
* RATES[9] is d/dt r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
* RATES[10] is d/dt s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
* RATES[11] is d/dt s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
* RATES[12] is d/dt r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
* RATES[13] is d/dt s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
* RATES[14] is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless).
* RATES[19] is d/dt z_1 in component CaRU_reduced_states (dimensionless).
* RATES[20] is d/dt z_2 in component CaRU_reduced_states (dimensionless).
* RATES[21] is d/dt z_3 in component CaRU_reduced_states (dimensionless).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -80.50146;
CONSTANTS[0] = 8314.5;
CONSTANTS[1] = 295;
CONSTANTS[2] = 96487;
CONSTANTS[3] = 0.0001;
CONSTANTS[4] = 1e3;
CONSTANTS[5] = 10;
CONSTANTS[6] = -0.6e-3;
CONSTANTS[7] = 25.85e3;
CONSTANTS[8] = 2.098e3;
CONSTANTS[9] = 25.85e-6;
CONSTANTS[10] = 2.098e-6;
CONSTANTS[11] = 0.95e-4;
CONSTANTS[12] = 1.5;
CONSTANTS[13] = 10;
CONSTANTS[14] = 5.4;
CONSTANTS[15] = 140;
STATES[1] = 10.73519;
STATES[2] = 0.00007901351;
STATES[3] = 700.0e-3;
STATES[4] = 139.2751;
STATES[5] = 0.067593139865;
CONSTANTS[16] = 70e-3;
CONSTANTS[17] = 0.8;
STATES[6] = 0.004164108;
STATES[7] = 0.6735613;
STATES[8] = 0.6729362;
CONSTANTS[18] = 0.035;
CONSTANTS[19] = 0.583;
CONSTANTS[20] = 0.417;
STATES[9] = 0.002191519;
STATES[10] = 0.9842542;
STATES[11] = 0.6421196;
CONSTANTS[21] = 0.007;
STATES[12] = 0.002907171;
STATES[13] = 0.3142767;
CONSTANTS[22] = 0.024;
CONSTANTS[23] = 0.00145;
CONSTANTS[24] = 0.2;
STATES[14] = 0.003578708;
CONSTANTS[25] = 0.00008015;
CONSTANTS[26] = 0.0000324;
CONSTANTS[27] = 0.000138;
CONSTANTS[28] = 65;
CONSTANTS[29] = 1.2;
CONSTANTS[30] = 0.065;
CONSTANTS[31] = 0.02;
CONSTANTS[32] = 9.13e-4;
CONSTANTS[33] = 50000;
CONSTANTS[34] = 87.5;
CONSTANTS[35] = 1.38;
CONSTANTS[36] = 0.35;
CONSTANTS[37] = 0.1;
CONSTANTS[38] = 38.5e-3;
CONSTANTS[39] = 0.45e-3;
CONSTANTS[40] = 0.5e-3;
CONSTANTS[41] = 0.0035e-3;
CONSTANTS[42] = 0.5e-3;
CONSTANTS[43] = 2.6875e-8;
CONSTANTS[44] = 1.8951e-5;
CONSTANTS[45] = 2.382e-3;
CONSTANTS[46] = 50e-3;
STATES[15] = 0.014417937837;
CONSTANTS[47] = 100;
CONSTANTS[48] = 0.2;
CONSTANTS[49] = 2;
CONSTANTS[50] = 8e-3;
CONSTANTS[51] = 2e-3;
CONSTANTS[52] = 1.75e-3;
CONSTANTS[53] = 3;
CONSTANTS[54] = 0.15;
CONSTANTS[55] = 3;
CONSTANTS[56] = 1.05e-3;
CONSTANTS[57] = 0.85;
CONSTANTS[58] = -4;
CONSTANTS[59] = 4.9;
CONSTANTS[60] = 56.2;
CONSTANTS[61] = 0.35;
STATES[16] = 0;
STATES[17] = 0;
STATES[18] = 0;
CONSTANTS[62] = -29;
CONSTANTS[63] = 138;
CONSTANTS[64] = 129;
CONSTANTS[65] = 0.03;
CONSTANTS[66] = 0.13;
CONSTANTS[67] = 0.625;
CONSTANTS[68] = -2;
CONSTANTS[69] = 7;
CONSTANTS[70] = 2.35;
CONSTANTS[71] = 1;
CONSTANTS[72] = 650;
CONSTANTS[73] = 2.43;
CONSTANTS[74] = 0.05;
CONSTANTS[75] = 0.012;
CONSTANTS[76] = 41e-3;
CONSTANTS[77] = 0.22e-3;
CONSTANTS[78] = 0.0625;
CONSTANTS[79] = 14;
CONSTANTS[80] = 0.01;
CONSTANTS[81] = 100;
STATES[19] = 0.98859;
STATES[20] = 0.0087302;
STATES[21] = 0.0026566;
CONSTANTS[82] = (exp(CONSTANTS[15]/67.3000) - 1.00000)/7.00000;
CONSTANTS[83] = 1.33000*CONSTANTS[17];
CONSTANTS[84] = 0.464700*CONSTANTS[18];
CONSTANTS[85] = 1.00000 - CONSTANTS[24];
CONSTANTS[86] = (( CONSTANTS[52]*pow(CONSTANTS[57], CONSTANTS[53]))/(pow(CONSTANTS[57], CONSTANTS[53])+pow(CONSTANTS[54], CONSTANTS[53])))*(1.00000 - ( CONSTANTS[53]*pow(CONSTANTS[54], CONSTANTS[53]))/(pow(CONSTANTS[57], CONSTANTS[53])+pow(CONSTANTS[54], CONSTANTS[53])));
CONSTANTS[87] = 0.00000;
CONSTANTS[88] = 2.10000;
CONSTANTS[89] = 1.17000*CONSTANTS[71];
CONSTANTS[90] = ( CONSTANTS[52]*pow(CONSTANTS[57], CONSTANTS[53] - 1.00000)*CONSTANTS[53]*pow(CONSTANTS[54], CONSTANTS[53]))/pow(pow(CONSTANTS[57], CONSTANTS[53])+pow(CONSTANTS[54], CONSTANTS[53]), 2.00000);
CONSTANTS[91] = CONSTANTS[70]/CONSTANTS[71];
CONSTANTS[92] = CONSTANTS[74]/CONSTANTS[89];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[16] = CONSTANTS[62]*CONSTANTS[87] - CONSTANTS[65]*STATES[16];
RATES[17] = CONSTANTS[63]*CONSTANTS[87] - CONSTANTS[66]*STATES[17];
RATES[18] = CONSTANTS[64]*CONSTANTS[87] - CONSTANTS[67]*STATES[18];
ALGEBRAIC[8] = 1.00000/(1.00000+exp((STATES[0]+87.5000)/10.3000));
RATES[13] = ((ALGEBRAIC[8] - STATES[13])/CONSTANTS[88])*0.00100000;
ALGEBRAIC[1] = 1.00000/(1.00000+exp((STATES[0]+45.0000)/- 6.50000));
ALGEBRAIC[12] = 0.00136000/(( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300)))+ 0.0800000*exp(- STATES[0]/11.0000));
RATES[6] = ((ALGEBRAIC[1] - STATES[6])/ALGEBRAIC[12])*0.00100000;
ALGEBRAIC[2] = 1.00000/(1.00000+exp((STATES[0]+76.1000)/6.07000));
ALGEBRAIC[13] = (STATES[0]>=- 40.0000 ? 0.000453700*(1.00000+exp(- (STATES[0]+10.6600)/11.1000)) : 0.00349000/( 0.135000*exp(- (STATES[0]+80.0000)/6.80000)+ 3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0])));
RATES[7] = ((ALGEBRAIC[2] - STATES[7])/ALGEBRAIC[13])*0.00100000;
ALGEBRAIC[3] = 1.00000/(1.00000+exp((STATES[0]+76.1000)/6.07000));
ALGEBRAIC[14] = (STATES[0]>=- 40.0000 ? ( 0.0116300*(1.00000+exp( - 0.100000*(STATES[0]+32.0000))))/exp( - 2.53500e-07*STATES[0]) : 0.00349000/( ((STATES[0]+37.7800)/(1.00000+exp( 0.311000*(STATES[0]+79.2300))))*( - 127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( - 0.0439100*STATES[0]))+( 0.121200*exp( - 0.0105200*STATES[0]))/(1.00000+exp( - 0.137800*(STATES[0]+40.1400)))));
RATES[8] = ((ALGEBRAIC[3] - STATES[8])/ALGEBRAIC[14])*0.00100000;
ALGEBRAIC[15] = 1.00000/( 45.1600*exp( 0.0357700*(STATES[0]+50.0000))+ 98.9000*exp( - 0.100000*(STATES[0]+38.0000)));
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+10.6000)/- 11.4200));
RATES[9] = ((ALGEBRAIC[4] - STATES[9])/ALGEBRAIC[15])*0.00100000;
ALGEBRAIC[16] = 0.550000*exp(- pow((STATES[0]+70.0000)/25.0000, 2.00000))+0.0490000;
ALGEBRAIC[5] = 1.00000/(1.00000+exp((STATES[0]+45.3000)/6.88410));
RATES[10] = ((ALGEBRAIC[5] - STATES[10])/ALGEBRAIC[16])*0.00100000;
ALGEBRAIC[17] = 3.30000*exp(( (- (STATES[0]+70.0000)/30.0000)*(STATES[0]+70.0000))/30.0000)+0.0490000;
ALGEBRAIC[6] = 1.00000/(1.00000+exp((STATES[0]+45.3000)/6.88410));
RATES[11] = ((ALGEBRAIC[6] - STATES[11])/ALGEBRAIC[17])*0.00100000;
ALGEBRAIC[18] = 10.0000/( 45.1600*exp( 0.0357700*(STATES[0]+50.0000))+ 98.9000*exp( - 0.100000*(STATES[0]+38.0000)));
ALGEBRAIC[7] = 1.00000/(1.00000+exp((STATES[0]+11.5000)/- 11.8200));
RATES[12] = ((ALGEBRAIC[7] - STATES[12])/ALGEBRAIC[18])*0.00100000;
ALGEBRAIC[19] = 1.00000/( 0.118850*exp((STATES[0]+80.0000)/28.3700)+ 0.562300*exp((STATES[0]+80.0000)/- 14.1900));
ALGEBRAIC[9] = 1.00000/(1.00000+exp((STATES[0]+138.600)/10.4800));
RATES[14] = ((ALGEBRAIC[9] - STATES[14])/ALGEBRAIC[19])*0.00100000;
ALGEBRAIC[0] = (VOI - floor(VOI/CONSTANTS[4])*CONSTANTS[4]>=0.00000&&VOI - floor(VOI/CONSTANTS[4])*CONSTANTS[4]<=CONSTANTS[5] ? CONSTANTS[6] : 0.00000);
ALGEBRAIC[26] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[14]/STATES[4]);
ALGEBRAIC[27] = CONSTANTS[84]*STATES[9]*( CONSTANTS[19]*STATES[10]+ CONSTANTS[20]*STATES[11])*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[28] = CONSTANTS[21]*STATES[12]*STATES[13]*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[29] = ( (48.0000/(exp((STATES[0]+37.0000)/25.0000)+exp((STATES[0]+37.0000)/- 25.0000))+10.0000)*0.00100000)/(1.00000+exp((STATES[0] - (ALGEBRAIC[26]+76.7700))/- 17.0000))+( CONSTANTS[22]*(STATES[0] - (ALGEBRAIC[26]+1.73000)))/( (1.00000+exp(( 1.61300*CONSTANTS[2]*(STATES[0] - (ALGEBRAIC[26]+1.73000)))/( CONSTANTS[0]*CONSTANTS[1])))*(1.00000+exp((CONSTANTS[14] - 0.998800)/- 0.124000)));
ALGEBRAIC[35] = CONSTANTS[27]*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[21] = ( (( (( CONSTANTS[11]*1.00000)/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ 0.0365000*CONSTANTS[82]*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*CONSTANTS[14])/(CONSTANTS[14]+CONSTANTS[12]))*1.00000)/(1.00000+pow(CONSTANTS[13]/STATES[1], 4.00000));
ALGEBRAIC[31] = CONSTANTS[23]*STATES[14]*CONSTANTS[85]*(STATES[0] - ALGEBRAIC[26]);
RATES[4] = ( - (ALGEBRAIC[0]+ALGEBRAIC[28]*0.00100000+ALGEBRAIC[35]*0.00100000+ALGEBRAIC[27]*0.00100000+ALGEBRAIC[29]*0.00100000+ALGEBRAIC[31]*0.00100000+ ALGEBRAIC[21]*- 2.00000)*1.00000)/( CONSTANTS[9]*CONSTANTS[2]);
ALGEBRAIC[24] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[15]/STATES[1]);
ALGEBRAIC[25] = CONSTANTS[83]*pow(STATES[6], 3.00000)*STATES[7]*STATES[8]*(STATES[0] - ALGEBRAIC[24]);
ALGEBRAIC[33] = CONSTANTS[25]*(STATES[0] - ALGEBRAIC[24]);
ALGEBRAIC[11] = ( CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]);
ALGEBRAIC[37] = ( CONSTANTS[38]*( exp( CONSTANTS[36]*ALGEBRAIC[11])*pow(STATES[1], 3.00000)*CONSTANTS[29] - exp( (CONSTANTS[36] - 1.00000)*ALGEBRAIC[11])*pow(CONSTANTS[15], 3.00000)*STATES[2]))/( (pow(CONSTANTS[15], 3.00000)+pow(CONSTANTS[34], 3.00000))*(CONSTANTS[29]+CONSTANTS[35])*(1.00000+ CONSTANTS[37]*exp( (CONSTANTS[36] - 1.00000)*ALGEBRAIC[11])));
ALGEBRAIC[38] = ALGEBRAIC[37]*CONSTANTS[9]*CONSTANTS[2];
ALGEBRAIC[30] = CONSTANTS[23]*STATES[14]*CONSTANTS[24]*(STATES[0] - ALGEBRAIC[24]);
RATES[1] = ( - (ALGEBRAIC[25]*0.00100000+ALGEBRAIC[33]*0.00100000+ ALGEBRAIC[38]*3.00000+ ALGEBRAIC[21]*3.00000+ALGEBRAIC[30]*0.00100000)*1.00000)/( CONSTANTS[9]*CONSTANTS[2]);
ALGEBRAIC[23] = CONSTANTS[16] - STATES[5];
ALGEBRAIC[48] = (VOI>300000. ? 1.00000 : 1.00000);
ALGEBRAIC[49] = (ALGEBRAIC[48]>0.800000&&ALGEBRAIC[48]<=1.15000 ? ALGEBRAIC[48] : ALGEBRAIC[48]>1.15000 ? 1.15000 : 0.800000);
ALGEBRAIC[51] = CONSTANTS[56]*(1.00000+ CONSTANTS[58]*(ALGEBRAIC[49] - 1.00000));
ALGEBRAIC[52] = ( ALGEBRAIC[51]*CONSTANTS[16])/(ALGEBRAIC[51]+ (CONSTANTS[48]/CONSTANTS[47])*(1.00000 - ( (1.00000+ CONSTANTS[59]*(ALGEBRAIC[49] - 1.00000))*0.500000)/CONSTANTS[49]));
ALGEBRAIC[53] = CONSTANTS[50]*pow(ALGEBRAIC[23]/ALGEBRAIC[52], CONSTANTS[55]);
ALGEBRAIC[47] = CONSTANTS[51]+( CONSTANTS[52]*pow(STATES[15], CONSTANTS[53] - 1.00000))/(pow(STATES[15], CONSTANTS[53])+pow(CONSTANTS[54], CONSTANTS[53]));
RATES[15] = ALGEBRAIC[53]*(1.00000 - STATES[15]) - ALGEBRAIC[47]*STATES[15];
ALGEBRAIC[55] = 1.00000+ CONSTANTS[59]*(ALGEBRAIC[49] - 1.00000);
ALGEBRAIC[54] = (CONSTANTS[50]/pow(ALGEBRAIC[52]/CONSTANTS[16], CONSTANTS[55]) - CONSTANTS[86])/(CONSTANTS[51]+CONSTANTS[90]+CONSTANTS[50]/pow(ALGEBRAIC[52]/CONSTANTS[16], CONSTANTS[55]));
ALGEBRAIC[56] = ( CONSTANTS[60]*STATES[15])/ALGEBRAIC[54];
ALGEBRAIC[57] = ALGEBRAIC[56]*ALGEBRAIC[55];
ALGEBRAIC[58] = STATES[16]+STATES[17]+STATES[18];
ALGEBRAIC[59] = (ALGEBRAIC[58]<0.00000 ? ( ALGEBRAIC[57]*( CONSTANTS[61]*ALGEBRAIC[58]+1.00000))/(1.00000 - ALGEBRAIC[58]) : ( ALGEBRAIC[57]*(1.00000+ (CONSTANTS[61]+2.00000)*ALGEBRAIC[58]))/(1.00000+ALGEBRAIC[58]));
ALGEBRAIC[60] = (1.00000 - ALGEBRAIC[59]/( CONSTANTS[49]*CONSTANTS[60])>0.100000 ? CONSTANTS[48]*(1.00000 - ALGEBRAIC[59]/( CONSTANTS[49]*CONSTANTS[60])) : CONSTANTS[48]*0.100000);
ALGEBRAIC[61] = (CONSTANTS[16] - STATES[5])*ALGEBRAIC[60] - STATES[2]*STATES[5]*CONSTANTS[47];
RATES[5] = ALGEBRAIC[61];
ALGEBRAIC[20] = 2.00000*ALGEBRAIC[11];
ALGEBRAIC[81] = (fabs(ALGEBRAIC[20])>1.00000e-09 ? (STATES[2]+( (CONSTANTS[32]/CONSTANTS[30])*CONSTANTS[29]*ALGEBRAIC[20]*exp(- ALGEBRAIC[20]))/(1.00000 - exp(- ALGEBRAIC[20])))/(1.00000+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20]))) : (STATES[2]+ (CONSTANTS[32]/CONSTANTS[30])*CONSTANTS[29])/(1.00000+CONSTANTS[32]/CONSTANTS[30]));
ALGEBRAIC[83] = (pow(ALGEBRAIC[81], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000))/( CONSTANTS[73]*(pow(ALGEBRAIC[81], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[67] = (pow(STATES[2], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000))/( CONSTANTS[73]*(pow(STATES[2], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[62] = exp((STATES[0] - CONSTANTS[68])/CONSTANTS[69]);
ALGEBRAIC[63] = ALGEBRAIC[62]/( CONSTANTS[71]*(ALGEBRAIC[62]+1.00000));
ALGEBRAIC[64] = pow(STATES[2], 2.00000)/( CONSTANTS[89]*(pow(STATES[2], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[82] = pow(ALGEBRAIC[81], 2.00000)/( CONSTANTS[89]*(pow(ALGEBRAIC[81], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[85] = (ALGEBRAIC[63]+CONSTANTS[91])*( (CONSTANTS[91]+CONSTANTS[92]+ALGEBRAIC[82])*(CONSTANTS[92]+ALGEBRAIC[64])+ ALGEBRAIC[63]*(CONSTANTS[92]+ALGEBRAIC[82]));
ALGEBRAIC[86] = ( ALGEBRAIC[63]*CONSTANTS[92]*(ALGEBRAIC[63]+CONSTANTS[91]+CONSTANTS[92]+ALGEBRAIC[64]))/ALGEBRAIC[85];
ALGEBRAIC[90] = ( CONSTANTS[91]*CONSTANTS[92]*(CONSTANTS[91]+ALGEBRAIC[63]+CONSTANTS[92]+ALGEBRAIC[82]))/ALGEBRAIC[85];
ALGEBRAIC[93] = ALGEBRAIC[86]*ALGEBRAIC[83]+ ALGEBRAIC[90]*ALGEBRAIC[67];
ALGEBRAIC[84] = ( CONSTANTS[75]*CONSTANTS[81]*(pow(ALGEBRAIC[81], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)))/( CONSTANTS[73]*( CONSTANTS[81]*pow(ALGEBRAIC[81], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[68] = ( CONSTANTS[75]*CONSTANTS[81]*(pow(STATES[2], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)))/( CONSTANTS[73]*( CONSTANTS[81]*pow(STATES[2], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[95] = ( ALGEBRAIC[63]*ALGEBRAIC[84]+ CONSTANTS[91]*ALGEBRAIC[68])/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[65] = ( STATES[2]*(ALGEBRAIC[62]+CONSTANTS[78]))/( CONSTANTS[72]*CONSTANTS[77]*(ALGEBRAIC[62]+1.00000));
ALGEBRAIC[76] = ( CONSTANTS[91]*ALGEBRAIC[65])/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[66] = ( CONSTANTS[79]*(ALGEBRAIC[62]+CONSTANTS[78]))/( CONSTANTS[72]*( CONSTANTS[79]*ALGEBRAIC[62]+CONSTANTS[78]));
ALGEBRAIC[77] = ALGEBRAIC[66];
ALGEBRAIC[78] = ((1.00000 - STATES[19]) - STATES[20]) - STATES[21];
RATES[20] = ( ALGEBRAIC[93]*STATES[19] - (ALGEBRAIC[95]+ALGEBRAIC[76])*STATES[20])+ ALGEBRAIC[77]*ALGEBRAIC[78];
ALGEBRAIC[79] = (STATES[2]+ (CONSTANTS[31]/CONSTANTS[30])*STATES[3])/(1.00000+CONSTANTS[31]/CONSTANTS[30]);
ALGEBRAIC[80] = ( ALGEBRAIC[79]*(ALGEBRAIC[62]+CONSTANTS[78]))/( CONSTANTS[72]*CONSTANTS[77]*(ALGEBRAIC[62]+1.00000));
ALGEBRAIC[87] = ( CONSTANTS[91]*( ALGEBRAIC[64]*(CONSTANTS[91]+CONSTANTS[92]+ALGEBRAIC[82])+ ALGEBRAIC[82]*ALGEBRAIC[63]))/ALGEBRAIC[85];
ALGEBRAIC[97] = ALGEBRAIC[87]*ALGEBRAIC[80]+ ALGEBRAIC[90]*ALGEBRAIC[65];
ALGEBRAIC[75] = ALGEBRAIC[66];
RATES[19] = - (ALGEBRAIC[93]+ALGEBRAIC[97])*STATES[19]+ ALGEBRAIC[95]*STATES[20]+ ALGEBRAIC[75]*STATES[21];
ALGEBRAIC[73] = ( CONSTANTS[92]*ALGEBRAIC[67])/(CONSTANTS[92]+ALGEBRAIC[64]);
ALGEBRAIC[74] = ALGEBRAIC[68];
RATES[21] = ( ALGEBRAIC[97]*STATES[19] - (ALGEBRAIC[75]+ALGEBRAIC[73])*STATES[21])+ ALGEBRAIC[74]*ALGEBRAIC[78];
ALGEBRAIC[88] = ( ALGEBRAIC[63]*( ALGEBRAIC[82]*(ALGEBRAIC[63]+CONSTANTS[92]+ALGEBRAIC[64])+ ALGEBRAIC[64]*CONSTANTS[91]))/ALGEBRAIC[85];
ALGEBRAIC[89] = ( CONSTANTS[31]*(STATES[3] - STATES[2]))/(1.00000+CONSTANTS[31]/CONSTANTS[30]);
ALGEBRAIC[91] = (fabs(ALGEBRAIC[20])>1.00000e-05 ? ( CONSTANTS[31]*((STATES[3] - STATES[2])+ (( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20])))*(STATES[3] - CONSTANTS[29]*exp(- ALGEBRAIC[20]))))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20]))) : ( CONSTANTS[31]*((STATES[3] - STATES[2])+ (( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*(STATES[3] - CONSTANTS[29]*exp(- 1.00000e-05))))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05))));
ALGEBRAIC[98] = ALGEBRAIC[88]*ALGEBRAIC[91]+ ALGEBRAIC[89]*ALGEBRAIC[87];
ALGEBRAIC[99] = ( ALGEBRAIC[89]*ALGEBRAIC[64])/(CONSTANTS[92]+ALGEBRAIC[64]);
ALGEBRAIC[100] = ( ( STATES[19]*ALGEBRAIC[98]+ STATES[21]*ALGEBRAIC[99])*CONSTANTS[33])/CONSTANTS[7];
ALGEBRAIC[101] = 1.50000*ALGEBRAIC[100];
ALGEBRAIC[39] = ( CONSTANTS[39]*pow(STATES[2], 2.00000))/(pow(CONSTANTS[40], 2.00000)+pow(STATES[2], 2.00000));
ALGEBRAIC[45] = CONSTANTS[44]*(STATES[3] - STATES[2]);
RATES[3] = (CONSTANTS[9]/CONSTANTS[10])*((- ALGEBRAIC[101]+ALGEBRAIC[39]) - ALGEBRAIC[45]);
ALGEBRAIC[32] = ALGEBRAIC[30]+ALGEBRAIC[31];
ALGEBRAIC[96] = (fabs(ALGEBRAIC[20])>1.00000e-05 ? ( (( CONSTANTS[32]*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20])))*(( CONSTANTS[29]*exp(- ALGEBRAIC[20]) - STATES[2])+ (CONSTANTS[31]/CONSTANTS[30])*( CONSTANTS[29]*exp(- ALGEBRAIC[20]) - STATES[3])))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(ALGEBRAIC[20]))) : ( (( CONSTANTS[32]*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*(( CONSTANTS[29]*exp(- 1.00000e-05) - STATES[2])+ (CONSTANTS[31]/CONSTANTS[30])*( CONSTANTS[29]*exp(- 1.00000e-05) - STATES[3])))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05))));
ALGEBRAIC[94] = (fabs(ALGEBRAIC[20])>1.00000e-05 ? ( (( CONSTANTS[32]*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20])))*( CONSTANTS[29]*exp(- ALGEBRAIC[20]) - STATES[2]))/(1.00000+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20]))) : ( (( CONSTANTS[32]*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*( CONSTANTS[29]*exp(- 1.00000e-05) - STATES[2]))/(1.00000+( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05))));
ALGEBRAIC[102] = ALGEBRAIC[96]*ALGEBRAIC[88]+ ALGEBRAIC[94]*ALGEBRAIC[86];
ALGEBRAIC[103] = ( ALGEBRAIC[94]*ALGEBRAIC[63])/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[104] = ( ( STATES[19]*ALGEBRAIC[102]+ STATES[20]*ALGEBRAIC[103])*CONSTANTS[33])/CONSTANTS[7];
ALGEBRAIC[105] = - 1.50000*ALGEBRAIC[104]*2.00000*CONSTANTS[9]*CONSTANTS[2];
ALGEBRAIC[40] = ( CONSTANTS[41]*STATES[2])/(CONSTANTS[42]+STATES[2]);
ALGEBRAIC[41] = ALGEBRAIC[40]*2.00000*CONSTANTS[9]*CONSTANTS[2];
ALGEBRAIC[42] = (( CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]))*log(CONSTANTS[29]/STATES[2]);
ALGEBRAIC[43] = CONSTANTS[43]*(ALGEBRAIC[42] - STATES[0]);
ALGEBRAIC[44] = - ALGEBRAIC[43]*2.00000*CONSTANTS[9]*CONSTANTS[2];
RATES[0] = - (ALGEBRAIC[25]*0.00100000+ALGEBRAIC[27]*0.00100000+ALGEBRAIC[28]*0.00100000+ALGEBRAIC[32]*0.00100000+ALGEBRAIC[29]*0.00100000+ALGEBRAIC[33]*0.00100000+ALGEBRAIC[35]*0.00100000+ALGEBRAIC[21]+ALGEBRAIC[0]+ALGEBRAIC[44]+ALGEBRAIC[38]+ALGEBRAIC[41]+ALGEBRAIC[105])/CONSTANTS[3];
ALGEBRAIC[46] = pow(1.00000+( CONSTANTS[45]*CONSTANTS[46])/pow(CONSTANTS[45]+STATES[2], 2.00000), - 1.00000);
RATES[2] = ALGEBRAIC[46]*(((ALGEBRAIC[101] - ALGEBRAIC[39])+ALGEBRAIC[45]+ALGEBRAIC[61]) - ( - 2.00000*ALGEBRAIC[38]+ALGEBRAIC[41]+ALGEBRAIC[44]+ALGEBRAIC[105])/( 2.00000*CONSTANTS[9]*CONSTANTS[2]));
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[8] = 1.00000/(1.00000+exp((STATES[0]+87.5000)/10.3000));
ALGEBRAIC[1] = 1.00000/(1.00000+exp((STATES[0]+45.0000)/- 6.50000));
ALGEBRAIC[12] = 0.00136000/(( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300)))+ 0.0800000*exp(- STATES[0]/11.0000));
ALGEBRAIC[2] = 1.00000/(1.00000+exp((STATES[0]+76.1000)/6.07000));
ALGEBRAIC[13] = (STATES[0]>=- 40.0000 ? 0.000453700*(1.00000+exp(- (STATES[0]+10.6600)/11.1000)) : 0.00349000/( 0.135000*exp(- (STATES[0]+80.0000)/6.80000)+ 3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0])));
ALGEBRAIC[3] = 1.00000/(1.00000+exp((STATES[0]+76.1000)/6.07000));
ALGEBRAIC[14] = (STATES[0]>=- 40.0000 ? ( 0.0116300*(1.00000+exp( - 0.100000*(STATES[0]+32.0000))))/exp( - 2.53500e-07*STATES[0]) : 0.00349000/( ((STATES[0]+37.7800)/(1.00000+exp( 0.311000*(STATES[0]+79.2300))))*( - 127140.*exp( 0.244400*STATES[0]) - 3.47400e-05*exp( - 0.0439100*STATES[0]))+( 0.121200*exp( - 0.0105200*STATES[0]))/(1.00000+exp( - 0.137800*(STATES[0]+40.1400)))));
ALGEBRAIC[15] = 1.00000/( 45.1600*exp( 0.0357700*(STATES[0]+50.0000))+ 98.9000*exp( - 0.100000*(STATES[0]+38.0000)));
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+10.6000)/- 11.4200));
ALGEBRAIC[16] = 0.550000*exp(- pow((STATES[0]+70.0000)/25.0000, 2.00000))+0.0490000;
ALGEBRAIC[5] = 1.00000/(1.00000+exp((STATES[0]+45.3000)/6.88410));
ALGEBRAIC[17] = 3.30000*exp(( (- (STATES[0]+70.0000)/30.0000)*(STATES[0]+70.0000))/30.0000)+0.0490000;
ALGEBRAIC[6] = 1.00000/(1.00000+exp((STATES[0]+45.3000)/6.88410));
ALGEBRAIC[18] = 10.0000/( 45.1600*exp( 0.0357700*(STATES[0]+50.0000))+ 98.9000*exp( - 0.100000*(STATES[0]+38.0000)));
ALGEBRAIC[7] = 1.00000/(1.00000+exp((STATES[0]+11.5000)/- 11.8200));
ALGEBRAIC[19] = 1.00000/( 0.118850*exp((STATES[0]+80.0000)/28.3700)+ 0.562300*exp((STATES[0]+80.0000)/- 14.1900));
ALGEBRAIC[9] = 1.00000/(1.00000+exp((STATES[0]+138.600)/10.4800));
ALGEBRAIC[0] = (VOI - floor(VOI/CONSTANTS[4])*CONSTANTS[4]>=0.00000&&VOI - floor(VOI/CONSTANTS[4])*CONSTANTS[4]<=CONSTANTS[5] ? CONSTANTS[6] : 0.00000);
ALGEBRAIC[26] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[14]/STATES[4]);
ALGEBRAIC[27] = CONSTANTS[84]*STATES[9]*( CONSTANTS[19]*STATES[10]+ CONSTANTS[20]*STATES[11])*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[28] = CONSTANTS[21]*STATES[12]*STATES[13]*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[29] = ( (48.0000/(exp((STATES[0]+37.0000)/25.0000)+exp((STATES[0]+37.0000)/- 25.0000))+10.0000)*0.00100000)/(1.00000+exp((STATES[0] - (ALGEBRAIC[26]+76.7700))/- 17.0000))+( CONSTANTS[22]*(STATES[0] - (ALGEBRAIC[26]+1.73000)))/( (1.00000+exp(( 1.61300*CONSTANTS[2]*(STATES[0] - (ALGEBRAIC[26]+1.73000)))/( CONSTANTS[0]*CONSTANTS[1])))*(1.00000+exp((CONSTANTS[14] - 0.998800)/- 0.124000)));
ALGEBRAIC[35] = CONSTANTS[27]*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[21] = ( (( (( CONSTANTS[11]*1.00000)/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ 0.0365000*CONSTANTS[82]*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*CONSTANTS[14])/(CONSTANTS[14]+CONSTANTS[12]))*1.00000)/(1.00000+pow(CONSTANTS[13]/STATES[1], 4.00000));
ALGEBRAIC[31] = CONSTANTS[23]*STATES[14]*CONSTANTS[85]*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[24] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[15]/STATES[1]);
ALGEBRAIC[25] = CONSTANTS[83]*pow(STATES[6], 3.00000)*STATES[7]*STATES[8]*(STATES[0] - ALGEBRAIC[24]);
ALGEBRAIC[33] = CONSTANTS[25]*(STATES[0] - ALGEBRAIC[24]);
ALGEBRAIC[11] = ( CONSTANTS[2]*STATES[0])/( CONSTANTS[0]*CONSTANTS[1]);
ALGEBRAIC[37] = ( CONSTANTS[38]*( exp( CONSTANTS[36]*ALGEBRAIC[11])*pow(STATES[1], 3.00000)*CONSTANTS[29] - exp( (CONSTANTS[36] - 1.00000)*ALGEBRAIC[11])*pow(CONSTANTS[15], 3.00000)*STATES[2]))/( (pow(CONSTANTS[15], 3.00000)+pow(CONSTANTS[34], 3.00000))*(CONSTANTS[29]+CONSTANTS[35])*(1.00000+ CONSTANTS[37]*exp( (CONSTANTS[36] - 1.00000)*ALGEBRAIC[11])));
ALGEBRAIC[38] = ALGEBRAIC[37]*CONSTANTS[9]*CONSTANTS[2];
ALGEBRAIC[30] = CONSTANTS[23]*STATES[14]*CONSTANTS[24]*(STATES[0] - ALGEBRAIC[24]);
ALGEBRAIC[23] = CONSTANTS[16] - STATES[5];
ALGEBRAIC[48] = (VOI>300000. ? 1.00000 : 1.00000);
ALGEBRAIC[49] = (ALGEBRAIC[48]>0.800000&&ALGEBRAIC[48]<=1.15000 ? ALGEBRAIC[48] : ALGEBRAIC[48]>1.15000 ? 1.15000 : 0.800000);
ALGEBRAIC[51] = CONSTANTS[56]*(1.00000+ CONSTANTS[58]*(ALGEBRAIC[49] - 1.00000));
ALGEBRAIC[52] = ( ALGEBRAIC[51]*CONSTANTS[16])/(ALGEBRAIC[51]+ (CONSTANTS[48]/CONSTANTS[47])*(1.00000 - ( (1.00000+ CONSTANTS[59]*(ALGEBRAIC[49] - 1.00000))*0.500000)/CONSTANTS[49]));
ALGEBRAIC[53] = CONSTANTS[50]*pow(ALGEBRAIC[23]/ALGEBRAIC[52], CONSTANTS[55]);
ALGEBRAIC[47] = CONSTANTS[51]+( CONSTANTS[52]*pow(STATES[15], CONSTANTS[53] - 1.00000))/(pow(STATES[15], CONSTANTS[53])+pow(CONSTANTS[54], CONSTANTS[53]));
ALGEBRAIC[55] = 1.00000+ CONSTANTS[59]*(ALGEBRAIC[49] - 1.00000);
ALGEBRAIC[54] = (CONSTANTS[50]/pow(ALGEBRAIC[52]/CONSTANTS[16], CONSTANTS[55]) - CONSTANTS[86])/(CONSTANTS[51]+CONSTANTS[90]+CONSTANTS[50]/pow(ALGEBRAIC[52]/CONSTANTS[16], CONSTANTS[55]));
ALGEBRAIC[56] = ( CONSTANTS[60]*STATES[15])/ALGEBRAIC[54];
ALGEBRAIC[57] = ALGEBRAIC[56]*ALGEBRAIC[55];
ALGEBRAIC[58] = STATES[16]+STATES[17]+STATES[18];
ALGEBRAIC[59] = (ALGEBRAIC[58]<0.00000 ? ( ALGEBRAIC[57]*( CONSTANTS[61]*ALGEBRAIC[58]+1.00000))/(1.00000 - ALGEBRAIC[58]) : ( ALGEBRAIC[57]*(1.00000+ (CONSTANTS[61]+2.00000)*ALGEBRAIC[58]))/(1.00000+ALGEBRAIC[58]));
ALGEBRAIC[60] = (1.00000 - ALGEBRAIC[59]/( CONSTANTS[49]*CONSTANTS[60])>0.100000 ? CONSTANTS[48]*(1.00000 - ALGEBRAIC[59]/( CONSTANTS[49]*CONSTANTS[60])) : CONSTANTS[48]*0.100000);
ALGEBRAIC[61] = (CONSTANTS[16] - STATES[5])*ALGEBRAIC[60] - STATES[2]*STATES[5]*CONSTANTS[47];
ALGEBRAIC[20] = 2.00000*ALGEBRAIC[11];
ALGEBRAIC[81] = (fabs(ALGEBRAIC[20])>1.00000e-09 ? (STATES[2]+( (CONSTANTS[32]/CONSTANTS[30])*CONSTANTS[29]*ALGEBRAIC[20]*exp(- ALGEBRAIC[20]))/(1.00000 - exp(- ALGEBRAIC[20])))/(1.00000+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20]))) : (STATES[2]+ (CONSTANTS[32]/CONSTANTS[30])*CONSTANTS[29])/(1.00000+CONSTANTS[32]/CONSTANTS[30]));
ALGEBRAIC[83] = (pow(ALGEBRAIC[81], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000))/( CONSTANTS[73]*(pow(ALGEBRAIC[81], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[67] = (pow(STATES[2], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000))/( CONSTANTS[73]*(pow(STATES[2], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[62] = exp((STATES[0] - CONSTANTS[68])/CONSTANTS[69]);
ALGEBRAIC[63] = ALGEBRAIC[62]/( CONSTANTS[71]*(ALGEBRAIC[62]+1.00000));
ALGEBRAIC[64] = pow(STATES[2], 2.00000)/( CONSTANTS[89]*(pow(STATES[2], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[82] = pow(ALGEBRAIC[81], 2.00000)/( CONSTANTS[89]*(pow(ALGEBRAIC[81], 2.00000)+pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[85] = (ALGEBRAIC[63]+CONSTANTS[91])*( (CONSTANTS[91]+CONSTANTS[92]+ALGEBRAIC[82])*(CONSTANTS[92]+ALGEBRAIC[64])+ ALGEBRAIC[63]*(CONSTANTS[92]+ALGEBRAIC[82]));
ALGEBRAIC[86] = ( ALGEBRAIC[63]*CONSTANTS[92]*(ALGEBRAIC[63]+CONSTANTS[91]+CONSTANTS[92]+ALGEBRAIC[64]))/ALGEBRAIC[85];
ALGEBRAIC[90] = ( CONSTANTS[91]*CONSTANTS[92]*(CONSTANTS[91]+ALGEBRAIC[63]+CONSTANTS[92]+ALGEBRAIC[82]))/ALGEBRAIC[85];
ALGEBRAIC[93] = ALGEBRAIC[86]*ALGEBRAIC[83]+ ALGEBRAIC[90]*ALGEBRAIC[67];
ALGEBRAIC[84] = ( CONSTANTS[75]*CONSTANTS[81]*(pow(ALGEBRAIC[81], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)))/( CONSTANTS[73]*( CONSTANTS[81]*pow(ALGEBRAIC[81], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[68] = ( CONSTANTS[75]*CONSTANTS[81]*(pow(STATES[2], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)))/( CONSTANTS[73]*( CONSTANTS[81]*pow(STATES[2], 2.00000)+ CONSTANTS[80]*pow(CONSTANTS[76], 2.00000)));
ALGEBRAIC[95] = ( ALGEBRAIC[63]*ALGEBRAIC[84]+ CONSTANTS[91]*ALGEBRAIC[68])/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[65] = ( STATES[2]*(ALGEBRAIC[62]+CONSTANTS[78]))/( CONSTANTS[72]*CONSTANTS[77]*(ALGEBRAIC[62]+1.00000));
ALGEBRAIC[76] = ( CONSTANTS[91]*ALGEBRAIC[65])/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[66] = ( CONSTANTS[79]*(ALGEBRAIC[62]+CONSTANTS[78]))/( CONSTANTS[72]*( CONSTANTS[79]*ALGEBRAIC[62]+CONSTANTS[78]));
ALGEBRAIC[77] = ALGEBRAIC[66];
ALGEBRAIC[78] = ((1.00000 - STATES[19]) - STATES[20]) - STATES[21];
ALGEBRAIC[79] = (STATES[2]+ (CONSTANTS[31]/CONSTANTS[30])*STATES[3])/(1.00000+CONSTANTS[31]/CONSTANTS[30]);
ALGEBRAIC[80] = ( ALGEBRAIC[79]*(ALGEBRAIC[62]+CONSTANTS[78]))/( CONSTANTS[72]*CONSTANTS[77]*(ALGEBRAIC[62]+1.00000));
ALGEBRAIC[87] = ( CONSTANTS[91]*( ALGEBRAIC[64]*(CONSTANTS[91]+CONSTANTS[92]+ALGEBRAIC[82])+ ALGEBRAIC[82]*ALGEBRAIC[63]))/ALGEBRAIC[85];
ALGEBRAIC[97] = ALGEBRAIC[87]*ALGEBRAIC[80]+ ALGEBRAIC[90]*ALGEBRAIC[65];
ALGEBRAIC[75] = ALGEBRAIC[66];
ALGEBRAIC[73] = ( CONSTANTS[92]*ALGEBRAIC[67])/(CONSTANTS[92]+ALGEBRAIC[64]);
ALGEBRAIC[74] = ALGEBRAIC[68];
ALGEBRAIC[88] = ( ALGEBRAIC[63]*( ALGEBRAIC[82]*(ALGEBRAIC[63]+CONSTANTS[92]+ALGEBRAIC[64])+ ALGEBRAIC[64]*CONSTANTS[91]))/ALGEBRAIC[85];
ALGEBRAIC[89] = ( CONSTANTS[31]*(STATES[3] - STATES[2]))/(1.00000+CONSTANTS[31]/CONSTANTS[30]);
ALGEBRAIC[91] = (fabs(ALGEBRAIC[20])>1.00000e-05 ? ( CONSTANTS[31]*((STATES[3] - STATES[2])+ (( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20])))*(STATES[3] - CONSTANTS[29]*exp(- ALGEBRAIC[20]))))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20]))) : ( CONSTANTS[31]*((STATES[3] - STATES[2])+ (( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*(STATES[3] - CONSTANTS[29]*exp(- 1.00000e-05))))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05))));
ALGEBRAIC[98] = ALGEBRAIC[88]*ALGEBRAIC[91]+ ALGEBRAIC[89]*ALGEBRAIC[87];
ALGEBRAIC[99] = ( ALGEBRAIC[89]*ALGEBRAIC[64])/(CONSTANTS[92]+ALGEBRAIC[64]);
ALGEBRAIC[100] = ( ( STATES[19]*ALGEBRAIC[98]+ STATES[21]*ALGEBRAIC[99])*CONSTANTS[33])/CONSTANTS[7];
ALGEBRAIC[101] = 1.50000*ALGEBRAIC[100];
ALGEBRAIC[39] = ( CONSTANTS[39]*pow(STATES[2], 2.00000))/(pow(CONSTANTS[40], 2.00000)+pow(STATES[2], 2.00000));
ALGEBRAIC[45] = CONSTANTS[44]*(STATES[3] - STATES[2]);
ALGEBRAIC[32] = ALGEBRAIC[30]+ALGEBRAIC[31];
ALGEBRAIC[96] = (fabs(ALGEBRAIC[20])>1.00000e-05 ? ( (( CONSTANTS[32]*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20])))*(( CONSTANTS[29]*exp(- ALGEBRAIC[20]) - STATES[2])+ (CONSTANTS[31]/CONSTANTS[30])*( CONSTANTS[29]*exp(- ALGEBRAIC[20]) - STATES[3])))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(ALGEBRAIC[20]))) : ( (( CONSTANTS[32]*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*(( CONSTANTS[29]*exp(- 1.00000e-05) - STATES[2])+ (CONSTANTS[31]/CONSTANTS[30])*( CONSTANTS[29]*exp(- 1.00000e-05) - STATES[3])))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05))));
ALGEBRAIC[94] = (fabs(ALGEBRAIC[20])>1.00000e-05 ? ( (( CONSTANTS[32]*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20])))*( CONSTANTS[29]*exp(- ALGEBRAIC[20]) - STATES[2]))/(1.00000+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20]))) : ( (( CONSTANTS[32]*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*( CONSTANTS[29]*exp(- 1.00000e-05) - STATES[2]))/(1.00000+( (CONSTANTS[32]/CONSTANTS[30])*1.00000e-05)/(1.00000 - exp(- 1.00000e-05))));
ALGEBRAIC[102] = ALGEBRAIC[96]*ALGEBRAIC[88]+ ALGEBRAIC[94]*ALGEBRAIC[86];
ALGEBRAIC[103] = ( ALGEBRAIC[94]*ALGEBRAIC[63])/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[104] = ( ( STATES[19]*ALGEBRAIC[102]+ STATES[20]*ALGEBRAIC[103])*CONSTANTS[33])/CONSTANTS[7];
ALGEBRAIC[105] = - 1.50000*ALGEBRAIC[104]*2.00000*CONSTANTS[9]*CONSTANTS[2];
ALGEBRAIC[40] = ( CONSTANTS[41]*STATES[2])/(CONSTANTS[42]+STATES[2]);
ALGEBRAIC[41] = ALGEBRAIC[40]*2.00000*CONSTANTS[9]*CONSTANTS[2];
ALGEBRAIC[42] = (( CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]))*log(CONSTANTS[29]/STATES[2]);
ALGEBRAIC[43] = CONSTANTS[43]*(ALGEBRAIC[42] - STATES[0]);
ALGEBRAIC[44] = - ALGEBRAIC[43]*2.00000*CONSTANTS[9]*CONSTANTS[2];
ALGEBRAIC[46] = pow(1.00000+( CONSTANTS[45]*CONSTANTS[46])/pow(CONSTANTS[45]+STATES[2], 2.00000), - 1.00000);
ALGEBRAIC[10] = STATES[2];
ALGEBRAIC[22] = (fabs(ALGEBRAIC[20])>1.00000e-09 ? (STATES[2]+ (CONSTANTS[31]/CONSTANTS[30])*STATES[3]+( (CONSTANTS[32]/CONSTANTS[30])*CONSTANTS[29]*ALGEBRAIC[20]*exp(- ALGEBRAIC[20]))/(1.00000 - exp(- ALGEBRAIC[20])))/(1.00000+CONSTANTS[31]/CONSTANTS[30]+( (CONSTANTS[32]/CONSTANTS[30])*ALGEBRAIC[20])/(1.00000 - exp(- ALGEBRAIC[20]))) : (STATES[2]+ (CONSTANTS[31]/CONSTANTS[30])*STATES[3]+ (CONSTANTS[32]/CONSTANTS[30])*CONSTANTS[29])/(1.00000+CONSTANTS[31]/CONSTANTS[30]+CONSTANTS[32]/CONSTANTS[30]));
ALGEBRAIC[34] = CONSTANTS[26]*(STATES[0] - CONSTANTS[28]);
ALGEBRAIC[36] = ALGEBRAIC[33]+ALGEBRAIC[34]+ALGEBRAIC[35];
ALGEBRAIC[50] = ALGEBRAIC[48];
ALGEBRAIC[69] = CONSTANTS[91]/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[70] = ALGEBRAIC[63]/(ALGEBRAIC[63]+CONSTANTS[91]);
ALGEBRAIC[71] = CONSTANTS[92]/(ALGEBRAIC[64]+CONSTANTS[92]);
ALGEBRAIC[72] = ALGEBRAIC[64]/(ALGEBRAIC[64]+CONSTANTS[92]);
ALGEBRAIC[92] = (((((((1.00000 - ALGEBRAIC[86]) - ALGEBRAIC[87]) - ALGEBRAIC[88]) - ALGEBRAIC[90]) - ALGEBRAIC[69]) - ALGEBRAIC[71]) - ALGEBRAIC[70]) - ALGEBRAIC[72];
}