/*
   There are a total of 41 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 65 entries in the constant variable array.
 */
/*
 * VOI is time in component Environment (second).
 * CONSTANTS[0] is RT in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[1] is F in component fixed_parameters (kilojoule_per_mole_per_millivolt).
 * CONSTANTS[2] is dG_C1o in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[3] is dG_C3o in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[4] is dG_C4o in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[3] is dG_C1op in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[6] is dG_C3op in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[7] is dG_C4op in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[5] is dG_F1o in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[0] is dG_H in component fixed_parameters (kilojoule_per_mole).
 * CONSTANTS[6] is k_dHPI in component fixed_parameters (molar).
 * CONSTANTS[7] is k_dHatp in component fixed_parameters (molar).
 * CONSTANTS[8] is k_dHadp in component fixed_parameters (molar).
 * CONSTANTS[9] is K_DT in component fixed_parameters (molar).
 * CONSTANTS[10] is K_DD in component fixed_parameters (molar).
 * CONSTANTS[11] is K_AK in component fixed_parameters (dimensionless).
 * CONSTANTS[12] is K_CK in component fixed_parameters (per_molar).
 * CONSTANTS[13] is W_c in component fixed_parameters (l_water_per_l_cyto).
 * CONSTANTS[14] is W_x in component fixed_parameters (l_water_per_l_mito).
 * CONSTANTS[15] is W_i in component fixed_parameters (l_water_per_l_mito).
 * CONSTANTS[16] is V_cyto in component fixed_parameters (cyto_per_cell).
 * CONSTANTS[17] is V_mito in component fixed_parameters (mito_per_cell).
 * CONSTANTS[61] is Rm_cyto in component fixed_parameters (mito_per_cyto).
 * CONSTANTS[58] is Rm_cell in component fixed_parameters (mito_per_cell).
 * CONSTANTS[63] is Rc_cell in component fixed_parameters (cyto_per_cell).
 * CONSTANTS[18] is n_A in component fixed_parameters (dimensionless).
 * CONSTANTS[19] is C_tot in component fixed_parameters (molar).
 * CONSTANTS[20] is CR_tot in component fixed_parameters (molar).
 * CONSTANTS[21] is Q_tot in component fixed_parameters (molar).
 * CONSTANTS[22] is NAD_tot in component fixed_parameters (molar).
 * CONSTANTS[23] is pH_C in component fixed_parameters (dimensionless).
 * CONSTANTS[59] is H_c in component fixed_parameters (molar).
 * CONSTANTS[24] is K_c in component fixed_parameters (molar).
 * CONSTANTS[25] is ATP_c in component fixed_parameters (molar).
 * CONSTANTS[26] is AMP_c in component fixed_parameters (molar).
 * CONSTANTS[62] is H_i in component fixed_parameters (molar).
 * CONSTANTS[60] is K_i in component fixed_parameters (molar).
 * STATES[0] is H_x in component dH_x_dt (molar).
 * STATES[1] is deltaPsi in component ddeltaPsi_dt (millivolt).
 * CONSTANTS[27] is k_PI1 in component adjustable_parameters (molar).
 * CONSTANTS[28] is k_PI2 in component adjustable_parameters (molar).
 * CONSTANTS[29] is k_PI3 in component adjustable_parameters (molar).
 * CONSTANTS[30] is k_PI4 in component adjustable_parameters (molar).
 * CONSTANTS[31] is k_PIH in component adjustable_parameters (molar).
 * CONSTANTS[32] is r in component adjustable_parameters (dimensionless).
 * CONSTANTS[33] is X_DH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
 * CONSTANTS[34] is X_C1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[35] is X_C3 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_half_molar).
 * CONSTANTS[36] is X_C4 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
 * CONSTANTS[37] is X_F1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[38] is X_ANT in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
 * CONSTANTS[39] is X_PI1 in component adjustable_parameters (mole_per_second_per_l_mito_per_molar).
 * CONSTANTS[40] is X_KH in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[41] is X_Hle in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
 * CONSTANTS[42] is X_K in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
 * CONSTANTS[43] is k_mADP in component adjustable_parameters (molar).
 * CONSTANTS[44] is X_AK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[45] is X_CK in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[46] is X_MgA in component adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
 * CONSTANTS[47] is k_O2 in component adjustable_parameters (molar).
 * CONSTANTS[48] is beta in component adjustable_parameters (molar).
 * CONSTANTS[49] is C_IM in component adjustable_parameters (mole_per_l_mito_per_millivolt).
 * CONSTANTS[50] is X_A in component adjustable_parameters (micron_per_second).
 * CONSTANTS[51] is X_PI2 in component adjustable_parameters (micron_per_second).
 * CONSTANTS[52] is gamma in component adjustable_parameters (per_micron).
 * ALGEBRAIC[8] is NAD_x in component dNAD_x_dt (molar).
 * STATES[2] is NADH_x in component dNADH_x_dt (molar).
 * ALGEBRAIC[9] is QH2_x in component dQH2_x_dt (molar).
 * STATES[3] is Q_x in component dQ_x_dt (molar).
 * ALGEBRAIC[10] is Cred_x in component dCred_x_dt (molar).
 * STATES[4] is Cox_x in component dCox_x_dt (molar).
 * ALGEBRAIC[11] is ATP_fx in component dATP_fx_dt (molar).
 * STATES[5] is ATP_x in component dATP_x_dt (molar).
 * STATES[6] is ATP_mx in component dATP_mx_dt (molar).
 * ALGEBRAIC[12] is ADP_fx in component dADP_fx_dt (molar).
 * STATES[7] is ADP_x in component dADP_x_dt (molar).
 * STATES[8] is ADP_mx in component dADP_mx_dt (molar).
 * ALGEBRAIC[19] is J_DH in component J_DH (mole_per_second_per_l_mito).
 * ALGEBRAIC[21] is J_C1 in component J_C1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[22] is J_C3 in component J_C3 (mole_per_second_per_l_mito).
 * ALGEBRAIC[23] is J_C4 in component J_C4 (mole_per_second_per_l_mito).
 * ALGEBRAIC[24] is J_F1 in component J_F1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[32] is J_PI1 in component J_PI1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[38] is J_KH in component J_KH (mole_per_second_per_l_mito).
 * ALGEBRAIC[34] 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[40] 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[14] is J_MgATP_x in component J_MgATP_x (mole_per_second_per_l_mito).
 * ALGEBRAIC[16] is J_MgADP_x in component J_MgADP_x (mole_per_second_per_l_mito).
 * ALGEBRAIC[27] 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[13] 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[15] 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[33] is J_ATP in component J_ATP (mole_per_second_per_l_mito).
 * ALGEBRAIC[35] is J_AKi in component J_AKi (mole_per_second_per_l_mito).
 * ALGEBRAIC[31] 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[29] is J_AMP in component J_AMP (mole_per_second_per_l_mito).
 * ALGEBRAIC[18] is J_MgATP_i in component J_MgATP_i (mole_per_second_per_l_mito).
 * ALGEBRAIC[20] 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[36] 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[37] is J_AKc in component J_AKc (mole_per_second_per_l_cyto).
 * CONSTANTS[64] is J_AtC in component J_AtC (mole_per_second_per_l_cell).
 * ALGEBRAIC[39] 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[4] 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[5] 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[2] is ADP_fc in component dADP_fc_dt (molar).
 * ALGEBRAIC[17] is Cr_c in component dCr_c_dt (molar).
 * CONSTANTS[53] is O2_x in component J_C4 (molar).
 * CONSTANTS[54] is mincon in component J_ANT (molar).
 * ALGEBRAIC[25] is Psi_x in component J_ANT (millivolt).
 * ALGEBRAIC[26] is Psi_i in component J_ANT (millivolt).
 * ALGEBRAIC[28] is H2PIi in component J_PI1 (molar).
 * ALGEBRAIC[30] is H2PIx in component J_PI1 (molar).
 * CONSTANTS[55] is mincond in component J_Hle (millivolt).
 * CONSTANTS[56] is mincond in component J_K (millivolt).
 * CONSTANTS[57] is X_AtC in component J_AtC (mole_per_second_per_l_cell).
 * RATES[0] is d/dt H_x in component dH_x_dt (molar).
 * 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[2] is d/dt NADH_x in component dNADH_x_dt (molar).
 * RATES[3] is d/dt Q_x in component dQ_x_dt (molar).
 * RATES[4] is d/dt Cox_x in component dCox_x_dt (molar).
 * RATES[5] is d/dt ATP_x in component dATP_x_dt (molar).
 * RATES[7] is d/dt ADP_x in component dADP_x_dt (molar).
 * RATES[6] is d/dt ATP_mx in component dATP_mx_dt (molar).
 * RATES[8] 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[1] is d/dt deltaPsi in component ddeltaPsi_dt (millivolt).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 2.5775;
CONSTANTS[1] = 0.096484;
CONSTANTS[2] = -69.37;
CONSTANTS[3] = -32.53;
CONSTANTS[4] = -122.94;
CONSTANTS[5] = 36.03;
CONSTANTS[6] = 1.7782754100389227e-7;
CONSTANTS[7] = 3.311311214825911e-7;
CONSTANTS[8] = 5.128613839913649e-7;
CONSTANTS[9] = 24e-6;
CONSTANTS[10] = 347e-6;
CONSTANTS[11] = 0.4331;
CONSTANTS[12] = 1.66e9;
CONSTANTS[13] = 0.8425;
CONSTANTS[14] = 0.6514;
CONSTANTS[15] = 0.0724;
CONSTANTS[16] = 0.894;
CONSTANTS[17] = 0.056;
CONSTANTS[18] = 3;
CONSTANTS[19] = 2.7e-3;
CONSTANTS[20] = 42.7e-3;
CONSTANTS[21] = 1.35e-3;
CONSTANTS[22] = 2.97e-3;
CONSTANTS[23] = 7;
CONSTANTS[24] = 0.15;
CONSTANTS[25] = 0;
CONSTANTS[26] = 0;
STATES[0] = 1e-7;
STATES[1] = 166.67;
CONSTANTS[27] = 0.1553e-3;
CONSTANTS[28] = 0.8222e-3;
CONSTANTS[29] = 0.3601e-3;
CONSTANTS[30] = 5.924e-3;
CONSTANTS[31] = 2.542e-4;
CONSTANTS[32] = 4.599;
CONSTANTS[33] = 0.0866;
CONSTANTS[34] = 4.405e3;
CONSTANTS[35] = 4.887;
CONSTANTS[36] = 6.766e-5;
CONSTANTS[37] = 1e3;
CONSTANTS[38] = 8.123e-3;
CONSTANTS[39] = 3.85e5;
CONSTANTS[40] = 5.651e7;
CONSTANTS[41] = 200;
CONSTANTS[42] = 0;
CONSTANTS[43] = 3.5e-6;
CONSTANTS[44] = 1e7;
CONSTANTS[45] = 1e7;
CONSTANTS[46] = 1e7;
CONSTANTS[47] = 1.2e-4;
CONSTANTS[48] = 0.01;
CONSTANTS[49] = 6.75e-6;
CONSTANTS[50] = 85;
CONSTANTS[51] = 327;
CONSTANTS[52] = 5.99;
STATES[2] = 0.0015723;
STATES[3] = 6.75e-4;
STATES[4] = 1.35e-3;
STATES[5] = 0.0026657;
STATES[6] = 0.0026046;
STATES[7] = 0.0073343;
STATES[8] = 0.0054765;
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[53] = 0;
CONSTANTS[54] = 1e-9;
CONSTANTS[55] = 1e-6;
CONSTANTS[56] = 1e-6;
CONSTANTS[57] = 0;
CONSTANTS[58] = CONSTANTS[17];
CONSTANTS[59] =  1.00000*pow(10.0000, - CONSTANTS[23]);
CONSTANTS[60] = CONSTANTS[24];
CONSTANTS[61] = CONSTANTS[17]/CONSTANTS[16];
CONSTANTS[62] = CONSTANTS[59];
CONSTANTS[63] = CONSTANTS[16];
CONSTANTS[64] = CONSTANTS[57];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = STATES[19] - STATES[22];
ALGEBRAIC[4] =  CONSTANTS[46]*( ALGEBRAIC[1]*STATES[25] -  CONSTANTS[9]*STATES[22]);
RATES[22] = ALGEBRAIC[4]/CONSTANTS[13];
ALGEBRAIC[2] = STATES[20] - STATES[23];
ALGEBRAIC[5] =  CONSTANTS[46]*( ALGEBRAIC[2]*STATES[25] -  CONSTANTS[10]*STATES[23]);
RATES[23] = ALGEBRAIC[5]/CONSTANTS[13];
RATES[25] = (- ALGEBRAIC[4] - ALGEBRAIC[4])/CONSTANTS[13];
ALGEBRAIC[11] = STATES[5] - STATES[6];
ALGEBRAIC[14] =  CONSTANTS[46]*( ALGEBRAIC[11]*STATES[10] -  CONSTANTS[9]*STATES[6]);
RATES[6] = ALGEBRAIC[14]/CONSTANTS[14];
ALGEBRAIC[12] = STATES[7] - STATES[8];
ALGEBRAIC[16] =  CONSTANTS[46]*( ALGEBRAIC[12]*STATES[10] -  CONSTANTS[10]*STATES[8]);
RATES[10] = (- ALGEBRAIC[16] - ALGEBRAIC[14])/CONSTANTS[14];
RATES[8] = ALGEBRAIC[16]/CONSTANTS[14];
ALGEBRAIC[13] = STATES[12] - STATES[13];
ALGEBRAIC[18] =  CONSTANTS[46]*( ALGEBRAIC[13]*STATES[18] -  CONSTANTS[9]*STATES[13]);
RATES[13] = ALGEBRAIC[18]/CONSTANTS[15];
ALGEBRAIC[15] = STATES[14] - STATES[15];
ALGEBRAIC[20] =  CONSTANTS[46]*( ALGEBRAIC[15]*STATES[18] -  CONSTANTS[10]*STATES[15]);
RATES[15] = ALGEBRAIC[20]/CONSTANTS[15];
RATES[18] = (- ALGEBRAIC[20] - ALGEBRAIC[18])/CONSTANTS[15];
ALGEBRAIC[8] = CONSTANTS[22] - STATES[2];
ALGEBRAIC[19] = ( CONSTANTS[33]*( CONSTANTS[32]*ALGEBRAIC[8] - STATES[2])*(1.00000+STATES[11]/CONSTANTS[27]))/(1.00000+STATES[11]/CONSTANTS[28]);
ALGEBRAIC[3] = CONSTANTS[2] -  1.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[0] =  CONSTANTS[1]*STATES[1]+ CONSTANTS[0]*arbitrary_log(CONSTANTS[62]/STATES[0], 10);
ALGEBRAIC[9] = CONSTANTS[21] - STATES[3];
ALGEBRAIC[21] =  CONSTANTS[34]*( exp(- (ALGEBRAIC[3]+ 4.00000*ALGEBRAIC[0])/CONSTANTS[0])*STATES[2]*STATES[3] -  ALGEBRAIC[8]*ALGEBRAIC[9]);
RATES[2] = (ALGEBRAIC[19] - ALGEBRAIC[21])/CONSTANTS[14];
ALGEBRAIC[6] = CONSTANTS[3]+ 2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[10] = CONSTANTS[19] - STATES[4];
ALGEBRAIC[22] =  (( CONSTANTS[35]*(1.00000+STATES[11]/CONSTANTS[29]))/(1.00000+STATES[11]/CONSTANTS[30]))*( exp(- ((ALGEBRAIC[6]+ 4.00000*ALGEBRAIC[0]) -  2.00000*CONSTANTS[1]*STATES[1])/( 2.00000*CONSTANTS[0]))*STATES[4]* pow(ALGEBRAIC[9], 1.0 / 2) -  ALGEBRAIC[10]* pow(STATES[3], 1.0 / 2));
RATES[3] = (ALGEBRAIC[21] - ALGEBRAIC[22])/CONSTANTS[14];
ALGEBRAIC[7] = CONSTANTS[4] -  2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[23] = (CONSTANTS[53]<1.00000e-12 ? 0.00000 :  (( (( CONSTANTS[36]*CONSTANTS[53])/(CONSTANTS[53]+CONSTANTS[47]))*ALGEBRAIC[10])/CONSTANTS[19])*( exp(- (ALGEBRAIC[7]+ 2.00000*ALGEBRAIC[0])/( 2.00000*CONSTANTS[0]))*ALGEBRAIC[10]*pow(CONSTANTS[53], 0.250000) -  STATES[4]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0])));
RATES[4] = ( 2.00000*ALGEBRAIC[22] -  2.00000*ALGEBRAIC[23])/CONSTANTS[14];
ALGEBRAIC[24] =  CONSTANTS[37]*( (( exp(- (CONSTANTS[5] -  CONSTANTS[18]*ALGEBRAIC[0])/CONSTANTS[0])*CONSTANTS[10])/CONSTANTS[9])*STATES[8]*STATES[11] -  STATES[6]*1.00000);
ALGEBRAIC[25] =  - 0.650000*STATES[1];
ALGEBRAIC[26] =  0.350000*STATES[1];
ALGEBRAIC[27] = (ALGEBRAIC[15]>CONSTANTS[54]&&ALGEBRAIC[13]>CONSTANTS[54] ? ( CONSTANTS[38]*(ALGEBRAIC[15]/(ALGEBRAIC[15]+ ALGEBRAIC[13]*exp(( - CONSTANTS[1]*ALGEBRAIC[26])/CONSTANTS[0])) - ALGEBRAIC[12]/(ALGEBRAIC[12]+ ALGEBRAIC[11]*exp(( - CONSTANTS[1]*ALGEBRAIC[25])/CONSTANTS[0])))*ALGEBRAIC[15])/(ALGEBRAIC[15]+CONSTANTS[43]) : 0.00000);
RATES[5] = (ALGEBRAIC[24] - ALGEBRAIC[27])/CONSTANTS[14];
RATES[7] = (- ALGEBRAIC[24]+ALGEBRAIC[27])/CONSTANTS[14];
ALGEBRAIC[28] = ( STATES[17]*CONSTANTS[62])/(CONSTANTS[62]+CONSTANTS[6]);
ALGEBRAIC[30] = ( STATES[11]*STATES[0])/(STATES[0]+CONSTANTS[6]);
ALGEBRAIC[32] = ( CONSTANTS[39]*( CONSTANTS[62]*ALGEBRAIC[28] -  STATES[0]*ALGEBRAIC[30]))/(ALGEBRAIC[28]+CONSTANTS[31]);
RATES[11] = (- ALGEBRAIC[24] - ALGEBRAIC[32])/CONSTANTS[14];
ALGEBRAIC[33] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[19] - STATES[12]);
ALGEBRAIC[35] =  CONSTANTS[44]*( CONSTANTS[11]*STATES[14]*STATES[14] -  STATES[16]*STATES[12]);
RATES[12] = (ALGEBRAIC[33]+ALGEBRAIC[27]+ALGEBRAIC[35])/CONSTANTS[15];
ALGEBRAIC[31] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[20] - STATES[14]);
RATES[14] = ((ALGEBRAIC[31] - ALGEBRAIC[27]) -  2.00000*ALGEBRAIC[35])/CONSTANTS[15];
ALGEBRAIC[29] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[21] - STATES[16]);
RATES[16] = (ALGEBRAIC[29]+ALGEBRAIC[35])/CONSTANTS[15];
ALGEBRAIC[36] =  CONSTANTS[52]*CONSTANTS[51]*(STATES[24] - STATES[17]);
RATES[17] = (- ALGEBRAIC[32]+ALGEBRAIC[36])/CONSTANTS[15];
ALGEBRAIC[37] =  CONSTANTS[44]*( CONSTANTS[11]*STATES[20]*STATES[20] -  STATES[21]*STATES[19]);
RATES[21] = ( - CONSTANTS[61]*ALGEBRAIC[29]+ALGEBRAIC[37])/CONSTANTS[13];
RATES[24] = ( - ALGEBRAIC[36]*CONSTANTS[61]+CONSTANTS[64]/CONSTANTS[63])/CONSTANTS[13];
ALGEBRAIC[38] =  CONSTANTS[40]*( CONSTANTS[60]*STATES[0] -  STATES[9]*CONSTANTS[62]);
ALGEBRAIC[34] = (fabs(STATES[1])>CONSTANTS[55] ? ( CONSTANTS[41]*STATES[1]*( CONSTANTS[62]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[0]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[41]*CONSTANTS[0]*(CONSTANTS[62] - STATES[0]))/CONSTANTS[1]);
RATES[0] = ( (STATES[0]/CONSTANTS[48])*((((((ALGEBRAIC[19] -  5.00000*ALGEBRAIC[21]) -  2.00000*ALGEBRAIC[22]) -  4.00000*ALGEBRAIC[23])+ (CONSTANTS[18] - 1.00000)*ALGEBRAIC[24]+ 2.00000*ALGEBRAIC[32]) - ALGEBRAIC[38])+ALGEBRAIC[34]))/CONSTANTS[14];
ALGEBRAIC[17] = CONSTANTS[20] - STATES[26];
ALGEBRAIC[39] =  CONSTANTS[45]*( CONSTANTS[12]*STATES[20]*STATES[26]*CONSTANTS[59] -  STATES[19]*ALGEBRAIC[17]);
RATES[19] = ((( - CONSTANTS[61]*ALGEBRAIC[33]+ALGEBRAIC[37]) - CONSTANTS[64]/CONSTANTS[63])+ALGEBRAIC[39])/CONSTANTS[13];
RATES[20] = ((( - CONSTANTS[61]*ALGEBRAIC[31] -  2.00000*ALGEBRAIC[37])+CONSTANTS[64]/CONSTANTS[63]) - ALGEBRAIC[39])/CONSTANTS[13];
RATES[26] = - ALGEBRAIC[39]/CONSTANTS[13];
ALGEBRAIC[40] = (fabs(STATES[1])>CONSTANTS[56] ? ( CONSTANTS[42]*STATES[1]*( CONSTANTS[60]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[9]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[42]*CONSTANTS[0]*(CONSTANTS[60] - STATES[9]))/CONSTANTS[1]);
RATES[9] = (ALGEBRAIC[38]+ALGEBRAIC[40])/CONSTANTS[14];
RATES[1] = ((((( 4.00000*ALGEBRAIC[21]+ 2.00000*ALGEBRAIC[22]+ 4.00000*ALGEBRAIC[23]) -  CONSTANTS[18]*ALGEBRAIC[24]) - ALGEBRAIC[27]) - ALGEBRAIC[34]) - ALGEBRAIC[40])/CONSTANTS[49];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = STATES[19] - STATES[22];
ALGEBRAIC[4] =  CONSTANTS[46]*( ALGEBRAIC[1]*STATES[25] -  CONSTANTS[9]*STATES[22]);
ALGEBRAIC[2] = STATES[20] - STATES[23];
ALGEBRAIC[5] =  CONSTANTS[46]*( ALGEBRAIC[2]*STATES[25] -  CONSTANTS[10]*STATES[23]);
ALGEBRAIC[11] = STATES[5] - STATES[6];
ALGEBRAIC[14] =  CONSTANTS[46]*( ALGEBRAIC[11]*STATES[10] -  CONSTANTS[9]*STATES[6]);
ALGEBRAIC[12] = STATES[7] - STATES[8];
ALGEBRAIC[16] =  CONSTANTS[46]*( ALGEBRAIC[12]*STATES[10] -  CONSTANTS[10]*STATES[8]);
ALGEBRAIC[13] = STATES[12] - STATES[13];
ALGEBRAIC[18] =  CONSTANTS[46]*( ALGEBRAIC[13]*STATES[18] -  CONSTANTS[9]*STATES[13]);
ALGEBRAIC[15] = STATES[14] - STATES[15];
ALGEBRAIC[20] =  CONSTANTS[46]*( ALGEBRAIC[15]*STATES[18] -  CONSTANTS[10]*STATES[15]);
ALGEBRAIC[8] = CONSTANTS[22] - STATES[2];
ALGEBRAIC[19] = ( CONSTANTS[33]*( CONSTANTS[32]*ALGEBRAIC[8] - STATES[2])*(1.00000+STATES[11]/CONSTANTS[27]))/(1.00000+STATES[11]/CONSTANTS[28]);
ALGEBRAIC[3] = CONSTANTS[2] -  1.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[0] =  CONSTANTS[1]*STATES[1]+ CONSTANTS[0]*arbitrary_log(CONSTANTS[62]/STATES[0], 10);
ALGEBRAIC[9] = CONSTANTS[21] - STATES[3];
ALGEBRAIC[21] =  CONSTANTS[34]*( exp(- (ALGEBRAIC[3]+ 4.00000*ALGEBRAIC[0])/CONSTANTS[0])*STATES[2]*STATES[3] -  ALGEBRAIC[8]*ALGEBRAIC[9]);
ALGEBRAIC[6] = CONSTANTS[3]+ 2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[10] = CONSTANTS[19] - STATES[4];
ALGEBRAIC[22] =  (( CONSTANTS[35]*(1.00000+STATES[11]/CONSTANTS[29]))/(1.00000+STATES[11]/CONSTANTS[30]))*( exp(- ((ALGEBRAIC[6]+ 4.00000*ALGEBRAIC[0]) -  2.00000*CONSTANTS[1]*STATES[1])/( 2.00000*CONSTANTS[0]))*STATES[4]* pow(ALGEBRAIC[9], 1.0 / 2) -  ALGEBRAIC[10]* pow(STATES[3], 1.0 / 2));
ALGEBRAIC[7] = CONSTANTS[4] -  2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[23] = (CONSTANTS[53]<1.00000e-12 ? 0.00000 :  (( (( CONSTANTS[36]*CONSTANTS[53])/(CONSTANTS[53]+CONSTANTS[47]))*ALGEBRAIC[10])/CONSTANTS[19])*( exp(- (ALGEBRAIC[7]+ 2.00000*ALGEBRAIC[0])/( 2.00000*CONSTANTS[0]))*ALGEBRAIC[10]*pow(CONSTANTS[53], 0.250000) -  STATES[4]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0])));
ALGEBRAIC[24] =  CONSTANTS[37]*( (( exp(- (CONSTANTS[5] -  CONSTANTS[18]*ALGEBRAIC[0])/CONSTANTS[0])*CONSTANTS[10])/CONSTANTS[9])*STATES[8]*STATES[11] -  STATES[6]*1.00000);
ALGEBRAIC[25] =  - 0.650000*STATES[1];
ALGEBRAIC[26] =  0.350000*STATES[1];
ALGEBRAIC[27] = (ALGEBRAIC[15]>CONSTANTS[54]&&ALGEBRAIC[13]>CONSTANTS[54] ? ( CONSTANTS[38]*(ALGEBRAIC[15]/(ALGEBRAIC[15]+ ALGEBRAIC[13]*exp(( - CONSTANTS[1]*ALGEBRAIC[26])/CONSTANTS[0])) - ALGEBRAIC[12]/(ALGEBRAIC[12]+ ALGEBRAIC[11]*exp(( - CONSTANTS[1]*ALGEBRAIC[25])/CONSTANTS[0])))*ALGEBRAIC[15])/(ALGEBRAIC[15]+CONSTANTS[43]) : 0.00000);
ALGEBRAIC[28] = ( STATES[17]*CONSTANTS[62])/(CONSTANTS[62]+CONSTANTS[6]);
ALGEBRAIC[30] = ( STATES[11]*STATES[0])/(STATES[0]+CONSTANTS[6]);
ALGEBRAIC[32] = ( CONSTANTS[39]*( CONSTANTS[62]*ALGEBRAIC[28] -  STATES[0]*ALGEBRAIC[30]))/(ALGEBRAIC[28]+CONSTANTS[31]);
ALGEBRAIC[33] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[19] - STATES[12]);
ALGEBRAIC[35] =  CONSTANTS[44]*( CONSTANTS[11]*STATES[14]*STATES[14] -  STATES[16]*STATES[12]);
ALGEBRAIC[31] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[20] - STATES[14]);
ALGEBRAIC[29] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[21] - STATES[16]);
ALGEBRAIC[36] =  CONSTANTS[52]*CONSTANTS[51]*(STATES[24] - STATES[17]);
ALGEBRAIC[37] =  CONSTANTS[44]*( CONSTANTS[11]*STATES[20]*STATES[20] -  STATES[21]*STATES[19]);
ALGEBRAIC[38] =  CONSTANTS[40]*( CONSTANTS[60]*STATES[0] -  STATES[9]*CONSTANTS[62]);
ALGEBRAIC[34] = (fabs(STATES[1])>CONSTANTS[55] ? ( CONSTANTS[41]*STATES[1]*( CONSTANTS[62]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[0]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[41]*CONSTANTS[0]*(CONSTANTS[62] - STATES[0]))/CONSTANTS[1]);
ALGEBRAIC[17] = CONSTANTS[20] - STATES[26];
ALGEBRAIC[39] =  CONSTANTS[45]*( CONSTANTS[12]*STATES[20]*STATES[26]*CONSTANTS[59] -  STATES[19]*ALGEBRAIC[17]);
ALGEBRAIC[40] = (fabs(STATES[1])>CONSTANTS[56] ? ( CONSTANTS[42]*STATES[1]*( CONSTANTS[60]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - STATES[9]))/(exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0]) - 1.00000) : ( CONSTANTS[42]*CONSTANTS[0]*(CONSTANTS[60] - STATES[9]))/CONSTANTS[1]);
}