/*
   There are a total of 53 entries in the algebraic variable array.
   There are a total of 22 entries in each of the rate and state variable arrays.
   There are a total of 83 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (hour).
 * CONSTANTS[0] is eps in component E (dimensionless).
 * CONSTANTS[1] is K1 in component A (per_hour).
 * CONSTANTS[2] is K1a in component A (per_hour).
 * CONSTANTS[3] is K2 in component A (per_hour).
 * CONSTANTS[4] is K2a in component A (per_hour).
 * CONSTANTS[5] is K2aa in component A (per_hour).
 * CONSTANTS[6] is K3 in component A (per_hour).
 * CONSTANTS[7] is K3a in component A (per_hour).
 * CONSTANTS[8] is K4 in component A (per_hour).
 * CONSTANTS[9] is J1 in component A (dimensionless).
 * CONSTANTS[10] is J3 in component A (dimensionless).
 * CONSTANTS[11] is J4 in component A (dimensionless).
 * ALGEBRAIC[1] is v_20 in component A (per_hour).
 * ALGEBRAIC[21] is v_21 in component A (per_hour).
 * ALGEBRAIC[0] is v_42 in component A (per_hour).
 * ALGEBRAIC[22] is v_19 in component A (per_hour).
 * CONSTANTS[12] is GA in component A (dimensionless).
 * CONSTANTS[13] is GB in component A (dimensionless).
 * CONSTANTS[14] is GE in component A (dimensionless).
 * ALGEBRAIC[12] is V_2 in component A (per_hour).
 * ALGEBRAIC[13] is V_4 in component A (per_hour).
 * STATES[0] is Cdh1 in component A (dimensionless).
 * CONSTANTS[15] is Cdc20 in component D (dimensionless).
 * STATES[1] is CYCA in component D (dimensionless).
 * STATES[2] is CYCB in component A (dimensionless).
 * STATES[3] is CYCE in component D (dimensionless).
 * CONSTANTS[16] is k15 in component B (per_hour).
 * CONSTANTS[17] is k16 in component B (per_hour).
 * CONSTANTS[18] is k17 in component B (per_hour).
 * CONSTANTS[19] is k17a in component B (per_hour).
 * CONSTANTS[20] is k18 in component B (per_hour).
 * CONSTANTS[21] is K9 in component B (per_hour).
 * CONSTANTS[22] is J17 in component B (dimensionless).
 * CONSTANTS[23] is J15 in component B (dimensionless).
 * ALGEBRAIC[14] is v_1 in component B (per_hour).
 * ALGEBRAIC[15] is v_2 in component B (per_hour).
 * ALGEBRAIC[2] is v_34 in component B (per_hour).
 * ALGEBRAIC[5] is v_39 in component B (per_hour).
 * ALGEBRAIC[3] is v_41 in component B (per_hour).
 * STATES[4] is DRG in component B (dimensionless).
 * STATES[5] is ERG in component B (dimensionless).
 * CONSTANTS[24] is K19 in component C (per_hour).
 * CONSTANTS[25] is K19a in component C (per_hour).
 * CONSTANTS[26] is K20 in component C (per_hour).
 * CONSTANTS[27] is K21 in component C (dimensionless).
 * CONSTANTS[28] is K22 in component C (per_hour).
 * CONSTANTS[29] is K23a in component C (per_hour).
 * CONSTANTS[30] is K23 in component C (per_hour).
 * CONSTANTS[31] is K26 in component C (per_hour).
 * CONSTANTS[32] is K26R in component C (per_hour).
 * ALGEBRAIC[16] is v_29 in component C (per_hour).
 * ALGEBRAIC[23] is v_30 in component C (per_hour).
 * ALGEBRAIC[25] is v_43 in component C (per_hour).
 * ALGEBRAIC[27] is v_44 in component C (per_hour).
 * ALGEBRAIC[29] is v_45 in component C (per_hour).
 * ALGEBRAIC[31] is v_46 in component C (per_hour).
 * ALGEBRAIC[33] is v_47 in component C (per_hour).
 * ALGEBRAIC[35] is v_48 in component C (per_hour).
 * ALGEBRAIC[37] is v_49 in component C (per_hour).
 * ALGEBRAIC[40] is v_50 in component C (per_hour).
 * ALGEBRAIC[42] is v_51 in component C (per_hour).
 * ALGEBRAIC[44] is v_52 in component C (per_hour).
 * CONSTANTS[33] is LA in component C (dimensionless).
 * CONSTANTS[34] is LB in component C (dimensionless).
 * CONSTANTS[35] is LD in component C (dimensionless).
 * CONSTANTS[36] is LE in component C (dimensionless).
 * CONSTANTS[37] is FE in component C (dimensionless).
 * CONSTANTS[38] is FB in component C (dimensionless).
 * STATES[6] is Rb in component C (dimensionless).
 * STATES[7] is PPRb in component C (dimensionless).
 * STATES[8] is E2F in component C (dimensionless).
 * STATES[9] is PE2F in component C (dimensionless).
 * STATES[10] is E2FRb in component C (dimensionless).
 * STATES[11] is PE2FRb in component C (dimensionless).
 * ALGEBRAIC[4] is PP1A in component C (dimensionless).
 * CONSTANTS[39] is PP1T in component C (dimensionless).
 * CONSTANTS[40] is CYCDT in component C (dimensionless).
 * CONSTANTS[41] is K5 in component D (per_hour).
 * CONSTANTS[42] is K6a in component D (per_hour).
 * CONSTANTS[43] is K6 in component D (per_hour).
 * CONSTANTS[44] is K7a in component D (per_hour).
 * CONSTANTS[45] is K7 in component D (per_hour).
 * CONSTANTS[46] is K8a in component D (per_hour).
 * CONSTANTS[47] is K8 in component D (per_hour).
 * CONSTANTS[48] is K10 in component D (per_hour).
 * CONSTANTS[49] is k24 in component D (per_hour).
 * CONSTANTS[50] is k24r in component D (per_hour).
 * CONSTANTS[51] is K25 in component D (per_hour).
 * CONSTANTS[52] is K25R in component D (per_hour).
 * CONSTANTS[53] is K29 in component D (per_hour).
 * CONSTANTS[54] is K30 in component D (per_hour).
 * CONSTANTS[55] is J8 in component D (dimensionless).
 * ALGEBRAIC[26] is v_3 in component D (per_hour).
 * ALGEBRAIC[28] is v_4 in component D (per_hour).
 * ALGEBRAIC[30] is v_5 in component D (per_hour).
 * ALGEBRAIC[32] is v_6 in component D (per_hour).
 * ALGEBRAIC[34] is v_7 in component D (per_hour).
 * ALGEBRAIC[36] is v_8 in component D (per_hour).
 * ALGEBRAIC[38] is v_9 in component D (per_hour).
 * ALGEBRAIC[41] is v_10 in component D (per_hour).
 * ALGEBRAIC[43] is v_11 in component D (per_hour).
 * ALGEBRAIC[45] is v_12 in component D (per_hour).
 * ALGEBRAIC[46] is v_13 in component D (per_hour).
 * ALGEBRAIC[49] is v_14 in component D (per_hour).
 * ALGEBRAIC[48] is v_15 in component D (per_hour).
 * ALGEBRAIC[51] is v_16 in component D (per_hour).
 * ALGEBRAIC[39] is v_17 in component D (per_hour).
 * ALGEBRAIC[47] is v_18 in component D (per_hour).
 * ALGEBRAIC[50] is v_36 in component D (per_hour).
 * ALGEBRAIC[52] is v_38 in component D (per_hour).
 * CONSTANTS[79] is v_40 in component D (per_hour).
 * CONSTANTS[56] is HA in component D (dimensionless).
 * CONSTANTS[57] is HB in component D (dimensionless).
 * CONSTANTS[58] is HE in component D (dimensionless).
 * CONSTANTS[59] is YE in component D (dimensionless).
 * CONSTANTS[60] is YB in component D (dimensionless).
 * STATES[12] is p27 in component D (dimensionless).
 * ALGEBRAIC[17] is V_6 in component D (per_hour).
 * ALGEBRAIC[24] is V_8 in component D (per_hour).
 * STATES[13] is MASS in component E (dimensionless).
 * STATES[14] is CA in component D (dimensionless).
 * STATES[15] is CD in component D (dimensionless).
 * STATES[16] is CE in component D (dimensionless).
 * STATES[17] is CYCD in component D (dimensionless).
 * CONSTANTS[61] is CYCET in component D (dimensionless).
 * CONSTANTS[62] is K27 in component E (per_hour).
 * CONSTANTS[63] is K28 in component E (per_hour).
 * ALGEBRAIC[6] is v_31 in component E (per_hour).
 * ALGEBRAIC[18] is v_32 in component E (per_hour).
 * ALGEBRAIC[7] is v_33 in component E (per_hour).
 * CONSTANTS[64] is r31switch in component E (dimensionless).
 * CONSTANTS[65] is MU in component E (per_hour).
 * STATES[18] is GM in component E (dimensionless).
 * CONSTANTS[66] is K11a in component F (per_hour).
 * CONSTANTS[67] is K11 in component F (per_hour).
 * CONSTANTS[68] is K12 in component F (per_hour).
 * CONSTANTS[69] is K13 in component F (per_hour).
 * CONSTANTS[70] is K14 in component F (per_hour).
 * CONSTANTS[71] is K31 in component F (per_hour).
 * CONSTANTS[72] is K32 in component F (per_hour).
 * CONSTANTS[73] is K33 in component F (per_hour).
 * CONSTANTS[74] is K34 in component F (per_hour).
 * CONSTANTS[75] is J13 in component F (dimensionless).
 * CONSTANTS[76] is J14 in component F (dimensionless).
 * CONSTANTS[77] is J31 in component F (dimensionless).
 * CONSTANTS[78] is J32 in component F (dimensionless).
 * ALGEBRAIC[8] is v_22 in component F (per_hour).
 * ALGEBRAIC[9] is v_23 in component F (per_hour).
 * ALGEBRAIC[19] is v_24 in component F (per_hour).
 * ALGEBRAIC[11] is v_25 in component F (per_hour).
 * ALGEBRAIC[10] is v_26 in component F (per_hour).
 * CONSTANTS[80] is v_27 in component F (per_hour).
 * CONSTANTS[81] is v_28 in component F (per_hour).
 * ALGEBRAIC[20] is v_35 in component F (per_hour).
 * CONSTANTS[82] is v_37 in component F (per_hour).
 * STATES[19] is IEP in component F (dimensionless).
 * STATES[20] is PPX in component F (dimensionless).
 * STATES[21] is Cdc20T in component F (dimensionless).
 * RATES[0] is d/dt Cdh1 in component A (dimensionless).
 * RATES[2] is d/dt CYCB in component A (dimensionless).
 * RATES[4] is d/dt DRG in component B (dimensionless).
 * RATES[5] is d/dt ERG in component B (dimensionless).
 * RATES[7] is d/dt PPRb in component C (dimensionless).
 * RATES[8] is d/dt E2F in component C (dimensionless).
 * RATES[9] is d/dt PE2F in component C (dimensionless).
 * RATES[6] is d/dt Rb in component C (dimensionless).
 * RATES[10] is d/dt E2FRb in component C (dimensionless).
 * RATES[11] is d/dt PE2FRb in component C (dimensionless).
 * RATES[14] is d/dt CA in component D (dimensionless).
 * RATES[15] is d/dt CD in component D (dimensionless).
 * RATES[16] is d/dt CE in component D (dimensionless).
 * RATES[1] is d/dt CYCA in component D (dimensionless).
 * RATES[17] is d/dt CYCD in component D (dimensionless).
 * RATES[3] is d/dt CYCE in component D (dimensionless).
 * RATES[12] is d/dt p27 in component D (dimensionless).
 * RATES[18] is d/dt GM in component E (dimensionless).
 * RATES[13] is d/dt MASS in component E (dimensionless).
 * RATES[21] is d/dt Cdc20T in component F (dimensionless).
 * RATES[19] is d/dt IEP in component F (dimensionless).
 * RATES[20] is d/dt PPX in component F (dimensionless).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1;
CONSTANTS[1] = 0.6;
CONSTANTS[2] = 0.1;
CONSTANTS[3] = 20;
CONSTANTS[4] = 0.05;
CONSTANTS[5] = 1;
CONSTANTS[6] = 140;
CONSTANTS[7] = 7.5;
CONSTANTS[8] = 40;
CONSTANTS[9] = 0.1;
CONSTANTS[10] = 0.01;
CONSTANTS[11] = 0.04;
CONSTANTS[12] = 0.3;
CONSTANTS[13] = 1;
CONSTANTS[14] = 0;
STATES[0] = 0.000653278;
CONSTANTS[15] = 0.00220177;
STATES[1] = 1.4094;
STATES[2] = 2.72898;
STATES[3] = 0.0229112;
CONSTANTS[16] = 0.25;
CONSTANTS[17] = 0.25;
CONSTANTS[18] = 10;
CONSTANTS[19] = 0.35;
CONSTANTS[20] = 10;
CONSTANTS[21] = 2.5;
CONSTANTS[22] = 0.3;
CONSTANTS[23] = 0.1;
STATES[4] = 0.900533;
STATES[5] = 0.0121809;
CONSTANTS[24] = 20;
CONSTANTS[25] = 0;
CONSTANTS[26] = 10;
CONSTANTS[27] = 1;
CONSTANTS[28] = 1;
CONSTANTS[29] = 0.005;
CONSTANTS[30] = 1;
CONSTANTS[31] = 10000;
CONSTANTS[32] = 200;
CONSTANTS[33] = 3;
CONSTANTS[34] = 5;
CONSTANTS[35] = 3.3;
CONSTANTS[36] = 5;
CONSTANTS[37] = 25;
CONSTANTS[38] = 2;
STATES[6] = 0.000190871;
STATES[7] = 9.97574;
STATES[8] = 0.989986;
STATES[9] = 3.98594;
STATES[10] = 0.00478911;
STATES[11] = 0.0192822;
CONSTANTS[39] = 1;
CONSTANTS[40] = 0.010976;
CONSTANTS[41] = 20;
CONSTANTS[42] = 10;
CONSTANTS[43] = 100;
CONSTANTS[44] = 0;
CONSTANTS[45] = 0.6;
CONSTANTS[46] = 0.1;
CONSTANTS[47] = 2;
CONSTANTS[48] = 5;
CONSTANTS[49] = 1000;
CONSTANTS[50] = 10;
CONSTANTS[51] = 1000;
CONSTANTS[52] = 10;
CONSTANTS[53] = 0.05;
CONSTANTS[54] = 20;
CONSTANTS[55] = 0.1;
CONSTANTS[56] = 0.5;
CONSTANTS[57] = 1;
CONSTANTS[58] = 0.5;
CONSTANTS[59] = 1;
CONSTANTS[60] = 0.05;
STATES[12] = 0.00922806;
STATES[13] = 1.68776;
STATES[14] = 0.0356927;
STATES[15] = 0.010976;
STATES[16] = 0.000542587;
STATES[17] = 0.43929;
CONSTANTS[61] = 0.000542587;
CONSTANTS[62] = 0.2;
CONSTANTS[63] = 0.2;
CONSTANTS[64] = 1;
CONSTANTS[65] = 0.061;
STATES[18] = 1.35565;
CONSTANTS[66] = 0;
CONSTANTS[67] = 1.5;
CONSTANTS[68] = 1.5;
CONSTANTS[69] = 5;
CONSTANTS[70] = 2.5;
CONSTANTS[71] = 0.7;
CONSTANTS[72] = 1.8;
CONSTANTS[73] = 0.05;
CONSTANTS[74] = 0.05;
CONSTANTS[75] = 0.005;
CONSTANTS[76] = 0.005;
CONSTANTS[77] = 0.01;
CONSTANTS[78] = 0.01;
STATES[19] = 0.154655;
STATES[20] = 1;
STATES[21] = 2.36733;
CONSTANTS[79] =  CONSTANTS[0]*CONSTANTS[41];
CONSTANTS[80] = ( CONSTANTS[70]*CONSTANTS[15])/(CONSTANTS[76]+CONSTANTS[15]);
CONSTANTS[81] =  CONSTANTS[68]*CONSTANTS[15];
CONSTANTS[82] =  CONSTANTS[0]*CONSTANTS[73];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[7] =  CONSTANTS[0]*CONSTANTS[65]*STATES[18];
RATES[13] = ALGEBRAIC[7];
ALGEBRAIC[8] =  CONSTANTS[74]*STATES[20];
RATES[20] = CONSTANTS[82] - ALGEBRAIC[8];
ALGEBRAIC[15] =  CONSTANTS[20]*STATES[4];
ALGEBRAIC[3] =  CONSTANTS[0]*(( CONSTANTS[18]*pow(STATES[4]/CONSTANTS[22], 2.00000))/(pow(STATES[4]/CONSTANTS[22], 2.00000)+1.00000)+ CONSTANTS[19]*STATES[5]);
RATES[4] = ALGEBRAIC[3] - ALGEBRAIC[15];
ALGEBRAIC[14] =  CONSTANTS[17]*STATES[5];
ALGEBRAIC[2] = ( CONSTANTS[0]*CONSTANTS[16])/(pow(STATES[4]/CONSTANTS[23], 2.00000)+1.00000);
RATES[5] = ALGEBRAIC[2] - ALGEBRAIC[14];
ALGEBRAIC[6] =  CONSTANTS[62]*STATES[13]*CONSTANTS[64];
ALGEBRAIC[18] =  CONSTANTS[63]*STATES[18];
RATES[18] = ALGEBRAIC[6] - ALGEBRAIC[18];
ALGEBRAIC[11] =  CONSTANTS[68]*STATES[21];
ALGEBRAIC[20] =  CONSTANTS[0]*CONSTANTS[66]+ CONSTANTS[67]*STATES[2];
RATES[21] = ALGEBRAIC[20] - ALGEBRAIC[11];
ALGEBRAIC[9] = ( CONSTANTS[71]*STATES[2]*(1.00000 - STATES[19]))/((CONSTANTS[77] - STATES[19])+1.00000);
ALGEBRAIC[19] = ( CONSTANTS[72]*STATES[20]*STATES[19])/(CONSTANTS[78]+STATES[19]);
RATES[19] = ALGEBRAIC[9] - ALGEBRAIC[19];
ALGEBRAIC[1] = ( (CONSTANTS[7]+ CONSTANTS[6]*CONSTANTS[15])*(1.00000 - STATES[0]))/((CONSTANTS[10] - STATES[0])+1.00000);
ALGEBRAIC[13] =  CONSTANTS[8]*( CONSTANTS[12]*STATES[1]+ CONSTANTS[13]*STATES[2]+ CONSTANTS[14]*STATES[3]);
ALGEBRAIC[21] = ( ALGEBRAIC[13]*STATES[0])/(CONSTANTS[11]+STATES[0]);
RATES[0] = ALGEBRAIC[1] - ALGEBRAIC[21];
ALGEBRAIC[0] =  CONSTANTS[0]*(( CONSTANTS[1]*pow(STATES[2]/CONSTANTS[9], 2.00000))/(pow(STATES[2]/CONSTANTS[9], 2.00000)+1.00000)+CONSTANTS[2]);
ALGEBRAIC[12] =  CONSTANTS[5]*CONSTANTS[15]+ CONSTANTS[4]*(1.00000 - STATES[0])+ CONSTANTS[3]*STATES[0];
ALGEBRAIC[22] =  ALGEBRAIC[12]*STATES[2];
RATES[2] = ALGEBRAIC[0] - ALGEBRAIC[22];
ALGEBRAIC[16] =  STATES[10]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[23] =  STATES[11]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[25] =  STATES[6]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[4] = CONSTANTS[39]/( CONSTANTS[27]*CONSTANTS[37]*(STATES[1]+STATES[3]+ CONSTANTS[38]*STATES[2])+1.00000);
ALGEBRAIC[27] =  STATES[7]*( CONSTANTS[25]*(CONSTANTS[39] - ALGEBRAIC[4])+ CONSTANTS[24]*ALGEBRAIC[4]);
RATES[7] = (ALGEBRAIC[16]+ALGEBRAIC[23]+ALGEBRAIC[25]) - ALGEBRAIC[27];
ALGEBRAIC[29] =  STATES[10]*CONSTANTS[32];
ALGEBRAIC[31] =  STATES[8]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[33] =  STATES[9]*CONSTANTS[28];
ALGEBRAIC[35] =  STATES[8]*STATES[6]*CONSTANTS[31];
RATES[8] = ((ALGEBRAIC[16]+ALGEBRAIC[29]+ALGEBRAIC[33]) - ALGEBRAIC[31]) - ALGEBRAIC[35];
ALGEBRAIC[26] =  CONSTANTS[48]*STATES[15];
ALGEBRAIC[34] =  CONSTANTS[49]*STATES[12]*STATES[17];
ALGEBRAIC[36] =  CONSTANTS[50]*STATES[15];
ALGEBRAIC[17] = CONSTANTS[42]+ CONSTANTS[43]*( CONSTANTS[56]*STATES[1]+ CONSTANTS[57]*STATES[2]+ CONSTANTS[58]*STATES[3]);
ALGEBRAIC[39] =  ALGEBRAIC[17]*STATES[15];
RATES[15] = ((ALGEBRAIC[34] - ALGEBRAIC[36]) - ALGEBRAIC[39]) - ALGEBRAIC[26];
ALGEBRAIC[5] =  CONSTANTS[0]*CONSTANTS[21]*STATES[4];
ALGEBRAIC[28] =  CONSTANTS[48]*STATES[17];
RATES[17] = ((ALGEBRAIC[5]+ALGEBRAIC[39]+ALGEBRAIC[36]) - ALGEBRAIC[34]) - ALGEBRAIC[28];
ALGEBRAIC[37] =  STATES[11]*CONSTANTS[32];
ALGEBRAIC[40] =  STATES[6]*STATES[9]*CONSTANTS[31];
RATES[9] = ((ALGEBRAIC[23]+ALGEBRAIC[37]+ALGEBRAIC[31]) - ALGEBRAIC[33]) - ALGEBRAIC[40];
RATES[6] = (((ALGEBRAIC[27]+ALGEBRAIC[29]+ALGEBRAIC[37]) - ALGEBRAIC[35]) - ALGEBRAIC[40]) - ALGEBRAIC[25];
ALGEBRAIC[42] =  STATES[11]*CONSTANTS[28];
ALGEBRAIC[44] =  STATES[10]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
RATES[10] = (((ALGEBRAIC[42]+ALGEBRAIC[35]) - ALGEBRAIC[44]) - ALGEBRAIC[16]) - ALGEBRAIC[29];
RATES[11] = (((ALGEBRAIC[44]+ALGEBRAIC[40]) - ALGEBRAIC[42]) - ALGEBRAIC[23]) - ALGEBRAIC[37];
ALGEBRAIC[32] =  CONSTANTS[51]*STATES[12]*STATES[1];
ALGEBRAIC[41] =  CONSTANTS[54]*CONSTANTS[15]*STATES[14];
ALGEBRAIC[45] =  CONSTANTS[52]*STATES[14];
ALGEBRAIC[47] =  ALGEBRAIC[17]*STATES[14];
RATES[14] = ((ALGEBRAIC[32] - ALGEBRAIC[45]) - ALGEBRAIC[47]) - ALGEBRAIC[41];
ALGEBRAIC[38] =  CONSTANTS[54]*CONSTANTS[15]*STATES[1];
ALGEBRAIC[50] =  CONSTANTS[0]*CONSTANTS[53]*STATES[8]*STATES[13];
RATES[1] = ((ALGEBRAIC[50] - ALGEBRAIC[38]) - ALGEBRAIC[32])+ALGEBRAIC[45]+ALGEBRAIC[47];
ALGEBRAIC[30] =  CONSTANTS[51]*STATES[12]*STATES[3];
ALGEBRAIC[43] =  CONSTANTS[52]*STATES[16];
ALGEBRAIC[24] = ( CONSTANTS[47]*( CONSTANTS[59]*(STATES[1]+STATES[3])+ CONSTANTS[60]*STATES[2]))/(CONSTANTS[61]+CONSTANTS[55])+CONSTANTS[46];
ALGEBRAIC[46] =  ALGEBRAIC[24]*STATES[16];
ALGEBRAIC[48] =  ALGEBRAIC[17]*STATES[12];
RATES[12] = (((((CONSTANTS[79]+ALGEBRAIC[26]+ALGEBRAIC[36]) - ALGEBRAIC[48]) - ALGEBRAIC[30]) - ALGEBRAIC[32]) - ALGEBRAIC[34])+ALGEBRAIC[43]+ALGEBRAIC[45]+ALGEBRAIC[46]+ALGEBRAIC[41];
ALGEBRAIC[51] =  ALGEBRAIC[17]*STATES[16];
RATES[16] = ((ALGEBRAIC[30] - ALGEBRAIC[43]) - ALGEBRAIC[46]) - ALGEBRAIC[51];
ALGEBRAIC[49] =  ALGEBRAIC[24]*STATES[3];
ALGEBRAIC[52] =  CONSTANTS[0]*(CONSTANTS[44]+ CONSTANTS[45]*STATES[8]);
RATES[3] = ((ALGEBRAIC[52] - ALGEBRAIC[49]) - ALGEBRAIC[30])+ALGEBRAIC[43]+ALGEBRAIC[51];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[7] =  CONSTANTS[0]*CONSTANTS[65]*STATES[18];
ALGEBRAIC[8] =  CONSTANTS[74]*STATES[20];
ALGEBRAIC[15] =  CONSTANTS[20]*STATES[4];
ALGEBRAIC[3] =  CONSTANTS[0]*(( CONSTANTS[18]*pow(STATES[4]/CONSTANTS[22], 2.00000))/(pow(STATES[4]/CONSTANTS[22], 2.00000)+1.00000)+ CONSTANTS[19]*STATES[5]);
ALGEBRAIC[14] =  CONSTANTS[17]*STATES[5];
ALGEBRAIC[2] = ( CONSTANTS[0]*CONSTANTS[16])/(pow(STATES[4]/CONSTANTS[23], 2.00000)+1.00000);
ALGEBRAIC[6] =  CONSTANTS[62]*STATES[13]*CONSTANTS[64];
ALGEBRAIC[18] =  CONSTANTS[63]*STATES[18];
ALGEBRAIC[11] =  CONSTANTS[68]*STATES[21];
ALGEBRAIC[20] =  CONSTANTS[0]*CONSTANTS[66]+ CONSTANTS[67]*STATES[2];
ALGEBRAIC[9] = ( CONSTANTS[71]*STATES[2]*(1.00000 - STATES[19]))/((CONSTANTS[77] - STATES[19])+1.00000);
ALGEBRAIC[19] = ( CONSTANTS[72]*STATES[20]*STATES[19])/(CONSTANTS[78]+STATES[19]);
ALGEBRAIC[1] = ( (CONSTANTS[7]+ CONSTANTS[6]*CONSTANTS[15])*(1.00000 - STATES[0]))/((CONSTANTS[10] - STATES[0])+1.00000);
ALGEBRAIC[13] =  CONSTANTS[8]*( CONSTANTS[12]*STATES[1]+ CONSTANTS[13]*STATES[2]+ CONSTANTS[14]*STATES[3]);
ALGEBRAIC[21] = ( ALGEBRAIC[13]*STATES[0])/(CONSTANTS[11]+STATES[0]);
ALGEBRAIC[0] =  CONSTANTS[0]*(( CONSTANTS[1]*pow(STATES[2]/CONSTANTS[9], 2.00000))/(pow(STATES[2]/CONSTANTS[9], 2.00000)+1.00000)+CONSTANTS[2]);
ALGEBRAIC[12] =  CONSTANTS[5]*CONSTANTS[15]+ CONSTANTS[4]*(1.00000 - STATES[0])+ CONSTANTS[3]*STATES[0];
ALGEBRAIC[22] =  ALGEBRAIC[12]*STATES[2];
ALGEBRAIC[16] =  STATES[10]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[23] =  STATES[11]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[25] =  STATES[6]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[4] = CONSTANTS[39]/( CONSTANTS[27]*CONSTANTS[37]*(STATES[1]+STATES[3]+ CONSTANTS[38]*STATES[2])+1.00000);
ALGEBRAIC[27] =  STATES[7]*( CONSTANTS[25]*(CONSTANTS[39] - ALGEBRAIC[4])+ CONSTANTS[24]*ALGEBRAIC[4]);
ALGEBRAIC[29] =  STATES[10]*CONSTANTS[32];
ALGEBRAIC[31] =  STATES[8]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[33] =  STATES[9]*CONSTANTS[28];
ALGEBRAIC[35] =  STATES[8]*STATES[6]*CONSTANTS[31];
ALGEBRAIC[26] =  CONSTANTS[48]*STATES[15];
ALGEBRAIC[34] =  CONSTANTS[49]*STATES[12]*STATES[17];
ALGEBRAIC[36] =  CONSTANTS[50]*STATES[15];
ALGEBRAIC[17] = CONSTANTS[42]+ CONSTANTS[43]*( CONSTANTS[56]*STATES[1]+ CONSTANTS[57]*STATES[2]+ CONSTANTS[58]*STATES[3]);
ALGEBRAIC[39] =  ALGEBRAIC[17]*STATES[15];
ALGEBRAIC[5] =  CONSTANTS[0]*CONSTANTS[21]*STATES[4];
ALGEBRAIC[28] =  CONSTANTS[48]*STATES[17];
ALGEBRAIC[37] =  STATES[11]*CONSTANTS[32];
ALGEBRAIC[40] =  STATES[6]*STATES[9]*CONSTANTS[31];
ALGEBRAIC[42] =  STATES[11]*CONSTANTS[28];
ALGEBRAIC[44] =  STATES[10]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[32] =  CONSTANTS[51]*STATES[12]*STATES[1];
ALGEBRAIC[41] =  CONSTANTS[54]*CONSTANTS[15]*STATES[14];
ALGEBRAIC[45] =  CONSTANTS[52]*STATES[14];
ALGEBRAIC[47] =  ALGEBRAIC[17]*STATES[14];
ALGEBRAIC[38] =  CONSTANTS[54]*CONSTANTS[15]*STATES[1];
ALGEBRAIC[50] =  CONSTANTS[0]*CONSTANTS[53]*STATES[8]*STATES[13];
ALGEBRAIC[30] =  CONSTANTS[51]*STATES[12]*STATES[3];
ALGEBRAIC[43] =  CONSTANTS[52]*STATES[16];
ALGEBRAIC[24] = ( CONSTANTS[47]*( CONSTANTS[59]*(STATES[1]+STATES[3])+ CONSTANTS[60]*STATES[2]))/(CONSTANTS[61]+CONSTANTS[55])+CONSTANTS[46];
ALGEBRAIC[46] =  ALGEBRAIC[24]*STATES[16];
ALGEBRAIC[48] =  ALGEBRAIC[17]*STATES[12];
ALGEBRAIC[51] =  ALGEBRAIC[17]*STATES[16];
ALGEBRAIC[49] =  ALGEBRAIC[24]*STATES[3];
ALGEBRAIC[52] =  CONSTANTS[0]*(CONSTANTS[44]+ CONSTANTS[45]*STATES[8]);
ALGEBRAIC[10] = ( CONSTANTS[69]*STATES[19]*(STATES[21] - CONSTANTS[15]))/((CONSTANTS[75] - CONSTANTS[15])+STATES[21]);
}