/* 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) { }