Generated Code

The following is python code generated by the CellML API from this CellML file. (Back to language selection)

The raw code is available.

# 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)