Generated Code

The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
   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[52] is v_21 in component A (per_hour).
 * ALGEBRAIC[0] is v_42 in component A (per_hour).
 * ALGEBRAIC[51] 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[2] is V_2 in component A (per_hour).
 * ALGEBRAIC[3] 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[8] is v_1 in component B (per_hour).
 * ALGEBRAIC[6] is v_2 in component B (per_hour).
 * ALGEBRAIC[7] is v_34 in component B (per_hour).
 * ALGEBRAIC[4] is v_39 in component B (per_hour).
 * ALGEBRAIC[5] 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[10] is v_29 in component C (per_hour).
 * ALGEBRAIC[11] is v_30 in component C (per_hour).
 * ALGEBRAIC[12] is v_43 in component C (per_hour).
 * ALGEBRAIC[13] is v_44 in component C (per_hour).
 * ALGEBRAIC[14] is v_45 in component C (per_hour).
 * ALGEBRAIC[15] is v_46 in component C (per_hour).
 * ALGEBRAIC[16] is v_47 in component C (per_hour).
 * ALGEBRAIC[17] is v_48 in component C (per_hour).
 * ALGEBRAIC[18] is v_49 in component C (per_hour).
 * ALGEBRAIC[19] is v_50 in component C (per_hour).
 * ALGEBRAIC[20] is v_51 in component C (per_hour).
 * ALGEBRAIC[21] 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[9] 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[24] is v_3 in component D (per_hour).
 * ALGEBRAIC[25] is v_4 in component D (per_hour).
 * ALGEBRAIC[26] is v_5 in component D (per_hour).
 * ALGEBRAIC[27] is v_6 in component D (per_hour).
 * ALGEBRAIC[28] is v_7 in component D (per_hour).
 * ALGEBRAIC[29] is v_8 in component D (per_hour).
 * ALGEBRAIC[30] is v_9 in component D (per_hour).
 * ALGEBRAIC[31] is v_10 in component D (per_hour).
 * ALGEBRAIC[32] is v_11 in component D (per_hour).
 * ALGEBRAIC[33] is v_12 in component D (per_hour).
 * ALGEBRAIC[34] is v_13 in component D (per_hour).
 * ALGEBRAIC[35] is v_14 in component D (per_hour).
 * ALGEBRAIC[36] is v_15 in component D (per_hour).
 * ALGEBRAIC[37] is v_16 in component D (per_hour).
 * ALGEBRAIC[38] is v_17 in component D (per_hour).
 * ALGEBRAIC[39] is v_18 in component D (per_hour).
 * ALGEBRAIC[40] is v_36 in component D (per_hour).
 * ALGEBRAIC[41] 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[22] is V_6 in component D (per_hour).
 * ALGEBRAIC[23] 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[42] is v_31 in component E (per_hour).
 * ALGEBRAIC[43] is v_32 in component E (per_hour).
 * ALGEBRAIC[44] 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[45] is v_22 in component F (per_hour).
 * ALGEBRAIC[46] is v_23 in component F (per_hour).
 * ALGEBRAIC[47] is v_24 in component F (per_hour).
 * ALGEBRAIC[48] is v_25 in component F (per_hour).
 * ALGEBRAIC[49] 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[50] 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).
 * There are a total of 0 condition variables.
 */
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];
RATES[0] = 0.1001;
RATES[2] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[6] = 0.1001;
RATES[10] = 0.1001;
RATES[11] = 0.1001;
RATES[14] = 0.1001;
RATES[15] = 0.1001;
RATES[16] = 0.1001;
RATES[1] = 0.1001;
RATES[17] = 0.1001;
RATES[3] = 0.1001;
RATES[12] = 0.1001;
RATES[18] = 0.1001;
RATES[13] = 0.1001;
RATES[21] = 0.1001;
RATES[19] = 0.1001;
RATES[20] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - ALGEBRAIC[1] - ALGEBRAIC[52];
resid[1] = RATES[2] - ALGEBRAIC[0] - ALGEBRAIC[51];
resid[2] = RATES[4] - ALGEBRAIC[5] - ALGEBRAIC[6];
resid[3] = RATES[5] - ALGEBRAIC[7] - ALGEBRAIC[8];
resid[4] = RATES[7] - (ALGEBRAIC[10]+ALGEBRAIC[11]+ALGEBRAIC[12]) - ALGEBRAIC[13];
resid[5] = RATES[8] - ((ALGEBRAIC[10]+ALGEBRAIC[14]+ALGEBRAIC[16]) - ALGEBRAIC[15]) - ALGEBRAIC[17];
resid[6] = RATES[9] - ((ALGEBRAIC[11]+ALGEBRAIC[18]+ALGEBRAIC[15]) - ALGEBRAIC[16]) - ALGEBRAIC[19];
resid[7] = RATES[6] - (((ALGEBRAIC[13]+ALGEBRAIC[14]+ALGEBRAIC[18]) - ALGEBRAIC[17]) - ALGEBRAIC[19]) - ALGEBRAIC[12];
resid[8] = RATES[10] - (((ALGEBRAIC[20]+ALGEBRAIC[17]) - ALGEBRAIC[21]) - ALGEBRAIC[10]) - ALGEBRAIC[14];
resid[9] = RATES[11] - (((ALGEBRAIC[21]+ALGEBRAIC[19]) - ALGEBRAIC[20]) - ALGEBRAIC[11]) - ALGEBRAIC[18];
resid[10] = RATES[14] - ((ALGEBRAIC[27] - ALGEBRAIC[33]) - ALGEBRAIC[39]) - ALGEBRAIC[31];
resid[11] = RATES[15] - ((ALGEBRAIC[28] - ALGEBRAIC[29]) - ALGEBRAIC[38]) - ALGEBRAIC[24];
resid[12] = RATES[16] - ((ALGEBRAIC[26] - ALGEBRAIC[32]) - ALGEBRAIC[34]) - ALGEBRAIC[37];
resid[13] = RATES[1] - ((ALGEBRAIC[40] - ALGEBRAIC[30]) - ALGEBRAIC[27])+ALGEBRAIC[33]+ALGEBRAIC[39];
resid[14] = RATES[17] - ((ALGEBRAIC[4]+ALGEBRAIC[38]+ALGEBRAIC[29]) - ALGEBRAIC[28]) - ALGEBRAIC[25];
resid[15] = RATES[3] - ((ALGEBRAIC[41] - ALGEBRAIC[35]) - ALGEBRAIC[26])+ALGEBRAIC[32]+ALGEBRAIC[37];
resid[16] = RATES[12] - (((((CONSTANTS[79]+ALGEBRAIC[24]+ALGEBRAIC[29]) - ALGEBRAIC[36]) - ALGEBRAIC[26]) - ALGEBRAIC[27]) - ALGEBRAIC[28])+ALGEBRAIC[32]+ALGEBRAIC[33]+ALGEBRAIC[34]+ALGEBRAIC[31];
resid[17] = RATES[18] - ALGEBRAIC[42] - ALGEBRAIC[43];
resid[18] = RATES[13] - ALGEBRAIC[44];
resid[19] = RATES[21] - ALGEBRAIC[50] - ALGEBRAIC[48];
resid[20] = RATES[19] - ALGEBRAIC[46] - ALGEBRAIC[47];
resid[21] = RATES[20] - CONSTANTS[82] - ALGEBRAIC[45];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[49] = ( CONSTANTS[69]*STATES[19]*(STATES[21] - CONSTANTS[15]))/((CONSTANTS[75] - CONSTANTS[15])+STATES[21]);
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
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[1] = ( (CONSTANTS[7]+ CONSTANTS[6]*CONSTANTS[15])*(1.00000 - STATES[0]))/((CONSTANTS[10] - STATES[0])+1.00000);
ALGEBRAIC[4] =  CONSTANTS[0]*CONSTANTS[21]*STATES[4];
ALGEBRAIC[5] =  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[6] =  CONSTANTS[20]*STATES[4];
ALGEBRAIC[7] = ( CONSTANTS[0]*CONSTANTS[16])/(pow(STATES[4]/CONSTANTS[23], 2.00000)+1.00000);
ALGEBRAIC[8] =  CONSTANTS[17]*STATES[5];
ALGEBRAIC[10] =  STATES[10]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[11] =  STATES[11]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[12] =  STATES[6]*CONSTANTS[26]*( CONSTANTS[40]*CONSTANTS[35]+ CONSTANTS[33]*STATES[1]+ CONSTANTS[34]*STATES[2]+ CONSTANTS[36]*STATES[3]);
ALGEBRAIC[9] = CONSTANTS[39]/( CONSTANTS[27]*CONSTANTS[37]*(STATES[1]+STATES[3]+ CONSTANTS[38]*STATES[2])+1.00000);
ALGEBRAIC[13] =  STATES[7]*( CONSTANTS[25]*(CONSTANTS[39] - ALGEBRAIC[9])+ CONSTANTS[24]*ALGEBRAIC[9]);
ALGEBRAIC[14] =  STATES[10]*CONSTANTS[32];
ALGEBRAIC[15] =  STATES[8]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[16] =  STATES[9]*CONSTANTS[28];
ALGEBRAIC[17] =  STATES[8]*STATES[6]*CONSTANTS[31];
ALGEBRAIC[18] =  STATES[11]*CONSTANTS[32];
ALGEBRAIC[19] =  STATES[6]*STATES[9]*CONSTANTS[31];
ALGEBRAIC[20] =  STATES[11]*CONSTANTS[28];
ALGEBRAIC[21] =  STATES[10]*(CONSTANTS[29]+ CONSTANTS[30]*(STATES[1]+STATES[2]));
ALGEBRAIC[24] =  CONSTANTS[48]*STATES[15];
ALGEBRAIC[25] =  CONSTANTS[48]*STATES[17];
ALGEBRAIC[26] =  CONSTANTS[51]*STATES[12]*STATES[3];
ALGEBRAIC[27] =  CONSTANTS[51]*STATES[12]*STATES[1];
ALGEBRAIC[28] =  CONSTANTS[49]*STATES[12]*STATES[17];
ALGEBRAIC[29] =  CONSTANTS[50]*STATES[15];
ALGEBRAIC[30] =  CONSTANTS[54]*CONSTANTS[15]*STATES[1];
ALGEBRAIC[31] =  CONSTANTS[54]*CONSTANTS[15]*STATES[14];
ALGEBRAIC[32] =  CONSTANTS[52]*STATES[16];
ALGEBRAIC[33] =  CONSTANTS[52]*STATES[14];
ALGEBRAIC[23] = ( CONSTANTS[47]*( CONSTANTS[59]*(STATES[1]+STATES[3])+ CONSTANTS[60]*STATES[2]))/(CONSTANTS[61]+CONSTANTS[55])+CONSTANTS[46];
ALGEBRAIC[34] =  ALGEBRAIC[23]*STATES[16];
ALGEBRAIC[35] =  ALGEBRAIC[23]*STATES[3];
ALGEBRAIC[22] = CONSTANTS[42]+ CONSTANTS[43]*( CONSTANTS[56]*STATES[1]+ CONSTANTS[57]*STATES[2]+ CONSTANTS[58]*STATES[3]);
ALGEBRAIC[36] =  ALGEBRAIC[22]*STATES[12];
ALGEBRAIC[37] =  ALGEBRAIC[22]*STATES[16];
ALGEBRAIC[38] =  ALGEBRAIC[22]*STATES[15];
ALGEBRAIC[39] =  ALGEBRAIC[22]*STATES[14];
ALGEBRAIC[40] =  CONSTANTS[0]*CONSTANTS[53]*STATES[8]*STATES[13];
ALGEBRAIC[41] =  CONSTANTS[0]*(CONSTANTS[44]+ CONSTANTS[45]*STATES[8]);
ALGEBRAIC[42] =  CONSTANTS[62]*STATES[13]*CONSTANTS[64];
ALGEBRAIC[43] =  CONSTANTS[63]*STATES[18];
ALGEBRAIC[44] =  CONSTANTS[0]*CONSTANTS[65]*STATES[18];
ALGEBRAIC[45] =  CONSTANTS[74]*STATES[20];
ALGEBRAIC[46] = ( CONSTANTS[71]*STATES[2]*(1.00000 - STATES[19]))/((CONSTANTS[77] - STATES[19])+1.00000);
ALGEBRAIC[47] = ( CONSTANTS[72]*STATES[20]*STATES[19])/(CONSTANTS[78]+STATES[19]);
ALGEBRAIC[48] =  CONSTANTS[68]*STATES[21];
ALGEBRAIC[50] =  CONSTANTS[0]*CONSTANTS[66]+ CONSTANTS[67]*STATES[2];
ALGEBRAIC[2] =  CONSTANTS[5]*CONSTANTS[15]+ CONSTANTS[4]*(1.00000 - STATES[0])+ CONSTANTS[3]*STATES[0];
ALGEBRAIC[51] =  ALGEBRAIC[2]*STATES[2];
ALGEBRAIC[3] =  CONSTANTS[8]*( CONSTANTS[12]*STATES[1]+ CONSTANTS[13]*STATES[2]+ CONSTANTS[14]*STATES[3]);
ALGEBRAIC[52] = ( ALGEBRAIC[3]*STATES[0])/(CONSTANTS[11]+STATES[0]);
}
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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}