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 32 entries in the algebraic variable array.
   There are a total of 24 entries in each of the rate and state variable arrays.
   There are a total of 142 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * CONSTANTS[0] is k_rel_TPI in component parameters (dimensionless).
 * CONSTANTS[1] is k_rel_GAPDH in component parameters (dimensionless).
 * CONSTANTS[2] is SUMAXP in component parameters (mM).
 * CONSTANTS[3] is CO2 in component extracellular (mM).
 * CONSTANTS[4] is ETOH in component extracellular (mM).
 * CONSTANTS[5] is SUCC in component extracellular (mM).
 * CONSTANTS[6] is GLY in component extracellular (mM).
 * CONSTANTS[7] is GLCo in component extracellular (mM).
 * CONSTANTS[8] is cytoplasm in component cytoplasm (liter).
 * STATES[0] is P in component cytoplasm (mM).
 * STATES[1] is G6P in component cytoplasm (mM).
 * STATES[2] is F6P in component cytoplasm (mM).
 * STATES[3] is F16P in component cytoplasm (mM).
 * STATES[4] is NADH in component cytoplasm (mM).
 * STATES[5] is NAD in component cytoplasm (mM).
 * STATES[6] is BPG in component cytoplasm (mM).
 * STATES[7] is P3G in component cytoplasm (mM).
 * STATES[8] is P2G in component cytoplasm (mM).
 * STATES[9] is PEP in component cytoplasm (mM).
 * STATES[10] is PYR in component cytoplasm (mM).
 * STATES[11] is ACE in component cytoplasm (mM).
 * CONSTANTS[9] is X in component cytoplasm (mM).
 * STATES[12] is GA3P in component cytoplasm (mM).
 * STATES[13] is DHAP in component cytoplasm (mM).
 * STATES[14] is D6PGluconoLactone in component cytoplasm (mM).
 * STATES[15] is D6PGluconate in component cytoplasm (mM).
 * STATES[16] is NADP in component cytoplasm (mM).
 * STATES[17] is NADPH in component cytoplasm (mM).
 * STATES[18] is Ribulose5P in component cytoplasm (mM).
 * STATES[19] is Ribose5P in component cytoplasm (mM).
 * STATES[20] is Xyl5P in component cytoplasm (mM).
 * STATES[21] is Seduhept7P in component cytoplasm (mM).
 * STATES[22] is Erythrose4P in component cytoplasm (mM).
 * STATES[23] is GLCi in component cytoplasm (mM).
 * CONSTANTS[10] is F26BP in component cytoplasm (mM).
 * ALGEBRAIC[1] is vGLK in component vGLK (mmol_per_min).
 * ALGEBRAIC[2] is vPGI in component vPGI (mmol_per_min).
 * ALGEBRAIC[31] is vPFK in component vPFK (mmol_per_min).
 * ALGEBRAIC[8] is vALD in component vALD (mmol_per_min).
 * ALGEBRAIC[9] is vG3PDH in component vG3PDH (mmol_per_min).
 * ALGEBRAIC[10] is vGAPDH in component vGAPDH (mmol_per_min).
 * ALGEBRAIC[11] is vPGK in component vPGK (mmol_per_min).
 * ALGEBRAIC[12] is vPGM in component vPGM (mmol_per_min).
 * ALGEBRAIC[13] is vENO in component vENO (mmol_per_min).
 * ALGEBRAIC[14] is vPYK in component vPYK (mmol_per_min).
 * ALGEBRAIC[15] is vPDC in component vPDC (mmol_per_min).
 * ALGEBRAIC[16] is vSUC in component vSUC (mmol_per_min).
 * ALGEBRAIC[17] is vADH in component vADH (mmol_per_min).
 * ALGEBRAIC[18] is vATP in component vATP (mmol_per_min).
 * ALGEBRAIC[19] is vTPI in component vTPI (mmol_per_min).
 * ALGEBRAIC[20] is vG6PDH in component vG6PDH (mmol_per_min).
 * ALGEBRAIC[21] is v6PGL in component v6PGL (mmol_per_min).
 * ALGEBRAIC[22] is vGluDH in component vGluDH (mmol_per_min).
 * ALGEBRAIC[23] is vPPI in component vPPI (mmol_per_min).
 * ALGEBRAIC[24] is vTransk1 in component vTransk1 (mmol_per_min).
 * ALGEBRAIC[25] is vR5PI in component vR5PI (mmol_per_min).
 * ALGEBRAIC[26] is vTransald in component vTransald (mmol_per_min).
 * ALGEBRAIC[27] is vTransk2 in component vTransk2 (mmol_per_min).
 * ALGEBRAIC[28] is vNADPH in component vNADPH (mmol_per_min).
 * ALGEBRAIC[29] is vGLT in component vGLT (mmol_per_min).
 * ALGEBRAIC[0] is ratio_NADPH_NADP in component rules (dimensionless).
 * CONSTANTS[11] is VmGLK in component vGLK (mM_per_min).
 * CONSTANTS[12] is KeqAK in component vGLK (dimensionless).
 * CONSTANTS[13] is KeqGLK in component vGLK (dimensionless).
 * CONSTANTS[14] is KmGLKATP in component vGLK (mM).
 * CONSTANTS[15] is KmGLKGLCi in component vGLK (mM).
 * CONSTANTS[16] is KmGLKG6P in component vGLK (mM).
 * CONSTANTS[17] is KmGLKADP in component vGLK (mM).
 * CONSTANTS[18] is VmPGI in component vPGI (mM_per_min).
 * CONSTANTS[19] is KmPGIG6P in component vPGI (mM).
 * CONSTANTS[20] is KeqPGI in component vPGI (dimensionless).
 * CONSTANTS[21] is KmPGIF6P in component vPGI (mM).
 * ALGEBRAIC[3] is numerator1 in component vPFK (mmol_mM2_per_minute).
 * ALGEBRAIC[4] is numerator2 in component vPFK (dimensionless).
 * ALGEBRAIC[30] is denominator in component vPFK (mM2).
 * ALGEBRAIC[5] is denom1 in component vPFK (dimensionless).
 * ALGEBRAIC[6] is denom2 in component vPFK (dimensionless).
 * ALGEBRAIC[7] is denom3 in component vPFK (dimensionless).
 * CONSTANTS[22] is gR in component vPFK (dimensionless).
 * CONSTANTS[23] is VmPFK in component vPFK (mM_per_min).
 * CONSTANTS[24] is KeqAK in component vPFK (dimensionless).
 * CONSTANTS[25] is KmPFKF6P in component vPFK (mM).
 * CONSTANTS[26] is KmPFKATP in component vPFK (mM).
 * CONSTANTS[27] is L0 in component vPFK (dimensionless).
 * CONSTANTS[28] is CPFKF26BP in component vPFK (dimensionless).
 * CONSTANTS[29] is KPFKF26BP in component vPFK (mM).
 * CONSTANTS[30] is CPFKF16BP in component vPFK (dimensionless).
 * CONSTANTS[31] is KPFKF16BP in component vPFK (mM).
 * CONSTANTS[32] is CPFKAMP in component vPFK (dimensionless).
 * CONSTANTS[33] is KPFKAMP in component vPFK (mM).
 * CONSTANTS[34] is CiPFKATP in component vPFK (dimensionless).
 * CONSTANTS[35] is KiPFKATP in component vPFK (mM).
 * CONSTANTS[36] is CPFKATP in component vPFK (dimensionless).
 * CONSTANTS[37] is VmALD in component vALD (mM_per_min).
 * CONSTANTS[38] is KeqTPI in component vALD (dimensionless).
 * CONSTANTS[39] is KeqALD in component vALD (mM).
 * CONSTANTS[40] is KmALDF16P in component vALD (mM).
 * CONSTANTS[41] is KmALDDHAP in component vALD (mM).
 * CONSTANTS[42] is KmALDGAP in component vALD (mM).
 * CONSTANTS[43] is KmALDGAPi in component vALD (mM).
 * CONSTANTS[44] is VmG3PDH in component vG3PDH (mM_per_min).
 * CONSTANTS[45] is KeqG3PDH in component vG3PDH (dimensionless).
 * CONSTANTS[46] is KeqTPI in component vG3PDH (dimensionless).
 * CONSTANTS[47] is KmG3PDHDHAP in component vG3PDH (mM).
 * CONSTANTS[48] is KmG3PDHNADH in component vG3PDH (mM).
 * CONSTANTS[49] is KmG3PDHNAD in component vG3PDH (mM).
 * CONSTANTS[50] is KmG3PDHGLY in component vG3PDH (mM).
 * CONSTANTS[51] is VmGAPDHr in component vGAPDH (mM_per_min).
 * CONSTANTS[52] is KmGAPDHBPG in component vGAPDH (mM).
 * CONSTANTS[53] is KmGAPDHNADH in component vGAPDH (mM).
 * CONSTANTS[54] is KeqTPI in component vGAPDH (dimensionless).
 * CONSTANTS[55] is VmGAPDHf in component vGAPDH (mM_per_min).
 * CONSTANTS[56] is KmGAPDHGAP in component vGAPDH (mM).
 * CONSTANTS[57] is KmGAPDHNAD in component vGAPDH (mM).
 * CONSTANTS[58] is KeqGAPDH in component vGAPDH (dimensionless).
 * CONSTANTS[59] is VmPGK in component vPGK (mM_per_min).
 * CONSTANTS[60] is KeqPGK in component vPGK (dimensionless).
 * CONSTANTS[61] is KeqAK in component vPGK (dimensionless).
 * CONSTANTS[62] is KmPGKATP in component vPGK (mM).
 * CONSTANTS[63] is KmPGKP3G in component vPGK (mM).
 * CONSTANTS[64] is KmPGKADP in component vPGK (mM).
 * CONSTANTS[65] is KmPGKBPG in component vPGK (mM).
 * CONSTANTS[66] is VmPGM in component vPGM (mM_per_min).
 * CONSTANTS[67] is KmPGMP3G in component vPGM (mM).
 * CONSTANTS[68] is KeqPGM in component vPGM (dimensionless).
 * CONSTANTS[69] is KmPGMP2G in component vPGM (mM).
 * CONSTANTS[70] is VmENO in component vENO (mM_per_min).
 * CONSTANTS[71] is KmENOP2G in component vENO (mM).
 * CONSTANTS[72] is KeqENO in component vENO (dimensionless).
 * CONSTANTS[73] is KmENOPEP in component vENO (mM).
 * CONSTANTS[74] is VmPYK in component vPYK (mM_per_min).
 * CONSTANTS[75] is KmPYKPEP in component vPYK (mM).
 * CONSTANTS[76] is KmPYKADP in component vPYK (mM).
 * CONSTANTS[77] is KeqAK in component vPYK (dimensionless).
 * CONSTANTS[78] is KeqPYK in component vPYK (dimensionless).
 * CONSTANTS[79] is KmPYKPYR in component vPYK (mM).
 * CONSTANTS[80] is KmPYKATP in component vPYK (mM).
 * CONSTANTS[81] is VmPDC in component vPDC (mM_per_min).
 * CONSTANTS[82] is nPDC in component vPDC (dimensionless).
 * CONSTANTS[83] is KmPDCPYR in component vPDC (mM).
 * CONSTANTS[84] is KSUCC in component vSUC (per_min).
 * CONSTANTS[85] is VmADH in component vADH (mM_per_min).
 * CONSTANTS[86] is KiADHNAD in component vADH (mM).
 * CONSTANTS[87] is KmADHETOH in component vADH (mM).
 * CONSTANTS[88] is KeqADH in component vADH (dimensionless).
 * CONSTANTS[89] is KmADHNAD in component vADH (mM).
 * CONSTANTS[90] is KmADHNADH in component vADH (mM).
 * CONSTANTS[91] is KiADHNADH in component vADH (mM).
 * CONSTANTS[92] is KmADHACE in component vADH (mM).
 * CONSTANTS[93] is KiADHACE in component vADH (mM).
 * CONSTANTS[94] is KiADHETOH in component vADH (mM).
 * CONSTANTS[95] is KATPASE in component vATP (per_min).
 * CONSTANTS[96] is KeqAK in component vATP (dimensionless).
 * CONSTANTS[97] is KmGA3P in component vTPI (mM).
 * CONSTANTS[98] is KmDHAP in component vTPI (mM).
 * CONSTANTS[99] is VmDHAP in component vTPI (mM_per_min).
 * CONSTANTS[100] is VmGA3P in component vTPI (mM_per_min).
 * CONSTANTS[101] is VmG6PDH in component vG6PDH (mM_per_min).
 * CONSTANTS[102] is KmG6P in component vG6PDH (mM).
 * CONSTANTS[103] is KmNADP in component vG6PDH (mM).
 * CONSTANTS[104] is KiNADPH in component vG6PDH (mM).
 * CONSTANTS[105] is Vm6PGL in component v6PGL (mM_per_min).
 * CONSTANTS[106] is Km6PGL in component v6PGL (mM).
 * CONSTANTS[107] is VmGluDH in component vGluDH (mM_per_min).
 * CONSTANTS[108] is KmGluconate in component vGluDH (mM).
 * CONSTANTS[109] is KmNADP in component vGluDH (mM).
 * CONSTANTS[110] is KiNADPH in component vGluDH (mM).
 * CONSTANTS[111] is VmPPIf in component vPPI (mM_per_min).
 * CONSTANTS[112] is VmPPIr in component vPPI (mM_per_min).
 * CONSTANTS[113] is KmRibu5P in component vPPI (mM).
 * CONSTANTS[114] is KmRibo5P in component vPPI (mM).
 * CONSTANTS[115] is VmTransk1f in component vTransk1 (mM_per_min).
 * CONSTANTS[116] is VmTransk1r in component vTransk1 (mM_per_min).
 * CONSTANTS[117] is KmRibose5P in component vTransk1 (mM).
 * CONSTANTS[118] is KmXyl5P in component vTransk1 (mM).
 * CONSTANTS[119] is KmGA3P in component vTransk1 (mM).
 * CONSTANTS[120] is KmSeduhept in component vTransk1 (mM).
 * CONSTANTS[121] is VmR5PIf in component vR5PI (mM_per_min).
 * CONSTANTS[122] is VmR5PIr in component vR5PI (mM_per_min).
 * CONSTANTS[123] is KmRibu5P in component vR5PI (mM).
 * CONSTANTS[124] is KmXyl in component vR5PI (mM).
 * CONSTANTS[125] is VmTransaldf in component vTransald (mM_per_min).
 * CONSTANTS[126] is VmTransaldr in component vTransald (mM_per_min).
 * CONSTANTS[127] is KmGA3P in component vTransald (mM).
 * CONSTANTS[128] is KmSeduhept in component vTransald (mM).
 * CONSTANTS[129] is KmF6P in component vTransald (mM).
 * CONSTANTS[130] is KmEry4P in component vTransald (mM).
 * CONSTANTS[131] is VmTransk2f in component vTransk2 (mM_per_min).
 * CONSTANTS[132] is VmTransk2r in component vTransk2 (mM_per_min).
 * CONSTANTS[133] is KmXyl5P in component vTransk2 (mM).
 * CONSTANTS[134] is KmEry4P in component vTransk2 (mM).
 * CONSTANTS[135] is KmF6P in component vTransk2 (mM).
 * CONSTANTS[136] is KmGA3P in component vTransk2 (mM).
 * CONSTANTS[137] is kNADPH in component vNADPH (per_min).
 * CONSTANTS[138] is VmGLT in component vGLT (mM_per_min).
 * CONSTANTS[139] is KeqGLT in component vGLT (dimensionless).
 * CONSTANTS[140] is KmGLTGLCo in component vGLT (mM).
 * CONSTANTS[141] is KmGLTGLCi in component vGLT (mM).
 * RATES[0] is d/dt P in component cytoplasm (mM).
 * RATES[1] is d/dt G6P in component cytoplasm (mM).
 * RATES[2] is d/dt F6P in component cytoplasm (mM).
 * RATES[3] is d/dt F16P in component cytoplasm (mM).
 * RATES[4] is d/dt NADH in component cytoplasm (mM).
 * RATES[5] is d/dt NAD in component cytoplasm (mM).
 * RATES[6] is d/dt BPG in component cytoplasm (mM).
 * RATES[7] is d/dt P3G in component cytoplasm (mM).
 * RATES[8] is d/dt P2G in component cytoplasm (mM).
 * RATES[9] is d/dt PEP in component cytoplasm (mM).
 * RATES[10] is d/dt PYR in component cytoplasm (mM).
 * RATES[11] is d/dt ACE in component cytoplasm (mM).
 * RATES[12] is d/dt GA3P in component cytoplasm (mM).
 * RATES[13] is d/dt DHAP in component cytoplasm (mM).
 * RATES[14] is d/dt D6PGluconoLactone in component cytoplasm (mM).
 * RATES[15] is d/dt D6PGluconate in component cytoplasm (mM).
 * RATES[16] is d/dt NADP in component cytoplasm (mM).
 * RATES[17] is d/dt NADPH in component cytoplasm (mM).
 * RATES[18] is d/dt Ribulose5P in component cytoplasm (mM).
 * RATES[19] is d/dt Ribose5P in component cytoplasm (mM).
 * RATES[20] is d/dt Xyl5P in component cytoplasm (mM).
 * RATES[21] is d/dt Seduhept7P in component cytoplasm (mM).
 * RATES[22] is d/dt Erythrose4P in component cytoplasm (mM).
 * RATES[23] is d/dt GLCi in component cytoplasm (mM).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1;
