# Size of variable arrays: sizeAlgebraic = 285 sizeStates = 32 sizeConstants = 230 from math import * from numpy import * def createLegends(): legend_states = [""] * sizeStates legend_rates = [""] * sizeStates legend_algebraic = [""] * sizeAlgebraic legend_voi = "" legend_constants = [""] * sizeConstants legend_voi = "time in component environment (minute)" legend_algebraic[0] = "pH_cy in component environment (dimensionless)" legend_states[0] = "pH_calc in component differential_equations (dimensionless)" legend_constants[1] = "CO2_in in component environment (molar)" legend_constants[48] = "HCO3_in in component environment (molar)" legend_constants[55] = "pH_in in component environment (dimensionless)" legend_constants[2] = "MOPSt in component environment (molar)" legend_constants[0] = "clamp in component environment (dimensionless)" legend_constants[46] = "Vfgly in component parameters (molar_per_minute)" legend_constants[3] = "fracA in component parameters (dimensionless)" legend_constants[4] = "KgpA_glyf in component parameters (molar)" legend_constants[51] = "KgpA_pi in component parameters (molar)" legend_constants[57] = "KgpA_igly in component parameters (molar)" legend_constants[5] = "KgpA_ipi in component parameters (molar)" legend_constants[58] = "KgpA_glyb in component parameters (molar)" legend_constants[6] = "KgpA_g1p in component parameters (molar)" legend_constants[59] = "KgpA_ig1p in component parameters (molar)" legend_constants[60] = "KgpB_pi in component parameters (molar)" legend_constants[7] = "KgpB_ipi in component parameters (molar)" legend_constants[61] = "KgpB_iglyf in component parameters (molar)" legend_constants[62] = "KgpB_g1p in component parameters (molar)" legend_constants[8] = "KgpB_ig1p in component parameters (molar)" legend_constants[63] = "KgpB_iglyb in component parameters (molar)" legend_constants[9] = "Kgp_amp in component parameters (molar)" legend_constants[10] = "nH in component parameters (dimensionless)" legend_constants[64] = "Vffpglm in component parameters (molar_per_minute)" legend_constants[65] = "Kpglm_g1p in component parameters (molar)" legend_constants[66] = "Kpglm_g6p in component parameters (molar)" legend_constants[67] = "Vbbpgi in component parameters (molar_per_minute)" legend_constants[68] = "Kpgi_g6p in component parameters (molar)" legend_constants[69] = "Kpgi_f6p in component parameters (molar)" legend_constants[70] = "Vffpfk in component parameters (molar_per_minute)" legend_constants[71] = "Kpfk_f6p in component parameters (molar)" legend_constants[72] = "Kpfk_f6pT in component parameters (molar)" legend_constants[73] = "Kpfk_atp in component parameters (molar)" legend_constants[74] = "Kpfk_atpT in component parameters (molar)" legend_constants[75] = "Kpfk_fbp in component parameters (molar)" legend_constants[11] = "Kpfk_fbpT in component parameters (molar)" legend_constants[77] = "Kpfk_adp in component parameters (molar)" legend_constants[12] = "Kpfk_adpT in component parameters (molar)" legend_constants[13] = "Kpfki in component parameters (molar)" legend_constants[14] = "Kmpfk in component parameters (molar)" legend_constants[15] = "d in component parameters (dimensionless)" legend_constants[16] = "e_ in component parameters (dimensionless)" legend_constants[17] = "Lo in component parameters (dimensionless)" legend_constants[78] = "Vffald in component parameters (molar_per_minute)" legend_constants[79] = "Kald_fbp in component parameters (molar)" legend_constants[80] = "Kald_dhap in component parameters (molar)" legend_constants[81] = "Kald_gap in component parameters (molar)" legend_constants[82] = "Vfftpi in component parameters (molar_per_minute)" legend_constants[83] = "Ktpi_gap in component parameters (molar)" legend_constants[84] = "Ktpi_dhap in component parameters (molar)" legend_constants[85] = "Vbbg3pdh in component parameters (molar_per_minute)" legend_constants[86] = "Kg3pdh_g3p in component parameters (molar)" legend_constants[87] = "Kg3pdh_nad in component parameters (molar)" legend_constants[88] = "Kg3pdh_dhap in component parameters (molar)" legend_constants[89] = "Kg3pdh_nadh in component parameters (molar)" legend_constants[90] = "Vffgad in component parameters (molar_per_minute)" legend_constants[91] = "Kgapdh_gap in component parameters (molar)" legend_constants[92] = "Kgapdh_nad in component parameters (molar)" legend_constants[93] = "Kgapdh_pi in component parameters (molar)" legend_constants[94] = "Kgapdh_bpg in component parameters (molar)" legend_constants[95] = "Kgapdh_nadh in component parameters (molar)" legend_constants[96] = "Vbbpgk in component parameters (molar_per_minute)" legend_constants[97] = "Kpgk_bpg in component parameters (molar)" legend_constants[98] = "Kpgk_adp in component parameters (molar)" legend_constants[99] = "Kpgk_3pg in component parameters (molar)" legend_constants[100] = "Kpgk_atp in component parameters (molar)" legend_constants[101] = "Vffpgm in component parameters (molar_per_minute)" legend_constants[102] = "Kpgm_3pg in component parameters (molar)" legend_constants[103] = "Kpgm_2pg in component parameters (molar)" legend_constants[104] = "Vffen in component parameters (molar_per_minute)" legend_constants[105] = "Ken_2pg in component parameters (molar)" legend_constants[106] = "Ken_pep in component parameters (molar)" legend_constants[107] = "Vffpk in component parameters (molar_per_minute)" legend_constants[108] = "Kpk_pep in component parameters (molar)" legend_constants[109] = "Kpk_adp in component parameters (molar)" legend_constants[110] = "Kpk_pyr in component parameters (molar)" legend_constants[111] = "Kpk_atp in component parameters (molar)" legend_constants[112] = "Vffldh in component parameters (molar_per_minute)" legend_constants[113] = "Kldh_pyr in component parameters (molar)" legend_constants[114] = "Kldh_nadh in component parameters (molar)" legend_constants[115] = "Kldh_lac in component parameters (molar)" legend_constants[116] = "Kldh_nad in component parameters (molar)" legend_constants[117] = "VforCK in component parameters (molar_per_minute)" legend_constants[118] = "Kck_pcr in component parameters (molar)" legend_constants[119] = "Kck_iatp in component parameters (molar)" legend_constants[120] = "Kck_iadp in component parameters (molar)" legend_constants[18] = "Kck_ipcr in component parameters (molar)" legend_constants[121] = "Kck_cr in component parameters (molar)" legend_constants[19] = "Kck_Iadp in component parameters (molar)" legend_constants[20] = "Kck_Ipcr in component parameters (molar)" legend_constants[21] = "Kck_Icr in component parameters (molar)" legend_constants[22] = "Kck_Iatp in component parameters (molar)" legend_constants[23] = "Kck_Ii in component parameters (molar)" legend_constants[122] = "Vfadk in component parameters (molar_per_minute)" legend_constants[123] = "Kadk_amp in component parameters (molar)" legend_constants[124] = "Kadk_atp in component parameters (molar)" legend_constants[125] = "Kadk_adp in component parameters (molar)" legend_constants[24] = "kut in component parameters (per_minute)" legend_constants[25] = "VmaxOxPhos in component parameters (molar_per_minute)" legend_constants[26] = "Kmadp in component parameters (molar)" legend_constants[27] = "nnHH in component parameters (dimensionless)" legend_constants[28] = "Vmax_MCT in component parameters (molar_per_minute)" legend_constants[29] = "Klac1 in component parameters (molar)" legend_constants[30] = "pKmct1 in component parameters (dimensionless)" legend_constants[31] = "PSVCO2 in component parameters (per_minute)" legend_constants[32] = "X_CA in component parameters (dimensionless)" legend_constants[33] = "kf_hydrCO2 in component parameters (per_minute)" legend_constants[49] = "Keq_CA in component parameters (dimensionless)" legend_constants[56] = "kb_hydrCO2 in component parameters (per_minute)" legend_constants[34] = "Flow in component parameters (liter_per_minute)" legend_constants[35] = "V_cell in component parameters (liter)" legend_constants[36] = "V_ec in component parameters (liter)" legend_constants[126] = "bufcapfixed in component parameters (molar)" legend_constants[127] = "T in component equilibrium_constants (kelvin)" legend_constants[128] = "I in component equilibrium_constants (molar)" legend_constants[129] = "k in component equilibrium_constants (molar)" legend_constants[130] = "c0 in component equilibrium_constants (molar)" legend_constants[131] = "R in component equilibrium_constants (kilojoule_per_mole_per_kelvin)" legend_constants[132] = "B in component equilibrium_constants (per_half_molar)" legend_constants[133] = "I1 in component equilibrium_constants (molar)" legend_constants[134] = "alphadebye in component equilibrium_constants (per_half_molar)" legend_constants[135] = "IcorrpKa in component equilibrium_constants (dimensionless)" legend_constants[137] = "RTalpha in component equilibrium_constants (kilojoule_half_liter_per_3_half_mole)" legend_constants[143] = "IcorrdeltaGpof in component equilibrium_constants (kilojoule_per_mole)" legend_constants[144] = "pKa1_Pi in component equilibrium_constants (dimensionless)" legend_constants[136] = "pKamg_Pi in component equilibrium_constants (dimensionless)" legend_constants[145] = "pKak_Pi in component equilibrium_constants (dimensionless)" legend_algebraic[1] = "P_Pi in component equilibrium_constants (dimensionless)" legend_algebraic[2] = "HPi2 in component equilibrium_constants (dimensionless)" legend_algebraic[3] = "H2Pi1 in component equilibrium_constants (dimensionless)" legend_algebraic[4] = "kPi in component equilibrium_constants (dimensionless)" legend_algebraic[5] = "mgPi in component equilibrium_constants (dimensionless)" legend_algebraic[6] = "Navg_Pi in component equilibrium_constants (dimensionless)" legend_algebraic[7] = "dNavgPidH in component equilibrium_constants (per_molar)" legend_constants[146] = "NH_HPi2 in component equilibrium_constants (dimensionless)" legend_constants[147] = "deltaGof_HPi2 in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[8] = "deltaGpof_HPi2 in component equilibrium_constants (kilojoule_per_mole)" legend_constants[148] = "pKa1_ATP in component equilibrium_constants (dimensionless)" legend_constants[138] = "pKamg_ATP in component equilibrium_constants (dimensionless)" legend_constants[149] = "pKak_ATP in component equilibrium_constants (dimensionless)" legend_algebraic[9] = "P_ATP in component equilibrium_constants (dimensionless)" legend_algebraic[10] = "ATP4 in component equilibrium_constants (dimensionless)" legend_algebraic[13] = "HATP3 in component equilibrium_constants (dimensionless)" legend_algebraic[11] = "mgATP2 in component equilibrium_constants (dimensionless)" legend_algebraic[12] = "kATP in component equilibrium_constants (dimensionless)" legend_algebraic[14] = "Navg_ATP in component equilibrium_constants (dimensionless)" legend_algebraic[15] = "dNavgATPdH in component equilibrium_constants (per_molar)" legend_constants[150] = "NH_ATP4 in component equilibrium_constants (dimensionless)" legend_constants[151] = "deltaGof_ATP4 in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[16] = "deltaGpof_ATP4 in component equilibrium_constants (kilojoule_per_mole)" legend_constants[152] = "pKak_ADP in component equilibrium_constants (dimensionless)" legend_constants[153] = "pKa1_ADP in component equilibrium_constants (dimensionless)" legend_constants[139] = "pKamg_ADP in component equilibrium_constants (dimensionless)" legend_algebraic[17] = "P_ADP in component equilibrium_constants (dimensionless)" legend_algebraic[18] = "ADP3 in component equilibrium_constants (dimensionless)" legend_algebraic[21] = "HADP2 in component equilibrium_constants (dimensionless)" legend_algebraic[19] = "mgADP in component equilibrium_constants (dimensionless)" legend_algebraic[20] = "kADP in component equilibrium_constants (dimensionless)" legend_algebraic[22] = "Navg_ADP in component equilibrium_constants (dimensionless)" legend_algebraic[23] = "dNavgADPdH in component equilibrium_constants (per_molar)" legend_constants[154] = "NH_ADP3 in component equilibrium_constants (dimensionless)" legend_constants[155] = "deltaGof_ADP3 in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[24] = "deltaGpof_ADP3 in component equilibrium_constants (kilojoule_per_mole)" legend_constants[156] = "pKa1_AMP in component equilibrium_constants (dimensionless)" legend_constants[140] = "pKamg_AMP in component equilibrium_constants (dimensionless)" legend_algebraic[25] = "P_AMP in component equilibrium_constants (dimensionless)" legend_algebraic[26] = "AMP2 in component equilibrium_constants (dimensionless)" legend_algebraic[28] = "HAMP1 in component equilibrium_constants (dimensionless)" legend_algebraic[27] = "mgAMP in component equilibrium_constants (dimensionless)" legend_algebraic[29] = "Navg_AMP in component equilibrium_constants (dimensionless)" legend_algebraic[30] = "dNavgAMPdH in component equilibrium_constants (per_molar)" legend_constants[157] = "NH_AMP2 in component equilibrium_constants (dimensionless)" legend_constants[158] = "deltaGof_AMP2 in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[31] = "deltaGpof_AMP2 in component equilibrium_constants (kilojoule_per_mole)" legend_constants[159] = "pKak_PCR in component equilibrium_constants (dimensionless)" legend_constants[160] = "pKa1_PCR in component equilibrium_constants (dimensionless)" legend_constants[141] = "pKamg_PCR in component equilibrium_constants (dimensionless)" legend_algebraic[32] = "P_PCR in component equilibrium_constants (dimensionless)" legend_algebraic[33] = "HPCR in component equilibrium_constants (dimensionless)" legend_algebraic[36] = "H2PCR in component equilibrium_constants (dimensionless)" legend_algebraic[34] = "kPCR in component equilibrium_constants (dimensionless)" legend_algebraic[35] = "mgPCR in component equilibrium_constants (dimensionless)" legend_algebraic[37] = "Navg_PCR in component equilibrium_constants (dimensionless)" legend_algebraic[38] = "dNavgPCRdH in component equilibrium_constants (per_molar)" legend_constants[161] = "NH_HPCR in component equilibrium_constants (dimensionless)" legend_constants[162] = "pKa1_CR in component equilibrium_constants (dimensionless)" legend_algebraic[39] = "P_CR in component equilibrium_constants (dimensionless)" legend_algebraic[40] = "HCR in component equilibrium_constants (dimensionless)" legend_algebraic[41] = "H2CR in component equilibrium_constants (dimensionless)" legend_algebraic[42] = "Navg_CR in component equilibrium_constants (dimensionless)" legend_algebraic[43] = "dNavgCRdH in component equilibrium_constants (per_molar)" legend_constants[163] = "NH_HCR in component equilibrium_constants (dimensionless)" legend_constants[164] = "pKa1_G1P in component equilibrium_constants (dimensionless)" legend_constants[142] = "pKamg_G1P in component equilibrium_constants (dimensionless)" legend_algebraic[44] = "P_G1P in component equilibrium_constants (dimensionless)" legend_algebraic[45] = "UG1P in component equilibrium_constants (dimensionless)" legend_algebraic[46] = "HG1P in component equilibrium_constants (dimensionless)" legend_algebraic[47] = "mgG1P in component equilibrium_constants (dimensionless)" legend_algebraic[48] = "Navg_G1P in component equilibrium_constants (dimensionless)" legend_algebraic[49] = "dNavgG1PdH in component equilibrium_constants (per_molar)" legend_constants[165] = "NH_UG1P in component equilibrium_constants (dimensionless)" legend_constants[166] = "deltaGof_UG1P in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[50] = "deltaGpof_UG1P in component equilibrium_constants (kilojoule_per_mole)" legend_constants[167] = "pKa1_G6P in component equilibrium_constants (dimensionless)" legend_algebraic[51] = "P_G6P in component equilibrium_constants (dimensionless)" legend_algebraic[52] = "UG6P in component equilibrium_constants (dimensionless)" legend_algebraic[53] = "HG6P in component equilibrium_constants (dimensionless)" legend_algebraic[54] = "Navg_G6P in component equilibrium_constants (dimensionless)" legend_algebraic[55] = "dNavgG6PdH in component equilibrium_constants (per_molar)" legend_constants[37] = "dNavgG6Pdmg in component equilibrium_constants (per_molar)" legend_constants[168] = "NH_UG6P in component equilibrium_constants (dimensionless)" legend_constants[169] = "deltaGof_UG6P in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[56] = "deltaGpof_UG6P in component equilibrium_constants (kilojoule_per_mole)" legend_constants[170] = "pKa1_F6P in component equilibrium_constants (dimensionless)" legend_algebraic[57] = "P_F6P in component equilibrium_constants (dimensionless)" legend_algebraic[58] = "UF6P in component equilibrium_constants (dimensionless)" legend_algebraic[59] = "HF6P in component equilibrium_constants (dimensionless)" legend_algebraic[60] = "Navg_F6P in component equilibrium_constants (dimensionless)" legend_algebraic[61] = "dNavgF6PdH in component equilibrium_constants (per_molar)" legend_constants[171] = "NH_UF6P in component equilibrium_constants (dimensionless)" legend_constants[172] = "deltaGof_UF6P in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[62] = "deltaGpof_UF6P in component equilibrium_constants (kilojoule_per_mole)" legend_constants[173] = "pKa1_FDP in component equilibrium_constants (dimensionless)" legend_constants[38] = "pKamg_FDP in component equilibrium_constants (dimensionless)" legend_algebraic[63] = "P_FDP in component equilibrium_constants (dimensionless)" legend_algebraic[64] = "UFDP in component equilibrium_constants (dimensionless)" legend_algebraic[66] = "HFDP in component equilibrium_constants (dimensionless)" legend_algebraic[65] = "mgFDP in component equilibrium_constants (dimensionless)" legend_algebraic[67] = "Navg_FDP in component equilibrium_constants (dimensionless)" legend_algebraic[68] = "dNavgFDPdH in component equilibrium_constants (per_molar)" legend_constants[174] = "NH_UFDP in component equilibrium_constants (dimensionless)" legend_constants[175] = "deltaGof_UFDP in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[69] = "deltaGpof_UFDP in component equilibrium_constants (kilojoule_per_mole)" legend_constants[176] = "pKa1_GAP in component equilibrium_constants (dimensionless)" legend_algebraic[70] = "P_GAP in component equilibrium_constants (dimensionless)" legend_algebraic[71] = "UGAP in component equilibrium_constants (dimensionless)" legend_algebraic[72] = "HGAP in component equilibrium_constants (dimensionless)" legend_algebraic[73] = "Navg_GAP in component equilibrium_constants (dimensionless)" legend_algebraic[74] = "dNavgGAPdH in component equilibrium_constants (per_molar)" legend_constants[39] = "dNavgGAPdmg in component equilibrium_constants (per_molar)" legend_constants[177] = "NH_UGAP in component equilibrium_constants (dimensionless)" legend_constants[178] = "deltaGof_UGAP in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[75] = "deltaGpof_UGAP in component equilibrium_constants (kilojoule_per_mole)" legend_constants[40] = "pKamg_G3P in component equilibrium_constants (dimensionless)" legend_constants[179] = "pKa1_G3P in component equilibrium_constants (dimensionless)" legend_algebraic[76] = "P_G3P in component equilibrium_constants (dimensionless)" legend_algebraic[77] = "UG3P in component equilibrium_constants (dimensionless)" legend_algebraic[79] = "HG3P in component equilibrium_constants (dimensionless)" legend_algebraic[78] = "mgG3P in component equilibrium_constants (dimensionless)" legend_algebraic[80] = "Navg_G3P in component equilibrium_constants (dimensionless)" legend_algebraic[81] = "dNavgG3PdH in component equilibrium_constants (per_molar)" legend_constants[180] = "NH_UG3P in component equilibrium_constants (dimensionless)" legend_constants[181] = "deltaGof_UG3P in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[82] = "deltaGpof_UG3P in component equilibrium_constants (kilojoule_per_mole)" legend_constants[182] = "pKa1_DHAP in component equilibrium_constants (dimensionless)" legend_constants[41] = "pKamg_DHAP in component equilibrium_constants (dimensionless)" legend_algebraic[83] = "P_DHAP in component equilibrium_constants (dimensionless)" legend_algebraic[84] = "UDHAP in component equilibrium_constants (dimensionless)" legend_algebraic[86] = "HDHAP in component equilibrium_constants (dimensionless)" legend_algebraic[85] = "mgDHAP in component equilibrium_constants (dimensionless)" legend_algebraic[87] = "Navg_DHAP in component equilibrium_constants (dimensionless)" legend_algebraic[88] = "dNavgDHAPdH in component equilibrium_constants (per_molar)" legend_constants[183] = "NH_UDHAP in component equilibrium_constants (dimensionless)" legend_constants[184] = "deltaGof_UDHAP in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[89] = "deltaGpof_UDHAP in component equilibrium_constants (kilojoule_per_mole)" legend_constants[185] = "pKa1_13DPG in component equilibrium_constants (dimensionless)" legend_algebraic[90] = "P_13DPG in component equilibrium_constants (dimensionless)" legend_algebraic[91] = "U13DPG in component equilibrium_constants (dimensionless)" legend_algebraic[92] = "H13DPG in component equilibrium_constants (dimensionless)" legend_algebraic[93] = "Navg_13DPG in component equilibrium_constants (dimensionless)" legend_algebraic[94] = "dNavg13DPGdH in component equilibrium_constants (per_molar)" legend_constants[186] = "NH_U13DPG in component equilibrium_constants (dimensionless)" legend_constants[187] = "deltaGof_U13DPG in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[95] = "deltaGpof_U13DPG in component equilibrium_constants (kilojoule_per_mole)" legend_constants[188] = "pKa1_3PG in component equilibrium_constants (dimensionless)" legend_algebraic[96] = "P_3PG in component equilibrium_constants (dimensionless)" legend_algebraic[97] = "U3PG in component equilibrium_constants (dimensionless)" legend_algebraic[98] = "H3PG in component equilibrium_constants (dimensionless)" legend_algebraic[99] = "Navg_3PG in component equilibrium_constants (dimensionless)" legend_algebraic[100] = "dNavg3PGdH in component equilibrium_constants (per_molar)" legend_constants[189] = "NH_U3PG in component equilibrium_constants (dimensionless)" legend_constants[190] = "deltaGof_U3PG in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[101] = "deltaGpof_U3PG in component equilibrium_constants (kilojoule_per_mole)" legend_constants[191] = "pKa1_2PG in component equilibrium_constants (dimensionless)" legend_constants[42] = "pKamg_2PG in component equilibrium_constants (dimensionless)" legend_constants[192] = "pKak_2PG in component equilibrium_constants (dimensionless)" legend_algebraic[102] = "P_2PG in component equilibrium_constants (dimensionless)" legend_algebraic[103] = "U2PG in component equilibrium_constants (dimensionless)" legend_algebraic[104] = "H2PG in component equilibrium_constants (dimensionless)" legend_algebraic[105] = "mg2PG in component equilibrium_constants (dimensionless)" legend_algebraic[106] = "k2PG in component equilibrium_constants (dimensionless)" legend_algebraic[107] = "Navg_2PG in component equilibrium_constants (dimensionless)" legend_algebraic[108] = "dNavg2PGdH in component equilibrium_constants (per_molar)" legend_constants[193] = "NH_U2PG in component equilibrium_constants (dimensionless)" legend_constants[194] = "deltaGof_U2PG in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[109] = "deltaGpof_U2PG in component equilibrium_constants (kilojoule_per_mole)" legend_constants[195] = "pKa1_PEP in component equilibrium_constants (dimensionless)" legend_constants[43] = "pKamg_PEP in component equilibrium_constants (dimensionless)" legend_constants[196] = "pKak_PEP in component equilibrium_constants (dimensionless)" legend_algebraic[110] = "P_PEP in component equilibrium_constants (dimensionless)" legend_algebraic[111] = "UPEP in component equilibrium_constants (dimensionless)" legend_algebraic[112] = "HPEP in component equilibrium_constants (dimensionless)" legend_algebraic[113] = "kPEP in component equilibrium_constants (dimensionless)" legend_algebraic[114] = "mgPEP in component equilibrium_constants (dimensionless)" legend_algebraic[115] = "Navg_PEP in component equilibrium_constants (dimensionless)" legend_algebraic[116] = "dNavgPEPdH in component equilibrium_constants (per_molar)" legend_constants[197] = "NH_UPEP in component equilibrium_constants (dimensionless)" legend_constants[198] = "deltaGof_UPEP in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[117] = "deltaGpof_UPEP in component equilibrium_constants (kilojoule_per_mole)" legend_constants[199] = "pKa1_PYR in component equilibrium_constants (dimensionless)" legend_algebraic[118] = "P_PYR in component equilibrium_constants (dimensionless)" legend_algebraic[119] = "UPYR in component equilibrium_constants (dimensionless)" legend_algebraic[120] = "HPYR in component equilibrium_constants (dimensionless)" legend_algebraic[121] = "Navg_PYR in component equilibrium_constants (dimensionless)" legend_algebraic[122] = "dNavgPYRdH in component equilibrium_constants (per_molar)" legend_constants[44] = "dNavgPYRdmg in component equilibrium_constants (per_molar)" legend_constants[200] = "NH_UPYR in component equilibrium_constants (dimensionless)" legend_constants[201] = "deltaGof_UPYR in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[123] = "deltaGpof_UPYR in component equilibrium_constants (kilojoule_per_mole)" legend_constants[45] = "pKamg_LAC in component equilibrium_constants (dimensionless)" legend_constants[202] = "pKa1_LAC in component equilibrium_constants (dimensionless)" legend_algebraic[124] = "P_LAC in component equilibrium_constants (dimensionless)" legend_algebraic[127] = "ULAC in component equilibrium_constants (dimensionless)" legend_algebraic[129] = "HLAC in component equilibrium_constants (dimensionless)" legend_algebraic[128] = "mgLAC in component equilibrium_constants (dimensionless)" legend_algebraic[130] = "Navg_LAC in component equilibrium_constants (dimensionless)" legend_algebraic[131] = "dNavgLACdH in component equilibrium_constants (per_molar)" legend_algebraic[125] = "dmgLACdmg in component equilibrium_constants (per_molar)" legend_algebraic[126] = "dmgLACdpH in component equilibrium_constants (dimensionless)" legend_constants[203] = "NH_ULAC in component equilibrium_constants (dimensionless)" legend_constants[204] = "deltaGof_ULAC in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[132] = "deltaGpof_ULAC in component equilibrium_constants (kilojoule_per_mole)" legend_constants[205] = "dNH_GLY in component equilibrium_constants (dimensionless)" legend_algebraic[133] = "deltaGpo_GLY in component equilibrium_constants (kilojoule_per_mole)" legend_constants[206] = "NH_NAD in component equilibrium_constants (dimensionless)" legend_constants[207] = "deltaGof_NAD in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[134] = "deltaGpof_NAD in component equilibrium_constants (kilojoule_per_mole)" legend_constants[208] = "NH_NADH in component equilibrium_constants (dimensionless)" legend_constants[209] = "deltaGof_NADH in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[135] = "deltaGpof_NADH in component equilibrium_constants (kilojoule_per_mole)" legend_constants[210] = "NH_H2O in component equilibrium_constants (dimensionless)" legend_constants[211] = "deltaGof_H2O in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[136] = "deltaGpof_H2O in component equilibrium_constants (kilojoule_per_mole)" legend_constants[212] = "NH_H in component equilibrium_constants (dimensionless)" legend_constants[213] = "deltaGof_H in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[137] = "deltaGpof_H in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[138] = "deltaH_CK in component equilibrium_constants (dimensionless)" legend_constants[214] = "slopeint_1 in component equilibrium_constants (kelvin)" legend_constants[215] = "slopeint_2 in component equilibrium_constants (dimensionless)" legend_constants[216] = "Kref_CKT in component equilibrium_constants (dimensionless)" legend_constants[217] = "I25 in component equilibrium_constants (molar)" legend_constants[218] = "alphadebye1 in component equilibrium_constants (dimensionless)" legend_constants[219] = "IcorrneglnKrefCK in component equilibrium_constants (dimensionless)" legend_constants[220] = "neglnKrefCK_Iexp in component equilibrium_constants (dimensionless)" legend_constants[221] = "Kref_CK in component equilibrium_constants (dimensionless)" legend_constants[222] = "deltaGpo_CK in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[140] = "Kapp_CK in component equilibrium_constants (dimensionless)" legend_algebraic[142] = "deltaH_ADK in component equilibrium_constants (dimensionless)" legend_algebraic[143] = "deltaGpo_ADK in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[144] = "Kapp_ADK in component equilibrium_constants (dimensionless)" legend_algebraic[145] = "deltaH_GP in component equilibrium_constants (dimensionless)" legend_algebraic[146] = "deltaGpo_GP in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[147] = "Kapp_GP in component equilibrium_constants (dimensionless)" legend_algebraic[148] = "deltaH_PGLM in component equilibrium_constants (dimensionless)" legend_algebraic[149] = "deltaGpo_PGLM in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[150] = "Kapp_PGLM in component equilibrium_constants (dimensionless)" legend_algebraic[151] = "deltaH_PGI in component equilibrium_constants (dimensionless)" legend_algebraic[152] = "deltaGpo_PGI in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[153] = "Kapp_PGI in component equilibrium_constants (dimensionless)" legend_algebraic[154] = "deltaH_PFK in component equilibrium_constants (dimensionless)" legend_algebraic[155] = "deltaGpo_PFK in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[156] = "Kapp_PFK in component equilibrium_constants (dimensionless)" legend_algebraic[157] = "deltaH_ALD in component equilibrium_constants (dimensionless)" legend_algebraic[158] = "deltaGpo_ALD in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[159] = "Kapp_ALD in component equilibrium_constants (dimensionless)" legend_algebraic[160] = "deltaH_TPI in component equilibrium_constants (dimensionless)" legend_algebraic[161] = "deltaGpo_TPI in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[162] = "Kapp_TPI in component equilibrium_constants (dimensionless)" legend_algebraic[163] = "deltaH_GAPDH in component equilibrium_constants (dimensionless)" legend_algebraic[164] = "deltaGpo_GAPDH in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[165] = "Kapp_GAPDH in component equilibrium_constants (dimensionless)" legend_algebraic[166] = "deltaH_G3PDH in component equilibrium_constants (dimensionless)" legend_algebraic[167] = "deltaGpo_G3PDH in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[168] = "Kapp_G3PDH in component equilibrium_constants (dimensionless)" legend_algebraic[169] = "deltaH_PGK in component equilibrium_constants (dimensionless)" legend_algebraic[170] = "deltaGpo_PGK in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[171] = "Kapp_PGK in component equilibrium_constants (dimensionless)" legend_algebraic[172] = "deltaH_PGM in component equilibrium_constants (dimensionless)" legend_algebraic[173] = "deltaGpo_PGM in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[174] = "Kapp_PGM in component equilibrium_constants (dimensionless)" legend_algebraic[175] = "deltaH_ENOL in component equilibrium_constants (dimensionless)" legend_algebraic[176] = "deltaGpo_ENOL in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[177] = "Kapp_ENOL in component equilibrium_constants (dimensionless)" legend_algebraic[178] = "deltaH_PK in component equilibrium_constants (dimensionless)" legend_algebraic[179] = "deltaGpo_PK in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[180] = "Kapp_PK in component equilibrium_constants (dimensionless)" legend_algebraic[181] = "deltaH_LDH in component equilibrium_constants (dimensionless)" legend_algebraic[182] = "deltaGpo_LDH in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[183] = "Kapp_LDH in component equilibrium_constants (dimensionless)" legend_algebraic[184] = "deltaH_ATPase in component equilibrium_constants (dimensionless)" legend_algebraic[139] = "deltaGpo_ATPase in component equilibrium_constants (kilojoule_per_mole)" legend_algebraic[141] = "Kapp_ATPase in component equilibrium_constants (dimensionless)" legend_states[1] = "mg in component differential_equations (molar)" legend_algebraic[185] = "Dglya in component glycogen_phosphorylase_A (dimensionless)" legend_algebraic[186] = "pa in component glycogen_phosphorylase_A (dimensionless)" legend_algebraic[187] = "VbglyA in component glycogen_phosphorylase_A (molar_per_minute)" legend_algebraic[188] = "glyAF in component glycogen_phosphorylase_A (per_minute)" legend_algebraic[189] = "glyAR in component glycogen_phosphorylase_A (per_minute)" legend_algebraic[190] = "flux_GPa in component glycogen_phosphorylase_A (molar_per_minute)" legend_states[2] = "Gly in component differential_equations (molar)" legend_states[3] = "Pi in component differential_equations (molar)" legend_states[4] = "G1P in component differential_equations (molar)" legend_constants[50] = "fracB in component glycogen_phosphorylase_B (dimensionless)" legend_algebraic[191] = "M in component glycogen_phosphorylase_B (dimensionless)" legend_algebraic[192] = "Dglyb in component glycogen_phosphorylase_B (dimensionless)" legend_algebraic[193] = "pb in component glycogen_phosphorylase_B (dimensionless)" legend_algebraic[194] = "VbglyB in component glycogen_phosphorylase_B (molar_per_minute)" legend_algebraic[195] = "glyBF in component glycogen_phosphorylase_B (per_minute)" legend_algebraic[196] = "glyBR in component glycogen_phosphorylase_B (per_minute)" legend_algebraic[197] = "flux_GPb in component glycogen_phosphorylase_B (molar_per_minute)" legend_states[5] = "AMP in component differential_equations (molar)" legend_algebraic[199] = "Vfpglm in component PGLM (molar_per_minute)" legend_algebraic[200] = "Vbpglm in component PGLM (molar_per_minute)" legend_algebraic[201] = "v_PGLM in component PGLM (molar_per_minute)" legend_states[6] = "G6P in component differential_equations (molar)" legend_algebraic[203] = "Vbpgi in component PGI (molar_per_minute)" legend_algebraic[204] = "Vfpgi in component PGI (molar_per_minute)" legend_algebraic[205] = "v_PGI in component PGI (molar_per_minute)" legend_states[7] = "F6P in component differential_equations (molar)" legend_algebraic[207] = "Vfpfk in component PFK (molar_per_minute)" legend_algebraic[208] = "Vbpfk in component PFK (molar_per_minute)" legend_algebraic[209] = "L in component PFK (dimensionless)" legend_constants[76] = "alpha in component PFK (dimensionless)" legend_algebraic[210] = "Delta in component PFK (dimensionless)" legend_algebraic[211] = "Deltap in component PFK (dimensionless)" legend_algebraic[212] = "v_PFK in component PFK (molar_per_minute)" legend_states[8] = "ATP in component differential_equations (molar)" legend_states[9] = "ADP in component differential_equations (molar)" legend_states[10] = "FBP in component differential_equations (molar)" legend_algebraic[214] = "Vfald in component ALD (molar_per_minute)" legend_algebraic[215] = "Vbald in component ALD (molar_per_minute)" legend_algebraic[216] = "v_ALD in component ALD (molar_per_minute)" legend_states[11] = "GAP in component differential_equations (molar)" legend_states[12] = "DHAP in component differential_equations (molar)" legend_constants[223] = "Vftpi in component TPI (molar_per_minute)" legend_algebraic[218] = "Vbtpi in component TPI (molar_per_minute)" legend_algebraic[219] = "v_TPI in component TPI (molar_per_minute)" legend_algebraic[220] = "Dg3pdh in component G3PDH (dimensionless)" legend_constants[224] = "Vbg3pdh in component G3PDH (molar_per_minute)" legend_algebraic[221] = "Vfg3pdh in component G3PDH (molar_per_minute)" legend_algebraic[222] = "v_G3PDH in component G3PDH (molar_per_minute)" legend_states[13] = "G3P in component differential_equations (molar)" legend_states[14] = "NADH in component differential_equations (molar)" legend_states[15] = "NAD in component differential_equations (molar)" legend_algebraic[225] = "Dgap in component GAPDH (dimensionless)" legend_algebraic[226] = "Vfgad in component GAPDH (molar_per_minute)" legend_algebraic[227] = "Vbgad in component GAPDH (molar_per_minute)" legend_algebraic[228] = "v_GAPDH in component GAPDH (molar_per_minute)" legend_states[16] = "BPG in component differential_equations (molar)" legend_constants[225] = "Vbpgk in component PGK (molar_per_minute)" legend_algebraic[230] = "Vfpgk in component PGK (molar_per_minute)" legend_algebraic[231] = "D_PGK in component PGK (dimensionless)" legend_algebraic[232] = "v_PGK in component PGK (molar_per_minute)" legend_states[17] = "P3G in component differential_equations (molar)" legend_algebraic[234] = "Vfpgm in component PGM (molar_per_minute)" legend_algebraic[235] = "Vbpgm in component PGM (molar_per_minute)" legend_algebraic[236] = "v_PGM in component PGM (molar_per_minute)" legend_states[18] = "P2G in component differential_equations (molar)" legend_constants[226] = "Vfen in component ENOL (molar_per_minute)" legend_algebraic[238] = "Vben in component ENOL (molar_per_minute)" legend_algebraic[239] = "v_ENOL in component ENOL (molar_per_minute)" legend_states[19] = "PEP in component differential_equations (molar)" legend_algebraic[241] = "Vfpk in component PK (molar_per_minute)" legend_algebraic[242] = "Vbpk in component PK (molar_per_minute)" legend_algebraic[243] = "v_PK in component PK (molar_per_minute)" legend_states[20] = "PYR in component differential_equations (molar)" legend_algebraic[245] = "Vfldh in component LDH (molar_per_minute)" legend_algebraic[246] = "Vbldh in component LDH (molar_per_minute)" legend_algebraic[247] = "v_LDH in component LDH (molar_per_minute)" legend_states[21] = "LAC in component differential_equations (molar)" legend_algebraic[251] = "VrevCK in component creatine_kinase (molar_per_minute)" legend_algebraic[252] = "CK in component creatine_kinase (molar_per_minute)" legend_states[22] = "Cr in component differential_equations (molar)" legend_states[23] = "PCr in component differential_equations (molar)" legend_algebraic[255] = "Vbadk in component adenylate_kinase (molar_per_minute)" legend_algebraic[256] = "ADK in component adenylate_kinase (molar_per_minute)" legend_algebraic[258] = "ATPase in component ATPase (molar_per_minute)" legend_algebraic[259] = "oxphos in component OXPHOS (molar_per_minute)" legend_constants[52] = "Kh in component lactate_transport (dimensionless)" legend_algebraic[260] = "v_MCT in component lactate_transport (molar_per_minute)" legend_states[24] = "pH_ec in component differential_equations (dimensionless)" legend_states[25] = "LAC_ec in component differential_equations (molar)" legend_constants[227] = "RQ in component CO2_generation_transport (dimensionless)" legend_constants[228] = "PO2ratio in component CO2_generation_transport (dimensionless)" legend_algebraic[264] = "fluxO2 in component CO2_generation_transport (molar_per_minute)" legend_algebraic[267] = "fluxCO2 in component CO2_generation_transport (molar_per_minute)" legend_algebraic[268] = "v_psCO2 in component Passive_permeation_flux (molar_per_minute)" legend_algebraic[269] = "v_hydrCO2cy in component Passive_permeation_flux (molar_per_minute)" legend_algebraic[270] = "v_hydrCO2ec in component Passive_permeation_flux (molar_per_minute)" legend_states[26] = "CO2_ec in component differential_equations (molar)" legend_states[27] = "HCO3_ec in component differential_equations (molar)" legend_states[28] = "CO2_cy in component differential_equations (molar)" legend_states[29] = "HCO3_cy in component differential_equations (molar)" legend_algebraic[271] = "bufcapmetab in component buffer_capacity (molar)" legend_algebraic[277] = "CKprtflux in component proton_flux (molar_per_minute)" legend_algebraic[281] = "glycprtflux in component proton_flux (molar_per_minute)" legend_algebraic[282] = "protons_consumed in component proton_flux (molar_per_minute)" legend_algebraic[283] = "pHODEterm1 in component proton_flux (per_minute)" legend_states[30] = "MOPS in component differential_equations (molar)" legend_states[31] = "MOPSH in component differential_equations (molar)" legend_algebraic[253] = "dPCrdt in component differential_equations (molar_per_minute)" legend_algebraic[254] = "dCrdt in component differential_equations (molar_per_minute)" legend_algebraic[248] = "dNADdt in component differential_equations (molar_per_minute)" legend_algebraic[249] = "dNADHdt in component differential_equations (molar_per_minute)" legend_algebraic[261] = "dATPdt in component differential_equations (molar_per_minute)" legend_algebraic[262] = "dADPdt in component differential_equations (molar_per_minute)" legend_algebraic[257] = "dAMPdt in component differential_equations (molar_per_minute)" legend_algebraic[263] = "dPidt in component differential_equations (molar_per_minute)" legend_algebraic[198] = "dglydt in component differential_equations (molar_per_minute)" legend_algebraic[202] = "dG1Pdt in component differential_equations (molar_per_minute)" legend_algebraic[206] = "dG6Pdt in component differential_equations (molar_per_minute)" legend_algebraic[213] = "dF6Pdt in component differential_equations (molar_per_minute)" legend_algebraic[217] = "dFBPdt in component differential_equations (molar_per_minute)" legend_algebraic[223] = "dDHAPdt in component differential_equations (molar_per_minute)" legend_algebraic[224] = "dG3Pdt in component differential_equations (molar_per_minute)" legend_algebraic[229] = "dGAPdt in component differential_equations (molar_per_minute)" legend_algebraic[233] = "dBPGdt in component differential_equations (molar_per_minute)" legend_algebraic[237] = "dP3Gdt in component differential_equations (molar_per_minute)" legend_algebraic[240] = "dP2Gdt in component differential_equations (molar_per_minute)" legend_algebraic[244] = "dPEPdt in component differential_equations (molar_per_minute)" legend_algebraic[250] = "dPYRdt in component differential_equations (molar_per_minute)" legend_algebraic[265] = "dLACdt in component differential_equations (molar_per_minute)" legend_constants[47] = "dmgdt in component differential_equations (molar_per_minute)" legend_algebraic[284] = "dpH_calcdt in component differential_equations (per_minute)" legend_algebraic[266] = "dLACecdt in component differential_equations (molar_per_minute)" legend_algebraic[272] = "dHCO3cydt in component differential_equations (molar_per_minute)" legend_algebraic[273] = "dCO2cydt in component differential_equations (molar_per_minute)" legend_algebraic[274] = "dHCO3ecdt in component differential_equations (molar_per_minute)" legend_algebraic[275] = "dCO2ecdt in component differential_equations (molar_per_minute)" legend_constants[53] = "MOPS_in in component differential_equations (molar)" legend_constants[54] = "MOPSH_in in component differential_equations (molar)" legend_algebraic[276] = "v_MOPSbuff in component differential_equations (molar_per_minute)" legend_algebraic[278] = "dMOPSdt in component differential_equations (molar_per_minute)" legend_algebraic[279] = "dMOPSHdt in component differential_equations (molar_per_minute)" legend_algebraic[280] = "dpHecdt in component differential_equations (per_minute)" legend_rates[23] = "d/dt PCr in component differential_equations (molar)" legend_rates[22] = "d/dt Cr in component differential_equations (molar)" legend_rates[15] = "d/dt NAD in component differential_equations (molar)" legend_rates[14] = "d/dt NADH in component differential_equations (molar)" legend_rates[8] = "d/dt ATP in component differential_equations (molar)" legend_rates[9] = "d/dt ADP in component differential_equations (molar)" legend_rates[5] = "d/dt AMP in component differential_equations (molar)" legend_rates[3] = "d/dt Pi in component differential_equations (molar)" legend_rates[2] = "d/dt Gly in component differential_equations (molar)" legend_rates[4] = "d/dt G1P in component differential_equations (molar)" legend_rates[6] = "d/dt G6P in component differential_equations (molar)" legend_rates[7] = "d/dt F6P in component differential_equations (molar)" legend_rates[10] = "d/dt FBP in component differential_equations (molar)" legend_rates[12] = "d/dt DHAP in component differential_equations (molar)" legend_rates[13] = "d/dt G3P in component differential_equations (molar)" legend_rates[11] = "d/dt GAP in component differential_equations (molar)" legend_rates[16] = "d/dt BPG in component differential_equations (molar)" legend_rates[17] = "d/dt P3G in component differential_equations (molar)" legend_rates[18] = "d/dt P2G in component differential_equations (molar)" legend_rates[19] = "d/dt PEP in component differential_equations (molar)" legend_rates[20] = "d/dt PYR in component differential_equations (molar)" legend_rates[21] = "d/dt LAC in component differential_equations (molar)" legend_rates[1] = "d/dt mg in component differential_equations (molar)" legend_rates[0] = "d/dt pH_calc in component differential_equations (dimensionless)" legend_rates[25] = "d/dt LAC_ec in component differential_equations (molar)" legend_rates[29] = "d/dt HCO3_cy in component differential_equations (molar)" legend_rates[28] = "d/dt CO2_cy in component differential_equations (molar)" legend_rates[27] = "d/dt HCO3_ec in component differential_equations (molar)" legend_rates[26] = "d/dt CO2_ec in component differential_equations (molar)" legend_rates[30] = "d/dt MOPS in component differential_equations (molar)" legend_rates[31] = "d/dt MOPSH in component differential_equations (molar)" legend_rates[24] = "d/dt pH_ec in component differential_equations (dimensionless)" return (legend_states, legend_algebraic, legend_voi, legend_constants) def initConsts(): constants = [0.0] * sizeConstants; states = [0.0] * sizeStates; states[0] = 7.14 constants[0] = 1 states[1] = 1e-3 states[2] = 0.02 states[3] = 0.0054 states[4] = 1e-9 states[5] = 1e-8 states[6] = 1e-9 states[7] = 1e-9 states[8] = 0.0055 states[9] = 1e-5 states[10] = 1e-9 states[11] = 1e-9 states[12] = 1e-9 states[13] = 1e-9 states[14] = 0.0000000406654344 states[15] = 0.000406654344 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.0075 states[23] = 14.5e-3 states[24] = 7.439766507222684 states[25] = 0 states[26] = 0.001564425 states[27] = 30e-3 states[28] = 0.001564425 states[29] = 15e-3 states[30] = 0 states[31] = 0 constants[1] = 0.0500000*0.0425000*(760.000-23.8000)*0.00100000 constants[2] = 0.00000 constants[3] = 0.0100000 constants[4] = 0.00170000 constants[5] = 0.00470000 constants[6] = 0.00270000 constants[7] = 0.00460000 constants[8] = 0.00740000 constants[9] = 6.46203e-06 constants[10] = 1.61590 constants[11] = 0.00402000 constants[12] = 0.00270000 constants[13] = 0.000870000*10.0000 constants[14] = 6.00000e-05 constants[15] = 0.0100000 constants[16] = 0.0100000 constants[17] = 13.0000 constants[18] = 0.00390000 constants[19] = 1.50000e-05 constants[20] = 0.00390000 constants[21] = 0.0580000 constants[22] = 0.00434600 constants[23] = 0.00100000 constants[24] = 0.124547 constants[25] = 0.00378440 constants[26] = 4.19392e-05 constants[27] = 2.21796 constants[28] = 0.000711494 constants[29] = 0.00400000 constants[30] = 8.89000 constants[31] = 2.30840 constants[32] = 740.000 constants[33] = 0.0469000*60.0000 constants[34] = 0.00400000 constants[35] = 0.0550000*0.670000*0.00100000 constants[36] = 0.000149000-(0.0550000/1.06000)*0.00100000 constants[37] = 0.00000 constants[38] = 2.70000 constants[39] = 0.00000 constants[40] = 1.63000 constants[41] = 1.57000 constants[42] = 2.45000 constants[43] = 2.26000 constants[44] = 0.00000 constants[45] = 0.980000 constants[46] = 0.00370000 constants[47] = constants[0]*0.00000 constants[48] = 0.0300000 constants[49] = power(10.0000, -6.15700) constants[50] = 1.00000-constants[3] constants[51] = 0.00400000 constants[52] = power(10.0000, -constants[30]) constants[53] = constants[2]/2.00000 constants[54] = constants[2]/2.00000 constants[229] = constants[47] constants[55] = 6.15700+log(constants[48]/constants[1], 10) constants[56] = constants[33]/constants[49] constants[57] = 0.00200000 constants[58] = 0.000150000 constants[59] = 0.0101000 constants[60] = 0.000200000 constants[61] = 0.0150000 constants[62] = 0.00150000 constants[63] = 0.00440000 constants[64] = 0.00520000 constants[65] = 6.30000e-05 constants[66] = 3.00000e-05 constants[67] = 0.383300 constants[68] = 0.000480000 constants[69] = 0.000119000 constants[70] = 0.00370000 constants[71] = 0.000180000 constants[72] = 0.0200000 constants[73] = 8.00000e-05 constants[74] = 0.000250000 constants[75] = 0.00402000 constants[76] = (constants[71]*constants[73])/(constants[72]*constants[74]) constants[77] = 0.00270000 constants[78] = 0.0120000 constants[79] = 5.00000e-05 constants[80] = 0.00200000 constants[81] = 0.00100000 constants[82] = 0.140600 constants[83] = 0.000320000 constants[84] = 0.000610000 constants[85] = 0.0557000 constants[86] = 0.000180000 constants[87] = 1.20000e-05 constants[88] = 0.000220000 constants[89] = 8.00000e-06 constants[90] = 0.0765000 constants[91] = 2.50000e-06 constants[92] = 9.00000e-05 constants[93] = 0.000290000 constants[94] = 8.00000e-07 constants[95] = 3.30000e-06 constants[96] = 0.00870000 constants[97] = 0.00200000 constants[98] = 8.00000e-06 constants[99] = 0.00120000 constants[100] = 0.000350000 constants[101] = 0.0839000 constants[102] = 0.000200000 constants[103] = 1.40000e-05 constants[104] = 0.796600 constants[105] = 0.000100000 constants[106] = 0.000370000 constants[107] = 0.0408000 constants[108] = 8.00000e-05 constants[109] = 0.000300000 constants[110] = 0.00705000 constants[111] = 0.00113000 constants[112] = 0.0956000 constants[113] = 0.000335000 constants[114] = 2.00000e-06 constants[115] = 0.0170000 constants[116] = 0.000849000 constants[117] = 8.82150 constants[118] = 0.00111000 constants[119] = 0.00350000 constants[120] = 0.000135000 constants[121] = 0.00380000 constants[122] = 2.02990 constants[123] = 0.000320000 constants[124] = 0.000270000 constants[125] = 0.000350000 constants[126] = 0.00413880 constants[127] = 298.150 constants[128] = 0.170000 constants[129] = 0.140000 constants[130] = 1.00000 constants[131] = 0.00831400 constants[132] = 1.60000 constants[133] = 0.100000 constants[134] = 1.17582 constants[135] = (constants[134]*((power(constants[133], 1.0/2))/(1.00000+constants[132]*(power(constants[133], 1.0/2)))-(power(constants[128], 1.0/2))/(1.00000+constants[132]*(power(constants[128], 1.0/2)))))/log(10.0000) constants[136] = 1.65000+constants[135]*((power(2.00000, 2.00000)+power(2.00000, 2.00000))-power(0.00000, 2.00000)) constants[137] = 2.91482 constants[138] = 4.19000+constants[135]*((power(4.00000, 2.00000)+power(2.00000, 2.00000))-power(2.00000, 2.00000)) constants[139] = 3.25000+constants[135]*((power(3.00000, 2.00000)+power(2.00000, 2.00000))-power(1.00000, 2.00000)) constants[140] = 1.92000+constants[135]*((power(2.00000, 2.00000)+power(2.00000, 2.00000))-power(0.00000, 2.00000)) constants[141] = 1.60000+constants[135]*((power(2.00000, 2.00000)+power(2.00000, 2.00000))-power(0.00000, 2.00000)) constants[142] = 2.48000+constants[135]*((power(2.00000, 2.00000)+power(2.00000, 2.00000))-power(0.00000, 2.00000)) constants[143] = (constants[137]*(power(constants[128], 1.0/2)))/(1.00000+constants[132]*(power(constants[128], 1.0/2))) constants[144] = 6.75000+constants[135]*((power(2.00000, 2.00000)+power(1.00000, 2.00000))-power(1.00000, 2.00000)) constants[145] = 0.500000 constants[146] = 1.00000 constants[147] = -1096.10 constants[148] = 6.48000+constants[135]*((power(4.00000, 2.00000)+power(1.00000, 2.00000))-power(3.00000, 2.00000)) constants[149] = 1.17000+constants[135]*((power(4.00000, 2.00000)+power(1.00000, 2.00000))-power(3.00000, 2.00000)) constants[150] = 12.0000 constants[151] = -2768.10 constants[152] = 1.00000 constants[153] = 6.38000+constants[135]*((power(3.00000, 2.00000)+power(1.00000, 2.00000))-power(2.00000, 2.00000)) constants[154] = 12.0000 constants[155] = -1906.13 constants[156] = 6.29000+constants[135]*((power(2.00000, 2.00000)+power(1.00000, 2.00000))-power(1.00000, 2.00000)) constants[157] = 12.0000 constants[158] = -1040.45 constants[159] = 0.310000 constants[160] = 4.50000+constants[135]*((power(2.00000, 2.00000)+power(1.00000, 2.00000))-power(1.00000, 2.00000)) constants[161] = 8.00000 constants[162] = 2.30000 constants[163] = 9.00000 constants[164] = 6.09000+constants[135]*((power(2.00000, 2.00000)+power(1.00000, 2.00000))-power(1.00000, 2.00000)) constants[165] = 11.0000 constants[166] = -1756.87 constants[167] = 6.11000+constants[135]*((power(2.00000, 2.00000)+power(1.00000, 2.00000))-power(1.00000, 2.00000)) constants[168] = 11.0000 constants[169] = -1763.94 constants[170] = 5.89000+constants[135]*((power(2.00000, 2.00000)+power(1.00000, 2.00000))-power(1.00000, 2.00000)) constants[171] = 11.0000 constants[172] = -1760.80 constants[173] = 6.40000 constants[174] = 10.0000 constants[175] = -2601.40 constants[176] = 6.45000 constants[177] = 5.00000 constants[178] = -1288.60 constants[179] = 6.22000+constants[135]*((power(2.00000, 2.00000)+power(1.00000, 2.00000))-power(1.00000, 2.00000)) constants[180] = 7.00000 constants[181] = -1339.25 constants[182] = 5.90000 constants[183] = 5.00000 constants[184] = -1296.26 constants[185] = 7.50000 constants[186] = 4.00000 constants[187] = -2356.14 constants[188] = 6.21000 constants[189] = 4.00000 constants[190] = -1502.54 constants[191] = 7.00000 constants[192] = 1.18000 constants[193] = 4.00000 constants[194] = -1496.38 constants[195] = 6.35000 constants[196] = 1.08000 constants[197] = 2.00000 constants[198] = -1263.65 constants[199] = 2.49000 constants[200] = 3.00000 constants[201] = -472.270 constants[202] = 3.67000+constants[135]*((power(1.00000, 2.00000)+power(1.00000, 2.00000))-power(0.00000, 2.00000)) constants[203] = 5.00000 constants[204] = -516.720 constants[205] = -10.0000 constants[206] = 26.0000 constants[207] = 0.00000 constants[208] = 27.0000 constants[209] = 22.6500 constants[210] = 2.00000 constants[211] = -237.190 constants[212] = 1.00000 constants[213] = 0.00000 constants[214] = 1000.00*2.04742 constants[215] = 1000.00*0.0138384 constants[216] = exp(constants[214]/constants[127]+constants[215]) constants[217] = 0.250000 constants[218] = (1.10708-0.00154508*constants[127])+5.95584e-06*(power(constants[127], 2.00000)) constants[219] = constants[218]*((power(constants[217], 1.0/2))/(1.00000+1.60000*(power(constants[217], 1.0/2)))-(power(constants[128], 1.0/2))/(1.00000+1.60000*(power(constants[128], 1.0/2)))) constants[220] = -log(constants[216])+constants[219]*((((power(4.00000, 2.00000)+power(0.00000, 2.00000))-power(2.00000, 2.00000))-power(3.00000, 2.00000))-power(1.00000, 2.00000)) constants[221] = exp(-constants[220]) constants[222] = -constants[131]*constants[127]*log(constants[221]) constants[223] = constants[82] constants[224] = constants[85] constants[225] = constants[96] constants[226] = constants[104] constants[227] = 1.00000 constants[228] = 4.20000 return (states, constants) def computeRates(voi, states, constants): rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic rates[1] = constants[229] algebraic[185] = 1.00000+states[2]/constants[4]+states[3]/constants[51]+(states[2]*states[3])/(constants[4]*constants[5])+states[2]/constants[58]+states[4]/constants[6]+(states[2]*states[4])/(constants[59]*constants[58]) algebraic[0] = states[0] algebraic[186] = 1.40400/(1.00000+power(10.0000, 5.94000-algebraic[0])+power(10.0000, algebraic[0]-7.29000)) algebraic[188] = ((algebraic[186]*constants[46]*states[3])/(constants[57]*constants[51]))/algebraic[185] algebraic[1] = 1.00000+power(10.0000, -algebraic[0]+constants[144])+(states[1]/constants[130])*(power(10.0000, constants[136]))+(constants[129]/constants[130])*(power(10.0000, constants[145])) algebraic[44] = 1.00000+(power(10.0000, -algebraic[0]))*(power(10.0000, constants[164]))+(states[1]/constants[130])*(power(10.0000, constants[142])) algebraic[8] = (constants[147]+constants[146]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[146]) algebraic[50] = (constants[166]+constants[165]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[165]) algebraic[133] = 655.700+constants[205]*log(10.0000)*constants[131]*constants[127]*algebraic[0] algebraic[146] = (algebraic[133]+algebraic[50])-algebraic[8] algebraic[147] = (exp(-algebraic[146]/(constants[131]*constants[127]))*algebraic[44])/algebraic[1] algebraic[187] = (algebraic[186]*constants[46]*constants[58]*constants[59])/(constants[57]*constants[51]*algebraic[147]) algebraic[189] = ((algebraic[187]*states[2])/(constants[58]*constants[59]))/algebraic[185] algebraic[190] = constants[3]*(states[2]*algebraic[188]-states[4]*algebraic[189]) algebraic[191] = (power(states[5]/constants[9], constants[10]))/(1.00000+power(states[5]/constants[9], constants[10])) algebraic[192] = 1.00000+states[2]/constants[61]+states[3]/constants[7]+states[2]/constants[63]+states[4]/constants[8]+(states[2]*states[3])/(constants[61]*constants[60])+(states[2]*states[4])/(constants[62]*constants[63]) algebraic[193] = 1.75000/(1.00000+power(10.0000, 6.12000-algebraic[0])+power(10.0000, algebraic[0]-7.03000)) algebraic[195] = ((algebraic[193]*algebraic[191]*constants[46]*states[3])/(constants[61]*constants[60]))/algebraic[192] algebraic[194] = (algebraic[193]*constants[46]*constants[62]*constants[63])/(constants[61]*constants[60]*algebraic[147]) algebraic[196] = ((algebraic[191]*algebraic[194]*states[2])/(constants[62]*constants[63]))/algebraic[192] algebraic[197] = constants[50]*(states[2]*algebraic[195]-states[4]*algebraic[196]) algebraic[198] = constants[0]*-(algebraic[190]+algebraic[197]) rates[2] = algebraic[198] algebraic[199] = (constants[64]*1.32900)/(1.00000+power(10.0000, -algebraic[0]+6.64000)+power(10.0000, algebraic[0]-8.36000)) algebraic[51] = 1.00000+power(10.0000, -algebraic[0]+constants[167]) algebraic[56] = (constants[169]+constants[168]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[168]) algebraic[149] = algebraic[56]-algebraic[50] algebraic[150] = (exp(-algebraic[149]/(constants[131]*constants[127]))*algebraic[51])/algebraic[44] algebraic[200] = (algebraic[199]*constants[66])/(constants[65]*algebraic[150]) algebraic[201] = ((algebraic[199]*states[4])/constants[65]-(algebraic[200]*states[6])/constants[66])/(1.00000+states[4]/constants[65]+states[6]/constants[66]) algebraic[202] = constants[0]*((algebraic[190]+algebraic[197])-algebraic[201]) rates[4] = algebraic[202] algebraic[203] = constants[67]/(1.00000+power(10.0000, -algebraic[0]+6.94000)+power(10.0000, algebraic[0]-9.35000)) algebraic[57] = 1.00000+power(10.0000, -algebraic[0]+constants[170]) algebraic[62] = (constants[172]+constants[171]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[171]) algebraic[152] = algebraic[62]-algebraic[56] algebraic[153] = (exp(-algebraic[152]/(constants[131]*constants[127]))*algebraic[57])/algebraic[51] algebraic[204] = ((algebraic[203]*constants[68])/constants[69])*algebraic[153] algebraic[205] = ((algebraic[204]*states[6])/constants[68]-(algebraic[203]*states[7])/constants[69])/(1.00000+states[7]/constants[69]+states[6]/constants[68]) algebraic[206] = constants[0]*(algebraic[201]-algebraic[205]) rates[6] = algebraic[206] algebraic[207] = constants[70]/(1.00000+power(algebraic[0]/6.80000, -30.0000)) algebraic[9] = 1.00000+power(10.0000, -algebraic[0]+constants[148])+(states[1]/constants[130])*(power(10.0000, constants[138]))+(constants[129]/constants[130])*(power(10.0000, constants[149])) algebraic[17] = 1.00000+power(10.0000, -algebraic[0]+constants[153])+(states[1]/constants[130])*(power(10.0000, constants[139]))+(constants[129]/constants[130])*(power(10.0000, constants[152])) algebraic[63] = 1.00000+power(10.0000, -algebraic[0]+constants[173])+(states[1]/constants[130])*(power(10.0000, constants[38])) algebraic[16] = (constants[151]+constants[150]*constants[131]*constants[127]*log(10.0000)*algebraic[0])-constants[143]*(16.0000-constants[150]) algebraic[24] = (constants[155]+constants[154]*constants[131]*constants[127]*log(10.0000)*algebraic[0])-constants[143]*(9.00000-constants[154]) algebraic[69] = (constants[175]+constants[174]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(16.0000-constants[174]) algebraic[137] = constants[213]+constants[212]*log(10.0000)*constants[131]*constants[127]*algebraic[0] algebraic[155] = ((algebraic[69]+algebraic[24]+algebraic[137])-algebraic[62])-algebraic[16] algebraic[156] = (exp(-algebraic[155]/(constants[131]*constants[127]))*algebraic[63]*algebraic[17])/(algebraic[57]*algebraic[9]*(power(10.0000, -algebraic[0]))) algebraic[208] = (algebraic[207]*constants[75]*constants[77])/(constants[71]*constants[73]*algebraic[156]) algebraic[209] = constants[17]*(power((((1.00000+states[8]/constants[13])/(1.00000+(constants[15]*states[8])/constants[13]))*(1.00000+(2.71828182846*states[5])/constants[14]))/(1.00000+states[5]/constants[14]), 4.00000)) algebraic[210] = (1.00000+states[7]/constants[71])*(1.00000+states[8]/constants[73])+states[10]/constants[75]+(states[9]/constants[77])*(1.00000+states[10]/constants[75]) algebraic[211] = (1.00000+states[7]/constants[72])*(1.00000+states[8]/constants[74])+states[10]/constants[11]+(states[9]/constants[12])*(1.00000+states[10]/constants[11]) algebraic[212] = ((((algebraic[207]*states[7]*states[8])/(constants[71]*constants[73]))/algebraic[210]-((algebraic[208]*states[9]*states[10])/(constants[77]*constants[75]))/algebraic[210])*(1.00000+constants[76]*algebraic[209]*(power(algebraic[211]/algebraic[210], 3.00000))))/(1.00000+algebraic[209]*(power(algebraic[211]/algebraic[210], 4.00000))) algebraic[213] = constants[0]*(algebraic[205]-algebraic[212]) rates[7] = algebraic[213] algebraic[214] = (constants[78]*1.01300)/(1.00000+power(10.0000, -algebraic[0]+5.32000)+power(10.0000, algebraic[0]-9.15000)) algebraic[70] = 1.00000+power(10.0000, -algebraic[0]+constants[176]) algebraic[83] = 1.00000+power(10.0000, -algebraic[0]+constants[182])+(states[1]/constants[130])*(power(10.0000, constants[41])) algebraic[75] = (constants[178]+constants[177]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[177]) algebraic[89] = (constants[184]+constants[183]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[183]) algebraic[158] = (algebraic[89]+algebraic[75])-algebraic[69] algebraic[159] = (exp(-algebraic[158]/(constants[131]*constants[127]))*algebraic[70]*algebraic[63])/algebraic[83] algebraic[215] = (algebraic[214]*constants[81]*constants[80])/(constants[79]*algebraic[159]*1.00000) algebraic[216] = ((algebraic[214]*states[10])/constants[79]-(algebraic[215]*states[11]*states[12])/(constants[81]*constants[80]))/(1.00000+states[10]/constants[79]+states[11]/constants[81]+states[12]/constants[80]) algebraic[217] = constants[0]*(algebraic[212]-algebraic[216]) rates[10] = algebraic[217] algebraic[161] = algebraic[89]-algebraic[75] algebraic[162] = (exp(-algebraic[161]/(constants[131]*constants[127]))*algebraic[83])/algebraic[70] algebraic[218] = (constants[223]*constants[84])/(constants[83]*algebraic[162]) algebraic[219] = ((constants[223]*states[11])/constants[83]-(algebraic[218]*states[12])/constants[84])/(1.00000+states[11]/constants[83]+states[12]/constants[84]) algebraic[220] = (1.00000+states[13]/constants[86]+states[14]/constants[89])*(1.00000+states[12]/constants[88]+states[15]/constants[89]) algebraic[76] = 1.00000+power(10.0000, -algebraic[0]+constants[179])+(states[1]/constants[130])*(power(10.0000, constants[40])) algebraic[82] = (constants[181]+constants[180]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[180]) algebraic[134] = (constants[207]+constants[206]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(1.00000-constants[206]) algebraic[135] = (constants[209]+constants[208]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[208]) algebraic[167] = ((algebraic[137]+algebraic[135]+algebraic[89])-algebraic[134])-algebraic[82] algebraic[168] = (exp(-algebraic[167]/(constants[131]*constants[127]))*algebraic[83])/(algebraic[76]*(power(10.0000, -algebraic[0]))) algebraic[221] = (constants[224]*constants[86]*constants[87]*algebraic[168])/(constants[88]*constants[89]) algebraic[222] = ((algebraic[221]*states[13]*states[15])/(constants[86]*constants[87])-(constants[224]*states[12]*states[14])/(constants[88]*constants[89]))/algebraic[220] algebraic[223] = constants[0]*(algebraic[216]+algebraic[219]+algebraic[222]) rates[12] = algebraic[223] algebraic[224] = constants[0]*-algebraic[222] rates[13] = algebraic[224] algebraic[225] = 1.00000+states[3]/constants[93]+states[11]/constants[91]+states[15]/constants[92]+(states[11]*states[15])/(constants[91]*constants[92])+(states[11]*states[15]*states[3])/(constants[91]*constants[92]*constants[93])+states[16]/constants[94]+states[14]/constants[95]+(states[16]*states[14])/(constants[95]*constants[94]) algebraic[226] = constants[90]*0.000700000*exp(algebraic[0]*0.897900) algebraic[90] = 1.00000+power(10.0000, -algebraic[0]+constants[185]) algebraic[95] = (constants[187]+constants[186]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(16.0000-constants[186]) algebraic[164] = (((algebraic[95]+algebraic[135]+algebraic[137])-algebraic[8])-algebraic[75])-algebraic[134] algebraic[165] = (exp(-algebraic[164]/(constants[131]*constants[127]))*algebraic[90])/(algebraic[1]*algebraic[70]*(power(10.0000, -algebraic[0]))) algebraic[227] = (algebraic[226]*constants[94]*constants[95]*1.00000)/(constants[91]*constants[93]*constants[92]*algebraic[165]) algebraic[228] = ((algebraic[226]*states[11]*states[15]*states[3])/(constants[92]*constants[91]*constants[93])-(algebraic[227]*states[16]*states[14])/(constants[94]*constants[95]))/algebraic[225] algebraic[229] = constants[0]*((algebraic[216]-algebraic[219])-algebraic[228]) rates[11] = algebraic[229] algebraic[96] = 1.00000+power(10.0000, -algebraic[0]+6.21000) algebraic[101] = (constants[190]+constants[189]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(9.00000-constants[189]) algebraic[170] = ((algebraic[16]+algebraic[101])-algebraic[95])-algebraic[24] algebraic[171] = (exp(-algebraic[170]/(constants[131]*constants[127]))*algebraic[9]*algebraic[96])/(algebraic[90]*algebraic[17]) algebraic[230] = ((constants[225]*constants[97]*constants[98])/(constants[99]*constants[100]))*algebraic[171] algebraic[231] = 1.00000+states[9]/constants[98]+states[16]/constants[97]+(states[16]*states[9])/(constants[97]*constants[98])+states[17]/constants[99]+states[8]/constants[100]+(states[17]*states[8])/(constants[99]*constants[100]) algebraic[232] = ((algebraic[230]*states[16]*states[9])/(constants[98]*constants[97])-(constants[225]*states[8]*states[17])/(constants[100]*constants[99]))/algebraic[231] algebraic[233] = constants[0]*(algebraic[228]-algebraic[232]) rates[16] = algebraic[233] algebraic[234] = (constants[101]*0.989000)/(1.00000+power(10.0000, -algebraic[0]+5.62000)+power(10.0000, algebraic[0]-8.74000)) algebraic[102] = 1.00000+power(10.0000, -algebraic[0]+constants[191])+(states[1]/constants[130])*(power(10.0000, constants[42]))+(constants[129]/constants[130])*(power(10.0000, constants[192])) algebraic[109] = (constants[194]+constants[193]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(9.00000-constants[193]) algebraic[173] = algebraic[109]-algebraic[101] algebraic[174] = (exp(-algebraic[173]/(constants[131]*constants[127]))*algebraic[102])/algebraic[96] algebraic[235] = (algebraic[234]*constants[103])/(constants[102]*algebraic[174]) algebraic[236] = ((algebraic[234]*states[17])/constants[102]-(algebraic[235]*states[18])/constants[103])/(1.00000+states[17]/constants[102]+states[18]/constants[103]) algebraic[237] = constants[0]*(algebraic[232]-algebraic[236]) rates[17] = algebraic[237] algebraic[110] = 1.00000+power(10.0000, constants[195]-algebraic[0])+(states[1]/constants[130])*(power(10.0000, constants[43]))+(constants[129]/constants[130])*(power(10.0000, constants[196])) algebraic[117] = (constants[198]+constants[197]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(9.00000-constants[197]) algebraic[136] = (constants[211]+constants[210]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(0.00000-constants[210]) algebraic[176] = (algebraic[136]+algebraic[117])-algebraic[109] algebraic[177] = (exp(-algebraic[176]/(constants[131]*constants[127]))*algebraic[110])/algebraic[102] algebraic[238] = (constants[226]*constants[106])/(constants[105]*algebraic[177]) algebraic[239] = ((constants[226]*states[18])/constants[105]-(algebraic[238]*states[19])/constants[106])/(1.00000+states[19]/constants[106]+states[18]/constants[105]) algebraic[240] = constants[0]*(algebraic[236]-algebraic[239]) rates[18] = algebraic[240] algebraic[241] = (constants[107]*1.05000)/(1.00000+power(10.0000, -algebraic[0]+5.58000)+power(10.0000, algebraic[0]-8.79000)) algebraic[118] = 1.00000+power(10.0000, -algebraic[0]+constants[199]) algebraic[123] = (constants[201]+constants[200]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(1.00000-constants[200]) algebraic[179] = (((algebraic[123]+algebraic[16])-algebraic[137])-algebraic[117])-algebraic[24] algebraic[180] = (exp(-algebraic[179]/(constants[131]*constants[127]))*algebraic[118]*algebraic[9]*(power(10.0000, -algebraic[0])))/(algebraic[110]*algebraic[17]) algebraic[242] = (algebraic[241]*constants[110]*constants[111])/(constants[108]*constants[109]*algebraic[180]) algebraic[243] = ((algebraic[241]*states[19]*states[9])/(constants[108]*constants[109])-(algebraic[242]*states[20]*states[8])/(constants[110]*constants[111]))/(1.00000+states[19]/constants[108]+states[9]/constants[109]+(states[19]*states[9])/(constants[108]*constants[109])+states[8]/constants[111]+states[20]/constants[110]+(states[20]*states[8])/(constants[110]*constants[111])) algebraic[244] = constants[0]*(algebraic[239]-algebraic[243]) rates[19] = algebraic[244] algebraic[245] = constants[112]*(-0.113400*algebraic[0]+1.60690) algebraic[124] = 1.00000+power(10.0000, -algebraic[0]+constants[202])+(states[1]/constants[130])*(power(10.0000, constants[45])) algebraic[132] = (constants[204]+constants[203]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(1.00000-constants[203]) algebraic[182] = (((algebraic[132]+algebraic[134])-algebraic[123])-algebraic[135])-algebraic[137] algebraic[183] = (exp(-algebraic[182]/(constants[131]*constants[127]))*algebraic[124]*(power(10.0000, -algebraic[0])))/algebraic[118] algebraic[246] = (algebraic[245]*constants[115]*constants[116])/(constants[113]*constants[114]*algebraic[183]) algebraic[247] = ((algebraic[245]*states[20]*states[14])/(constants[113]*constants[114])-(algebraic[246]*states[21]*states[15])/(constants[115]*constants[116]))/(1.00000+states[20]/constants[113]+states[14]/constants[114]+(states[20]*states[14])/(constants[113]*constants[114])+states[21]/constants[115]+states[15]/constants[116]+(states[21]*states[15])/(constants[115]*constants[116])) algebraic[248] = constants[0]*((-algebraic[228]-algebraic[222])+algebraic[247]) rates[15] = algebraic[248] algebraic[249] = constants[0]*((algebraic[228]+algebraic[222])-algebraic[247]) rates[14] = algebraic[249] algebraic[250] = constants[0]*(algebraic[243]-algebraic[247]) rates[20] = algebraic[250] algebraic[32] = 1.00000+power(10.0000, -algebraic[0]+constants[160])+(states[1]/constants[130])*(power(10.0000, constants[141]))+(constants[129]/constants[130])*(power(10.0000, constants[159])) algebraic[39] = 1.00000+power(10.0000, -algebraic[0]+constants[162]) algebraic[140] = (exp(-constants[222]/(constants[131]*constants[127]))*(power(10.0000, -algebraic[0]))*algebraic[9]*algebraic[39])/(algebraic[32]*algebraic[17]) algebraic[251] = ((constants[117]/algebraic[140])*constants[119]*constants[121])/(constants[120]*constants[118]) algebraic[252] = ((algebraic[251]*states[8]*states[22])/(constants[119]*constants[121])-(constants[117]*states[9]*states[23])/(constants[120]*constants[118]))/(1.00000+states[9]/constants[120]+states[23]/constants[18]+(states[23]*states[9])/(constants[120]*constants[118])+states[8]/constants[119]+(states[22]*states[8])/(constants[121]*constants[119])) algebraic[253] = constants[0]*algebraic[252] rates[23] = algebraic[253] algebraic[254] = constants[0]*-algebraic[252] rates[22] = algebraic[254] algebraic[25] = 1.00000+power(10.0000, constants[156]-algebraic[0])+(states[1]/constants[130])*(power(10.0000, constants[140])) algebraic[31] = (constants[158]+constants[157]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[157]) algebraic[143] = (2.00000*algebraic[24]-algebraic[16])-algebraic[31] algebraic[144] = (exp(-algebraic[143]/(constants[131]*constants[127]))*(power(algebraic[17], 2.00000)))/(algebraic[9]*algebraic[25]) algebraic[255] = (constants[122]*(power(constants[125], 2.00000)))/(constants[123]*constants[124]*algebraic[144]) algebraic[256] = ((constants[122]*states[8]*states[5])/(constants[124]*constants[123])-algebraic[255]*(power(states[9]/constants[125], 2.00000)))/(1.00000+states[8]/constants[124]+states[5]/constants[123]+(states[8]*states[5])/(constants[124]*constants[123])+(2.00000*states[9])/constants[125]+(power(states[9], 2.00000))/(power(constants[125], 2.00000))) algebraic[257] = constants[0]*-algebraic[256] rates[5] = algebraic[257] algebraic[258] = constants[24]*states[8] algebraic[259] = (constants[25]*(power(states[9]/constants[26], constants[27])))/(1.00000+power(states[9]/constants[26], constants[27])) algebraic[261] = constants[0]*(((((-algebraic[252]-algebraic[256])-algebraic[212])+algebraic[232]+algebraic[243])-algebraic[258])+algebraic[259]) rates[8] = algebraic[261] algebraic[262] = constants[0]*(((((algebraic[252]+2.00000*algebraic[256]+algebraic[212])-algebraic[232])-algebraic[243])+algebraic[258])-algebraic[259]) rates[9] = algebraic[262] algebraic[263] = constants[0]*(((-(algebraic[190]+algebraic[197])-algebraic[228])+algebraic[258])-algebraic[259]) rates[3] = algebraic[263] algebraic[260] = (constants[28]*((states[21]*(power(10.0000, -algebraic[0])))/(constants[52]*constants[29])-(states[25]*(power(10.0000, -states[24])))/(constants[52]*constants[29])))/((states[21]*(power(10.0000, -algebraic[0])))/(constants[52]*constants[29])+2.00000+(power(10.0000, -algebraic[0]))/constants[52]+states[21]/constants[29]+states[25]/constants[29]+(power(10.0000, -states[24]))/constants[52]+(states[25]*(power(10.0000, -states[24])))/(constants[29]*constants[52])) algebraic[265] = constants[0]*(algebraic[247]-algebraic[260]) rates[21] = algebraic[265] algebraic[266] = constants[0]*((constants[35]*algebraic[260])/constants[36]-(constants[34]*states[25])/constants[36]) rates[25] = algebraic[266] algebraic[269] = constants[32]*(constants[33]*states[28]-constants[56]*(power(10.0000, -algebraic[0]))*states[29]) algebraic[272] = constants[0]*algebraic[269] rates[29] = algebraic[272] algebraic[264] = algebraic[259]/constants[228] algebraic[267] = constants[227]*algebraic[264] algebraic[268] = constants[31]*(states[28]-states[26]) algebraic[273] = constants[0]*((algebraic[267]-algebraic[269])-algebraic[268]) rates[28] = algebraic[273] algebraic[270] = 1.00000*(constants[33]*states[26]-constants[56]*(power(10.0000, -states[24]))*states[27]) algebraic[274] = constants[0]*((constants[34]*(constants[48]-states[27]))/constants[36]+algebraic[270]) rates[27] = algebraic[274] algebraic[275] = constants[0]*(((constants[34]*(constants[1]-states[26]))/constants[36]+(algebraic[268]*constants[35])/constants[36])-algebraic[270]) rates[26] = algebraic[275] algebraic[276] = 1.00000e+07*((power(10.0000, -7.18400))*states[30]-states[31]*(power(10.0000, -states[24]))) algebraic[278] = constants[0]*((constants[34]*(constants[53]-states[30]))/constants[36]-algebraic[276]) rates[30] = algebraic[278] algebraic[279] = constants[0]*((constants[34]*(constants[54]-states[31]))/constants[36]+algebraic[276]) rates[31] = algebraic[279] algebraic[280] = (constants[0]*((constants[34]*(power(10.0000, -constants[55])-power(10.0000, -states[24])))/constants[36]+algebraic[270]/1.00000+algebraic[276]/1.00000+(constants[35]*algebraic[260])/(constants[36]*1.00000)))/(-log(10.0000)*(power(10.0000, -states[24]))) rates[24] = algebraic[280] algebraic[2] = 1.00000/algebraic[1] algebraic[3] = (power(10.0000, -algebraic[0]+constants[144]))*algebraic[2] algebraic[7] = ((power(10.0000, constants[144]))*algebraic[1]*(1.00000-algebraic[3]))/(constants[130]*(power(algebraic[1], 2.00000))) algebraic[10] = 1.00000/algebraic[9] algebraic[13] = (power(10.0000, -algebraic[0]+constants[148]))*algebraic[10] algebraic[15] = ((power(10.0000, constants[148]))*algebraic[9]*(1.00000-algebraic[13]))/(constants[130]*(power(algebraic[9], 2.00000))) algebraic[18] = 1.00000/algebraic[17] algebraic[21] = (power(10.0000, -algebraic[0]+constants[153]))*algebraic[18] algebraic[23] = ((power(10.0000, constants[153]))*algebraic[17]*(1.00000-algebraic[21]))/(constants[130]*(power(algebraic[17], 2.00000))) algebraic[30] = ((power(10.0000, constants[156]))*(algebraic[25]-power(10.0000, -algebraic[0]+constants[156])))/(constants[130]*(power(algebraic[25], 2.00000))) algebraic[38] = ((power(10.0000, constants[160]))*(algebraic[32]-power(10.0000, -algebraic[0]+constants[160])))/(constants[130]*(power(algebraic[32], 2.00000))) algebraic[43] = ((power(10.0000, constants[162]))*(algebraic[39]-power(10.0000, -algebraic[0]+constants[162])))/(constants[130]*(power(algebraic[39], 2.00000))) algebraic[49] = ((power(10.0000, constants[164]))*(algebraic[44]-power(10.0000, -algebraic[0]+constants[164])))/(constants[130]*(power(algebraic[44], 2.00000))) algebraic[55] = ((power(10.0000, constants[167]))*(algebraic[51]-power(10.0000, -algebraic[0]+constants[167])))/(constants[130]*(power(algebraic[51], 2.00000))) algebraic[61] = ((power(10.0000, constants[170]))*(algebraic[57]-power(10.0000, -algebraic[0]+constants[170])))/(constants[130]*(power(algebraic[57], 2.00000))) algebraic[64] = 1.00000/algebraic[63] algebraic[66] = algebraic[64]*(power(10.0000, -algebraic[0]+constants[173])) algebraic[68] = ((power(10.0000, constants[173]))*algebraic[63]*(1.00000-algebraic[66]))/(constants[130]*(power(algebraic[63], 2.00000))) algebraic[74] = ((power(10.0000, constants[176]))*(algebraic[70]-power(10.0000, -algebraic[0]+constants[176])))/(constants[130]*(power(algebraic[70], 2.00000))) algebraic[81] = ((power(10.0000, constants[179]))*(algebraic[76]-power(10.0000, -algebraic[0]+constants[179])))/(constants[130]*(power(algebraic[76], 2.00000))) algebraic[88] = ((power(10.0000, constants[182]))*(algebraic[83]-power(10.0000, -algebraic[0]+constants[182])))/(constants[130]*(power(algebraic[83], 2.00000))) algebraic[94] = ((power(10.0000, constants[185]))*(algebraic[90]-power(10.0000, -algebraic[0]+constants[185])))/(constants[130]*(power(algebraic[90], 2.00000))) algebraic[100] = ((power(10.0000, constants[188]))*(algebraic[96]-power(10.0000, -algebraic[0]+constants[188])))/(constants[130]*(power(algebraic[96], 2.00000))) algebraic[108] = ((power(10.0000, constants[191]))*(algebraic[102]-power(10.0000, -algebraic[0]+constants[191])))/(constants[130]*(power(algebraic[102], 2.00000))) algebraic[116] = ((power(10.0000, constants[195]))*(algebraic[110]-power(10.0000, -algebraic[0]+constants[195])))/(constants[130]*(power(algebraic[110], 2.00000))) algebraic[122] = ((power(10.0000, constants[199]))*(algebraic[118]-power(10.0000, -algebraic[0]+constants[199])))/(constants[130]*(power(algebraic[118], 2.00000))) algebraic[131] = ((power(10.0000, constants[202]))*(algebraic[124]-power(10.0000, -algebraic[0]+constants[202])))/(constants[130]*(power(algebraic[124], 2.00000))) algebraic[271] = log(10.0000)*(power(10.0000, -algebraic[0]))*constants[130]*(1.00000+algebraic[7]*states[3]+algebraic[15]*states[8]+algebraic[23]*states[9]+algebraic[30]*states[5]+algebraic[38]*states[23]+algebraic[43]*states[22]+algebraic[49]*states[4]+algebraic[55]*states[6]+algebraic[61]*states[7]+algebraic[68]*states[10]+algebraic[74]*states[11]+algebraic[88]*states[12]+algebraic[81]*states[13]+algebraic[94]*states[16]+algebraic[100]*states[17]+algebraic[108]*states[18]+algebraic[116]*states[19]+algebraic[122]*states[20]+algebraic[131]*states[21]) algebraic[14] = algebraic[13] algebraic[22] = algebraic[21] algebraic[26] = 1.00000/algebraic[25] algebraic[28] = algebraic[26]*(power(10.0000, constants[156]-algebraic[0])) algebraic[29] = algebraic[28] algebraic[142] = ((2.00000*algebraic[22]-algebraic[14])-algebraic[29])+((2.00000*constants[154]-constants[150])-constants[157]) algebraic[6] = algebraic[3] algebraic[184] = ((algebraic[22]+algebraic[6])-algebraic[14])+(((constants[154]+constants[146])-constants[150])-constants[210]) algebraic[33] = 1.00000/algebraic[32] algebraic[36] = (power(10.0000, -algebraic[0]+constants[160]))*algebraic[33] algebraic[37] = algebraic[36] algebraic[40] = 1.00000/algebraic[39] algebraic[41] = algebraic[40]*(power(10.0000, -algebraic[0]+constants[162])) algebraic[42] = 0.00000*algebraic[40]+1.00000*algebraic[41] algebraic[138] = (((algebraic[14]+algebraic[42])-algebraic[37])-algebraic[22])+(((constants[150]+constants[163])-constants[161])-constants[154]) algebraic[277] = algebraic[138]*-algebraic[252] algebraic[45] = 1.00000/algebraic[44] algebraic[46] = algebraic[45]*(power(10.0000, -algebraic[0]+constants[164])) algebraic[48] = algebraic[46] algebraic[145] = (algebraic[48]-algebraic[6])+(1.00000-constants[146]) algebraic[52] = 1.00000/algebraic[51] algebraic[53] = algebraic[52]*(power(10.0000, -algebraic[0]+constants[167])) algebraic[54] = algebraic[53] algebraic[148] = (algebraic[54]-algebraic[48])+(constants[168]-constants[165]) algebraic[58] = 1.00000/algebraic[57] algebraic[59] = algebraic[58]*(power(10.0000, -algebraic[0]+constants[170])) algebraic[60] = algebraic[59] algebraic[151] = (algebraic[60]-algebraic[54])+(constants[171]-constants[168]) algebraic[67] = algebraic[66] algebraic[154] = (((((algebraic[22]+algebraic[67])-algebraic[60])-algebraic[14])+constants[154]+constants[174])-constants[171])-constants[150] algebraic[71] = 1.00000/algebraic[70] algebraic[72] = algebraic[71]*(power(10.0000, -algebraic[0]+constants[176])) algebraic[73] = algebraic[72] algebraic[84] = 1.00000/algebraic[83] algebraic[86] = algebraic[84]*(power(10.0000, -algebraic[0]+constants[182])) algebraic[87] = algebraic[86] algebraic[157] = ((algebraic[87]+algebraic[73])-algebraic[67])+((constants[183]+constants[177])-constants[174]) algebraic[160] = (algebraic[87]-algebraic[73])+(constants[183]-constants[177]) algebraic[91] = 1.00000/algebraic[90] algebraic[92] = algebraic[91]*(power(10.0000, -algebraic[0]+constants[185])) algebraic[93] = algebraic[92] algebraic[163] = ((algebraic[93]-algebraic[73])-algebraic[6])+((((constants[186]+constants[208])-constants[177])-constants[146])-constants[206]) algebraic[77] = 1.00000/algebraic[76] algebraic[79] = algebraic[77]*(power(10.0000, -algebraic[0]+constants[179])) algebraic[80] = algebraic[79] algebraic[166] = (algebraic[87]-algebraic[80])+(((constants[183]+constants[208])-constants[206])-constants[180]) algebraic[97] = 1.00000/algebraic[96] algebraic[98] = algebraic[97]*(power(10.0000, -algebraic[0]+6.21000)) algebraic[99] = algebraic[98] algebraic[169] = (((algebraic[99]+algebraic[14])-algebraic[93])-algebraic[22])+(((constants[189]+constants[150])-constants[186])-constants[154]) algebraic[103] = 1.00000/algebraic[102] algebraic[104] = algebraic[103]*(power(10.0000, -algebraic[0]+constants[191])) algebraic[107] = algebraic[104] algebraic[172] = (algebraic[107]-algebraic[99])+(constants[193]-constants[189]) algebraic[111] = 1.00000/algebraic[110] algebraic[112] = algebraic[111]*(power(10.0000, constants[195]-algebraic[0])) algebraic[115] = algebraic[112] algebraic[175] = (algebraic[115]-algebraic[107])+((constants[210]+constants[197])-constants[193]) algebraic[119] = 1.00000/algebraic[118] algebraic[120] = algebraic[119]*(power(10.0000, -algebraic[0]+constants[199])) algebraic[121] = algebraic[120] algebraic[178] = (((algebraic[121]+algebraic[14])-algebraic[115])-algebraic[22])+(((constants[200]+constants[150])-constants[197])-constants[154]) algebraic[127] = 1.00000/algebraic[124] algebraic[129] = algebraic[127]*(power(10.0000, -algebraic[0]+constants[202])) algebraic[130] = algebraic[129] algebraic[181] = (algebraic[130]-algebraic[121])+(((constants[203]+constants[206])-constants[200])-constants[208]) algebraic[281] = algebraic[148]*algebraic[201]+algebraic[145]*(algebraic[190]+algebraic[197])+algebraic[151]*algebraic[205]+algebraic[154]*algebraic[212]+algebraic[157]*algebraic[216]+algebraic[160]*algebraic[219]+algebraic[163]*algebraic[228]+algebraic[169]*algebraic[232]+algebraic[172]*algebraic[236]+algebraic[175]*algebraic[239]+algebraic[178]*algebraic[243]+algebraic[181]*algebraic[247]+algebraic[166]*algebraic[222] algebraic[282] = ((algebraic[277]+algebraic[281]+algebraic[184]*algebraic[258]+-algebraic[184]*algebraic[259]+algebraic[142]*algebraic[256])-algebraic[269])+algebraic[260] algebraic[283] = algebraic[282]/(constants[126]+algebraic[271]) algebraic[284] = constants[0]*algebraic[283] rates[0] = algebraic[284] return(rates) def computeAlgebraic(constants, states, voi): algebraic = array([[0.0] * len(voi)] * sizeAlgebraic) states = array(states) voi = array(voi) algebraic[185] = 1.00000+states[2]/constants[4]+states[3]/constants[51]+(states[2]*states[3])/(constants[4]*constants[5])+states[2]/constants[58]+states[4]/constants[6]+(states[2]*states[4])/(constants[59]*constants[58]) algebraic[0] = states[0] algebraic[186] = 1.40400/(1.00000+power(10.0000, 5.94000-algebraic[0])+power(10.0000, algebraic[0]-7.29000)) algebraic[188] = ((algebraic[186]*constants[46]*states[3])/(constants[57]*constants[51]))/algebraic[185] algebraic[1] = 1.00000+power(10.0000, -algebraic[0]+constants[144])+(states[1]/constants[130])*(power(10.0000, constants[136]))+(constants[129]/constants[130])*(power(10.0000, constants[145])) algebraic[44] = 1.00000+(power(10.0000, -algebraic[0]))*(power(10.0000, constants[164]))+(states[1]/constants[130])*(power(10.0000, constants[142])) algebraic[8] = (constants[147]+constants[146]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[146]) algebraic[50] = (constants[166]+constants[165]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[165]) algebraic[133] = 655.700+constants[205]*log(10.0000)*constants[131]*constants[127]*algebraic[0] algebraic[146] = (algebraic[133]+algebraic[50])-algebraic[8] algebraic[147] = (exp(-algebraic[146]/(constants[131]*constants[127]))*algebraic[44])/algebraic[1] algebraic[187] = (algebraic[186]*constants[46]*constants[58]*constants[59])/(constants[57]*constants[51]*algebraic[147]) algebraic[189] = ((algebraic[187]*states[2])/(constants[58]*constants[59]))/algebraic[185] algebraic[190] = constants[3]*(states[2]*algebraic[188]-states[4]*algebraic[189]) algebraic[191] = (power(states[5]/constants[9], constants[10]))/(1.00000+power(states[5]/constants[9], constants[10])) algebraic[192] = 1.00000+states[2]/constants[61]+states[3]/constants[7]+states[2]/constants[63]+states[4]/constants[8]+(states[2]*states[3])/(constants[61]*constants[60])+(states[2]*states[4])/(constants[62]*constants[63]) algebraic[193] = 1.75000/(1.00000+power(10.0000, 6.12000-algebraic[0])+power(10.0000, algebraic[0]-7.03000)) algebraic[195] = ((algebraic[193]*algebraic[191]*constants[46]*states[3])/(constants[61]*constants[60]))/algebraic[192] algebraic[194] = (algebraic[193]*constants[46]*constants[62]*constants[63])/(constants[61]*constants[60]*algebraic[147]) algebraic[196] = ((algebraic[191]*algebraic[194]*states[2])/(constants[62]*constants[63]))/algebraic[192] algebraic[197] = constants[50]*(states[2]*algebraic[195]-states[4]*algebraic[196]) algebraic[198] = constants[0]*-(algebraic[190]+algebraic[197]) algebraic[199] = (constants[64]*1.32900)/(1.00000+power(10.0000, -algebraic[0]+6.64000)+power(10.0000, algebraic[0]-8.36000)) algebraic[51] = 1.00000+power(10.0000, -algebraic[0]+constants[167]) algebraic[56] = (constants[169]+constants[168]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[168]) algebraic[149] = algebraic[56]-algebraic[50] algebraic[150] = (exp(-algebraic[149]/(constants[131]*constants[127]))*algebraic[51])/algebraic[44] algebraic[200] = (algebraic[199]*constants[66])/(constants[65]*algebraic[150]) algebraic[201] = ((algebraic[199]*states[4])/constants[65]-(algebraic[200]*states[6])/constants[66])/(1.00000+states[4]/constants[65]+states[6]/constants[66]) algebraic[202] = constants[0]*((algebraic[190]+algebraic[197])-algebraic[201]) algebraic[203] = constants[67]/(1.00000+power(10.0000, -algebraic[0]+6.94000)+power(10.0000, algebraic[0]-9.35000)) algebraic[57] = 1.00000+power(10.0000, -algebraic[0]+constants[170]) algebraic[62] = (constants[172]+constants[171]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[171]) algebraic[152] = algebraic[62]-algebraic[56] algebraic[153] = (exp(-algebraic[152]/(constants[131]*constants[127]))*algebraic[57])/algebraic[51] algebraic[204] = ((algebraic[203]*constants[68])/constants[69])*algebraic[153] algebraic[205] = ((algebraic[204]*states[6])/constants[68]-(algebraic[203]*states[7])/constants[69])/(1.00000+states[7]/constants[69]+states[6]/constants[68]) algebraic[206] = constants[0]*(algebraic[201]-algebraic[205]) algebraic[207] = constants[70]/(1.00000+power(algebraic[0]/6.80000, -30.0000)) algebraic[9] = 1.00000+power(10.0000, -algebraic[0]+constants[148])+(states[1]/constants[130])*(power(10.0000, constants[138]))+(constants[129]/constants[130])*(power(10.0000, constants[149])) algebraic[17] = 1.00000+power(10.0000, -algebraic[0]+constants[153])+(states[1]/constants[130])*(power(10.0000, constants[139]))+(constants[129]/constants[130])*(power(10.0000, constants[152])) algebraic[63] = 1.00000+power(10.0000, -algebraic[0]+constants[173])+(states[1]/constants[130])*(power(10.0000, constants[38])) algebraic[16] = (constants[151]+constants[150]*constants[131]*constants[127]*log(10.0000)*algebraic[0])-constants[143]*(16.0000-constants[150]) algebraic[24] = (constants[155]+constants[154]*constants[131]*constants[127]*log(10.0000)*algebraic[0])-constants[143]*(9.00000-constants[154]) algebraic[69] = (constants[175]+constants[174]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(16.0000-constants[174]) algebraic[137] = constants[213]+constants[212]*log(10.0000)*constants[131]*constants[127]*algebraic[0] algebraic[155] = ((algebraic[69]+algebraic[24]+algebraic[137])-algebraic[62])-algebraic[16] algebraic[156] = (exp(-algebraic[155]/(constants[131]*constants[127]))*algebraic[63]*algebraic[17])/(algebraic[57]*algebraic[9]*(power(10.0000, -algebraic[0]))) algebraic[208] = (algebraic[207]*constants[75]*constants[77])/(constants[71]*constants[73]*algebraic[156]) algebraic[209] = constants[17]*(power((((1.00000+states[8]/constants[13])/(1.00000+(constants[15]*states[8])/constants[13]))*(1.00000+(2.71828182846*states[5])/constants[14]))/(1.00000+states[5]/constants[14]), 4.00000)) algebraic[210] = (1.00000+states[7]/constants[71])*(1.00000+states[8]/constants[73])+states[10]/constants[75]+(states[9]/constants[77])*(1.00000+states[10]/constants[75]) algebraic[211] = (1.00000+states[7]/constants[72])*(1.00000+states[8]/constants[74])+states[10]/constants[11]+(states[9]/constants[12])*(1.00000+states[10]/constants[11]) algebraic[212] = ((((algebraic[207]*states[7]*states[8])/(constants[71]*constants[73]))/algebraic[210]-((algebraic[208]*states[9]*states[10])/(constants[77]*constants[75]))/algebraic[210])*(1.00000+constants[76]*algebraic[209]*(power(algebraic[211]/algebraic[210], 3.00000))))/(1.00000+algebraic[209]*(power(algebraic[211]/algebraic[210], 4.00000))) algebraic[213] = constants[0]*(algebraic[205]-algebraic[212]) algebraic[214] = (constants[78]*1.01300)/(1.00000+power(10.0000, -algebraic[0]+5.32000)+power(10.0000, algebraic[0]-9.15000)) algebraic[70] = 1.00000+power(10.0000, -algebraic[0]+constants[176]) algebraic[83] = 1.00000+power(10.0000, -algebraic[0]+constants[182])+(states[1]/constants[130])*(power(10.0000, constants[41])) algebraic[75] = (constants[178]+constants[177]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[177]) algebraic[89] = (constants[184]+constants[183]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[183]) algebraic[158] = (algebraic[89]+algebraic[75])-algebraic[69] algebraic[159] = (exp(-algebraic[158]/(constants[131]*constants[127]))*algebraic[70]*algebraic[63])/algebraic[83] algebraic[215] = (algebraic[214]*constants[81]*constants[80])/(constants[79]*algebraic[159]*1.00000) algebraic[216] = ((algebraic[214]*states[10])/constants[79]-(algebraic[215]*states[11]*states[12])/(constants[81]*constants[80]))/(1.00000+states[10]/constants[79]+states[11]/constants[81]+states[12]/constants[80]) algebraic[217] = constants[0]*(algebraic[212]-algebraic[216]) algebraic[161] = algebraic[89]-algebraic[75] algebraic[162] = (exp(-algebraic[161]/(constants[131]*constants[127]))*algebraic[83])/algebraic[70] algebraic[218] = (constants[223]*constants[84])/(constants[83]*algebraic[162]) algebraic[219] = ((constants[223]*states[11])/constants[83]-(algebraic[218]*states[12])/constants[84])/(1.00000+states[11]/constants[83]+states[12]/constants[84]) algebraic[220] = (1.00000+states[13]/constants[86]+states[14]/constants[89])*(1.00000+states[12]/constants[88]+states[15]/constants[89]) algebraic[76] = 1.00000+power(10.0000, -algebraic[0]+constants[179])+(states[1]/constants[130])*(power(10.0000, constants[40])) algebraic[82] = (constants[181]+constants[180]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[180]) algebraic[134] = (constants[207]+constants[206]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(1.00000-constants[206]) algebraic[135] = (constants[209]+constants[208]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[208]) algebraic[167] = ((algebraic[137]+algebraic[135]+algebraic[89])-algebraic[134])-algebraic[82] algebraic[168] = (exp(-algebraic[167]/(constants[131]*constants[127]))*algebraic[83])/(algebraic[76]*(power(10.0000, -algebraic[0]))) algebraic[221] = (constants[224]*constants[86]*constants[87]*algebraic[168])/(constants[88]*constants[89]) algebraic[222] = ((algebraic[221]*states[13]*states[15])/(constants[86]*constants[87])-(constants[224]*states[12]*states[14])/(constants[88]*constants[89]))/algebraic[220] algebraic[223] = constants[0]*(algebraic[216]+algebraic[219]+algebraic[222]) algebraic[224] = constants[0]*-algebraic[222] algebraic[225] = 1.00000+states[3]/constants[93]+states[11]/constants[91]+states[15]/constants[92]+(states[11]*states[15])/(constants[91]*constants[92])+(states[11]*states[15]*states[3])/(constants[91]*constants[92]*constants[93])+states[16]/constants[94]+states[14]/constants[95]+(states[16]*states[14])/(constants[95]*constants[94]) algebraic[226] = constants[90]*0.000700000*exp(algebraic[0]*0.897900) algebraic[90] = 1.00000+power(10.0000, -algebraic[0]+constants[185]) algebraic[95] = (constants[187]+constants[186]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(16.0000-constants[186]) algebraic[164] = (((algebraic[95]+algebraic[135]+algebraic[137])-algebraic[8])-algebraic[75])-algebraic[134] algebraic[165] = (exp(-algebraic[164]/(constants[131]*constants[127]))*algebraic[90])/(algebraic[1]*algebraic[70]*(power(10.0000, -algebraic[0]))) algebraic[227] = (algebraic[226]*constants[94]*constants[95]*1.00000)/(constants[91]*constants[93]*constants[92]*algebraic[165]) algebraic[228] = ((algebraic[226]*states[11]*states[15]*states[3])/(constants[92]*constants[91]*constants[93])-(algebraic[227]*states[16]*states[14])/(constants[94]*constants[95]))/algebraic[225] algebraic[229] = constants[0]*((algebraic[216]-algebraic[219])-algebraic[228]) algebraic[96] = 1.00000+power(10.0000, -algebraic[0]+6.21000) algebraic[101] = (constants[190]+constants[189]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(9.00000-constants[189]) algebraic[170] = ((algebraic[16]+algebraic[101])-algebraic[95])-algebraic[24] algebraic[171] = (exp(-algebraic[170]/(constants[131]*constants[127]))*algebraic[9]*algebraic[96])/(algebraic[90]*algebraic[17]) algebraic[230] = ((constants[225]*constants[97]*constants[98])/(constants[99]*constants[100]))*algebraic[171] algebraic[231] = 1.00000+states[9]/constants[98]+states[16]/constants[97]+(states[16]*states[9])/(constants[97]*constants[98])+states[17]/constants[99]+states[8]/constants[100]+(states[17]*states[8])/(constants[99]*constants[100]) algebraic[232] = ((algebraic[230]*states[16]*states[9])/(constants[98]*constants[97])-(constants[225]*states[8]*states[17])/(constants[100]*constants[99]))/algebraic[231] algebraic[233] = constants[0]*(algebraic[228]-algebraic[232]) algebraic[234] = (constants[101]*0.989000)/(1.00000+power(10.0000, -algebraic[0]+5.62000)+power(10.0000, algebraic[0]-8.74000)) algebraic[102] = 1.00000+power(10.0000, -algebraic[0]+constants[191])+(states[1]/constants[130])*(power(10.0000, constants[42]))+(constants[129]/constants[130])*(power(10.0000, constants[192])) algebraic[109] = (constants[194]+constants[193]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(9.00000-constants[193]) algebraic[173] = algebraic[109]-algebraic[101] algebraic[174] = (exp(-algebraic[173]/(constants[131]*constants[127]))*algebraic[102])/algebraic[96] algebraic[235] = (algebraic[234]*constants[103])/(constants[102]*algebraic[174]) algebraic[236] = ((algebraic[234]*states[17])/constants[102]-(algebraic[235]*states[18])/constants[103])/(1.00000+states[17]/constants[102]+states[18]/constants[103]) algebraic[237] = constants[0]*(algebraic[232]-algebraic[236]) algebraic[110] = 1.00000+power(10.0000, constants[195]-algebraic[0])+(states[1]/constants[130])*(power(10.0000, constants[43]))+(constants[129]/constants[130])*(power(10.0000, constants[196])) algebraic[117] = (constants[198]+constants[197]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(9.00000-constants[197]) algebraic[136] = (constants[211]+constants[210]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(0.00000-constants[210]) algebraic[176] = (algebraic[136]+algebraic[117])-algebraic[109] algebraic[177] = (exp(-algebraic[176]/(constants[131]*constants[127]))*algebraic[110])/algebraic[102] algebraic[238] = (constants[226]*constants[106])/(constants[105]*algebraic[177]) algebraic[239] = ((constants[226]*states[18])/constants[105]-(algebraic[238]*states[19])/constants[106])/(1.00000+states[19]/constants[106]+states[18]/constants[105]) algebraic[240] = constants[0]*(algebraic[236]-algebraic[239]) algebraic[241] = (constants[107]*1.05000)/(1.00000+power(10.0000, -algebraic[0]+5.58000)+power(10.0000, algebraic[0]-8.79000)) algebraic[118] = 1.00000+power(10.0000, -algebraic[0]+constants[199]) algebraic[123] = (constants[201]+constants[200]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(1.00000-constants[200]) algebraic[179] = (((algebraic[123]+algebraic[16])-algebraic[137])-algebraic[117])-algebraic[24] algebraic[180] = (exp(-algebraic[179]/(constants[131]*constants[127]))*algebraic[118]*algebraic[9]*(power(10.0000, -algebraic[0])))/(algebraic[110]*algebraic[17]) algebraic[242] = (algebraic[241]*constants[110]*constants[111])/(constants[108]*constants[109]*algebraic[180]) algebraic[243] = ((algebraic[241]*states[19]*states[9])/(constants[108]*constants[109])-(algebraic[242]*states[20]*states[8])/(constants[110]*constants[111]))/(1.00000+states[19]/constants[108]+states[9]/constants[109]+(states[19]*states[9])/(constants[108]*constants[109])+states[8]/constants[111]+states[20]/constants[110]+(states[20]*states[8])/(constants[110]*constants[111])) algebraic[244] = constants[0]*(algebraic[239]-algebraic[243]) algebraic[245] = constants[112]*(-0.113400*algebraic[0]+1.60690) algebraic[124] = 1.00000+power(10.0000, -algebraic[0]+constants[202])+(states[1]/constants[130])*(power(10.0000, constants[45])) algebraic[132] = (constants[204]+constants[203]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(1.00000-constants[203]) algebraic[182] = (((algebraic[132]+algebraic[134])-algebraic[123])-algebraic[135])-algebraic[137] algebraic[183] = (exp(-algebraic[182]/(constants[131]*constants[127]))*algebraic[124]*(power(10.0000, -algebraic[0])))/algebraic[118] algebraic[246] = (algebraic[245]*constants[115]*constants[116])/(constants[113]*constants[114]*algebraic[183]) algebraic[247] = ((algebraic[245]*states[20]*states[14])/(constants[113]*constants[114])-(algebraic[246]*states[21]*states[15])/(constants[115]*constants[116]))/(1.00000+states[20]/constants[113]+states[14]/constants[114]+(states[20]*states[14])/(constants[113]*constants[114])+states[21]/constants[115]+states[15]/constants[116]+(states[21]*states[15])/(constants[115]*constants[116])) algebraic[248] = constants[0]*((-algebraic[228]-algebraic[222])+algebraic[247]) algebraic[249] = constants[0]*((algebraic[228]+algebraic[222])-algebraic[247]) algebraic[250] = constants[0]*(algebraic[243]-algebraic[247]) algebraic[32] = 1.00000+power(10.0000, -algebraic[0]+constants[160])+(states[1]/constants[130])*(power(10.0000, constants[141]))+(constants[129]/constants[130])*(power(10.0000, constants[159])) algebraic[39] = 1.00000+power(10.0000, -algebraic[0]+constants[162]) algebraic[140] = (exp(-constants[222]/(constants[131]*constants[127]))*(power(10.0000, -algebraic[0]))*algebraic[9]*algebraic[39])/(algebraic[32]*algebraic[17]) algebraic[251] = ((constants[117]/algebraic[140])*constants[119]*constants[121])/(constants[120]*constants[118]) algebraic[252] = ((algebraic[251]*states[8]*states[22])/(constants[119]*constants[121])-(constants[117]*states[9]*states[23])/(constants[120]*constants[118]))/(1.00000+states[9]/constants[120]+states[23]/constants[18]+(states[23]*states[9])/(constants[120]*constants[118])+states[8]/constants[119]+(states[22]*states[8])/(constants[121]*constants[119])) algebraic[253] = constants[0]*algebraic[252] algebraic[254] = constants[0]*-algebraic[252] algebraic[25] = 1.00000+power(10.0000, constants[156]-algebraic[0])+(states[1]/constants[130])*(power(10.0000, constants[140])) algebraic[31] = (constants[158]+constants[157]*log(10.0000)*constants[131]*constants[127]*algebraic[0])-constants[143]*(4.00000-constants[157]) algebraic[143] = (2.00000*algebraic[24]-algebraic[16])-algebraic[31] algebraic[144] = (exp(-algebraic[143]/(constants[131]*constants[127]))*(power(algebraic[17], 2.00000)))/(algebraic[9]*algebraic[25]) algebraic[255] = (constants[122]*(power(constants[125], 2.00000)))/(constants[123]*constants[124]*algebraic[144]) algebraic[256] = ((constants[122]*states[8]*states[5])/(constants[124]*constants[123])-algebraic[255]*(power(states[9]/constants[125], 2.00000)))/(1.00000+states[8]/constants[124]+states[5]/constants[123]+(states[8]*states[5])/(constants[124]*constants[123])+(2.00000*states[9])/constants[125]+(power(states[9], 2.00000))/(power(constants[125], 2.00000))) algebraic[257] = constants[0]*-algebraic[256] algebraic[258] = constants[24]*states[8] algebraic[259] = (constants[25]*(power(states[9]/constants[26], constants[27])))/(1.00000+power(states[9]/constants[26], constants[27])) algebraic[261] = constants[0]*(((((-algebraic[252]-algebraic[256])-algebraic[212])+algebraic[232]+algebraic[243])-algebraic[258])+algebraic[259]) algebraic[262] = constants[0]*(((((algebraic[252]+2.00000*algebraic[256]+algebraic[212])-algebraic[232])-algebraic[243])+algebraic[258])-algebraic[259]) algebraic[263] = constants[0]*(((-(algebraic[190]+algebraic[197])-algebraic[228])+algebraic[258])-algebraic[259]) algebraic[260] = (constants[28]*((states[21]*(power(10.0000, -algebraic[0])))/(constants[52]*constants[29])-(states[25]*(power(10.0000, -states[24])))/(constants[52]*constants[29])))/((states[21]*(power(10.0000, -algebraic[0])))/(constants[52]*constants[29])+2.00000+(power(10.0000, -algebraic[0]))/constants[52]+states[21]/constants[29]+states[25]/constants[29]+(power(10.0000, -states[24]))/constants[52]+(states[25]*(power(10.0000, -states[24])))/(constants[29]*constants[52])) algebraic[265] = constants[0]*(algebraic[247]-algebraic[260]) algebraic[266] = constants[0]*((constants[35]*algebraic[260])/constants[36]-(constants[34]*states[25])/constants[36]) algebraic[269] = constants[32]*(constants[33]*states[28]-constants[56]*(power(10.0000, -algebraic[0]))*states[29]) algebraic[272] = constants[0]*algebraic[269] algebraic[264] = algebraic[259]/constants[228] algebraic[267] = constants[227]*algebraic[264] algebraic[268] = constants[31]*(states[28]-states[26]) algebraic[273] = constants[0]*((algebraic[267]-algebraic[269])-algebraic[268]) algebraic[270] = 1.00000*(constants[33]*states[26]-constants[56]*(power(10.0000, -states[24]))*states[27]) algebraic[274] = constants[0]*((constants[34]*(constants[48]-states[27]))/constants[36]+algebraic[270]) algebraic[275] = constants[0]*(((constants[34]*(constants[1]-states[26]))/constants[36]+(algebraic[268]*constants[35])/constants[36])-algebraic[270]) algebraic[276] = 1.00000e+07*((power(10.0000, -7.18400))*states[30]-states[31]*(power(10.0000, -states[24]))) algebraic[278] = constants[0]*((constants[34]*(constants[53]-states[30]))/constants[36]-algebraic[276]) algebraic[279] = constants[0]*((constants[34]*(constants[54]-states[31]))/constants[36]+algebraic[276]) algebraic[280] = (constants[0]*((constants[34]*(power(10.0000, -constants[55])-power(10.0000, -states[24])))/constants[36]+algebraic[270]/1.00000+algebraic[276]/1.00000+(constants[35]*algebraic[260])/(constants[36]*1.00000)))/(-log(10.0000)*(power(10.0000, -states[24]))) algebraic[2] = 1.00000/algebraic[1] algebraic[3] = (power(10.0000, -algebraic[0]+constants[144]))*algebraic[2] algebraic[7] = ((power(10.0000, constants[144]))*algebraic[1]*(1.00000-algebraic[3]))/(constants[130]*(power(algebraic[1], 2.00000))) algebraic[10] = 1.00000/algebraic[9] algebraic[13] = (power(10.0000, -algebraic[0]+constants[148]))*algebraic[10] algebraic[15] = ((power(10.0000, constants[148]))*algebraic[9]*(1.00000-algebraic[13]))/(constants[130]*(power(algebraic[9], 2.00000))) algebraic[18] = 1.00000/algebraic[17] algebraic[21] = (power(10.0000, -algebraic[0]+constants[153]))*algebraic[18] algebraic[23] = ((power(10.0000, constants[153]))*algebraic[17]*(1.00000-algebraic[21]))/(constants[130]*(power(algebraic[17], 2.00000))) algebraic[30] = ((power(10.0000, constants[156]))*(algebraic[25]-power(10.0000, -algebraic[0]+constants[156])))/(constants[130]*(power(algebraic[25], 2.00000))) algebraic[38] = ((power(10.0000, constants[160]))*(algebraic[32]-power(10.0000, -algebraic[0]+constants[160])))/(constants[130]*(power(algebraic[32], 2.00000))) algebraic[43] = ((power(10.0000, constants[162]))*(algebraic[39]-power(10.0000, -algebraic[0]+constants[162])))/(constants[130]*(power(algebraic[39], 2.00000))) algebraic[49] = ((power(10.0000, constants[164]))*(algebraic[44]-power(10.0000, -algebraic[0]+constants[164])))/(constants[130]*(power(algebraic[44], 2.00000))) algebraic[55] = ((power(10.0000, constants[167]))*(algebraic[51]-power(10.0000, -algebraic[0]+constants[167])))/(constants[130]*(power(algebraic[51], 2.00000))) algebraic[61] = ((power(10.0000, constants[170]))*(algebraic[57]-power(10.0000, -algebraic[0]+constants[170])))/(constants[130]*(power(algebraic[57], 2.00000))) algebraic[64] = 1.00000/algebraic[63] algebraic[66] = algebraic[64]*(power(10.0000, -algebraic[0]+constants[173])) algebraic[68] = ((power(10.0000, constants[173]))*algebraic[63]*(1.00000-algebraic[66]))/(constants[130]*(power(algebraic[63], 2.00000))) algebraic[74] = ((power(10.0000, constants[176]))*(algebraic[70]-power(10.0000, -algebraic[0]+constants[176])))/(constants[130]*(power(algebraic[70], 2.00000))) algebraic[81] = ((power(10.0000, constants[179]))*(algebraic[76]-power(10.0000, -algebraic[0]+constants[179])))/(constants[130]*(power(algebraic[76], 2.00000))) algebraic[88] = ((power(10.0000, constants[182]))*(algebraic[83]-power(10.0000, -algebraic[0]+constants[182])))/(constants[130]*(power(algebraic[83], 2.00000))) algebraic[94] = ((power(10.0000, constants[185]))*(algebraic[90]-power(10.0000, -algebraic[0]+constants[185])))/(constants[130]*(power(algebraic[90], 2.00000))) algebraic[100] = ((power(10.0000, constants[188]))*(algebraic[96]-power(10.0000, -algebraic[0]+constants[188])))/(constants[130]*(power(algebraic[96], 2.00000))) algebraic[108] = ((power(10.0000, constants[191]))*(algebraic[102]-power(10.0000, -algebraic[0]+constants[191])))/(constants[130]*(power(algebraic[102], 2.00000))) algebraic[116] = ((power(10.0000, constants[195]))*(algebraic[110]-power(10.0000, -algebraic[0]+constants[195])))/(constants[130]*(power(algebraic[110], 2.00000))) algebraic[122] = ((power(10.0000, constants[199]))*(algebraic[118]-power(10.0000, -algebraic[0]+constants[199])))/(constants[130]*(power(algebraic[118], 2.00000))) algebraic[131] = ((power(10.0000, constants[202]))*(algebraic[124]-power(10.0000, -algebraic[0]+constants[202])))/(constants[130]*(power(algebraic[124], 2.00000))) algebraic[271] = log(10.0000)*(power(10.0000, -algebraic[0]))*constants[130]*(1.00000+algebraic[7]*states[3]+algebraic[15]*states[8]+algebraic[23]*states[9]+algebraic[30]*states[5]+algebraic[38]*states[23]+algebraic[43]*states[22]+algebraic[49]*states[4]+algebraic[55]*states[6]+algebraic[61]*states[7]+algebraic[68]*states[10]+algebraic[74]*states[11]+algebraic[88]*states[12]+algebraic[81]*states[13]+algebraic[94]*states[16]+algebraic[100]*states[17]+algebraic[108]*states[18]+algebraic[116]*states[19]+algebraic[122]*states[20]+algebraic[131]*states[21]) algebraic[14] = algebraic[13] algebraic[22] = algebraic[21] algebraic[26] = 1.00000/algebraic[25] algebraic[28] = algebraic[26]*(power(10.0000, constants[156]-algebraic[0])) algebraic[29] = algebraic[28] algebraic[142] = ((2.00000*algebraic[22]-algebraic[14])-algebraic[29])+((2.00000*constants[154]-constants[150])-constants[157]) algebraic[6] = algebraic[3] algebraic[184] = ((algebraic[22]+algebraic[6])-algebraic[14])+(((constants[154]+constants[146])-constants[150])-constants[210]) algebraic[33] = 1.00000/algebraic[32] algebraic[36] = (power(10.0000, -algebraic[0]+constants[160]))*algebraic[33] algebraic[37] = algebraic[36] algebraic[40] = 1.00000/algebraic[39] algebraic[41] = algebraic[40]*(power(10.0000, -algebraic[0]+constants[162])) algebraic[42] = 0.00000*algebraic[40]+1.00000*algebraic[41] algebraic[138] = (((algebraic[14]+algebraic[42])-algebraic[37])-algebraic[22])+(((constants[150]+constants[163])-constants[161])-constants[154]) algebraic[277] = algebraic[138]*-algebraic[252] algebraic[45] = 1.00000/algebraic[44] algebraic[46] = algebraic[45]*(power(10.0000, -algebraic[0]+constants[164])) algebraic[48] = algebraic[46] algebraic[145] = (algebraic[48]-algebraic[6])+(1.00000-constants[146]) algebraic[52] = 1.00000/algebraic[51] algebraic[53] = algebraic[52]*(power(10.0000, -algebraic[0]+constants[167])) algebraic[54] = algebraic[53] algebraic[148] = (algebraic[54]-algebraic[48])+(constants[168]-constants[165]) algebraic[58] = 1.00000/algebraic[57] algebraic[59] = algebraic[58]*(power(10.0000, -algebraic[0]+constants[170])) algebraic[60] = algebraic[59] algebraic[151] = (algebraic[60]-algebraic[54])+(constants[171]-constants[168]) algebraic[67] = algebraic[66] algebraic[154] = (((((algebraic[22]+algebraic[67])-algebraic[60])-algebraic[14])+constants[154]+constants[174])-constants[171])-constants[150] algebraic[71] = 1.00000/algebraic[70] algebraic[72] = algebraic[71]*(power(10.0000, -algebraic[0]+constants[176])) algebraic[73] = algebraic[72] algebraic[84] = 1.00000/algebraic[83] algebraic[86] = algebraic[84]*(power(10.0000, -algebraic[0]+constants[182])) algebraic[87] = algebraic[86] algebraic[157] = ((algebraic[87]+algebraic[73])-algebraic[67])+((constants[183]+constants[177])-constants[174]) algebraic[160] = (algebraic[87]-algebraic[73])+(constants[183]-constants[177]) algebraic[91] = 1.00000/algebraic[90] algebraic[92] = algebraic[91]*(power(10.0000, -algebraic[0]+constants[185])) algebraic[93] = algebraic[92] algebraic[163] = ((algebraic[93]-algebraic[73])-algebraic[6])+((((constants[186]+constants[208])-constants[177])-constants[146])-constants[206]) algebraic[77] = 1.00000/algebraic[76] algebraic[79] = algebraic[77]*(power(10.0000, -algebraic[0]+constants[179])) algebraic[80] = algebraic[79] algebraic[166] = (algebraic[87]-algebraic[80])+(((constants[183]+constants[208])-constants[206])-constants[180]) algebraic[97] = 1.00000/algebraic[96] algebraic[98] = algebraic[97]*(power(10.0000, -algebraic[0]+6.21000)) algebraic[99] = algebraic[98] algebraic[169] = (((algebraic[99]+algebraic[14])-algebraic[93])-algebraic[22])+(((constants[189]+constants[150])-constants[186])-constants[154]) algebraic[103] = 1.00000/algebraic[102] algebraic[104] = algebraic[103]*(power(10.0000, -algebraic[0]+constants[191])) algebraic[107] = algebraic[104] algebraic[172] = (algebraic[107]-algebraic[99])+(constants[193]-constants[189]) algebraic[111] = 1.00000/algebraic[110] algebraic[112] = algebraic[111]*(power(10.0000, constants[195]-algebraic[0])) algebraic[115] = algebraic[112] algebraic[175] = (algebraic[115]-algebraic[107])+((constants[210]+constants[197])-constants[193]) algebraic[119] = 1.00000/algebraic[118] algebraic[120] = algebraic[119]*(power(10.0000, -algebraic[0]+constants[199])) algebraic[121] = algebraic[120] algebraic[178] = (((algebraic[121]+algebraic[14])-algebraic[115])-algebraic[22])+(((constants[200]+constants[150])-constants[197])-constants[154]) algebraic[127] = 1.00000/algebraic[124] algebraic[129] = algebraic[127]*(power(10.0000, -algebraic[0]+constants[202])) algebraic[130] = algebraic[129] algebraic[181] = (algebraic[130]-algebraic[121])+(((constants[203]+constants[206])-constants[200])-constants[208]) algebraic[281] = algebraic[148]*algebraic[201]+algebraic[145]*(algebraic[190]+algebraic[197])+algebraic[151]*algebraic[205]+algebraic[154]*algebraic[212]+algebraic[157]*algebraic[216]+algebraic[160]*algebraic[219]+algebraic[163]*algebraic[228]+algebraic[169]*algebraic[232]+algebraic[172]*algebraic[236]+algebraic[175]*algebraic[239]+algebraic[178]*algebraic[243]+algebraic[181]*algebraic[247]+algebraic[166]*algebraic[222] algebraic[282] = ((algebraic[277]+algebraic[281]+algebraic[184]*algebraic[258]+-algebraic[184]*algebraic[259]+algebraic[142]*algebraic[256])-algebraic[269])+algebraic[260] algebraic[283] = algebraic[282]/(constants[126]+algebraic[271]) algebraic[284] = constants[0]*algebraic[283] algebraic[4] = (constants[129]/constants[130])*(power(10.0000, constants[145]))*algebraic[2] algebraic[5] = (states[1]/constants[130])*(power(10.0000, constants[136]))*algebraic[2] algebraic[11] = (states[1]/constants[130])*(power(10.0000, constants[138]))*algebraic[10] algebraic[12] = (constants[129]/constants[130])*(power(10.0000, constants[149]))*algebraic[10] algebraic[19] = ((algebraic[18]*states[1])/constants[130])*(power(10.0000, constants[139])) algebraic[20] = ((algebraic[18]*constants[129])/constants[130])*(power(10.0000, constants[152])) algebraic[27] = (states[1]/constants[130])*(power(10.0000, constants[140]))*algebraic[26] algebraic[34] = (constants[129]/constants[130])*(power(10.0000, constants[159]))*algebraic[33] algebraic[35] = (states[1]/constants[130])*(power(10.0000, constants[141]))*algebraic[33] algebraic[47] = ((algebraic[45]*states[1])/constants[130])*(power(10.0000, constants[142])) algebraic[65] = ((algebraic[64]*states[1])/constants[130])*(power(10.0000, constants[38])) algebraic[78] = ((algebraic[77]*states[1])/constants[130])*(power(10.0000, constants[40])) algebraic[85] = ((algebraic[84]*states[1])/constants[130])*(power(10.0000, constants[41])) algebraic[105] = ((algebraic[103]*states[1])/constants[130])*(power(10.0000, constants[42])) algebraic[106] = ((algebraic[103]*constants[129])/constants[130])*(power(10.0000, constants[192])) algebraic[113] = ((algebraic[111]*constants[129])/constants[130])*(power(10.0000, constants[196])) algebraic[114] = ((algebraic[111]*states[1])/constants[130])*(power(10.0000, constants[43])) algebraic[125] = ((algebraic[124]*(power(10.0000, constants[45])))/constants[130]-((states[1]/constants[130])*(power(10.0000, 2.00000*constants[45])))/constants[130])/(power(algebraic[124], 2.00000)) algebraic[126] = ((states[1]/constants[130])*(power(10.0000, -algebraic[0]+constants[202]+constants[45]))*log(10.0000))/(power(algebraic[124], 2.00000)) algebraic[128] = (states[1]/constants[130])*(power(10.0000, constants[45]))*algebraic[127] algebraic[139] = ((algebraic[24]+algebraic[8]+algebraic[137])-algebraic[136])-algebraic[16] algebraic[141] = (exp(-algebraic[139]/(constants[131]*constants[127]))*algebraic[17]*algebraic[1])/(algebraic[9]*(power(10.0000, -algebraic[0]))) return algebraic def solve_model(): """Solve model with ODE solver""" from scipy.integrate import ode # Initialise constants and state variables (init_states, constants) = initConsts() # Set timespan to solve over voi = linspace(0, 10, 500) # Construct ODE object to solve r = ode(computeRates) r.set_integrator('vode', method='bdf', atol=1e-06, rtol=1e-06, max_step=1) r.set_initial_value(init_states, voi[0]) r.set_f_params(constants) # Solve model states = array([[0.0] * len(voi)] * sizeStates) states[:,0] = init_states for (i,t) in enumerate(voi[1:]): if r.successful(): r.integrate(t) states[:,i+1] = r.y else: break # Compute algebraic variables algebraic = computeAlgebraic(constants, states, voi) return (voi, states, algebraic) def plot_model(voi, states, algebraic): """Plot variables against variable of integration""" import pylab (legend_states, legend_algebraic, legend_voi, legend_constants) = createLegends() pylab.figure(1) pylab.plot(voi,vstack((states,algebraic)).T) pylab.xlabel(legend_voi) pylab.legend(legend_states + legend_algebraic, loc='best') pylab.show() if __name__ == "__main__": (voi, states, algebraic) = solve_model() plot_model(voi, states, algebraic)