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 62 entries in the algebraic variable array.
There are a total of 31 entries in each of the rate and state variable arrays.
There are a total of 70 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 (gas_constant_units).
* CONSTANTS[1] is T in component membrane (kelvin).
* CONSTANTS[2] is F in component membrane (faradays_constant_units).
* ALGEBRAIC[15] is i_Na in component fast_sodium_current (uA_per_mm2).
* ALGEBRAIC[19] is i_Ca_L_Ca in component L_type_Ca_channel (uA_per_mm2).
* ALGEBRAIC[21] is i_Ca_L_K in component L_type_Ca_channel (uA_per_mm2).
* ALGEBRAIC[25] is i_K in component time_dependent_potassium_current (uA_per_mm2).
* ALGEBRAIC[31] is i_K1 in component time_independent_potassium_current (uA_per_mm2).
* ALGEBRAIC[35] is i_NaCa in component Na_Ca_exchanger (uA_per_mm2).
* ALGEBRAIC[34] is i_Kp in component plateau_potassium_current (uA_per_mm2).
* ALGEBRAIC[36] is i_p_Ca in component sarcolemmal_calcium_pump (uA_per_mm2).
* ALGEBRAIC[38] is i_Na_b in component sodium_background_current (uA_per_mm2).
* ALGEBRAIC[40] is i_Ca_b in component calcium_background_current (uA_per_mm2).
* ALGEBRAIC[43] is i_NaK in component sodium_potassium_pump (uA_per_mm2).
* ALGEBRAIC[60] is i_ns_Ca in component non_specific_calcium_activated_current (uA_per_mm2).
* CONSTANTS[3] is Cm in component membrane (uF_per_mm2).
* ALGEBRAIC[6] is I_stim in component membrane (uA_per_mm2).
* CONSTANTS[4] is stim_start in component membrane (ms).
* CONSTANTS[5] is stim_end in component membrane (ms).
* CONSTANTS[6] is stim_period in component membrane (ms).
* CONSTANTS[7] is stim_duration in component membrane (ms).
* CONSTANTS[8] is stim_amplitude in component membrane (uA_per_mm2).
* ALGEBRAIC[13] is E_Na in component fast_sodium_current (mV).
* CONSTANTS[9] is g_Na in component fast_sodium_current (mS_per_mm2).
* STATES[1] is Nai in component ionic_concentrations (mM).
* CONSTANTS[10] 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[7] 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[8] 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[9] is beta_j in component fast_sodium_current_j_gate (per_ms).
* CONSTANTS[11] is P_Ca in component L_type_Ca_channel (mm_per_ms).
* CONSTANTS[12] is P_K in component L_type_Ca_channel (mm_per_ms).
* ALGEBRAIC[20] is p_k in component L_type_Ca_channel (mm_per_ms).
* CONSTANTS[13] is i_Ca_L_Ca_half in component L_type_Ca_channel (uA_per_mm2).
* ALGEBRAIC[17] 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[10] is beta in component L_type_Ca_channel (per_ms).
* ALGEBRAIC[18] is gamma in component L_type_Ca_channel (per_ms).
* ALGEBRAIC[14] is alpha_a in component L_type_Ca_channel (per_ms).
* ALGEBRAIC[16] is beta_b in component L_type_Ca_channel (per_ms).
* CONSTANTS[14] is a in component L_type_Ca_channel (dimensionless).
* CONSTANTS[15] is b in component L_type_Ca_channel (dimensionless).
* CONSTANTS[16] is g in component L_type_Ca_channel (per_ms).
* CONSTANTS[17] is f in component L_type_Ca_channel (per_ms).
* 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 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[21] 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[11] is tau_y in component L_type_Ca_channel_y_gate (ms).
* ALGEBRAIC[22] is g_K in component time_dependent_potassium_current (mS_per_mm2).
* CONSTANTS[22] is g_K_max in component time_dependent_potassium_current (mS_per_mm2).
* ALGEBRAIC[23] is E_K in component time_dependent_potassium_current (mV).
* CONSTANTS[23] is P_NaK in component time_dependent_potassium_current (dimensionless).
* STATES[21] is X in component time_dependent_potassium_current_X_gate (dimensionless).
* ALGEBRAIC[24] 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[12] is beta_X in component time_dependent_potassium_current_X_gate (per_ms).
* ALGEBRAIC[27] is E_K1 in component time_independent_potassium_current (mV).
* ALGEBRAIC[26] is g_K1 in component time_independent_potassium_current (mS_per_mm2).
* CONSTANTS[24] is g_K1_max in component time_independent_potassium_current (mS_per_mm2).
* ALGEBRAIC[30] is K1_infinity in component time_independent_potassium_current_K1_gate (dimensionless).
* ALGEBRAIC[28] is alpha_K1 in component time_independent_potassium_current_K1_gate (per_ms).
* ALGEBRAIC[29] is beta_K1 in component time_independent_potassium_current_K1_gate (per_ms).
* ALGEBRAIC[32] is E_Kp in component plateau_potassium_current (mV).
* CONSTANTS[25] is g_Kp in component plateau_potassium_current (mS_per_mm2).
* ALGEBRAIC[33] is Kp in component plateau_potassium_current (dimensionless).
* CONSTANTS[26] is k_NaCa in component Na_Ca_exchanger (uA_per_mm2).
* CONSTANTS[27] is K_mNa in component Na_Ca_exchanger (mM).
* CONSTANTS[28] is K_mCa in component Na_Ca_exchanger (mM).
* CONSTANTS[29] is k_sat in component Na_Ca_exchanger (dimensionless).
* CONSTANTS[30] is eta in component Na_Ca_exchanger (dimensionless).
* STATES[22] is Cai in component calcium_subsystem (mM).
* CONSTANTS[31] is K_mpCa in component sarcolemmal_calcium_pump (mM).
* CONSTANTS[32] is I_pCa in component sarcolemmal_calcium_pump (uA_per_mm2).
* CONSTANTS[33] is g_Nab in component sodium_background_current (mS_per_mm2).
* ALGEBRAIC[37] is E_NaN in component sodium_background_current (mV).
* CONSTANTS[34] is g_Cab in component calcium_background_current (mS_per_mm2).
* ALGEBRAIC[39] is E_CaN in component calcium_background_current (mV).
* CONSTANTS[35] is I_NaK in component sodium_potassium_pump (uA_per_mm2).
* ALGEBRAIC[41] is f_NaK in component sodium_potassium_pump (dimensionless).
* CONSTANTS[36] is K_mNai in component sodium_potassium_pump (mM).
* CONSTANTS[37] is K_mKo in component sodium_potassium_pump (mM).
* CONSTANTS[66] is sigma in component sodium_potassium_pump (dimensionless).
* ALGEBRAIC[52] is i_ns_Na in component non_specific_calcium_activated_current (uA_per_mm2).
* ALGEBRAIC[58] is i_ns_K in component non_specific_calcium_activated_current (uA_per_mm2).
* ALGEBRAIC[49] is I_ns_Na in component non_specific_calcium_activated_current (uA_per_mm2).
* ALGEBRAIC[56] is I_ns_K in component non_specific_calcium_activated_current (uA_per_mm2).
* CONSTANTS[38] is K_m_ns_Ca in component non_specific_calcium_activated_current (mM).
* CONSTANTS[39] is P_ns_Ca in component non_specific_calcium_activated_current (mm_per_ms).
* ALGEBRAIC[45] is EnsCa in component non_specific_calcium_activated_current (mV).
* ALGEBRAIC[47] is VnsCa in component non_specific_calcium_activated_current (mV).
* CONSTANTS[40] is Am in component calcium_subsystem (per_mm).
* CONSTANTS[41] is V_myo in component calcium_subsystem (dimensionless).
* ALGEBRAIC[42] 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[42] is v1 in component calcium_subsystem (per_ms).
* CONSTANTS[43] is v2 in component calcium_subsystem (per_ms).
* CONSTANTS[44] is v3 in component calcium_subsystem (mM_per_ms).
* CONSTANTS[45] is nCa in component calcium_subsystem (dimensionless).
* CONSTANTS[46] is mCa in component calcium_subsystem (dimensionless).
* CONSTANTS[47] is k_a_plus in component calcium_subsystem (per_mM4_per_ms).
* CONSTANTS[48] is k_a_minus in component calcium_subsystem (per_ms).
* CONSTANTS[49] is k_b_plus in component calcium_subsystem (per_mM3_per_ms).
* CONSTANTS[50] is k_b_minus in component calcium_subsystem (per_ms).
* CONSTANTS[51] is k_c_plus in component calcium_subsystem (per_ms).
* CONSTANTS[52] is k_c_minus in component calcium_subsystem (per_ms).
* CONSTANTS[53] is k_htrpn_plus in component calcium_subsystem (per_mM_per_ms).
* CONSTANTS[54] is k_htrpn_minus in component calcium_subsystem (per_ms).
* CONSTANTS[55] is k_ltrpn_plus in component calcium_subsystem (per_mM_per_ms).
* CONSTANTS[56] is k_ltrpn_minus in component calcium_subsystem (per_ms).
* CONSTANTS[57] 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[69] is V_JSR in component calcium_subsystem (dimensionless).
* CONSTANTS[68] is V_NSR in component calcium_subsystem (dimensionless).
* CONSTANTS[67] is V_SS in component calcium_subsystem (dimensionless).
* CONSTANTS[58] is K_mup in component calcium_subsystem (mM).
* CONSTANTS[59] is K_mCMDN in component calcium_subsystem (mM).
* CONSTANTS[60] is K_mCSQN in component calcium_subsystem (mM).
* CONSTANTS[61] is tau_xfer in component calcium_subsystem (ms).
* CONSTANTS[62] is HTRPN_tot in component calcium_subsystem (mM).
* CONSTANTS[63] 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[64] is CSQN_tot in component calcium_subsystem (mM).
* CONSTANTS[65] is CMDN_tot in component calcium_subsystem (mM).
* ALGEBRAIC[61] is Bi in component calcium_subsystem (dimensionless).
* ALGEBRAIC[53] is B_SS in component calcium_subsystem (dimensionless).
* ALGEBRAIC[54] is B_JSR in component calcium_subsystem (dimensionless).
* ALGEBRAIC[44] is J_rel in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[46] is J_leak in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[48] is J_up in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[50] is J_tr in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[51] is J_xfer in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[59] is J_trpn in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[55] is J_htrpn in component calcium_subsystem (mM_per_ms).
* ALGEBRAIC[57] is J_ltrpn in component calcium_subsystem (mM_per_ms).
* 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).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -84.1638;
CONSTANTS[0] = 8.3145e3;
CONSTANTS[1] = 310;
CONSTANTS[2] = 9.6845e4;
CONSTANTS[3] = 0.01;
CONSTANTS[4] = 100;
CONSTANTS[5] = 10100;
CONSTANTS[6] = 500;
CONSTANTS[7] = 1;
CONSTANTS[8] = 0.516289;
CONSTANTS[9] = 0.128;
STATES[1] = 10.2042;
CONSTANTS[10] = 140;
STATES[2] = 0.0328302;
STATES[3] = 0.988354;
STATES[4] = 0.99254;
CONSTANTS[11] = 33.75e-6;
CONSTANTS[12] = 1e-9;
CONSTANTS[13] = -4.58e-3;
STATES[5] = 9.84546e-21;
STATES[6] = 0;
CONSTANTS[14] = 2;
CONSTANTS[15] = 2;
CONSTANTS[16] = 2;
CONSTANTS[17] = 0.3;
CONSTANTS[18] = 0;
CONSTANTS[19] = 0;
CONSTANTS[20] = 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[21] = 1.8;
STATES[18] = 5.4;
STATES[19] = 143.727;
STATES[20] = 0.998983;
CONSTANTS[22] = 0.001128;
CONSTANTS[23] = 0.01833;
STATES[21] = 0.000928836;
CONSTANTS[24] = 7.5e-3;
CONSTANTS[25] = 8.28e-5;
CONSTANTS[26] = 50;
CONSTANTS[27] = 87.5;
CONSTANTS[28] = 1.38;
CONSTANTS[29] = 0.1;
CONSTANTS[30] = 0.35;
STATES[22] = 9.94893e-11;
CONSTANTS[31] = 0.5e-3;
CONSTANTS[32] = 1.15e-2;
CONSTANTS[33] = 1.41e-5;
CONSTANTS[34] = 6.032e-5;
CONSTANTS[35] = 0.013;
CONSTANTS[36] = 10;
CONSTANTS[37] = 1.5;
CONSTANTS[38] = 1.2e-3;
CONSTANTS[39] = 1.75e-9;
CONSTANTS[40] = 546.69;
CONSTANTS[41] = 0.92;
STATES[23] = 1.19168e-3;
STATES[24] = 6.30613e-9;
STATES[25] = 0.762527;
STATES[26] = 0.236283;
CONSTANTS[42] = 1.8;
CONSTANTS[43] = 0.58e-4;
CONSTANTS[44] = 1.8e-3;
CONSTANTS[45] = 4;
CONSTANTS[46] = 3;
CONSTANTS[47] = 1.215e10;
CONSTANTS[48] = 0.1425;
CONSTANTS[49] = 4.05e7;
CONSTANTS[50] = 1.93;
CONSTANTS[51] = 0.018;
CONSTANTS[52] = 0.0008;
CONSTANTS[53] = 20;
CONSTANTS[54] = 0.066e-3;
CONSTANTS[55] = 40;
CONSTANTS[56] = 0.04;
CONSTANTS[57] = 34.48;
STATES[27] = 1.17504;
STATES[28] = 1.243891;
CONSTANTS[58] = 0.5e-3;
CONSTANTS[59] = 2.38e-3;
CONSTANTS[60] = 0.8;
CONSTANTS[61] = 3.125;
CONSTANTS[62] = 0.14;
CONSTANTS[63] = 0.07;
STATES[29] = 0.13598;
STATES[30] = 0.00635;
CONSTANTS[64] = 15;
CONSTANTS[65] = 0.05;
CONSTANTS[66] = (1.00000/7.00000)*(exp(CONSTANTS[10]/67.3000) - 1.00000);
CONSTANTS[67] = 5.82800e-05*CONSTANTS[41];
CONSTANTS[68] = 0.0810000*CONSTANTS[41];
CONSTANTS[69] = 0.00464000*CONSTANTS[41];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[5] = CONSTANTS[17]*STATES[11] - CONSTANTS[16]*STATES[5];
RATES[6] = CONSTANTS[19]*STATES[16] - CONSTANTS[18]*STATES[6];
RATES[25] = - CONSTANTS[47]*pow(STATES[17], CONSTANTS[45])*STATES[25]+ CONSTANTS[48]*STATES[23];
RATES[23] = ( CONSTANTS[47]*pow(STATES[17], CONSTANTS[45])*STATES[25] - ( CONSTANTS[48]*STATES[23]+ CONSTANTS[49]*pow(STATES[17], CONSTANTS[46])*STATES[23]+ CONSTANTS[51]*STATES[23]))+ CONSTANTS[50]*STATES[24]+ CONSTANTS[52]*STATES[26];
RATES[24] = CONSTANTS[49]*pow(STATES[17], CONSTANTS[46])*STATES[23] - CONSTANTS[50]*STATES[24];
RATES[26] = CONSTANTS[51]*STATES[23] - CONSTANTS[52]*STATES[26];
ALGEBRAIC[0] = ( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300)));
ALGEBRAIC[7] = 0.0800000*exp(- STATES[0]/11.0000);
RATES[2] = ALGEBRAIC[0]*(1.00000 - STATES[2]) - ALGEBRAIC[7]*STATES[2];
ALGEBRAIC[1] = (STATES[0]<- 40.0000 ? 0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[8] = (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[8]*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[9] = (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[9]*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[11] = 20.0000+600.000/(1.00000+exp((STATES[0]+30.0000)/9.50000));
RATES[20] = (ALGEBRAIC[4] - STATES[20])/ALGEBRAIC[11];
ALGEBRAIC[5] = ( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)));
ALGEBRAIC[12] = ( 0.000131000*(STATES[0]+30.0000))/(- 1.00000+exp( 0.0687000*(STATES[0]+30.0000)));
RATES[21] = ALGEBRAIC[5]*(1.00000 - STATES[21]) - ALGEBRAIC[12]*STATES[21];
ALGEBRAIC[3] = 0.400000*exp((STATES[0]+12.0000)/10.0000);
ALGEBRAIC[10] = 0.0500000*exp((STATES[0]+12.0000)/- 13.0000);
ALGEBRAIC[18] = ( 0.187500*STATES[17])/1.00000;
RATES[7] = ( ALGEBRAIC[10]*STATES[8]+ CONSTANTS[20]*STATES[12]) - ( 4.00000*ALGEBRAIC[3]+ALGEBRAIC[18])*STATES[7];
RATES[8] = ( 4.00000*ALGEBRAIC[3]*STATES[7]+ 2.00000*ALGEBRAIC[10]*STATES[9]+ (CONSTANTS[20]/CONSTANTS[15])*STATES[13]) - (ALGEBRAIC[10]+ 3.00000*ALGEBRAIC[3]+ ALGEBRAIC[18]*CONSTANTS[14])*STATES[8];
RATES[9] = ( 3.00000*ALGEBRAIC[3]*STATES[8]+ 3.00000*ALGEBRAIC[10]*STATES[10]+ (CONSTANTS[20]/pow(CONSTANTS[15], 2.00000))*STATES[14]) - ( ALGEBRAIC[10]*2.00000+ 2.00000*ALGEBRAIC[3]+ ALGEBRAIC[18]*pow(CONSTANTS[14], 2.00000))*STATES[9];
RATES[10] = ( 2.00000*ALGEBRAIC[3]*STATES[9]+ 4.00000*ALGEBRAIC[10]*STATES[11]+ (CONSTANTS[20]/pow(CONSTANTS[15], 3.00000))*STATES[15]) - ( ALGEBRAIC[10]*3.00000+ALGEBRAIC[3]+ ALGEBRAIC[18]*pow(CONSTANTS[14], 3.00000))*STATES[10];
RATES[11] = ( ALGEBRAIC[3]*STATES[10]+ CONSTANTS[16]*STATES[5]+ (CONSTANTS[20]/pow(CONSTANTS[15], 4.00000))*STATES[16]) - ( ALGEBRAIC[10]*4.00000+CONSTANTS[17]+ ALGEBRAIC[18]*pow(CONSTANTS[14], 4.00000))*STATES[11];
ALGEBRAIC[14] = ALGEBRAIC[3]*CONSTANTS[14];
ALGEBRAIC[16] = ALGEBRAIC[10]/CONSTANTS[15];
RATES[12] = ( ALGEBRAIC[16]*STATES[13]+ ALGEBRAIC[18]*STATES[7]) - ( 4.00000*ALGEBRAIC[14]+CONSTANTS[20])*STATES[12];
RATES[13] = ( 4.00000*ALGEBRAIC[14]*STATES[12]+ 2.00000*ALGEBRAIC[16]*STATES[14]+ ALGEBRAIC[18]*CONSTANTS[14]*STATES[8]) - (ALGEBRAIC[16]+ 3.00000*ALGEBRAIC[14]+CONSTANTS[20]/CONSTANTS[15])*STATES[13];
RATES[14] = ( 3.00000*ALGEBRAIC[14]*STATES[13]+ 3.00000*ALGEBRAIC[16]*STATES[15]+ ALGEBRAIC[18]*pow(CONSTANTS[14], 2.00000)*STATES[9]) - ( ALGEBRAIC[16]*2.00000+ 2.00000*ALGEBRAIC[14]+CONSTANTS[20]/pow(CONSTANTS[15], 2.00000))*STATES[14];
RATES[15] = ( 2.00000*ALGEBRAIC[14]*STATES[14]+ 4.00000*ALGEBRAIC[16]*STATES[16]+ ALGEBRAIC[18]*pow(CONSTANTS[14], 3.00000)*STATES[10]) - ( ALGEBRAIC[16]*3.00000+ALGEBRAIC[14]+CONSTANTS[20]/pow(CONSTANTS[15], 3.00000))*STATES[15];
RATES[16] = ( ALGEBRAIC[14]*STATES[15]+ CONSTANTS[18]*STATES[6]+ ALGEBRAIC[18]*pow(CONSTANTS[14], 4.00000)*STATES[11]) - ( ALGEBRAIC[16]*4.00000+CONSTANTS[19]+CONSTANTS[20]/pow(CONSTANTS[15], 4.00000))*STATES[16];
ALGEBRAIC[46] = CONSTANTS[43]*(STATES[28] - STATES[22]);
ALGEBRAIC[48] = ( CONSTANTS[44]*pow(STATES[22], 2.00000))/(pow(CONSTANTS[58], 2.00000)+pow(STATES[22], 2.00000));
ALGEBRAIC[50] = (STATES[28] - STATES[27])/CONSTANTS[57];
RATES[28] = ( (ALGEBRAIC[48] - ALGEBRAIC[46])*CONSTANTS[41])/CONSTANTS[68] - ( ALGEBRAIC[50]*CONSTANTS[69])/CONSTANTS[68];
ALGEBRAIC[55] = CONSTANTS[53]*STATES[22]*(CONSTANTS[62] - STATES[29]) - CONSTANTS[54]*STATES[29];
RATES[29] = ALGEBRAIC[55];
ALGEBRAIC[17] = ( (( CONSTANTS[11]*4.00000*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.00100000*exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.341000*CONSTANTS[21]))/(exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[19] = ALGEBRAIC[17]*STATES[20]*(STATES[5]+STATES[6]);
ALGEBRAIC[53] = 1.00000/(1.00000+( CONSTANTS[65]*CONSTANTS[59])/pow(CONSTANTS[59]+STATES[17], 2.00000));
ALGEBRAIC[42] = STATES[23]+STATES[24];
ALGEBRAIC[44] = CONSTANTS[42]*ALGEBRAIC[42]*(STATES[27] - STATES[17]);
ALGEBRAIC[51] = (STATES[17] - STATES[22])/CONSTANTS[61];
RATES[17] = ALGEBRAIC[53]*((( ALGEBRAIC[44]*CONSTANTS[69])/CONSTANTS[67] - ( ALGEBRAIC[51]*CONSTANTS[41])/CONSTANTS[67]) - ( ALGEBRAIC[19]*CONSTANTS[40])/( 2.00000*CONSTANTS[67]*CONSTANTS[2]));
ALGEBRAIC[54] = 1.00000/(1.00000+( CONSTANTS[64]*CONSTANTS[60])/pow(CONSTANTS[60]+STATES[27], 2.00000));
RATES[27] = ALGEBRAIC[54]*(ALGEBRAIC[50] - ALGEBRAIC[44]);
ALGEBRAIC[13] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[10]/STATES[1]);
ALGEBRAIC[15] = CONSTANTS[9]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*(STATES[0] - ALGEBRAIC[13]);
ALGEBRAIC[35] = (( (( (( CONSTANTS[26]*1.00000)/(pow(CONSTANTS[27], 3.00000)+pow(CONSTANTS[10], 3.00000)))*1.00000)/(CONSTANTS[28]+CONSTANTS[21]))*1.00000)/(1.00000+ CONSTANTS[29]*exp(( (CONSTANTS[30] - 1.00000)*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*( exp(( CONSTANTS[30]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[1], 3.00000)*CONSTANTS[21] - exp(( (CONSTANTS[30] - 1.00000)*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(CONSTANTS[10], 3.00000)*STATES[22]);
ALGEBRAIC[37] = ALGEBRAIC[13];
ALGEBRAIC[38] = CONSTANTS[33]*(STATES[0] - ALGEBRAIC[37]);
ALGEBRAIC[41] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ 0.0365000*CONSTANTS[66]*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])));
ALGEBRAIC[43] = ( (( CONSTANTS[35]*ALGEBRAIC[41]*1.00000)/(1.00000+pow(CONSTANTS[36]/STATES[1], 1.50000)))*STATES[18])/(STATES[18]+CONSTANTS[37]);
ALGEBRAIC[45] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log((STATES[18]+CONSTANTS[10])/(STATES[19]+STATES[1]));
ALGEBRAIC[47] = STATES[0] - ALGEBRAIC[45];
ALGEBRAIC[49] = ( (( CONSTANTS[39]*pow(1.00000, 2.00000)*ALGEBRAIC[47]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[1]*exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*CONSTANTS[10]))/(exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[52] = ( ALGEBRAIC[49]*1.00000)/(1.00000+pow(CONSTANTS[38]/STATES[22], 3.00000));
RATES[1] = ( - (ALGEBRAIC[15]+ALGEBRAIC[38]+ALGEBRAIC[52]+ ALGEBRAIC[35]*3.00000+ ALGEBRAIC[43]*3.00000)*CONSTANTS[40])/( CONSTANTS[41]*CONSTANTS[2]);
ALGEBRAIC[57] = CONSTANTS[55]*STATES[22]*(CONSTANTS[63] - STATES[30]) - CONSTANTS[56]*STATES[30];
RATES[30] = ALGEBRAIC[57];
ALGEBRAIC[20] = CONSTANTS[12]/(1.00000+ALGEBRAIC[17]/CONSTANTS[13]);
ALGEBRAIC[21] = ( (( ALGEBRAIC[20]*STATES[20]*(STATES[5]+STATES[6])*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( STATES[19]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - STATES[18]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[22] = CONSTANTS[22]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[23] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log((STATES[18]+ CONSTANTS[23]*CONSTANTS[10])/(STATES[19]+ CONSTANTS[23]*STATES[1]));
ALGEBRAIC[24] = 1.00000/(1.00000+exp((STATES[0] - 56.2600)/32.1000));
ALGEBRAIC[25] = ALGEBRAIC[22]*ALGEBRAIC[24]*pow(STATES[21], 2.00000)*(STATES[0] - ALGEBRAIC[23]);
ALGEBRAIC[27] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(STATES[18]/STATES[19]);
ALGEBRAIC[26] = CONSTANTS[24]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[28] = 1.02000/(1.00000+exp( 0.238500*((STATES[0] - ALGEBRAIC[27]) - 59.2150)));
ALGEBRAIC[29] = ( 0.491240*(exp( 0.0803200*((STATES[0]+5.47600) - ALGEBRAIC[27]))+exp( 0.0617500*(STATES[0] - (ALGEBRAIC[27]+594.310)))))/(1.00000+exp( - 0.514300*((STATES[0] - ALGEBRAIC[27])+4.75300)));
ALGEBRAIC[30] = ALGEBRAIC[28]/(ALGEBRAIC[28]+ALGEBRAIC[29]);
ALGEBRAIC[31] = ALGEBRAIC[26]*ALGEBRAIC[30]*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[32] = ALGEBRAIC[27];
ALGEBRAIC[33] = 1.00000/(1.00000+exp((7.48800 - STATES[0])/5.98000));
ALGEBRAIC[34] = CONSTANTS[25]*ALGEBRAIC[33]*(STATES[0] - ALGEBRAIC[32]);
ALGEBRAIC[56] = ( (( CONSTANTS[39]*pow(1.00000, 2.00000)*ALGEBRAIC[47]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[19]*exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*STATES[18]))/(exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[58] = ( ALGEBRAIC[56]*1.00000)/(1.00000+pow(CONSTANTS[38]/STATES[22], 3.00000));
RATES[19] = ( - (ALGEBRAIC[21]+ALGEBRAIC[25]+ALGEBRAIC[31]+ALGEBRAIC[34]+ALGEBRAIC[58]+ - ALGEBRAIC[43]*2.00000)*CONSTANTS[40])/( CONSTANTS[41]*CONSTANTS[2]);
RATES[18] = ( (ALGEBRAIC[21]+ALGEBRAIC[25]+ALGEBRAIC[31]+ALGEBRAIC[34]+ALGEBRAIC[58]+ - ALGEBRAIC[43]*2.00000)*CONSTANTS[40])/( CONSTANTS[41]*CONSTANTS[2]);
ALGEBRAIC[36] = ( CONSTANTS[32]*STATES[22])/(CONSTANTS[31]+STATES[22]);
ALGEBRAIC[39] = (( CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]))*log(CONSTANTS[21]/STATES[22]);
ALGEBRAIC[40] = CONSTANTS[34]*(STATES[0] - ALGEBRAIC[39]);
ALGEBRAIC[60] = ALGEBRAIC[52]+ALGEBRAIC[58];
ALGEBRAIC[6] = (VOI>=CONSTANTS[4]&&VOI<=CONSTANTS[5]&&(VOI - CONSTANTS[4]) - floor((VOI - CONSTANTS[4])/CONSTANTS[6])*CONSTANTS[6]<=CONSTANTS[7] ? CONSTANTS[8] : 0.00000);
RATES[0] = (ALGEBRAIC[6] - (ALGEBRAIC[15]+ALGEBRAIC[19]+ALGEBRAIC[21]+ALGEBRAIC[25]+ALGEBRAIC[35]+ALGEBRAIC[31]+ALGEBRAIC[34]+ALGEBRAIC[36]+ALGEBRAIC[38]+ALGEBRAIC[40]+ALGEBRAIC[43]+ALGEBRAIC[60]))/CONSTANTS[3];
ALGEBRAIC[61] = 1.00000/(1.00000+( CONSTANTS[65]*CONSTANTS[59])/pow(CONSTANTS[59]+STATES[22], 2.00000));
ALGEBRAIC[59] = ALGEBRAIC[55]+ALGEBRAIC[57];
RATES[22] = ALGEBRAIC[61]*((ALGEBRAIC[46]+ALGEBRAIC[51]) - (ALGEBRAIC[48]+ALGEBRAIC[59]+( ((ALGEBRAIC[40] - 2.00000*ALGEBRAIC[35])+ALGEBRAIC[36])*CONSTANTS[40])/( 2.00000*CONSTANTS[41]*CONSTANTS[2])));
}
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[7] = 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[8] = (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[9] = (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[11] = 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[12] = ( 0.000131000*(STATES[0]+30.0000))/(- 1.00000+exp( 0.0687000*(STATES[0]+30.0000)));
ALGEBRAIC[3] = 0.400000*exp((STATES[0]+12.0000)/10.0000);
ALGEBRAIC[10] = 0.0500000*exp((STATES[0]+12.0000)/- 13.0000);
ALGEBRAIC[18] = ( 0.187500*STATES[17])/1.00000;
ALGEBRAIC[14] = ALGEBRAIC[3]*CONSTANTS[14];
ALGEBRAIC[16] = ALGEBRAIC[10]/CONSTANTS[15];
ALGEBRAIC[46] = CONSTANTS[43]*(STATES[28] - STATES[22]);
ALGEBRAIC[48] = ( CONSTANTS[44]*pow(STATES[22], 2.00000))/(pow(CONSTANTS[58], 2.00000)+pow(STATES[22], 2.00000));
ALGEBRAIC[50] = (STATES[28] - STATES[27])/CONSTANTS[57];
ALGEBRAIC[55] = CONSTANTS[53]*STATES[22]*(CONSTANTS[62] - STATES[29]) - CONSTANTS[54]*STATES[29];
ALGEBRAIC[17] = ( (( CONSTANTS[11]*4.00000*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.00100000*exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.341000*CONSTANTS[21]))/(exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[19] = ALGEBRAIC[17]*STATES[20]*(STATES[5]+STATES[6]);
ALGEBRAIC[53] = 1.00000/(1.00000+( CONSTANTS[65]*CONSTANTS[59])/pow(CONSTANTS[59]+STATES[17], 2.00000));
ALGEBRAIC[42] = STATES[23]+STATES[24];
ALGEBRAIC[44] = CONSTANTS[42]*ALGEBRAIC[42]*(STATES[27] - STATES[17]);
ALGEBRAIC[51] = (STATES[17] - STATES[22])/CONSTANTS[61];
ALGEBRAIC[54] = 1.00000/(1.00000+( CONSTANTS[64]*CONSTANTS[60])/pow(CONSTANTS[60]+STATES[27], 2.00000));
ALGEBRAIC[13] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[10]/STATES[1]);
ALGEBRAIC[15] = CONSTANTS[9]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*(STATES[0] - ALGEBRAIC[13]);
ALGEBRAIC[35] = (( (( (( CONSTANTS[26]*1.00000)/(pow(CONSTANTS[27], 3.00000)+pow(CONSTANTS[10], 3.00000)))*1.00000)/(CONSTANTS[28]+CONSTANTS[21]))*1.00000)/(1.00000+ CONSTANTS[29]*exp(( (CONSTANTS[30] - 1.00000)*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*( exp(( CONSTANTS[30]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[1], 3.00000)*CONSTANTS[21] - exp(( (CONSTANTS[30] - 1.00000)*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(CONSTANTS[10], 3.00000)*STATES[22]);
ALGEBRAIC[37] = ALGEBRAIC[13];
ALGEBRAIC[38] = CONSTANTS[33]*(STATES[0] - ALGEBRAIC[37]);
ALGEBRAIC[41] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ 0.0365000*CONSTANTS[66]*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])));
ALGEBRAIC[43] = ( (( CONSTANTS[35]*ALGEBRAIC[41]*1.00000)/(1.00000+pow(CONSTANTS[36]/STATES[1], 1.50000)))*STATES[18])/(STATES[18]+CONSTANTS[37]);
ALGEBRAIC[45] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log((STATES[18]+CONSTANTS[10])/(STATES[19]+STATES[1]));
ALGEBRAIC[47] = STATES[0] - ALGEBRAIC[45];
ALGEBRAIC[49] = ( (( CONSTANTS[39]*pow(1.00000, 2.00000)*ALGEBRAIC[47]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[1]*exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*CONSTANTS[10]))/(exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[52] = ( ALGEBRAIC[49]*1.00000)/(1.00000+pow(CONSTANTS[38]/STATES[22], 3.00000));
ALGEBRAIC[57] = CONSTANTS[55]*STATES[22]*(CONSTANTS[63] - STATES[30]) - CONSTANTS[56]*STATES[30];
ALGEBRAIC[20] = CONSTANTS[12]/(1.00000+ALGEBRAIC[17]/CONSTANTS[13]);
ALGEBRAIC[21] = ( (( ALGEBRAIC[20]*STATES[20]*(STATES[5]+STATES[6])*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( STATES[19]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - STATES[18]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[22] = CONSTANTS[22]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[23] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log((STATES[18]+ CONSTANTS[23]*CONSTANTS[10])/(STATES[19]+ CONSTANTS[23]*STATES[1]));
ALGEBRAIC[24] = 1.00000/(1.00000+exp((STATES[0] - 56.2600)/32.1000));
ALGEBRAIC[25] = ALGEBRAIC[22]*ALGEBRAIC[24]*pow(STATES[21], 2.00000)*(STATES[0] - ALGEBRAIC[23]);
ALGEBRAIC[27] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(STATES[18]/STATES[19]);
ALGEBRAIC[26] = CONSTANTS[24]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[28] = 1.02000/(1.00000+exp( 0.238500*((STATES[0] - ALGEBRAIC[27]) - 59.2150)));
ALGEBRAIC[29] = ( 0.491240*(exp( 0.0803200*((STATES[0]+5.47600) - ALGEBRAIC[27]))+exp( 0.0617500*(STATES[0] - (ALGEBRAIC[27]+594.310)))))/(1.00000+exp( - 0.514300*((STATES[0] - ALGEBRAIC[27])+4.75300)));
ALGEBRAIC[30] = ALGEBRAIC[28]/(ALGEBRAIC[28]+ALGEBRAIC[29]);
ALGEBRAIC[31] = ALGEBRAIC[26]*ALGEBRAIC[30]*(STATES[0] - ALGEBRAIC[27]);
ALGEBRAIC[32] = ALGEBRAIC[27];
ALGEBRAIC[33] = 1.00000/(1.00000+exp((7.48800 - STATES[0])/5.98000));
ALGEBRAIC[34] = CONSTANTS[25]*ALGEBRAIC[33]*(STATES[0] - ALGEBRAIC[32]);
ALGEBRAIC[56] = ( (( CONSTANTS[39]*pow(1.00000, 2.00000)*ALGEBRAIC[47]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[19]*exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*STATES[18]))/(exp(( ALGEBRAIC[47]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000);
ALGEBRAIC[58] = ( ALGEBRAIC[56]*1.00000)/(1.00000+pow(CONSTANTS[38]/STATES[22], 3.00000));
ALGEBRAIC[36] = ( CONSTANTS[32]*STATES[22])/(CONSTANTS[31]+STATES[22]);
ALGEBRAIC[39] = (( CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]))*log(CONSTANTS[21]/STATES[22]);
ALGEBRAIC[40] = CONSTANTS[34]*(STATES[0] - ALGEBRAIC[39]);
ALGEBRAIC[60] = ALGEBRAIC[52]+ALGEBRAIC[58];
ALGEBRAIC[6] = (VOI>=CONSTANTS[4]&&VOI<=CONSTANTS[5]&&(VOI - CONSTANTS[4]) - floor((VOI - CONSTANTS[4])/CONSTANTS[6])*CONSTANTS[6]<=CONSTANTS[7] ? CONSTANTS[8] : 0.00000);
ALGEBRAIC[61] = 1.00000/(1.00000+( CONSTANTS[65]*CONSTANTS[59])/pow(CONSTANTS[59]+STATES[22], 2.00000));
ALGEBRAIC[59] = ALGEBRAIC[55]+ALGEBRAIC[57];
}