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 312 entries in the algebraic variable array.
There are a total of 25 entries in each of the rate and state variable arrays.
There are a total of 342 entries in the constant variable array.
*/
/*
* VOI is time in component environment (minute).
* ALGEBRAIC[0] is pH_cy in component environment (dimensionless).
* CONSTANTS[34] is addbuffer in component environment (dimensionless).
* CONSTANTS[35] is pHstat in component environment (dimensionless).
* CONSTANTS[2] is Par_97 in component parameters (dimensionless).
* CONSTANTS[3] is Par_98 in component parameters (dimensionless).
* STATES[0] is pH_calc in component differential_equations (dimensionless).
* CONSTANTS[0] is R in component global_parameters (kilojoule_per_kelvin_per_mole).
* CONSTANTS[1] is T1 in component global_parameters (kelvin).
* CONSTANTS[130] is T in component global_parameters (kelvin).
* CONSTANTS[125] is I in component global_parameters (molar).
* CONSTANTS[124] is Par_90 in component parameters (molar).
* CONSTANTS[129] is Par_94 in component parameters (kelvin).
* CONSTANTS[4] is Par_1 in component parameters (molar_per_minute).
* CONSTANTS[5] is Par_2 in component parameters (dimensionless).
* CONSTANTS[6] is Par_3 in component parameters (molar).
* CONSTANTS[36] is Par_4 in component parameters (molar).
* CONSTANTS[56] is Par_5 in component parameters (molar).
* CONSTANTS[7] is Par_6 in component parameters (molar).
* CONSTANTS[57] is Par_7 in component parameters (molar).
* CONSTANTS[8] is Par_8 in component parameters (molar).
* CONSTANTS[58] is Par_9 in component parameters (molar).
* CONSTANTS[59] is Par_10 in component parameters (molar).
* CONSTANTS[9] is Par_11 in component parameters (molar).
* CONSTANTS[60] is Par_12 in component parameters (molar).
* CONSTANTS[61] is Par_13 in component parameters (molar).
* CONSTANTS[10] is Par_14 in component parameters (molar).
* CONSTANTS[62] is Par_15 in component parameters (molar).
* CONSTANTS[11] is Par_16 in component parameters (molar).
* CONSTANTS[12] is Par_17 in component parameters (dimensionless).
* CONSTANTS[13] is Par_18 in component parameters (dimensionless).
* CONSTANTS[63] is Par_19 in component parameters (molar_per_minute).
* CONSTANTS[64] is Par_20 in component parameters (molar).
* CONSTANTS[65] is Par_21 in component parameters (molar).
* CONSTANTS[66] is Par_22 in component parameters (molar_per_minute).
* CONSTANTS[67] is Par_23 in component parameters (molar).
* CONSTANTS[68] is Par_24 in component parameters (molar).
* CONSTANTS[69] is Par_25 in component parameters (molar_per_minute).
* CONSTANTS[70] is Par_26 in component parameters (molar).
* CONSTANTS[71] is Par_27 in component parameters (molar).
* CONSTANTS[72] is Par_28 in component parameters (molar).
* CONSTANTS[73] is Par_29 in component parameters (molar).
* CONSTANTS[74] is Par_30 in component parameters (molar).
* CONSTANTS[14] is Par_31 in component parameters (molar).
* CONSTANTS[75] is Par_32 in component parameters (molar).
* CONSTANTS[15] is Par_33 in component parameters (molar).
* CONSTANTS[16] is Par_34 in component parameters (molar).
* CONSTANTS[17] is Par_35 in component parameters (molar).
* CONSTANTS[18] is Par_36 in component parameters (dimensionless).
* CONSTANTS[19] is Par_37 in component parameters (dimensionless).
* CONSTANTS[20] is Par_38 in component parameters (dimensionless).
* CONSTANTS[76] is Par_39 in component parameters (molar_per_minute).
* CONSTANTS[77] is Par_40 in component parameters (molar).
* CONSTANTS[78] is Par_41 in component parameters (molar).
* CONSTANTS[79] is Par_42 in component parameters (molar).
* CONSTANTS[80] is Par_43 in component parameters (molar_per_minute).
* CONSTANTS[81] is Par_44 in component parameters (molar).
* CONSTANTS[82] is Par_45 in component parameters (molar).
* CONSTANTS[83] is Par_46 in component parameters (molar_per_minute).
* CONSTANTS[84] is Par_47 in component parameters (molar).
* CONSTANTS[85] is Par_48 in component parameters (molar).
* CONSTANTS[86] is Par_49 in component parameters (molar).
* CONSTANTS[87] is Par_50 in component parameters (molar).
* CONSTANTS[88] is Par_51 in component parameters (molar_per_minute).
* CONSTANTS[89] is Par_52 in component parameters (molar).
* CONSTANTS[90] is Par_53 in component parameters (molar).
* CONSTANTS[91] is Par_54 in component parameters (molar).
* CONSTANTS[92] is Par_55 in component parameters (molar).
* CONSTANTS[93] is Par_56 in component parameters (molar).
* CONSTANTS[94] is Par_57 in component parameters (molar_per_minute).
* CONSTANTS[95] is Par_58 in component parameters (molar).
* CONSTANTS[96] is Par_59 in component parameters (molar).
* CONSTANTS[97] is Par_60 in component parameters (molar).
* CONSTANTS[98] is Par_61 in component parameters (molar).
* CONSTANTS[99] is Par_62 in component parameters (molar_per_minute).
* CONSTANTS[100] is Par_63 in component parameters (molar).
* CONSTANTS[101] is Par_64 in component parameters (molar).
* CONSTANTS[102] is Par_65 in component parameters (molar_per_minute).
* CONSTANTS[103] is Par_66 in component parameters (molar).
* CONSTANTS[104] is Par_67 in component parameters (molar).
* CONSTANTS[105] is Par_68 in component parameters (molar_per_minute).
* CONSTANTS[106] is Par_69 in component parameters (molar).
* CONSTANTS[107] is Par_70 in component parameters (molar).
* CONSTANTS[108] is Par_71 in component parameters (molar).
* CONSTANTS[109] is Par_72 in component parameters (molar).
* CONSTANTS[110] is Par_73 in component parameters (molar_per_minute).
* CONSTANTS[111] is Par_74 in component parameters (molar).
* CONSTANTS[112] is Par_75 in component parameters (molar).
* CONSTANTS[113] is Par_76 in component parameters (molar).
* CONSTANTS[114] is Par_77 in component parameters (molar).
* CONSTANTS[115] is Par_78 in component parameters (molar_per_minute).
* CONSTANTS[116] is Par_79 in component parameters (molar).
* CONSTANTS[117] is Par_80 in component parameters (molar).
* CONSTANTS[118] is Par_81 in component parameters (molar).
* CONSTANTS[21] is Par_82 in component parameters (molar).
* CONSTANTS[119] is Par_83 in component parameters (molar).
* CONSTANTS[120] is Par_84 in component parameters (molar_per_minute).
* CONSTANTS[121] is Par_85 in component parameters (molar).
* CONSTANTS[122] is Par_86 in component parameters (molar).
* CONSTANTS[123] is Par_87 in component parameters (molar).
* CONSTANTS[22] is Par_88 in component parameters (molar_per_minute).
* CONSTANTS[23] is Par_89 in component parameters (molar).
* CONSTANTS[126] is Par_91 in component parameters (molar).
* CONSTANTS[127] is Par_92 in component parameters (molar).
* CONSTANTS[128] is Par_93 in component parameters (molar).
* CONSTANTS[24] is Par_95 in component parameters (molar).
* CONSTANTS[131] is Par_96 in component parameters (molar).
* CONSTANTS[132] is Par_99 in component parameters (dimensionless).
* CONSTANTS[25] is Par_100 in component parameters (dimensionless).
* CONSTANTS[37] is mgT in component equilibrium_constants (molar).
* CONSTANTS[133] is k in component equilibrium_constants (molar).
* CONSTANTS[134] is c0 in component equilibrium_constants (molar).
* CONSTANTS[135] is RT2dadT in component correction_factors (kilojoule_half_liter_per_3_half_mole).
* CONSTANTS[136] is B in component correction_factors (per_half_molar).
* CONSTANTS[137] is Icorr in component correction_factors (kilojoule_per_mole).
* CONSTANTS[138] is I1 in component correction_factors (molar).
* CONSTANTS[139] is alphadebye in component correction_factors (per_half_molar).
* CONSTANTS[140] is IcorrpKa in component correction_factors (kilojoule_per_mole).
* CONSTANTS[141] is TcorrpKa in component correction_factors (mole_per_kilojoule).
* CONSTANTS[142] is RTalpha in component correction_factors (kilojoule_half_liter_per_3_half_mole).
* CONSTANTS[143] is IcorrdeltaGpof in component correction_factors (kilojoule_per_mole).
* CONSTANTS[144] is pKak_Pi in component correction_factors (dimensionless).
* CONSTANTS[145] is deltaH1o_Pi in component correction_factors (kilojoule_per_mole).
* CONSTANTS[146] is deltaHmgo_Pi in component correction_factors (kilojoule_per_mole).
* CONSTANTS[147] is deltaH1_Pi in component correction_factors (kilojoule_per_mole).
* CONSTANTS[148] is deltaHmg_Pi in component correction_factors (kilojoule_per_mole).
* CONSTANTS[149] is pKa1_Pi in component correction_factors (dimensionless).
* CONSTANTS[150] is pKamg_Pi in component correction_factors (dimensionless).
* ALGEBRAIC[1] is P_Pi in component correction_factors (dimensionless).
* ALGEBRAIC[2] is HPi2 in component correction_factors (dimensionless).
* ALGEBRAIC[3] is H2Pi1 in component correction_factors (dimensionless).
* ALGEBRAIC[4] is kPi in component correction_factors (dimensionless).
* ALGEBRAIC[5] is mgPi in component correction_factors (dimensionless).
* ALGEBRAIC[6] is Navg_Pi in component correction_factors (dimensionless).
* ALGEBRAIC[7] is dNavgPidH in component correction_factors (per_molar).
* ALGEBRAIC[8] is dNavgPidmg in component correction_factors (per_molar).
* ALGEBRAIC[9] is dmgPidmg in component correction_factors (per_molar).
* ALGEBRAIC[10] is dmgPidpH in component correction_factors (dimensionless).
* CONSTANTS[151] is NH_HPi2 in component correction_factors (dimensionless).
* CONSTANTS[152] is deltaGof_HPi2 in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[11] is deltaGpof_HPi2 in component correction_factors (kilojoule_per_mole).
* CONSTANTS[153] is deltaH1o_ATP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[154] is deltaHmgo_ATP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[155] is deltaHko_ATP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[156] is deltaH1_ATP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[157] is deltaHmg_ATP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[158] is deltaHk_ATP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[159] is pKa1_ATP in component correction_factors (dimensionless).
* CONSTANTS[160] is pKamg_ATP in component correction_factors (dimensionless).
* CONSTANTS[161] is pKak_ATP in component correction_factors (dimensionless).
* ALGEBRAIC[12] is P_ATP in component correction_factors (dimensionless).
* ALGEBRAIC[13] is ATP4 in component correction_factors (dimensionless).
* ALGEBRAIC[14] is HATP3 in component correction_factors (dimensionless).
* ALGEBRAIC[15] is mgATP2 in component correction_factors (dimensionless).
* ALGEBRAIC[16] is kATP in component correction_factors (dimensionless).
* ALGEBRAIC[17] is Navg_ATP in component correction_factors (dimensionless).
* ALGEBRAIC[18] is dNavgATPdH in component correction_factors (per_molar).
* ALGEBRAIC[19] is dNavgATPdmg in component correction_factors (per_molar).
* ALGEBRAIC[20] is dmgATP2dmg in component correction_factors (per_molar).
* ALGEBRAIC[21] is dmgATP2dpH in component correction_factors (dimensionless).
* CONSTANTS[162] is NH_ATP4 in component correction_factors (dimensionless).
* CONSTANTS[163] is deltaGof_ATP4 in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[22] is deltaGpof_ATP4 in component correction_factors (kilojoule_per_mole).
* CONSTANTS[164] is pKak_ADP in component correction_factors (dimensionless).
* CONSTANTS[165] is deltaH1o_ADP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[166] is deltaHmgo_ADP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[167] is deltaH1_ADP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[168] is deltaHmg_ADP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[169] is pKa1_ADP in component correction_factors (dimensionless).
* CONSTANTS[170] is pKamg_ADP in component correction_factors (dimensionless).
* ALGEBRAIC[23] is P_ADP in component correction_factors (dimensionless).
* ALGEBRAIC[24] is ADP3 in component correction_factors (dimensionless).
* ALGEBRAIC[25] is HADP2 in component correction_factors (dimensionless).
* ALGEBRAIC[26] is mgADP in component correction_factors (dimensionless).
* ALGEBRAIC[27] is kADP in component correction_factors (dimensionless).
* ALGEBRAIC[28] is Navg_ADP in component correction_factors (dimensionless).
* ALGEBRAIC[29] is dNavgADPdH in component correction_factors (per_molar).
* ALGEBRAIC[30] is dNavgADPdmg in component correction_factors (per_molar).
* ALGEBRAIC[31] is dmgADPdmg in component correction_factors (per_molar).
* ALGEBRAIC[32] is dmgADPdpH in component correction_factors (dimensionless).
* CONSTANTS[171] is NH_ADP3 in component correction_factors (dimensionless).
* CONSTANTS[172] is deltaGof_ADP3 in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[33] is deltaGpof_ADP3 in component correction_factors (kilojoule_per_mole).
* CONSTANTS[173] is deltaH1o_AMP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[174] is deltaHmgo_AMP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[175] is deltaH1_AMP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[176] is deltaHmg_AMP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[177] is pKa1_AMP in component correction_factors (dimensionless).
* CONSTANTS[178] is pKamg_AMP in component correction_factors (dimensionless).
* ALGEBRAIC[34] is P_AMP in component correction_factors (dimensionless).
* ALGEBRAIC[35] is AMP2 in component correction_factors (dimensionless).
* ALGEBRAIC[36] is HAMP1 in component correction_factors (dimensionless).
* ALGEBRAIC[37] is mgAMP in component correction_factors (dimensionless).
* ALGEBRAIC[38] is Navg_AMP in component correction_factors (dimensionless).
* ALGEBRAIC[39] is dNavgAMPdH in component correction_factors (per_molar).
* ALGEBRAIC[40] is dNavgAMPdmg in component correction_factors (per_molar).
* ALGEBRAIC[41] is dmgAMPdmg in component correction_factors (per_molar).
* ALGEBRAIC[42] is dmgAMPdpH in component correction_factors (dimensionless).
* CONSTANTS[179] is NH_AMP2 in component correction_factors (dimensionless).
* CONSTANTS[180] is deltaGof_AMP2 in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[43] is deltaGpof_AMP2 in component correction_factors (kilojoule_per_mole).
* CONSTANTS[181] is pKak_PCR in component correction_factors (dimensionless).
* CONSTANTS[182] is deltaH1o_PCR in component correction_factors (kilojoule_per_mole).
* CONSTANTS[183] is deltaHmgo_PCR in component correction_factors (kilojoule_per_mole).
* CONSTANTS[184] is deltaH1_PCR in component correction_factors (kilojoule_per_mole).
* CONSTANTS[185] is deltaHmg_PCR in component correction_factors (kilojoule_per_mole).
* CONSTANTS[186] is pKa1_PCR in component correction_factors (dimensionless).
* CONSTANTS[187] is pKamg_PCR in component correction_factors (dimensionless).
* ALGEBRAIC[44] is P_PCR in component correction_factors (dimensionless).
* ALGEBRAIC[45] is HPCR in component correction_factors (dimensionless).
* ALGEBRAIC[46] is H2PCR in component correction_factors (dimensionless).
* ALGEBRAIC[47] is kPCR in component correction_factors (dimensionless).
* ALGEBRAIC[48] is mgPCR in component correction_factors (dimensionless).
* ALGEBRAIC[49] is Navg_PCR in component correction_factors (dimensionless).
* ALGEBRAIC[50] is dNavgPCRdH in component correction_factors (per_molar).
* ALGEBRAIC[51] is dNavgPCRdmg in component correction_factors (per_molar).
* ALGEBRAIC[52] is dmgPCRdmg in component correction_factors (per_molar).
* ALGEBRAIC[53] is dmgPCRdpH in component correction_factors (dimensionless).
* CONSTANTS[188] is NH_HPCR in component correction_factors (dimensionless).
* CONSTANTS[189] is pKa1_CR in component correction_factors (dimensionless).
* ALGEBRAIC[54] is P_CR in component correction_factors (dimensionless).
* ALGEBRAIC[55] is HCR in component correction_factors (dimensionless).
* ALGEBRAIC[56] is H2CR in component correction_factors (dimensionless).
* ALGEBRAIC[57] is Navg_CR in component correction_factors (dimensionless).
* ALGEBRAIC[58] is dNavgCRdH in component correction_factors (per_molar).
* CONSTANTS[26] is dNavgCRdmg in component correction_factors (per_molar).
* CONSTANTS[190] is NH_HCR in component correction_factors (dimensionless).
* CONSTANTS[191] is deltaH1o_G1P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[192] is deltaHmgo_G1P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[193] is deltaH1_G1P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[194] is deltaHmg_G1P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[195] is pKa1_G1P in component correction_factors (dimensionless).
* CONSTANTS[196] is pKamg_G1P in component correction_factors (dimensionless).
* ALGEBRAIC[59] is P_G1P in component correction_factors (dimensionless).
* ALGEBRAIC[60] is UG1P in component correction_factors (dimensionless).
* ALGEBRAIC[61] is HG1P in component correction_factors (dimensionless).
* ALGEBRAIC[62] is mgG1P in component correction_factors (dimensionless).
* ALGEBRAIC[63] is Navg_G1P in component correction_factors (dimensionless).
* ALGEBRAIC[64] is dNavgG1PdH in component correction_factors (per_molar).
* ALGEBRAIC[65] is dNavgG1Pdmg in component correction_factors (per_molar).
* ALGEBRAIC[66] is dmgG1Pdmg in component correction_factors (per_molar).
* ALGEBRAIC[67] is dmgG1PdpH in component correction_factors (dimensionless).
* CONSTANTS[197] is NH_UG1P in component correction_factors (dimensionless).
* CONSTANTS[198] is deltaGof_UG1P in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[68] is deltaGpof_UG1P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[199] is pKa1_G6P in component correction_factors (dimensionless).
* ALGEBRAIC[69] is P_G6P in component correction_factors (dimensionless).
* ALGEBRAIC[70] is UG6P in component correction_factors (dimensionless).
* ALGEBRAIC[71] is HG6P in component correction_factors (dimensionless).
* ALGEBRAIC[72] is Navg_G6P in component correction_factors (dimensionless).
* ALGEBRAIC[73] is dNavgG6PdH in component correction_factors (per_molar).
* CONSTANTS[27] is dNavgG6Pdmg in component correction_factors (per_molar).
* CONSTANTS[200] is NH_UG6P in component correction_factors (dimensionless).
* CONSTANTS[201] is deltaGof_UG6P in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[74] is deltaGpof_UG6P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[202] is pKa1_F6P in component correction_factors (dimensionless).
* ALGEBRAIC[75] is P_F6P in component correction_factors (dimensionless).
* ALGEBRAIC[76] is UF6P in component correction_factors (dimensionless).
* ALGEBRAIC[77] is HF6P in component correction_factors (dimensionless).
* ALGEBRAIC[78] is Navg_F6P in component correction_factors (dimensionless).
* ALGEBRAIC[79] is dNavgF6PdH in component correction_factors (per_molar).
* CONSTANTS[28] is dNavgF6Pdmg in component correction_factors (per_molar).
* CONSTANTS[203] is NH_UF6P in component correction_factors (dimensionless).
* CONSTANTS[204] is deltaGof_UF6P in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[80] is deltaGpof_UF6P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[205] is pKa1_FDP in component correction_factors (dimensionless).
* CONSTANTS[206] is pKa2_FDP in component correction_factors (dimensionless).
* CONSTANTS[207] is pKamg_FDP in component correction_factors (dimensionless).
* ALGEBRAIC[81] is P_FDP in component correction_factors (dimensionless).
* ALGEBRAIC[82] is UFDP in component correction_factors (dimensionless).
* ALGEBRAIC[83] is HFDP in component correction_factors (dimensionless).
* ALGEBRAIC[84] is H2FDP in component correction_factors (dimensionless).
* ALGEBRAIC[85] is mgFDP in component correction_factors (dimensionless).
* ALGEBRAIC[86] is Navg_FDP in component correction_factors (dimensionless).
* ALGEBRAIC[87] is dNavgFDPdH in component correction_factors (per_molar).
* ALGEBRAIC[88] is dNavgFDPdmg in component correction_factors (per_molar).
* ALGEBRAIC[89] is dmgFDPdmg in component correction_factors (per_molar).
* ALGEBRAIC[90] is dmgFDPdpH in component correction_factors (dimensionless).
* CONSTANTS[208] is NH_UFDP in component correction_factors (dimensionless).
* CONSTANTS[209] is deltaGof_UFDP in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[91] is deltaGpof_UFDP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[210] is pKa1_GAP in component correction_factors (dimensionless).
* ALGEBRAIC[92] is P_GAP in component correction_factors (dimensionless).
* ALGEBRAIC[93] is UGAP in component correction_factors (dimensionless).
* ALGEBRAIC[94] is HGAP in component correction_factors (dimensionless).
* ALGEBRAIC[95] is Navg_GAP in component correction_factors (dimensionless).
* ALGEBRAIC[96] is dNavgGAPdH in component correction_factors (per_molar).
* CONSTANTS[29] is dNavgGAPdmg in component correction_factors (per_molar).
* CONSTANTS[211] is NH_UGAP in component correction_factors (dimensionless).
* CONSTANTS[212] is deltaGof_UGAP in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[97] is deltaGpof_UGAP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[213] is pKamg_G3P in component correction_factors (dimensionless).
* CONSTANTS[214] is deltaH1o_G3P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[215] is deltaH1_G3P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[216] is pKa1_G3P in component correction_factors (dimensionless).
* ALGEBRAIC[98] is P_G3P in component correction_factors (dimensionless).
* ALGEBRAIC[99] is UG3P in component correction_factors (dimensionless).
* ALGEBRAIC[100] is HG3P in component correction_factors (dimensionless).
* ALGEBRAIC[101] is mgG3P in component correction_factors (dimensionless).
* ALGEBRAIC[102] is Navg_G3P in component correction_factors (dimensionless).
* ALGEBRAIC[103] is dNavgG3PdH in component correction_factors (per_molar).
* ALGEBRAIC[104] is dNavgG3Pdmg in component correction_factors (per_molar).
* ALGEBRAIC[105] is dmgG3Pdmg in component correction_factors (per_molar).
* ALGEBRAIC[106] is dmgG3PdpH in component correction_factors (dimensionless).
* CONSTANTS[217] is NH_UG3P in component correction_factors (dimensionless).
* CONSTANTS[218] is deltaGof_UG3P in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[107] is deltaGpof_UG3P in component correction_factors (kilojoule_per_mole).
* CONSTANTS[219] is pKa1_DHAP in component correction_factors (dimensionless).
* CONSTANTS[220] is pKamg_DHAP in component correction_factors (dimensionless).
* ALGEBRAIC[108] is P_DHAP in component correction_factors (dimensionless).
* ALGEBRAIC[109] is UDHAP in component correction_factors (dimensionless).
* ALGEBRAIC[110] is HDHAP in component correction_factors (dimensionless).
* ALGEBRAIC[111] is mgDHAP in component correction_factors (dimensionless).
* ALGEBRAIC[112] is Navg_DHAP in component correction_factors (dimensionless).
* ALGEBRAIC[113] is dNavgDHAPdH in component correction_factors (per_molar).
* ALGEBRAIC[114] is dNavgDHAPdmg in component correction_factors (per_molar).
* ALGEBRAIC[115] is dmgDHAPdmg in component correction_factors (per_molar).
* ALGEBRAIC[116] is dmgDHAPdpH in component correction_factors (dimensionless).
* CONSTANTS[221] is NH_UDHAP in component correction_factors (dimensionless).
* CONSTANTS[222] is deltaGof_UDHAP in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[117] is deltaGpof_UDHAP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[223] is pKa1_13DPG in component correction_factors (dimensionless).
* ALGEBRAIC[118] is P_13DPG in component correction_factors (dimensionless).
* ALGEBRAIC[119] is U13DPG in component correction_factors (dimensionless).
* ALGEBRAIC[120] is H13DPG in component correction_factors (dimensionless).
* ALGEBRAIC[121] is Navg_13DPG in component correction_factors (dimensionless).
* ALGEBRAIC[122] is dNavg13DPGdH in component correction_factors (per_molar).
* CONSTANTS[30] is dNavg13DPGdmg in component correction_factors (per_molar).
* CONSTANTS[224] is NH_U13DPG in component correction_factors (dimensionless).
* CONSTANTS[225] is deltaGof_U13DPG in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[123] is deltaGpof_U13DPG in component correction_factors (kilojoule_per_mole).
* CONSTANTS[226] is pKa1_3PG in component correction_factors (dimensionless).
* ALGEBRAIC[124] is P_3PG in component correction_factors (dimensionless).
* ALGEBRAIC[125] is U3PG in component correction_factors (dimensionless).
* ALGEBRAIC[126] is H3PG in component correction_factors (dimensionless).
* ALGEBRAIC[127] is Navg_3PG in component correction_factors (dimensionless).
* ALGEBRAIC[128] is dNavg3PGdH in component correction_factors (per_molar).
* CONSTANTS[31] is dNavg3PGdmg in component correction_factors (per_molar).
* CONSTANTS[227] is NH_U3PG in component correction_factors (dimensionless).
* CONSTANTS[228] is deltaGof_U3PG in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[129] is deltaGpof_U3PG in component correction_factors (kilojoule_per_mole).
* CONSTANTS[229] is pKa1_2PG in component correction_factors (dimensionless).
* CONSTANTS[230] is pKamg_2PG in component correction_factors (dimensionless).
* CONSTANTS[231] is pKak_2PG in component correction_factors (dimensionless).
* ALGEBRAIC[130] is P_2PG in component correction_factors (dimensionless).
* ALGEBRAIC[131] is U2PG in component correction_factors (dimensionless).
* ALGEBRAIC[132] is H2PG in component correction_factors (dimensionless).
* ALGEBRAIC[133] is mg2PG in component correction_factors (dimensionless).
* ALGEBRAIC[134] is k2PG in component correction_factors (dimensionless).
* ALGEBRAIC[135] is Navg_2PG in component correction_factors (dimensionless).
* ALGEBRAIC[136] is dNavg2PGdH in component correction_factors (per_molar).
* ALGEBRAIC[137] is dNavg2PGdmg in component correction_factors (per_molar).
* ALGEBRAIC[138] is dmg2PGdmg in component correction_factors (per_molar).
* ALGEBRAIC[139] is dmg2PGdpH in component correction_factors (dimensionless).
* CONSTANTS[232] is NH_U2PG in component correction_factors (dimensionless).
* CONSTANTS[233] is deltaGof_U2PG in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[140] is deltaGpof_U2PG in component correction_factors (kilojoule_per_mole).
* CONSTANTS[234] is pKa1_PEP in component correction_factors (dimensionless).
* CONSTANTS[235] is pKamg_PEP in component correction_factors (dimensionless).
* CONSTANTS[236] is pKak_PEP in component correction_factors (dimensionless).
* ALGEBRAIC[141] is P_PEP in component correction_factors (dimensionless).
* ALGEBRAIC[142] is UPEP in component correction_factors (dimensionless).
* ALGEBRAIC[143] is HPEP in component correction_factors (dimensionless).
* ALGEBRAIC[144] is kPEP in component correction_factors (dimensionless).
* ALGEBRAIC[145] is mgPEP in component correction_factors (dimensionless).
* ALGEBRAIC[146] is Navg_PEP in component correction_factors (dimensionless).
* ALGEBRAIC[147] is dNavgPEPdH in component correction_factors (per_molar).
* ALGEBRAIC[148] is dNavgPEPdmg in component correction_factors (per_molar).
* ALGEBRAIC[149] is dmgPEPdmg in component correction_factors (per_molar).
* ALGEBRAIC[150] is dmgPEPdpH in component correction_factors (dimensionless).
* CONSTANTS[237] is NH_UPEP in component correction_factors (dimensionless).
* CONSTANTS[238] is deltaGof_UPEP in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[151] is deltaGpof_UPEP in component correction_factors (kilojoule_per_mole).
* CONSTANTS[239] is pKa1_PYR in component correction_factors (dimensionless).
* ALGEBRAIC[152] is P_PYR in component correction_factors (dimensionless).
* ALGEBRAIC[153] is UPYR in component correction_factors (dimensionless).
* ALGEBRAIC[154] is HPYR in component correction_factors (dimensionless).
* ALGEBRAIC[155] is Navg_PYR in component correction_factors (dimensionless).
* ALGEBRAIC[156] is dNavgPYRdH in component correction_factors (per_molar).
* CONSTANTS[32] is dNavgPYRdmg in component correction_factors (per_molar).
* CONSTANTS[240] is NH_UPYR in component correction_factors (dimensionless).
* CONSTANTS[241] is deltaGof_UPYR in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[157] is deltaGpof_UPYR in component correction_factors (kilojoule_per_mole).
* CONSTANTS[242] is pKamg_LAC in component correction_factors (dimensionless).
* CONSTANTS[243] is deltaH1o_LAC in component correction_factors (kilojoule_per_mole).
* CONSTANTS[244] is deltaH1_LAC in component correction_factors (kilojoule_per_mole).
* CONSTANTS[245] is pKa1_LAC in component correction_factors (dimensionless).
* ALGEBRAIC[158] is P_LAC in component correction_factors (dimensionless).
* ALGEBRAIC[159] is ULAC in component correction_factors (dimensionless).
* ALGEBRAIC[160] is HLAC in component correction_factors (dimensionless).
* ALGEBRAIC[161] is mgLAC in component correction_factors (dimensionless).
* ALGEBRAIC[162] is Navg_LAC in component correction_factors (dimensionless).
* ALGEBRAIC[163] is dNavgLACdH in component correction_factors (per_molar).
* ALGEBRAIC[164] is dNavgLACdmg in component correction_factors (per_molar).
* ALGEBRAIC[165] is dmgLACdmg in component correction_factors (per_molar).
* ALGEBRAIC[166] is dmgLACdpH in component correction_factors (dimensionless).
* CONSTANTS[246] is NH_ULAC in component correction_factors (dimensionless).
* CONSTANTS[247] is deltaGof_ULAC in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[167] is deltaGpof_ULAC in component correction_factors (kilojoule_per_mole).
* CONSTANTS[248] is dNH_GLY in component correction_factors (dimensionless).
* ALGEBRAIC[168] is deltaGpo_GLY in component correction_factors (kilojoule_per_mole).
* CONSTANTS[249] is NH_NAD in component correction_factors (dimensionless).
* CONSTANTS[250] is deltaGof_NAD in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[169] is deltaGpof_NAD in component correction_factors (kilojoule_per_mole).
* CONSTANTS[251] is NH_NADH in component correction_factors (dimensionless).
* CONSTANTS[252] is deltaGof_NADH in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[170] is deltaGpof_NADH in component correction_factors (kilojoule_per_mole).
* CONSTANTS[253] is NH_H2O in component correction_factors (dimensionless).
* CONSTANTS[254] is deltaGof_H2O in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[171] is deltaGpof_H2O in component correction_factors (kilojoule_per_mole).
* CONSTANTS[255] is NH_H in component correction_factors (dimensionless).
* CONSTANTS[256] is deltaGof_H in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[172] is deltaGpof_H in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[173] is deltaH_CK in component correction_factors (dimensionless).
* CONSTANTS[257] is Kref_CK in component correction_factors (dimensionless).
* CONSTANTS[258] is deltaHo_CKo in component correction_factors (kilojoule_per_mole).
* CONSTANTS[259] is deltaH1_CK in component correction_factors (kilojoule_per_mole).
* CONSTANTS[260] is Kref_CKI in component correction_factors (dimensionless).
* CONSTANTS[261] is Kref_CKT in component correction_factors (dimensionless).
* CONSTANTS[262] is deltaGpo_CK in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[174] is Kapp_CK in component correction_factors (dimensionless).
* ALGEBRAIC[175] is deltaH_ADK in component correction_factors (dimensionless).
* ALGEBRAIC[176] is deltaGpo_ADK in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[177] is Kapp_ADK in component correction_factors (dimensionless).
* ALGEBRAIC[178] is deltaH_GP in component correction_factors (dimensionless).
* ALGEBRAIC[179] is deltaGpo_GP in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[180] is Kapp_GP in component correction_factors (dimensionless).
* ALGEBRAIC[181] is deltaH_PGLM in component correction_factors (dimensionless).
* ALGEBRAIC[182] is deltaGpo_PGLM in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[183] is Kapp_PGLM in component correction_factors (dimensionless).
* ALGEBRAIC[184] is deltaH_PGI in component correction_factors (dimensionless).
* ALGEBRAIC[185] is deltaGpo_PGI in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[186] is Kapp_PGI in component correction_factors (dimensionless).
* ALGEBRAIC[187] is deltaH_PFK in component correction_factors (dimensionless).
* ALGEBRAIC[188] is deltaGpo_PFK in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[189] is Kapp_PFK in component correction_factors (dimensionless).
* ALGEBRAIC[190] is deltaH_ALD in component correction_factors (dimensionless).
* ALGEBRAIC[191] is deltaGpo_ALD in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[192] is Kapp_ALD in component correction_factors (molar).
* ALGEBRAIC[193] is deltaH_TPI in component correction_factors (dimensionless).
* ALGEBRAIC[194] is deltaGpo_TPI in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[195] is Kapp_TPI in component correction_factors (dimensionless).
* ALGEBRAIC[196] is deltaH_GAPDH in component correction_factors (dimensionless).
* ALGEBRAIC[197] is deltaGpo_GAPDH in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[198] is Kapp_GAPDH in component correction_factors (per_molar).
* ALGEBRAIC[199] is deltaH_G3PDH in component correction_factors (dimensionless).
* ALGEBRAIC[200] is deltaGpo_G3PDH in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[201] is Kapp_G3PDH in component correction_factors (dimensionless).
* ALGEBRAIC[202] is deltaH_PGK in component correction_factors (dimensionless).
* ALGEBRAIC[203] is deltaGpo_PGK in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[204] is Kapp_PGK in component correction_factors (dimensionless).
* ALGEBRAIC[205] is deltaH_PGM in component correction_factors (dimensionless).
* ALGEBRAIC[206] is deltaGpo_PGM in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[207] is Kapp_PGM in component correction_factors (dimensionless).
* ALGEBRAIC[208] is deltaH_ENOL in component correction_factors (dimensionless).
* ALGEBRAIC[209] is deltaGpo_ENOL in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[210] is Kapp_ENOL in component correction_factors (dimensionless).
* ALGEBRAIC[211] is deltaH_PK in component correction_factors (dimensionless).
* ALGEBRAIC[212] is deltaGpo_PK in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[213] is Kapp_PK in component correction_factors (dimensionless).
* ALGEBRAIC[214] is deltaH_LDH in component correction_factors (dimensionless).
* ALGEBRAIC[215] is deltaGpo_LDH in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[216] is Kapp_LDH in component correction_factors (dimensionless).
* ALGEBRAIC[217] is deltaH_ATPase in component correction_factors (dimensionless).
* ALGEBRAIC[218] is deltaGpo_ATPase in component correction_factors (kilojoule_per_mole).
* ALGEBRAIC[219] is Kapp_ATPase in component correction_factors (dimensionless).
* STATES[1] is Mg in component differential_equations (molar).
* CONSTANTS[263] is Vfgly in component glycogen_phosphorylase (molar_per_minute).
* CONSTANTS[38] is expno in component glycogen_phosphorylase (dimensionless).
* ALGEBRAIC[220] is fracA in component glycogen_phosphorylase (dimensionless).
* CONSTANTS[39] is KgpA_glyf in component glycogen_phosphorylase (molar).
* CONSTANTS[264] is KgpA_pi in component glycogen_phosphorylase (molar).
* CONSTANTS[265] is KgpA_igly in component glycogen_phosphorylase (molar).
* CONSTANTS[40] is KgpA_ipi in component glycogen_phosphorylase (molar).
* CONSTANTS[266] is KgpA_glyb in component glycogen_phosphorylase (molar).
* CONSTANTS[41] is KgpA_g1p in component glycogen_phosphorylase (molar).
* CONSTANTS[267] is KgpA_ig1p in component glycogen_phosphorylase (molar).
* ALGEBRAIC[221] is Dglya in component glycogen_phosphorylase (dimensionless).
* ALGEBRAIC[222] is pa in component glycogen_phosphorylase (dimensionless).
* ALGEBRAIC[223] is VbglyA in component glycogen_phosphorylase (molar_per_minute).
* ALGEBRAIC[224] is glyAF in component glycogen_phosphorylase (per_minute).
* ALGEBRAIC[225] is glyAR in component glycogen_phosphorylase (per_minute).
* ALGEBRAIC[226] is flux_GPa in component glycogen_phosphorylase (molar_per_minute).
* STATES[2] is G1P in component differential_equations (molar).
* STATES[3] is Pi in component differential_equations (molar).
* STATES[4] is Gly in component differential_equations (molar).
* ALGEBRAIC[227] is fracB in component glycogen_phosphorylase_B (dimensionless).
* CONSTANTS[268] is KgpB_pi in component glycogen_phosphorylase_B (molar).
* CONSTANTS[42] is KgpB_ipi in component glycogen_phosphorylase_B (molar).
* CONSTANTS[269] is KgpB_iglyf in component glycogen_phosphorylase_B (molar).
* CONSTANTS[270] is KgpB_g1p in component glycogen_phosphorylase_B (molar).
* CONSTANTS[43] is KgpB_ig1p in component glycogen_phosphorylase_B (molar).
* CONSTANTS[271] is KgpB_iglyb in component glycogen_phosphorylase_B (molar).
* CONSTANTS[44] is Kgp_amp in component glycogen_phosphorylase_B (molar).
* CONSTANTS[45] is interactioncoeff in component glycogen_phosphorylase_B (dimensionless).
* CONSTANTS[46] is nH in component glycogen_phosphorylase_B (dimensionless).
* ALGEBRAIC[228] is M in component glycogen_phosphorylase_B (dimensionless).
* ALGEBRAIC[229] is Dglyb in component glycogen_phosphorylase_B (dimensionless).
* ALGEBRAIC[230] is pb in component glycogen_phosphorylase_B (dimensionless).
* ALGEBRAIC[231] is VbglyB in component glycogen_phosphorylase_B (molar_per_minute).
* ALGEBRAIC[232] is glyBF in component glycogen_phosphorylase_B (per_minute).
* ALGEBRAIC[233] is glyBR in component glycogen_phosphorylase_B (per_minute).
* ALGEBRAIC[234] is flux_GPb in component glycogen_phosphorylase_B (molar_per_minute).
* STATES[5] is AMP in component differential_equations (molar).
* CONSTANTS[272] is Vffpglm in component PGLM (molar_per_minute).
* CONSTANTS[273] is Kpglm_g1p in component PGLM (molar).
* CONSTANTS[274] is Kpglm_g6p in component PGLM (molar).
* ALGEBRAIC[235] is Vfpglm in component PGLM (molar_per_minute).
* ALGEBRAIC[236] is Vbpglm in component PGLM (molar_per_minute).
* ALGEBRAIC[237] is v_PGLM in component PGLM (molar_per_minute).
* STATES[6] is G6P in component differential_equations (molar).
* CONSTANTS[275] is Vbbpgi in component PGI (molar_per_minute).
* CONSTANTS[276] is Kpgi_g6p in component PGI (molar).
* CONSTANTS[277] is Kpgi_f6p in component PGI (molar).
* ALGEBRAIC[238] is Vbpgi in component PGI (molar_per_minute).
* ALGEBRAIC[239] is Vfpgi in component PGI (molar_per_minute).
* ALGEBRAIC[240] is v_PGI in component PGI (molar_per_minute).
* STATES[7] is F6P in component differential_equations (molar).
* CONSTANTS[278] is Vffpfk in component PFK (molar_per_minute).
* CONSTANTS[279] is Kpfk_f6p in component PFK (molar).
* CONSTANTS[280] is Kpfk_f6pT in component PFK (molar).
* CONSTANTS[281] is Kpfk_atp in component PFK (molar).
* CONSTANTS[283] is Kpfk_atpT in component PFK (molar).
* CONSTANTS[282] is Kpfk_fbp in component PFK (molar).
* CONSTANTS[47] is Kpfk_fbpT in component PFK (molar).
* CONSTANTS[284] is Kpfk_adp in component PFK (molar).
* CONSTANTS[48] is Kpfk_adpT in component PFK (molar).
* CONSTANTS[49] is Kpfki in component PFK (molar).
* CONSTANTS[50] is Kmpfk in component PFK (molar).
* CONSTANTS[51] is d in component PFK (dimensionless).
* CONSTANTS[52] is e_ in component PFK (dimensionless).
* CONSTANTS[53] is Lo in component PFK (dimensionless).
* ALGEBRAIC[241] is Vfpfk in component PFK (molar_per_minute).
* ALGEBRAIC[242] is Vbpfk in component PFK (molar_per_minute).
* ALGEBRAIC[243] is L in component PFK (dimensionless).
* CONSTANTS[285] is alpha in component PFK (dimensionless).
* ALGEBRAIC[244] is Delta in component PFK (dimensionless).
* ALGEBRAIC[245] is Deltap in component PFK (dimensionless).
* ALGEBRAIC[246] is v_PFK in component PFK (molar_per_minute).
* STATES[8] is FBP in component differential_equations (molar).
* STATES[9] is ADP in component differential_equations (molar).
* STATES[10] is ATP in component differential_equations (molar).
* CONSTANTS[286] is Vffald in component ALD (molar_per_minute).
* CONSTANTS[287] is Kald_fbp in component ALD (molar).
* CONSTANTS[288] is Kald_dhap in component ALD (molar).
* CONSTANTS[289] is Kald_gap in component ALD (molar).
* ALGEBRAIC[247] is Vfald in component ALD (molar_per_minute).
* ALGEBRAIC[248] is Vbald in component ALD (molar_per_minute).
* ALGEBRAIC[249] is v_ALD in component ALD (molar_per_minute).
* STATES[11] is DHAP in component differential_equations (molar).
* STATES[12] is GAP in component differential_equations (molar).
* CONSTANTS[290] is Vfftpi in component TPI (molar_per_minute).
* CONSTANTS[291] is Ktpi_gap in component TPI (molar).
* CONSTANTS[292] is Ktpi_dhap in component TPI (molar).
* CONSTANTS[293] is Vftpi in component TPI (molar_per_minute).
* ALGEBRAIC[250] is Vbtpi in component TPI (molar_per_minute).
* ALGEBRAIC[251] is v_TPI in component TPI (molar_per_minute).
* CONSTANTS[294] is Vbbg3pdh in component G3PDH (molar_per_minute).
* CONSTANTS[295] is Kg3pdh_g3p in component G3PDH (molar).
* CONSTANTS[296] is Kg3pdh_nad in component G3PDH (molar).
* CONSTANTS[297] is Kg3pdh_dhap in component G3PDH (molar).
* CONSTANTS[298] is Kg3pdh_nadh in component G3PDH (molar).
* ALGEBRAIC[252] is Dg3pdh in component G3PDH (dimensionless).
* CONSTANTS[299] is Vbg3pdh in component G3PDH (molar_per_minute).
* ALGEBRAIC[253] is Vfg3pdh in component G3PDH (molar_per_minute).
* ALGEBRAIC[254] is v_G3PDH in component G3PDH (molar_per_minute).
* STATES[13] is G3P in component differential_equations (molar).
* STATES[14] is NAD in component differential_equations (molar).
* STATES[15] is NADH in component differential_equations (molar).
* CONSTANTS[300] is Vffgad in component GAPDH (molar_per_minute).
* CONSTANTS[301] is Kgapdh_gap in component GAPDH (molar).
* CONSTANTS[302] is Kgapdh_nad in component GAPDH (molar).
* CONSTANTS[303] is Kgapdh_pi in component GAPDH (molar).
* CONSTANTS[304] is Kgapdh_bpg in component GAPDH (molar).
* CONSTANTS[305] is Kgapdh_nadh in component GAPDH (molar).
* ALGEBRAIC[255] is Dgap in component GAPDH (dimensionless).
* ALGEBRAIC[256] is Vfgad in component GAPDH (molar_per_minute).
* ALGEBRAIC[257] is Vbgad in component GAPDH (molar_per_minute).
* ALGEBRAIC[258] is v_GAPDH in component GAPDH (molar_per_minute).
* STATES[16] is BPG in component differential_equations (molar).
* CONSTANTS[306] is Vbbpgk in component PGK (molar_per_minute).
* CONSTANTS[307] is Kpgk_bpg in component PGK (molar).
* CONSTANTS[308] is Kpgk_adp in component PGK (molar).
* CONSTANTS[309] is Kpgk_3pg in component PGK (molar).
* CONSTANTS[310] is Kpgk_atp in component PGK (molar).
* CONSTANTS[311] is Vbpgk in component PGK (molar_per_minute).
* ALGEBRAIC[259] is Vfpgk in component PGK (molar_per_minute).
* ALGEBRAIC[260] is D_PGK in component PGK (dimensionless).
* ALGEBRAIC[261] is v_PGK in component PGK (molar_per_minute).
* STATES[17] is P3G in component differential_equations (molar).
* CONSTANTS[312] is Vffpgm in component PGM (molar_per_minute).
* CONSTANTS[313] is Kpgm_3pg in component PGM (molar).
* CONSTANTS[314] is Kpgm_2pg in component PGM (molar).
* ALGEBRAIC[262] is Vfpgm in component PGM (molar_per_minute).
* ALGEBRAIC[263] is Vbpgm in component PGM (molar_per_minute).
* ALGEBRAIC[264] is v_PGM in component PGM (molar_per_minute).
* STATES[18] is P2G in component differential_equations (molar).
* CONSTANTS[315] is Vffen in component ENOL (molar_per_minute).
* CONSTANTS[316] is Ken_2pg in component ENOL (molar).
* CONSTANTS[317] is Ken_pep in component ENOL (molar).
* CONSTANTS[318] is Vfen in component ENOL (molar_per_minute).
* ALGEBRAIC[265] is Vben in component ENOL (molar_per_minute).
* ALGEBRAIC[266] is v_ENOL in component ENOL (molar_per_minute).
* STATES[19] is PEP in component differential_equations (molar).
* CONSTANTS[319] is Vffpk in component PK (molar_per_minute).
* CONSTANTS[320] is Kpk_pep in component PK (molar).
* CONSTANTS[321] is Kpk_adp in component PK (molar).
* CONSTANTS[322] is Kpk_pyr in component PK (molar).
* CONSTANTS[323] is Kpk_atp in component PK (molar).
* ALGEBRAIC[267] is Vfpk in component PK (molar_per_minute).
* ALGEBRAIC[268] is Vbpk in component PK (molar_per_minute).
* ALGEBRAIC[269] is v_PK in component PK (molar_per_minute).
* STATES[20] is PYR in component differential_equations (molar).
* CONSTANTS[324] is Vffldh in component LDH (molar_per_minute).
* CONSTANTS[325] is Kldh_pyr in component LDH (molar).
* CONSTANTS[326] is Kldh_nadh in component LDH (molar).
* CONSTANTS[327] is Kldh_lac in component LDH (molar).
* CONSTANTS[328] is Kldh_nad in component LDH (molar).
* ALGEBRAIC[270] is Vfldh in component LDH (molar_per_minute).
* ALGEBRAIC[271] is Vbldh in component LDH (molar_per_minute).
* ALGEBRAIC[272] is v_LDH in component LDH (molar_per_minute).
* STATES[21] is LAC in component differential_equations (molar).
* ALGEBRAIC[273] is VmaxATPase in component ATPase (molar_per_minute).
* CONSTANTS[54] is Katp_ATPase in component ATPase (molar).
* ALGEBRAIC[274] is ATPase in component ATPase (molar_per_minute).
* CONSTANTS[329] is VforCK in component creatine_kinase (molar_per_minute).
* CONSTANTS[330] is Kck_pcr in component creatine_kinase (molar).
* CONSTANTS[331] is Kck_iatp in component creatine_kinase (molar).
* CONSTANTS[332] is Kck_iadp in component creatine_kinase (molar).
* CONSTANTS[55] is Kck_ipcr in component creatine_kinase (molar).
* CONSTANTS[333] is Kck_cr in component creatine_kinase (molar).
* ALGEBRAIC[275] is VrevCK in component creatine_kinase (molar_per_minute).
* ALGEBRAIC[276] is CK in component creatine_kinase (molar_per_minute).
* STATES[22] is Cr in component differential_equations (molar).
* STATES[23] is PCr in component differential_equations (molar).
* CONSTANTS[334] is Vfadk in component adenylate_kinase (molar_per_minute).
* CONSTANTS[335] is Kadk_amp in component adenylate_kinase (molar).
* CONSTANTS[336] is Kadk_atp in component adenylate_kinase (molar).
* CONSTANTS[337] is Kadk_adp in component adenylate_kinase (molar).
* ALGEBRAIC[277] is Vbadk in component adenylate_kinase (molar_per_minute).
* ALGEBRAIC[278] is ADK in component adenylate_kinase (molar_per_minute).
* CONSTANTS[338] is carnosine in component buffer_capacity (molar).
* CONSTANTS[339] is tris in component buffer_capacity (molar).
* CONSTANTS[340] is acetate in component buffer_capacity (molar).
* ALGEBRAIC[279] is bufcapfixed in component buffer_capacity (molar).
* ALGEBRAIC[280] is bufcapmetab in component buffer_capacity (molar).
* ALGEBRAIC[281] is protons_consumed in component buffer_capacity (molar_per_minute).
* ALGEBRAIC[282] is CKprtflux in component buffer_capacity (molar_per_minute).
* ALGEBRAIC[283] is glycprtflux in component buffer_capacity (molar_per_minute).
* ALGEBRAIC[284] is pHODEterm1 in component buffer_capacity (per_minute).
* ALGEBRAIC[285] is pHODEterm2 in component buffer_capacity (per_molar).
* ALGEBRAIC[286] is denom_mgODE in component buffer_capacity (dimensionless).
* ALGEBRAIC[287] is RHSterm1_mgODE in component buffer_capacity (molar).
* ALGEBRAIC[288] is denomMgpHODE in component buffer_capacity (dimensionless).
* ALGEBRAIC[311] is RHSterm2_mgODE in component differential_equations (molar_per_minute).
* CONSTANTS[33] is fixmg in component differential_equations (dimensionless).
* CONSTANTS[341] is fixpH in component differential_equations (dimensionless).
* STATES[24] is protonload in component differential_equations (molar).
* ALGEBRAIC[293] is dATPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[294] is dADPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[295] is dAMPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[302] is dDHAPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[301] is dFBPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[298] is dG1Pdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[303] is dG3Pdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[310] is dLACdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[307] is dP2Gdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[289] is dPCrdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[308] is dPEPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[296] is dPidt in component differential_equations (molar_per_minute).
* ALGEBRAIC[290] is dCrdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[291] is dNADdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[292] is dNADHdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[297] is dGlydt in component differential_equations (molar_per_minute).
* ALGEBRAIC[299] is dG6Pdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[300] is dF6Pdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[304] is dGAPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[305] is dBGPdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[306] is dP3Gdt in component differential_equations (molar_per_minute).
* ALGEBRAIC[309] is dPYRdt in component differential_equations (molar_per_minute).
* RATES[23] is d/dt PCr in component differential_equations (molar).
* RATES[22] is d/dt Cr in component differential_equations (molar).
* RATES[14] is d/dt NAD in component differential_equations (molar).
* RATES[15] is d/dt NADH in component differential_equations (molar).
* RATES[10] is d/dt ATP in component differential_equations (molar).
* RATES[9] is d/dt ADP in component differential_equations (molar).
* RATES[5] is d/dt AMP in component differential_equations (molar).
* RATES[3] is d/dt Pi in component differential_equations (molar).
* RATES[4] is d/dt Gly in component differential_equations (molar).
* RATES[2] is d/dt G1P in component differential_equations (molar).
* RATES[6] is d/dt G6P in component differential_equations (molar).
* RATES[7] is d/dt F6P in component differential_equations (molar).
* RATES[8] is d/dt FBP in component differential_equations (molar).
* RATES[11] is d/dt DHAP in component differential_equations (molar).
* RATES[13] is d/dt G3P in component differential_equations (molar).
* RATES[12] is d/dt GAP in component differential_equations (molar).
* RATES[16] is d/dt BPG in component differential_equations (molar).
* RATES[17] is d/dt P3G in component differential_equations (molar).
* RATES[18] is d/dt P2G in component differential_equations (molar).
* RATES[19] is d/dt PEP in component differential_equations (molar).
* RATES[20] is d/dt PYR in component differential_equations (molar).
* RATES[21] is d/dt LAC in component differential_equations (molar).
* RATES[1] is d/dt Mg in component differential_equations (molar).
* RATES[0] is d/dt pH_calc in component differential_equations (dimensionless).
* RATES[24] is d/dt protonload in component differential_equations (molar).
* There are a total of 11 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 7.8;
CONSTANTS[0] = 8.314e-3;
CONSTANTS[1] = 298.15;
STATES[1] = 5.132658807e-4;
STATES[2] = 1e-9;
STATES[3] = 0.03;
STATES[4] = 0.04;
STATES[5] = 1e-9;
STATES[6] = 1e-9;
STATES[7] = 1e-9;
STATES[8] = 1e-9;
STATES[9] = 1e-9;
STATES[10] = 0.005;
STATES[11] = 1e-9;
STATES[12] = 1e-9;
STATES[13] = 1e-9;
STATES[14] = 0.0005;
STATES[15] = 1e-9;
STATES[16] = 1e-9;
STATES[17] = 1e-9;
STATES[18] = 1e-9;
STATES[19] = 1e-9;
STATES[20] = 1e-9;
STATES[21] = 1e-9;
STATES[22] = 0.029999999;
STATES[23] = 1e-9;
STATES[24] = 0;
CONSTANTS[2] = 1.00000;
CONSTANTS[3] = 7.40000;
CONSTANTS[4] = 0.0500000;
CONSTANTS[5] = 0.400000;
CONSTANTS[6] = 0.00170000;
CONSTANTS[7] = 0.00470000;
CONSTANTS[8] = 0.00270000;
CONSTANTS[9] = 0.00460000;
CONSTANTS[10] = 0.00740000;
CONSTANTS[11] = 0.00266055;
CONSTANTS[12] = 0.0200000;
CONSTANTS[13] = 1.75052;
CONSTANTS[14] = 0.00402000;
CONSTANTS[15] = 0.00270000;
CONSTANTS[16] = 0.000870000;
CONSTANTS[17] = 6.00000e-05;
CONSTANTS[18] = 0.0100000;
CONSTANTS[19] = 0.0100000;
CONSTANTS[20] = 13.0000;
CONSTANTS[21] = 0.00390000;
CONSTANTS[22] = 0.00000;
CONSTANTS[23] = 0.000100000;
CONSTANTS[24] = 0.00500000;
CONSTANTS[25] = 29.0000;
CONSTANTS[26] = 0.00000;
CONSTANTS[27] = 0.00000;
CONSTANTS[28] = 0.00000;
CONSTANTS[29] = 0.00000;
CONSTANTS[30] = 0.00000;
CONSTANTS[31] = 0.00000;
CONSTANTS[32] = 0.00000;
CONSTANTS[33] = 1.00000;
CONSTANTS[34] = CONSTANTS[2];
CONSTANTS[35] = CONSTANTS[3];
CONSTANTS[36] = 0.00400000;
CONSTANTS[37] = CONSTANTS[24];
CONSTANTS[38] = CONSTANTS[25];
CONSTANTS[39] = CONSTANTS[6];
CONSTANTS[40] = CONSTANTS[7];
CONSTANTS[41] = CONSTANTS[8];
CONSTANTS[42] = CONSTANTS[9];
CONSTANTS[43] = CONSTANTS[10];
CONSTANTS[44] = CONSTANTS[11];
CONSTANTS[45] = CONSTANTS[12];
CONSTANTS[46] = CONSTANTS[13];
CONSTANTS[47] = CONSTANTS[14];
CONSTANTS[48] = CONSTANTS[15];
CONSTANTS[49] = CONSTANTS[16];
CONSTANTS[50] = CONSTANTS[17];
CONSTANTS[51] = CONSTANTS[18];
CONSTANTS[52] = CONSTANTS[19];
CONSTANTS[53] = CONSTANTS[20];
CONSTANTS[54] = CONSTANTS[23];
CONSTANTS[55] = CONSTANTS[21];
CONSTANTS[56] = 0.00200000/1.50000;
CONSTANTS[57] = 0.000150000;
CONSTANTS[58] = 0.0101000;
CONSTANTS[59] = 0.000200000;
CONSTANTS[60] = 0.0150000;
CONSTANTS[61] = 0.00150000;
CONSTANTS[62] = 0.00440000;
CONSTANTS[63] = 0.480000;
CONSTANTS[64] = 6.30000e-05;
CONSTANTS[65] = 3.00000e-05;
CONSTANTS[66] = 0.880000;
CONSTANTS[67] = 0.000480000;
CONSTANTS[68] = 0.000119000;
CONSTANTS[69] = 0.0560000;
CONSTANTS[70] = 0.000180000;
CONSTANTS[71] = 0.0200000;
CONSTANTS[72] = 8.00000e-05;
CONSTANTS[73] = 0.000250000;
CONSTANTS[74] = 0.00402000;
CONSTANTS[75] = 0.00270000;
CONSTANTS[76] = 0.0106591;
CONSTANTS[77] = 5.00000e-05;
CONSTANTS[78] = 0.00200000;
CONSTANTS[79] = 0.00100000;
CONSTANTS[80] = 12.0000;
CONSTANTS[81] = 0.000320000;
CONSTANTS[82] = 0.000610000;
CONSTANTS[83] = 0.0825000;
CONSTANTS[84] = 0.000180000;
CONSTANTS[85] = 1.20000e-05;
CONSTANTS[86] = 0.000220000;
CONSTANTS[87] = 8.00000e-06;
CONSTANTS[88] = 1.26500;
CONSTANTS[89] = 2.50000e-06;
CONSTANTS[90] = 9.00000e-05;
CONSTANTS[91] = 0.000290000;
CONSTANTS[92] = 8.00000e-07;
CONSTANTS[93] = 3.30000e-06;
CONSTANTS[94] = 1.12000;
CONSTANTS[95] = 0.00200000;
CONSTANTS[96] = 8.00000e-06;
CONSTANTS[97] = 0.00120000;
CONSTANTS[98] = 0.000350000;
CONSTANTS[99] = 1.12000;
CONSTANTS[100] = 0.000200000;
CONSTANTS[101] = 1.40000e-05;
CONSTANTS[102] = 0.192000;
CONSTANTS[103] = 0.000100000;
CONSTANTS[104] = 0.000370000;
CONSTANTS[105] = 1.44000;
CONSTANTS[106] = 8.00000e-05;
CONSTANTS[107] = 0.000300000;
CONSTANTS[108] = 0.00705000;
CONSTANTS[109] = 0.00113000;
CONSTANTS[110] = 1.92000;
CONSTANTS[111] = 0.000335000;
CONSTANTS[112] = 2.00000e-06;
CONSTANTS[113] = 0.0170000;
CONSTANTS[114] = 0.000849000;
CONSTANTS[115] = 0.500000;
CONSTANTS[116] = 0.00111000;
CONSTANTS[117] = 0.00350000;
CONSTANTS[118] = 0.000135000;
CONSTANTS[119] = 0.00380000;
CONSTANTS[120] = 0.880000;
CONSTANTS[121] = 0.000320000;
CONSTANTS[122] = 0.000270000;
CONSTANTS[123] = 0.000350000;
CONSTANTS[124] = 0.100000;
CONSTANTS[125] = CONSTANTS[124];
CONSTANTS[126] = 0.0150000;
CONSTANTS[127] = 0.0250000;
CONSTANTS[128] = 0.0100000;
CONSTANTS[129] = 303.150;
CONSTANTS[130] = CONSTANTS[129];
CONSTANTS[131] = 0.0800000;
CONSTANTS[132] = 1.00000;
CONSTANTS[133] = CONSTANTS[131];
CONSTANTS[134] = 1.00000;
CONSTANTS[135] = 1.47750;
CONSTANTS[136] = 1.60000;
CONSTANTS[137] = ( CONSTANTS[135]* pow(CONSTANTS[125], 1.0 / 2))/(1.00000+ CONSTANTS[136]* pow(CONSTANTS[125], 1.0 / 2));
CONSTANTS[138] = 0.100000;
CONSTANTS[139] = 1.17582;
CONSTANTS[140] = ( 1.00000*CONSTANTS[139]*( pow(CONSTANTS[138], 1.0 / 2)/(1.00000+ CONSTANTS[136]* pow(CONSTANTS[138], 1.0 / 2)) - pow(CONSTANTS[125], 1.0 / 2)/(1.00000+ CONSTANTS[136]* pow(CONSTANTS[125], 1.0 / 2))))/log(10.0000);
CONSTANTS[141] = (1.00000/CONSTANTS[130] - 1.00000/CONSTANTS[1])/( log(10.0000)*CONSTANTS[0]);
CONSTANTS[142] = 2.91482;
CONSTANTS[143] = ( CONSTANTS[142]* pow(CONSTANTS[125], 1.0 / 2))/(1.00000+ CONSTANTS[136]* pow(CONSTANTS[125], 1.0 / 2));
CONSTANTS[144] = 0.500000;
CONSTANTS[145] = 3.00000;
CONSTANTS[146] = - 2.90000;
CONSTANTS[147] = CONSTANTS[145]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000));
CONSTANTS[148] = CONSTANTS[146]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(0.00000, 2.00000));
CONSTANTS[149] = 6.75000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[147];
CONSTANTS[150] = 1.65000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[148];
CONSTANTS[151] = 1.00000;
CONSTANTS[152] = - 1096.10;
CONSTANTS[153] = - 5.00000;
CONSTANTS[154] = - 18.0000;
CONSTANTS[155] = - 1.00000;
CONSTANTS[156] = CONSTANTS[153]+ CONSTANTS[137]*((pow(4.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(3.00000, 2.00000));
CONSTANTS[157] = CONSTANTS[154]+ CONSTANTS[137]*((pow(4.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(2.00000, 2.00000));
CONSTANTS[158] = CONSTANTS[155]+ CONSTANTS[137]*((pow(4.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(3.00000, 2.00000));
CONSTANTS[159] = 6.48000+ (CONSTANTS[140]/1.00000)*((pow(4.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(3.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[156];
CONSTANTS[160] = 4.19000+ (CONSTANTS[140]/1.00000)*((pow(4.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(2.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[157];
CONSTANTS[161] = 1.17000+ (CONSTANTS[140]/1.00000)*((pow(4.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(3.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[158];
CONSTANTS[162] = 12.0000;
CONSTANTS[163] = - 2768.10;
CONSTANTS[164] = 1.00000;
CONSTANTS[165] = - 3.00000;
CONSTANTS[166] = - 15.0000;
CONSTANTS[167] = CONSTANTS[165]+ CONSTANTS[137]*((pow(3.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(2.00000, 2.00000));
CONSTANTS[168] = CONSTANTS[166]+ CONSTANTS[137]*((pow(3.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(1.00000, 2.00000));
CONSTANTS[169] = 6.38000+ (CONSTANTS[140]/1.00000)*((pow(3.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(2.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[167];
CONSTANTS[170] = 3.25000+ (CONSTANTS[140]/1.00000)*((pow(3.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(1.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[168];
CONSTANTS[171] = 12.0000;
CONSTANTS[172] = - 1906.13;
CONSTANTS[173] = - 3.00000;
CONSTANTS[174] = - 7.50000;
CONSTANTS[175] = CONSTANTS[173]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000));
CONSTANTS[176] = CONSTANTS[174]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(0.00000, 2.00000));
CONSTANTS[177] = 6.29000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[175];
CONSTANTS[178] = 1.92000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(0.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[176];
CONSTANTS[179] = 12.0000;
CONSTANTS[180] = - 1040.45;
CONSTANTS[181] = 0.310000;
CONSTANTS[182] = 2.66000;
CONSTANTS[183] = 8.19000;
CONSTANTS[184] = CONSTANTS[182]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000));
CONSTANTS[185] = CONSTANTS[183]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(0.00000, 2.00000));
CONSTANTS[186] = 4.50000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[184];
CONSTANTS[187] = 1.60000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(0.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[185];
CONSTANTS[188] = 8.00000;
CONSTANTS[189] = 2.30000;
CONSTANTS[190] = 9.00000;
CONSTANTS[191] = - 1.70000;
CONSTANTS[192] = - 12.0000;
CONSTANTS[193] = CONSTANTS[191]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000));
CONSTANTS[194] = CONSTANTS[192]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(0.00000, 2.00000));
CONSTANTS[195] = 6.09000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[193];
CONSTANTS[196] = 2.48000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(2.00000, 2.00000)) - pow(0.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[194];
CONSTANTS[197] = 11.0000;
CONSTANTS[198] = - 1756.87;
CONSTANTS[199] = 6.11000;
CONSTANTS[200] = 11.0000;
CONSTANTS[201] = - 1763.94;
CONSTANTS[202] = 5.89000;
CONSTANTS[203] = 11.0000;
CONSTANTS[204] = - 1760.80;
CONSTANTS[205] = 6.40000;
CONSTANTS[206] = 5.92000;
CONSTANTS[207] = 2.70000;
CONSTANTS[208] = 10.0000;
CONSTANTS[209] = - 2601.40;
CONSTANTS[210] = 6.45000;
CONSTANTS[211] = 5.00000;
CONSTANTS[212] = - 1288.60;
CONSTANTS[213] = 1.63000;
CONSTANTS[214] = - 3.10000;
CONSTANTS[215] = CONSTANTS[214]+ CONSTANTS[137]*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000));
CONSTANTS[216] = 6.22000+ (CONSTANTS[140]/1.00000)*((pow(2.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(1.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[215];
CONSTANTS[217] = 7.00000;
CONSTANTS[218] = - 1339.25;
CONSTANTS[219] = 5.90000;
CONSTANTS[220] = 1.57000;
CONSTANTS[221] = 5.00000;
CONSTANTS[222] = - 1296.26;
CONSTANTS[223] = 7.50000;
CONSTANTS[224] = 4.00000;
CONSTANTS[225] = - 2356.14;
CONSTANTS[226] = 6.21000;
CONSTANTS[227] = 4.00000;
CONSTANTS[228] = - 1502.54;
CONSTANTS[229] = 7.00000;
CONSTANTS[230] = 2.45000;
CONSTANTS[231] = 1.18000;
CONSTANTS[232] = 4.00000;
CONSTANTS[233] = - 1496.38;
CONSTANTS[234] = 6.35000;
CONSTANTS[235] = 2.26000;
CONSTANTS[236] = 1.08000;
CONSTANTS[237] = 2.00000;
CONSTANTS[238] = - 1263.65;
CONSTANTS[239] = 2.49000;
CONSTANTS[240] = 3.00000;
CONSTANTS[241] = - 472.270;
CONSTANTS[242] = 0.980000;
CONSTANTS[243] = - 0.330000;
CONSTANTS[244] = CONSTANTS[243]+ CONSTANTS[137]*((pow(1.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(0.00000, 2.00000));
CONSTANTS[245] = 3.67000+ (CONSTANTS[140]/1.00000)*((pow(1.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(0.00000, 2.00000))+ CONSTANTS[141]*CONSTANTS[244];
CONSTANTS[246] = 5.00000;
CONSTANTS[247] = - 516.720;
CONSTANTS[248] = - 10.0000;
CONSTANTS[249] = 26.0000;
CONSTANTS[250] = 0.00000;
CONSTANTS[251] = 27.0000;
CONSTANTS[252] = 22.6500;
CONSTANTS[253] = 2.00000;
CONSTANTS[254] = - 237.190;
CONSTANTS[255] = 1.00000;
CONSTANTS[256] = 0.00000;
CONSTANTS[257] = 2.58000e+08;
CONSTANTS[258] = - 17.5500;
CONSTANTS[259] = CONSTANTS[258]+ CONSTANTS[137]*(((pow(2.00000, 2.00000)+pow(3.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(4.00000, 2.00000)) - pow(0.00000, 2.00000));
CONSTANTS[260] = exp(log(CONSTANTS[257])+( CONSTANTS[139]* pow(CONSTANTS[125], 1.0 / 2)*(((pow(2.00000, 2.00000)+pow(3.00000, 2.00000)+pow(1.00000, 2.00000)) - pow(4.00000, 2.00000)) - pow(0.00000, 2.00000)))/(1.00000+ CONSTANTS[136]* pow(CONSTANTS[125], 1.0 / 2)));
CONSTANTS[261] = pow(10.0000, arbitrary_log(CONSTANTS[260], 10) - CONSTANTS[141]*CONSTANTS[259]);
CONSTANTS[262] = - CONSTANTS[0]*CONSTANTS[130]*log(CONSTANTS[261]);
CONSTANTS[263] = CONSTANTS[4];
CONSTANTS[264] = CONSTANTS[36];
CONSTANTS[265] = CONSTANTS[56];
CONSTANTS[266] = CONSTANTS[57];
CONSTANTS[267] = CONSTANTS[58];
CONSTANTS[268] = CONSTANTS[59];
CONSTANTS[269] = CONSTANTS[60];
CONSTANTS[270] = CONSTANTS[61];
CONSTANTS[271] = CONSTANTS[62];
CONSTANTS[272] = CONSTANTS[63];
CONSTANTS[273] = CONSTANTS[64];
CONSTANTS[274] = CONSTANTS[65];
CONSTANTS[275] = CONSTANTS[66];
CONSTANTS[276] = CONSTANTS[67];
CONSTANTS[277] = CONSTANTS[68];
CONSTANTS[278] = CONSTANTS[69];
CONSTANTS[279] = CONSTANTS[70];
CONSTANTS[280] = CONSTANTS[71];
CONSTANTS[281] = CONSTANTS[72];
CONSTANTS[282] = CONSTANTS[74];
CONSTANTS[283] = CONSTANTS[73];
CONSTANTS[284] = CONSTANTS[75];
CONSTANTS[285] = ( CONSTANTS[279]*CONSTANTS[281])/( CONSTANTS[280]*CONSTANTS[283]);
CONSTANTS[286] = CONSTANTS[76];
CONSTANTS[287] = CONSTANTS[77];
CONSTANTS[288] = CONSTANTS[78];
CONSTANTS[289] = CONSTANTS[79];
CONSTANTS[290] = CONSTANTS[80];
CONSTANTS[291] = CONSTANTS[81];
CONSTANTS[292] = CONSTANTS[82];
CONSTANTS[293] = CONSTANTS[290];
CONSTANTS[294] = CONSTANTS[83];
CONSTANTS[295] = CONSTANTS[84];
CONSTANTS[296] = CONSTANTS[85];
CONSTANTS[297] = CONSTANTS[86];
CONSTANTS[298] = CONSTANTS[87];
CONSTANTS[299] = CONSTANTS[294];
CONSTANTS[300] = CONSTANTS[88];
CONSTANTS[301] = CONSTANTS[89];
CONSTANTS[302] = CONSTANTS[90];
CONSTANTS[303] = CONSTANTS[91];
CONSTANTS[304] = CONSTANTS[92];
CONSTANTS[305] = CONSTANTS[93];
CONSTANTS[306] = CONSTANTS[94];
CONSTANTS[307] = CONSTANTS[95];
CONSTANTS[308] = CONSTANTS[96];
CONSTANTS[309] = CONSTANTS[97];
CONSTANTS[310] = CONSTANTS[98];
CONSTANTS[311] = CONSTANTS[306];
CONSTANTS[312] = CONSTANTS[99];
CONSTANTS[313] = CONSTANTS[100];
CONSTANTS[314] = CONSTANTS[101];
CONSTANTS[315] = CONSTANTS[102];
CONSTANTS[316] = CONSTANTS[103];
CONSTANTS[317] = CONSTANTS[104];
CONSTANTS[318] = CONSTANTS[315];
CONSTANTS[319] = CONSTANTS[105];
CONSTANTS[320] = CONSTANTS[106];
CONSTANTS[321] = CONSTANTS[107];
CONSTANTS[322] = CONSTANTS[108];
CONSTANTS[323] = CONSTANTS[109];
CONSTANTS[324] = CONSTANTS[110];
CONSTANTS[325] = CONSTANTS[111];
CONSTANTS[326] = CONSTANTS[112];
CONSTANTS[327] = CONSTANTS[113];
CONSTANTS[328] = CONSTANTS[114];
CONSTANTS[329] = CONSTANTS[115];
CONSTANTS[330] = CONSTANTS[116];
CONSTANTS[331] = CONSTANTS[117];
CONSTANTS[332] = CONSTANTS[118];
CONSTANTS[333] = CONSTANTS[119];
CONSTANTS[334] = CONSTANTS[120];
CONSTANTS[335] = CONSTANTS[121];
CONSTANTS[336] = CONSTANTS[122];
CONSTANTS[337] = CONSTANTS[123];
CONSTANTS[338] = CONSTANTS[127];
CONSTANTS[339] = CONSTANTS[126];
CONSTANTS[340] = CONSTANTS[128];
CONSTANTS[341] = CONSTANTS[132];
RATES[23] = 0.1001;
RATES[22] = 0.1001;
RATES[14] = 0.1001;
RATES[15] = 0.1001;
RATES[10] = 0.1001;
RATES[9] = 0.1001;
RATES[5] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[2] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[11] = 0.1001;
RATES[13] = 0.1001;
RATES[12] = 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[1] = 0.1001;
RATES[0] = 0.1001;
RATES[24] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[23] - ALGEBRAIC[289];
resid[1] = RATES[22] - ALGEBRAIC[290];
resid[2] = RATES[14] - ALGEBRAIC[291];
resid[3] = RATES[15] - ALGEBRAIC[292];
resid[4] = RATES[10] - ALGEBRAIC[293];
resid[5] = RATES[9] - ALGEBRAIC[294];
resid[6] = RATES[5] - ALGEBRAIC[295];
resid[7] = RATES[3] - ALGEBRAIC[296];
resid[8] = RATES[4] - ALGEBRAIC[297];
resid[9] = RATES[2] - ALGEBRAIC[298];
resid[10] = RATES[6] - ALGEBRAIC[299];
resid[11] = RATES[7] - ALGEBRAIC[300];
resid[12] = RATES[8] - ALGEBRAIC[301];
resid[13] = RATES[11] - ALGEBRAIC[302];
resid[14] = RATES[13] - ALGEBRAIC[303];
resid[15] = RATES[12] - ALGEBRAIC[304];
resid[16] = RATES[16] - ALGEBRAIC[305];
resid[17] = RATES[17] - ALGEBRAIC[306];
resid[18] = RATES[18] - ALGEBRAIC[307];
resid[19] = RATES[19] - ALGEBRAIC[308];
resid[20] = RATES[20] - ALGEBRAIC[309];
resid[21] = RATES[21] - ALGEBRAIC[310];
resid[22] = RATES[1] - (CONDVAR[10]<=0.00000 ? ( CONSTANTS[33]*(ALGEBRAIC[311]+ ALGEBRAIC[287]*ALGEBRAIC[284]))/(1.00000 - ALGEBRAIC[287]*ALGEBRAIC[285]) : CONSTANTS[33]*ALGEBRAIC[311]);
resid[23] = RATES[0] - ( CONSTANTS[341]*(ALGEBRAIC[284]+ ALGEBRAIC[311]*ALGEBRAIC[285]))/(1.00000 - ALGEBRAIC[287]*ALGEBRAIC[285]);
resid[24] = RATES[24] - - ALGEBRAIC[281];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[4] = (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[144])*ALGEBRAIC[2];
ALGEBRAIC[47] = (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[181])*ALGEBRAIC[45];
ALGEBRAIC[134] = (( ALGEBRAIC[131]*CONSTANTS[133])/CONSTANTS[134])*pow(10.0000, CONSTANTS[231]);
ALGEBRAIC[144] = (( ALGEBRAIC[142]*CONSTANTS[133])/CONSTANTS[134])*pow(10.0000, CONSTANTS[236]);
ALGEBRAIC[217] = ((ALGEBRAIC[28]+ALGEBRAIC[6]) - ALGEBRAIC[17])+(((CONSTANTS[171]+CONSTANTS[151]) - CONSTANTS[162]) - CONSTANTS[253]);
ALGEBRAIC[218] = ((ALGEBRAIC[33]+ALGEBRAIC[11]+ALGEBRAIC[172]) - ALGEBRAIC[171]) - ALGEBRAIC[22];
ALGEBRAIC[219] = ( exp(- ALGEBRAIC[218]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[23]*ALGEBRAIC[1])/( ALGEBRAIC[12]*pow(10.0000, - ALGEBRAIC[0]));
ALGEBRAIC[282] = ALGEBRAIC[173]*- ALGEBRAIC[276];
ALGEBRAIC[283] = ALGEBRAIC[181]*ALGEBRAIC[237]+ ALGEBRAIC[178]*(ALGEBRAIC[226]+ALGEBRAIC[234])+ ALGEBRAIC[184]*ALGEBRAIC[240]+ ALGEBRAIC[187]*ALGEBRAIC[246]+ ALGEBRAIC[190]*ALGEBRAIC[249]+ ALGEBRAIC[193]*ALGEBRAIC[251]+ ALGEBRAIC[196]*ALGEBRAIC[258]+ ALGEBRAIC[202]*ALGEBRAIC[261]+ ALGEBRAIC[205]*ALGEBRAIC[264]+ ALGEBRAIC[208]*ALGEBRAIC[266]+ ALGEBRAIC[211]*ALGEBRAIC[269]+ ALGEBRAIC[214]*ALGEBRAIC[272]+ ALGEBRAIC[199]*ALGEBRAIC[254];
ALGEBRAIC[288] = 1.00000 - ALGEBRAIC[287]*ALGEBRAIC[285];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (CONDVAR[0]<=0.00000||CONDVAR[1]>0.00000&&CONSTANTS[34]==0.00000 ? STATES[0] : CONSTANTS[35]);
ALGEBRAIC[12] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[159])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[160])+ (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[161]);
ALGEBRAIC[13] = 1.00000/ALGEBRAIC[12];
ALGEBRAIC[14] = pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[159])*ALGEBRAIC[13];
ALGEBRAIC[15] = (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[160])*ALGEBRAIC[13];
ALGEBRAIC[16] = (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[161])*ALGEBRAIC[13];
ALGEBRAIC[17] = 0.00000*ALGEBRAIC[13]+ 1.00000*ALGEBRAIC[14]+ 0.00000*ALGEBRAIC[15]+ 0.00000*ALGEBRAIC[16];
ALGEBRAIC[23] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[169])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[170])+ (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[164]);
ALGEBRAIC[24] = 1.00000/ALGEBRAIC[23];
ALGEBRAIC[25] = pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[169])*ALGEBRAIC[24];
ALGEBRAIC[26] = (( ALGEBRAIC[24]*STATES[1])/CONSTANTS[134])*pow(10.0000, CONSTANTS[170]);
ALGEBRAIC[27] = (( ALGEBRAIC[24]*CONSTANTS[133])/CONSTANTS[134])*pow(10.0000, CONSTANTS[164]);
ALGEBRAIC[28] = 0.00000*ALGEBRAIC[24]+ 1.00000*ALGEBRAIC[25]+ 0.00000*ALGEBRAIC[26]+ 0.00000*ALGEBRAIC[27];
ALGEBRAIC[44] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[186])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[187])+ (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[181]);
ALGEBRAIC[45] = 1.00000/ALGEBRAIC[44];
ALGEBRAIC[46] = pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[186])*ALGEBRAIC[45];
ALGEBRAIC[49] = ALGEBRAIC[46];
ALGEBRAIC[54] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[189]);
ALGEBRAIC[55] = 1.00000/ALGEBRAIC[54];
ALGEBRAIC[56] = ALGEBRAIC[55]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[189]);
ALGEBRAIC[57] = 0.00000*ALGEBRAIC[55]+ 1.00000*ALGEBRAIC[56];
ALGEBRAIC[173] = (((ALGEBRAIC[17]+ALGEBRAIC[57]) - ALGEBRAIC[49]) - ALGEBRAIC[28])+(((CONSTANTS[162]+CONSTANTS[190]) - CONSTANTS[188]) - CONSTANTS[171]);
ALGEBRAIC[34] = 1.00000+pow(10.0000, CONSTANTS[177] - ALGEBRAIC[0])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[178]);
ALGEBRAIC[35] = 1.00000/ALGEBRAIC[34];
ALGEBRAIC[36] = ALGEBRAIC[35]*pow(10.0000, CONSTANTS[177] - ALGEBRAIC[0]);
ALGEBRAIC[37] = (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[178])*ALGEBRAIC[35];
ALGEBRAIC[38] = 0.00000*ALGEBRAIC[35]+ALGEBRAIC[36]+ 0.00000*ALGEBRAIC[37];
ALGEBRAIC[175] = (( 2.00000*ALGEBRAIC[28] - ALGEBRAIC[17]) - ALGEBRAIC[38])+(( 2.00000*CONSTANTS[171] - CONSTANTS[162]) - CONSTANTS[179]);
ALGEBRAIC[1] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[149])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[150])+ (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[144]);
ALGEBRAIC[2] = 1.00000/ALGEBRAIC[1];
ALGEBRAIC[3] = pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[149])*ALGEBRAIC[2];
ALGEBRAIC[6] = 1.00000*ALGEBRAIC[3];
ALGEBRAIC[59] = 1.00000+ pow(10.0000, - ALGEBRAIC[0])*pow(10.0000, CONSTANTS[195])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[196]);
ALGEBRAIC[60] = 1.00000/ALGEBRAIC[59];
ALGEBRAIC[61] = ALGEBRAIC[60]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[195]);
ALGEBRAIC[63] = ALGEBRAIC[61];
ALGEBRAIC[178] = (ALGEBRAIC[63] - ALGEBRAIC[6])+(1.00000 - CONSTANTS[151]);
ALGEBRAIC[69] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[199]);
ALGEBRAIC[70] = 1.00000/ALGEBRAIC[69];
ALGEBRAIC[71] = ALGEBRAIC[70]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[199]);
ALGEBRAIC[72] = ALGEBRAIC[71];
ALGEBRAIC[181] = (ALGEBRAIC[72] - ALGEBRAIC[63])+(CONSTANTS[200] - CONSTANTS[197]);
ALGEBRAIC[75] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[202]);
ALGEBRAIC[76] = 1.00000/ALGEBRAIC[75];
ALGEBRAIC[77] = ALGEBRAIC[76]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[202]);
ALGEBRAIC[78] = ALGEBRAIC[77];
ALGEBRAIC[184] = (ALGEBRAIC[78] - ALGEBRAIC[72])+(CONSTANTS[203] - CONSTANTS[200]);
ALGEBRAIC[81] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[205])+pow(10.0000, - 2.00000*ALGEBRAIC[0]+CONSTANTS[205]+CONSTANTS[206])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[207]);
ALGEBRAIC[82] = 1.00000/ALGEBRAIC[81];
ALGEBRAIC[83] = ALGEBRAIC[82]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[205]);
ALGEBRAIC[84] = ALGEBRAIC[82]*pow(10.0000, - 2.00000*ALGEBRAIC[0]+CONSTANTS[205]+CONSTANTS[206]);
ALGEBRAIC[86] = ALGEBRAIC[83]+ 2.00000*ALGEBRAIC[84];
ALGEBRAIC[187] = (((((ALGEBRAIC[28]+ALGEBRAIC[86]) - ALGEBRAIC[78]) - ALGEBRAIC[17])+CONSTANTS[171]+CONSTANTS[208]) - CONSTANTS[203]) - CONSTANTS[162];
ALGEBRAIC[92] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[210]);
ALGEBRAIC[93] = 1.00000/ALGEBRAIC[92];
ALGEBRAIC[94] = ALGEBRAIC[93]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[210]);
ALGEBRAIC[95] = ALGEBRAIC[94];
ALGEBRAIC[108] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[219])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[220]);
ALGEBRAIC[109] = 1.00000/ALGEBRAIC[108];
ALGEBRAIC[110] = ALGEBRAIC[109]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[219]);
ALGEBRAIC[112] = ALGEBRAIC[110];
ALGEBRAIC[190] = ((ALGEBRAIC[112]+ALGEBRAIC[95]) - ALGEBRAIC[86])+((CONSTANTS[221]+CONSTANTS[211]) - CONSTANTS[208]);
ALGEBRAIC[193] = (ALGEBRAIC[112] - ALGEBRAIC[95])+(CONSTANTS[221] - CONSTANTS[211]);
ALGEBRAIC[118] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[223]);
ALGEBRAIC[119] = 1.00000/ALGEBRAIC[118];
ALGEBRAIC[120] = ALGEBRAIC[119]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[223]);
ALGEBRAIC[121] = ALGEBRAIC[120];
ALGEBRAIC[196] = ((ALGEBRAIC[121] - ALGEBRAIC[95]) - ALGEBRAIC[6])+((((CONSTANTS[224]+CONSTANTS[251]) - CONSTANTS[211]) - CONSTANTS[151]) - CONSTANTS[249]);
ALGEBRAIC[98] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[216])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[213]);
ALGEBRAIC[99] = 1.00000/ALGEBRAIC[98];
ALGEBRAIC[100] = ALGEBRAIC[99]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[216]);
ALGEBRAIC[102] = ALGEBRAIC[100];
ALGEBRAIC[199] = (ALGEBRAIC[112] - ALGEBRAIC[102])+(((CONSTANTS[221]+CONSTANTS[251]) - CONSTANTS[249]) - CONSTANTS[217]);
ALGEBRAIC[124] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+6.21000);
ALGEBRAIC[125] = 1.00000/ALGEBRAIC[124];
ALGEBRAIC[126] = ALGEBRAIC[125]*pow(10.0000, - ALGEBRAIC[0]+6.21000);
ALGEBRAIC[127] = ALGEBRAIC[126];
ALGEBRAIC[202] = (((ALGEBRAIC[127]+ALGEBRAIC[17]) - ALGEBRAIC[121]) - ALGEBRAIC[28])+(((CONSTANTS[227]+CONSTANTS[162]) - CONSTANTS[224]) - CONSTANTS[171]);
ALGEBRAIC[130] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[229])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[230])+ (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[231]);
ALGEBRAIC[131] = 1.00000/ALGEBRAIC[130];
ALGEBRAIC[132] = ALGEBRAIC[131]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[229]);
ALGEBRAIC[135] = ALGEBRAIC[132];
ALGEBRAIC[205] = (ALGEBRAIC[135] - ALGEBRAIC[127])+(CONSTANTS[232] - CONSTANTS[227]);
ALGEBRAIC[141] = 1.00000+pow(10.0000, CONSTANTS[234] - ALGEBRAIC[0])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[235])+ (CONSTANTS[133]/CONSTANTS[134])*pow(10.0000, CONSTANTS[236]);
ALGEBRAIC[142] = 1.00000/ALGEBRAIC[141];
ALGEBRAIC[143] = ALGEBRAIC[142]*pow(10.0000, CONSTANTS[234] - ALGEBRAIC[0]);
ALGEBRAIC[146] = ALGEBRAIC[143];
ALGEBRAIC[208] = (ALGEBRAIC[146] - ALGEBRAIC[135])+((CONSTANTS[253]+CONSTANTS[237]) - CONSTANTS[232]);
ALGEBRAIC[152] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[239]);
ALGEBRAIC[153] = 1.00000/ALGEBRAIC[152];
ALGEBRAIC[154] = ALGEBRAIC[153]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[239]);
ALGEBRAIC[155] = ALGEBRAIC[154];
ALGEBRAIC[211] = (((ALGEBRAIC[155]+ALGEBRAIC[17]) - ALGEBRAIC[146]) - ALGEBRAIC[28])+(((CONSTANTS[240]+CONSTANTS[162]) - CONSTANTS[237]) - CONSTANTS[171]);
ALGEBRAIC[158] = 1.00000+pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[245])+ (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[242]);
ALGEBRAIC[159] = 1.00000/ALGEBRAIC[158];
ALGEBRAIC[160] = ALGEBRAIC[159]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[245]);
ALGEBRAIC[162] = ALGEBRAIC[160];
ALGEBRAIC[214] = (ALGEBRAIC[162] - ALGEBRAIC[155])+(((CONSTANTS[246]+CONSTANTS[249]) - CONSTANTS[240]) - CONSTANTS[251]);
ALGEBRAIC[220] = (CONSTANTS[38] != 45.0000 ? CONSTANTS[5] : CONDVAR[2]<0.00000 ? 0.00100000 : CONDVAR[3]>=0.00000&&CONDVAR[4]<0.00000 ? 0.00400000 : CONDVAR[5]>=0.00000&&CONDVAR[6]<0.00000 ? 0.0100000 : CONDVAR[7]>=0.00000 ? 0.0400000 : 0.0/0.0);
ALGEBRAIC[221] = 1.00000+STATES[4]/CONSTANTS[39]+STATES[3]/CONSTANTS[264]+( STATES[4]*STATES[3])/( CONSTANTS[39]*CONSTANTS[40])+STATES[4]/CONSTANTS[266]+STATES[2]/CONSTANTS[41]+( STATES[4]*STATES[2])/( CONSTANTS[267]*CONSTANTS[266]);
ALGEBRAIC[222] = 1.40400/(1.00000+pow(10.0000, 5.94000 - ALGEBRAIC[0])+pow(10.0000, ALGEBRAIC[0] - 7.29000));
ALGEBRAIC[224] = (( ALGEBRAIC[222]*CONSTANTS[263]*STATES[3])/( CONSTANTS[265]*CONSTANTS[264]))/ALGEBRAIC[221];
ALGEBRAIC[11] = (CONSTANTS[152]+ CONSTANTS[151]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[151]);
ALGEBRAIC[68] = (CONSTANTS[198]+ CONSTANTS[197]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[197]);
ALGEBRAIC[168] = 655.700+ CONSTANTS[248]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0];
ALGEBRAIC[179] = (ALGEBRAIC[168]+ALGEBRAIC[68]) - ALGEBRAIC[11];
ALGEBRAIC[180] = ( exp(- ALGEBRAIC[179]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[59])/ALGEBRAIC[1];
ALGEBRAIC[223] = ( ALGEBRAIC[222]*CONSTANTS[263]*CONSTANTS[266]*CONSTANTS[267])/( CONSTANTS[265]*CONSTANTS[264]*ALGEBRAIC[180]);
ALGEBRAIC[225] = (( ALGEBRAIC[223]*STATES[4])/( CONSTANTS[266]*CONSTANTS[267]))/ALGEBRAIC[221];
ALGEBRAIC[226] = ALGEBRAIC[220]*( STATES[4]*ALGEBRAIC[224] - STATES[2]*ALGEBRAIC[225]);
ALGEBRAIC[227] = 1.00000 - ALGEBRAIC[220];
ALGEBRAIC[228] = (pow(STATES[5]/CONSTANTS[44], CONSTANTS[46])/CONSTANTS[45])/(1.00000+pow(STATES[5]/CONSTANTS[44], CONSTANTS[46])/CONSTANTS[45]);
ALGEBRAIC[229] = 1.00000+STATES[4]/CONSTANTS[269]+STATES[3]/CONSTANTS[42]+STATES[4]/CONSTANTS[271]+STATES[2]/CONSTANTS[43]+( STATES[4]*STATES[3])/( CONSTANTS[269]*CONSTANTS[268])+( STATES[4]*STATES[2])/( CONSTANTS[270]*CONSTANTS[271]);
ALGEBRAIC[230] = 1.75000/(1.00000+pow(10.0000, 6.12000 - ALGEBRAIC[0])+pow(10.0000, ALGEBRAIC[0] - 7.03000));
ALGEBRAIC[232] = (( ALGEBRAIC[230]*ALGEBRAIC[228]*CONSTANTS[263]*STATES[3])/( CONSTANTS[269]*CONSTANTS[268]))/ALGEBRAIC[229];
ALGEBRAIC[231] = ( ALGEBRAIC[230]*CONSTANTS[263]*CONSTANTS[270]*CONSTANTS[271])/( CONSTANTS[269]*CONSTANTS[268]*ALGEBRAIC[180]);
ALGEBRAIC[233] = (( ALGEBRAIC[228]*ALGEBRAIC[231]*STATES[4])/( CONSTANTS[270]*CONSTANTS[271]))/ALGEBRAIC[229];
ALGEBRAIC[234] = ALGEBRAIC[227]*( STATES[4]*ALGEBRAIC[232] - STATES[2]*ALGEBRAIC[233]);
ALGEBRAIC[235] = ( CONSTANTS[272]*1.32900)/(1.00000+pow(10.0000, - ALGEBRAIC[0]+6.64000)+pow(10.0000, ALGEBRAIC[0] - 8.36000));
ALGEBRAIC[74] = (CONSTANTS[201]+ CONSTANTS[200]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[200]);
ALGEBRAIC[182] = ALGEBRAIC[74] - ALGEBRAIC[68];
ALGEBRAIC[183] = ( exp(- ALGEBRAIC[182]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[69])/ALGEBRAIC[59];
ALGEBRAIC[236] = ( ALGEBRAIC[235]*CONSTANTS[274])/( CONSTANTS[273]*ALGEBRAIC[183]);
ALGEBRAIC[237] = (( ALGEBRAIC[235]*STATES[2])/CONSTANTS[273] - ( ALGEBRAIC[236]*STATES[6])/CONSTANTS[274])/(1.00000+STATES[2]/CONSTANTS[273]+STATES[6]/CONSTANTS[274]);
ALGEBRAIC[238] = CONSTANTS[275]/(1.00000+pow(10.0000, - ALGEBRAIC[0]+6.94000)+pow(10.0000, ALGEBRAIC[0] - 9.35000));
ALGEBRAIC[80] = (CONSTANTS[204]+ CONSTANTS[203]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[203]);
ALGEBRAIC[185] = ALGEBRAIC[80] - ALGEBRAIC[74];
ALGEBRAIC[186] = ( exp(- ALGEBRAIC[185]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[75])/ALGEBRAIC[69];
ALGEBRAIC[239] = (( ALGEBRAIC[238]*CONSTANTS[276])/CONSTANTS[277])*ALGEBRAIC[186];
ALGEBRAIC[240] = (( ALGEBRAIC[239]*STATES[6])/CONSTANTS[276] - ( ALGEBRAIC[238]*STATES[7])/CONSTANTS[277])/(1.00000+STATES[7]/CONSTANTS[277]+STATES[6]/CONSTANTS[276]);
ALGEBRAIC[241] = CONSTANTS[278]/(1.00000+pow(ALGEBRAIC[0]/6.80000, - 30.0000));
ALGEBRAIC[22] = (CONSTANTS[163]+ CONSTANTS[162]*CONSTANTS[0]*CONSTANTS[130]*log(10.0000)*ALGEBRAIC[0]) - CONSTANTS[143]*(16.0000 - CONSTANTS[162]);
ALGEBRAIC[33] = (CONSTANTS[172]+ CONSTANTS[171]*CONSTANTS[0]*CONSTANTS[130]*log(10.0000)*ALGEBRAIC[0]) - CONSTANTS[143]*(9.00000 - CONSTANTS[171]);
ALGEBRAIC[91] = (CONSTANTS[209]+ CONSTANTS[208]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(16.0000 - CONSTANTS[208]);
ALGEBRAIC[172] = CONSTANTS[256]+ CONSTANTS[255]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0];
ALGEBRAIC[188] = ((ALGEBRAIC[91]+ALGEBRAIC[33]+ALGEBRAIC[172]) - ALGEBRAIC[80]) - ALGEBRAIC[22];
ALGEBRAIC[189] = ( exp(- ALGEBRAIC[188]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[81]*ALGEBRAIC[23])/( ALGEBRAIC[75]*ALGEBRAIC[12]*pow(10.0000, - ALGEBRAIC[0]));
ALGEBRAIC[242] = ( ALGEBRAIC[241]*CONSTANTS[282]*CONSTANTS[284])/( CONSTANTS[279]*CONSTANTS[281]*ALGEBRAIC[189]);
ALGEBRAIC[243] = CONSTANTS[53]*pow(( ((1.00000+STATES[10]/CONSTANTS[49])/(1.00000+( CONSTANTS[51]*STATES[10])/CONSTANTS[49]))*(1.00000+( CONSTANTS[52]*STATES[5])/CONSTANTS[50]))/(1.00000+STATES[5]/CONSTANTS[50]), 4.00000);
ALGEBRAIC[244] = (1.00000+STATES[7]/CONSTANTS[279])*(1.00000+STATES[10]/CONSTANTS[281])+STATES[8]/CONSTANTS[282]+ (STATES[9]/CONSTANTS[284])*(1.00000+STATES[8]/CONSTANTS[282]);
ALGEBRAIC[245] = (1.00000+STATES[7]/CONSTANTS[280])*(1.00000+STATES[10]/CONSTANTS[283])+STATES[8]/CONSTANTS[47]+ (STATES[9]/CONSTANTS[48])*(1.00000+STATES[8]/CONSTANTS[47]);
ALGEBRAIC[246] = ( ((( ALGEBRAIC[241]*STATES[7]*STATES[10])/( CONSTANTS[279]*CONSTANTS[281]))/ALGEBRAIC[244] - (( ALGEBRAIC[242]*STATES[9]*STATES[8])/( CONSTANTS[284]*CONSTANTS[282]))/ALGEBRAIC[244])*(1.00000+ CONSTANTS[285]*ALGEBRAIC[243]*pow(ALGEBRAIC[245]/ALGEBRAIC[244], 3.00000)))/(1.00000+ ALGEBRAIC[243]*pow(ALGEBRAIC[245]/ALGEBRAIC[244], 4.00000));
ALGEBRAIC[247] = ( CONSTANTS[286]*1.01300)/(1.00000+pow(10.0000, - ALGEBRAIC[0]+5.32000)+pow(10.0000, ALGEBRAIC[0] - 9.15000));
ALGEBRAIC[97] = (CONSTANTS[212]+ CONSTANTS[211]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[211]);
ALGEBRAIC[117] = (CONSTANTS[222]+ CONSTANTS[221]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[221]);
ALGEBRAIC[191] = (ALGEBRAIC[117]+ALGEBRAIC[97]) - ALGEBRAIC[91];
ALGEBRAIC[192] = ( 1.00000*exp(- ALGEBRAIC[191]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[92]*ALGEBRAIC[81])/ALGEBRAIC[108];
ALGEBRAIC[248] = ( ALGEBRAIC[247]*CONSTANTS[289]*CONSTANTS[288])/( CONSTANTS[287]*ALGEBRAIC[192]);
ALGEBRAIC[249] = (( ALGEBRAIC[247]*STATES[8])/CONSTANTS[287] - ( ALGEBRAIC[248]*STATES[12]*STATES[11])/( CONSTANTS[289]*CONSTANTS[288]))/(1.00000+STATES[8]/CONSTANTS[287]+STATES[12]/CONSTANTS[289]+STATES[11]/CONSTANTS[288]);
ALGEBRAIC[194] = ALGEBRAIC[117] - ALGEBRAIC[97];
ALGEBRAIC[195] = ( exp(- ALGEBRAIC[194]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[108])/ALGEBRAIC[92];
ALGEBRAIC[250] = ( CONSTANTS[293]*CONSTANTS[292])/( CONSTANTS[291]*ALGEBRAIC[195]);
ALGEBRAIC[251] = (( CONSTANTS[293]*STATES[12])/CONSTANTS[291] - ( ALGEBRAIC[250]*STATES[11])/CONSTANTS[292])/(1.00000+STATES[12]/CONSTANTS[291]+STATES[11]/CONSTANTS[292]);
ALGEBRAIC[252] = (1.00000+STATES[13]/CONSTANTS[295]+STATES[15]/CONSTANTS[298])*(1.00000+STATES[11]/CONSTANTS[297]+STATES[14]/CONSTANTS[298]);
ALGEBRAIC[107] = (CONSTANTS[218]+ CONSTANTS[217]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[217]);
ALGEBRAIC[169] = (CONSTANTS[250]+ CONSTANTS[249]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(1.00000 - CONSTANTS[249]);
ALGEBRAIC[170] = (CONSTANTS[252]+ CONSTANTS[251]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[251]);
ALGEBRAIC[200] = ((ALGEBRAIC[172]+ALGEBRAIC[170]+ALGEBRAIC[117]) - ALGEBRAIC[169]) - ALGEBRAIC[107];
ALGEBRAIC[201] = ( exp(- ALGEBRAIC[200]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[108])/( ALGEBRAIC[98]*pow(10.0000, - ALGEBRAIC[0]));
ALGEBRAIC[253] = ( CONSTANTS[299]*CONSTANTS[295]*CONSTANTS[296]*ALGEBRAIC[201])/( CONSTANTS[297]*CONSTANTS[298]);
ALGEBRAIC[254] = (( ALGEBRAIC[253]*STATES[13]*STATES[14])/( CONSTANTS[295]*CONSTANTS[296]) - ( CONSTANTS[299]*STATES[11]*STATES[15])/( CONSTANTS[297]*CONSTANTS[298]))/ALGEBRAIC[252];
ALGEBRAIC[255] = 1.00000+STATES[3]/CONSTANTS[303]+STATES[12]/CONSTANTS[301]+STATES[14]/CONSTANTS[302]+( STATES[12]*STATES[14])/( CONSTANTS[301]*CONSTANTS[302])+( STATES[12]*STATES[14]*STATES[3])/( CONSTANTS[301]*CONSTANTS[302]*CONSTANTS[303])+STATES[16]/CONSTANTS[304]+STATES[15]/CONSTANTS[305]+( STATES[16]*STATES[15])/( CONSTANTS[305]*CONSTANTS[304]);
ALGEBRAIC[256] = CONSTANTS[300]*0.000700000*exp( ALGEBRAIC[0]*0.897900);
ALGEBRAIC[123] = (CONSTANTS[225]+ CONSTANTS[224]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(16.0000 - CONSTANTS[224]);
ALGEBRAIC[197] = (((ALGEBRAIC[123]+ALGEBRAIC[170]+ALGEBRAIC[172]) - ALGEBRAIC[11]) - ALGEBRAIC[97]) - ALGEBRAIC[169];
ALGEBRAIC[198] = ( exp(- ALGEBRAIC[197]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[118])/( ALGEBRAIC[1]*ALGEBRAIC[92]*pow(10.0000, - ALGEBRAIC[0])*1.00000);
ALGEBRAIC[257] = ( ALGEBRAIC[256]*CONSTANTS[304]*CONSTANTS[305])/( CONSTANTS[301]*CONSTANTS[303]*CONSTANTS[302]*ALGEBRAIC[198]);
ALGEBRAIC[258] = (( ALGEBRAIC[256]*STATES[12]*STATES[14]*STATES[3])/( CONSTANTS[302]*CONSTANTS[301]*CONSTANTS[303]) - ( ALGEBRAIC[257]*STATES[16]*STATES[15])/( CONSTANTS[304]*CONSTANTS[305]))/ALGEBRAIC[255];
ALGEBRAIC[129] = (CONSTANTS[228]+ CONSTANTS[227]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(9.00000 - CONSTANTS[227]);
ALGEBRAIC[203] = ((ALGEBRAIC[22]+ALGEBRAIC[129]) - ALGEBRAIC[123]) - ALGEBRAIC[33];
ALGEBRAIC[204] = ( exp(- ALGEBRAIC[203]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[12]*ALGEBRAIC[124])/( ALGEBRAIC[118]*ALGEBRAIC[23]);
ALGEBRAIC[259] = (( CONSTANTS[311]*CONSTANTS[307]*CONSTANTS[308])/( CONSTANTS[309]*CONSTANTS[310]))*ALGEBRAIC[204];
ALGEBRAIC[260] = 1.00000+STATES[9]/CONSTANTS[308]+STATES[16]/CONSTANTS[307]+( STATES[16]*STATES[9])/( CONSTANTS[307]*CONSTANTS[308])+STATES[17]/CONSTANTS[309]+STATES[10]/CONSTANTS[310]+( STATES[17]*STATES[10])/( CONSTANTS[309]*CONSTANTS[310]);
ALGEBRAIC[261] = (( ALGEBRAIC[259]*STATES[16]*STATES[9])/( CONSTANTS[308]*CONSTANTS[307]) - ( CONSTANTS[311]*STATES[10]*STATES[17])/( CONSTANTS[310]*CONSTANTS[309]))/ALGEBRAIC[260];
ALGEBRAIC[262] = ( CONSTANTS[312]*0.989000)/(1.00000+pow(10.0000, - ALGEBRAIC[0]+5.62000)+pow(10.0000, ALGEBRAIC[0] - 8.74000));
ALGEBRAIC[140] = (CONSTANTS[233]+ CONSTANTS[232]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(9.00000 - CONSTANTS[232]);
ALGEBRAIC[206] = ALGEBRAIC[140] - ALGEBRAIC[129];
ALGEBRAIC[207] = ( exp(- ALGEBRAIC[206]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[130])/ALGEBRAIC[124];
ALGEBRAIC[263] = ( ALGEBRAIC[262]*CONSTANTS[314])/( CONSTANTS[313]*ALGEBRAIC[207]);
ALGEBRAIC[264] = (( ALGEBRAIC[262]*STATES[17])/CONSTANTS[313] - ( ALGEBRAIC[263]*STATES[18])/CONSTANTS[314])/(1.00000+STATES[17]/CONSTANTS[313]+STATES[18]/CONSTANTS[314]);
ALGEBRAIC[151] = (CONSTANTS[238]+ CONSTANTS[237]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(9.00000 - CONSTANTS[237]);
ALGEBRAIC[171] = (CONSTANTS[254]+ CONSTANTS[253]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(0.00000 - CONSTANTS[253]);
ALGEBRAIC[209] = (ALGEBRAIC[171]+ALGEBRAIC[151]) - ALGEBRAIC[140];
ALGEBRAIC[210] = ( exp(- ALGEBRAIC[209]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[141])/ALGEBRAIC[130];
ALGEBRAIC[265] = ( CONSTANTS[318]*CONSTANTS[317])/( CONSTANTS[316]*ALGEBRAIC[210]);
ALGEBRAIC[266] = (( CONSTANTS[318]*STATES[18])/CONSTANTS[316] - ( ALGEBRAIC[265]*STATES[19])/CONSTANTS[317])/(1.00000+STATES[19]/CONSTANTS[317]+STATES[18]/CONSTANTS[316]);
ALGEBRAIC[267] = ( CONSTANTS[319]*1.05000)/(1.00000+pow(10.0000, - ALGEBRAIC[0]+5.58000)+pow(10.0000, ALGEBRAIC[0] - 8.79000));
ALGEBRAIC[157] = (CONSTANTS[241]+ CONSTANTS[240]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(1.00000 - CONSTANTS[240]);
ALGEBRAIC[212] = (((ALGEBRAIC[157]+ALGEBRAIC[22]) - ALGEBRAIC[172]) - ALGEBRAIC[151]) - ALGEBRAIC[33];
ALGEBRAIC[213] = ( exp(- ALGEBRAIC[212]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[152]*ALGEBRAIC[12]*pow(10.0000, - ALGEBRAIC[0]))/( ALGEBRAIC[141]*ALGEBRAIC[23]);
ALGEBRAIC[268] = ( ALGEBRAIC[267]*CONSTANTS[322]*CONSTANTS[323])/( CONSTANTS[320]*CONSTANTS[321]*ALGEBRAIC[213]);
ALGEBRAIC[269] = (( ALGEBRAIC[267]*STATES[19]*STATES[9])/( CONSTANTS[320]*CONSTANTS[321]) - ( ALGEBRAIC[268]*STATES[20]*STATES[10])/( CONSTANTS[322]*CONSTANTS[323]))/(1.00000+STATES[19]/CONSTANTS[320]+STATES[9]/CONSTANTS[321]+( STATES[19]*STATES[9])/( CONSTANTS[320]*CONSTANTS[321])+STATES[10]/CONSTANTS[323]+STATES[20]/CONSTANTS[322]+( STATES[20]*STATES[10])/( CONSTANTS[322]*CONSTANTS[323]));
ALGEBRAIC[270] = CONSTANTS[324]*( - 0.113400*ALGEBRAIC[0]+1.60690);
ALGEBRAIC[167] = (CONSTANTS[247]+ CONSTANTS[246]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(1.00000 - CONSTANTS[246]);
ALGEBRAIC[215] = (((ALGEBRAIC[167]+ALGEBRAIC[169]) - ALGEBRAIC[157]) - ALGEBRAIC[170]) - ALGEBRAIC[172];
ALGEBRAIC[216] = ( exp(- ALGEBRAIC[215]/( CONSTANTS[0]*CONSTANTS[130]))*ALGEBRAIC[158]*pow(10.0000, - ALGEBRAIC[0]))/ALGEBRAIC[152];
ALGEBRAIC[271] = ( ALGEBRAIC[270]*CONSTANTS[327]*CONSTANTS[328])/( CONSTANTS[325]*CONSTANTS[326]*ALGEBRAIC[216]);
ALGEBRAIC[272] = (( ALGEBRAIC[270]*STATES[20]*STATES[15])/( CONSTANTS[325]*CONSTANTS[326]) - ( ALGEBRAIC[271]*STATES[21]*STATES[14])/( CONSTANTS[327]*CONSTANTS[328]))/(1.00000+STATES[20]/CONSTANTS[325]+STATES[15]/CONSTANTS[326]+( STATES[20]*STATES[15])/( CONSTANTS[325]*CONSTANTS[326])+STATES[21]/CONSTANTS[327]+STATES[14]/CONSTANTS[328]+( STATES[21]*STATES[14])/( CONSTANTS[327]*CONSTANTS[328]));
ALGEBRAIC[174] = ( exp(- CONSTANTS[262]/( CONSTANTS[0]*CONSTANTS[130]))*pow(10.0000, - ALGEBRAIC[0])*ALGEBRAIC[12]*ALGEBRAIC[54])/( ALGEBRAIC[44]*ALGEBRAIC[23]);
ALGEBRAIC[275] = ( (CONSTANTS[329]/ALGEBRAIC[174])*CONSTANTS[331]*CONSTANTS[333])/( CONSTANTS[332]*CONSTANTS[330]);
ALGEBRAIC[276] = (( ALGEBRAIC[275]*STATES[10]*STATES[22])/( CONSTANTS[331]*CONSTANTS[333]) - ( CONSTANTS[329]*STATES[9]*STATES[23])/( CONSTANTS[332]*CONSTANTS[330]))/(1.00000+STATES[9]/CONSTANTS[332]+STATES[23]/CONSTANTS[55]+( STATES[23]*STATES[9])/( CONSTANTS[332]*CONSTANTS[330])+STATES[10]/CONSTANTS[331]+( STATES[22]*STATES[10])/( CONSTANTS[333]*CONSTANTS[331]));
ALGEBRAIC[43] = (CONSTANTS[180]+ CONSTANTS[179]*log(10.0000)*CONSTANTS[0]*CONSTANTS[130]*ALGEBRAIC[0]) - CONSTANTS[143]*(4.00000 - CONSTANTS[179]);
ALGEBRAIC[176] = ( 2.00000*ALGEBRAIC[33] - ALGEBRAIC[22]) - ALGEBRAIC[43];
ALGEBRAIC[177] = ( exp(- ALGEBRAIC[176]/( CONSTANTS[0]*CONSTANTS[130]))*pow(ALGEBRAIC[23], 2.00000))/( ALGEBRAIC[12]*ALGEBRAIC[34]);
ALGEBRAIC[277] = ( CONSTANTS[334]*pow(CONSTANTS[337], 2.00000))/( CONSTANTS[335]*CONSTANTS[336]*ALGEBRAIC[177]);
ALGEBRAIC[278] = (( CONSTANTS[334]*STATES[10]*STATES[5])/( CONSTANTS[336]*CONSTANTS[335]) - ALGEBRAIC[277]*pow(STATES[9]/CONSTANTS[337], 2.00000))/(1.00000+STATES[10]/CONSTANTS[336]+STATES[5]/CONSTANTS[335]+( STATES[10]*STATES[5])/( CONSTANTS[336]*CONSTANTS[335])+( 2.00000*STATES[9])/CONSTANTS[337]+pow(STATES[9], 2.00000)/pow(CONSTANTS[337], 2.00000));
ALGEBRAIC[281] = ALGEBRAIC[173]*- ALGEBRAIC[276]+ ALGEBRAIC[175]*ALGEBRAIC[278]+ ALGEBRAIC[181]*ALGEBRAIC[237]+ ALGEBRAIC[178]*(ALGEBRAIC[226]+ALGEBRAIC[234])+ ALGEBRAIC[184]*ALGEBRAIC[240]+ ALGEBRAIC[187]*ALGEBRAIC[246]+ ALGEBRAIC[190]*ALGEBRAIC[249]+ ALGEBRAIC[193]*ALGEBRAIC[251]+ ALGEBRAIC[196]*ALGEBRAIC[258]+ ALGEBRAIC[202]*ALGEBRAIC[261]+ ALGEBRAIC[205]*ALGEBRAIC[264]+ ALGEBRAIC[208]*ALGEBRAIC[266]+ ALGEBRAIC[211]*ALGEBRAIC[269]+ ALGEBRAIC[214]*ALGEBRAIC[272]+ ALGEBRAIC[199]*ALGEBRAIC[254];
ALGEBRAIC[279] = ( log(10.0000)*CONSTANTS[338]*pow(10.0000, - ALGEBRAIC[0] - 6.87000))/pow(pow(10.0000, - ALGEBRAIC[0])+pow(10.0000, - 6.87000), 2.00000)+( log(10.0000)*CONSTANTS[339]*pow(10.0000, - ALGEBRAIC[0] - 8.30000))/pow(pow(10.0000, - ALGEBRAIC[0])+pow(10.0000, - 8.30000), 2.00000)+( log(10.0000)*CONSTANTS[340]*pow(10.0000, - ALGEBRAIC[0] - 4.80000))/pow(pow(10.0000, - ALGEBRAIC[0])+pow(10.0000, - 4.80000), 2.00000);
ALGEBRAIC[7] = ( pow(10.0000, CONSTANTS[149])*(ALGEBRAIC[1] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[149])))/( CONSTANTS[134]*pow(ALGEBRAIC[1], 2.00000));
ALGEBRAIC[18] = ( pow(10.0000, CONSTANTS[159])*(ALGEBRAIC[12] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[159])))/( CONSTANTS[134]*pow(ALGEBRAIC[12], 2.00000));
ALGEBRAIC[29] = ( pow(10.0000, CONSTANTS[169])*(ALGEBRAIC[23] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[169])))/( CONSTANTS[134]*pow(ALGEBRAIC[23], 2.00000));
ALGEBRAIC[39] = ( pow(10.0000, CONSTANTS[177])*(ALGEBRAIC[34] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[177])))/( CONSTANTS[134]*pow(ALGEBRAIC[34], 2.00000));
ALGEBRAIC[50] = ( pow(10.0000, CONSTANTS[186])*(ALGEBRAIC[44] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[186])))/( CONSTANTS[134]*pow(ALGEBRAIC[44], 2.00000));
ALGEBRAIC[58] = ( pow(10.0000, CONSTANTS[189])*(ALGEBRAIC[54] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[189])))/( CONSTANTS[134]*pow(ALGEBRAIC[54], 2.00000));
ALGEBRAIC[64] = ( pow(10.0000, CONSTANTS[195])*(ALGEBRAIC[59] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[195])))/( CONSTANTS[134]*pow(ALGEBRAIC[59], 2.00000));
ALGEBRAIC[73] = ( pow(10.0000, CONSTANTS[199])*(ALGEBRAIC[69] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[199])))/( CONSTANTS[134]*pow(ALGEBRAIC[69], 2.00000));
ALGEBRAIC[79] = ( pow(10.0000, CONSTANTS[202])*(ALGEBRAIC[75] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[202])))/( CONSTANTS[134]*pow(ALGEBRAIC[75], 2.00000));
ALGEBRAIC[87] = ( ALGEBRAIC[81]*(pow(10.0000, CONSTANTS[205])+ 2.00000*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[205]+CONSTANTS[206])) - (pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[205])+ 2.00000*pow(10.0000, - 2.00000*ALGEBRAIC[0]+CONSTANTS[205]+CONSTANTS[206]))*(pow(10.0000, CONSTANTS[205])+ 2.00000*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[205]+CONSTANTS[206])))/( CONSTANTS[134]*pow(ALGEBRAIC[81], 2.00000));
ALGEBRAIC[96] = ( pow(10.0000, CONSTANTS[210])*(ALGEBRAIC[92] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[210])))/( CONSTANTS[134]*pow(ALGEBRAIC[92], 2.00000));
ALGEBRAIC[103] = ( pow(10.0000, CONSTANTS[216])*(ALGEBRAIC[98] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[216])))/( CONSTANTS[134]*pow(ALGEBRAIC[98], 2.00000));
ALGEBRAIC[113] = ( pow(10.0000, CONSTANTS[219])*(ALGEBRAIC[108] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[219])))/( CONSTANTS[134]*pow(ALGEBRAIC[108], 2.00000));
ALGEBRAIC[122] = ( pow(10.0000, CONSTANTS[223])*(ALGEBRAIC[118] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[223])))/( CONSTANTS[134]*pow(ALGEBRAIC[118], 2.00000));
ALGEBRAIC[128] = ( pow(10.0000, CONSTANTS[226])*(ALGEBRAIC[124] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[226])))/( CONSTANTS[134]*pow(ALGEBRAIC[124], 2.00000));
ALGEBRAIC[136] = ( pow(10.0000, CONSTANTS[229])*(ALGEBRAIC[130] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[229])))/( CONSTANTS[134]*pow(ALGEBRAIC[130], 2.00000));
ALGEBRAIC[147] = ( pow(10.0000, CONSTANTS[234])*(ALGEBRAIC[141] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[234])))/( CONSTANTS[134]*pow(ALGEBRAIC[141], 2.00000));
ALGEBRAIC[156] = ( pow(10.0000, CONSTANTS[239])*(ALGEBRAIC[152] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[239])))/( CONSTANTS[134]*pow(ALGEBRAIC[152], 2.00000));
ALGEBRAIC[163] = ( pow(10.0000, CONSTANTS[245])*(ALGEBRAIC[158] - pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[245])))/( CONSTANTS[134]*pow(ALGEBRAIC[158], 2.00000));
ALGEBRAIC[280] = log(10.0000)*pow(10.0000, - ALGEBRAIC[0])*CONSTANTS[134]*(1.00000+ ALGEBRAIC[7]*STATES[3]+ ALGEBRAIC[18]*STATES[10]+ ALGEBRAIC[29]*STATES[9]+ ALGEBRAIC[39]*STATES[5]+ ALGEBRAIC[50]*STATES[23]+ ALGEBRAIC[58]*STATES[22]+ ALGEBRAIC[64]*STATES[2]+ ALGEBRAIC[73]*STATES[6]+ ALGEBRAIC[79]*STATES[7]+ ALGEBRAIC[87]*STATES[8]+ ALGEBRAIC[96]*STATES[12]+ ALGEBRAIC[113]*STATES[11]+ ALGEBRAIC[103]*STATES[13]+ ALGEBRAIC[122]*STATES[16]+ ALGEBRAIC[128]*STATES[17]+ ALGEBRAIC[136]*STATES[18]+ ALGEBRAIC[147]*STATES[19]+ ALGEBRAIC[156]*STATES[20]+ ALGEBRAIC[163]*STATES[21]);
ALGEBRAIC[284] = ALGEBRAIC[281]/(ALGEBRAIC[279]+ALGEBRAIC[280]);
ALGEBRAIC[8] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[149]+CONSTANTS[150]))/pow(ALGEBRAIC[1], 2.00000);
ALGEBRAIC[19] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[159]+CONSTANTS[160]))/pow(ALGEBRAIC[12], 2.00000);
ALGEBRAIC[30] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[169]+CONSTANTS[170]))/pow(ALGEBRAIC[23], 2.00000);
ALGEBRAIC[40] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[177]+CONSTANTS[178]))/pow(ALGEBRAIC[34], 2.00000);
ALGEBRAIC[51] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[186]+CONSTANTS[187]))/pow(ALGEBRAIC[44], 2.00000);
ALGEBRAIC[65] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[195]+CONSTANTS[196]))/pow(ALGEBRAIC[59], 2.00000);
ALGEBRAIC[88] = ( - (pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[205])+ 2.00000*pow(10.0000, - 2.00000*ALGEBRAIC[0]+CONSTANTS[205]+CONSTANTS[206]))*pow(10.0000, CONSTANTS[207]))/( CONSTANTS[134]*pow(ALGEBRAIC[81], 2.00000));
ALGEBRAIC[104] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[216]+CONSTANTS[213]))/pow(ALGEBRAIC[98], 2.00000);
ALGEBRAIC[114] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[219]+CONSTANTS[220]))/pow(ALGEBRAIC[108], 2.00000);
ALGEBRAIC[137] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[229]+CONSTANTS[230]))/pow(ALGEBRAIC[130], 2.00000);
ALGEBRAIC[148] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[234]+CONSTANTS[235]))/pow(ALGEBRAIC[141], 2.00000);
ALGEBRAIC[164] = ( (- pow(10.0000, - ALGEBRAIC[0])/CONSTANTS[134])*pow(10.0000, CONSTANTS[245]+CONSTANTS[242]))/pow(ALGEBRAIC[158], 2.00000);
ALGEBRAIC[285] = ( ALGEBRAIC[8]*STATES[3]+ ALGEBRAIC[19]*STATES[10]+ ALGEBRAIC[30]*STATES[9]+ ALGEBRAIC[40]*STATES[5]+ ALGEBRAIC[51]*STATES[23]+ CONSTANTS[26]*STATES[22]+ ALGEBRAIC[65]*STATES[2]+ CONSTANTS[27]*STATES[6]+ CONSTANTS[28]*STATES[7]+ ALGEBRAIC[88]*STATES[8]+ CONSTANTS[29]*STATES[12]+ ALGEBRAIC[114]*STATES[11]+ ALGEBRAIC[104]*STATES[13]+ CONSTANTS[30]*STATES[16]+ CONSTANTS[31]*STATES[17]+ ALGEBRAIC[137]*STATES[18]+ ALGEBRAIC[148]*STATES[19]+ CONSTANTS[32]*STATES[20]+ ALGEBRAIC[164]*STATES[21])/(ALGEBRAIC[279]+ALGEBRAIC[280]);
ALGEBRAIC[10] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[149]+CONSTANTS[150])*log(10.0000))/pow(ALGEBRAIC[1], 2.00000);
ALGEBRAIC[21] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[159]+CONSTANTS[160])*log(10.0000))/pow(ALGEBRAIC[12], 2.00000);
ALGEBRAIC[32] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[169]+CONSTANTS[170])*log(10.0000))/pow(ALGEBRAIC[23], 2.00000);
ALGEBRAIC[42] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[177]+CONSTANTS[178])*log(10.0000))/pow(ALGEBRAIC[34], 2.00000);
ALGEBRAIC[53] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[186]+CONSTANTS[187])*log(10.0000))/pow(ALGEBRAIC[44], 2.00000);
ALGEBRAIC[67] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[195]+CONSTANTS[196])*log(10.0000))/pow(ALGEBRAIC[59], 2.00000);
ALGEBRAIC[90] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[207])*( ALGEBRAIC[81]*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[205])*log(10.0000) - pow(10.0000, - 2.00000*ALGEBRAIC[0]+CONSTANTS[205]+CONSTANTS[206])*2.00000*log(10.0000)))/pow(ALGEBRAIC[81], 2.00000);
ALGEBRAIC[106] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[216]+CONSTANTS[213])*log(10.0000))/pow(ALGEBRAIC[98], 2.00000);
ALGEBRAIC[116] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[219]+CONSTANTS[220])*log(10.0000))/pow(ALGEBRAIC[108], 2.00000);
ALGEBRAIC[139] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[229]+CONSTANTS[230])*log(10.0000))/pow(ALGEBRAIC[130], 2.00000);
ALGEBRAIC[150] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[234]+CONSTANTS[235])*log(10.0000))/pow(ALGEBRAIC[141], 2.00000);
ALGEBRAIC[166] = ( (STATES[1]/CONSTANTS[134])*pow(10.0000, - ALGEBRAIC[0]+CONSTANTS[245]+CONSTANTS[242])*log(10.0000))/pow(ALGEBRAIC[158], 2.00000);
ALGEBRAIC[9] = (( ALGEBRAIC[1]*pow(10.0000, CONSTANTS[150]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[150]))/CONSTANTS[134])/pow(ALGEBRAIC[1], 2.00000);
ALGEBRAIC[20] = (( ALGEBRAIC[12]*pow(10.0000, CONSTANTS[160]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[160]))/CONSTANTS[134])/pow(ALGEBRAIC[12], 2.00000);
ALGEBRAIC[31] = (( ALGEBRAIC[23]*pow(10.0000, CONSTANTS[170]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[170]))/CONSTANTS[134])/pow(ALGEBRAIC[23], 2.00000);
ALGEBRAIC[41] = (( ALGEBRAIC[34]*pow(10.0000, CONSTANTS[178]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[178]))/CONSTANTS[134])/pow(ALGEBRAIC[34], 2.00000);
ALGEBRAIC[52] = (( ALGEBRAIC[44]*pow(10.0000, CONSTANTS[187]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[187]))/CONSTANTS[134])/pow(ALGEBRAIC[44], 2.00000);
ALGEBRAIC[66] = (( ALGEBRAIC[59]*pow(10.0000, CONSTANTS[196]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[196]))/CONSTANTS[134])/pow(ALGEBRAIC[59], 2.00000);
ALGEBRAIC[89] = (( ALGEBRAIC[81]*pow(10.0000, CONSTANTS[207]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[207]))/CONSTANTS[134])/pow(ALGEBRAIC[81], 2.00000);
ALGEBRAIC[105] = (( ALGEBRAIC[98]*pow(10.0000, CONSTANTS[213]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[213]))/CONSTANTS[134])/pow(ALGEBRAIC[98], 2.00000);
ALGEBRAIC[115] = (( ALGEBRAIC[108]*pow(10.0000, CONSTANTS[220]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[220]))/CONSTANTS[134])/pow(ALGEBRAIC[108], 2.00000);
ALGEBRAIC[138] = (( ALGEBRAIC[130]*pow(10.0000, CONSTANTS[230]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[230]))/CONSTANTS[134])/pow(ALGEBRAIC[130], 2.00000);
ALGEBRAIC[149] = (( ALGEBRAIC[141]*pow(10.0000, CONSTANTS[235]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[235]))/CONSTANTS[134])/pow(ALGEBRAIC[141], 2.00000);
ALGEBRAIC[165] = (( ALGEBRAIC[158]*pow(10.0000, CONSTANTS[242]))/CONSTANTS[134] - ( (STATES[1]/CONSTANTS[134])*pow(10.0000, 2.00000*CONSTANTS[242]))/CONSTANTS[134])/pow(ALGEBRAIC[158], 2.00000);
ALGEBRAIC[286] = - 1.00000 - ( ALGEBRAIC[20]*STATES[10]+ ALGEBRAIC[31]*STATES[9]+ ALGEBRAIC[41]*STATES[5]+ ALGEBRAIC[9]*STATES[3]+ ALGEBRAIC[52]*STATES[23]+ ALGEBRAIC[66]*STATES[2]+ ALGEBRAIC[89]*STATES[8]+ ALGEBRAIC[105]*STATES[13]+ ALGEBRAIC[115]*STATES[11]+ ALGEBRAIC[149]*STATES[19]+ ALGEBRAIC[138]*STATES[18]+ ALGEBRAIC[165]*STATES[21]);
ALGEBRAIC[287] = ( ALGEBRAIC[21]*STATES[10]+ ALGEBRAIC[32]*STATES[9]+ ALGEBRAIC[42]*STATES[5]+ ALGEBRAIC[10]*STATES[3]+ ALGEBRAIC[53]*STATES[23]+ ALGEBRAIC[67]*STATES[2]+ ALGEBRAIC[90]*STATES[8]+ ALGEBRAIC[106]*STATES[13]+ ALGEBRAIC[116]*STATES[11]+ ALGEBRAIC[150]*STATES[19]+ ALGEBRAIC[139]*STATES[18]+ ALGEBRAIC[166]*STATES[21])/ALGEBRAIC[286];
ALGEBRAIC[289] = ALGEBRAIC[276];
ALGEBRAIC[290] = - ALGEBRAIC[276];
ALGEBRAIC[291] = (- ALGEBRAIC[258] - ALGEBRAIC[254])+ALGEBRAIC[272];
ALGEBRAIC[292] = (ALGEBRAIC[258]+ALGEBRAIC[254]) - ALGEBRAIC[272];
ALGEBRAIC[273] = (CONSTANTS[38] != 5.00000&&CONSTANTS[38] != 45.0000 ? CONSTANTS[22] : CONSTANTS[38] != 5.00000&&CONSTANTS[38]==45.0000&&CONDVAR[8]>=0.00000 ? CONSTANTS[22] : CONSTANTS[38]==5.00000&&CONDVAR[9]>0.00000 ? CONSTANTS[22] : 0.00000);
ALGEBRAIC[274] = ( ALGEBRAIC[273]*STATES[10])/(CONSTANTS[54]+STATES[10]);
ALGEBRAIC[293] = (((- ALGEBRAIC[276] - ALGEBRAIC[278]) - ALGEBRAIC[246])+ALGEBRAIC[261]+ALGEBRAIC[269]) - ALGEBRAIC[274];
ALGEBRAIC[294] = (((ALGEBRAIC[276]+ 2.00000*ALGEBRAIC[278]+ALGEBRAIC[246]) - ALGEBRAIC[261]) - ALGEBRAIC[269])+ALGEBRAIC[274];
ALGEBRAIC[295] = - ALGEBRAIC[278];
ALGEBRAIC[296] = (- (ALGEBRAIC[226]+ALGEBRAIC[234]) - ALGEBRAIC[258])+ALGEBRAIC[274];
ALGEBRAIC[297] = - (ALGEBRAIC[226]+ALGEBRAIC[234]);
ALGEBRAIC[298] = (ALGEBRAIC[226]+ALGEBRAIC[234]) - ALGEBRAIC[237];
ALGEBRAIC[299] = ALGEBRAIC[237] - ALGEBRAIC[240];
ALGEBRAIC[300] = ALGEBRAIC[240] - ALGEBRAIC[246];
ALGEBRAIC[301] = ALGEBRAIC[246] - ALGEBRAIC[249];
ALGEBRAIC[302] = ALGEBRAIC[249]+ALGEBRAIC[251]+ALGEBRAIC[254];
ALGEBRAIC[303] = - ALGEBRAIC[254];
ALGEBRAIC[304] = (ALGEBRAIC[249] - ALGEBRAIC[251]) - ALGEBRAIC[258];
ALGEBRAIC[305] = ALGEBRAIC[258] - ALGEBRAIC[261];
ALGEBRAIC[306] = ALGEBRAIC[261] - ALGEBRAIC[264];
ALGEBRAIC[307] = ALGEBRAIC[264] - ALGEBRAIC[266];
ALGEBRAIC[308] = ALGEBRAIC[266] - ALGEBRAIC[269];
ALGEBRAIC[309] = ALGEBRAIC[269] - ALGEBRAIC[272];
ALGEBRAIC[310] = ALGEBRAIC[272];
ALGEBRAIC[5] = (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[150])*ALGEBRAIC[2];
ALGEBRAIC[48] = (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[187])*ALGEBRAIC[45];
ALGEBRAIC[62] = (( ALGEBRAIC[60]*STATES[1])/CONSTANTS[134])*pow(10.0000, CONSTANTS[196]);
ALGEBRAIC[85] = (( ALGEBRAIC[82]*STATES[1])/CONSTANTS[134])*pow(10.0000, CONSTANTS[207]);
ALGEBRAIC[101] = (( ALGEBRAIC[99]*STATES[1])/CONSTANTS[134])*pow(10.0000, CONSTANTS[213]);
ALGEBRAIC[111] = (( ALGEBRAIC[109]*STATES[1])/CONSTANTS[134])*pow(10.0000, CONSTANTS[220]);
ALGEBRAIC[133] = (( ALGEBRAIC[131]*STATES[1])/CONSTANTS[134])*pow(10.0000, CONSTANTS[230]);
ALGEBRAIC[145] = (( ALGEBRAIC[142]*STATES[1])/CONSTANTS[134])*pow(10.0000, CONSTANTS[235]);
ALGEBRAIC[161] = (STATES[1]/CONSTANTS[134])*pow(10.0000, CONSTANTS[242])*ALGEBRAIC[159];
ALGEBRAIC[311] = ( ALGEBRAIC[15]*ALGEBRAIC[293]+ ALGEBRAIC[26]*ALGEBRAIC[294]+ ALGEBRAIC[37]*ALGEBRAIC[295]+ ALGEBRAIC[5]*ALGEBRAIC[296]+ ALGEBRAIC[48]*ALGEBRAIC[289]+ ALGEBRAIC[62]*ALGEBRAIC[298]+ ALGEBRAIC[85]*ALGEBRAIC[301]+ ALGEBRAIC[101]*ALGEBRAIC[303]+ ALGEBRAIC[111]*ALGEBRAIC[302]+ ALGEBRAIC[145]*ALGEBRAIC[308]+ ALGEBRAIC[133]*ALGEBRAIC[307]+ ALGEBRAIC[161]*ALGEBRAIC[310])/ALGEBRAIC[286];
}
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;
SI[24] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - 1.00000;
CONDVAR[1] = VOI - 1.00000;
CONDVAR[2] = VOI - 40.0000;
CONDVAR[3] = VOI - 40.0000;
CONDVAR[4] = VOI - 80.0000;
CONDVAR[5] = VOI - 80.0000;
CONDVAR[6] = VOI - 100.000;
CONDVAR[7] = VOI - 100.000;
CONDVAR[8] = VOI - 110.000;
CONDVAR[9] = VOI - 30.0000;
CONDVAR[10] = VOI - 1.00000;
}