/*
   There are a total of 33 entries in the algebraic variable array.
   There are a total of 10 entries in each of the rate and state variable arrays.
   There are a total of 61 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * CONSTANTS[0] is Vc_Vg in component volume_ratio (dimensionless).
 * CONSTANTS[1] is Glc_o in component Glc_o (millimolar).
 * STATES[0] is Glc_i in component Glc_i (millimolar).
 * ALGEBRAIC[22] is V_HK in component V_HK (flux).
 * ALGEBRAIC[19] is V_glucose_transport in component V_glucose_transport (flux).
 * ALGEBRAIC[0] is Glc_6_P_g in component Glc_6_P_g (millimolar).
 * CONSTANTS[2] is Keq_PGI in component Glc_6_P_g (dimensionless).
 * ALGEBRAIC[1] is Fru_6_P_g in component Fru_6_P_g (millimolar).
 * STATES[1] is hexose_P_g in component hexose_P_g (millimolar).
 * ALGEBRAIC[29] is V_PFK in component V_PFK (flux).
 * STATES[2] is Fru_1_6_BP_g in component Fru_1_6_BP_g (millimolar).
 * ALGEBRAIC[32] is V_ALD in component V_ALD (flux).
 * ALGEBRAIC[8] is GA_3_P_g in component GA_3_P_g (millimolar).
 * CONSTANTS[3] is Keq_TIM in component DHAP (dimensionless).
 * ALGEBRAIC[6] is DHAP_g in component DHAP_g (millimolar).
 * STATES[3] is triose_P in component triose_P (millimolar).
 * ALGEBRAIC[23] is V_GAPDH in component V_GAPDH (flux).
 * ALGEBRAIC[27] is V_GDH in component V_GDH (flux).
 * ALGEBRAIC[21] is V_GPO in component V_GPO (flux).
 * STATES[4] is one_three_BPGA_g in component one_three_BPGA_g (millimolar).
 * ALGEBRAIC[24] is V_PGK in component V_PGK (flux).
 * ALGEBRAIC[2] is three_PGA in component three_PGA (millimolar).
 * STATES[5] is N in component N (millimolar).
 * CONSTANTS[4] is Keq_PGM in component three_PGA (dimensionless).
 * CONSTANTS[5] is Keq_ENO in component three_PGA (dimensionless).
 * ALGEBRAIC[3] is two_PGA_c in component two_PGA_c (millimolar).
 * ALGEBRAIC[28] is V_PYK in component V_PYK (flux).
 * ALGEBRAIC[4] is PEP_c in component PEP_c (millimolar).
 * STATES[6] is PYR_c in component PYR_c (millimolar).
 * ALGEBRAIC[20] is V_pyruvate_transport in component V_pyruvate_transport (flux).
 * CONSTANTS[6] is glycerol_g in component glycerol_g (millimolar).
 * ALGEBRAIC[5] is DHAP in component DHAP (millimolar).
 * ALGEBRAIC[7] is DHAP_c in component DHAP_c (millimolar).
 * ALGEBRAIC[13] is Gly_3_P in component Gly_3_P (millimolar).
 * ALGEBRAIC[10] is ATP_g in component ATP_g (millimolar).
 * ALGEBRAIC[12] is ADP_g in component ADP_g (millimolar).
 * CONSTANTS[7] is C4 in component C4 (millimolar).
 * ALGEBRAIC[14] is Gly_3_P_c in component Gly_3_P_c (millimolar).
 * ALGEBRAIC[15] is Gly_3_P_g in component Gly_3_P_g (millimolar).
 * ALGEBRAIC[9] is NAD_g in component NAD_g (millimolar).
 * STATES[7] is NADH_g in component NADH_g (millimolar).
 * CONSTANTS[8] is C3 in component C3 (millimolar).
 * STATES[8] is P_g in component P_g (millimolar).
 * ALGEBRAIC[25] is V_GK in component V_GK (flux).
 * STATES[9] is P_c in component P_c (millimolar).
 * ALGEBRAIC[30] is V_ATP_utilisation in component V_ATP_utilisation (flux).
 * CONSTANTS[9] is Keq_AK in component ATP_g (dimensionless).
 * CONSTANTS[10] is C1 in component C1 (millimolar).
 * ALGEBRAIC[11] is ATP_c in component ATP_c (millimolar).
 * CONSTANTS[11] is Keq_AK in component ATP_c (dimensionless).
 * CONSTANTS[12] is C2 in component C2 (millimolar).
 * ALGEBRAIC[16] is ADP_c in component ADP_c (millimolar).
 * ALGEBRAIC[17] is AMP_g in component AMP_g (millimolar).
 * ALGEBRAIC[18] is AMP_c in component AMP_c (millimolar).
 * CONSTANTS[13] is K_Glc in component V_glucose_transport (millimolar).
 * CONSTANTS[14] is alpha in component V_glucose_transport (dimensionless).
 * CONSTANTS[15] is V_glucose_transport_max in component V_glucose_transport (flux).
 * CONSTANTS[16] is K_pyruvate in component V_pyruvate_transport (millimolar).
 * CONSTANTS[17] is V_pyruvate_transport_max in component V_pyruvate_transport (flux).
 * CONSTANTS[18] is K_Gly_3_P in component V_GPO (millimolar).
 * CONSTANTS[19] is V_GPO_max in component V_GPO (flux).
 * CONSTANTS[20] is K_Glc_i in component V_HK (millimolar).
 * CONSTANTS[21] is K_ATP in component V_HK (millimolar).
 * CONSTANTS[22] is K_ADP in component V_HK (millimolar).
 * CONSTANTS[23] is V_HK_max in component V_HK (flux).
 * CONSTANTS[24] is K_NAD in component V_GAPDH (millimolar).
 * CONSTANTS[25] is K_GA_3_P in component V_GAPDH (millimolar).
 * CONSTANTS[26] is K_1_3_BPGA in component V_GAPDH (millimolar).
 * CONSTANTS[27] is K_NADH in component V_GAPDH (millimolar).
 * CONSTANTS[28] is V_GAPDH_max_plus in component V_GAPDH (flux).
 * CONSTANTS[29] is V_GAPDH_max_ratio in component V_GAPDH (dimensionless).
 * CONSTANTS[30] is K_ADP in component V_PGK (millimolar).
 * CONSTANTS[31] is K_1_3_BPGA in component V_PGK (millimolar).
 * CONSTANTS[32] is K_3_PGA in component V_PGK (millimolar).
 * CONSTANTS[33] is K_ATP in component V_PGK (millimolar).
 * CONSTANTS[34] is V_PGK_max_plus in component V_PGK (flux).
 * CONSTANTS[35] is V_PGK_max_ratio in component V_PGK (dimensionless).
 * CONSTANTS[36] is K_ADP in component V_GK (millimolar).
 * CONSTANTS[37] is K_Gly_3_P in component V_GK (millimolar).
 * CONSTANTS[38] is K_glycerol in component V_GK (millimolar).
 * CONSTANTS[39] is K_ATP in component V_GK (millimolar).
 * CONSTANTS[40] is V_GK_max_plus in component V_GK (flux).
 * CONSTANTS[41] is V_GK_max_ratio in component V_GK (dimensionless).
 * CONSTANTS[42] is K_NADH in component V_GDH (millimolar).
 * CONSTANTS[43] is K_Gly_3_P in component V_GDH (millimolar).
 * CONSTANTS[44] is K_DHAP in component V_GDH (millimolar).
 * CONSTANTS[45] is K_NAD in component V_GDH (millimolar).
 * CONSTANTS[46] is V_GDH_max_plus in component V_GDH (flux).
 * CONSTANTS[47] is V_GDH_max_ratio in component V_GDH (dimensionless).
 * CONSTANTS[48] is n in component V_PFK (dimensionless).
 * CONSTANTS[49] is Km_Fru_6_P in component V_PFK (millimolar).
 * CONSTANTS[50] is Km_ATP in component V_PFK (millimolar).
 * CONSTANTS[51] is V_PFK_max in component V_PFK (flux).
 * ALGEBRAIC[26] is Km_PEP in component V_PYK (millimolar).
 * CONSTANTS[52] is Km_ADP in component V_PYK (millimolar).
 * CONSTANTS[53] is n in component V_PYK (dimensionless).
 * CONSTANTS[54] is V_PYK_max in component V_PYK (flux).
 * ALGEBRAIC[31] is Km_Fru_1_6_BP in component V_ALD (millimolar).
 * CONSTANTS[55] is Km_GA_3_P in component V_ALD (millimolar).
 * CONSTANTS[56] is Ki_GA_3_P in component V_ALD (millimolar).
 * CONSTANTS[57] is Km_DHAP in component V_ALD (millimolar).
 * CONSTANTS[58] is V_ALD_max_plus in component V_ALD (flux).
 * CONSTANTS[59] is V_ALD_max_ratio in component V_ALD (dimensionless).
 * CONSTANTS[60] is k in component V_ATP_utilisation (flux).
 * RATES[0] is d/dt Glc_i in component Glc_i (millimolar).
 * RATES[1] is d/dt hexose_P_g in component hexose_P_g (millimolar).
 * RATES[2] is d/dt Fru_1_6_BP_g in component Fru_1_6_BP_g (millimolar).
 * RATES[3] is d/dt triose_P in component triose_P (millimolar).
 * RATES[4] is d/dt one_three_BPGA_g in component one_three_BPGA_g (millimolar).
 * RATES[5] is d/dt N in component N (millimolar).
 * RATES[6] is d/dt PYR_c in component PYR_c (millimolar).
 * RATES[7] is d/dt NADH_g in component NADH_g (millimolar).
 * RATES[8] is d/dt P_g in component P_g (millimolar).
 * RATES[9] is d/dt P_c in component P_c (millimolar).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 22.3;
CONSTANTS[1] = 5.0;
STATES[0] = 0.01;
CONSTANTS[2] = 0.29;
STATES[1] = 0.01;
STATES[2] = 0.01;
CONSTANTS[3] = 0.045;
STATES[3] = 0.01;
STATES[4] = 0.01;
STATES[5] = 0.01;
CONSTANTS[4] = 0.187;
CONSTANTS[5] = 6.7;
STATES[6] = 0.01;
CONSTANTS[6] = 0.00;
CONSTANTS[7] = 120.0;
STATES[7] = 0.01;
CONSTANTS[8] = 4.0;
STATES[8] = 0.01;
STATES[9] = 0.01;
CONSTANTS[9] = 0.442;
CONSTANTS[10] = 3.9;
CONSTANTS[11] = 0.442;
CONSTANTS[12] = 3.9;
CONSTANTS[13] = 2.0;
CONSTANTS[14] = 0.75;
CONSTANTS[15] = 106.2;
CONSTANTS[16] = 1.96;
CONSTANTS[17] = 160.0;
CONSTANTS[18] = 1.7;
CONSTANTS[19] = 200.0;
CONSTANTS[20] = 0.1;
CONSTANTS[21] = 0.116;
CONSTANTS[22] = 0.126;
CONSTANTS[23] = 625.0;
CONSTANTS[24] = 0.45;
CONSTANTS[25] = 0.15;
CONSTANTS[26] = 0.1;
CONSTANTS[27] = 0.02;
CONSTANTS[28] = 1470.0;
CONSTANTS[29] = 0.67;
CONSTANTS[30] = 0.1;
CONSTANTS[31] = 0.05;
CONSTANTS[32] = 1.62;
CONSTANTS[33] = 0.29;
CONSTANTS[34] = 640.0;
CONSTANTS[35] = 0.029;
CONSTANTS[36] = 0.12;
CONSTANTS[37] = 5.1;
CONSTANTS[38] = 0.12;
CONSTANTS[39] = 0.19;
CONSTANTS[40] = 0.0;
CONSTANTS[41] = 167.0;
CONSTANTS[42] = 0.01;
CONSTANTS[43] = 2.0;
CONSTANTS[44] = 0.1;
CONSTANTS[45] = 0.4;
CONSTANTS[46] = 533.0;
CONSTANTS[47] = 0.28;
CONSTANTS[48] = 1.2;
CONSTANTS[49] = 0.82;
CONSTANTS[50] = 2.6E-2;
CONSTANTS[51] = 780.0;
CONSTANTS[52] = 0.114;
CONSTANTS[53] = 2.5;
CONSTANTS[54] = 2.6E3;
CONSTANTS[55] = 6.7E-2;
CONSTANTS[56] = 9.8E-2;
CONSTANTS[57] = 1.5E-2;
CONSTANTS[58] = 780.0;
CONSTANTS[59] = 1.19;
CONSTANTS[60] = 50;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[10] = ((CONSTANTS[10] - - ( STATES[8]*(1.00000 -  4.00000*CONSTANTS[9])))+pow(pow(CONSTANTS[10] - - ( STATES[8]*(1.00000 -  4.00000*CONSTANTS[9])), 2.00000) -  4.00000*(1.00000 -  4.00000*CONSTANTS[9])*( - CONSTANTS[9]*pow(STATES[8], 2.00000)), 0.500000))/( 2.00000*(1.00000 -  4.00000*CONSTANTS[9]));
ALGEBRAIC[12] = STATES[8] -  2.00000*ALGEBRAIC[10];
ALGEBRAIC[22] =  CONSTANTS[23]*(( (ALGEBRAIC[10]/CONSTANTS[21])*(STATES[0]/CONSTANTS[20]))/( (1.00000+ALGEBRAIC[10]/CONSTANTS[21]+ALGEBRAIC[12]/CONSTANTS[22])*(1.00000+STATES[0]/CONSTANTS[20])));
ALGEBRAIC[19] =  CONSTANTS[15]*((CONSTANTS[1] - STATES[0])/(CONSTANTS[13]+CONSTANTS[1]+STATES[0]+ CONSTANTS[14]*CONSTANTS[1]*(STATES[0]/CONSTANTS[13])));
RATES[0] = ALGEBRAIC[19] - ALGEBRAIC[22];
ALGEBRAIC[5] = ( STATES[3]*(1.00000+CONSTANTS[0]))/(1.00000+CONSTANTS[0]+CONSTANTS[3]);
ALGEBRAIC[6] = ALGEBRAIC[5];
ALGEBRAIC[8] =  CONSTANTS[3]*ALGEBRAIC[6];
ALGEBRAIC[9] = CONSTANTS[8] - STATES[7];
ALGEBRAIC[23] =  CONSTANTS[28]*(( (ALGEBRAIC[8]/CONSTANTS[25])*(ALGEBRAIC[9]/CONSTANTS[24] - CONSTANTS[29])*(STATES[4]/CONSTANTS[26])*(STATES[7]/CONSTANTS[27]))/( (1.00000+ALGEBRAIC[8]/CONSTANTS[25]+STATES[4]/CONSTANTS[26])*(1.00000+ALGEBRAIC[9]/CONSTANTS[24]+STATES[7]/CONSTANTS[27])));
ALGEBRAIC[2] = ( STATES[5]*(1.00000+CONSTANTS[0]))/(1.00000+ (1.00000+CONSTANTS[4]+ CONSTANTS[4]*CONSTANTS[5])*CONSTANTS[0]);
ALGEBRAIC[24] =  CONSTANTS[34]*(( (STATES[4]/CONSTANTS[31])*(ALGEBRAIC[12]/CONSTANTS[30] - CONSTANTS[35])*(ALGEBRAIC[2]/CONSTANTS[32])*(ALGEBRAIC[10]/CONSTANTS[33]))/( (1.00000+STATES[4]/CONSTANTS[31]+ALGEBRAIC[2]/CONSTANTS[32])*(1.00000+ALGEBRAIC[12]/CONSTANTS[30]+ALGEBRAIC[10]/CONSTANTS[33])));
RATES[4] = ALGEBRAIC[23] - ALGEBRAIC[24];
ALGEBRAIC[3] =  CONSTANTS[4]*ALGEBRAIC[2];
ALGEBRAIC[4] =  CONSTANTS[5]*ALGEBRAIC[3];
ALGEBRAIC[11] = ((CONSTANTS[12] - - ( STATES[9]*(1.00000 -  4.00000*CONSTANTS[11])))+pow(pow(CONSTANTS[12] - - ( STATES[9]*(1.00000 -  4.00000*CONSTANTS[11])), 2.00000) -  4.00000*(1.00000 -  4.00000*CONSTANTS[11])*( - CONSTANTS[11]*pow(STATES[9], 2.00000)), 0.500000))/( 2.00000*(1.00000 -  4.00000*CONSTANTS[11]));
ALGEBRAIC[16] = STATES[9] -  2.00000*ALGEBRAIC[11];
ALGEBRAIC[26] =  0.340000*(1.00000+ALGEBRAIC[11]/0.570000+ALGEBRAIC[16]/0.640000);
ALGEBRAIC[28] =  CONSTANTS[54]*(( pow(ALGEBRAIC[4]/ALGEBRAIC[26], CONSTANTS[53])*(ALGEBRAIC[16]/CONSTANTS[52]))/( (1.00000+pow(ALGEBRAIC[4]/ALGEBRAIC[26], CONSTANTS[53]))*(1.00000+ALGEBRAIC[16]/CONSTANTS[52])));
RATES[5] = ALGEBRAIC[24] - ALGEBRAIC[28];
ALGEBRAIC[20] =  CONSTANTS[17]*(STATES[6]/CONSTANTS[16])*(1.00000+STATES[6]/CONSTANTS[16]);
RATES[6] = ALGEBRAIC[28] - ALGEBRAIC[20];
ALGEBRAIC[0] = STATES[1]/CONSTANTS[2];
ALGEBRAIC[1] = STATES[1] - ALGEBRAIC[0];
ALGEBRAIC[13] = (CONSTANTS[7] - (ALGEBRAIC[0]+ALGEBRAIC[1]+ 2.00000*STATES[2]+ALGEBRAIC[8]+STATES[4]+ 2.00000*ALGEBRAIC[10]+ALGEBRAIC[12]))/(1.00000+CONSTANTS[0]) - ALGEBRAIC[5];
ALGEBRAIC[15] = ALGEBRAIC[13];
ALGEBRAIC[27] =  CONSTANTS[46]*(( (ALGEBRAIC[6]/CONSTANTS[44])*(STATES[7]/CONSTANTS[42] - CONSTANTS[47])*(ALGEBRAIC[15]/CONSTANTS[43])*(ALGEBRAIC[9]/CONSTANTS[45]))/( (1.00000+ALGEBRAIC[6]/CONSTANTS[44]+ALGEBRAIC[15]/CONSTANTS[43])*(1.00000+STATES[7]/CONSTANTS[42]+ALGEBRAIC[9]/CONSTANTS[45])));
RATES[7] = ALGEBRAIC[23] - ALGEBRAIC[27];
ALGEBRAIC[29] =  CONSTANTS[51]*(( pow(ALGEBRAIC[1]/CONSTANTS[49], CONSTANTS[48])*(ALGEBRAIC[10]/CONSTANTS[50]))/( (1.00000+pow(ALGEBRAIC[1]/CONSTANTS[49], CONSTANTS[48]))*(1.00000+ALGEBRAIC[10]/CONSTANTS[50])));
RATES[1] = ALGEBRAIC[22] - ALGEBRAIC[29];
ALGEBRAIC[25] =  CONSTANTS[40]*(( (ALGEBRAIC[15]/CONSTANTS[37])*(ALGEBRAIC[12]/CONSTANTS[36] - CONSTANTS[41])*(CONSTANTS[6]/CONSTANTS[38])*(ALGEBRAIC[10]/CONSTANTS[39]))/( (1.00000+ALGEBRAIC[15]/CONSTANTS[37]+CONSTANTS[6]/CONSTANTS[38])*(1.00000+ALGEBRAIC[12]/CONSTANTS[36]+ALGEBRAIC[10]/CONSTANTS[39])));
RATES[8] = (ALGEBRAIC[24]+ALGEBRAIC[25]) - (ALGEBRAIC[22]+ALGEBRAIC[29]);
ALGEBRAIC[30] =  CONSTANTS[60]*(ALGEBRAIC[11]/ALGEBRAIC[16]);
RATES[9] = ALGEBRAIC[28] - ALGEBRAIC[30];
ALGEBRAIC[17] = CONSTANTS[10] - (ALGEBRAIC[10]+ALGEBRAIC[12]);
ALGEBRAIC[31] =  0.00900000*(1.00000+ALGEBRAIC[10]/0.680000+ALGEBRAIC[12]/1.51000+ALGEBRAIC[17]/3.65000);
ALGEBRAIC[32] =  CONSTANTS[58]*((STATES[2]/ALGEBRAIC[31] -  CONSTANTS[59]*(( ALGEBRAIC[8]*ALGEBRAIC[6])/( CONSTANTS[55]*CONSTANTS[57])))/(1.00000+STATES[2]/ALGEBRAIC[31]+ALGEBRAIC[8]/CONSTANTS[55]+ALGEBRAIC[6]/CONSTANTS[57]+( STATES[2]*ALGEBRAIC[8])/( ALGEBRAIC[31]*CONSTANTS[56])+( ALGEBRAIC[6]*ALGEBRAIC[8])/( CONSTANTS[57]*CONSTANTS[55])));
RATES[2] = ALGEBRAIC[29] - ALGEBRAIC[32];
ALGEBRAIC[14] = ALGEBRAIC[13];
ALGEBRAIC[21] =  CONSTANTS[19]*(ALGEBRAIC[14]/CONSTANTS[18])*(1.00000+ALGEBRAIC[14]/CONSTANTS[18]);
RATES[3] = ( 2.00000*ALGEBRAIC[32]+ALGEBRAIC[21]) - (ALGEBRAIC[23]+ALGEBRAIC[27]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[10] = ((CONSTANTS[10] - - ( STATES[8]*(1.00000 -  4.00000*CONSTANTS[9])))+pow(pow(CONSTANTS[10] - - ( STATES[8]*(1.00000 -  4.00000*CONSTANTS[9])), 2.00000) -  4.00000*(1.00000 -  4.00000*CONSTANTS[9])*( - CONSTANTS[9]*pow(STATES[8], 2.00000)), 0.500000))/( 2.00000*(1.00000 -  4.00000*CONSTANTS[9]));
ALGEBRAIC[12] = STATES[8] -  2.00000*ALGEBRAIC[10];
ALGEBRAIC[22] =  CONSTANTS[23]*(( (ALGEBRAIC[10]/CONSTANTS[21])*(STATES[0]/CONSTANTS[20]))/( (1.00000+ALGEBRAIC[10]/CONSTANTS[21]+ALGEBRAIC[12]/CONSTANTS[22])*(1.00000+STATES[0]/CONSTANTS[20])));
ALGEBRAIC[19] =  CONSTANTS[15]*((CONSTANTS[1] - STATES[0])/(CONSTANTS[13]+CONSTANTS[1]+STATES[0]+ CONSTANTS[14]*CONSTANTS[1]*(STATES[0]/CONSTANTS[13])));
ALGEBRAIC[5] = ( STATES[3]*(1.00000+CONSTANTS[0]))/(1.00000+CONSTANTS[0]+CONSTANTS[3]);
ALGEBRAIC[6] = ALGEBRAIC[5];
ALGEBRAIC[8] =  CONSTANTS[3]*ALGEBRAIC[6];
ALGEBRAIC[9] = CONSTANTS[8] - STATES[7];
ALGEBRAIC[23] =  CONSTANTS[28]*(( (ALGEBRAIC[8]/CONSTANTS[25])*(ALGEBRAIC[9]/CONSTANTS[24] - CONSTANTS[29])*(STATES[4]/CONSTANTS[26])*(STATES[7]/CONSTANTS[27]))/( (1.00000+ALGEBRAIC[8]/CONSTANTS[25]+STATES[4]/CONSTANTS[26])*(1.00000+ALGEBRAIC[9]/CONSTANTS[24]+STATES[7]/CONSTANTS[27])));
ALGEBRAIC[2] = ( STATES[5]*(1.00000+CONSTANTS[0]))/(1.00000+ (1.00000+CONSTANTS[4]+ CONSTANTS[4]*CONSTANTS[5])*CONSTANTS[0]);
ALGEBRAIC[24] =  CONSTANTS[34]*(( (STATES[4]/CONSTANTS[31])*(ALGEBRAIC[12]/CONSTANTS[30] - CONSTANTS[35])*(ALGEBRAIC[2]/CONSTANTS[32])*(ALGEBRAIC[10]/CONSTANTS[33]))/( (1.00000+STATES[4]/CONSTANTS[31]+ALGEBRAIC[2]/CONSTANTS[32])*(1.00000+ALGEBRAIC[12]/CONSTANTS[30]+ALGEBRAIC[10]/CONSTANTS[33])));
ALGEBRAIC[3] =  CONSTANTS[4]*ALGEBRAIC[2];
ALGEBRAIC[4] =  CONSTANTS[5]*ALGEBRAIC[3];
ALGEBRAIC[11] = ((CONSTANTS[12] - - ( STATES[9]*(1.00000 -  4.00000*CONSTANTS[11])))+pow(pow(CONSTANTS[12] - - ( STATES[9]*(1.00000 -  4.00000*CONSTANTS[11])), 2.00000) -  4.00000*(1.00000 -  4.00000*CONSTANTS[11])*( - CONSTANTS[11]*pow(STATES[9], 2.00000)), 0.500000))/( 2.00000*(1.00000 -  4.00000*CONSTANTS[11]));
ALGEBRAIC[16] = STATES[9] -  2.00000*ALGEBRAIC[11];
ALGEBRAIC[26] =  0.340000*(1.00000+ALGEBRAIC[11]/0.570000+ALGEBRAIC[16]/0.640000);
ALGEBRAIC[28] =  CONSTANTS[54]*(( pow(ALGEBRAIC[4]/ALGEBRAIC[26], CONSTANTS[53])*(ALGEBRAIC[16]/CONSTANTS[52]))/( (1.00000+pow(ALGEBRAIC[4]/ALGEBRAIC[26], CONSTANTS[53]))*(1.00000+ALGEBRAIC[16]/CONSTANTS[52])));
ALGEBRAIC[20] =  CONSTANTS[17]*(STATES[6]/CONSTANTS[16])*(1.00000+STATES[6]/CONSTANTS[16]);
ALGEBRAIC[0] = STATES[1]/CONSTANTS[2];
ALGEBRAIC[1] = STATES[1] - ALGEBRAIC[0];
ALGEBRAIC[13] = (CONSTANTS[7] - (ALGEBRAIC[0]+ALGEBRAIC[1]+ 2.00000*STATES[2]+ALGEBRAIC[8]+STATES[4]+ 2.00000*ALGEBRAIC[10]+ALGEBRAIC[12]))/(1.00000+CONSTANTS[0]) - ALGEBRAIC[5];
ALGEBRAIC[15] = ALGEBRAIC[13];
ALGEBRAIC[27] =  CONSTANTS[46]*(( (ALGEBRAIC[6]/CONSTANTS[44])*(STATES[7]/CONSTANTS[42] - CONSTANTS[47])*(ALGEBRAIC[15]/CONSTANTS[43])*(ALGEBRAIC[9]/CONSTANTS[45]))/( (1.00000+ALGEBRAIC[6]/CONSTANTS[44]+ALGEBRAIC[15]/CONSTANTS[43])*(1.00000+STATES[7]/CONSTANTS[42]+ALGEBRAIC[9]/CONSTANTS[45])));
ALGEBRAIC[29] =  CONSTANTS[51]*(( pow(ALGEBRAIC[1]/CONSTANTS[49], CONSTANTS[48])*(ALGEBRAIC[10]/CONSTANTS[50]))/( (1.00000+pow(ALGEBRAIC[1]/CONSTANTS[49], CONSTANTS[48]))*(1.00000+ALGEBRAIC[10]/CONSTANTS[50])));
ALGEBRAIC[25] =  CONSTANTS[40]*(( (ALGEBRAIC[15]/CONSTANTS[37])*(ALGEBRAIC[12]/CONSTANTS[36] - CONSTANTS[41])*(CONSTANTS[6]/CONSTANTS[38])*(ALGEBRAIC[10]/CONSTANTS[39]))/( (1.00000+ALGEBRAIC[15]/CONSTANTS[37]+CONSTANTS[6]/CONSTANTS[38])*(1.00000+ALGEBRAIC[12]/CONSTANTS[36]+ALGEBRAIC[10]/CONSTANTS[39])));
ALGEBRAIC[30] =  CONSTANTS[60]*(ALGEBRAIC[11]/ALGEBRAIC[16]);
ALGEBRAIC[17] = CONSTANTS[10] - (ALGEBRAIC[10]+ALGEBRAIC[12]);
ALGEBRAIC[31] =  0.00900000*(1.00000+ALGEBRAIC[10]/0.680000+ALGEBRAIC[12]/1.51000+ALGEBRAIC[17]/3.65000);
ALGEBRAIC[32] =  CONSTANTS[58]*((STATES[2]/ALGEBRAIC[31] -  CONSTANTS[59]*(( ALGEBRAIC[8]*ALGEBRAIC[6])/( CONSTANTS[55]*CONSTANTS[57])))/(1.00000+STATES[2]/ALGEBRAIC[31]+ALGEBRAIC[8]/CONSTANTS[55]+ALGEBRAIC[6]/CONSTANTS[57]+( STATES[2]*ALGEBRAIC[8])/( ALGEBRAIC[31]*CONSTANTS[56])+( ALGEBRAIC[6]*ALGEBRAIC[8])/( CONSTANTS[57]*CONSTANTS[55])));
ALGEBRAIC[14] = ALGEBRAIC[13];
ALGEBRAIC[21] =  CONSTANTS[19]*(ALGEBRAIC[14]/CONSTANTS[18])*(1.00000+ALGEBRAIC[14]/CONSTANTS[18]);
ALGEBRAIC[7] = ALGEBRAIC[5];
ALGEBRAIC[18] = CONSTANTS[12] - (ALGEBRAIC[11]+ALGEBRAIC[16]);
}