/*
   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[1] is dG_C1op in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[2] is dG_C3op in component fixed_parameters (kilojoule_per_mole).
 * ALGEBRAIC[3] 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[4] is NAD_x in component dNAD_x_dt (molar).
 * STATES[2] is NADH_x in component dNADH_x_dt (molar).
 * ALGEBRAIC[5] is QH2_x in component dQH2_x_dt (molar).
 * STATES[3] is Q_x in component dQ_x_dt (molar).
 * ALGEBRAIC[6] is Cred_x in component dCred_x_dt (molar).
 * STATES[4] is Cox_x in component dCox_x_dt (molar).
 * ALGEBRAIC[7] 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[8] 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[14] is J_DH in component J_DH (mole_per_second_per_l_mito).
 * ALGEBRAIC[15] is J_C1 in component J_C1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[16] is J_C3 in component J_C3 (mole_per_second_per_l_mito).
 * ALGEBRAIC[17] is J_C4 in component J_C4 (mole_per_second_per_l_mito).
 * ALGEBRAIC[18] is J_F1 in component J_F1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[24] is J_PI1 in component J_PI1 (mole_per_second_per_l_mito).
 * ALGEBRAIC[26] is J_KH in component J_KH (mole_per_second_per_l_mito).
 * ALGEBRAIC[25] 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[27] 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[32] is J_MgATP_x in component J_MgATP_x (mole_per_second_per_l_mito).
 * ALGEBRAIC[33] is J_MgADP_x in component J_MgADP_x (mole_per_second_per_l_mito).
 * ALGEBRAIC[21] 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[9] 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[10] 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[30] is J_ATP in component J_ATP (mole_per_second_per_l_mito).
 * ALGEBRAIC[38] is J_AKi in component J_AKi (mole_per_second_per_l_mito).
 * ALGEBRAIC[29] 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[28] is J_AMP in component J_AMP (mole_per_second_per_l_mito).
 * ALGEBRAIC[34] is J_MgATP_i in component J_MgATP_i (mole_per_second_per_l_mito).
 * ALGEBRAIC[35] 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[31] 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[39] 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[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[36] 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[37] 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[11] is ATP_fc in component dATP_fc_dt (molar).
 * ALGEBRAIC[12] is ADP_fc in component dADP_fc_dt (molar).
 * ALGEBRAIC[13] 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[19] is Psi_x in component J_ANT (millivolt).
 * ALGEBRAIC[20] is Psi_i in component J_ANT (millivolt).
 * ALGEBRAIC[22] is H2PIi in component J_PI1 (molar).
 * ALGEBRAIC[23] 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).
 * There are a total of 5 condition variables.
 */
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];
RATES[0] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[7] = 0.1001;
RATES[6] = 0.1001;
RATES[8] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[14] = 0.1001;
RATES[16] = 0.1001;
RATES[13] = 0.1001;
RATES[15] = 0.1001;
RATES[17] = 0.1001;
RATES[18] = 0.1001;
RATES[19] = 0.1001;
RATES[20] = 0.1001;
RATES[21] = 0.1001;
RATES[22] = 0.1001;
RATES[23] = 0.1001;
RATES[24] = 0.1001;
RATES[25] = 0.1001;
RATES[26] = 0.1001;
RATES[1] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - ( (STATES[0]/CONSTANTS[48])*((((((ALGEBRAIC[14] -  5.00000*ALGEBRAIC[15]) -  2.00000*ALGEBRAIC[16]) -  4.00000*ALGEBRAIC[17])+ (CONSTANTS[18] - 1.00000)*ALGEBRAIC[18]+ 2.00000*ALGEBRAIC[24]) - ALGEBRAIC[26])+ALGEBRAIC[25]))/CONSTANTS[14];
resid[1] = RATES[9] - (ALGEBRAIC[26]+ALGEBRAIC[27])/CONSTANTS[14];
resid[2] = RATES[10] - (- ALGEBRAIC[33] - ALGEBRAIC[32])/CONSTANTS[14];
resid[3] = RATES[2] - (ALGEBRAIC[14] - ALGEBRAIC[15])/CONSTANTS[14];
resid[4] = RATES[3] - (ALGEBRAIC[15] - ALGEBRAIC[16])/CONSTANTS[14];
resid[5] = RATES[4] - ( 2.00000*ALGEBRAIC[16] -  2.00000*ALGEBRAIC[17])/CONSTANTS[14];
resid[6] = RATES[5] - (ALGEBRAIC[18] - ALGEBRAIC[21])/CONSTANTS[14];
resid[7] = RATES[7] - (- ALGEBRAIC[18]+ALGEBRAIC[21])/CONSTANTS[14];
resid[8] = RATES[6] - ALGEBRAIC[32]/CONSTANTS[14];
resid[9] = RATES[8] - ALGEBRAIC[33]/CONSTANTS[14];
resid[10] = RATES[11] - (- ALGEBRAIC[18] - ALGEBRAIC[24])/CONSTANTS[14];
resid[11] = RATES[12] - (ALGEBRAIC[30]+ALGEBRAIC[21]+ALGEBRAIC[38])/CONSTANTS[15];
resid[12] = RATES[14] - ((ALGEBRAIC[29] - ALGEBRAIC[21]) -  2.00000*ALGEBRAIC[38])/CONSTANTS[15];
resid[13] = RATES[16] - (ALGEBRAIC[28]+ALGEBRAIC[38])/CONSTANTS[15];
resid[14] = RATES[13] - ALGEBRAIC[34]/CONSTANTS[15];
resid[15] = RATES[15] - ALGEBRAIC[35]/CONSTANTS[15];
resid[16] = RATES[17] - (- ALGEBRAIC[24]+ALGEBRAIC[31])/CONSTANTS[15];
resid[17] = RATES[18] - (- ALGEBRAIC[35] - ALGEBRAIC[34])/CONSTANTS[15];
resid[18] = RATES[19] - ((( - CONSTANTS[61]*ALGEBRAIC[30]+ALGEBRAIC[39]) - CONSTANTS[64]/CONSTANTS[63])+ALGEBRAIC[40])/CONSTANTS[13];
resid[19] = RATES[20] - ((( - CONSTANTS[61]*ALGEBRAIC[29] -  2.00000*ALGEBRAIC[39])+CONSTANTS[64]/CONSTANTS[63]) - ALGEBRAIC[40])/CONSTANTS[13];
resid[20] = RATES[21] - ( - CONSTANTS[61]*ALGEBRAIC[28]+ALGEBRAIC[39])/CONSTANTS[13];
resid[21] = RATES[22] - ALGEBRAIC[36]/CONSTANTS[13];
resid[22] = RATES[23] - ALGEBRAIC[37]/CONSTANTS[13];
resid[23] = RATES[24] - ( - ALGEBRAIC[31]*CONSTANTS[61]+CONSTANTS[64]/CONSTANTS[63])/CONSTANTS[13];
resid[24] = RATES[25] - (- ALGEBRAIC[36] - ALGEBRAIC[36])/CONSTANTS[13];
resid[25] = RATES[26] - - ALGEBRAIC[40]/CONSTANTS[13];
resid[26] = RATES[1] - ((((( 4.00000*ALGEBRAIC[15]+ 2.00000*ALGEBRAIC[16]+ 4.00000*ALGEBRAIC[17]) -  CONSTANTS[18]*ALGEBRAIC[18]) - ALGEBRAIC[21]) - ALGEBRAIC[25]) - ALGEBRAIC[27])/CONSTANTS[49];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = CONSTANTS[22] - STATES[2];
ALGEBRAIC[14] = ( CONSTANTS[33]*( CONSTANTS[32]*ALGEBRAIC[4] - STATES[2])*(1.00000+STATES[11]/CONSTANTS[27]))/(1.00000+STATES[11]/CONSTANTS[28]);
ALGEBRAIC[1] = 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[5] = CONSTANTS[21] - STATES[3];
ALGEBRAIC[15] =  CONSTANTS[34]*( exp(- (ALGEBRAIC[1]+ 4.00000*ALGEBRAIC[0])/CONSTANTS[0])*STATES[2]*STATES[3] -  ALGEBRAIC[4]*ALGEBRAIC[5]);
ALGEBRAIC[2] = CONSTANTS[3]+ 2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[6] = CONSTANTS[19] - STATES[4];
ALGEBRAIC[16] =  (( CONSTANTS[35]*(1.00000+STATES[11]/CONSTANTS[29]))/(1.00000+STATES[11]/CONSTANTS[30]))*( exp(- ((ALGEBRAIC[2]+ 4.00000*ALGEBRAIC[0]) -  2.00000*CONSTANTS[1]*STATES[1])/( 2.00000*CONSTANTS[0]))*STATES[4]* pow(ALGEBRAIC[5], 1.0 / 2) -  ALGEBRAIC[6]* pow(STATES[3], 1.0 / 2));
ALGEBRAIC[3] = CONSTANTS[4] -  2.00000*CONSTANTS[0]*arbitrary_log( (STATES[0]/1.00000)*pow(10.0000, 7.00000), 10);
ALGEBRAIC[17] = (CONDVAR[0]<0.00000 ? 0.00000 :  (( (( CONSTANTS[36]*CONSTANTS[53])/(CONSTANTS[53]+CONSTANTS[47]))*ALGEBRAIC[6])/CONSTANTS[19])*( exp(- (ALGEBRAIC[3]+ 2.00000*ALGEBRAIC[0])/( 2.00000*CONSTANTS[0]))*ALGEBRAIC[6]*pow(CONSTANTS[53], 0.250000) -  STATES[4]*exp(( CONSTANTS[1]*STATES[1])/CONSTANTS[0])));
ALGEBRAIC[18] =  CONSTANTS[37]*( (( exp(- (CONSTANTS[5] -  CONSTANTS[18]*ALGEBRAIC[0])/CONSTANTS[0])*CONSTANTS[10])/CONSTANTS[9])*STATES[8]*STATES[11] -  STATES[6]*1.00000);
ALGEBRAIC[7] = STATES[5] - STATES[6];
ALGEBRAIC[8] = STATES[7] - STATES[8];
ALGEBRAIC[9] = STATES[12] - STATES[13];
ALGEBRAIC[10] = STATES[14] - STATES[15];
ALGEBRAIC[19] =  - 0.650000*STATES[1];
ALGEBRAIC[20] =  0.350000*STATES[1];
ALGEBRAIC[21] = (CONDVAR[1]>0.00000&&CONDVAR[2]>0.00000 ? ( CONSTANTS[38]*(ALGEBRAIC[10]/(ALGEBRAIC[10]+ ALGEBRAIC[9]*exp(( - CONSTANTS[1]*ALGEBRAIC[20])/CONSTANTS[0])) - ALGEBRAIC[8]/(ALGEBRAIC[8]+ ALGEBRAIC[7]*exp(( - CONSTANTS[1]*ALGEBRAIC[19])/CONSTANTS[0])))*ALGEBRAIC[10])/(ALGEBRAIC[10]+CONSTANTS[43]) : 0.00000);
ALGEBRAIC[22] = ( STATES[17]*CONSTANTS[62])/(CONSTANTS[62]+CONSTANTS[6]);
ALGEBRAIC[23] = ( STATES[11]*STATES[0])/(STATES[0]+CONSTANTS[6]);
ALGEBRAIC[24] = ( CONSTANTS[39]*( CONSTANTS[62]*ALGEBRAIC[22] -  STATES[0]*ALGEBRAIC[23]))/(ALGEBRAIC[22]+CONSTANTS[31]);
ALGEBRAIC[25] = (CONDVAR[3]>0.00000 ? ( 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[26] =  CONSTANTS[40]*( CONSTANTS[60]*STATES[0] -  STATES[9]*CONSTANTS[62]);
ALGEBRAIC[27] = (CONDVAR[4]>0.00000 ? ( 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]);
ALGEBRAIC[28] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[21] - STATES[16]);
ALGEBRAIC[29] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[20] - STATES[14]);
ALGEBRAIC[30] =  CONSTANTS[52]*CONSTANTS[50]*(STATES[19] - STATES[12]);
ALGEBRAIC[31] =  CONSTANTS[52]*CONSTANTS[51]*(STATES[24] - STATES[17]);
ALGEBRAIC[32] =  CONSTANTS[46]*( ALGEBRAIC[7]*STATES[10] -  CONSTANTS[9]*STATES[6]);
ALGEBRAIC[33] =  CONSTANTS[46]*( ALGEBRAIC[8]*STATES[10] -  CONSTANTS[10]*STATES[8]);
ALGEBRAIC[34] =  CONSTANTS[46]*( ALGEBRAIC[9]*STATES[18] -  CONSTANTS[9]*STATES[13]);
ALGEBRAIC[35] =  CONSTANTS[46]*( ALGEBRAIC[10]*STATES[18] -  CONSTANTS[10]*STATES[15]);
ALGEBRAIC[11] = STATES[19] - STATES[22];
ALGEBRAIC[36] =  CONSTANTS[46]*( ALGEBRAIC[11]*STATES[25] -  CONSTANTS[9]*STATES[22]);
ALGEBRAIC[12] = STATES[20] - STATES[23];
ALGEBRAIC[37] =  CONSTANTS[46]*( ALGEBRAIC[12]*STATES[25] -  CONSTANTS[10]*STATES[23]);
ALGEBRAIC[38] =  CONSTANTS[44]*( CONSTANTS[11]*STATES[14]*STATES[14] -  STATES[16]*STATES[12]);
ALGEBRAIC[39] =  CONSTANTS[44]*( CONSTANTS[11]*STATES[20]*STATES[20] -  STATES[21]*STATES[19]);
ALGEBRAIC[13] = CONSTANTS[20] - STATES[26];
ALGEBRAIC[40] =  CONSTANTS[45]*( CONSTANTS[12]*STATES[20]*STATES[26]*CONSTANTS[59] -  STATES[19]*ALGEBRAIC[13]);
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
SI[8] = 1.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[13] = 1.0;
SI[14] = 1.0;
SI[15] = 1.0;
SI[16] = 1.0;
SI[17] = 1.0;
SI[18] = 1.0;
SI[19] = 1.0;
SI[20] = 1.0;
SI[21] = 1.0;
SI[22] = 1.0;
SI[23] = 1.0;
SI[24] = 1.0;
SI[25] = 1.0;
SI[26] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = CONSTANTS[53] - 1.00000e-12;
CONDVAR[1] = ALGEBRAIC[10] - CONSTANTS[54];
CONDVAR[2] = ALGEBRAIC[9] - CONSTANTS[54];
CONDVAR[3] = fabs(STATES[1]) - CONSTANTS[55];
CONDVAR[4] = fabs(STATES[1]) - CONSTANTS[56];
}