/*
   There are a total of 42 entries in the algebraic variable array.
   There are a total of 27 entries in each of the rate and state variable arrays.
   There are a total of 64 entries in the constant variable array.
 */
/*
 * VOI is time in component Environment (second).
 * CONSTANTS[0] is O2_x in component Environment (molar).
 * CONSTANTS[1] is X_AtC in component Environment (mole_per_second_per_l_cell).
 * CONSTANTS[2] is RT in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[3] is F in component fixed_parameters (kilojoule_per_mole_per_millivolt).
 * CONSTANTS[4] is dG_C1o in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[5] is dG_C3o in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[6] is dG_C4o in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[13] is dG_C1op in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[14] is dG_C3op in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[15] is dG_C4op in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[7] is dG_F1o in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[11] is dG_H in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[50] is k_dHPI in component fixed_parameters (molar).
 * CONSTANTS[51] is k_dHatp in component fixed_parameters (molar).
 * CONSTANTS[52] is k_dHadp in component fixed_parameters (molar).
 * CONSTANTS[8] is K_DT in component fixed_parameters (molar).
 * CONSTANTS[9] is K_DD in component fixed_parameters (molar).
 * CONSTANTS[10] is K_AK in component fixed_parameters (dimensionless).
 * CONSTANTS[11] is K_CK in component fixed_parameters (per_molar).
 * CONSTANTS[56] is W_m in component fixed_parameters (l_water_per_l_mito).
 * CONSTANTS[59] is W_c in component fixed_parameters (l_water_per_l_cyto).
 * CONSTANTS[57] is W_x in component fixed_parameters (l_water_per_l_mito).
 * CONSTANTS[58] is W_i in component fixed_parameters (l_water_per_l_mito).
 * CONSTANTS[12] is V_cyto in component fixed_parameters (cyto_per_cell).
 * CONSTANTS[13] is V_mito in component fixed_parameters (mito_per_cell).
 * CONSTANTS[61] is Rm_cyto in component fixed_parameters (mito_per_cyto).
 * CONSTANTS[54] is Rm_cell in component fixed_parameters (mito_per_cell).
 * CONSTANTS[62] is Rc_cell in component fixed_parameters (cyto_per_cell).
 * CONSTANTS[14] is n_A in component fixed_parameters (dimensionless).
 * CONSTANTS[15] is C_tot in component fixed_parameters (molar).
 * CONSTANTS[16] is CR_tot in component fixed_parameters (molar).
 * CONSTANTS[17] is Q_tot in component fixed_parameters (molar).
 * CONSTANTS[18] is NAD_tot in component fixed_parameters (molar).
 * CONSTANTS[19] is pH_C in component fixed_parameters (dimensionless).
 * CONSTANTS[53] is H_c in component fixed_parameters (molar).
 * CONSTANTS[20] is K_c in component fixed_parameters (molar).
 * CONSTANTS[60] is H_i in component fixed_parameters (molar).
 * CONSTANTS[55] is K_i in component fixed_parameters (molar).
 * ALGEBRAIC[9] is H_x in component dH_x_dt (molar).
 * STATES[0] is deltaPsi in component ddeltaPsi_dt (millivolt).
 * CONSTANTS[21] is k_PI1 in component adjustable_parameters (molar).
 * CONSTANTS[22] is k_PI2 in component adjustable_parameters (molar).
 * CONSTANTS[23] is k_PI3 in component adjustable_parameters (molar).
 * CONSTANTS[24] is k_PI4 in component adjustable_parameters (molar).
 * CONSTANTS[25] is k_PIH in component adjustable_parameters (molar).
 * CONSTANTS[26] is r in component adjustable_parameters (dimensionless).
 * CONSTANTS[27] is X_DH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
 * CONSTANTS[28] is X_C1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[29] is X_C3 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_half_molar).
 * CONSTANTS[30] is X_C4 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
 * CONSTANTS[31] is X_F1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[32] is X_ANT in component adjustable_parameters (mole_per_second_per_l_mito).
 * CONSTANTS[33] is X_PI1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
 * CONSTANTS[34] is X_KH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[35] is X_Hle in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
 * CONSTANTS[36] is X_K in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
 * CONSTANTS[37] is k_mADP in component adjustable_parameters (molar).
 * CONSTANTS[38] is X_AK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[39] is X_CK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[40] is X_MgA in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[41] is k_O2 in component adjustable_parameters (molar).
 * CONSTANTS[42] is x_buff in component adjustable_parameters (per_molar).
 * CONSTANTS[43] is C_IM in component adjustable_parameters (mole_per_l_mito_per_millivolt).
 * CONSTANTS[44] is X_A in component adjustable_parameters (micron_per_second).
 * CONSTANTS[45] is X_PI2 in component adjustable_parameters (micron_per_second).
 * CONSTANTS[46] is gamma in component adjustable_parameters (per_micron).
 * ALGEBRAIC[0] is NAD_x in component dNAD_x_dt (molar).
 * STATES[1] is NADH_x in component dNADH_x_dt (molar).
 * ALGEBRAIC[2] is QH2_x in component dQH2_x_dt (molar).
 * STATES[2] is Q_x in component dQ_x_dt (molar).
 * ALGEBRAIC[4] is Cred_x in component dCred_x_dt (molar).
 * STATES[3] is Cox_x in component dCox_x_dt (molar).
 * ALGEBRAIC[6] is ATP_fx in component dATP_fx_dt (molar).
 * STATES[4] is ATP_x in component dATP_x_dt (molar).
 * STATES[5] is ATP_mx in component dATP_mx_dt (molar).
 * ALGEBRAIC[8] is ADP_fx in component dADP_fx_dt (molar).
 * STATES[6] is ADP_x in component dADP_x_dt (molar).
 * STATES[7] is ADP_mx in component dADP_mx_dt (molar).
 * STATES[8] is pH_x in component dH_x_dt (dimensionless).
 * ALGEBRAIC[20] is J_DH in component J_DH (mole_per_second_per_l_mito).
 * ALGEBRAIC[22] is J_C1 in component J_C1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[23] is J_C3 in component J_C3 (mole_per_second_per_l_mito).
 * ALGEBRAIC[24] is J_C4 in component J_C4 (mole_per_second_per_l_mito).
 * ALGEBRAIC[25] is J_F1 in component J_F1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[33] is J_PI1 in component J_PI1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[39] is J_KH in component J_KH (mole_per_second_per_l_mito).
 * ALGEBRAIC[35] is J_Hle in component J_Hle (mole_per_second_per_l_mito).
 * STATES[9] is K_x in component dK_x_dt (molar).
 * ALGEBRAIC[41] is J_K in component J_K (mole_per_second_per_l_mito).
 * STATES[10] is Mg_x in component dMg_x_dt (molar).
 * ALGEBRAIC[10] is J_MgATP_x in component J_MgATP_x (mole_per_second_per_l_mito).
 * ALGEBRAIC[12] is J_MgADP_x in component J_MgADP_x (mole_per_second_per_l_mito).
 * ALGEBRAIC[28] is J_ANT in component J_ANT (mole_per_second_per_l_mito).
 * STATES[11] is PI_x in component dPI_x_dt (molar).
 * ALGEBRAIC[16] is ATP_fi in component dATP_fi_dt (molar).
 * STATES[12] is ATP_i in component dATP_i_dt (molar).
 * STATES[13] is ATP_mi in component dATP_mi_dt (molar).
 * ALGEBRAIC[17] is ADP_fi in component dADP_fi_dt (molar).
 * STATES[14] is ADP_i in component dADP_i_dt (molar).
 * STATES[15] is ADP_mi in component dADP_mi_dt (molar).
 * ALGEBRAIC[34] is J_ATP in component J_ATP (mole_per_second_per_l_mito).
 * ALGEBRAIC[36] is J_AKi in component J_AKi (mole_per_second_per_l_mito).
 * ALGEBRAIC[32] is J_ADP in component J_ADP (mole_per_second_per_l_mito).
 * STATES[16] is AMP_i in component dAMP_i_dt (molar).
 * ALGEBRAIC[30] is J_AMP in component J_AMP (mole_per_second_per_l_mito).
 * ALGEBRAIC[19] is J_MgATP_i in component J_MgATP_i (mole_per_second_per_l_mito).
 * ALGEBRAIC[21] is J_MgADP_i in component J_MgADP_i (mole_per_second_per_l_mito).
 * STATES[17] is PI_i in component dPI_i_dt (molar).
 * ALGEBRAIC[37] is J_PI2 in component J_PI2 (mole_per_second_per_l_mito).
 * STATES[18] is Mg_i in component dMg_i_dt (molar).
 * STATES[19] is ATP_c in component dATP_c_dt (molar).
 * ALGEBRAIC[38] is J_AKc in component J_AKc (mole_per_second_per_l_cyto).
 * CONSTANTS[63] is J_AtC in component J_AtC (mole_per_second_per_l_cell).
 * ALGEBRAIC[40] is J_CKc in component J_CKc (mole_per_second_per_l_cyto).
 * STATES[20] is ADP_c in component dADP_c_dt (molar).
 * STATES[21] is AMP_c in component dAMP_c_dt (molar).
 * STATES[22] is ATP_mc in component dATP_mc_dt (molar).
 * ALGEBRAIC[5] is J_MgATP_c in component J_MgATP_c (mole_per_second_per_l_cyto).
 * STATES[23] is ADP_mc in component dADP_mc_dt (molar).
 * ALGEBRAIC[7] is J_MgADP_c in component J_MgADP_c (mole_per_second_per_l_cyto).
 * STATES[24] is PI_c in component dPI_c_dt (molar).
 * STATES[25] is Mg_c in component dMg_c_dt (molar).
 * STATES[26] is PCr_c in component dPCr_c_dt (molar).
 * ALGEBRAIC[1] is ATP_fc in component dATP_fc_dt (molar).
 * ALGEBRAIC[3] is ADP_fc in component dADP_fc_dt (molar).
 * ALGEBRAIC[18] is Cr_c in component dCr_c_dt (molar).
 * CONSTANTS[47] is mincon in component J_ANT (molar).
 * ALGEBRAIC[26] is Psi_x in component J_ANT (millivolt).
 * ALGEBRAIC[27] is Psi_i in component J_ANT (millivolt).
 * ALGEBRAIC[29] is H2PIi in component J_PI1 (molar).
 * ALGEBRAIC[31] is H2PIx in component J_PI1 (molar).
 * CONSTANTS[48] is mincond in component J_Hle (millivolt).
 * CONSTANTS[49] is mincond in component J_K (millivolt).
 * RATES[8] is d/dt pH_x in component dH_x_dt (dimensionless).
 * RATES[9] is d/dt K_x in component dK_x_dt (molar).
 * RATES[10] is d/dt Mg_x in component dMg_x_dt (molar).
 * RATES[1] is d/dt NADH_x in component dNADH_x_dt (molar).
 * RATES[2] is d/dt Q_x in component dQ_x_dt (molar).
 * RATES[3] is d/dt Cox_x in component dCox_x_dt (molar).
 * RATES[4] is d/dt ATP_x in component dATP_x_dt (molar).
 * RATES[6] is d/dt ADP_x in component dADP_x_dt (molar).
 * RATES[5] is d/dt ATP_mx in component dATP_mx_dt (molar).
 * RATES[7] is d/dt ADP_mx in component dADP_mx_dt (molar).
 * RATES[11] is d/dt PI_x in component dPI_x_dt (molar).
 * RATES[12] is d/dt ATP_i in component dATP_i_dt (molar).
 * RATES[14] is d/dt ADP_i in component dADP_i_dt (molar).
 * RATES[16] is d/dt AMP_i in component dAMP_i_dt (molar).
 * RATES[13] is d/dt ATP_mi in component dATP_mi_dt (molar).
 * RATES[15] is d/dt ADP_mi in component dADP_mi_dt (molar).
 * RATES[17] is d/dt PI_i in component dPI_i_dt (molar).
 * RATES[18] is d/dt Mg_i in component dMg_i_dt (molar).
 * RATES[19] is d/dt ATP_c in component dATP_c_dt (molar).
 * RATES[20] is d/dt ADP_c in component dADP_c_dt (molar).
 * RATES[21] is d/dt AMP_c in component dAMP_c_dt (molar).
 * RATES[22] is d/dt ATP_mc in component dATP_mc_dt (molar).
 * RATES[23] is d/dt ADP_mc in component dADP_mc_dt (molar).
 * RATES[24] is d/dt PI_c in component dPI_c_dt (molar).
 * RATES[25] is d/dt Mg_c in component dMg_c_dt (molar).
 * RATES[26] is d/dt PCr_c in component dPCr_c_dt (molar).
 * RATES[0] is d/dt deltaPsi in component ddeltaPsi_dt (millivolt).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 3.48e-5;
CONSTANTS[1] = 0e-4;
CONSTANTS[2] = 2.5775;
CONSTANTS[3] = 0.096484;
CONSTANTS[4] = -69.37;
CONSTANTS[5] = -32.53;
CONSTANTS[6] = -122.94;
CONSTANTS[7] = 36.03;
CONSTANTS[8] = 24e-6;
CONSTANTS[9] = 347e-6;
CONSTANTS[10] = 0.4331;
CONSTANTS[11] = 1.66e9;
CONSTANTS[12] = 0.894;
CONSTANTS[13] = 0.056;
CONSTANTS[14] = 3;
CONSTANTS[15] = 2.7e-3;
CONSTANTS[16] = 42.7e-3;
CONSTANTS[17] = 1.35e-3;
CONSTANTS[18] = 2.97e-3;
CONSTANTS[19] = 7;
CONSTANTS[20] = 0.15;
STATES[0] = 166.67;
CONSTANTS[21] = 0.1553e-3;
CONSTANTS[22] = 0.8222e-3;
CONSTANTS[23] = 0.3601e-3;
CONSTANTS[24] = 5.651e-3;
CONSTANTS[25] = 2.542e-4;
CONSTANTS[26] = 4.559;
CONSTANTS[27] = 0.0866;
CONSTANTS[28] = 4.405e3;
CONSTANTS[29] = 4.887;
CONSTANTS[30] = 6.766e-5;
CONSTANTS[31] = 1e3;
CONSTANTS[32] = 8.123e-3;
CONSTANTS[33] = 3.85e5;
CONSTANTS[34] = 5.651e7;
CONSTANTS[35] = 200;
CONSTANTS[36] = 0;
CONSTANTS[37] = 3.5e-6;
CONSTANTS[38] = 1e7;
CONSTANTS[39] = 1e7;
CONSTANTS[40] = 1e7;
CONSTANTS[41] = 1.2e-4;
CONSTANTS[42] = 100;
CONSTANTS[43] = 6.75e-6;
CONSTANTS[44] = 85;
CONSTANTS[45] = 327;
CONSTANTS[46] = 5.99;
STATES[1] = 0.0015723;
STATES[2] = 6.75e-4;
STATES[3] = 1.35e-3;
STATES[4] = 0.0026657;
STATES[5] = 0.0026046;
STATES[6] = 0.0073343;
STATES[7] = 0.0054765;
STATES[8] = 7;
STATES[9] = 0.14661;
STATES[10] = 1.0229e-3;
STATES[11] = 1.72e-4;
STATES[12] = 0.0065339;
STATES[13] = 0.0063812;
STATES[14] = 6.5773e-5;
STATES[15] = 4.8866e-5;
STATES[16] = 2.8837e-7;
STATES[17] = 1.72e-4;
STATES[18] = 0.0010029;
STATES[19] = 0.0081312;
STATES[20] = 6.85e-5;
STATES[21] = 3.0911e-7;
STATES[22] = 0.0079786;
STATES[23] = 5.1958e-5;
STATES[24] = 8.7702e-3;
STATES[25] = 0.001003;
STATES[26] = 23.589e-3;
CONSTANTS[47] = 1e-9;
CONSTANTS[48] = 1e-6;
CONSTANTS[49] = 1e-6;
CONSTANTS[50] =  1.00000*pow(10.0000, - 6.75000);
CONSTANTS[51] =  1.00000*pow(10.0000, - 6.48000);
CONSTANTS[52] =  1.00000*pow(10.0000, - 6.29000);
CONSTANTS[53] =  1.00000*pow(10.0000, - CONSTANTS[19]);
CONSTANTS[54] = CONSTANTS[13];
CONSTANTS[55] = CONSTANTS[20];
CONSTANTS[56] =  0.664000*1.09000;
CONSTANTS[57] =  0.900000*CONSTANTS[56];
CONSTANTS[58] =  0.100000*CONSTANTS[56];
CONSTANTS[59] =  0.807000*1.04400;
CONSTANTS[60] = CONSTANTS[53];
CONSTANTS[61] = CONSTANTS[13]/CONSTANTS[12];
CONSTANTS[62] = CONSTANTS[12];
CONSTANTS[63] = CONSTANTS[1];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = STATES[19] - STATES[22];
ALGEBRAIC[5] =  CONSTANTS[40]*1.00000*( ALGEBRAIC[1]*STATES[25] -  CONSTANTS[8]*STATES[22]);
RATES[22] = ALGEBRAIC[5]/CONSTANTS[59];
ALGEBRAIC[3] = STATES[20] - STATES[23];
ALGEBRAIC[7] =  CONSTANTS[40]*1.00000*( ALGEBRAIC[3]*STATES[25] -  CONSTANTS[9]*STATES[23]);
RATES[23] = ALGEBRAIC[7]/CONSTANTS[59];
RATES[25] = (- ALGEBRAIC[5] - ALGEBRAIC[7])/CONSTANTS[59];
ALGEBRAIC[6] = STATES[4] - STATES[5];
ALGEBRAIC[10] =  CONSTANTS[40]*( ALGEBRAIC[6]*STATES[10] -  CONSTANTS[8]*STATES[5]);
RATES[5] = ALGEBRAIC[10]/CONSTANTS[57];
ALGEBRAIC[8] = STATES[6] - STATES[7];
ALGEBRAIC[12] =  CONSTANTS[40]*( ALGEBRAIC[8]*STATES[10] -  CONSTANTS[9]*STATES[7]);
RATES[10] = (- ALGEBRAIC[12] - ALGEBRAIC[10])/CONSTANTS[57];
RATES[7] = ALGEBRAIC[12]/CONSTANTS[57];
ALGEBRAIC[16] = STATES[12] - STATES[13];
ALGEBRAIC[19] =  CONSTANTS[40]*( ALGEBRAIC[16]*STATES[18] -  CONSTANTS[8]*STATES[13]);
RATES[13] = ALGEBRAIC[19]/CONSTANTS[58];
ALGEBRAIC[17] = STATES[14] - STATES[15];
ALGEBRAIC[21] =  CONSTANTS[40]*( ALGEBRAIC[17]*STATES[18] -  CONSTANTS[9]*STATES[15]);
RATES[15] = ALGEBRAIC[21]/CONSTANTS[58];
RATES[18] = (- ALGEBRAIC[21] - ALGEBRAIC[19])/CONSTANTS[58];
ALGEBRAIC[0] = CONSTANTS[18] - STATES[1];
ALGEBRAIC[20] = ( CONSTANTS[27]*( CONSTANTS[26]*ALGEBRAIC[0] - STATES[1])*(1.00000+STATES[11]/CONSTANTS[21]))/(1.00000+STATES[11]/CONSTANTS[22]);
ALGEBRAIC[9] =  1.00000*pow(10.0000, - STATES[8]);
ALGEBRAIC[13] = CONSTANTS[4] -  1.00000*CONSTANTS[2]*log(ALGEBRAIC[9]/1.00000e-07);
ALGEBRAIC[11] =  CONSTANTS[3]*STATES[0]+ CONSTANTS[2]*log(CONSTANTS[60]/ALGEBRAIC[9]);
ALGEBRAIC[2] = CONSTANTS[17] - STATES[2];
ALGEBRAIC[22] =  CONSTANTS[28]*( exp(- (ALGEBRAIC[13]+ 4.00000*ALGEBRAIC[11])/CONSTANTS[2])*STATES[1]*STATES[2] -  ALGEBRAIC[0]*ALGEBRAIC[2]);
RATES[1] = (ALGEBRAIC[20] - ALGEBRAIC[22])/CONSTANTS[57];
ALGEBRAIC[14] = CONSTANTS[5]+ 2.00000*CONSTANTS[2]*log(ALGEBRAIC[9]/1.00000e-07);
ALGEBRAIC[4] = CONSTANTS[15] - STATES[3];
ALGEBRAIC[23] =  (( CONSTANTS[29]*(1.00000+STATES[11]/CONSTANTS[23]))/(1.00000+STATES[11]/CONSTANTS[24]))*( exp(- ((ALGEBRAIC[14]+ 4.00000*ALGEBRAIC[11]) -  2.00000*CONSTANTS[3]*STATES[0])/( 2.00000*CONSTANTS[2]))*STATES[3]* pow(ALGEBRAIC[2], 1.0 / 2) -  ALGEBRAIC[4]* pow(STATES[2], 1.0 / 2));
RATES[2] = - (ALGEBRAIC[22] - ALGEBRAIC[23])/CONSTANTS[57];
ALGEBRAIC[15] = CONSTANTS[6] -  2.00000*CONSTANTS[2]*log(ALGEBRAIC[9]/1.00000e-07);
ALGEBRAIC[24] = (CONSTANTS[0]<1.00000e-12 ? 0.00000 :  (( CONSTANTS[30]*CONSTANTS[0])/(CONSTANTS[0]+CONSTANTS[41]))*(ALGEBRAIC[4]/CONSTANTS[15])*( exp(- (ALGEBRAIC[15]+ 2.00000*ALGEBRAIC[11])/( 2.00000*CONSTANTS[2]))*ALGEBRAIC[4]*pow(CONSTANTS[0]/1.00000, 0.250000) -  STATES[3]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2])));
RATES[3] = - ( 2.00000*ALGEBRAIC[23] -  2.00000*ALGEBRAIC[24])/CONSTANTS[57];
ALGEBRAIC[25] =  CONSTANTS[31]*( (( exp(- (CONSTANTS[7] -  CONSTANTS[14]*ALGEBRAIC[11])/CONSTANTS[2])*CONSTANTS[9])/CONSTANTS[8])*STATES[7]*STATES[11] -  STATES[5]*1.00000);
ALGEBRAIC[26] =  - 0.650000*STATES[0];
ALGEBRAIC[27] =  0.350000*STATES[0];
ALGEBRAIC[28] = (ALGEBRAIC[17]>CONSTANTS[47]||ALGEBRAIC[16]>CONSTANTS[47] ?  CONSTANTS[32]*(ALGEBRAIC[17]/(ALGEBRAIC[17]+ ALGEBRAIC[16]*exp(( - CONSTANTS[3]*ALGEBRAIC[27])/CONSTANTS[2])) - ALGEBRAIC[8]/(ALGEBRAIC[8]+ ALGEBRAIC[6]*exp(( - CONSTANTS[3]*ALGEBRAIC[26])/CONSTANTS[2])))*(ALGEBRAIC[17]/(ALGEBRAIC[17]+CONSTANTS[37])) : 0.00000);
RATES[4] = (ALGEBRAIC[25] - ALGEBRAIC[28])/CONSTANTS[57];
RATES[6] = (- ALGEBRAIC[25]+ALGEBRAIC[28])/CONSTANTS[57];
ALGEBRAIC[29] = ( STATES[17]*CONSTANTS[60])/(CONSTANTS[60]+CONSTANTS[50]);
ALGEBRAIC[31] = ( STATES[11]*ALGEBRAIC[9])/(ALGEBRAIC[9]+CONSTANTS[50]);
ALGEBRAIC[33] = ( CONSTANTS[33]*( CONSTANTS[60]*ALGEBRAIC[29] -  ALGEBRAIC[9]*ALGEBRAIC[31]))/(ALGEBRAIC[29]+CONSTANTS[25]);
RATES[11] = (- ALGEBRAIC[25]+ALGEBRAIC[33])/CONSTANTS[57];
ALGEBRAIC[34] =  CONSTANTS[46]*CONSTANTS[44]*(STATES[19] - STATES[12]);
ALGEBRAIC[36] =  CONSTANTS[38]*( CONSTANTS[10]*STATES[14]*STATES[14] -  STATES[16]*STATES[12]);
RATES[12] = (ALGEBRAIC[34]+ALGEBRAIC[28]+ALGEBRAIC[36])/CONSTANTS[58];
ALGEBRAIC[32] =  CONSTANTS[46]*CONSTANTS[44]*(STATES[20] - STATES[14]);
RATES[14] = ((ALGEBRAIC[32] - ALGEBRAIC[28]) -  2.00000*ALGEBRAIC[36])/CONSTANTS[58];
ALGEBRAIC[30] =  CONSTANTS[46]*CONSTANTS[44]*(STATES[21] - STATES[16]);
RATES[16] = (ALGEBRAIC[30]+ALGEBRAIC[36])/CONSTANTS[58];
ALGEBRAIC[37] =  CONSTANTS[46]*CONSTANTS[45]*(STATES[24] - STATES[17]);
RATES[17] = (- ALGEBRAIC[33]+ALGEBRAIC[37])/CONSTANTS[58];
ALGEBRAIC[38] =  CONSTANTS[38]*1.00000*( CONSTANTS[10]*STATES[20]*STATES[20] -  STATES[21]*STATES[19]);
RATES[21] = ( - CONSTANTS[61]*ALGEBRAIC[30]+ALGEBRAIC[38])/CONSTANTS[59];
RATES[24] = ( - ALGEBRAIC[37]*CONSTANTS[61]+CONSTANTS[63]/CONSTANTS[62])/CONSTANTS[59];
ALGEBRAIC[39] =  CONSTANTS[34]*( CONSTANTS[55]*ALGEBRAIC[9] -  STATES[9]*CONSTANTS[60]);
ALGEBRAIC[35] = (fabs(STATES[0])>CONSTANTS[48] ? ( CONSTANTS[35]*STATES[0]*( CONSTANTS[60]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - ALGEBRAIC[9]))/(exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - 1.00000) : ( CONSTANTS[35]*CONSTANTS[2]*(CONSTANTS[60] - ALGEBRAIC[9]))/CONSTANTS[3]);
RATES[8] = ( (- 1.00000/2.30300)*CONSTANTS[42]*((((((ALGEBRAIC[20] -  5.00000*ALGEBRAIC[22]) -  2.00000*ALGEBRAIC[23]) -  4.00000*ALGEBRAIC[24])+ (CONSTANTS[14] - 1.00000)*ALGEBRAIC[25]+ 2.00000*ALGEBRAIC[33]) - ALGEBRAIC[39])+ALGEBRAIC[35]))/CONSTANTS[57];
ALGEBRAIC[18] = CONSTANTS[16] - STATES[26];
ALGEBRAIC[40] =  CONSTANTS[39]*1.00000*( CONSTANTS[11]*STATES[20]*STATES[26]*CONSTANTS[53] -  STATES[19]*ALGEBRAIC[18]);
RATES[19] = ((( - CONSTANTS[61]*ALGEBRAIC[34]+ALGEBRAIC[38]) - CONSTANTS[63]/CONSTANTS[62])+ALGEBRAIC[40])/CONSTANTS[59];
RATES[20] = ((( - CONSTANTS[61]*ALGEBRAIC[32] -  2.00000*ALGEBRAIC[38])+CONSTANTS[63]/CONSTANTS[62]) - ALGEBRAIC[40])/CONSTANTS[59];
RATES[26] = - ALGEBRAIC[40]/CONSTANTS[59];
ALGEBRAIC[41] = (fabs(STATES[0])>CONSTANTS[49] ? ( CONSTANTS[36]*STATES[0]*( CONSTANTS[55]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - STATES[9]))/(exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - 1.00000) : ( CONSTANTS[36]*CONSTANTS[2]*(CONSTANTS[55] - STATES[9]))/CONSTANTS[3]);
RATES[9] = (ALGEBRAIC[39]+ALGEBRAIC[41])/CONSTANTS[57];
RATES[0] = ((((( 4.00000*ALGEBRAIC[22]+ 2.00000*ALGEBRAIC[23]+ 4.00000*ALGEBRAIC[24]) -  CONSTANTS[14]*ALGEBRAIC[25]) - ALGEBRAIC[28]) - ALGEBRAIC[35]) - ALGEBRAIC[41])/CONSTANTS[43];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = STATES[19] - STATES[22];
ALGEBRAIC[5] =  CONSTANTS[40]*1.00000*( ALGEBRAIC[1]*STATES[25] -  CONSTANTS[8]*STATES[22]);
ALGEBRAIC[3] = STATES[20] - STATES[23];
ALGEBRAIC[7] =  CONSTANTS[40]*1.00000*( ALGEBRAIC[3]*STATES[25] -  CONSTANTS[9]*STATES[23]);
ALGEBRAIC[6] = STATES[4] - STATES[5];
ALGEBRAIC[10] =  CONSTANTS[40]*( ALGEBRAIC[6]*STATES[10] -  CONSTANTS[8]*STATES[5]);
ALGEBRAIC[8] = STATES[6] - STATES[7];
ALGEBRAIC[12] =  CONSTANTS[40]*( ALGEBRAIC[8]*STATES[10] -  CONSTANTS[9]*STATES[7]);
ALGEBRAIC[16] = STATES[12] - STATES[13];
ALGEBRAIC[19] =  CONSTANTS[40]*( ALGEBRAIC[16]*STATES[18] -  CONSTANTS[8]*STATES[13]);
ALGEBRAIC[17] = STATES[14] - STATES[15];
ALGEBRAIC[21] =  CONSTANTS[40]*( ALGEBRAIC[17]*STATES[18] -  CONSTANTS[9]*STATES[15]);
ALGEBRAIC[0] = CONSTANTS[18] - STATES[1];
ALGEBRAIC[20] = ( CONSTANTS[27]*( CONSTANTS[26]*ALGEBRAIC[0] - STATES[1])*(1.00000+STATES[11]/CONSTANTS[21]))/(1.00000+STATES[11]/CONSTANTS[22]);
ALGEBRAIC[9] =  1.00000*pow(10.0000, - STATES[8]);
ALGEBRAIC[13] = CONSTANTS[4] -  1.00000*CONSTANTS[2]*log(ALGEBRAIC[9]/1.00000e-07);
ALGEBRAIC[11] =  CONSTANTS[3]*STATES[0]+ CONSTANTS[2]*log(CONSTANTS[60]/ALGEBRAIC[9]);
ALGEBRAIC[2] = CONSTANTS[17] - STATES[2];
ALGEBRAIC[22] =  CONSTANTS[28]*( exp(- (ALGEBRAIC[13]+ 4.00000*ALGEBRAIC[11])/CONSTANTS[2])*STATES[1]*STATES[2] -  ALGEBRAIC[0]*ALGEBRAIC[2]);
ALGEBRAIC[14] = CONSTANTS[5]+ 2.00000*CONSTANTS[2]*log(ALGEBRAIC[9]/1.00000e-07);
ALGEBRAIC[4] = CONSTANTS[15] - STATES[3];
ALGEBRAIC[23] =  (( CONSTANTS[29]*(1.00000+STATES[11]/CONSTANTS[23]))/(1.00000+STATES[11]/CONSTANTS[24]))*( exp(- ((ALGEBRAIC[14]+ 4.00000*ALGEBRAIC[11]) -  2.00000*CONSTANTS[3]*STATES[0])/( 2.00000*CONSTANTS[2]))*STATES[3]* pow(ALGEBRAIC[2], 1.0 / 2) -  ALGEBRAIC[4]* pow(STATES[2], 1.0 / 2));
ALGEBRAIC[15] = CONSTANTS[6] -  2.00000*CONSTANTS[2]*log(ALGEBRAIC[9]/1.00000e-07);
ALGEBRAIC[24] = (CONSTANTS[0]<1.00000e-12 ? 0.00000 :  (( CONSTANTS[30]*CONSTANTS[0])/(CONSTANTS[0]+CONSTANTS[41]))*(ALGEBRAIC[4]/CONSTANTS[15])*( exp(- (ALGEBRAIC[15]+ 2.00000*ALGEBRAIC[11])/( 2.00000*CONSTANTS[2]))*ALGEBRAIC[4]*pow(CONSTANTS[0]/1.00000, 0.250000) -  STATES[3]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2])));
ALGEBRAIC[25] =  CONSTANTS[31]*( (( exp(- (CONSTANTS[7] -  CONSTANTS[14]*ALGEBRAIC[11])/CONSTANTS[2])*CONSTANTS[9])/CONSTANTS[8])*STATES[7]*STATES[11] -  STATES[5]*1.00000);
ALGEBRAIC[26] =  - 0.650000*STATES[0];
ALGEBRAIC[27] =  0.350000*STATES[0];
ALGEBRAIC[28] = (ALGEBRAIC[17]>CONSTANTS[47]||ALGEBRAIC[16]>CONSTANTS[47] ?  CONSTANTS[32]*(ALGEBRAIC[17]/(ALGEBRAIC[17]+ ALGEBRAIC[16]*exp(( - CONSTANTS[3]*ALGEBRAIC[27])/CONSTANTS[2])) - ALGEBRAIC[8]/(ALGEBRAIC[8]+ ALGEBRAIC[6]*exp(( - CONSTANTS[3]*ALGEBRAIC[26])/CONSTANTS[2])))*(ALGEBRAIC[17]/(ALGEBRAIC[17]+CONSTANTS[37])) : 0.00000);
ALGEBRAIC[29] = ( STATES[17]*CONSTANTS[60])/(CONSTANTS[60]+CONSTANTS[50]);
ALGEBRAIC[31] = ( STATES[11]*ALGEBRAIC[9])/(ALGEBRAIC[9]+CONSTANTS[50]);
ALGEBRAIC[33] = ( CONSTANTS[33]*( CONSTANTS[60]*ALGEBRAIC[29] -  ALGEBRAIC[9]*ALGEBRAIC[31]))/(ALGEBRAIC[29]+CONSTANTS[25]);
ALGEBRAIC[34] =  CONSTANTS[46]*CONSTANTS[44]*(STATES[19] - STATES[12]);
ALGEBRAIC[36] =  CONSTANTS[38]*( CONSTANTS[10]*STATES[14]*STATES[14] -  STATES[16]*STATES[12]);
ALGEBRAIC[32] =  CONSTANTS[46]*CONSTANTS[44]*(STATES[20] - STATES[14]);
ALGEBRAIC[30] =  CONSTANTS[46]*CONSTANTS[44]*(STATES[21] - STATES[16]);
ALGEBRAIC[37] =  CONSTANTS[46]*CONSTANTS[45]*(STATES[24] - STATES[17]);
ALGEBRAIC[38] =  CONSTANTS[38]*1.00000*( CONSTANTS[10]*STATES[20]*STATES[20] -  STATES[21]*STATES[19]);
ALGEBRAIC[39] =  CONSTANTS[34]*( CONSTANTS[55]*ALGEBRAIC[9] -  STATES[9]*CONSTANTS[60]);
ALGEBRAIC[35] = (fabs(STATES[0])>CONSTANTS[48] ? ( CONSTANTS[35]*STATES[0]*( CONSTANTS[60]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - ALGEBRAIC[9]))/(exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - 1.00000) : ( CONSTANTS[35]*CONSTANTS[2]*(CONSTANTS[60] - ALGEBRAIC[9]))/CONSTANTS[3]);
ALGEBRAIC[18] = CONSTANTS[16] - STATES[26];
ALGEBRAIC[40] =  CONSTANTS[39]*1.00000*( CONSTANTS[11]*STATES[20]*STATES[26]*CONSTANTS[53] -  STATES[19]*ALGEBRAIC[18]);
ALGEBRAIC[41] = (fabs(STATES[0])>CONSTANTS[49] ? ( CONSTANTS[36]*STATES[0]*( CONSTANTS[55]*exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - STATES[9]))/(exp(( CONSTANTS[3]*STATES[0])/CONSTANTS[2]) - 1.00000) : ( CONSTANTS[36]*CONSTANTS[2]*(CONSTANTS[55] - STATES[9]))/CONSTANTS[3]);
}