CONSTANTS[1] = 1;
CONSTANTS[2] = 4.1;
CONSTANTS[3] = 1;
CONSTANTS[4] = 50;
CONSTANTS[5] = 0.1;
CONSTANTS[6] = 0.15;
CONSTANTS[7] = 50;
CONSTANTS[8] = 1;
STATES[0] = 5;
STATES[1] = 1.39;
STATES[2] = 0.28;
STATES[3] = 0.1;
STATES[4] = 0.39;
STATES[5] = 1.2;
STATES[6] = 0.1;
STATES[7] = 0.1;
STATES[8] = 0.1;
STATES[9] = 0.1;
STATES[10] = 3.36;
STATES[11] = 0.04;
CONSTANTS[9] = 0.1;
STATES[12] = 0.05;
STATES[13] = 1;
STATES[14] = 0.1;
STATES[15] = 0.1;
STATES[16] = 0.4;
STATES[17] = 1.6;
STATES[18] = 0.1;
STATES[19] = 0.1;
STATES[20] = 0.1;
STATES[21] = 0.1;
STATES[22] = 0;
STATES[23] = 0.087;
CONSTANTS[10] = 0.02;
CONSTANTS[11] = 226.452;
CONSTANTS[12] = 0.45;
CONSTANTS[13] = 3800;
CONSTANTS[14] = 0.15;
CONSTANTS[15] = 0.08;
CONSTANTS[16] = 30;
CONSTANTS[17] = 0.23;
CONSTANTS[18] = 339.677;
CONSTANTS[19] = 1.4;
CONSTANTS[20] = 0.314;
CONSTANTS[21] = 0.3;
CONSTANTS[22] = 5.12;
CONSTANTS[23] = 182.903;
CONSTANTS[24] = 0.45;
CONSTANTS[25] = 0.1;
CONSTANTS[26] = 0.71;
CONSTANTS[27] = 0.66;
CONSTANTS[28] = 0.0174;
CONSTANTS[29] = 0.000682;
CONSTANTS[30] = 0.397;
CONSTANTS[31] = 0.111;
CONSTANTS[32] = 0.0845;
CONSTANTS[33] = 0.0995;
CONSTANTS[34] = 100;
CONSTANTS[35] = 0.65;
CONSTANTS[36] = 3;
CONSTANTS[37] = 322.258;
CONSTANTS[38] = 0.045;
CONSTANTS[39] = 0.069;
CONSTANTS[40] = 0.3;
CONSTANTS[41] = 2.4;
CONSTANTS[42] = 2;
CONSTANTS[43] = 10;
CONSTANTS[44] = 70.15;
CONSTANTS[45] = 4300;
CONSTANTS[46] = 0.045;
CONSTANTS[47] = 0.4;
CONSTANTS[48] = 0.023;
CONSTANTS[49] = 0.93;
CONSTANTS[50] = 1;
CONSTANTS[51] = 6549.68;
CONSTANTS[52] = 0.0098;
CONSTANTS[53] = 0.06;
CONSTANTS[54] = 0.045;
CONSTANTS[55] = 1184.52;
CONSTANTS[56] = 0.21;
CONSTANTS[57] = 0.09;
CONSTANTS[58] = 0.005;
CONSTANTS[59] = 1306.45;
CONSTANTS[60] = 3200;
CONSTANTS[61] = 0.45;
CONSTANTS[62] = 0.3;
CONSTANTS[63] = 0.53;
CONSTANTS[64] = 0.2;
CONSTANTS[65] = 0.003;
CONSTANTS[66] = 2525.81;
CONSTANTS[67] = 1.2;
CONSTANTS[68] = 0.19;
CONSTANTS[69] = 0.08;
CONSTANTS[70] = 365.806;
CONSTANTS[71] = 0.04;
CONSTANTS[72] = 6.7;
CONSTANTS[73] = 0.5;
CONSTANTS[74] = 1088.71;
CONSTANTS[75] = 0.14;
CONSTANTS[76] = 0.53;
CONSTANTS[77] = 0.45;
CONSTANTS[78] = 6500;
CONSTANTS[79] = 21;
CONSTANTS[80] = 1.5;
CONSTANTS[81] = 174.194;
CONSTANTS[82] = 1.9;
CONSTANTS[83] = 4.33;
CONSTANTS[84] = 21.4;
CONSTANTS[85] = 810;
CONSTANTS[86] = 0.92;
CONSTANTS[87] = 17;
CONSTANTS[88] = 6.9e-5;
CONSTANTS[89] = 0.17;
CONSTANTS[90] = 0.11;
CONSTANTS[91] = 0.031;
CONSTANTS[92] = 1.11;
CONSTANTS[93] = 1.1;
CONSTANTS[94] = 90;
CONSTANTS[95] = 39.5;
CONSTANTS[96] = 0.45;
CONSTANTS[97] = 1.27;
CONSTANTS[98] = 1.23;
CONSTANTS[99] = 10900;
CONSTANTS[100] = 555;
CONSTANTS[101] = 4;
CONSTANTS[102] = 0.04;
CONSTANTS[103] = 0.02;
CONSTANTS[104] = 0.017;
CONSTANTS[105] = 4;
CONSTANTS[106] = 0.8;
CONSTANTS[107] = 4;
CONSTANTS[108] = 0.02;
CONSTANTS[109] = 0.03;
CONSTANTS[110] = 0.03;
CONSTANTS[111] = 3458;
CONSTANTS[112] = 3458;
CONSTANTS[113] = 1.6;
CONSTANTS[114] = 1.6;
CONSTANTS[115] = 4;
CONSTANTS[116] = 2;
CONSTANTS[117] = 0.1;
CONSTANTS[118] = 0.15;
CONSTANTS[119] = 0.1;
CONSTANTS[120] = 0.15;
CONSTANTS[121] = 1039;
CONSTANTS[122] = 1039;
CONSTANTS[123] = 1.5;
CONSTANTS[124] = 1.5;
CONSTANTS[125] = 55;
CONSTANTS[126] = 10;
CONSTANTS[127] = 0.22;
CONSTANTS[128] = 0.18;
CONSTANTS[129] = 0.32;
CONSTANTS[130] = 0.018;
CONSTANTS[131] = 3.2;
CONSTANTS[132] = 43;
CONSTANTS[133] = 0.16;
CONSTANTS[134] = 0.09;
CONSTANTS[135] = 1.1;
CONSTANTS[136] = 2.1;
CONSTANTS[137] = 2;
CONSTANTS[138] = 97.264;
CONSTANTS[139] = 1;
CONSTANTS[140] = 1.1918;
CONSTANTS[141] = 1.1918;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[13] = 0.1001;
RATES[14] = 0.1001;
RATES[15] = 0.1001;
RATES[16] = 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;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - (((ALGEBRAIC[11] - (ALGEBRAIC[1]+ALGEBRAIC[31]))+ALGEBRAIC[14]) - ALGEBRAIC[18])/CONSTANTS[8];
resid[1] = RATES[1] - ((ALGEBRAIC[1] - ALGEBRAIC[2]) - ALGEBRAIC[20])/CONSTANTS[8];
resid[2] = RATES[2] - ((ALGEBRAIC[2] - ALGEBRAIC[31])+ALGEBRAIC[26]+ALGEBRAIC[27])/CONSTANTS[8];
resid[3] = RATES[3] - (ALGEBRAIC[31] - ALGEBRAIC[8])/CONSTANTS[8];
resid[4] = RATES[4] - (((ALGEBRAIC[10] - ALGEBRAIC[9])+ 3.00000*ALGEBRAIC[16]) - ALGEBRAIC[17])/CONSTANTS[8];
resid[5] = RATES[5] - (((ALGEBRAIC[9] - ALGEBRAIC[10]) -  3.00000*ALGEBRAIC[16])+ALGEBRAIC[17])/CONSTANTS[8];
resid[6] = RATES[6] - (ALGEBRAIC[10] - ALGEBRAIC[11])/CONSTANTS[8];
resid[7] = RATES[7] - (ALGEBRAIC[11] - ALGEBRAIC[12])/CONSTANTS[8];
resid[8] = RATES[8] - (ALGEBRAIC[12] - ALGEBRAIC[13])/CONSTANTS[8];
resid[9] = RATES[9] - (ALGEBRAIC[13] - ALGEBRAIC[14])/CONSTANTS[8];
resid[10] = RATES[10] - (ALGEBRAIC[14] - ALGEBRAIC[15])/CONSTANTS[8];
resid[11] = RATES[11] - ((ALGEBRAIC[15] -  2.00000*ALGEBRAIC[16]) - ALGEBRAIC[17])/CONSTANTS[8];
resid[12] = RATES[12] - (((((ALGEBRAIC[8] - ALGEBRAIC[10]) - ALGEBRAIC[19])+ALGEBRAIC[24]) - ALGEBRAIC[26])+ALGEBRAIC[27])/CONSTANTS[8];
resid[13] = RATES[13] - ((ALGEBRAIC[8] - ALGEBRAIC[9])+ALGEBRAIC[19])/CONSTANTS[8];
resid[14] = RATES[14] - (ALGEBRAIC[20] - ALGEBRAIC[21])/CONSTANTS[8];
resid[15] = RATES[15] - (ALGEBRAIC[21] - ALGEBRAIC[22])/CONSTANTS[8];
resid[16] = RATES[16] - (ALGEBRAIC[28] - (ALGEBRAIC[20]+ALGEBRAIC[22]))/CONSTANTS[8];
resid[17] = RATES[17] - ((ALGEBRAIC[20]+ALGEBRAIC[22]) - ALGEBRAIC[28])/CONSTANTS[8];
resid[18] = RATES[18] - ((ALGEBRAIC[22] - ALGEBRAIC[23]) - ALGEBRAIC[25])/CONSTANTS[8];
resid[19] = RATES[19] - (ALGEBRAIC[23] - ALGEBRAIC[24])/CONSTANTS[8];
resid[20] = RATES[20] - ((ALGEBRAIC[25] - ALGEBRAIC[24]) - ALGEBRAIC[27])/CONSTANTS[8];
resid[21] = RATES[21] - (ALGEBRAIC[24] - ALGEBRAIC[26])/CONSTANTS[8];
resid[22] = RATES[22] - (ALGEBRAIC[26] - ALGEBRAIC[27])/CONSTANTS[8];
resid[23] = RATES[23] - (ALGEBRAIC[29] - ALGEBRAIC[1])/CONSTANTS[8];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = STATES[17]/STATES[16];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = ( CONSTANTS[8]*CONSTANTS[11]*(( - STATES[1]*(CONSTANTS[2] -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[12]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[12]*pow(STATES[0], 2.00000)), 1.0 / 2)))/( (1.00000 -  4.00000*CONSTANTS[12])*CONSTANTS[13])+( STATES[23]*(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[12]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[12]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[12]*pow(STATES[0], 2.00000)), 1.0 / 2)))/(2.00000 -  8.00000*CONSTANTS[12])))/( CONSTANTS[14]*CONSTANTS[15]*(1.00000+STATES[1]/CONSTANTS[16]+STATES[23]/CONSTANTS[15])*(1.00000+(CONSTANTS[2] -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[12]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[12]*pow(STATES[0], 2.00000)), 1.0 / 2))/( (1.00000 -  4.00000*CONSTANTS[12])*CONSTANTS[17])+(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[12]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[12]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[12]*pow(STATES[0], 2.00000)), 1.0 / 2))/( (2.00000 -  8.00000*CONSTANTS[12])*CONSTANTS[14])));
ALGEBRAIC[2] = ( (( CONSTANTS[8]*CONSTANTS[18])/CONSTANTS[19])*(STATES[1] - STATES[2]/CONSTANTS[20]))/(1.00000+STATES[1]/CONSTANTS[19]+STATES[2]/CONSTANTS[21]);
ALGEBRAIC[8] = ( (( CONSTANTS[8]*CONSTANTS[37]*STATES[3])/CONSTANTS[40])*(1.00000 - ( STATES[13]*STATES[12])/( STATES[3]*CONSTANTS[39])))/(1.00000+STATES[3]/CONSTANTS[40]+STATES[13]/CONSTANTS[41]+STATES[12]/CONSTANTS[42]+( STATES[3]*STATES[12])/( CONSTANTS[40]*CONSTANTS[43])+( STATES[13]*STATES[12])/( CONSTANTS[41]*CONSTANTS[42]));
ALGEBRAIC[9] = ( CONSTANTS[8]*CONSTANTS[44]*(( - CONSTANTS[6]*STATES[5])/CONSTANTS[45]+( STATES[4]*STATES[13])/(1.00000+CONSTANTS[46])))/( CONSTANTS[47]*CONSTANTS[48]*(1.00000+STATES[5]/CONSTANTS[49]+STATES[4]/CONSTANTS[48])*(1.00000+CONSTANTS[6]/CONSTANTS[50]+STATES[13]/( (1.00000+CONSTANTS[46])*CONSTANTS[47])));
ALGEBRAIC[10] = ( (( CONSTANTS[8]*CONSTANTS[1]*CONSTANTS[55]*STATES[12]*STATES[5])/( CONSTANTS[56]*CONSTANTS[57]))*(1.00000 - ( STATES[6]*STATES[4])/( STATES[12]*STATES[5]*CONSTANTS[58])))/( (1.00000+STATES[12]/CONSTANTS[56]+STATES[6]/CONSTANTS[52])*(1.00000+STATES[5]/CONSTANTS[57]+STATES[4]/CONSTANTS[53]));
ALGEBRAIC[11] = ( 1.00000*CONSTANTS[8]*CONSTANTS[59]*(( CONSTANTS[60]*STATES[6]*(CONSTANTS[2] -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[61]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[61]*pow(STATES[0], 2.00000)), 1.0 / 2)))/(1.00000 -  4.00000*CONSTANTS[61]) - ( (((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[61]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[61]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[61]*pow(STATES[0], 2.00000)), 1.0 / 2))*STATES[7])/(2.00000 -  8.00000*CONSTANTS[61])))/( CONSTANTS[62]*CONSTANTS[63]*(1.00000+(CONSTANTS[2] -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[61]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[61]*pow(STATES[0], 2.00000)), 1.0 / 2))/( (1.00000 -  4.00000*CONSTANTS[61])*CONSTANTS[64])+(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[61]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[61]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[61]*pow(STATES[0], 2.00000)), 1.0 / 2))/( (2.00000 -  8.00000*CONSTANTS[61])*CONSTANTS[62]))*(1.00000+STATES[6]/CONSTANTS[65]+STATES[7]/CONSTANTS[63]));
ALGEBRAIC[12] = ( (( CONSTANTS[8]*CONSTANTS[66])/CONSTANTS[67])*(STATES[7] - STATES[8]/CONSTANTS[68]))/(1.00000+STATES[7]/CONSTANTS[67]+STATES[8]/CONSTANTS[69]);
ALGEBRAIC[13] = ( (( CONSTANTS[8]*CONSTANTS[70])/CONSTANTS[71])*(STATES[8] - STATES[9]/CONSTANTS[72]))/(1.00000+STATES[8]/CONSTANTS[71]+STATES[9]/CONSTANTS[73]);
ALGEBRAIC[14] = ( (( CONSTANTS[8]*CONSTANTS[74])/( CONSTANTS[75]*CONSTANTS[76]))*(( STATES[9]*(CONSTANTS[2] -  pow((((pow(STATES[0], 2.00000) -  4.00000*CONSTANTS[77]*pow(STATES[0], 2.00000)) -  2.00000*STATES[0]*CONSTANTS[2])+ 8.00000*CONSTANTS[77]*STATES[0]*CONSTANTS[2]+pow(CONSTANTS[2], 2.00000)), 1.0 / 2)))/(1.00000 -  4.00000*CONSTANTS[77]) - (( STATES[10]*(((STATES[0] -  4.00000*CONSTANTS[77]*STATES[0]) - CONSTANTS[2])+ pow((((pow(STATES[0], 2.00000) -  4.00000*CONSTANTS[77]*pow(STATES[0], 2.00000)) -  2.00000*STATES[0]*CONSTANTS[2])+ 8.00000*CONSTANTS[77]*STATES[0]*CONSTANTS[2]+pow(CONSTANTS[2], 2.00000)), 1.0 / 2)))/(2.00000 -  8.00000*CONSTANTS[77]))/CONSTANTS[78]))/( (1.00000+STATES[9]/CONSTANTS[75]+STATES[10]/CONSTANTS[79])*(1.00000+((((STATES[0] -  4.00000*CONSTANTS[77]*STATES[0]) - CONSTANTS[2])+ pow((((pow(STATES[0], 2.00000) -  4.00000*CONSTANTS[77]*pow(STATES[0], 2.00000)) -  2.00000*STATES[0]*CONSTANTS[2])+ 8.00000*CONSTANTS[77]*STATES[0]*CONSTANTS[2]+pow(CONSTANTS[2], 2.00000)), 1.0 / 2))/(2.00000 -  8.00000*CONSTANTS[77]))/CONSTANTS[80]+((CONSTANTS[2] -  pow((((pow(STATES[0], 2.00000) -  4.00000*CONSTANTS[77]*pow(STATES[0], 2.00000)) -  2.00000*STATES[0]*CONSTANTS[2])+ 8.00000*CONSTANTS[77]*STATES[0]*CONSTANTS[2]+pow(CONSTANTS[2], 2.00000)), 1.0 / 2))/(1.00000 -  4.00000*CONSTANTS[77]))/CONSTANTS[76]));
ALGEBRAIC[15] = (( CONSTANTS[8]*CONSTANTS[81]*pow(STATES[10], CONSTANTS[82]))/pow(CONSTANTS[83], CONSTANTS[82]))/(1.00000+pow(STATES[10], CONSTANTS[82])/pow(CONSTANTS[83], CONSTANTS[82]));
ALGEBRAIC[16] =  CONSTANTS[8]*CONSTANTS[84]*STATES[11];
ALGEBRAIC[17] = ( (( CONSTANTS[8]*- CONSTANTS[85])/( CONSTANTS[86]*CONSTANTS[87]))*( STATES[5]*CONSTANTS[4] - ( STATES[4]*STATES[11])/CONSTANTS[88]))/(1.00000+STATES[5]/CONSTANTS[86]+( CONSTANTS[89]*CONSTANTS[4])/( CONSTANTS[86]*CONSTANTS[87])+( CONSTANTS[90]*STATES[11])/( CONSTANTS[91]*CONSTANTS[92])+STATES[4]/CONSTANTS[91]+( STATES[5]*CONSTANTS[4])/( CONSTANTS[86]*CONSTANTS[87])+( CONSTANTS[90]*STATES[5]*STATES[11])/( CONSTANTS[86]*CONSTANTS[91]*CONSTANTS[92])+( CONSTANTS[89]*CONSTANTS[4]*STATES[4])/( CONSTANTS[86]*CONSTANTS[87]*CONSTANTS[91])+( STATES[4]*STATES[11])/( CONSTANTS[91]*CONSTANTS[92])+( STATES[5]*CONSTANTS[4]*STATES[11])/( CONSTANTS[86]*CONSTANTS[87]*CONSTANTS[93])+( CONSTANTS[4]*STATES[4]*STATES[11])/( CONSTANTS[94]*CONSTANTS[91]*CONSTANTS[92]));
ALGEBRAIC[18] = ( CONSTANTS[8]*CONSTANTS[95]*(((STATES[0] -  4.00000*CONSTANTS[96]*STATES[0]) - CONSTANTS[2])+ pow((((pow(STATES[0], 2.00000) -  4.00000*CONSTANTS[96]*pow(STATES[0], 2.00000)) -  2.00000*STATES[0]*CONSTANTS[2])+ 8.00000*CONSTANTS[96]*STATES[0]*CONSTANTS[2]+pow(CONSTANTS[2], 2.00000)), 1.0 / 2)))/(2.00000 -  8.00000*CONSTANTS[96]);
ALGEBRAIC[19] = ( CONSTANTS[8]*CONSTANTS[0]*(( CONSTANTS[99]*STATES[12])/CONSTANTS[97] - ( CONSTANTS[100]*STATES[13])/CONSTANTS[98]))/(1.00000+STATES[12]/CONSTANTS[97]+STATES[13]/CONSTANTS[98]);
ALGEBRAIC[20] = (( CONSTANTS[8]*CONSTANTS[101]*STATES[1]*STATES[16])/( CONSTANTS[102]*CONSTANTS[103]))/( (1.00000+STATES[1]/CONSTANTS[102]+STATES[17]/CONSTANTS[104])*(1.00000+STATES[16]/CONSTANTS[103]));
ALGEBRAIC[21] = ( CONSTANTS[8]*CONSTANTS[105]*STATES[14])/(CONSTANTS[106]+STATES[14]);
ALGEBRAIC[22] = (( CONSTANTS[8]*CONSTANTS[107]*STATES[15]*STATES[16])/( CONSTANTS[108]*CONSTANTS[109]))/( (1.00000+STATES[15]/CONSTANTS[108]+STATES[17]/CONSTANTS[110])*(1.00000+STATES[16]/CONSTANTS[109]));
ALGEBRAIC[23] = ( CONSTANTS[8]*(( CONSTANTS[111]*STATES[18])/CONSTANTS[113] - ( CONSTANTS[112]*STATES[19])/CONSTANTS[114]))/(1.00000+STATES[18]/CONSTANTS[113]+STATES[19]/CONSTANTS[114]);
ALGEBRAIC[24] = ( CONSTANTS[8]*(( CONSTANTS[115]*STATES[19]*STATES[20])/( CONSTANTS[117]*CONSTANTS[118]) - ( CONSTANTS[116]*STATES[12]*STATES[21])/( CONSTANTS[119]*CONSTANTS[120])))/( (1.00000+STATES[19]/CONSTANTS[117]+STATES[12]/CONSTANTS[119])*(1.00000+STATES[20]/CONSTANTS[118]+STATES[21]/CONSTANTS[120]));
ALGEBRAIC[25] = ( CONSTANTS[8]*(( CONSTANTS[121]*STATES[18])/CONSTANTS[123] - ( CONSTANTS[122]*STATES[20])/CONSTANTS[124]))/(1.00000+STATES[18]/CONSTANTS[123]+STATES[20]/CONSTANTS[124]);
ALGEBRAIC[26] = ( CONSTANTS[8]*(( CONSTANTS[125]*STATES[12]*STATES[21])/( CONSTANTS[127]*CONSTANTS[128]) - ( CONSTANTS[126]*STATES[2]*STATES[22])/( CONSTANTS[129]*CONSTANTS[130])))/( (1.00000+STATES[12]/CONSTANTS[127]+STATES[2]/CONSTANTS[129])*(1.00000+STATES[21]/CONSTANTS[128]+STATES[22]/CONSTANTS[130]));
ALGEBRAIC[27] = ( CONSTANTS[8]*(( CONSTANTS[131]*STATES[22]*STATES[20])/( CONSTANTS[134]*CONSTANTS[133]) - ( CONSTANTS[132]*STATES[2]*STATES[12])/( CONSTANTS[135]*CONSTANTS[136])))/( (1.00000+STATES[20]/CONSTANTS[133]+STATES[12]/CONSTANTS[136])*(1.00000+STATES[22]/CONSTANTS[134]+STATES[2]/CONSTANTS[135]));
ALGEBRAIC[28] =  CONSTANTS[8]*CONSTANTS[137]*STATES[17];
ALGEBRAIC[29] = ( CONSTANTS[8]*CONSTANTS[138]*(CONSTANTS[7] - STATES[23]/CONSTANTS[139]))/( CONSTANTS[140]*(1.00000+CONSTANTS[7]/CONSTANTS[140]+STATES[23]/CONSTANTS[141]+( 0.910000*CONSTANTS[7]*STATES[23])/( CONSTANTS[141]*CONSTANTS[140])));
ALGEBRAIC[3] =  CONSTANTS[8]*CONSTANTS[22]*CONSTANTS[23]*STATES[2]*(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2));
ALGEBRAIC[4] = 1.00000+STATES[2]/CONSTANTS[25]+(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2))/( (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[26])+( CONSTANTS[22]*STATES[2]*(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2)))/( (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[26]*CONSTANTS[25]);
ALGEBRAIC[5] =  pow(1.00000+( CONSTANTS[28]*CONSTANTS[10])/CONSTANTS[29]+( CONSTANTS[30]*STATES[3])/CONSTANTS[31], 2.00000)*pow(1.00000+( 2.00000*CONSTANTS[32]*CONSTANTS[24]*pow(CONSTANTS[2] -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2), 2.00000))/( (- 1.00000+ 4.00000*CONSTANTS[24])*CONSTANTS[33]*(((CONSTANTS[2] - STATES[0])+ 4.00000*CONSTANTS[24]*STATES[0]) -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2))), 2.00000);
ALGEBRAIC[6] =  pow(1.00000+( CONSTANTS[34]*(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2)))/( (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[35]), 2.00000)*pow(1.00000+( CONSTANTS[36]*(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2)))/( (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[26]), 2.00000);
ALGEBRAIC[7] =  pow(1.00000+CONSTANTS[10]/CONSTANTS[29]+STATES[3]/CONSTANTS[31], 2.00000)*pow(1.00000+( 2.00000*CONSTANTS[24]*pow(CONSTANTS[2] -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2), 2.00000))/( (- 1.00000+ 4.00000*CONSTANTS[24])*CONSTANTS[33]*(((CONSTANTS[2] - STATES[0])+ 4.00000*CONSTANTS[24]*STATES[0]) -  pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2))), 2.00000)*pow(1.00000+(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2))/( (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[35]), 2.00000);
ALGEBRAIC[30] =  (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[26]*CONSTANTS[25]*(( CONSTANTS[27]*ALGEBRAIC[5]*ALGEBRAIC[6])/ALGEBRAIC[7]+pow(1.00000+STATES[2]/CONSTANTS[25]+(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2))/( (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[26])+( CONSTANTS[22]*STATES[2]*(((- CONSTANTS[2]+STATES[0]) -  4.00000*CONSTANTS[24]*STATES[0])+ pow((((pow(CONSTANTS[2], 2.00000) -  2.00000*CONSTANTS[2]*STATES[0])+ 8.00000*CONSTANTS[24]*CONSTANTS[2]*STATES[0]+pow(STATES[0], 2.00000)) -  4.00000*CONSTANTS[24]*pow(STATES[0], 2.00000)), 1.0 / 2)))/( (2.00000 -  8.00000*CONSTANTS[24])*CONSTANTS[26]*CONSTANTS[25]), 2.00000));
ALGEBRAIC[31] = ( ALGEBRAIC[3]*ALGEBRAIC[4])/ALGEBRAIC[30];
}
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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}