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 = 207
sizeStates = 60
sizeConstants = 388
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 (second)"
legend_constants[0] = "ExpType in component Environment (dimensionless)"
legend_constants[1] = "StateType in component Environment (dimensionless)"
legend_constants[2] = "F in component parameters (kilojoule_per_mole_per_millivolt)"
legend_constants[323] = "V_cyto in component parameters (cyto_per_cell)"
legend_constants[307] = "V_mito in component parameters (mito_per_cell)"
legend_constants[3] = "W_x in component parameters (l_water_per_l_mito)"
legend_constants[4] = "W_i in component parameters (l_water_per_l_mito)"
legend_constants[286] = "W_c in component parameters (l_water_per_l_cyto)"
legend_constants[5] = "gamma in component parameters (per_micrometer)"
legend_constants[6] = "rho_m in component parameters (mg_per_mitochondria)"
legend_constants[7] = "n_A in component parameters (dimensionless)"
legend_constants[8] = "p_PI in component parameters (micrometer_per_second)"
legend_constants[9] = "p_A in component parameters (micrometer_per_second)"
legend_constants[10] = "p_TI in component parameters (micrometer_per_second)"
legend_constants[11] = "k_mADP in component parameters (molar)"
legend_constants[12] = "theta in component parameters (dimensionless)"
legend_constants[13] = "k_O2 in component parameters (molar)"
legend_constants[14] = "NAD_tot in component parameters (molar)"
legend_constants[15] = "Q_tot in component parameters (molar)"
legend_constants[16] = "cytC_tot in component parameters (molar)"
legend_constants[17] = "A_tot in component parameters (molar)"
legend_constants[18] = "CR_tot in component parameters (molar)"
legend_constants[19] = "CO2_tot_x in component parameters (molar)"
legend_constants[20] = "FAD_tot in component parameters (molar)"
legend_constants[21] = "C_IM in component parameters (mole_per_l_mito_volume_per_mv)"
legend_constants[22] = "X_pdh in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[23] = "X_cits in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[24] = "X_acon in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[25] = "X_isod in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[26] = "X_akgd in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[27] = "K_ir_akgd in component parameters (molar)"
legend_constants[28] = "X_scoas in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[29] = "X_sdh in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[30] = "X_fum in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[31] = "X_mdh in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[32] = "X_ndk in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[33] = "X_got in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[34] = "X_PYRH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[35] = "X_GLUH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[36] = "X_CITMAL in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[37] = "X_AKGMAL in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[38] = "X_SUCMAL in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[39] = "X_MALPI in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[40] = "X_ASPGLU in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[41] = "X_C1 in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[42] = "X_C3 in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_half_molar)"
legend_constants[43] = "k_PI_1 in component parameters (molar)"
legend_constants[44] = "k_PI_2 in component parameters (molar)"
legend_constants[45] = "X_C4 in component parameters (mole_per_second_per_l_mito_volume_per_molar)"
legend_constants[46] = "X_F1 in component parameters (mole_per_second_per_l_mito_volume_per_molar)"
legend_constants[47] = "X_ANT in component parameters (mole_per_second_per_l_mito_volume)"
legend_constants[48] = "X_PiHt in component parameters (mole_per_second_per_l_mito_volume_per_molar)"
legend_constants[49] = "k_PiHt in component parameters (molar)"
legend_constants[50] = "X_KH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[51] = "X_Hle in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_mv)"
legend_constants[273] = "X_HK in component parameters (mole_per_second_per_l_cyto_volume)"
legend_constants[52] = "J_AtC in component parameters (mole_per_second_per_l_cyto_volume)"
legend_constants[274] = "RT in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[53] = "dGf1_H2O in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[54] = "dGf1_O2 in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[55] = "dGf1_NADH in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[56] = "dGf1_NAD in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[57] = "dGf1_QH2 in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[58] = "dGf1_COQ in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[59] = "dGf1_ATP in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[60] = "dGf1_ADP in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[61] = "dGf1_AMP in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[62] = "dGf1_GTP in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[63] = "dGf1_GDP in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[64] = "dGf1_Cred in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[65] = "dGf1_Cox in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[66] = "dGf1_Pi in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[67] = "dGf1_PCr in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[68] = "dGf1_Cr in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[69] = "dGf1_FADH2 in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[70] = "dGf1_FAD in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[71] = "dGf1_COASH in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[72] = "dGf1_ACCOA in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[73] = "dGf1_OAA in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[74] = "dGf1_CIT in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[75] = "dGf1_ICIT in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[76] = "dGf1_AKG in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[77] = "dGf1_SCOA in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[78] = "dGf1_SUC in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[79] = "dGf1_FUM in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[80] = "dGf1_MAL in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[81] = "dGf1_ASP in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[82] = "dGf1_GLU in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[83] = "dGf1_CO2_tot in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[84] = "dGf1_PYR in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[85] = "dGf1_GLC in component thermodynamic_data (kilojoule_per_mole)"
legend_constants[86] = "dGf1_G6P in component thermodynamic_data (kilojoule_per_mole)"
legend_algebraic[0] = "P_ATP_x in component P_ATP (dimensionless)"
legend_algebraic[163] = "P_ATP_i in component P_ATP (dimensionless)"
legend_algebraic[17] = "P_ATP_c in component P_ATP (dimensionless)"
legend_algebraic[18] = "uATP_x in component P_ATP (molar)"
legend_algebraic[20] = "mATP_c in component P_ATP (molar)"
legend_constants[296] = "K_H in component P_ATP (molar)"
legend_constants[311] = "K_Mg in component P_ATP (molar)"
legend_constants[325] = "K_K in component P_ATP (molar)"
legend_states[0] = "ATP_x in component dATP_x_dt (molar)"
legend_states[1] = "H_x in component dH_x_dt (molar)"
legend_states[2] = "Mg_x in component dMg_x_dt (molar)"
legend_states[3] = "K_x in component dK_x_dt (molar)"
legend_algebraic[132] = "H_i in component dH_i_dt (molar)"
legend_algebraic[149] = "Mg_i in component dMg_i_dt (molar)"
legend_algebraic[162] = "K_i in component dK_i_dt (molar)"
legend_states[4] = "ATP_c in component dATP_c_dt (molar)"
legend_states[5] = "H_c in component dH_c_dt (molar)"
legend_states[6] = "Mg_c in component dMg_c_dt (molar)"
legend_states[7] = "K_c in component dK_c_dt (molar)"
legend_algebraic[21] = "P_ADP_x in component P_ADP (dimensionless)"
legend_algebraic[165] = "P_ADP_i in component P_ADP (dimensionless)"
legend_algebraic[22] = "P_ADP_c in component P_ADP (dimensionless)"
legend_algebraic[23] = "uADP_x in component P_ADP (molar)"
legend_algebraic[25] = "mADP_c in component P_ADP (molar)"
legend_constants[332] = "K_H in component P_ADP (molar)"
legend_constants[336] = "K_Mg in component P_ADP (molar)"
legend_constants[340] = "K_K in component P_ADP (molar)"
legend_states[8] = "ADP_x in component dADP_x_dt (molar)"
legend_states[9] = "ADP_c in component dADP_c_dt (molar)"
legend_algebraic[26] = "P_AMP_x in component P_AMP (dimensionless)"
legend_algebraic[167] = "P_AMP_i in component P_AMP (dimensionless)"
legend_algebraic[27] = "P_AMP_c in component P_AMP (dimensionless)"
legend_algebraic[28] = "uAMP_x in component P_AMP (molar)"
legend_constants[344] = "K_H in component P_AMP (molar)"
legend_constants[347] = "K_Mg in component P_AMP (molar)"
legend_constants[349] = "K_K in component P_AMP (molar)"
legend_states[10] = "AMP_x in component dAMP_x_dt (molar)"
legend_algebraic[29] = "P_GTP_x in component P_GTP (dimensionless)"
legend_algebraic[164] = "P_GTP_i in component P_GTP (dimensionless)"
legend_algebraic[19] = "P_GTP_c in component P_GTP (dimensionless)"
legend_algebraic[30] = "uGTP_x in component P_GTP (molar)"
legend_constants[313] = "K_H in component P_GTP (molar)"
legend_constants[326] = "K_Mg in component P_GTP (molar)"
legend_constants[333] = "K_K in component P_GTP (molar)"
legend_states[11] = "GTP_x in component dGTP_x_dt (molar)"
legend_algebraic[31] = "P_GDP_x in component P_GDP (dimensionless)"
legend_algebraic[166] = "P_GDP_i in component P_GDP (dimensionless)"
legend_algebraic[24] = "P_GDP_c in component P_GDP (dimensionless)"
legend_algebraic[32] = "uGDP_x in component P_GDP (molar)"
legend_constants[337] = "K_H in component P_GDP (molar)"
legend_constants[341] = "K_Mg in component P_GDP (molar)"
legend_constants[345] = "K_K in component P_GDP (molar)"
legend_states[12] = "GDP_x in component dGDP_x_dt (molar)"
legend_algebraic[33] = "P_PI_x in component P_PI (dimensionless)"
legend_algebraic[168] = "P_PI_i in component P_PI (dimensionless)"
legend_algebraic[1] = "P_PI_c in component P_PI (dimensionless)"
legend_algebraic[34] = "hPI_x in component P_PI (molar)"
legend_algebraic[169] = "hPI_i in component P_PI (molar)"
legend_states[13] = "PI_x in component dPI_x_dt (molar)"
legend_states[14] = "PI_i in component dPI_i_dt (molar)"
legend_constants[351] = "K_H in component P_PI (molar)"
legend_constants[353] = "K_Mg in component P_PI (molar)"
legend_constants[355] = "K_K in component P_PI (molar)"
legend_algebraic[35] = "P_COASH_x in component P_COASH (dimensionless)"
legend_algebraic[133] = "P_COASH_i in component P_COASH (dimensionless)"
legend_algebraic[2] = "P_COASH_c in component P_COASH (dimensionless)"
legend_constants[357] = "K_H in component P_COASH (molar)"
legend_constants[87] = "K_Mg in component P_COASH (molar)"
legend_constants[88] = "K_K in component P_COASH (molar)"
legend_algebraic[36] = "P_OAA_x in component P_OAA (dimensionless)"
legend_algebraic[150] = "P_OAA_i in component P_OAA (dimensionless)"
legend_algebraic[3] = "P_OAA_c in component P_OAA (dimensionless)"
legend_constants[359] = "K_Mg in component P_OAA (molar)"
legend_constants[89] = "K_H in component P_OAA (molar)"
legend_constants[90] = "K_K in component P_OAA (molar)"
legend_algebraic[37] = "P_CIT_x in component P_CIT (dimensionless)"
legend_algebraic[170] = "P_CIT_i in component P_CIT (dimensionless)"
legend_algebraic[4] = "P_CIT_c in component P_CIT (dimensionless)"
legend_algebraic[38] = "uCIT_x in component P_CIT (molar)"
legend_algebraic[171] = "uCIT_i in component P_CIT (molar)"
legend_algebraic[39] = "hCIT_x in component P_CIT (molar)"
legend_algebraic[172] = "hCIT_i in component P_CIT (molar)"
legend_states[15] = "CIT_x in component dCIT_x_dt (molar)"
legend_states[16] = "CIT_i in component dCIT_i_dt (molar)"
legend_constants[360] = "K_H in component P_CIT (molar)"
legend_constants[361] = "K_Mg in component P_CIT (molar)"
legend_constants[362] = "K_K in component P_CIT (molar)"
legend_algebraic[40] = "P_ICIT_x in component P_ICIT (dimensionless)"
legend_algebraic[151] = "P_ICIT_i in component P_ICIT (dimensionless)"
legend_algebraic[5] = "P_ICIT_c in component P_ICIT (dimensionless)"
legend_constants[363] = "K_H in component P_ICIT (molar)"
legend_constants[364] = "K_Mg in component P_ICIT (molar)"
legend_constants[91] = "K_K in component P_ICIT (molar)"
legend_algebraic[41] = "P_SCOA_x in component P_SCOA (dimensionless)"
legend_algebraic[134] = "P_SCOA_i in component P_SCOA (dimensionless)"
legend_algebraic[6] = "P_SCOA_c in component P_SCOA (dimensionless)"
legend_constants[365] = "K_H in component P_SCOA (molar)"
legend_constants[92] = "K_Mg in component P_SCOA (molar)"
legend_constants[93] = "K_K in component P_SCOA (molar)"
legend_algebraic[42] = "P_SUC_x in component P_SUC (dimensionless)"
legend_algebraic[173] = "P_SUC_i in component P_SUC (dimensionless)"
legend_algebraic[7] = "P_SUC_c in component P_SUC (dimensionless)"
legend_constants[366] = "K_H in component P_SUC (molar)"
legend_constants[367] = "K_Mg in component P_SUC (molar)"
legend_constants[368] = "K_K in component P_SUC (molar)"
legend_algebraic[43] = "P_FUM_x in component P_FUM (dimensionless)"
legend_algebraic[135] = "P_FUM_i in component P_FUM (dimensionless)"
legend_algebraic[8] = "P_FUM_c in component P_FUM (dimensionless)"
legend_constants[369] = "K_H in component P_FUM (molar)"
legend_constants[94] = "K_Mg in component P_FUM (molar)"
legend_constants[95] = "K_K in component P_FUM (molar)"
legend_algebraic[44] = "P_MAL_x in component P_MAL (dimensionless)"
legend_algebraic[174] = "P_MAL_i in component P_MAL (dimensionless)"
legend_algebraic[9] = "P_MAL_c in component P_MAL (dimensionless)"
legend_constants[370] = "K_H in component P_MAL (molar)"
legend_constants[371] = "K_Mg in component P_MAL (molar)"
legend_constants[372] = "K_K in component P_MAL (molar)"
legend_algebraic[45] = "P_CO2_tot_x in component P_CO2_tot (dimensionless)"
legend_algebraic[136] = "P_CO2_tot_i in component P_CO2_tot (dimensionless)"
legend_algebraic[10] = "P_CO2_tot_c in component P_CO2_tot (dimensionless)"
legend_constants[373] = "K_H in component P_CO2_tot (molar)"
legend_constants[96] = "K_Mg in component P_CO2_tot (molar)"
legend_constants[97] = "K_K in component P_CO2_tot (molar)"
legend_algebraic[46] = "P_PYR_x in component P_PYR (dimensionless)"
legend_algebraic[152] = "P_PYR_i in component P_PYR (dimensionless)"
legend_algebraic[11] = "P_PYR_c in component P_PYR (dimensionless)"
legend_constants[374] = "K_Mg in component P_PYR (molar)"
legend_constants[98] = "K_H in component P_PYR (molar)"
legend_constants[99] = "K_K in component P_PYR (molar)"
legend_algebraic[12] = "P_G6P_x in component P_G6P (dimensionless)"
legend_algebraic[137] = "P_G6P_i in component P_G6P (dimensionless)"
legend_algebraic[47] = "P_G6P_c in component P_G6P (dimensionless)"
legend_constants[375] = "K_H in component P_G6P (molar)"
legend_constants[100] = "K_Mg in component P_G6P (molar)"
legend_constants[101] = "K_K in component P_G6P (molar)"
legend_algebraic[48] = "P_GLU_x in component P_GLU (dimensionless)"
legend_algebraic[153] = "P_GLU_i in component P_GLU (dimensionless)"
legend_algebraic[13] = "P_GLU_c in component P_GLU (dimensionless)"
legend_constants[376] = "K_H in component P_GLU (molar)"
legend_constants[377] = "K_Mg in component P_GLU (molar)"
legend_constants[102] = "K_K in component P_GLU (molar)"
legend_algebraic[49] = "P_ASP_x in component P_ASP (dimensionless)"
legend_algebraic[154] = "P_ASP_i in component P_ASP (dimensionless)"
legend_algebraic[14] = "P_ASP_c in component P_ASP (dimensionless)"
legend_constants[378] = "K_H in component P_ASP (molar)"
legend_constants[379] = "K_Mg in component P_ASP (molar)"
legend_constants[103] = "K_K in component P_ASP (molar)"
legend_constants[275] = "P_GLC_c in component P_GLC (dimensionless)"
legend_constants[104] = "K_H in component P_GLC (molar)"
legend_constants[105] = "K_Mg in component P_GLC (molar)"
legend_constants[106] = "K_K in component P_GLC (molar)"
legend_constants[276] = "P_NAD_x in component P_NAD (dimensionless)"
legend_constants[107] = "K_H in component P_NAD (molar)"
legend_constants[108] = "K_Mg in component P_NAD (molar)"
legend_constants[109] = "K_K in component P_NAD (molar)"
legend_constants[277] = "P_NADH_x in component P_NADH (dimensionless)"
legend_constants[110] = "K_H in component P_NADH (molar)"
legend_constants[111] = "K_Mg in component P_NADH (molar)"
legend_constants[112] = "K_K in component P_NADH (molar)"
legend_constants[304] = "P_QH2_x in component P_QH2 (dimensionless)"
legend_constants[113] = "K_H in component P_QH2 (molar)"
legend_constants[114] = "K_Mg in component P_QH2 (molar)"
legend_constants[115] = "K_K in component P_QH2 (molar)"
legend_constants[320] = "P_COQ_x in component P_COQ (dimensionless)"
legend_constants[116] = "K_H in component P_COQ (molar)"
legend_constants[117] = "K_Mg in component P_COQ (molar)"
legend_constants[118] = "K_K in component P_COQ (molar)"
legend_constants[278] = "P_ACCOA_x in component P_ACCOA (dimensionless)"
legend_constants[119] = "K_H in component P_ACCOA (molar)"
legend_constants[120] = "K_Mg in component P_ACCOA (molar)"
legend_constants[121] = "K_K in component P_ACCOA (molar)"
legend_constants[306] = "P_AKG_x in component P_AKG (dimensionless)"
legend_constants[122] = "K_H in component P_AKG (molar)"
legend_constants[123] = "K_Mg in component P_AKG (molar)"
legend_constants[124] = "K_K in component P_AKG (molar)"
legend_constants[279] = "P_O2_x in component P_O2 (dimensionless)"
legend_constants[125] = "K_H in component P_O2 (molar)"
legend_constants[126] = "K_Mg in component P_O2 (molar)"
legend_constants[127] = "K_K in component P_O2 (molar)"
legend_constants[280] = "P_FADH2_x in component P_FADH2 (dimensionless)"
legend_constants[128] = "K_H in component P_FADH2 (molar)"
legend_constants[129] = "K_Mg in component P_FADH2 (molar)"
legend_constants[130] = "K_K in component P_FADH2 (molar)"
legend_constants[281] = "P_FAD_x in component P_FAD (dimensionless)"
legend_constants[131] = "K_H in component P_FAD (molar)"
legend_constants[132] = "K_Mg in component P_FAD (molar)"
legend_constants[133] = "K_K in component P_FAD (molar)"
legend_constants[282] = "P_H_x in component P_H (dimensionless)"
legend_constants[134] = "K_H in component P_H (molar)"
legend_constants[135] = "K_Mg in component P_H (molar)"
legend_constants[136] = "K_K in component P_H (molar)"
legend_constants[283] = "P_Mg_x in component P_Mg (dimensionless)"
legend_constants[137] = "K_H in component P_Mg (molar)"
legend_constants[138] = "K_Mg in component P_Mg (molar)"
legend_constants[139] = "K_K in component P_Mg (molar)"
legend_constants[284] = "P_K_x in component P_K (dimensionless)"
legend_constants[140] = "K_H in component P_K (molar)"
legend_constants[141] = "K_Mg in component P_K (molar)"
legend_constants[142] = "K_K in component P_K (molar)"
legend_algebraic[114] = "J_pdh in component J_pdh (mole_per_second_per_l_mito_volume)"
legend_states[17] = "PYR_x in component dPYR_x_dt (molar)"
legend_states[18] = "COASH_x in component dCOASH_x_dt (molar)"
legend_algebraic[113] = "NAD_x in component dNAD_x_dt (molar)"
legend_states[19] = "CO2_tot_x in component dCO2_tot_x_dt (molar)"
legend_states[20] = "ACCOA_x in component dACCOA_x_dt (molar)"
legend_states[21] = "NADH_x in component dNADH_x_dt (molar)"
legend_algebraic[50] = "K_eq_pdh in component J_pdh (dimensionless)"
legend_constants[314] = "K0_eq_pdh in component J_pdh (molar)"
legend_constants[298] = "dGr0_pdh in component J_pdh (kilojoule_per_mole)"
legend_constants[143] = "K_mA in component J_pdh (molar)"
legend_constants[144] = "K_mB in component J_pdh (molar)"
legend_constants[145] = "K_mC in component J_pdh (molar)"
legend_algebraic[51] = "alpha_i1 in component J_pdh (dimensionless)"
legend_algebraic[52] = "alpha_i2 in component J_pdh (dimensionless)"
legend_constants[146] = "K_iACCOA in component J_pdh (molar)"
legend_constants[147] = "K_iNADH in component J_pdh (molar)"
legend_constants[148] = "minCon in component J_pdh (molar)"
legend_algebraic[56] = "J_cits in component J_cits (mole_per_second_per_l_mito_volume)"
legend_states[22] = "OAA_x in component dOAA_x_dt (molar)"
legend_states[23] = "SCOA_x in component dSCOA_x_dt (molar)"
legend_algebraic[53] = "K_eq_cits in component J_cits (dimensionless)"
legend_constants[315] = "K0_eq_cits in component J_cits (molar_squared)"
legend_constants[299] = "dGr0_cits in component J_cits (kilojoule_per_mole)"
legend_constants[149] = "K_mA in component J_cits (molar)"
legend_constants[150] = "K_mB in component J_cits (molar)"
legend_constants[151] = "K_ia in component J_cits (molar)"
legend_constants[152] = "K_iCIT in component J_cits (molar)"
legend_constants[153] = "K_iATP in component J_cits (molar)"
legend_constants[154] = "K_iADP in component J_cits (molar)"
legend_constants[155] = "K_iAMP in component J_cits (molar)"
legend_constants[156] = "K_iCOASH in component J_cits (molar)"
legend_constants[157] = "K_iSCOA in component J_cits (molar)"
legend_algebraic[54] = "alpha_i1 in component J_cits (dimensionless)"
legend_algebraic[55] = "alpha_i2 in component J_cits (dimensionless)"
legend_algebraic[59] = "J_acon in component J_acon (mole_per_second_per_l_mito_volume)"
legend_states[24] = "ICIT_x in component dICIT_x_dt (molar)"
legend_algebraic[58] = "V_mr in component J_acon (mole_per_second_per_l_mito_volume)"
legend_algebraic[57] = "K_eq_acon in component J_acon (dimensionless)"
legend_constants[316] = "K0_eq_acon in component J_acon (dimensionless)"
legend_constants[300] = "dGr0_acon in component J_acon (kilojoule_per_mole)"
legend_constants[158] = "K_mA in component J_acon (molar)"
legend_constants[159] = "K_mP in component J_acon (molar)"
legend_algebraic[115] = "J_isod in component J_isod (mole_per_second_per_l_mito_volume)"
legend_states[25] = "AKG_x in component dAKG_x_dt (molar)"
legend_algebraic[60] = "K_eq_isod in component J_isod (molar)"
legend_constants[317] = "K0_eq_isod in component J_isod (molar_cubed)"
legend_constants[301] = "dGr0_isod in component J_isod (kilojoule_per_mole)"
legend_constants[160] = "K_mA in component J_isod (molar)"
legend_constants[161] = "K_mB in component J_isod (molar)"
legend_constants[162] = "K_ib in component J_isod (molar)"
legend_constants[163] = "K_iq in component J_isod (molar)"
legend_algebraic[61] = "alpha_i in component J_isod (dimensionless)"
legend_constants[164] = "n_H in component J_isod (dimensionless)"
legend_constants[165] = "K_iATP in component J_isod (molar)"
legend_constants[166] = "K_aADP in component J_isod (molar)"
legend_constants[167] = "minCon in component J_isod (molar)"
legend_algebraic[116] = "J_akgd in component J_akgd (mole_per_second_per_l_mito_volume)"
legend_algebraic[62] = "K_eq_akgd in component J_akgd (dimensionless)"
legend_constants[318] = "K0_eq_akgd in component J_akgd (molar)"
legend_constants[302] = "dGr0_akgd in component J_akgd (kilojoule_per_mole)"
legend_constants[168] = "K_mA in component J_akgd (molar)"
legend_constants[169] = "K_mB in component J_akgd (molar)"
legend_constants[170] = "K_mC in component J_akgd (molar)"
legend_constants[171] = "K_iq in component J_akgd (molar)"
legend_constants[172] = "K_ir1 in component J_akgd (molar)"
legend_constants[173] = "K_ir2 in component J_akgd (molar)"
legend_constants[174] = "K_iATP in component J_akgd (molar)"
legend_constants[175] = "K_aADP in component J_akgd (molar)"
legend_algebraic[63] = "alpha_i in component J_akgd (dimensionless)"
legend_constants[176] = "minCon in component J_akgd (molar)"
legend_algebraic[66] = "J_scoas in component J_scoas (mole_per_second_per_l_mito_volume)"
legend_states[26] = "SUC_x in component dSUC_x_dt (molar)"
legend_algebraic[64] = "K_eq_scoas in component J_scoas (dimensionless)"
legend_constants[319] = "K0_eq_scoas in component J_scoas (molar)"
legend_constants[303] = "dGr0_scoas in component J_scoas (kilojoule_per_mole)"
legend_algebraic[65] = "V_mr in component J_scoas (mole_per_second_per_l_mito_volume)"
legend_constants[177] = "K_ia in component J_scoas (molar)"
legend_constants[178] = "K_ib in component J_scoas (molar)"
legend_constants[179] = "K_ic in component J_scoas (molar)"
legend_constants[180] = "K_ip in component J_scoas (molar)"
legend_constants[181] = "K_iq in component J_scoas (molar)"
legend_constants[182] = "K_ir in component J_scoas (molar)"
legend_constants[183] = "K_mA in component J_scoas (molar)"
legend_constants[184] = "K_mB in component J_scoas (molar)"
legend_constants[185] = "K_mC in component J_scoas (molar)"
legend_constants[186] = "K_mP in component J_scoas (molar)"
legend_constants[187] = "K_mQ in component J_scoas (molar)"
legend_constants[188] = "K_mR in component J_scoas (molar)"
legend_algebraic[119] = "J_sdh in component J_sdh (mole_per_second_per_l_mito_volume)"
legend_algebraic[118] = "COQ_x in component dCOQ_x_dt (molar)"
legend_states[27] = "QH2_x in component dQH2_x_dt (molar)"
legend_states[28] = "FUM_x in component dFUM_x_dt (molar)"
legend_algebraic[67] = "K_eq_sdh in component J_sdh (dimensionless)"
legend_constants[334] = "K0_eq_sdh in component J_sdh (dimensionless)"
legend_constants[327] = "dGr0_sdh in component J_sdh (kilojoule_per_mole)"
legend_algebraic[68] = "V_mr in component J_sdh (mole_per_second_per_l_mito_volume)"
legend_constants[189] = "K_ia in component J_sdh (molar)"
legend_constants[190] = "K_iq in component J_sdh (molar)"
legend_constants[191] = "K_mSUC in component J_sdh (molar)"
legend_constants[192] = "K_mCOQ in component J_sdh (molar)"
legend_constants[193] = "K_mQH2 in component J_sdh (molar)"
legend_constants[194] = "K_mFUM in component J_sdh (molar)"
legend_algebraic[69] = "alpha_i in component J_sdh (dimensionless)"
legend_constants[195] = "K_iOAA in component J_sdh (molar)"
legend_constants[196] = "K_aSUC in component J_sdh (molar)"
legend_constants[197] = "K_aFUM in component J_sdh (molar)"
legend_algebraic[73] = "J_fum in component J_fum (mole_per_second_per_l_mito_volume)"
legend_states[29] = "MAL_x in component dMAL_x_dt (molar)"
legend_algebraic[70] = "K_eq_fum in component J_fum (dimensionless)"
legend_constants[342] = "K0_eq_fum in component J_fum (dimensionless)"
legend_constants[338] = "dGr0_fum in component J_fum (kilojoule_per_mole)"
legend_algebraic[71] = "V_mr in component J_fum (mole_per_second_per_l_mito_volume)"
legend_constants[198] = "K_mFUM in component J_fum (molar)"
legend_constants[199] = "K_mMAL in component J_fum (molar)"
legend_algebraic[72] = "alpha_i in component J_fum (dimensionless)"
legend_constants[200] = "K_iCIT in component J_fum (molar)"
legend_constants[201] = "K_iATP in component J_fum (molar)"
legend_constants[202] = "K_iADP in component J_fum (molar)"
legend_constants[203] = "K_iGTP in component J_fum (molar)"
legend_constants[204] = "K_iGDP in component J_fum (molar)"
legend_algebraic[117] = "J_mdh in component J_mdh (mole_per_second_per_l_mito_volume)"
legend_algebraic[74] = "K_eq_mdh in component J_mdh (dimensionless)"
legend_constants[321] = "K0_eq_mdh in component J_mdh (molar)"
legend_constants[305] = "dGr0_mdh in component J_mdh (kilojoule_per_mole)"
legend_algebraic[75] = "V_mr in component J_mdh (mole_per_second_per_l_mito_volume)"
legend_constants[205] = "K_ia in component J_mdh (molar)"
legend_constants[206] = "K_ib in component J_mdh (molar)"
legend_constants[207] = "K_ip in component J_mdh (molar)"
legend_constants[208] = "K_iq in component J_mdh (molar)"
legend_constants[209] = "K_mNAD in component J_mdh (molar)"
legend_constants[210] = "K_mMAL in component J_mdh (molar)"
legend_constants[211] = "K_mOAA in component J_mdh (molar)"
legend_constants[212] = "K_mNADH in component J_mdh (molar)"
legend_algebraic[76] = "alpha_i in component J_mdh (dimensionless)"
legend_constants[213] = "K_iATP in component J_mdh (molar)"
legend_constants[214] = "K_iADP in component J_mdh (molar)"
legend_constants[215] = "K_iAMP in component J_mdh (molar)"
legend_algebraic[80] = "J_ndk_f in component J_ndk (mole_per_second_per_l_mito_volume)"
legend_algebraic[81] = "J_ndk in component J_ndk (mole_per_second_per_l_mito_volume)"
legend_algebraic[77] = "K_eq_ndk in component J_ndk (dimensionless)"
legend_constants[312] = "K0_eq_ndk in component J_ndk (dimensionless)"
legend_constants[297] = "dGr0_ndk in component J_ndk (kilojoule_per_mole)"
legend_algebraic[78] = "V_mr in component J_ndk (mole_per_second_per_l_mito_volume)"
legend_constants[216] = "K_ia in component J_ndk (molar)"
legend_constants[217] = "K_ib in component J_ndk (molar)"
legend_constants[218] = "K_ip in component J_ndk (molar)"
legend_constants[219] = "K_iq in component J_ndk (molar)"
legend_constants[220] = "K_mA in component J_ndk (molar)"
legend_constants[221] = "K_mB in component J_ndk (molar)"
legend_constants[222] = "K_mP in component J_ndk (molar)"
legend_constants[223] = "K_mQ in component J_ndk (molar)"
legend_algebraic[79] = "alpha_i in component J_ndk (dimensionless)"
legend_constants[224] = "K_iAMP in component J_ndk (molar)"
legend_constants[225] = "minCon in component J_ndk (molar)"
legend_algebraic[85] = "J_got_f in component J_got (mole_per_second_per_l_mito_volume)"
legend_algebraic[86] = "J_got in component J_got (mole_per_second_per_l_mito_volume)"
legend_states[30] = "ASP_x in component dASP_x_dt (molar)"
legend_states[31] = "GLU_x in component dGLU_x_dt (molar)"
legend_algebraic[82] = "K_eq_got in component J_got (dimensionless)"
legend_constants[328] = "K0_eq_got in component J_got (dimensionless)"
legend_constants[322] = "dGr0_got in component J_got (kilojoule_per_mole)"
legend_algebraic[83] = "V_mr in component J_got (mole_per_second_per_l_mito_volume)"
legend_constants[226] = "K_ia in component J_got (molar)"
legend_constants[227] = "K_ib in component J_got (molar)"
legend_constants[228] = "K_ip in component J_got (molar)"
legend_constants[229] = "K_iq in component J_got (molar)"
legend_constants[230] = "K_mA in component J_got (molar)"
legend_constants[231] = "K_mB in component J_got (molar)"
legend_constants[232] = "K_mP in component J_got (molar)"
legend_constants[233] = "K_mQ in component J_got (molar)"
legend_algebraic[84] = "alpha_i in component J_got (dimensionless)"
legend_constants[234] = "K_iAKG in component J_got (molar)"
legend_constants[235] = "minCon in component J_got (molar)"
legend_algebraic[155] = "J_PYRH in component J_PYRH (mole_per_second_per_l_mito_volume)"
legend_states[32] = "PYR_i in component dPYR_i_dt (molar)"
legend_algebraic[156] = "J_GLUH in component J_GLUH (mole_per_second_per_l_mito_volume)"
legend_states[33] = "GLU_i in component dGLU_i_dt (molar)"
legend_algebraic[175] = "J_CITMAL in component J_CITMAL (mole_per_second_per_l_mito_volume)"
legend_states[34] = "MAL_i in component dMAL_i_dt (molar)"
legend_algebraic[87] = "J_AKGMAL in component J_AKGMAL (mole_per_second_per_l_mito_volume)"
legend_constants[236] = "K_mAKGi in component J_AKGMAL (molar)"
legend_constants[237] = "K_mAKGx in component J_AKGMAL (molar)"
legend_constants[238] = "K_mMALi in component J_AKGMAL (molar)"
legend_constants[239] = "K_mMALx in component J_AKGMAL (molar)"
legend_states[35] = "AKG_i in component dAKG_i_dt (molar)"
legend_algebraic[176] = "J_MALPI in component J_MALPI (mole_per_second_per_l_mito_volume)"
legend_algebraic[158] = "J_ASPGLU in component J_ASPGLU (mole_per_second_per_l_mito_volume)"
legend_algebraic[157] = "K_eqASPGLU in component J_ASPGLU (dimensionless)"
legend_constants[285] = "K_hASPGLU in component J_ASPGLU (molar)"
legend_constants[240] = "K_iASPi in component J_ASPGLU (molar)"
legend_constants[241] = "K_iASPx in component J_ASPGLU (molar)"
legend_constants[242] = "K_iGLUi in component J_ASPGLU (molar)"
legend_constants[243] = "K_iGLUx in component J_ASPGLU (molar)"
legend_states[36] = "ASP_i in component dASP_i_dt (molar)"
legend_constants[244] = "m in component J_ASPGLU (dimensionless)"
legend_constants[245] = "pKa_ASPGLU in component J_ASPGLU (dimensionless)"
legend_states[37] = "dPsi in component ddPsi_dt (millivolt)"
legend_algebraic[177] = "J_SUCMAL in component J_SUCMAL (mole_per_second_per_l_mito_volume)"
legend_states[38] = "SUC_i in component dSUC_i_dt (molar)"
legend_algebraic[88] = "J_PYRt in component J_PYRt (mole_per_second_per_l_mito_volume)"
legend_states[39] = "PYR_c in component dPYR_c_dt (molar)"
legend_algebraic[90] = "J_CITt in component J_CITt (mole_per_second_per_l_mito_volume)"
legend_states[40] = "CIT_c in component dCIT_c_dt (molar)"
legend_algebraic[91] = "J_MALt in component J_MALt (mole_per_second_per_l_mito_volume)"
legend_states[41] = "MAL_c in component dMAL_c_dt (molar)"
legend_algebraic[89] = "J_AKGt in component J_AKGt (mole_per_second_per_l_mito_volume)"
legend_states[42] = "AKG_c in component dAKG_c_dt (molar)"
legend_algebraic[92] = "J_SUCt in component J_SUCt (mole_per_second_per_l_mito_volume)"
legend_states[43] = "SUC_c in component dSUC_c_dt (molar)"
legend_algebraic[93] = "J_GLUt in component J_GLUt (mole_per_second_per_l_mito_volume)"
legend_states[44] = "GLU_c in component dGLU_c_dt (molar)"
legend_algebraic[94] = "J_ASPt in component J_ASPt (mole_per_second_per_l_mito_volume)"
legend_states[45] = "ASP_c in component dASP_c_dt (molar)"
legend_algebraic[15] = "J_FUMt in component J_FUMt (mole_per_second_per_l_mito_volume)"
legend_constants[246] = "zero in component J_FUMt (micrometer_per_second)"
legend_states[46] = "FUM_c in component dFUM_c_dt (molar)"
legend_states[47] = "FUM_i in component dFUM_i_dt (molar)"
legend_algebraic[16] = "J_ICITt in component J_ICITt (mole_per_second_per_l_mito_volume)"
legend_states[48] = "ICIT_c in component dICIT_c_dt (molar)"
legend_states[49] = "ICIT_i in component dICIT_i_dt (molar)"
legend_algebraic[95] = "J_ADPt in component J_ADPt (mole_per_second_per_l_mito_volume)"
legend_states[50] = "ADP_i in component dADP_i_dt (molar)"
legend_algebraic[96] = "J_ATPt in component J_ATPt (mole_per_second_per_l_mito_volume)"
legend_states[51] = "ATP_i in component dATP_i_dt (molar)"
legend_algebraic[97] = "J_AMPt in component J_AMPt (mole_per_second_per_l_mito_volume)"
legend_states[52] = "AMP_c in component dAMP_c_dt (molar)"
legend_states[53] = "AMP_i in component dAMP_i_dt (molar)"
legend_algebraic[98] = "J_PIt in component J_PIt (mole_per_second_per_l_mito_volume)"
legend_states[54] = "PI_c in component dPI_c_dt (molar)"
legend_algebraic[139] = "J_C1 in component J_C1 (mole_per_second_per_l_mito_volume)"
legend_constants[287] = "dGr0_C1 in component J_C1 (kilojoule_per_mole)"
legend_algebraic[100] = "K_eq_C1 in component J_C1 (per_molar)"
legend_algebraic[138] = "K_app_C1 in component J_C1 (dimensionless)"
legend_algebraic[122] = "J_C3 in component J_C3 (mole_per_second_per_l_mito_volume)"
legend_constants[288] = "dGr0_C3 in component J_C3 (kilojoule_per_mole)"
legend_algebraic[102] = "K_eq_C3 in component J_C3 (molar_squared)"
legend_algebraic[104] = "K_app_C3 in component J_C3 (dimensionless)"
legend_algebraic[106] = "QH2_x in component J_C3 (molar)"
legend_algebraic[120] = "COQ_x in component J_C3 (molar)"
legend_states[55] = "Cred_i in component dCred_i_dt (molar)"
legend_algebraic[121] = "Cox_i in component dCox_i_dt (molar)"
legend_constants[247] = "minCon in component J_C3 (molar)"
legend_algebraic[141] = "J_C4 in component J_C4 (mole_per_second_per_l_mito_volume)"
legend_constants[289] = "dGr0_C4 in component J_C4 (kilojoule_per_mole)"
legend_algebraic[108] = "K_eq_C4 in component J_C4 (dimensionless)"
legend_algebraic[140] = "K_app_C4 in component J_C4 (dimensionless)"
legend_algebraic[110] = "O2_x in component J_C4 (molar)"
legend_states[56] = "O2_x in component dO2_x_dt (molar)"
legend_constants[248] = "minCon in component J_C4 (molar)"
legend_algebraic[143] = "J_F1 in component J_F1 (mole_per_second_per_l_mito_volume)"
legend_constants[290] = "dGr0_F1 in component J_F1 (kilojoule_per_mole)"
legend_algebraic[112] = "K_eq_F1 in component J_F1 (per_molar_squared)"
legend_algebraic[142] = "K_app_F1 in component J_F1 (per_molar)"
legend_algebraic[178] = "J_ANT in component J_ANT (mole_per_second_per_l_mito_volume)"
legend_constants[249] = "minCond in component J_ANT (molar)"
legend_algebraic[179] = "J_PIHt in component J_PIHt (mole_per_second_per_l_mito_volume)"
legend_algebraic[144] = "J_Hle in component J_Hle (mole_per_second_per_l_mito_volume)"
legend_constants[250] = "minCond in component J_Hle (millivolt)"
legend_algebraic[181] = "J_KH in component J_KH (mole_per_second_per_l_mito_volume)"
legend_algebraic[103] = "J_HK in component J_HK (mole_per_second_per_l_cyto_volume)"
legend_constants[308] = "dGr0_HK in component J_HK (kilojoule_per_mole)"
legend_states[57] = "GLC_c in component dGLC_c_dt (molar)"
legend_states[58] = "G6P_c in component dG6P_c_dt (molar)"
legend_constants[251] = "K_ia in component J_HK (molar)"
legend_constants[252] = "K_ib in component J_HK (molar)"
legend_constants[253] = "K_ip in component J_HK (molar)"
legend_constants[254] = "K_iq in component J_HK (molar)"
legend_constants[255] = "K_mA in component J_HK (molar)"
legend_constants[256] = "K_mB in component J_HK (molar)"
legend_constants[257] = "K_mP in component J_HK (molar)"
legend_constants[258] = "K_mQ in component J_HK (molar)"
legend_constants[259] = "Ki_G6P in component J_HK (molar)"
legend_constants[324] = "K_eq_HK in component J_HK (molar)"
legend_algebraic[99] = "Kapp_HK in component J_HK (dimensionless)"
legend_algebraic[101] = "Kapp_HK_m in component J_HK (dimensionless)"
legend_algebraic[107] = "J_CKc in component J_CKc (mole_per_second_per_l_cyto_volume)"
legend_constants[260] = "X_CK in component J_CKc (mole_per_second_per_l_cyto_volume_per_molar_per_molar)"
legend_constants[309] = "K_CK in component J_CKc (per_molar)"
legend_states[59] = "PCr_c in component dPCr_c_dt (molar)"
legend_algebraic[105] = "Cr_c in component J_CKc (molar)"
legend_algebraic[182] = "J_AKi in component J_AK (mole_per_second_per_l_mito_volume)"
legend_algebraic[111] = "J_AKc in component J_AK (mole_per_second_per_l_mito_volume)"
legend_constants[261] = "X_AK in component J_AK (mole_per_second_per_l_mito_volume_per_molar_per_molar)"
legend_constants[380] = "dGr0_AK in component J_AK (kilojoule_per_mole)"
legend_constants[381] = "K_eq_AK in component J_AK (dimensionless)"
legend_algebraic[180] = "Kapp_AKi in component J_AK (dimensionless)"
legend_algebraic[109] = "Kapp_AKc in component J_AK (dimensionless)"
legend_algebraic[146] = "dNAD_x_dt in component dNAD_x_dt (molar_per_second)"
legend_algebraic[145] = "dNADH_x_dt in component dNADH_x_dt (molar_per_second)"
legend_algebraic[148] = "dCOQ_x_dt in component dCOQ_x_dt (molar_per_second)"
legend_algebraic[147] = "dQH2_x_dt in component dQH2_x_dt (molar_per_second)"
legend_algebraic[183] = "dATP_x_dt in component dATP_x_dt (molar_per_second)"
legend_algebraic[184] = "dADP_x_dt in component dADP_x_dt (molar_per_second)"
legend_constants[262] = "dAMP_x_dt in component dAMP_x_dt (molar_per_second)"
legend_algebraic[123] = "dGTP_x_dt in component dGTP_x_dt (molar_per_second)"
legend_algebraic[124] = "dGDP_x_dt in component dGDP_x_dt (molar_per_second)"
legend_algebraic[185] = "dPI_x_dt in component dPI_x_dt (molar_per_second)"
legend_algebraic[159] = "dPYR_x_dt in component dPYR_x_dt (molar_per_second)"
legend_algebraic[125] = "dACCOA_x_dt in component dACCOA_x_dt (molar_per_second)"
legend_algebraic[186] = "dCIT_x_dt in component dCIT_x_dt (molar_per_second)"
legend_algebraic[126] = "dICIT_x_dt in component dICIT_x_dt (molar_per_second)"
legend_algebraic[127] = "dAKG_x_dt in component dAKG_x_dt (molar_per_second)"
legend_algebraic[128] = "dSCOA_x_dt in component dSCOA_x_dt (molar_per_second)"
legend_algebraic[129] = "dCOASH_x_dt in component dCOASH_x_dt (molar_per_second)"
legend_algebraic[187] = "dSUC_x_dt in component dSUC_x_dt (molar_per_second)"
legend_algebraic[130] = "dFUM_x_dt in component dFUM_x_dt (molar_per_second)"
legend_algebraic[188] = "dMAL_x_dt in component dMAL_x_dt (molar_per_second)"
legend_algebraic[131] = "dOAA_x_dt in component dOAA_x_dt (molar_per_second)"
legend_algebraic[160] = "dGLU_x_dt in component dGLU_x_dt (molar_per_second)"
legend_algebraic[161] = "dASP_x_dt in component dASP_x_dt (molar_per_second)"
legend_constants[263] = "dO2_x_dt in component dO2_x_dt (molar_per_second)"
legend_constants[264] = "dCO2_tot_x_dt in component dCO2_tot_x_dt (molar_per_second)"
legend_constants[331] = "Rm_cyto in component dPYR_c_dt (mito_per_cyto)"
legend_constants[335] = "Rm_cyto in component dCIT_c_dt (mito_per_cyto)"
legend_constants[330] = "Rm_cyto in component dICIT_c_dt (mito_per_cyto)"
legend_constants[339] = "Rm_cyto in component dAKG_c_dt (mito_per_cyto)"
legend_constants[343] = "Rm_cyto in component dSUC_c_dt (mito_per_cyto)"
legend_constants[329] = "Rm_cyto in component dFUM_c_dt (mito_per_cyto)"
legend_constants[346] = "Rm_cyto in component dMAL_c_dt (mito_per_cyto)"
legend_constants[348] = "Rm_cyto in component dGLU_c_dt (mito_per_cyto)"
legend_constants[350] = "Rm_cyto in component dASP_c_dt (mito_per_cyto)"
legend_constants[352] = "Rm_cyto in component dPI_c_dt (mito_per_cyto)"
legend_constants[265] = "conversion_factor in component dPI_c_dt (mito_per_cyto)"
legend_constants[354] = "Rm_cyto in component dATP_c_dt (mito_per_cyto)"
legend_constants[266] = "conversion_factor in component dATP_c_dt (mito_per_cyto)"
legend_constants[356] = "Rm_cyto in component dADP_c_dt (mito_per_cyto)"
legend_constants[267] = "conversion_factor in component dADP_c_dt (mito_per_cyto)"
legend_constants[358] = "Rm_cyto in component dAMP_c_dt (mito_per_cyto)"
legend_constants[268] = "conversion_factor in component dAMP_c_dt (mito_per_cyto)"
legend_constants[310] = "FAD_x in component dFAD_x_dt (molar)"
legend_constants[269] = "dFAD_x_dt in component dFAD_x_dt (molar_per_second)"
legend_constants[270] = "dFADH2_x_dt in component dFAD_x_dt (molar_per_second)"
legend_constants[291] = "FADH2_x in component dFAD_x_dt (molar)"
legend_constants[292] = "dH_x_dt in component dH_x_dt (molar_per_second)"
legend_algebraic[206] = "D in component BINDING_IONS (dimensionless)"
legend_algebraic[200] = "phi_H in component BINDING_IONS (molar_per_second)"
legend_algebraic[201] = "phi_Mg in component BINDING_IONS (molar_per_second)"
legend_algebraic[202] = "phi_K in component BINDING_IONS (molar_per_second)"
legend_algebraic[189] = "dHB_dMg in component BINDING_IONS (dimensionless)"
legend_algebraic[190] = "dHB_dK in component BINDING_IONS (dimensionless)"
legend_algebraic[193] = "dMgB_dK in component BINDING_IONS (dimensionless)"
legend_algebraic[196] = "dKB_dMg in component BINDING_IONS (dimensionless)"
legend_algebraic[205] = "alpha_K in component BINDING_IONS (dimensionless)"
legend_algebraic[204] = "alpha_Mg in component BINDING_IONS (dimensionless)"
legend_constants[293] = "dMg_x_dt in component dMg_x_dt (molar_per_second)"
legend_algebraic[203] = "alpha_H in component BINDING_IONS (dimensionless)"
legend_algebraic[195] = "dKB_dH in component BINDING_IONS (dimensionless)"
legend_algebraic[192] = "dMgB_dH in component BINDING_IONS (dimensionless)"
legend_constants[294] = "dK_x_dt in component dK_x_dt (molar_per_second)"
legend_algebraic[191] = "dHB_dH in component BINDING_IONS (dimensionless)"
legend_algebraic[194] = "dMgB_dMg in component BINDING_IONS (dimensionless)"
legend_algebraic[197] = "dKB_dK in component BINDING_IONS (dimensionless)"
legend_algebraic[198] = "J_H_t in component BINDING_IONS (molar_per_second)"
legend_constants[295] = "J_Mg_t in component BINDING_IONS (molar_per_second)"
legend_algebraic[199] = "J_K_t in component BINDING_IONS (molar_per_second)"
legend_constants[271] = "B_x in component BINDING_IONS (molar)"
legend_constants[272] = "K_BX in component BINDING_IONS (molar)"
legend_rates[37] = "d/dt dPsi in component ddPsi_dt (millivolt)"
legend_rates[0] = "d/dt ATP_x in component dATP_x_dt (molar)"
legend_rates[8] = "d/dt ADP_x in component dADP_x_dt (molar)"
legend_rates[10] = "d/dt AMP_x in component dAMP_x_dt (molar)"
legend_rates[11] = "d/dt GTP_x in component dGTP_x_dt (molar)"
legend_rates[12] = "d/dt GDP_x in component dGDP_x_dt (molar)"
legend_rates[13] = "d/dt PI_x in component dPI_x_dt (molar)"
legend_rates[21] = "d/dt NADH_x in component dNADH_x_dt (molar)"
legend_rates[27] = "d/dt QH2_x in component dQH2_x_dt (molar)"
legend_rates[17] = "d/dt PYR_x in component dPYR_x_dt (molar)"
legend_rates[20] = "d/dt ACCOA_x in component dACCOA_x_dt (molar)"
legend_rates[15] = "d/dt CIT_x in component dCIT_x_dt (molar)"
legend_rates[24] = "d/dt ICIT_x in component dICIT_x_dt (molar)"
legend_rates[25] = "d/dt AKG_x in component dAKG_x_dt (molar)"
legend_rates[23] = "d/dt SCOA_x in component dSCOA_x_dt (molar)"
legend_rates[18] = "d/dt COASH_x in component dCOASH_x_dt (molar)"
legend_rates[26] = "d/dt SUC_x in component dSUC_x_dt (molar)"
legend_rates[28] = "d/dt FUM_x in component dFUM_x_dt (molar)"
legend_rates[29] = "d/dt MAL_x in component dMAL_x_dt (molar)"
legend_rates[22] = "d/dt OAA_x in component dOAA_x_dt (molar)"
legend_rates[31] = "d/dt GLU_x in component dGLU_x_dt (molar)"
legend_rates[30] = "d/dt ASP_x in component dASP_x_dt (molar)"
legend_rates[56] = "d/dt O2_x in component dO2_x_dt (molar)"
legend_rates[19] = "d/dt CO2_tot_x in component dCO2_tot_x_dt (molar)"
legend_rates[55] = "d/dt Cred_i in component dCred_i_dt (molar)"
legend_rates[51] = "d/dt ATP_i in component dATP_i_dt (molar)"
legend_rates[50] = "d/dt ADP_i in component dADP_i_dt (molar)"
legend_rates[53] = "d/dt AMP_i in component dAMP_i_dt (molar)"
legend_rates[14] = "d/dt PI_i in component dPI_i_dt (molar)"
legend_rates[32] = "d/dt PYR_i in component dPYR_i_dt (molar)"
legend_rates[16] = "d/dt CIT_i in component dCIT_i_dt (molar)"
legend_rates[49] = "d/dt ICIT_i in component dICIT_i_dt (molar)"
legend_rates[35] = "d/dt AKG_i in component dAKG_i_dt (molar)"
legend_rates[38] = "d/dt SUC_i in component dSUC_i_dt (molar)"
legend_rates[47] = "d/dt FUM_i in component dFUM_i_dt (molar)"
legend_rates[34] = "d/dt MAL_i in component dMAL_i_dt (molar)"
legend_rates[33] = "d/dt GLU_i in component dGLU_i_dt (molar)"
legend_rates[36] = "d/dt ASP_i in component dASP_i_dt (molar)"
legend_rates[39] = "d/dt PYR_c in component dPYR_c_dt (molar)"
legend_rates[40] = "d/dt CIT_c in component dCIT_c_dt (molar)"
legend_rates[48] = "d/dt ICIT_c in component dICIT_c_dt (molar)"
legend_rates[42] = "d/dt AKG_c in component dAKG_c_dt (molar)"
legend_rates[43] = "d/dt SUC_c in component dSUC_c_dt (molar)"
legend_rates[46] = "d/dt FUM_c in component dFUM_c_dt (molar)"
legend_rates[41] = "d/dt MAL_c in component dMAL_c_dt (molar)"
legend_rates[44] = "d/dt GLU_c in component dGLU_c_dt (molar)"
legend_rates[45] = "d/dt ASP_c in component dASP_c_dt (molar)"
legend_rates[5] = "d/dt H_c in component dH_c_dt (molar)"
legend_rates[6] = "d/dt Mg_c in component dMg_c_dt (molar)"
legend_rates[7] = "d/dt K_c in component dK_c_dt (molar)"
legend_rates[54] = "d/dt PI_c in component dPI_c_dt (molar)"
legend_rates[4] = "d/dt ATP_c in component dATP_c_dt (molar)"
legend_rates[9] = "d/dt ADP_c in component dADP_c_dt (molar)"
legend_rates[57] = "d/dt GLC_c in component dGLC_c_dt (molar)"
legend_rates[58] = "d/dt G6P_c in component dG6P_c_dt (molar)"
legend_rates[59] = "d/dt PCr_c in component dPCr_c_dt (molar)"
legend_rates[52] = "d/dt AMP_c in component dAMP_c_dt (molar)"
legend_rates[1] = "d/dt H_x in component dH_x_dt (molar)"
legend_rates[2] = "d/dt Mg_x in component dMg_x_dt (molar)"
legend_rates[3] = "d/dt K_x in component dK_x_dt (molar)"
return (legend_states, legend_algebraic, legend_voi, legend_constants)
def initConsts():
constants = [0.0] * sizeConstants; states = [0.0] * sizeStates;
constants[0] = 1
constants[1] = 2
constants[2] = 0.096484
constants[3] = 0.651384
constants[4] = 0.072376
constants[5] = 5.99
constants[6] = 2.725e5
constants[7] = 3
constants[8] = 327
constants[9] = 85
constants[10] = 85
constants[11] = 3.5e-6
constants[12] = 0.6
constants[13] = 1.2e-4
constants[14] = 2.97e-3
constants[15] = 1.35e-3
constants[16] = 2.7e-3
constants[17] = 10e-3
constants[18] = 42.7e-3
constants[19] = 21.4e-3
constants[20] = 1e-4
constants[21] = 6.75e-6
constants[22] = 0.205093643103884
constants[23] = 9.82530750761898
constants[24] = 0.027665411180519
constants[25] = 0.492047553306503
constants[26] = 0.087630003425524
constants[27] = 6.04e-7
constants[28] = 0.448520452008943
constants[29] = 0.08578867843174
constants[30] = 0.00707728025683
constants[31] = 0.077060004692409
constants[32] = 0.025543515147654
constants[33] = 4.551391837332316
constants[34] = 2.996779273769102e8
constants[35] = 2.725114215177993e8
constants[36] = 80.432110322866379
constants[37] = 0.282091450119444
constants[38] = 96.995618247389871
constants[39] = 16.062441848715917
constants[40] = 5.865621634754698e-5
constants[41] = 3.236481869934019e4
constants[42] = 0.790813746051711
constants[43] = 3.248913595933763e-5
constants[44] = 0.003474542589709
constants[45] = 1.076088893810943e-4
constants[46] = 8.120393248498723e3
constants[47] = 0.007522124149343
constants[48] = 3.335568500871044e7
constants[49] = 0.00161251709168
constants[50] = 4.758023126099895e6
constants[51] = 3.02092849850504e2
constants[52] = 0
constants[53] = -235.74
constants[54] = 16.4
constants[55] = 39.31
constants[56] = 18.1
constants[57] = -23.3
constants[58] = 65.17
constants[59] = -2771
constants[60] = -1903.96
constants[61] = -1034.66
constants[62] = -2771
constants[63] = -1903.96
constants[64] = -27.41
constants[65] = -6.52
constants[66] = -1098.27
constants[67] = 0
constants[68] = -252.68
constants[69] = -67.6
constants[70] = 19.55
constants[71] = -0.72
constants[72] = -178.19
constants[73] = -794.74
constants[74] = -1165.59
constants[75] = -1158.94
constants[76] = -793.41
constants[77] = -507.55
constants[78] = -690.44
constants[79] = -603.32
constants[80] = -842.66
constants[81] = -692.26
constants[82] = -692.4
constants[83] = -530.71
constants[84] = -470.82
constants[85] = -907.21
constants[86] = -1758.87
states[0] = 1.025e-11
states[1] = 5.7919e-8
states[2] = 0.0046388
states[3] = 0.086878
states[4] = 1e-9
states[5] = 1e-7
states[6] = 0.005
states[7] = 0.15
states[8] = 0.01
states[9] = 1e-9
states[10] = 1e-6
states[11] = 2.7157e-10
states[12] = 0.005
states[13] = 2.9918e-4
states[14] = 1.1574e-4
constants[87] = 1e20
constants[88] = 1e20
constants[89] = 1e20
constants[90] = 1e20
states[15] = 6.815100000000001e-5
states[16] = 1e-9
constants[91] = 1e20
constants[92] = 1e20
constants[93] = 1e20
constants[94] = 1e20
constants[95] = 1e20
constants[96] = 1e20
constants[97] = 1e20
constants[98] = 1e20
constants[99] = 1e20
constants[100] = 1e20
constants[101] = 1e20
constants[102] = 1e20
constants[103] = 1e20
constants[104] = 1e20
constants[105] = 1e20
constants[106] = 1e20
constants[107] = 1e20
constants[108] = 1e20
constants[109] = 1e20
constants[110] = 1e20
constants[111] = 1e20
constants[112] = 1e20
constants[113] = 1e20
constants[114] = 1e20
constants[115] = 1e20
constants[116] = 1e20
constants[117] = 1e20
constants[118] = 1e20
constants[119] = 1e20
constants[120] = 1e20
constants[121] = 1e20
constants[122] = 1e20
constants[123] = 1e20
constants[124] = 1e20
constants[125] = 1e20
constants[126] = 1e20
constants[127] = 1e20
constants[128] = 1e20
constants[129] = 1e20
constants[130] = 1e20
constants[131] = 1e20
constants[132] = 1e20
constants[133] = 1e20
constants[134] = 1e20
constants[135] = 1e20
constants[136] = 1e20
constants[137] = 1e20
constants[138] = 1e20
constants[139] = 1e20
constants[140] = 1e20
constants[141] = 1e20
constants[142] = 1e20
states[17] = 1.2076e-11
states[18] = 0.0030107
states[19] = 0.0214
states[20] = 2.5507e-12
states[21] = 8.484200000000001e-21
constants[143] = 38.3e-6
constants[144] = 9.9e-6
constants[145] = 60.7e-6
constants[146] = 40.2e-6
constants[147] = 40e-6
constants[148] = 1e-32
states[22] = 8.9236e-4
states[23] = 6.9627e-10
constants[149] = 4e-6
constants[150] = 14e-6
constants[151] = 3.33e-6
constants[152] = 1600e-6
constants[153] = 900e-6
constants[154] = 1800e-6
constants[155] = 6000e-6
constants[156] = 67e-6
constants[157] = 140e-6
states[24] = 9.3558e-7
constants[158] = 1161e-6
constants[159] = 434e-6
states[25] = 4.7979e-11
constants[160] = 74e-6
constants[161] = 183e-6
constants[162] = 23.8e-6
constants[163] = 29e-6
constants[164] = 3
constants[165] = 91e-6
constants[166] = 50e-6
constants[167] = 1e-32
constants[168] = 80e-6
constants[169] = 55e-6
constants[170] = 21e-6
constants[171] = 6.9e-6
constants[172] = 6.010527470838273e-7
constants[173] = 1e3
constants[174] = 50e-6
constants[175] = 100e-6
constants[176] = 1e-32
states[26] = 2.1681e-7
constants[177] = 5.5e-6
constants[178] = 100e-6
constants[179] = 2000e-6
constants[180] = 20e-6
constants[181] = 3000e-6
constants[182] = 11.1e-6
constants[183] = 16e-6
constants[184] = 55e-6
constants[185] = 660e-6
constants[186] = 20e-6
constants[187] = 880e-6
constants[188] = 11.1e-6
states[27] = 0
states[28] = 7.8476e-9
constants[189] = 120e-6
constants[190] = 1275e-6
constants[191] = 467e-6
constants[192] = 480e-6
constants[193] = 2.45e-6
constants[194] = 1200e-6
constants[195] = 1.5e-6
constants[196] = 450e-6
constants[197] = 375e-6
states[29] = 3.2191e-8
constants[198] = 44.7e-6
constants[199] = 197.7e-6
constants[200] = 3500e-6
constants[201] = 40e-6
constants[202] = 400e-6
constants[203] = 80e-6
constants[204] = 330e-6
constants[205] = 279e-6
constants[206] = 360e-6
constants[207] = 5.5e-6
constants[208] = 3.18e-6
constants[209] = 90.55e-6
constants[210] = 250e-6
constants[211] = 6.128e-6
constants[212] = 2.58e-6
constants[213] = 183.2e-6
constants[214] = 394.4e-6
constants[215] = 420e-6
constants[216] = 170e-6
constants[217] = 143.6e-6
constants[218] = 146.6e-6
constants[219] = 156.5e-6
constants[220] = 111e-6
constants[221] = 100e-6
constants[222] = 260e-6
constants[223] = 278e-6
constants[224] = 650e-6
constants[225] = 1e-32
states[30] = 1.0217e-4
states[31] = 9.9311e-12
constants[226] = 3480e-6
constants[227] = 710e-6
constants[228] = 50e-6
constants[229] = 8400e-6
constants[230] = 3900e-6
constants[231] = 430e-6
constants[232] = 88e-6
constants[233] = 8900e-6
constants[234] = 16.6e-3
constants[235] = 1e-32
states[32] = 9.1804e-10
states[33] = 5.7539e-12
states[34] = 1e-9
constants[236] = 0.3e-3
constants[237] = 0.17e-3
constants[238] = 1.4e-3
constants[239] = 0.7e-3
states[35] = 9.999500000000001e-10
constants[240] = 28e-6
constants[241] = 2.8e-3
constants[242] = 180e-6
constants[243] = 1.6e-3
states[36] = 4.4404e-9
constants[244] = 1.8
constants[245] = 6.5
states[37] = -14.087
states[38] = 1e-9
states[39] = 1e-9
states[40] = 1e-9
states[41] = 1e-9
states[42] = 1e-9
states[43] = 1e-9
states[44] = 5.754e-12
states[45] = 4.4404e-9
constants[246] = 0
states[46] = 1e-9
states[47] = 1e-9
states[48] = 1e-9
states[49] = 1e-9
states[50] = 1e-9
states[51] = 1e-9
states[52] = 1e-9
states[53] = 1e-9
states[54] = 1.1574e-4
states[55] = 9.6996e-8
constants[247] = 1e-32
states[56] = 6.499999999999999e-5
constants[248] = 1e-32
constants[249] = 1e-9
constants[250] = 1e-9
states[57] = 1e-9
states[58] = 1e-9
constants[251] = 1e-3
constants[252] = 47e-6
constants[253] = 47e-6
constants[254] = 1e-3
constants[255] = 1e-3
constants[256] = 47e-6
constants[257] = 47e-6
constants[258] = 1e-3
constants[259] = 10e-6
constants[260] = 1e7
states[59] = 1e-9
constants[261] = 1e7
constants[262] = 0
constants[263] = 0
constants[264] = 0
constants[265] = 1
constants[266] = 1
constants[267] = 1
constants[268] = 1
constants[269] = 0
constants[270] = 0
constants[271] = 0.02
constants[272] = 1e-7
constants[273] = custom_piecewise([equal(constants[0] , 1.00000) & equal(constants[1] , 3.00000), 0.107108 , True, 0.00000])
constants[274] = custom_piecewise([equal(constants[0] , 1.00000), (8.31400*(28.0000+273.150))/1000.00 , True, (8.31400*(37.0000+273.150))/1000.00])
constants[275] = 1.00000
constants[276] = 1.00000
constants[277] = 1.00000
constants[278] = 1.00000
constants[279] = 1.00000
constants[280] = 1.00000
constants[281] = 1.00000
constants[282] = 1.00000
constants[283] = 1.00000
constants[284] = 1.00000
constants[285] = 1.00000*(power(10.0000, -constants[245]))
constants[286] = custom_piecewise([equal(constants[0] , 1.00000), 80.0000 , equal(constants[0] , 2.00000), 800.000 , equal(constants[0] , 3.00000) | equal(constants[0] , 4.00000), 0.807000*1.04400 , True, 80.0000])
constants[287] = ((constants[56]+constants[57])-constants[55])-constants[58]
constants[288] = ((constants[58]+2.00000*constants[64])-constants[57])-2.00000*constants[65]
constants[289] = ((2.00000*constants[65]+constants[53])-2.00000*constants[64])-0.500000*constants[54]
constants[290] = ((constants[59]+constants[53])-constants[60])-constants[66]
constants[291] = constants[20]/2.00000
constants[292] = 0.00000
constants[293] = 0.00000
constants[294] = 0.00000
constants[295] = 0.00000
constants[382] = 0.00000
constants[383] = 0.00000
constants[384] = 0.00000
constants[385] = 0.00000
constants[386] = 0.00000
constants[387] = 0.00000
constants[296] = 1.00000*(power(10.0000, -6.59000))
constants[297] = ((constants[63]+constants[59])-constants[62])-constants[60]
constants[298] = ((((constants[83]+constants[72]+constants[55])-constants[84])-constants[71])-constants[56])-constants[53]
constants[299] = (((constants[71]+constants[74])-constants[72])-constants[73])-constants[53]
constants[300] = constants[75]-constants[74]
constants[301] = (((constants[76]+constants[55]+constants[83])-constants[75])-constants[56])-constants[53]
constants[302] = ((((constants[83]+constants[77]+constants[55])-constants[76])-constants[71])-constants[56])-constants[53]
constants[303] = (((constants[71]+constants[78]+constants[62])-constants[63])-constants[77])-constants[66]
constants[304] = 1.00000
constants[305] = ((constants[73]+constants[55])-constants[56])-constants[80]
constants[306] = 1.00000
constants[307] = custom_piecewise([equal(constants[0] , 3.00000), 0.288200 , equal(constants[0] , 4.00000), 0.0560000 , True, 0.00000])
constants[308] = ((constants[86]+constants[60])-constants[85])-constants[59]
constants[309] = exp(50.7800/constants[274])/1.00000
constants[310] = constants[20]-constants[291]
constants[311] = 1.00000*(power(10.0000, -3.82000))
constants[312] = exp(-constants[297]/constants[274])
constants[313] = constants[296]
constants[314] = 1.00000*exp(-constants[298]/constants[274])
constants[315] = 1.00000*exp(-constants[299]/constants[274])
constants[316] = exp(-constants[300]/constants[274])
constants[317] = 1.00000*exp(-constants[301]/constants[274])
constants[318] = 1.00000*exp(-constants[302]/constants[274])
constants[319] = 1.00000*exp(-constants[303]/constants[274])
constants[320] = 1.00000
constants[321] = 1.00000*exp(-constants[305]/constants[274])
constants[322] = ((constants[73]+constants[82])-constants[81])-constants[76]
constants[323] = custom_piecewise([equal(constants[0] , 3.00000), 0.680100 , equal(constants[0] , 4.00000), 0.894000 , True, 1.00000])
constants[324] = 1.00000*exp(-constants[308]/constants[274])
constants[325] = 1.00000*(power(10.0000, -1.01300))
constants[326] = constants[311]
constants[327] = ((constants[57]+constants[79])-constants[78])-constants[58]
constants[328] = exp(-constants[322]/constants[274])
constants[329] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[330] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[331] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 3.00000) | equal(constants[0] , 4.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[332] = 1.00000*(power(10.0000, -6.42000))
constants[333] = constants[325]
constants[334] = exp(-constants[327]/constants[274])
constants[335] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[336] = 1.00000*(power(10.0000, -2.79000))
constants[337] = constants[332]
constants[338] = (constants[80]-constants[79])-constants[53]
constants[339] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[340] = 1.00000*(power(10.0000, -0.882000))
constants[341] = constants[336]
constants[342] = exp(-constants[338]/constants[274])
constants[343] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[344] = 1.00000*(power(10.0000, -6.22000))
constants[345] = constants[340]
constants[346] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[347] = 1.00000*(power(10.0000, -1.86000))
constants[348] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[349] = 1.00000*(power(10.0000, -0.621500))
constants[350] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), 1.00000 , True, constants[307]/constants[323]])
constants[351] = 1.00000*(power(10.0000, -6.71000))
constants[352] = constants[307]/constants[323]
constants[353] = 1.00000*(power(10.0000, -1.69000))
constants[354] = constants[307]/constants[323]
constants[355] = 1.00000*(power(10.0000, 0.00740000))
constants[356] = constants[307]/constants[323]
constants[357] = 1.00000*(power(10.0000, -8.13000))
constants[358] = constants[307]/constants[323]
constants[359] = 1.00000*(power(10.0000, -0.862900))
constants[360] = 1.00000*(power(10.0000, -5.63000))
constants[361] = 1.00000*(power(10.0000, -3.37000))
constants[362] = 1.00000*(power(10.0000, -0.339000))
constants[363] = 1.00000*(power(10.0000, -5.64000))
constants[364] = 1.00000*(power(10.0000, -2.46000))
constants[365] = 1.00000*(power(10.0000, -3.96000))
constants[366] = 1.00000*(power(10.0000, -5.13000))
constants[367] = 1.00000*(power(10.0000, -1.17000))
constants[368] = 1.00000*(power(10.0000, -0.352500))
constants[369] = 1.00000*(power(10.0000, -4.10000))
constants[370] = 1.00000*(power(10.0000, -4.75000))
constants[371] = 1.00000*(power(10.0000, -1.55000))
constants[372] = 1.00000*(power(10.0000, 0.107000))
constants[373] = 1.00000*(power(10.0000, -9.82000))
constants[374] = 1.00000*(power(10.0000, -1.02000))
constants[375] = 1.00000*(power(10.0000, -5.91000))
constants[376] = 1.00000*(power(10.0000, -4.06000))
constants[377] = 1.00000*(power(10.0000, -1.82000))
constants[378] = 1.00000*(power(10.0000, -3.65000))
constants[379] = 1.00000*(power(10.0000, -2.32000))
constants[380] = (constants[59]+constants[61])-2.00000*constants[60]
constants[381] = exp(-constants[380]/constants[274])
return (states, constants)
def computeRates(voi, states, constants):
rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic
rates[10] = constants[382]
rates[56] = constants[383]
rates[19] = constants[384]
rates[5] = constants[385]
rates[6] = constants[386]
rates[7] = constants[387]
algebraic[16] = constants[5]*constants[10]*(states[48]-states[49])
rates[49] = algebraic[16]/constants[4]
algebraic[15] = constants[5]*constants[246]*(states[46]-states[47])
rates[47] = algebraic[15]/constants[4]
rates[48] = (-constants[330]*algebraic[16])/constants[286]
rates[46] = (-constants[329]*algebraic[15])/constants[286]
algebraic[0] = 1.00000+states[1]/constants[296]+states[2]/constants[311]+states[3]/constants[325]
algebraic[29] = algebraic[0]
algebraic[21] = 1.00000+states[1]/constants[332]+states[2]/constants[336]+states[3]/constants[340]
algebraic[31] = algebraic[21]
algebraic[33] = 1.00000+states[1]/constants[351]+states[2]/constants[353]+states[3]/constants[355]
algebraic[35] = 1.00000+states[1]/constants[357]
algebraic[41] = 1.00000+states[1]/constants[365]
algebraic[42] = 1.00000+states[1]/constants[366]+states[2]/constants[367]+states[3]/constants[368]
algebraic[64] = (constants[319]*algebraic[35]*algebraic[42]*algebraic[29])/(states[1]*algebraic[31]*algebraic[41]*algebraic[33])
algebraic[65] = (constants[28]*constants[186]*constants[181]*constants[182])/(algebraic[64]*constants[177]*constants[178]*constants[185])
algebraic[66] = (constants[28]*algebraic[65]*(states[12]*states[23]*states[13]-(states[18]*states[26]*states[11])/algebraic[64]))/(algebraic[65]*constants[177]*constants[178]*constants[185]+algebraic[65]*constants[178]*constants[185]*states[12]+algebraic[65]*constants[177]*constants[184]*states[13]+algebraic[65]*constants[185]*states[12]*states[23]+algebraic[65]*constants[184]*states[12]*states[13]+algebraic[65]*constants[183]*states[23]*states[13]+algebraic[65]*states[12]*states[23]*states[13]+(constants[28]*constants[182]*constants[187]*states[18])/algebraic[64]+(constants[28]*constants[181]*constants[186]*states[11])/algebraic[64]+(constants[28]*constants[188]*states[18]*states[26])/algebraic[64]+(constants[28]*constants[187]*states[18]*states[11])/algebraic[64]+(constants[28]*constants[186]*states[26]*states[11])/algebraic[64]+(constants[28]*states[18]*states[26]*states[11])/algebraic[64]+(constants[28]*constants[187]*constants[182]*states[12]*states[18])/(constants[177]*algebraic[64])+(algebraic[65]*constants[177]*constants[184]*states[13]*states[11])/constants[182]+(constants[28]*constants[187]*constants[182]*states[12]*states[23]*states[18])/(constants[177]*constants[178]*algebraic[64])+(algebraic[65]*constants[183]*states[23]*states[13]*states[11])/constants[182]+(constants[28]*constants[188]*states[12]*states[18]*states[26])/(constants[177]*algebraic[64])+(algebraic[65]*constants[177]*constants[184]*states[13]*states[26]*states[11])/(constants[181]*constants[182])+(constants[28]*constants[182]*constants[187]*states[12]*states[23]*states[13]*states[18])/(constants[177]*constants[178]*constants[179]*algebraic[64])+(constants[28]*constants[180]*constants[188]*states[12]*states[23]*states[13]*states[26])/(constants[177]*constants[178]*constants[179]*algebraic[64])+(constants[28]*constants[188]*states[12]*states[23]*states[18]*states[26])/(constants[177]*constants[178]*algebraic[64])+(algebraic[65]*constants[183]*states[23]*states[13]*states[26]*states[11])/(constants[181]*constants[182])+(algebraic[65]*constants[183]*constants[179]*states[23]*states[18]*states[26]*states[11])/(constants[180]*constants[181]*constants[182])+(algebraic[65]*constants[177]*constants[184]*states[13]*states[18]*states[26]*states[11])/(constants[180]*constants[181]*constants[182])+(constants[28]*constants[188]*states[12]*states[23]*states[13]*states[18]*states[26])/(constants[177]*constants[178]*constants[179]*algebraic[64])+(algebraic[65]*constants[183]*states[23]*states[13]*states[18]*states[26]*states[11])/(constants[180]*constants[181]*constants[182]))
algebraic[77] = (constants[312]*algebraic[31]*algebraic[0])/(algebraic[29]*algebraic[21])
algebraic[78] = (constants[32]*constants[223]*constants[218])/(algebraic[77]*constants[216]*constants[221])
algebraic[26] = 1.00000+states[1]/constants[344]+states[2]/constants[347]+states[3]/constants[349]
algebraic[28] = (states[10]*(1.00000+states[1]/constants[344]))/algebraic[26]
algebraic[79] = 1.00000+algebraic[28]/constants[224]
algebraic[80] = custom_piecewise([greater(states[11] , constants[225]) & greater(states[8] , constants[225]), (((constants[32]*algebraic[78])/algebraic[79])*states[11]*states[8])/(algebraic[78]*constants[221]*states[11]+algebraic[78]*constants[220]*states[8]+(constants[32]*constants[223]*states[12])/algebraic[77]+(constants[32]*constants[222]*states[0])/algebraic[77]+algebraic[78]*states[11]*states[8]+(constants[32]*constants[223]*states[11]*states[12])/(algebraic[77]*constants[216])+(constants[32]*states[12]*states[0])/algebraic[77]+(algebraic[78]*constants[220]*states[8]*states[0])/constants[219]) , True, 0.00000])
algebraic[81] = custom_piecewise([greater(states[12] , constants[225]) & greater(states[0] , constants[225]), algebraic[80]-((((constants[32]*algebraic[78])/algebraic[79])*states[12]*states[0])/algebraic[77])/(algebraic[78]*constants[221]*states[11]+algebraic[78]*constants[220]*states[8]+(constants[32]*constants[223]*states[12])/algebraic[77]+(constants[32]*constants[222]*states[0])/algebraic[77]+algebraic[78]*states[11]*states[8]+(constants[32]*constants[223]*states[11]*states[12])/(algebraic[77]*constants[216])+(constants[32]*states[12]*states[0])/algebraic[77]+(algebraic[78]*constants[220]*states[8]*states[0])/constants[219]) , True, algebraic[80]])
rates[11] = (algebraic[66]-algebraic[81])/constants[3]
rates[12] = (-algebraic[66]+algebraic[81])/constants[3]
algebraic[87] = (constants[37]*(states[35]*states[29]-states[25]*states[34]))/(constants[236]*constants[239]*(2.00000+states[34]/constants[238]+states[29]/constants[239]+states[35]/constants[236]+states[25]/constants[237]+(states[34]*states[25])/(constants[238]*constants[237])+(states[29]*states[35])/(constants[239]*constants[236])))
algebraic[89] = constants[5]*constants[10]*(states[42]-states[35])
rates[35] = (-algebraic[87]+algebraic[89])/constants[4]
algebraic[88] = constants[5]*constants[10]*(states[39]-states[32])
rates[39] = (-constants[331]*algebraic[88])/constants[286]
rates[42] = (-constants[339]*algebraic[89])/constants[286]
algebraic[90] = constants[5]*constants[10]*(states[40]-states[16])
rates[40] = (-constants[335]*algebraic[90])/constants[286]
algebraic[91] = constants[5]*constants[10]*(states[41]-states[34])
rates[41] = (-constants[346]*algebraic[91])/constants[286]
algebraic[92] = constants[5]*constants[10]*(states[43]-states[38])
rates[43] = (-constants[343]*algebraic[92])/constants[286]
algebraic[93] = constants[5]*constants[10]*(states[44]-states[33])
rates[44] = (-constants[348]*algebraic[93])/constants[286]
algebraic[94] = constants[5]*constants[10]*(states[45]-states[36])
rates[45] = (-constants[350]*algebraic[94])/constants[286]
algebraic[98] = constants[5]*constants[8]*(states[54]-states[14])
rates[54] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), (-algebraic[98]*constants[265])/constants[286] , True, (constants[352]*algebraic[98]+constants[52])/constants[286]])
algebraic[17] = 1.00000+states[5]/constants[296]+states[6]/constants[311]+states[7]/constants[325]
algebraic[20] = ((states[4]*states[6])/constants[311])/algebraic[17]
algebraic[22] = 1.00000+states[5]/constants[332]+states[6]/constants[336]+states[7]/constants[340]
algebraic[25] = ((states[9]*states[6])/constants[336])/algebraic[22]
algebraic[47] = 1.00000+states[5]/constants[375]
algebraic[99] = (constants[324]*algebraic[47]*algebraic[22])/(states[5]*algebraic[17]*constants[275])
algebraic[101] = (algebraic[99]*constants[311]*algebraic[17])/(constants[336]*algebraic[22])
algebraic[103] = ((constants[273]/(constants[252]*constants[255]))*(algebraic[20]*states[57]-(states[58]*algebraic[25])/algebraic[101]))/(1.00000+algebraic[20]/constants[251]+states[57]/constants[252]+(algebraic[20]*states[57])/(constants[252]*constants[255])+states[58]/constants[253]+algebraic[25]/constants[254]+(states[58]*algebraic[25])/(constants[254]*constants[257])+(states[58]*states[57])/(constants[259]*constants[252]))
rates[57] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), (-(1.00000+constants[286]/1.00000)*algebraic[103])/constants[286] , True, 0.00000])
rates[58] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), ((1.00000+constants[286]/1.00000)*algebraic[103])/constants[286] , True, 0.00000])
algebraic[105] = constants[18]-states[59]
algebraic[107] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[260]*(((constants[309]*states[9])/algebraic[22])*states[59]*states[5]-(states[4]/algebraic[17])*algebraic[105])])
rates[59] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), 0.00000 , True, -algebraic[107]/constants[286]])
algebraic[96] = constants[5]*constants[9]*(states[4]-states[51])
algebraic[27] = 1.00000+states[5]/constants[344]+states[6]/constants[347]+states[7]/constants[349]
algebraic[109] = (constants[381]*algebraic[17]*algebraic[27])/(algebraic[22]*algebraic[22])
algebraic[111] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[261]*(algebraic[109]*states[9]*states[9]-states[52]*states[4])])
rates[4] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), (-(1.00000+constants[286]/1.00000)*algebraic[103]-constants[266]*algebraic[96])/constants[286] , True, ((-constants[354]*algebraic[96]+algebraic[107]+constants[266]*algebraic[111])-constants[52])/constants[286]])
algebraic[95] = constants[5]*constants[9]*(states[9]-states[50])
rates[9] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), ((1.00000+constants[286]/1.00000)*algebraic[103]-constants[267]*algebraic[95])/constants[286] , True, (((-constants[356]*algebraic[95]-algebraic[107])-2.00000*constants[267]*algebraic[111])+constants[52])/constants[286]])
algebraic[97] = constants[5]*constants[9]*(states[52]-states[53])
rates[52] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), 0.00000 , True, (-constants[358]*algebraic[97]+constants[268]*algebraic[111])/constants[286]])
algebraic[113] = constants[14]-states[21]
algebraic[45] = 1.00000+states[1]/constants[373]
algebraic[46] = 1.00000+states[2]/constants[374]
algebraic[50] = (constants[314]*algebraic[45]*constants[278]*constants[277])/(states[1]*algebraic[46]*algebraic[35]*constants[276])
algebraic[51] = 1.00000+states[20]/constants[146]
algebraic[52] = 1.00000+states[21]/constants[147]
algebraic[114] = custom_piecewise([greater(states[17] , constants[148]) & greater(states[18] , constants[148]) & greater(algebraic[113] , constants[148]), (constants[22]*(states[17]*states[18]*algebraic[113]-(states[19]*states[20]*states[21])/algebraic[50]))/(constants[145]*algebraic[52]*states[17]*states[18]+constants[144]*algebraic[51]*states[17]*algebraic[113]+constants[143]*states[18]*algebraic[113]+states[17]*states[18]*algebraic[113]) , True, 0.00000])
algebraic[36] = 1.00000+states[2]/constants[359]
algebraic[37] = 1.00000+states[1]/constants[360]+states[2]/constants[361]+states[3]/constants[362]
algebraic[53] = (constants[315]*algebraic[37]*algebraic[35])/((power(states[1], 2.00000))*constants[278]*algebraic[36])
algebraic[38] = (states[15]*(1.00000+states[1]/constants[360]))/algebraic[37]
algebraic[54] = 1.00000+algebraic[38]/constants[152]
algebraic[18] = (states[0]*(1.00000+states[1]/constants[296]))/algebraic[0]
algebraic[23] = (states[8]*(1.00000+states[1]/constants[332]))/algebraic[21]
algebraic[55] = 1.00000+algebraic[18]/constants[153]+algebraic[23]/constants[154]+algebraic[28]/constants[155]+states[18]/constants[156]+states[23]/constants[157]
algebraic[56] = (constants[23]*(states[20]*states[22]-(states[18]*states[15])/algebraic[53]))/(constants[151]*constants[150]*algebraic[54]+constants[149]*algebraic[54]*states[20]+constants[150]*algebraic[55]*states[22]+states[22]*states[20])
rates[20] = (algebraic[114]-algebraic[56])/constants[3]
algebraic[40] = 1.00000+states[1]/constants[363]+states[2]/constants[364]
algebraic[57] = (constants[316]*algebraic[40])/algebraic[37]
algebraic[58] = (constants[24]*constants[159])/(constants[158]*algebraic[57])
algebraic[59] = (constants[24]*algebraic[58]*(states[15]-states[24]/algebraic[57]))/(constants[158]*algebraic[58]+algebraic[58]*states[15]+(constants[24]*states[24])/algebraic[57])
algebraic[60] = (constants[317]*algebraic[45]*constants[306]*constants[277])/((power(states[1], 2.00000))*constants[276]*algebraic[40])
algebraic[61] = 1.00000+(constants[166]/algebraic[23])*(1.00000+algebraic[18]/constants[165])
algebraic[115] = custom_piecewise([greater(algebraic[113] , constants[167]) & greater(states[24] , constants[167]), (constants[25]*(1.00000-(states[19]*states[25]*states[21])/(algebraic[60]*algebraic[113]*states[24])))/(1.00000+(power(constants[161]/states[24], constants[164]))*algebraic[61]+(constants[160]/algebraic[113])*(1.00000+(power(constants[162]/states[24], constants[164]))*algebraic[61]+(states[21]*algebraic[61])/constants[163])) , True, 0.00000])
rates[24] = (algebraic[59]-algebraic[115])/constants[3]
algebraic[62] = (constants[318]*algebraic[45]*algebraic[41]*constants[277])/(states[1]*algebraic[35]*constants[306]*constants[276])
algebraic[63] = 1.00000+(constants[175]/algebraic[23])*(1.00000+algebraic[18]/constants[174])
algebraic[116] = custom_piecewise([greater(states[25] , constants[176]) & greater(states[18] , constants[176]) & greater(algebraic[113] , constants[176]), ((constants[26]*(1.00000-(states[19]*states[23]*states[21])/(algebraic[62]*states[25]*states[18]*algebraic[113])))/(1.00000+(constants[168]/states[25])*algebraic[63]+(constants[169]/states[18])*(1.00000+states[23]/constants[171])+(constants[170]/algebraic[113])*(1.00000+states[21]/constants[172])))/(1.00000+states[21]/constants[173]) , True, 0.00000])
algebraic[48] = 1.00000+states[1]/constants[376]+states[2]/constants[377]
algebraic[49] = 1.00000+states[1]/constants[378]+states[2]/constants[379]
algebraic[82] = (constants[328]*algebraic[48]*algebraic[36])/(algebraic[49]*constants[306])
algebraic[83] = (constants[33]*constants[233]*constants[228])/(algebraic[82]*constants[226]*constants[231])
algebraic[84] = 1.00000+states[25]/constants[234]
algebraic[85] = custom_piecewise([greater(states[30] , constants[235]) & greater(states[25] , constants[235]), (constants[33]*algebraic[83]*states[30]*states[25])/(algebraic[83]*constants[231]*states[30]+algebraic[83]*constants[230]*algebraic[84]*states[25]+(constants[33]*constants[233]*algebraic[84]*states[22])/algebraic[82]+(constants[33]*constants[232]*states[31])/algebraic[82]+algebraic[83]*states[30]*states[25]+(constants[33]*constants[233]*states[30]*states[22])/(algebraic[82]*constants[226])+(constants[33]*states[22]*states[31])/algebraic[82]+(algebraic[83]*constants[230]*states[25]*states[31])/constants[229]) , True, 0.00000])
algebraic[86] = custom_piecewise([greater(states[22] , constants[235]) & greater(states[31] , constants[235]), algebraic[85]-((constants[33]*algebraic[83]*states[22]*states[31])/algebraic[82])/(algebraic[83]*constants[231]*states[30]+algebraic[83]*constants[230]*algebraic[84]*states[25]+(constants[33]*constants[233]*algebraic[84]*states[22])/algebraic[82]+(constants[33]*constants[232]*states[31])/algebraic[82]+algebraic[83]*states[30]*states[25]+(constants[33]*constants[233]*states[30]*states[22])/(algebraic[82]*constants[226])+(constants[33]*states[22]*states[31])/algebraic[82]+(algebraic[83]*constants[230]*states[25]*states[31])/constants[229]) , True, algebraic[85]])
rates[25] = (((algebraic[115]-algebraic[116])-algebraic[86])+algebraic[87])/constants[3]
rates[23] = (algebraic[116]-algebraic[66])/constants[3]
rates[18] = ((-algebraic[114]-algebraic[116])+algebraic[66]+algebraic[56])/constants[3]
algebraic[44] = 1.00000+states[1]/constants[370]+states[2]/constants[371]+states[3]/constants[372]
algebraic[74] = (constants[321]*algebraic[36]*constants[277])/(states[1]*algebraic[44]*constants[276])
algebraic[75] = (constants[31]*constants[208]*constants[211])/(algebraic[74]*constants[205]*constants[210])
algebraic[76] = 1.00000+algebraic[18]/constants[213]+algebraic[23]/constants[214]+algebraic[28]/constants[215]
algebraic[117] = (constants[31]*algebraic[75]*(algebraic[113]*states[29]-(states[22]*states[21])/algebraic[74]))/(algebraic[75]*constants[205]*constants[210]*algebraic[76]+algebraic[75]*constants[210]*algebraic[113]+algebraic[75]*constants[209]*algebraic[76]*states[29]+(constants[31]*constants[212]*algebraic[76]*states[22])/algebraic[74]+(constants[31]*constants[211]*states[21])/algebraic[74]+algebraic[75]*algebraic[113]*states[29]+(constants[31]*constants[212]*algebraic[113]*states[22])/(algebraic[74]*constants[205])+(constants[31]*states[22]*states[21])/algebraic[74]+(algebraic[75]*constants[209]*states[29]*states[21])/constants[208]+(algebraic[75]*algebraic[113]*states[29]*states[22])/constants[207]+(constants[31]*states[29]*states[22]*states[21])/(constants[206]*algebraic[74]))
rates[22] = (-algebraic[56]+algebraic[117]+algebraic[86])/constants[3]
algebraic[118] = constants[15]-states[27]
algebraic[43] = 1.00000+states[1]/constants[369]
algebraic[67] = (constants[334]*algebraic[43]*constants[304])/(algebraic[42]*constants[320])
algebraic[68] = (constants[29]*constants[193]*constants[190])/(algebraic[67]*constants[189]*constants[192])
algebraic[69] = (1.00000+states[22]/constants[195]+states[26]/constants[196]+states[28]/constants[197])/(1.00000+states[26]/constants[196]+states[28]/constants[197])
algebraic[119] = (constants[29]*algebraic[68]*(states[26]*algebraic[118]-(states[27]*states[28])/algebraic[67]))/(algebraic[68]*constants[189]*constants[192]*algebraic[69]+algebraic[68]*constants[192]*states[26]+algebraic[68]*constants[191]*algebraic[69]*algebraic[118]+(constants[29]*constants[194]*algebraic[69]*states[27])/algebraic[67]+(constants[29]*constants[193]*states[28])/algebraic[67]+algebraic[68]*states[26]*algebraic[118]+(constants[29]*constants[194]*states[26]*states[27])/(algebraic[67]*constants[189])+(algebraic[68]*constants[191]*algebraic[118]*states[28])/constants[190]+(constants[29]*states[27]*states[28])/algebraic[67])
algebraic[70] = (constants[342]*algebraic[44])/algebraic[43]
algebraic[71] = (constants[30]*constants[199])/(algebraic[70]*constants[198])
algebraic[30] = (states[11]*(1.00000+states[1]/constants[313]))/algebraic[29]
algebraic[32] = (states[12]*(1.00000+states[1]/constants[337]))/algebraic[31]
algebraic[72] = 1.00000+states[15]/constants[200]+algebraic[18]/constants[201]+algebraic[23]/constants[202]+algebraic[30]/constants[203]+algebraic[32]/constants[204]
algebraic[73] = (algebraic[71]*constants[30]*(states[28]-states[29]/algebraic[70]))/(constants[198]*algebraic[71]*algebraic[72]+algebraic[71]*states[28]+(constants[30]*states[29])/algebraic[70])
rates[28] = (algebraic[119]-algebraic[73])/constants[3]
algebraic[132] = states[5]
algebraic[100] = 1.00000*exp(-(constants[287]+4.00000*constants[2]*states[37])/constants[274])
algebraic[138] = (algebraic[100]*(power(states[1], 5.00000)))/(power(algebraic[132], 4.00000))
algebraic[139] = constants[41]*(algebraic[138]*states[21]*algebraic[118]-algebraic[113]*states[27])
rates[21] = ((algebraic[114]+algebraic[115]+algebraic[116]+algebraic[117])-algebraic[139])/constants[3]
algebraic[102] = 1.00000*exp(-(constants[288]+2.00000*constants[2]*states[37])/constants[274])
algebraic[104] = (algebraic[102]*(power(states[1], 2.00000)))/(power(states[1], 4.00000))
algebraic[106] = custom_piecewise([greater(states[27] , constants[247]), states[27] , True, constants[247]])
algebraic[120] = custom_piecewise([greater(algebraic[118] , constants[247]), algebraic[118] , True, constants[247]])
algebraic[121] = constants[16]-states[55]
algebraic[122] = ((constants[42]*(1.00000+states[13]/constants[43]))/(1.00000+states[13]/constants[44]))*((power(algebraic[104], 1.0/2))*algebraic[121]*(power(algebraic[106], 1.0/2))-states[55]*(power(algebraic[120], 1.0/2)))
rates[27] = ((algebraic[119]+algebraic[139])-algebraic[122])/constants[3]
algebraic[108] = exp(-(constants[289]+4.00000*constants[2]*states[37])/constants[274])
algebraic[140] = (algebraic[108]*(power(states[1], 4.00000)))/((power(algebraic[132], 2.00000))*(power(1.00000, 2.00000)))
algebraic[110] = custom_piecewise([greater(states[56] , constants[248]), states[56] , True, constants[248]])
algebraic[141] = constants[45]*((algebraic[110]/(algebraic[110]+constants[13]))*exp((constants[2]*states[37])/constants[274])*(states[55]/constants[16]))*((power(algebraic[140], 1.00000/2.00000))*states[55]*(power(algebraic[110]/1.00000, 1.00000/4.00000))-algebraic[121])
rates[55] = (2.00000*algebraic[122]-2.00000*algebraic[141])/constants[4]
algebraic[149] = states[6]
algebraic[152] = 1.00000+algebraic[149]/constants[374]
algebraic[155] = constants[34]*((states[32]/algebraic[152])*algebraic[132]-(states[17]/algebraic[46])*states[1])
rates[17] = (-algebraic[114]+algebraic[155])/constants[3]
rates[32] = (-algebraic[155]+algebraic[88])/constants[4]
algebraic[153] = 1.00000+algebraic[132]/constants[376]+algebraic[149]/constants[377]
algebraic[156] = constants[35]*((states[33]/algebraic[153])*algebraic[132]-(states[31]/algebraic[48])*states[1])
algebraic[154] = 1.00000+algebraic[132]/constants[378]+algebraic[149]/constants[379]
algebraic[157] = (exp((-constants[2]*states[37])/constants[274])*algebraic[49]*algebraic[153])/(algebraic[154]*algebraic[48])
algebraic[158] = (constants[40]*(algebraic[157]*states[36]*states[31]*states[1]-states[30]*states[33]*algebraic[132]))/(algebraic[157]*constants[240]*constants[243]*constants[285]*(2.00000*constants[244]+(constants[244]*states[36])/constants[240]+(states[36]*states[31]*states[1])/(constants[240]*constants[243]*constants[285])+(constants[244]*states[30]*algebraic[132])/(constants[241]*constants[285])+(states[30]*states[33]*algebraic[132])/(constants[241]*constants[242]*constants[285])+(constants[244]*states[30])/constants[241]+(constants[244]*states[36]*states[1])/(constants[240]*constants[285])+(constants[244]*states[1])/constants[285]+(constants[244]*states[33]*algebraic[132])/(constants[242]*constants[285])+(constants[244]*algebraic[132])/constants[285]+(constants[244]*states[31]*states[1])/(constants[243]*constants[285])))
rates[31] = ((algebraic[86]+algebraic[156])-algebraic[158])/constants[3]
rates[30] = (-algebraic[86]+algebraic[158])/constants[3]
rates[33] = (-algebraic[156]+algebraic[158]+algebraic[93])/constants[4]
rates[36] = (-algebraic[158]+algebraic[94])/constants[4]
algebraic[39] = ((states[15]*states[1])/constants[360])/algebraic[37]
algebraic[162] = states[7]
algebraic[170] = 1.00000+algebraic[132]/constants[360]+algebraic[149]/constants[361]+algebraic[162]/constants[362]
algebraic[172] = ((states[16]*algebraic[132])/constants[360])/algebraic[170]
algebraic[174] = 1.00000+algebraic[132]/constants[370]+algebraic[149]/constants[371]+algebraic[162]/constants[372]
algebraic[175] = constants[36]*((algebraic[172]*states[29])/algebraic[44]-(algebraic[39]*states[34])/algebraic[174])
rates[15] = (-algebraic[59]+algebraic[175]+algebraic[56])/constants[3]
rates[16] = (-algebraic[175]+algebraic[90])/constants[4]
algebraic[173] = 1.00000+algebraic[132]/constants[366]+algebraic[149]/constants[367]+algebraic[162]/constants[368]
algebraic[177] = constants[38]*(((states[38]/algebraic[173])*states[29])/algebraic[44]-((states[26]/algebraic[42])*states[34])/algebraic[174])
rates[26] = ((algebraic[66]-algebraic[119])+algebraic[177])/constants[3]
algebraic[168] = 1.00000+algebraic[132]/constants[351]+algebraic[149]/constants[353]+algebraic[162]/constants[355]
algebraic[176] = constants[39]*(((states[34]/algebraic[174])*states[13])/algebraic[33]-((states[29]/algebraic[44])*states[14])/algebraic[168])
rates[29] = (((((algebraic[73]-algebraic[117])+algebraic[176])-algebraic[87])-algebraic[175])-algebraic[177])/constants[3]
rates[38] = (algebraic[92]-algebraic[177])/constants[4]
rates[34] = (-algebraic[176]+algebraic[91]+algebraic[87]+algebraic[175]+algebraic[177])/constants[4]
algebraic[112] = exp(-(constants[290]-constants[7]*constants[2]*states[37])/constants[274])/1.00000
algebraic[142] = (((algebraic[112]*(power(algebraic[132], constants[7])))/(power(states[1], constants[7]-1.00000)))*algebraic[0])/(algebraic[21]*algebraic[33])
algebraic[143] = constants[46]*(algebraic[142]*states[8]*states[13]-states[0])
algebraic[163] = 1.00000+algebraic[132]/constants[296]+algebraic[149]/constants[311]+algebraic[162]/constants[325]
algebraic[165] = 1.00000+algebraic[132]/constants[332]+algebraic[149]/constants[336]+algebraic[162]/constants[340]
algebraic[178] = custom_piecewise([greater(states[50] , constants[249]) & greater(states[0] , constants[249]), constants[47]*((states[50]/algebraic[165])/(states[50]/algebraic[165]+(states[51]/algebraic[163])*exp((constants[12]*constants[2]*states[37])/constants[274]))-(states[8]/algebraic[21])/(states[8]/algebraic[21]+(states[0]/algebraic[0])*exp((constants[12]*(constants[2]-1.00000)*states[37])/constants[274])))*((states[50]/algebraic[165])/(states[50]/algebraic[165]+constants[11])) , True, 0.00000])
algebraic[144] = custom_piecewise([greater(fabs(states[37]) , constants[250]), (constants[51]*states[37]*(algebraic[132]*exp((constants[2]*states[37])/constants[274])-states[1]))/(exp((constants[2]*states[37])/constants[274])-1.00000) , True, (constants[51]*constants[274]*(algebraic[132]-states[1]))/constants[2]])
rates[37] = (((((4.00000*algebraic[139]+2.00000*algebraic[122]+4.00000*algebraic[141])-constants[7]*algebraic[143])-algebraic[178])-algebraic[144])+algebraic[158])/constants[21]
rates[0] = ((algebraic[81]+algebraic[143])-algebraic[178])/constants[3]
rates[8] = ((-algebraic[81]-algebraic[143])+algebraic[178])/constants[3]
algebraic[34] = ((states[13]*states[1])/constants[351])/algebraic[33]
algebraic[169] = ((states[14]*algebraic[132])/constants[351])/algebraic[168]
algebraic[179] = (constants[48]*(algebraic[169]*algebraic[132]-algebraic[34]*states[1]))/(constants[49]*(1.00000+algebraic[169]/constants[49])*(1.00000+algebraic[34]/constants[49]))
rates[13] = (((-algebraic[66]-algebraic[143])+algebraic[179])-algebraic[176])/constants[3]
rates[14] = (-algebraic[179]+algebraic[98]+algebraic[176])/constants[4]
algebraic[167] = 1.00000+algebraic[132]/constants[344]+algebraic[149]/constants[347]+algebraic[162]/constants[349]
algebraic[180] = (constants[381]*algebraic[163]*algebraic[167])/(algebraic[165]*algebraic[165])
algebraic[182] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[261]*(algebraic[180]*states[50]*states[50]-states[53]*states[51])])
rates[51] = (algebraic[96]+algebraic[178]+algebraic[182])/constants[4]
rates[50] = ((algebraic[95]-algebraic[178])-2.00000*algebraic[182])/constants[4]
rates[53] = (algebraic[97]+algebraic[182])/constants[4]
algebraic[189] = -((states[1]*states[1])/(constants[134]*constants[135]*(power(constants[282], 2.00000)))+(states[1]*states[0])/(constants[296]*constants[311]*(power(algebraic[0], 2.00000)))+(states[1]*states[8])/(constants[332]*constants[336]*(power(algebraic[21], 2.00000)))+(states[1]*states[10])/(constants[344]*constants[347]*(power(algebraic[26], 2.00000)))+(states[1]*states[11])/(constants[296]*constants[311]*(power(algebraic[29], 2.00000)))+(states[1]*states[12])/(constants[332]*constants[336]*(power(algebraic[31], 2.00000)))+(states[1]*states[13])/(constants[351]*constants[353]*(power(algebraic[33], 2.00000)))+(states[1]*states[21])/(constants[110]*constants[111]*(power(constants[277], 2.00000)))+(states[1]*algebraic[113])/(constants[107]*constants[108]*(power(constants[276], 2.00000)))+(states[1]*states[27])/(constants[113]*constants[114]*(power(constants[304], 2.00000)))+(states[1]*algebraic[118])/(constants[116]*constants[117]*(power(constants[320], 2.00000)))+(states[1]*states[17])/(constants[98]*constants[374]*(power(algebraic[46], 2.00000)))+(states[1]*states[22])/(constants[89]*constants[359]*(power(algebraic[36], 2.00000)))+(states[1]*states[20])/(constants[119]*constants[120]*(power(constants[278], 2.00000)))+(states[1]*states[15])/(constants[360]*constants[361]*(power(algebraic[37], 2.00000)))+(states[1]*states[24])/(constants[363]*constants[364]*(power(algebraic[40], 2.00000)))+(states[1]*states[25])/(constants[122]*constants[123]*(power(constants[306], 2.00000)))+(states[1]*states[23])/(constants[365]*constants[92]*(power(algebraic[41], 2.00000)))+(states[1]*states[18])/(constants[357]*constants[87]*(power(algebraic[35], 2.00000)))+(states[1]*states[26])/(constants[366]*constants[367]*(power(algebraic[42], 2.00000)))+(states[1]*states[28])/(constants[369]*constants[94]*(power(algebraic[43], 2.00000)))+(states[1]*states[29])/(constants[370]*constants[371]*(power(algebraic[44], 2.00000)))+(states[1]*states[31])/(constants[376]*constants[377]*(power(algebraic[48], 2.00000)))+(states[1]*states[30])/(constants[378]*constants[379]*(power(algebraic[49], 2.00000)))+(states[1]*states[3])/(constants[140]*constants[141]*(power(constants[284], 2.00000)))+(states[1]*states[2])/(constants[137]*constants[138]*(power(constants[283], 2.00000)))+(states[1]*states[56])/(constants[125]*constants[126]*(power(constants[279], 2.00000)))+(states[1]*constants[291])/(constants[128]*constants[129]*(power(constants[280], 2.00000)))+(states[1]*constants[310])/(constants[131]*constants[132]*(power(constants[281], 2.00000)))+(states[1]*states[19])/(constants[373]*constants[96]*(power(algebraic[45], 2.00000))))
algebraic[190] = -((states[1]*states[1])/(constants[134]*constants[136]*(power(constants[282], 2.00000)))+(states[1]*states[0])/(constants[296]*constants[325]*(power(algebraic[0], 2.00000)))+(states[1]*states[8])/(constants[332]*constants[340]*(power(algebraic[21], 2.00000)))+(states[1]*states[10])/(constants[344]*constants[349]*(power(algebraic[26], 2.00000)))+(states[1]*states[11])/(constants[296]*constants[325]*(power(algebraic[29], 2.00000)))+(states[1]*states[12])/(constants[332]*constants[340]*(power(algebraic[31], 2.00000)))+(states[1]*states[13])/(constants[351]*constants[355]*(power(algebraic[33], 2.00000)))+(states[1]*states[21])/(constants[110]*constants[112]*(power(constants[277], 2.00000)))+(states[1]*algebraic[113])/(constants[107]*constants[109]*(power(constants[276], 2.00000)))+(states[1]*states[27])/(constants[113]*constants[115]*(power(constants[304], 2.00000)))+(states[1]*algebraic[118])/(constants[116]*constants[118]*(power(constants[320], 2.00000)))+(states[1]*states[17])/(constants[98]*constants[99]*(power(algebraic[46], 2.00000)))+(states[1]*states[22])/(constants[89]*constants[90]*(power(algebraic[36], 2.00000)))+(states[1]*states[20])/(constants[119]*constants[121]*(power(constants[278], 2.00000)))+(states[1]*states[15])/(constants[360]*constants[362]*(power(algebraic[37], 2.00000)))+(states[1]*states[24])/(constants[363]*constants[91]*(power(algebraic[40], 2.00000)))+(states[1]*states[25])/(constants[122]*constants[124]*(power(constants[306], 2.00000)))+(states[1]*states[23])/(constants[365]*constants[93]*(power(algebraic[41], 2.00000)))+(states[1]*states[18])/(constants[357]*constants[88]*(power(algebraic[35], 2.00000)))+(states[1]*states[26])/(constants[366]*constants[368]*(power(algebraic[42], 2.00000)))+(states[1]*states[28])/(constants[369]*constants[95]*(power(algebraic[43], 2.00000)))+(states[1]*states[29])/(constants[370]*constants[372]*(power(algebraic[44], 2.00000)))+(states[1]*states[31])/(constants[376]*constants[102]*(power(algebraic[48], 2.00000)))+(states[1]*states[30])/(constants[378]*constants[103]*(power(algebraic[49], 2.00000)))+(states[1]*states[3])/(constants[140]*constants[142]*(power(constants[284], 2.00000)))+(states[1]*states[2])/(constants[137]*constants[139]*(power(constants[283], 2.00000)))+(states[1]*states[56])/(constants[125]*constants[127]*(power(constants[279], 2.00000)))+(states[1]*constants[291])/(constants[128]*constants[130]*(power(constants[280], 2.00000)))+(states[1]*constants[310])/(constants[131]*constants[133]*(power(constants[281], 2.00000)))+(states[1]*states[19])/(constants[373]*constants[97]*(power(algebraic[45], 2.00000))))
algebraic[193] = -((states[2]*states[1])/(constants[135]*constants[136]*(power(constants[282], 2.00000)))+(states[2]*states[0])/(constants[311]*constants[325]*(power(algebraic[0], 2.00000)))+(states[2]*states[8])/(constants[336]*constants[340]*(power(algebraic[21], 2.00000)))+(states[2]*states[10])/(constants[347]*constants[349]*(power(algebraic[26], 2.00000)))+(states[2]*states[11])/(constants[311]*constants[325]*(power(algebraic[29], 2.00000)))+(states[2]*states[12])/(constants[336]*constants[340]*(power(algebraic[31], 2.00000)))+(states[2]*states[13])/(constants[353]*constants[355]*(power(algebraic[33], 2.00000)))+(states[2]*states[21])/(constants[111]*constants[112]*(power(constants[277], 2.00000)))+(states[2]*algebraic[113])/(constants[108]*constants[109]*(power(constants[276], 2.00000)))+(states[2]*states[27])/(constants[114]*constants[115]*(power(constants[304], 2.00000)))+(states[2]*algebraic[118])/(constants[117]*constants[118]*(power(constants[320], 2.00000)))+(states[2]*states[17])/(constants[374]*constants[99]*(power(algebraic[46], 2.00000)))+(states[2]*states[22])/(constants[359]*constants[90]*(power(algebraic[36], 2.00000)))+(states[2]*states[20])/(constants[120]*constants[121]*(power(constants[278], 2.00000)))+(states[2]*states[15])/(constants[361]*constants[362]*(power(algebraic[37], 2.00000)))+(states[2]*states[24])/(constants[364]*constants[91]*(power(algebraic[40], 2.00000)))+(states[2]*states[25])/(constants[123]*constants[124]*(power(constants[306], 2.00000)))+(states[2]*states[23])/(constants[92]*constants[93]*(power(algebraic[41], 2.00000)))+(states[2]*states[18])/(constants[87]*constants[88]*(power(algebraic[35], 2.00000)))+(states[2]*states[26])/(constants[367]*constants[368]*(power(algebraic[42], 2.00000)))+(states[2]*states[28])/(constants[94]*constants[95]*(power(algebraic[43], 2.00000)))+(states[2]*states[29])/(constants[371]*constants[372]*(power(algebraic[44], 2.00000)))+(states[2]*states[31])/(constants[377]*constants[102]*(power(algebraic[48], 2.00000)))+(states[2]*states[30])/(constants[379]*constants[103]*(power(algebraic[49], 2.00000)))+(states[2]*states[3])/(constants[141]*constants[142]*(power(constants[284], 2.00000)))+(states[2]*states[2])/(constants[138]*constants[139]*(power(constants[283], 2.00000)))+(states[2]*states[56])/(constants[126]*constants[127]*(power(constants[279], 2.00000)))+(states[2]*constants[291])/(constants[129]*constants[130]*(power(constants[280], 2.00000)))+(states[2]*constants[310])/(constants[132]*constants[133]*(power(constants[281], 2.00000)))+(states[2]*states[19])/(constants[96]*constants[97]*(power(algebraic[45], 2.00000))))
algebraic[196] = -((states[3]*states[1])/(constants[136]*constants[135]*(power(constants[282], 2.00000)))+(states[3]*states[0])/(constants[325]*constants[311]*(power(algebraic[0], 2.00000)))+(states[3]*states[8])/(constants[340]*constants[336]*(power(algebraic[21], 2.00000)))+(states[3]*states[10])/(constants[349]*constants[347]*(power(algebraic[26], 2.00000)))+(states[3]*states[11])/(constants[325]*constants[311]*(power(algebraic[29], 2.00000)))+(states[3]*states[12])/(constants[340]*constants[336]*(power(algebraic[31], 2.00000)))+(states[3]*states[13])/(constants[355]*constants[353]*(power(algebraic[33], 2.00000)))+(states[3]*states[21])/(constants[112]*constants[111]*(power(constants[277], 2.00000)))+(states[3]*algebraic[113])/(constants[109]*constants[108]*(power(constants[276], 2.00000)))+(states[3]*states[27])/(constants[115]*constants[114]*(power(constants[304], 2.00000)))+(states[3]*algebraic[118])/(constants[118]*constants[117]*(power(constants[320], 2.00000)))+(states[3]*states[17])/(constants[99]*constants[374]*(power(algebraic[46], 2.00000)))+(states[3]*states[22])/(constants[90]*constants[359]*(power(algebraic[36], 2.00000)))+(states[3]*states[20])/(constants[121]*constants[120]*(power(constants[278], 2.00000)))+(states[3]*states[15])/(constants[362]*constants[361]*(power(algebraic[37], 2.00000)))+(states[3]*states[24])/(constants[91]*constants[364]*(power(algebraic[40], 2.00000)))+(states[3]*states[25])/(constants[124]*constants[123]*(power(constants[306], 2.00000)))+(states[3]*states[23])/(constants[93]*constants[92]*(power(algebraic[41], 2.00000)))+(states[3]*states[18])/(constants[88]*constants[87]*(power(algebraic[35], 2.00000)))+(states[3]*states[26])/(constants[368]*constants[367]*(power(algebraic[42], 2.00000)))+(states[3]*states[28])/(constants[95]*constants[94]*(power(algebraic[43], 2.00000)))+(states[3]*states[29])/(constants[372]*constants[371]*(power(algebraic[44], 2.00000)))+(states[3]*states[31])/(constants[102]*constants[377]*(power(algebraic[48], 2.00000)))+(states[3]*states[30])/(constants[103]*constants[379]*(power(algebraic[49], 2.00000)))+(states[3]*states[3])/(constants[142]*constants[141]*(power(constants[284], 2.00000)))+(states[3]*states[2])/(constants[139]*constants[138]*(power(constants[283], 2.00000)))+(states[3]*states[56])/(constants[127]*constants[126]*(power(constants[279], 2.00000)))+(states[3]*constants[291])/(constants[130]*constants[129]*(power(constants[280], 2.00000)))+(states[3]*constants[310])/(constants[133]*constants[132]*(power(constants[281], 2.00000)))+(states[3]*states[19])/(constants[97]*constants[96]*(power(algebraic[45], 2.00000))))
algebraic[197] = ((1.00000+states[1]/constants[134]+states[2]/constants[135])*states[1])/(constants[136]*(power(constants[282], 2.00000)))+((1.00000+states[1]/constants[296]+states[2]/constants[311])*states[0])/(constants[325]*(power(algebraic[0], 2.00000)))+((1.00000+states[1]/constants[332]+states[2]/constants[336])*states[8])/(constants[340]*(power(algebraic[21], 2.00000)))+((1.00000+states[1]/constants[344]+states[2]/constants[347])*states[10])/(constants[349]*(power(algebraic[26], 2.00000)))+((1.00000+states[1]/constants[296]+states[2]/constants[311])*states[11])/(constants[325]*(power(algebraic[29], 2.00000)))+((1.00000+states[1]/constants[332]+states[2]/constants[336])*states[12])/(constants[340]*(power(algebraic[31], 2.00000)))+((1.00000+states[1]/constants[351]+states[2]/constants[353])*states[13])/(constants[355]*(power(algebraic[33], 2.00000)))+((1.00000+states[1]/constants[110]+states[2]/constants[111])*states[21])/(constants[112]*(power(constants[277], 2.00000)))+((1.00000+states[1]/constants[107]+states[2]/constants[108])*algebraic[113])/(constants[109]*(power(constants[276], 2.00000)))+((1.00000+states[1]/constants[113]+states[2]/constants[114])*states[27])/(constants[115]*(power(constants[304], 2.00000)))+((1.00000+states[1]/constants[116]+states[2]/constants[117])*algebraic[118])/(constants[118]*(power(constants[320], 2.00000)))+((1.00000+states[1]/constants[98]+states[2]/constants[374])*states[17])/(constants[99]*(power(algebraic[46], 2.00000)))+((1.00000+states[1]/constants[89]+states[2]/constants[359])*states[22])/(constants[90]*(power(algebraic[36], 2.00000)))+((1.00000+states[1]/constants[119]+states[2]/constants[120])*states[20])/(constants[121]*(power(constants[278], 2.00000)))+((1.00000+states[1]/constants[360]+states[2]/constants[361])*states[15])/(constants[362]*(power(algebraic[37], 2.00000)))+((1.00000+states[1]/constants[363]+states[2]/constants[364])*states[24])/(constants[91]*(power(algebraic[40], 2.00000)))+((1.00000+states[1]/constants[122]+states[2]/constants[123])*states[25])/(constants[124]*(power(constants[306], 2.00000)))+((1.00000+states[1]/constants[365]+states[2]/constants[92])*states[23])/(constants[93]*(power(algebraic[41], 2.00000)))+((1.00000+states[1]/constants[357]+states[2]/constants[87])*states[18])/(constants[88]*(power(algebraic[35], 2.00000)))+((1.00000+states[1]/constants[366]+states[2]/constants[367])*states[26])/(constants[368]*(power(algebraic[42], 2.00000)))+((1.00000+states[1]/constants[369]+states[2]/constants[94])*states[28])/(constants[95]*(power(algebraic[43], 2.00000)))+((1.00000+states[1]/constants[370]+states[2]/constants[371])*states[29])/(constants[372]*(power(algebraic[44], 2.00000)))+((1.00000+states[1]/constants[376]+states[2]/constants[377])*states[31])/(constants[102]*(power(algebraic[48], 2.00000)))+((1.00000+states[1]/constants[378]+states[2]/constants[379])*states[30])/(constants[103]*(power(algebraic[49], 2.00000)))+((1.00000+states[1]/constants[140]+states[2]/constants[141])*states[3])/(constants[142]*(power(constants[284], 2.00000)))+((1.00000+states[1]/constants[137]+states[2]/constants[138])*states[2])/(constants[139]*(power(constants[283], 2.00000)))+((1.00000+states[1]/constants[125]+states[2]/constants[126])*states[56])/(constants[127]*(power(constants[279], 2.00000)))+((1.00000+states[1]/constants[128]+states[2]/constants[129])*constants[291])/(constants[130]*(power(constants[280], 2.00000)))+((1.00000+states[1]/constants[131]+states[2]/constants[132])*constants[310])/(constants[133]*(power(constants[281], 2.00000)))+((1.00000+states[1]/constants[373]+states[2]/constants[96])*states[19])/(constants[97]*(power(algebraic[45], 2.00000)))
algebraic[205] = 1.00000+algebraic[197]
algebraic[194] = ((1.00000+states[1]/constants[134]+states[3]/constants[136])*states[1])/(constants[135]*(power(constants[282], 2.00000)))+((1.00000+states[1]/constants[296]+states[3]/constants[325])*states[0])/(constants[311]*(power(algebraic[0], 2.00000)))+((1.00000+states[1]/constants[332]+states[3]/constants[340])*states[8])/(constants[336]*(power(algebraic[21], 2.00000)))+((1.00000+states[1]/constants[344]+states[3]/constants[349])*states[10])/(constants[347]*(power(algebraic[26], 2.00000)))+((1.00000+states[1]/constants[296]+states[3]/constants[325])*states[11])/(constants[311]*(power(algebraic[29], 2.00000)))+((1.00000+states[1]/constants[332]+states[3]/constants[340])*states[12])/(constants[336]*(power(algebraic[31], 2.00000)))+((1.00000+states[1]/constants[351]+states[3]/constants[355])*states[13])/(constants[353]*(power(algebraic[33], 2.00000)))+((1.00000+states[1]/constants[110]+states[3]/constants[112])*states[21])/(constants[111]*(power(constants[277], 2.00000)))+((1.00000+states[1]/constants[107]+states[3]/constants[109])*algebraic[113])/(constants[108]*(power(constants[276], 2.00000)))+((1.00000+states[1]/constants[113]+states[3]/constants[115])*states[27])/(constants[114]*(power(constants[304], 2.00000)))+((1.00000+states[1]/constants[116]+states[3]/constants[118])*algebraic[118])/(constants[117]*(power(constants[320], 2.00000)))+((1.00000+states[1]/constants[98]+states[3]/constants[99])*states[17])/(constants[374]*(power(algebraic[46], 2.00000)))+((1.00000+states[1]/constants[89]+states[3]/constants[90])*states[22])/(constants[359]*(power(algebraic[36], 2.00000)))+((1.00000+states[1]/constants[119]+states[3]/constants[121])*states[20])/(constants[120]*(power(constants[278], 2.00000)))+((1.00000+states[1]/constants[360]+states[3]/constants[362])*states[15])/(constants[361]*(power(algebraic[37], 2.00000)))+((1.00000+states[1]/constants[363]+states[3]/constants[91])*states[24])/(constants[364]*(power(algebraic[40], 2.00000)))+((1.00000+states[1]/constants[122]+states[3]/constants[124])*states[25])/(constants[123]*(power(constants[306], 2.00000)))+((1.00000+states[1]/constants[365]+states[3]/constants[93])*states[23])/(constants[92]*(power(algebraic[41], 2.00000)))+((1.00000+states[1]/constants[357]+states[3]/constants[88])*states[18])/(constants[87]*(power(algebraic[35], 2.00000)))+((1.00000+states[1]/constants[366]+states[3]/constants[368])*states[26])/(constants[367]*(power(algebraic[42], 2.00000)))+((1.00000+states[1]/constants[369]+states[3]/constants[95])*states[28])/(constants[94]*(power(algebraic[43], 2.00000)))+((1.00000+states[1]/constants[370]+states[3]/constants[372])*states[29])/(constants[371]*(power(algebraic[44], 2.00000)))+((1.00000+states[1]/constants[376]+states[3]/constants[102])*states[31])/(constants[377]*(power(algebraic[48], 2.00000)))+((1.00000+states[1]/constants[378]+states[3]/constants[103])*states[30])/(constants[379]*(power(algebraic[49], 2.00000)))+((1.00000+states[1]/constants[140]+states[3]/constants[142])*states[3])/(constants[141]*(power(constants[284], 2.00000)))+((1.00000+states[1]/constants[137]+states[3]/constants[139])*states[2])/(constants[138]*(power(constants[283], 2.00000)))+((1.00000+states[1]/constants[125]+states[3]/constants[127])*states[56])/(constants[126]*(power(constants[279], 2.00000)))+((1.00000+states[1]/constants[128]+states[3]/constants[130])*constants[291])/(constants[129]*(power(constants[280], 2.00000)))+((1.00000+states[1]/constants[131]+states[3]/constants[133])*constants[310])/(constants[132]*(power(constants[281], 2.00000)))+((1.00000+states[1]/constants[373]+states[3]/constants[97])*states[19])/(constants[96]*(power(algebraic[45], 2.00000)))
algebraic[204] = 1.00000+algebraic[194]
algebraic[191] = ((1.00000+states[2]/constants[135]+states[3]/constants[136])*states[1])/(constants[134]*(power(constants[282], 2.00000)))+((1.00000+states[2]/constants[311]+states[3]/constants[325])*states[0])/(constants[296]*(power(algebraic[0], 2.00000)))+((1.00000+states[2]/constants[336]+states[3]/constants[340])*states[8])/(constants[332]*(power(algebraic[21], 2.00000)))+((1.00000+states[2]/constants[347]+states[3]/constants[349])*states[10])/(constants[344]*(power(algebraic[26], 2.00000)))+((1.00000+states[2]/constants[311]+states[3]/constants[325])*states[11])/(constants[296]*(power(algebraic[29], 2.00000)))+((1.00000+states[2]/constants[336]+states[3]/constants[340])*states[12])/(constants[332]*(power(algebraic[31], 2.00000)))+((1.00000+states[2]/constants[353]+states[3]/constants[355])*states[13])/(constants[351]*(power(algebraic[33], 2.00000)))+((1.00000+states[2]/constants[111]+states[3]/constants[112])*states[21])/(constants[110]*(power(constants[277], 2.00000)))+((1.00000+states[2]/constants[108]+states[3]/constants[109])*algebraic[113])/(constants[107]*(power(constants[276], 2.00000)))+((1.00000+states[2]/constants[114]+states[3]/constants[115])*states[27])/(constants[113]*(power(constants[304], 2.00000)))+((1.00000+states[2]/constants[117]+states[3]/constants[118])*algebraic[118])/(constants[116]*(power(constants[320], 2.00000)))+((1.00000+states[2]/constants[374]+states[3]/constants[99])*states[17])/(constants[98]*(power(algebraic[46], 2.00000)))+((1.00000+states[2]/constants[359]+states[3]/constants[90])*states[22])/(constants[89]*(power(algebraic[36], 2.00000)))+((1.00000+states[2]/constants[120]+states[3]/constants[121])*states[20])/(constants[119]*(power(constants[278], 2.00000)))+((1.00000+states[2]/constants[361]+states[3]/constants[362])*states[15])/(constants[360]*(power(algebraic[37], 2.00000)))+((1.00000+states[2]/constants[364]+states[3]/constants[91])*states[24])/(constants[363]*(power(algebraic[40], 2.00000)))+((1.00000+states[2]/constants[123]+states[3]/constants[124])*states[25])/(constants[122]*(power(constants[306], 2.00000)))+((1.00000+states[2]/constants[92]+states[3]/constants[93])*states[23])/(constants[365]*(power(algebraic[41], 2.00000)))+((1.00000+states[2]/constants[87]+states[3]/constants[88])*states[18])/(constants[357]*(power(algebraic[35], 2.00000)))+((1.00000+states[2]/constants[367]+states[3]/constants[368])*states[26])/(constants[366]*(power(algebraic[42], 2.00000)))+((1.00000+states[2]/constants[94]+states[3]/constants[95])*states[28])/(constants[369]*(power(algebraic[43], 2.00000)))+((1.00000+states[2]/constants[371]+states[3]/constants[372])*states[29])/(constants[370]*(power(algebraic[44], 2.00000)))+((1.00000+states[2]/constants[377]+states[3]/constants[102])*states[31])/(constants[376]*(power(algebraic[48], 2.00000)))+((1.00000+states[2]/constants[379]+states[3]/constants[103])*states[30])/(constants[378]*(power(algebraic[49], 2.00000)))+((1.00000+states[2]/constants[141]+states[3]/constants[142])*states[3])/(constants[140]*(power(constants[284], 2.00000)))+((1.00000+states[2]/constants[138]+states[3]/constants[139])*states[2])/(constants[137]*(power(constants[283], 2.00000)))+((1.00000+states[2]/constants[126]+states[3]/constants[127])*states[56])/(constants[125]*(power(constants[279], 2.00000)))+((1.00000+states[2]/constants[129]+states[3]/constants[130])*constants[291])/(constants[128]*(power(constants[280], 2.00000)))+((1.00000+states[2]/constants[132]+states[3]/constants[133])*constants[310])/(constants[131]*(power(constants[281], 2.00000)))+((1.00000+states[2]/constants[96]+states[3]/constants[97])*states[19])/(constants[373]*(power(algebraic[45], 2.00000)))
algebraic[203] = 1.00000+algebraic[191]+constants[271]/(constants[272]*(power(1.00000+states[1]/constants[272], 2.00000)))
algebraic[195] = -((states[3]*states[1])/(constants[136]*constants[134]*(power(constants[282], 2.00000)))+(states[3]*states[0])/(constants[325]*constants[296]*(power(algebraic[0], 2.00000)))+(states[3]*states[8])/(constants[340]*constants[332]*(power(algebraic[21], 2.00000)))+(states[3]*states[10])/(constants[349]*constants[344]*(power(algebraic[26], 2.00000)))+(states[3]*states[11])/(constants[325]*constants[296]*(power(algebraic[29], 2.00000)))+(states[3]*states[12])/(constants[340]*constants[332]*(power(algebraic[31], 2.00000)))+(states[3]*states[13])/(constants[355]*constants[351]*(power(algebraic[33], 2.00000)))+(states[3]*states[21])/(constants[112]*constants[110]*(power(constants[277], 2.00000)))+(states[3]*algebraic[113])/(constants[109]*constants[107]*(power(constants[276], 2.00000)))+(states[3]*states[27])/(constants[115]*constants[113]*(power(constants[304], 2.00000)))+(states[3]*algebraic[118])/(constants[118]*constants[116]*(power(constants[320], 2.00000)))+(states[3]*states[17])/(constants[99]*constants[98]*(power(algebraic[46], 2.00000)))+(states[3]*states[22])/(constants[90]*constants[89]*(power(algebraic[36], 2.00000)))+(states[3]*states[20])/(constants[121]*constants[119]*(power(constants[278], 2.00000)))+(states[3]*states[15])/(constants[362]*constants[360]*(power(algebraic[37], 2.00000)))+(states[3]*states[24])/(constants[91]*constants[363]*(power(algebraic[40], 2.00000)))+(states[3]*states[25])/(constants[124]*constants[122]*(power(constants[306], 2.00000)))+(states[3]*states[23])/(constants[93]*constants[365]*(power(algebraic[41], 2.00000)))+(states[3]*states[18])/(constants[88]*constants[357]*(power(algebraic[35], 2.00000)))+(states[3]*states[26])/(constants[368]*constants[366]*(power(algebraic[42], 2.00000)))+(states[3]*states[28])/(constants[95]*constants[369]*(power(algebraic[43], 2.00000)))+(states[3]*states[29])/(constants[372]*constants[370]*(power(algebraic[44], 2.00000)))+(states[3]*states[31])/(constants[102]*constants[376]*(power(algebraic[48], 2.00000)))+(states[3]*states[30])/(constants[103]*constants[378]*(power(algebraic[49], 2.00000)))+(states[3]*states[3])/(constants[142]*constants[140]*(power(constants[284], 2.00000)))+(states[3]*states[2])/(constants[139]*constants[137]*(power(constants[283], 2.00000)))+(states[3]*states[56])/(constants[127]*constants[125]*(power(constants[279], 2.00000)))+(states[3]*constants[291])/(constants[130]*constants[128]*(power(constants[280], 2.00000)))+(states[3]*constants[310])/(constants[133]*constants[131]*(power(constants[281], 2.00000)))+(states[3]*states[19])/(constants[97]*constants[373]*(power(algebraic[45], 2.00000))))
algebraic[192] = -((states[2]*states[1])/(constants[135]*constants[134]*(power(constants[282], 2.00000)))+(states[2]*states[0])/(constants[311]*constants[296]*(power(algebraic[0], 2.00000)))+(states[2]*states[8])/(constants[336]*constants[332]*(power(algebraic[21], 2.00000)))+(states[2]*states[10])/(constants[347]*constants[344]*(power(algebraic[26], 2.00000)))+(states[2]*states[11])/(constants[311]*constants[296]*(power(algebraic[29], 2.00000)))+(states[2]*states[12])/(constants[336]*constants[332]*(power(algebraic[31], 2.00000)))+(states[2]*states[13])/(constants[353]*constants[351]*(power(algebraic[33], 2.00000)))+(states[2]*states[21])/(constants[111]*constants[110]*(power(constants[277], 2.00000)))+(states[2]*algebraic[113])/(constants[108]*constants[107]*(power(constants[276], 2.00000)))+(states[2]*states[27])/(constants[114]*constants[113]*(power(constants[304], 2.00000)))+(states[2]*algebraic[118])/(constants[117]*constants[116]*(power(constants[320], 2.00000)))+(states[2]*states[17])/(constants[374]*constants[98]*(power(algebraic[46], 2.00000)))+(states[2]*states[22])/(constants[359]*constants[89]*(power(algebraic[36], 2.00000)))+(states[2]*states[20])/(constants[120]*constants[119]*(power(constants[278], 2.00000)))+(states[2]*states[15])/(constants[361]*constants[360]*(power(algebraic[37], 2.00000)))+(states[2]*states[24])/(constants[364]*constants[363]*(power(algebraic[40], 2.00000)))+(states[2]*states[25])/(constants[123]*constants[122]*(power(constants[306], 2.00000)))+(states[2]*states[23])/(constants[92]*constants[365]*(power(algebraic[41], 2.00000)))+(states[2]*states[18])/(constants[87]*constants[357]*(power(algebraic[35], 2.00000)))+(states[2]*states[26])/(constants[367]*constants[366]*(power(algebraic[42], 2.00000)))+(states[2]*states[28])/(constants[94]*constants[369]*(power(algebraic[43], 2.00000)))+(states[2]*states[29])/(constants[371]*constants[370]*(power(algebraic[44], 2.00000)))+(states[2]*states[31])/(constants[377]*constants[376]*(power(algebraic[48], 2.00000)))+(states[2]*states[30])/(constants[379]*constants[378]*(power(algebraic[49], 2.00000)))+(states[2]*states[3])/(constants[141]*constants[140]*(power(constants[284], 2.00000)))+(states[2]*states[2])/(constants[138]*constants[137]*(power(constants[283], 2.00000)))+(states[2]*states[56])/(constants[126]*constants[125]*(power(constants[279], 2.00000)))+(states[2]*constants[291])/(constants[129]*constants[128]*(power(constants[280], 2.00000)))+(states[2]*constants[310])/(constants[132]*constants[131]*(power(constants[281], 2.00000)))+(states[2]*states[19])/(constants[96]*constants[373]*(power(algebraic[45], 2.00000))))
algebraic[206] = (((algebraic[203]*algebraic[196]*algebraic[193]+algebraic[205]*algebraic[189]*algebraic[192]+algebraic[204]*algebraic[190]*algebraic[195])-algebraic[204]*algebraic[205]*algebraic[203])-algebraic[190]*algebraic[196]*algebraic[192])-algebraic[189]*algebraic[193]*algebraic[195]
algebraic[145] = ((algebraic[114]+algebraic[115]+algebraic[116]+algebraic[117])-algebraic[139])/constants[3]
algebraic[146] = -algebraic[145]
algebraic[147] = ((algebraic[119]+algebraic[139])-algebraic[122])/constants[3]
algebraic[148] = -algebraic[147]
algebraic[183] = ((algebraic[81]+algebraic[143])-algebraic[178])/constants[3]
algebraic[184] = ((-algebraic[81]-algebraic[143])+algebraic[178])/constants[3]
algebraic[123] = (algebraic[66]-algebraic[81])/constants[3]
algebraic[124] = (-algebraic[66]+algebraic[81])/constants[3]
algebraic[185] = (((-algebraic[66]-algebraic[143])+algebraic[179])-algebraic[176])/constants[3]
algebraic[159] = (-algebraic[114]+algebraic[155])/constants[3]
algebraic[125] = (algebraic[114]-algebraic[56])/constants[3]
algebraic[186] = (-algebraic[59]+algebraic[175]+algebraic[56])/constants[3]
algebraic[126] = (algebraic[59]-algebraic[115])/constants[3]
algebraic[127] = (((algebraic[115]-algebraic[116])-algebraic[86])+algebraic[87])/constants[3]
algebraic[128] = (algebraic[116]-algebraic[66])/constants[3]
algebraic[129] = ((-algebraic[114]-algebraic[116])+algebraic[66]+algebraic[56])/constants[3]
algebraic[187] = ((algebraic[66]-algebraic[119])+algebraic[177])/constants[3]
algebraic[130] = (algebraic[119]-algebraic[73])/constants[3]
algebraic[188] = (((((algebraic[73]-algebraic[117])+algebraic[176])-algebraic[87])-algebraic[175])-algebraic[177])/constants[3]
algebraic[131] = (-algebraic[56]+algebraic[117]+algebraic[86])/constants[3]
algebraic[160] = ((algebraic[86]+algebraic[156])-algebraic[158])/constants[3]
algebraic[161] = (-algebraic[86]+algebraic[158])/constants[3]
algebraic[181] = constants[50]*(algebraic[162]*states[1]-states[3]*algebraic[132])
algebraic[198] = (((((((((-algebraic[114]+2.00000*algebraic[56])-algebraic[116])+algebraic[66]+algebraic[117]+algebraic[155]+algebraic[156]+algebraic[175])-algebraic[158])-5.00000*algebraic[139])-2.00000*algebraic[122])-4.00000*algebraic[141])+(constants[7]-1.00000)*algebraic[143]+2.00000*algebraic[179]+algebraic[144])-algebraic[181])/constants[3]
algebraic[200] = -((states[1]*constants[292])/(constants[134]*constants[282])+(states[1]*algebraic[183])/(constants[296]*algebraic[0])+(states[1]*algebraic[184])/(constants[332]*algebraic[21])+(states[1]*constants[262])/(constants[344]*algebraic[26])+(states[1]*algebraic[123])/(constants[296]*algebraic[29])+(states[1]*algebraic[124])/(constants[332]*algebraic[31])+(states[1]*algebraic[185])/(constants[351]*algebraic[33])+(states[1]*algebraic[145])/(constants[110]*constants[277])+(states[1]*algebraic[146])/(constants[107]*constants[276])+(states[1]*algebraic[147])/(constants[113]*constants[304])+(states[1]*algebraic[148])/(constants[116]*constants[320])+(states[1]*algebraic[159])/(constants[98]*algebraic[46])+(states[1]*algebraic[131])/(constants[89]*algebraic[36])+(states[1]*algebraic[125])/(constants[119]*constants[278])+(states[1]*algebraic[186])/(constants[360]*algebraic[37])+(states[1]*algebraic[126])/(constants[363]*algebraic[40])+(states[1]*algebraic[127])/(constants[122]*constants[306])+(states[1]*algebraic[128])/(constants[365]*algebraic[41])+(states[1]*algebraic[129])/(constants[357]*algebraic[35])+(states[1]*algebraic[187])/(constants[366]*algebraic[42])+(states[1]*algebraic[130])/(constants[369]*algebraic[43])+(states[1]*algebraic[188])/(constants[370]*algebraic[44])+(states[1]*algebraic[160])/(constants[376]*algebraic[48])+(states[1]*algebraic[161])/(constants[378]*algebraic[49])+(states[1]*constants[294])/(constants[140]*constants[284])+(states[1]*constants[293])/(constants[137]*constants[283])+(states[1]*constants[263])/(constants[125]*constants[279])+(states[1]*constants[270])/(constants[128]*constants[280])+(states[1]*constants[269])/(constants[131]*constants[281])+(states[1]*constants[264])/(constants[373]*algebraic[45]))+algebraic[198]
algebraic[201] = -((states[2]*constants[292])/(constants[135]*constants[282])+(states[2]*algebraic[183])/(constants[311]*algebraic[0])+(states[2]*algebraic[184])/(constants[336]*algebraic[21])+(states[2]*constants[262])/(constants[347]*algebraic[26])+(states[2]*algebraic[123])/(constants[311]*algebraic[29])+(states[2]*algebraic[124])/(constants[336]*algebraic[31])+(states[2]*algebraic[185])/(constants[353]*algebraic[33])+(states[2]*algebraic[145])/(constants[111]*constants[277])+(states[2]*algebraic[146])/(constants[108]*constants[276])+(states[2]*algebraic[147])/(constants[114]*constants[304])+(states[2]*algebraic[148])/(constants[117]*constants[320])+(states[2]*algebraic[159])/(constants[374]*algebraic[46])+(states[2]*algebraic[131])/(constants[359]*algebraic[36])+(states[2]*algebraic[125])/(constants[120]*constants[278])+(states[2]*algebraic[186])/(constants[361]*algebraic[37])+(states[2]*algebraic[126])/(constants[364]*algebraic[40])+(states[2]*algebraic[127])/(constants[123]*constants[306])+(states[2]*algebraic[128])/(constants[92]*algebraic[41])+(states[2]*algebraic[129])/(constants[87]*algebraic[35])+(states[2]*algebraic[187])/(constants[367]*algebraic[42])+(states[2]*algebraic[130])/(constants[94]*algebraic[43])+(states[2]*algebraic[188])/(constants[371]*algebraic[44])+(states[2]*algebraic[160])/(constants[377]*algebraic[48])+(states[2]*algebraic[161])/(constants[379]*algebraic[49])+(states[2]*constants[294])/(constants[141]*constants[284])+(states[2]*constants[293])/(constants[138]*constants[283])+(states[2]*constants[263])/(constants[126]*constants[279])+(states[2]*constants[270])/(constants[129]*constants[280])+(states[2]*constants[269])/(constants[132]*constants[281])+(states[2]*constants[264])/(constants[96]*algebraic[45]))+constants[295]
algebraic[199] = algebraic[181]/constants[3]
algebraic[202] = -((states[3]*constants[292])/(constants[136]*constants[282])+(states[3]*algebraic[183])/(constants[325]*algebraic[0])+(states[3]*algebraic[184])/(constants[340]*algebraic[21])+(states[3]*constants[262])/(constants[349]*algebraic[26])+(states[3]*algebraic[123])/(constants[325]*algebraic[29])+(states[3]*algebraic[124])/(constants[340]*algebraic[31])+(states[3]*algebraic[185])/(constants[355]*algebraic[33])+(states[3]*algebraic[145])/(constants[112]*constants[277])+(states[3]*algebraic[146])/(constants[109]*constants[276])+(states[3]*algebraic[147])/(constants[115]*constants[304])+(states[3]*algebraic[148])/(constants[118]*constants[320])+(states[3]*algebraic[159])/(constants[99]*algebraic[46])+(states[3]*algebraic[131])/(constants[90]*algebraic[36])+(states[3]*algebraic[125])/(constants[121]*constants[278])+(states[3]*algebraic[186])/(constants[362]*algebraic[37])+(states[3]*algebraic[126])/(constants[91]*algebraic[40])+(states[3]*algebraic[127])/(constants[124]*constants[306])+(states[3]*algebraic[128])/(constants[93]*algebraic[41])+(states[3]*algebraic[129])/(constants[88]*algebraic[35])+(states[3]*algebraic[187])/(constants[368]*algebraic[42])+(states[3]*algebraic[130])/(constants[95]*algebraic[43])+(states[3]*algebraic[188])/(constants[372]*algebraic[44])+(states[3]*algebraic[160])/(constants[102]*algebraic[48])+(states[3]*algebraic[161])/(constants[103]*algebraic[49])+(states[3]*constants[294])/(constants[142]*constants[284])+(states[3]*constants[293])/(constants[139]*constants[283])+(states[3]*constants[263])/(constants[127]*constants[279])+(states[3]*constants[270])/(constants[130]*constants[280])+(states[3]*constants[269])/(constants[133]*constants[281])+(states[3]*constants[264])/(constants[97]*algebraic[45]))+algebraic[199]
rates[1] = ((algebraic[196]*algebraic[193]-algebraic[204]*algebraic[205])*algebraic[200]+(algebraic[205]*algebraic[189]-algebraic[190]*algebraic[196])*algebraic[201]+(algebraic[204]*algebraic[190]-algebraic[189]*algebraic[193])*algebraic[202])/algebraic[206]
rates[2] = ((algebraic[205]*algebraic[192]-algebraic[195]*algebraic[193])*algebraic[200]+(algebraic[195]*algebraic[190]-algebraic[203]*algebraic[205])*algebraic[201]+(algebraic[203]*algebraic[193]-algebraic[190]*algebraic[192])*algebraic[202])/algebraic[206]
rates[3] = ((algebraic[204]*algebraic[195]-algebraic[196]*algebraic[192])*algebraic[200]+(algebraic[203]*algebraic[196]-algebraic[195]*algebraic[189])*algebraic[201]+(algebraic[192]*algebraic[189]-algebraic[203]*algebraic[204])*algebraic[202])/algebraic[206]
return(rates)
def computeAlgebraic(constants, states, voi):
algebraic = array([[0.0] * len(voi)] * sizeAlgebraic)
states = array(states)
voi = array(voi)
algebraic[16] = constants[5]*constants[10]*(states[48]-states[49])
algebraic[15] = constants[5]*constants[246]*(states[46]-states[47])
algebraic[0] = 1.00000+states[1]/constants[296]+states[2]/constants[311]+states[3]/constants[325]
algebraic[29] = algebraic[0]
algebraic[21] = 1.00000+states[1]/constants[332]+states[2]/constants[336]+states[3]/constants[340]
algebraic[31] = algebraic[21]
algebraic[33] = 1.00000+states[1]/constants[351]+states[2]/constants[353]+states[3]/constants[355]
algebraic[35] = 1.00000+states[1]/constants[357]
algebraic[41] = 1.00000+states[1]/constants[365]
algebraic[42] = 1.00000+states[1]/constants[366]+states[2]/constants[367]+states[3]/constants[368]
algebraic[64] = (constants[319]*algebraic[35]*algebraic[42]*algebraic[29])/(states[1]*algebraic[31]*algebraic[41]*algebraic[33])
algebraic[65] = (constants[28]*constants[186]*constants[181]*constants[182])/(algebraic[64]*constants[177]*constants[178]*constants[185])
algebraic[66] = (constants[28]*algebraic[65]*(states[12]*states[23]*states[13]-(states[18]*states[26]*states[11])/algebraic[64]))/(algebraic[65]*constants[177]*constants[178]*constants[185]+algebraic[65]*constants[178]*constants[185]*states[12]+algebraic[65]*constants[177]*constants[184]*states[13]+algebraic[65]*constants[185]*states[12]*states[23]+algebraic[65]*constants[184]*states[12]*states[13]+algebraic[65]*constants[183]*states[23]*states[13]+algebraic[65]*states[12]*states[23]*states[13]+(constants[28]*constants[182]*constants[187]*states[18])/algebraic[64]+(constants[28]*constants[181]*constants[186]*states[11])/algebraic[64]+(constants[28]*constants[188]*states[18]*states[26])/algebraic[64]+(constants[28]*constants[187]*states[18]*states[11])/algebraic[64]+(constants[28]*constants[186]*states[26]*states[11])/algebraic[64]+(constants[28]*states[18]*states[26]*states[11])/algebraic[64]+(constants[28]*constants[187]*constants[182]*states[12]*states[18])/(constants[177]*algebraic[64])+(algebraic[65]*constants[177]*constants[184]*states[13]*states[11])/constants[182]+(constants[28]*constants[187]*constants[182]*states[12]*states[23]*states[18])/(constants[177]*constants[178]*algebraic[64])+(algebraic[65]*constants[183]*states[23]*states[13]*states[11])/constants[182]+(constants[28]*constants[188]*states[12]*states[18]*states[26])/(constants[177]*algebraic[64])+(algebraic[65]*constants[177]*constants[184]*states[13]*states[26]*states[11])/(constants[181]*constants[182])+(constants[28]*constants[182]*constants[187]*states[12]*states[23]*states[13]*states[18])/(constants[177]*constants[178]*constants[179]*algebraic[64])+(constants[28]*constants[180]*constants[188]*states[12]*states[23]*states[13]*states[26])/(constants[177]*constants[178]*constants[179]*algebraic[64])+(constants[28]*constants[188]*states[12]*states[23]*states[18]*states[26])/(constants[177]*constants[178]*algebraic[64])+(algebraic[65]*constants[183]*states[23]*states[13]*states[26]*states[11])/(constants[181]*constants[182])+(algebraic[65]*constants[183]*constants[179]*states[23]*states[18]*states[26]*states[11])/(constants[180]*constants[181]*constants[182])+(algebraic[65]*constants[177]*constants[184]*states[13]*states[18]*states[26]*states[11])/(constants[180]*constants[181]*constants[182])+(constants[28]*constants[188]*states[12]*states[23]*states[13]*states[18]*states[26])/(constants[177]*constants[178]*constants[179]*algebraic[64])+(algebraic[65]*constants[183]*states[23]*states[13]*states[18]*states[26]*states[11])/(constants[180]*constants[181]*constants[182]))
algebraic[77] = (constants[312]*algebraic[31]*algebraic[0])/(algebraic[29]*algebraic[21])
algebraic[78] = (constants[32]*constants[223]*constants[218])/(algebraic[77]*constants[216]*constants[221])
algebraic[26] = 1.00000+states[1]/constants[344]+states[2]/constants[347]+states[3]/constants[349]
algebraic[28] = (states[10]*(1.00000+states[1]/constants[344]))/algebraic[26]
algebraic[79] = 1.00000+algebraic[28]/constants[224]
algebraic[80] = custom_piecewise([greater(states[11] , constants[225]) & greater(states[8] , constants[225]), (((constants[32]*algebraic[78])/algebraic[79])*states[11]*states[8])/(algebraic[78]*constants[221]*states[11]+algebraic[78]*constants[220]*states[8]+(constants[32]*constants[223]*states[12])/algebraic[77]+(constants[32]*constants[222]*states[0])/algebraic[77]+algebraic[78]*states[11]*states[8]+(constants[32]*constants[223]*states[11]*states[12])/(algebraic[77]*constants[216])+(constants[32]*states[12]*states[0])/algebraic[77]+(algebraic[78]*constants[220]*states[8]*states[0])/constants[219]) , True, 0.00000])
algebraic[81] = custom_piecewise([greater(states[12] , constants[225]) & greater(states[0] , constants[225]), algebraic[80]-((((constants[32]*algebraic[78])/algebraic[79])*states[12]*states[0])/algebraic[77])/(algebraic[78]*constants[221]*states[11]+algebraic[78]*constants[220]*states[8]+(constants[32]*constants[223]*states[12])/algebraic[77]+(constants[32]*constants[222]*states[0])/algebraic[77]+algebraic[78]*states[11]*states[8]+(constants[32]*constants[223]*states[11]*states[12])/(algebraic[77]*constants[216])+(constants[32]*states[12]*states[0])/algebraic[77]+(algebraic[78]*constants[220]*states[8]*states[0])/constants[219]) , True, algebraic[80]])
algebraic[87] = (constants[37]*(states[35]*states[29]-states[25]*states[34]))/(constants[236]*constants[239]*(2.00000+states[34]/constants[238]+states[29]/constants[239]+states[35]/constants[236]+states[25]/constants[237]+(states[34]*states[25])/(constants[238]*constants[237])+(states[29]*states[35])/(constants[239]*constants[236])))
algebraic[89] = constants[5]*constants[10]*(states[42]-states[35])
algebraic[88] = constants[5]*constants[10]*(states[39]-states[32])
algebraic[90] = constants[5]*constants[10]*(states[40]-states[16])
algebraic[91] = constants[5]*constants[10]*(states[41]-states[34])
algebraic[92] = constants[5]*constants[10]*(states[43]-states[38])
algebraic[93] = constants[5]*constants[10]*(states[44]-states[33])
algebraic[94] = constants[5]*constants[10]*(states[45]-states[36])
algebraic[98] = constants[5]*constants[8]*(states[54]-states[14])
algebraic[17] = 1.00000+states[5]/constants[296]+states[6]/constants[311]+states[7]/constants[325]
algebraic[20] = ((states[4]*states[6])/constants[311])/algebraic[17]
algebraic[22] = 1.00000+states[5]/constants[332]+states[6]/constants[336]+states[7]/constants[340]
algebraic[25] = ((states[9]*states[6])/constants[336])/algebraic[22]
algebraic[47] = 1.00000+states[5]/constants[375]
algebraic[99] = (constants[324]*algebraic[47]*algebraic[22])/(states[5]*algebraic[17]*constants[275])
algebraic[101] = (algebraic[99]*constants[311]*algebraic[17])/(constants[336]*algebraic[22])
algebraic[103] = ((constants[273]/(constants[252]*constants[255]))*(algebraic[20]*states[57]-(states[58]*algebraic[25])/algebraic[101]))/(1.00000+algebraic[20]/constants[251]+states[57]/constants[252]+(algebraic[20]*states[57])/(constants[252]*constants[255])+states[58]/constants[253]+algebraic[25]/constants[254]+(states[58]*algebraic[25])/(constants[254]*constants[257])+(states[58]*states[57])/(constants[259]*constants[252]))
algebraic[105] = constants[18]-states[59]
algebraic[107] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[260]*(((constants[309]*states[9])/algebraic[22])*states[59]*states[5]-(states[4]/algebraic[17])*algebraic[105])])
algebraic[96] = constants[5]*constants[9]*(states[4]-states[51])
algebraic[27] = 1.00000+states[5]/constants[344]+states[6]/constants[347]+states[7]/constants[349]
algebraic[109] = (constants[381]*algebraic[17]*algebraic[27])/(algebraic[22]*algebraic[22])
algebraic[111] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[261]*(algebraic[109]*states[9]*states[9]-states[52]*states[4])])
algebraic[95] = constants[5]*constants[9]*(states[9]-states[50])
algebraic[97] = constants[5]*constants[9]*(states[52]-states[53])
algebraic[113] = constants[14]-states[21]
algebraic[45] = 1.00000+states[1]/constants[373]
algebraic[46] = 1.00000+states[2]/constants[374]
algebraic[50] = (constants[314]*algebraic[45]*constants[278]*constants[277])/(states[1]*algebraic[46]*algebraic[35]*constants[276])
algebraic[51] = 1.00000+states[20]/constants[146]
algebraic[52] = 1.00000+states[21]/constants[147]
algebraic[114] = custom_piecewise([greater(states[17] , constants[148]) & greater(states[18] , constants[148]) & greater(algebraic[113] , constants[148]), (constants[22]*(states[17]*states[18]*algebraic[113]-(states[19]*states[20]*states[21])/algebraic[50]))/(constants[145]*algebraic[52]*states[17]*states[18]+constants[144]*algebraic[51]*states[17]*algebraic[113]+constants[143]*states[18]*algebraic[113]+states[17]*states[18]*algebraic[113]) , True, 0.00000])
algebraic[36] = 1.00000+states[2]/constants[359]
algebraic[37] = 1.00000+states[1]/constants[360]+states[2]/constants[361]+states[3]/constants[362]
algebraic[53] = (constants[315]*algebraic[37]*algebraic[35])/((power(states[1], 2.00000))*constants[278]*algebraic[36])
algebraic[38] = (states[15]*(1.00000+states[1]/constants[360]))/algebraic[37]
algebraic[54] = 1.00000+algebraic[38]/constants[152]
algebraic[18] = (states[0]*(1.00000+states[1]/constants[296]))/algebraic[0]
algebraic[23] = (states[8]*(1.00000+states[1]/constants[332]))/algebraic[21]
algebraic[55] = 1.00000+algebraic[18]/constants[153]+algebraic[23]/constants[154]+algebraic[28]/constants[155]+states[18]/constants[156]+states[23]/constants[157]
algebraic[56] = (constants[23]*(states[20]*states[22]-(states[18]*states[15])/algebraic[53]))/(constants[151]*constants[150]*algebraic[54]+constants[149]*algebraic[54]*states[20]+constants[150]*algebraic[55]*states[22]+states[22]*states[20])
algebraic[40] = 1.00000+states[1]/constants[363]+states[2]/constants[364]
algebraic[57] = (constants[316]*algebraic[40])/algebraic[37]
algebraic[58] = (constants[24]*constants[159])/(constants[158]*algebraic[57])
algebraic[59] = (constants[24]*algebraic[58]*(states[15]-states[24]/algebraic[57]))/(constants[158]*algebraic[58]+algebraic[58]*states[15]+(constants[24]*states[24])/algebraic[57])
algebraic[60] = (constants[317]*algebraic[45]*constants[306]*constants[277])/((power(states[1], 2.00000))*constants[276]*algebraic[40])
algebraic[61] = 1.00000+(constants[166]/algebraic[23])*(1.00000+algebraic[18]/constants[165])
algebraic[115] = custom_piecewise([greater(algebraic[113] , constants[167]) & greater(states[24] , constants[167]), (constants[25]*(1.00000-(states[19]*states[25]*states[21])/(algebraic[60]*algebraic[113]*states[24])))/(1.00000+(power(constants[161]/states[24], constants[164]))*algebraic[61]+(constants[160]/algebraic[113])*(1.00000+(power(constants[162]/states[24], constants[164]))*algebraic[61]+(states[21]*algebraic[61])/constants[163])) , True, 0.00000])
algebraic[62] = (constants[318]*algebraic[45]*algebraic[41]*constants[277])/(states[1]*algebraic[35]*constants[306]*constants[276])
algebraic[63] = 1.00000+(constants[175]/algebraic[23])*(1.00000+algebraic[18]/constants[174])
algebraic[116] = custom_piecewise([greater(states[25] , constants[176]) & greater(states[18] , constants[176]) & greater(algebraic[113] , constants[176]), ((constants[26]*(1.00000-(states[19]*states[23]*states[21])/(algebraic[62]*states[25]*states[18]*algebraic[113])))/(1.00000+(constants[168]/states[25])*algebraic[63]+(constants[169]/states[18])*(1.00000+states[23]/constants[171])+(constants[170]/algebraic[113])*(1.00000+states[21]/constants[172])))/(1.00000+states[21]/constants[173]) , True, 0.00000])
algebraic[48] = 1.00000+states[1]/constants[376]+states[2]/constants[377]
algebraic[49] = 1.00000+states[1]/constants[378]+states[2]/constants[379]
algebraic[82] = (constants[328]*algebraic[48]*algebraic[36])/(algebraic[49]*constants[306])
algebraic[83] = (constants[33]*constants[233]*constants[228])/(algebraic[82]*constants[226]*constants[231])
algebraic[84] = 1.00000+states[25]/constants[234]
algebraic[85] = custom_piecewise([greater(states[30] , constants[235]) & greater(states[25] , constants[235]), (constants[33]*algebraic[83]*states[30]*states[25])/(algebraic[83]*constants[231]*states[30]+algebraic[83]*constants[230]*algebraic[84]*states[25]+(constants[33]*constants[233]*algebraic[84]*states[22])/algebraic[82]+(constants[33]*constants[232]*states[31])/algebraic[82]+algebraic[83]*states[30]*states[25]+(constants[33]*constants[233]*states[30]*states[22])/(algebraic[82]*constants[226])+(constants[33]*states[22]*states[31])/algebraic[82]+(algebraic[83]*constants[230]*states[25]*states[31])/constants[229]) , True, 0.00000])
algebraic[86] = custom_piecewise([greater(states[22] , constants[235]) & greater(states[31] , constants[235]), algebraic[85]-((constants[33]*algebraic[83]*states[22]*states[31])/algebraic[82])/(algebraic[83]*constants[231]*states[30]+algebraic[83]*constants[230]*algebraic[84]*states[25]+(constants[33]*constants[233]*algebraic[84]*states[22])/algebraic[82]+(constants[33]*constants[232]*states[31])/algebraic[82]+algebraic[83]*states[30]*states[25]+(constants[33]*constants[233]*states[30]*states[22])/(algebraic[82]*constants[226])+(constants[33]*states[22]*states[31])/algebraic[82]+(algebraic[83]*constants[230]*states[25]*states[31])/constants[229]) , True, algebraic[85]])
algebraic[44] = 1.00000+states[1]/constants[370]+states[2]/constants[371]+states[3]/constants[372]
algebraic[74] = (constants[321]*algebraic[36]*constants[277])/(states[1]*algebraic[44]*constants[276])
algebraic[75] = (constants[31]*constants[208]*constants[211])/(algebraic[74]*constants[205]*constants[210])
algebraic[76] = 1.00000+algebraic[18]/constants[213]+algebraic[23]/constants[214]+algebraic[28]/constants[215]
algebraic[117] = (constants[31]*algebraic[75]*(algebraic[113]*states[29]-(states[22]*states[21])/algebraic[74]))/(algebraic[75]*constants[205]*constants[210]*algebraic[76]+algebraic[75]*constants[210]*algebraic[113]+algebraic[75]*constants[209]*algebraic[76]*states[29]+(constants[31]*constants[212]*algebraic[76]*states[22])/algebraic[74]+(constants[31]*constants[211]*states[21])/algebraic[74]+algebraic[75]*algebraic[113]*states[29]+(constants[31]*constants[212]*algebraic[113]*states[22])/(algebraic[74]*constants[205])+(constants[31]*states[22]*states[21])/algebraic[74]+(algebraic[75]*constants[209]*states[29]*states[21])/constants[208]+(algebraic[75]*algebraic[113]*states[29]*states[22])/constants[207]+(constants[31]*states[29]*states[22]*states[21])/(constants[206]*algebraic[74]))
algebraic[118] = constants[15]-states[27]
algebraic[43] = 1.00000+states[1]/constants[369]
algebraic[67] = (constants[334]*algebraic[43]*constants[304])/(algebraic[42]*constants[320])
algebraic[68] = (constants[29]*constants[193]*constants[190])/(algebraic[67]*constants[189]*constants[192])
algebraic[69] = (1.00000+states[22]/constants[195]+states[26]/constants[196]+states[28]/constants[197])/(1.00000+states[26]/constants[196]+states[28]/constants[197])
algebraic[119] = (constants[29]*algebraic[68]*(states[26]*algebraic[118]-(states[27]*states[28])/algebraic[67]))/(algebraic[68]*constants[189]*constants[192]*algebraic[69]+algebraic[68]*constants[192]*states[26]+algebraic[68]*constants[191]*algebraic[69]*algebraic[118]+(constants[29]*constants[194]*algebraic[69]*states[27])/algebraic[67]+(constants[29]*constants[193]*states[28])/algebraic[67]+algebraic[68]*states[26]*algebraic[118]+(constants[29]*constants[194]*states[26]*states[27])/(algebraic[67]*constants[189])+(algebraic[68]*constants[191]*algebraic[118]*states[28])/constants[190]+(constants[29]*states[27]*states[28])/algebraic[67])
algebraic[70] = (constants[342]*algebraic[44])/algebraic[43]
algebraic[71] = (constants[30]*constants[199])/(algebraic[70]*constants[198])
algebraic[30] = (states[11]*(1.00000+states[1]/constants[313]))/algebraic[29]
algebraic[32] = (states[12]*(1.00000+states[1]/constants[337]))/algebraic[31]
algebraic[72] = 1.00000+states[15]/constants[200]+algebraic[18]/constants[201]+algebraic[23]/constants[202]+algebraic[30]/constants[203]+algebraic[32]/constants[204]
algebraic[73] = (algebraic[71]*constants[30]*(states[28]-states[29]/algebraic[70]))/(constants[198]*algebraic[71]*algebraic[72]+algebraic[71]*states[28]+(constants[30]*states[29])/algebraic[70])
algebraic[132] = states[5]
algebraic[100] = 1.00000*exp(-(constants[287]+4.00000*constants[2]*states[37])/constants[274])
algebraic[138] = (algebraic[100]*(power(states[1], 5.00000)))/(power(algebraic[132], 4.00000))
algebraic[139] = constants[41]*(algebraic[138]*states[21]*algebraic[118]-algebraic[113]*states[27])
algebraic[102] = 1.00000*exp(-(constants[288]+2.00000*constants[2]*states[37])/constants[274])
algebraic[104] = (algebraic[102]*(power(states[1], 2.00000)))/(power(states[1], 4.00000))
algebraic[106] = custom_piecewise([greater(states[27] , constants[247]), states[27] , True, constants[247]])
algebraic[120] = custom_piecewise([greater(algebraic[118] , constants[247]), algebraic[118] , True, constants[247]])
algebraic[121] = constants[16]-states[55]
algebraic[122] = ((constants[42]*(1.00000+states[13]/constants[43]))/(1.00000+states[13]/constants[44]))*((power(algebraic[104], 1.0/2))*algebraic[121]*(power(algebraic[106], 1.0/2))-states[55]*(power(algebraic[120], 1.0/2)))
algebraic[108] = exp(-(constants[289]+4.00000*constants[2]*states[37])/constants[274])
algebraic[140] = (algebraic[108]*(power(states[1], 4.00000)))/((power(algebraic[132], 2.00000))*(power(1.00000, 2.00000)))
algebraic[110] = custom_piecewise([greater(states[56] , constants[248]), states[56] , True, constants[248]])
algebraic[141] = constants[45]*((algebraic[110]/(algebraic[110]+constants[13]))*exp((constants[2]*states[37])/constants[274])*(states[55]/constants[16]))*((power(algebraic[140], 1.00000/2.00000))*states[55]*(power(algebraic[110]/1.00000, 1.00000/4.00000))-algebraic[121])
algebraic[149] = states[6]
algebraic[152] = 1.00000+algebraic[149]/constants[374]
algebraic[155] = constants[34]*((states[32]/algebraic[152])*algebraic[132]-(states[17]/algebraic[46])*states[1])
algebraic[153] = 1.00000+algebraic[132]/constants[376]+algebraic[149]/constants[377]
algebraic[156] = constants[35]*((states[33]/algebraic[153])*algebraic[132]-(states[31]/algebraic[48])*states[1])
algebraic[154] = 1.00000+algebraic[132]/constants[378]+algebraic[149]/constants[379]
algebraic[157] = (exp((-constants[2]*states[37])/constants[274])*algebraic[49]*algebraic[153])/(algebraic[154]*algebraic[48])
algebraic[158] = (constants[40]*(algebraic[157]*states[36]*states[31]*states[1]-states[30]*states[33]*algebraic[132]))/(algebraic[157]*constants[240]*constants[243]*constants[285]*(2.00000*constants[244]+(constants[244]*states[36])/constants[240]+(states[36]*states[31]*states[1])/(constants[240]*constants[243]*constants[285])+(constants[244]*states[30]*algebraic[132])/(constants[241]*constants[285])+(states[30]*states[33]*algebraic[132])/(constants[241]*constants[242]*constants[285])+(constants[244]*states[30])/constants[241]+(constants[244]*states[36]*states[1])/(constants[240]*constants[285])+(constants[244]*states[1])/constants[285]+(constants[244]*states[33]*algebraic[132])/(constants[242]*constants[285])+(constants[244]*algebraic[132])/constants[285]+(constants[244]*states[31]*states[1])/(constants[243]*constants[285])))
algebraic[39] = ((states[15]*states[1])/constants[360])/algebraic[37]
algebraic[162] = states[7]
algebraic[170] = 1.00000+algebraic[132]/constants[360]+algebraic[149]/constants[361]+algebraic[162]/constants[362]
algebraic[172] = ((states[16]*algebraic[132])/constants[360])/algebraic[170]
algebraic[174] = 1.00000+algebraic[132]/constants[370]+algebraic[149]/constants[371]+algebraic[162]/constants[372]
algebraic[175] = constants[36]*((algebraic[172]*states[29])/algebraic[44]-(algebraic[39]*states[34])/algebraic[174])
algebraic[173] = 1.00000+algebraic[132]/constants[366]+algebraic[149]/constants[367]+algebraic[162]/constants[368]
algebraic[177] = constants[38]*(((states[38]/algebraic[173])*states[29])/algebraic[44]-((states[26]/algebraic[42])*states[34])/algebraic[174])
algebraic[168] = 1.00000+algebraic[132]/constants[351]+algebraic[149]/constants[353]+algebraic[162]/constants[355]
algebraic[176] = constants[39]*(((states[34]/algebraic[174])*states[13])/algebraic[33]-((states[29]/algebraic[44])*states[14])/algebraic[168])
algebraic[112] = exp(-(constants[290]-constants[7]*constants[2]*states[37])/constants[274])/1.00000
algebraic[142] = (((algebraic[112]*(power(algebraic[132], constants[7])))/(power(states[1], constants[7]-1.00000)))*algebraic[0])/(algebraic[21]*algebraic[33])
algebraic[143] = constants[46]*(algebraic[142]*states[8]*states[13]-states[0])
algebraic[163] = 1.00000+algebraic[132]/constants[296]+algebraic[149]/constants[311]+algebraic[162]/constants[325]
algebraic[165] = 1.00000+algebraic[132]/constants[332]+algebraic[149]/constants[336]+algebraic[162]/constants[340]
algebraic[178] = custom_piecewise([greater(states[50] , constants[249]) & greater(states[0] , constants[249]), constants[47]*((states[50]/algebraic[165])/(states[50]/algebraic[165]+(states[51]/algebraic[163])*exp((constants[12]*constants[2]*states[37])/constants[274]))-(states[8]/algebraic[21])/(states[8]/algebraic[21]+(states[0]/algebraic[0])*exp((constants[12]*(constants[2]-1.00000)*states[37])/constants[274])))*((states[50]/algebraic[165])/(states[50]/algebraic[165]+constants[11])) , True, 0.00000])
algebraic[144] = custom_piecewise([greater(fabs(states[37]) , constants[250]), (constants[51]*states[37]*(algebraic[132]*exp((constants[2]*states[37])/constants[274])-states[1]))/(exp((constants[2]*states[37])/constants[274])-1.00000) , True, (constants[51]*constants[274]*(algebraic[132]-states[1]))/constants[2]])
algebraic[34] = ((states[13]*states[1])/constants[351])/algebraic[33]
algebraic[169] = ((states[14]*algebraic[132])/constants[351])/algebraic[168]
algebraic[179] = (constants[48]*(algebraic[169]*algebraic[132]-algebraic[34]*states[1]))/(constants[49]*(1.00000+algebraic[169]/constants[49])*(1.00000+algebraic[34]/constants[49]))
algebraic[167] = 1.00000+algebraic[132]/constants[344]+algebraic[149]/constants[347]+algebraic[162]/constants[349]
algebraic[180] = (constants[381]*algebraic[163]*algebraic[167])/(algebraic[165]*algebraic[165])
algebraic[182] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[261]*(algebraic[180]*states[50]*states[50]-states[53]*states[51])])
algebraic[189] = -((states[1]*states[1])/(constants[134]*constants[135]*(power(constants[282], 2.00000)))+(states[1]*states[0])/(constants[296]*constants[311]*(power(algebraic[0], 2.00000)))+(states[1]*states[8])/(constants[332]*constants[336]*(power(algebraic[21], 2.00000)))+(states[1]*states[10])/(constants[344]*constants[347]*(power(algebraic[26], 2.00000)))+(states[1]*states[11])/(constants[296]*constants[311]*(power(algebraic[29], 2.00000)))+(states[1]*states[12])/(constants[332]*constants[336]*(power(algebraic[31], 2.00000)))+(states[1]*states[13])/(constants[351]*constants[353]*(power(algebraic[33], 2.00000)))+(states[1]*states[21])/(constants[110]*constants[111]*(power(constants[277], 2.00000)))+(states[1]*algebraic[113])/(constants[107]*constants[108]*(power(constants[276], 2.00000)))+(states[1]*states[27])/(constants[113]*constants[114]*(power(constants[304], 2.00000)))+(states[1]*algebraic[118])/(constants[116]*constants[117]*(power(constants[320], 2.00000)))+(states[1]*states[17])/(constants[98]*constants[374]*(power(algebraic[46], 2.00000)))+(states[1]*states[22])/(constants[89]*constants[359]*(power(algebraic[36], 2.00000)))+(states[1]*states[20])/(constants[119]*constants[120]*(power(constants[278], 2.00000)))+(states[1]*states[15])/(constants[360]*constants[361]*(power(algebraic[37], 2.00000)))+(states[1]*states[24])/(constants[363]*constants[364]*(power(algebraic[40], 2.00000)))+(states[1]*states[25])/(constants[122]*constants[123]*(power(constants[306], 2.00000)))+(states[1]*states[23])/(constants[365]*constants[92]*(power(algebraic[41], 2.00000)))+(states[1]*states[18])/(constants[357]*constants[87]*(power(algebraic[35], 2.00000)))+(states[1]*states[26])/(constants[366]*constants[367]*(power(algebraic[42], 2.00000)))+(states[1]*states[28])/(constants[369]*constants[94]*(power(algebraic[43], 2.00000)))+(states[1]*states[29])/(constants[370]*constants[371]*(power(algebraic[44], 2.00000)))+(states[1]*states[31])/(constants[376]*constants[377]*(power(algebraic[48], 2.00000)))+(states[1]*states[30])/(constants[378]*constants[379]*(power(algebraic[49], 2.00000)))+(states[1]*states[3])/(constants[140]*constants[141]*(power(constants[284], 2.00000)))+(states[1]*states[2])/(constants[137]*constants[138]*(power(constants[283], 2.00000)))+(states[1]*states[56])/(constants[125]*constants[126]*(power(constants[279], 2.00000)))+(states[1]*constants[291])/(constants[128]*constants[129]*(power(constants[280], 2.00000)))+(states[1]*constants[310])/(constants[131]*constants[132]*(power(constants[281], 2.00000)))+(states[1]*states[19])/(constants[373]*constants[96]*(power(algebraic[45], 2.00000))))
algebraic[190] = -((states[1]*states[1])/(constants[134]*constants[136]*(power(constants[282], 2.00000)))+(states[1]*states[0])/(constants[296]*constants[325]*(power(algebraic[0], 2.00000)))+(states[1]*states[8])/(constants[332]*constants[340]*(power(algebraic[21], 2.00000)))+(states[1]*states[10])/(constants[344]*constants[349]*(power(algebraic[26], 2.00000)))+(states[1]*states[11])/(constants[296]*constants[325]*(power(algebraic[29], 2.00000)))+(states[1]*states[12])/(constants[332]*constants[340]*(power(algebraic[31], 2.00000)))+(states[1]*states[13])/(constants[351]*constants[355]*(power(algebraic[33], 2.00000)))+(states[1]*states[21])/(constants[110]*constants[112]*(power(constants[277], 2.00000)))+(states[1]*algebraic[113])/(constants[107]*constants[109]*(power(constants[276], 2.00000)))+(states[1]*states[27])/(constants[113]*constants[115]*(power(constants[304], 2.00000)))+(states[1]*algebraic[118])/(constants[116]*constants[118]*(power(constants[320], 2.00000)))+(states[1]*states[17])/(constants[98]*constants[99]*(power(algebraic[46], 2.00000)))+(states[1]*states[22])/(constants[89]*constants[90]*(power(algebraic[36], 2.00000)))+(states[1]*states[20])/(constants[119]*constants[121]*(power(constants[278], 2.00000)))+(states[1]*states[15])/(constants[360]*constants[362]*(power(algebraic[37], 2.00000)))+(states[1]*states[24])/(constants[363]*constants[91]*(power(algebraic[40], 2.00000)))+(states[1]*states[25])/(constants[122]*constants[124]*(power(constants[306], 2.00000)))+(states[1]*states[23])/(constants[365]*constants[93]*(power(algebraic[41], 2.00000)))+(states[1]*states[18])/(constants[357]*constants[88]*(power(algebraic[35], 2.00000)))+(states[1]*states[26])/(constants[366]*constants[368]*(power(algebraic[42], 2.00000)))+(states[1]*states[28])/(constants[369]*constants[95]*(power(algebraic[43], 2.00000)))+(states[1]*states[29])/(constants[370]*constants[372]*(power(algebraic[44], 2.00000)))+(states[1]*states[31])/(constants[376]*constants[102]*(power(algebraic[48], 2.00000)))+(states[1]*states[30])/(constants[378]*constants[103]*(power(algebraic[49], 2.00000)))+(states[1]*states[3])/(constants[140]*constants[142]*(power(constants[284], 2.00000)))+(states[1]*states[2])/(constants[137]*constants[139]*(power(constants[283], 2.00000)))+(states[1]*states[56])/(constants[125]*constants[127]*(power(constants[279], 2.00000)))+(states[1]*constants[291])/(constants[128]*constants[130]*(power(constants[280], 2.00000)))+(states[1]*constants[310])/(constants[131]*constants[133]*(power(constants[281], 2.00000)))+(states[1]*states[19])/(constants[373]*constants[97]*(power(algebraic[45], 2.00000))))
algebraic[193] = -((states[2]*states[1])/(constants[135]*constants[136]*(power(constants[282], 2.00000)))+(states[2]*states[0])/(constants[311]*constants[325]*(power(algebraic[0], 2.00000)))+(states[2]*states[8])/(constants[336]*constants[340]*(power(algebraic[21], 2.00000)))+(states[2]*states[10])/(constants[347]*constants[349]*(power(algebraic[26], 2.00000)))+(states[2]*states[11])/(constants[311]*constants[325]*(power(algebraic[29], 2.00000)))+(states[2]*states[12])/(constants[336]*constants[340]*(power(algebraic[31], 2.00000)))+(states[2]*states[13])/(constants[353]*constants[355]*(power(algebraic[33], 2.00000)))+(states[2]*states[21])/(constants[111]*constants[112]*(power(constants[277], 2.00000)))+(states[2]*algebraic[113])/(constants[108]*constants[109]*(power(constants[276], 2.00000)))+(states[2]*states[27])/(constants[114]*constants[115]*(power(constants[304], 2.00000)))+(states[2]*algebraic[118])/(constants[117]*constants[118]*(power(constants[320], 2.00000)))+(states[2]*states[17])/(constants[374]*constants[99]*(power(algebraic[46], 2.00000)))+(states[2]*states[22])/(constants[359]*constants[90]*(power(algebraic[36], 2.00000)))+(states[2]*states[20])/(constants[120]*constants[121]*(power(constants[278], 2.00000)))+(states[2]*states[15])/(constants[361]*constants[362]*(power(algebraic[37], 2.00000)))+(states[2]*states[24])/(constants[364]*constants[91]*(power(algebraic[40], 2.00000)))+(states[2]*states[25])/(constants[123]*constants[124]*(power(constants[306], 2.00000)))+(states[2]*states[23])/(constants[92]*constants[93]*(power(algebraic[41], 2.00000)))+(states[2]*states[18])/(constants[87]*constants[88]*(power(algebraic[35], 2.00000)))+(states[2]*states[26])/(constants[367]*constants[368]*(power(algebraic[42], 2.00000)))+(states[2]*states[28])/(constants[94]*constants[95]*(power(algebraic[43], 2.00000)))+(states[2]*states[29])/(constants[371]*constants[372]*(power(algebraic[44], 2.00000)))+(states[2]*states[31])/(constants[377]*constants[102]*(power(algebraic[48], 2.00000)))+(states[2]*states[30])/(constants[379]*constants[103]*(power(algebraic[49], 2.00000)))+(states[2]*states[3])/(constants[141]*constants[142]*(power(constants[284], 2.00000)))+(states[2]*states[2])/(constants[138]*constants[139]*(power(constants[283], 2.00000)))+(states[2]*states[56])/(constants[126]*constants[127]*(power(constants[279], 2.00000)))+(states[2]*constants[291])/(constants[129]*constants[130]*(power(constants[280], 2.00000)))+(states[2]*constants[310])/(constants[132]*constants[133]*(power(constants[281], 2.00000)))+(states[2]*states[19])/(constants[96]*constants[97]*(power(algebraic[45], 2.00000))))
algebraic[196] = -((states[3]*states[1])/(constants[136]*constants[135]*(power(constants[282], 2.00000)))+(states[3]*states[0])/(constants[325]*constants[311]*(power(algebraic[0], 2.00000)))+(states[3]*states[8])/(constants[340]*constants[336]*(power(algebraic[21], 2.00000)))+(states[3]*states[10])/(constants[349]*constants[347]*(power(algebraic[26], 2.00000)))+(states[3]*states[11])/(constants[325]*constants[311]*(power(algebraic[29], 2.00000)))+(states[3]*states[12])/(constants[340]*constants[336]*(power(algebraic[31], 2.00000)))+(states[3]*states[13])/(constants[355]*constants[353]*(power(algebraic[33], 2.00000)))+(states[3]*states[21])/(constants[112]*constants[111]*(power(constants[277], 2.00000)))+(states[3]*algebraic[113])/(constants[109]*constants[108]*(power(constants[276], 2.00000)))+(states[3]*states[27])/(constants[115]*constants[114]*(power(constants[304], 2.00000)))+(states[3]*algebraic[118])/(constants[118]*constants[117]*(power(constants[320], 2.00000)))+(states[3]*states[17])/(constants[99]*constants[374]*(power(algebraic[46], 2.00000)))+(states[3]*states[22])/(constants[90]*constants[359]*(power(algebraic[36], 2.00000)))+(states[3]*states[20])/(constants[121]*constants[120]*(power(constants[278], 2.00000)))+(states[3]*states[15])/(constants[362]*constants[361]*(power(algebraic[37], 2.00000)))+(states[3]*states[24])/(constants[91]*constants[364]*(power(algebraic[40], 2.00000)))+(states[3]*states[25])/(constants[124]*constants[123]*(power(constants[306], 2.00000)))+(states[3]*states[23])/(constants[93]*constants[92]*(power(algebraic[41], 2.00000)))+(states[3]*states[18])/(constants[88]*constants[87]*(power(algebraic[35], 2.00000)))+(states[3]*states[26])/(constants[368]*constants[367]*(power(algebraic[42], 2.00000)))+(states[3]*states[28])/(constants[95]*constants[94]*(power(algebraic[43], 2.00000)))+(states[3]*states[29])/(constants[372]*constants[371]*(power(algebraic[44], 2.00000)))+(states[3]*states[31])/(constants[102]*constants[377]*(power(algebraic[48], 2.00000)))+(states[3]*states[30])/(constants[103]*constants[379]*(power(algebraic[49], 2.00000)))+(states[3]*states[3])/(constants[142]*constants[141]*(power(constants[284], 2.00000)))+(states[3]*states[2])/(constants[139]*constants[138]*(power(constants[283], 2.00000)))+(states[3]*states[56])/(constants[127]*constants[126]*(power(constants[279], 2.00000)))+(states[3]*constants[291])/(constants[130]*constants[129]*(power(constants[280], 2.00000)))+(states[3]*constants[310])/(constants[133]*constants[132]*(power(constants[281], 2.00000)))+(states[3]*states[19])/(constants[97]*constants[96]*(power(algebraic[45], 2.00000))))
algebraic[197] = ((1.00000+states[1]/constants[134]+states[2]/constants[135])*states[1])/(constants[136]*(power(constants[282], 2.00000)))+((1.00000+states[1]/constants[296]+states[2]/constants[311])*states[0])/(constants[325]*(power(algebraic[0], 2.00000)))+((1.00000+states[1]/constants[332]+states[2]/constants[336])*states[8])/(constants[340]*(power(algebraic[21], 2.00000)))+((1.00000+states[1]/constants[344]+states[2]/constants[347])*states[10])/(constants[349]*(power(algebraic[26], 2.00000)))+((1.00000+states[1]/constants[296]+states[2]/constants[311])*states[11])/(constants[325]*(power(algebraic[29], 2.00000)))+((1.00000+states[1]/constants[332]+states[2]/constants[336])*states[12])/(constants[340]*(power(algebraic[31], 2.00000)))+((1.00000+states[1]/constants[351]+states[2]/constants[353])*states[13])/(constants[355]*(power(algebraic[33], 2.00000)))+((1.00000+states[1]/constants[110]+states[2]/constants[111])*states[21])/(constants[112]*(power(constants[277], 2.00000)))+((1.00000+states[1]/constants[107]+states[2]/constants[108])*algebraic[113])/(constants[109]*(power(constants[276], 2.00000)))+((1.00000+states[1]/constants[113]+states[2]/constants[114])*states[27])/(constants[115]*(power(constants[304], 2.00000)))+((1.00000+states[1]/constants[116]+states[2]/constants[117])*algebraic[118])/(constants[118]*(power(constants[320], 2.00000)))+((1.00000+states[1]/constants[98]+states[2]/constants[374])*states[17])/(constants[99]*(power(algebraic[46], 2.00000)))+((1.00000+states[1]/constants[89]+states[2]/constants[359])*states[22])/(constants[90]*(power(algebraic[36], 2.00000)))+((1.00000+states[1]/constants[119]+states[2]/constants[120])*states[20])/(constants[121]*(power(constants[278], 2.00000)))+((1.00000+states[1]/constants[360]+states[2]/constants[361])*states[15])/(constants[362]*(power(algebraic[37], 2.00000)))+((1.00000+states[1]/constants[363]+states[2]/constants[364])*states[24])/(constants[91]*(power(algebraic[40], 2.00000)))+((1.00000+states[1]/constants[122]+states[2]/constants[123])*states[25])/(constants[124]*(power(constants[306], 2.00000)))+((1.00000+states[1]/constants[365]+states[2]/constants[92])*states[23])/(constants[93]*(power(algebraic[41], 2.00000)))+((1.00000+states[1]/constants[357]+states[2]/constants[87])*states[18])/(constants[88]*(power(algebraic[35], 2.00000)))+((1.00000+states[1]/constants[366]+states[2]/constants[367])*states[26])/(constants[368]*(power(algebraic[42], 2.00000)))+((1.00000+states[1]/constants[369]+states[2]/constants[94])*states[28])/(constants[95]*(power(algebraic[43], 2.00000)))+((1.00000+states[1]/constants[370]+states[2]/constants[371])*states[29])/(constants[372]*(power(algebraic[44], 2.00000)))+((1.00000+states[1]/constants[376]+states[2]/constants[377])*states[31])/(constants[102]*(power(algebraic[48], 2.00000)))+((1.00000+states[1]/constants[378]+states[2]/constants[379])*states[30])/(constants[103]*(power(algebraic[49], 2.00000)))+((1.00000+states[1]/constants[140]+states[2]/constants[141])*states[3])/(constants[142]*(power(constants[284], 2.00000)))+((1.00000+states[1]/constants[137]+states[2]/constants[138])*states[2])/(constants[139]*(power(constants[283], 2.00000)))+((1.00000+states[1]/constants[125]+states[2]/constants[126])*states[56])/(constants[127]*(power(constants[279], 2.00000)))+((1.00000+states[1]/constants[128]+states[2]/constants[129])*constants[291])/(constants[130]*(power(constants[280], 2.00000)))+((1.00000+states[1]/constants[131]+states[2]/constants[132])*constants[310])/(constants[133]*(power(constants[281], 2.00000)))+((1.00000+states[1]/constants[373]+states[2]/constants[96])*states[19])/(constants[97]*(power(algebraic[45], 2.00000)))
algebraic[205] = 1.00000+algebraic[197]
algebraic[194] = ((1.00000+states[1]/constants[134]+states[3]/constants[136])*states[1])/(constants[135]*(power(constants[282], 2.00000)))+((1.00000+states[1]/constants[296]+states[3]/constants[325])*states[0])/(constants[311]*(power(algebraic[0], 2.00000)))+((1.00000+states[1]/constants[332]+states[3]/constants[340])*states[8])/(constants[336]*(power(algebraic[21], 2.00000)))+((1.00000+states[1]/constants[344]+states[3]/constants[349])*states[10])/(constants[347]*(power(algebraic[26], 2.00000)))+((1.00000+states[1]/constants[296]+states[3]/constants[325])*states[11])/(constants[311]*(power(algebraic[29], 2.00000)))+((1.00000+states[1]/constants[332]+states[3]/constants[340])*states[12])/(constants[336]*(power(algebraic[31], 2.00000)))+((1.00000+states[1]/constants[351]+states[3]/constants[355])*states[13])/(constants[353]*(power(algebraic[33], 2.00000)))+((1.00000+states[1]/constants[110]+states[3]/constants[112])*states[21])/(constants[111]*(power(constants[277], 2.00000)))+((1.00000+states[1]/constants[107]+states[3]/constants[109])*algebraic[113])/(constants[108]*(power(constants[276], 2.00000)))+((1.00000+states[1]/constants[113]+states[3]/constants[115])*states[27])/(constants[114]*(power(constants[304], 2.00000)))+((1.00000+states[1]/constants[116]+states[3]/constants[118])*algebraic[118])/(constants[117]*(power(constants[320], 2.00000)))+((1.00000+states[1]/constants[98]+states[3]/constants[99])*states[17])/(constants[374]*(power(algebraic[46], 2.00000)))+((1.00000+states[1]/constants[89]+states[3]/constants[90])*states[22])/(constants[359]*(power(algebraic[36], 2.00000)))+((1.00000+states[1]/constants[119]+states[3]/constants[121])*states[20])/(constants[120]*(power(constants[278], 2.00000)))+((1.00000+states[1]/constants[360]+states[3]/constants[362])*states[15])/(constants[361]*(power(algebraic[37], 2.00000)))+((1.00000+states[1]/constants[363]+states[3]/constants[91])*states[24])/(constants[364]*(power(algebraic[40], 2.00000)))+((1.00000+states[1]/constants[122]+states[3]/constants[124])*states[25])/(constants[123]*(power(constants[306], 2.00000)))+((1.00000+states[1]/constants[365]+states[3]/constants[93])*states[23])/(constants[92]*(power(algebraic[41], 2.00000)))+((1.00000+states[1]/constants[357]+states[3]/constants[88])*states[18])/(constants[87]*(power(algebraic[35], 2.00000)))+((1.00000+states[1]/constants[366]+states[3]/constants[368])*states[26])/(constants[367]*(power(algebraic[42], 2.00000)))+((1.00000+states[1]/constants[369]+states[3]/constants[95])*states[28])/(constants[94]*(power(algebraic[43], 2.00000)))+((1.00000+states[1]/constants[370]+states[3]/constants[372])*states[29])/(constants[371]*(power(algebraic[44], 2.00000)))+((1.00000+states[1]/constants[376]+states[3]/constants[102])*states[31])/(constants[377]*(power(algebraic[48], 2.00000)))+((1.00000+states[1]/constants[378]+states[3]/constants[103])*states[30])/(constants[379]*(power(algebraic[49], 2.00000)))+((1.00000+states[1]/constants[140]+states[3]/constants[142])*states[3])/(constants[141]*(power(constants[284], 2.00000)))+((1.00000+states[1]/constants[137]+states[3]/constants[139])*states[2])/(constants[138]*(power(constants[283], 2.00000)))+((1.00000+states[1]/constants[125]+states[3]/constants[127])*states[56])/(constants[126]*(power(constants[279], 2.00000)))+((1.00000+states[1]/constants[128]+states[3]/constants[130])*constants[291])/(constants[129]*(power(constants[280], 2.00000)))+((1.00000+states[1]/constants[131]+states[3]/constants[133])*constants[310])/(constants[132]*(power(constants[281], 2.00000)))+((1.00000+states[1]/constants[373]+states[3]/constants[97])*states[19])/(constants[96]*(power(algebraic[45], 2.00000)))
algebraic[204] = 1.00000+algebraic[194]
algebraic[191] = ((1.00000+states[2]/constants[135]+states[3]/constants[136])*states[1])/(constants[134]*(power(constants[282], 2.00000)))+((1.00000+states[2]/constants[311]+states[3]/constants[325])*states[0])/(constants[296]*(power(algebraic[0], 2.00000)))+((1.00000+states[2]/constants[336]+states[3]/constants[340])*states[8])/(constants[332]*(power(algebraic[21], 2.00000)))+((1.00000+states[2]/constants[347]+states[3]/constants[349])*states[10])/(constants[344]*(power(algebraic[26], 2.00000)))+((1.00000+states[2]/constants[311]+states[3]/constants[325])*states[11])/(constants[296]*(power(algebraic[29], 2.00000)))+((1.00000+states[2]/constants[336]+states[3]/constants[340])*states[12])/(constants[332]*(power(algebraic[31], 2.00000)))+((1.00000+states[2]/constants[353]+states[3]/constants[355])*states[13])/(constants[351]*(power(algebraic[33], 2.00000)))+((1.00000+states[2]/constants[111]+states[3]/constants[112])*states[21])/(constants[110]*(power(constants[277], 2.00000)))+((1.00000+states[2]/constants[108]+states[3]/constants[109])*algebraic[113])/(constants[107]*(power(constants[276], 2.00000)))+((1.00000+states[2]/constants[114]+states[3]/constants[115])*states[27])/(constants[113]*(power(constants[304], 2.00000)))+((1.00000+states[2]/constants[117]+states[3]/constants[118])*algebraic[118])/(constants[116]*(power(constants[320], 2.00000)))+((1.00000+states[2]/constants[374]+states[3]/constants[99])*states[17])/(constants[98]*(power(algebraic[46], 2.00000)))+((1.00000+states[2]/constants[359]+states[3]/constants[90])*states[22])/(constants[89]*(power(algebraic[36], 2.00000)))+((1.00000+states[2]/constants[120]+states[3]/constants[121])*states[20])/(constants[119]*(power(constants[278], 2.00000)))+((1.00000+states[2]/constants[361]+states[3]/constants[362])*states[15])/(constants[360]*(power(algebraic[37], 2.00000)))+((1.00000+states[2]/constants[364]+states[3]/constants[91])*states[24])/(constants[363]*(power(algebraic[40], 2.00000)))+((1.00000+states[2]/constants[123]+states[3]/constants[124])*states[25])/(constants[122]*(power(constants[306], 2.00000)))+((1.00000+states[2]/constants[92]+states[3]/constants[93])*states[23])/(constants[365]*(power(algebraic[41], 2.00000)))+((1.00000+states[2]/constants[87]+states[3]/constants[88])*states[18])/(constants[357]*(power(algebraic[35], 2.00000)))+((1.00000+states[2]/constants[367]+states[3]/constants[368])*states[26])/(constants[366]*(power(algebraic[42], 2.00000)))+((1.00000+states[2]/constants[94]+states[3]/constants[95])*states[28])/(constants[369]*(power(algebraic[43], 2.00000)))+((1.00000+states[2]/constants[371]+states[3]/constants[372])*states[29])/(constants[370]*(power(algebraic[44], 2.00000)))+((1.00000+states[2]/constants[377]+states[3]/constants[102])*states[31])/(constants[376]*(power(algebraic[48], 2.00000)))+((1.00000+states[2]/constants[379]+states[3]/constants[103])*states[30])/(constants[378]*(power(algebraic[49], 2.00000)))+((1.00000+states[2]/constants[141]+states[3]/constants[142])*states[3])/(constants[140]*(power(constants[284], 2.00000)))+((1.00000+states[2]/constants[138]+states[3]/constants[139])*states[2])/(constants[137]*(power(constants[283], 2.00000)))+((1.00000+states[2]/constants[126]+states[3]/constants[127])*states[56])/(constants[125]*(power(constants[279], 2.00000)))+((1.00000+states[2]/constants[129]+states[3]/constants[130])*constants[291])/(constants[128]*(power(constants[280], 2.00000)))+((1.00000+states[2]/constants[132]+states[3]/constants[133])*constants[310])/(constants[131]*(power(constants[281], 2.00000)))+((1.00000+states[2]/constants[96]+states[3]/constants[97])*states[19])/(constants[373]*(power(algebraic[45], 2.00000)))
algebraic[203] = 1.00000+algebraic[191]+constants[271]/(constants[272]*(power(1.00000+states[1]/constants[272], 2.00000)))
algebraic[195] = -((states[3]*states[1])/(constants[136]*constants[134]*(power(constants[282], 2.00000)))+(states[3]*states[0])/(constants[325]*constants[296]*(power(algebraic[0], 2.00000)))+(states[3]*states[8])/(constants[340]*constants[332]*(power(algebraic[21], 2.00000)))+(states[3]*states[10])/(constants[349]*constants[344]*(power(algebraic[26], 2.00000)))+(states[3]*states[11])/(constants[325]*constants[296]*(power(algebraic[29], 2.00000)))+(states[3]*states[12])/(constants[340]*constants[332]*(power(algebraic[31], 2.00000)))+(states[3]*states[13])/(constants[355]*constants[351]*(power(algebraic[33], 2.00000)))+(states[3]*states[21])/(constants[112]*constants[110]*(power(constants[277], 2.00000)))+(states[3]*algebraic[113])/(constants[109]*constants[107]*(power(constants[276], 2.00000)))+(states[3]*states[27])/(constants[115]*constants[113]*(power(constants[304], 2.00000)))+(states[3]*algebraic[118])/(constants[118]*constants[116]*(power(constants[320], 2.00000)))+(states[3]*states[17])/(constants[99]*constants[98]*(power(algebraic[46], 2.00000)))+(states[3]*states[22])/(constants[90]*constants[89]*(power(algebraic[36], 2.00000)))+(states[3]*states[20])/(constants[121]*constants[119]*(power(constants[278], 2.00000)))+(states[3]*states[15])/(constants[362]*constants[360]*(power(algebraic[37], 2.00000)))+(states[3]*states[24])/(constants[91]*constants[363]*(power(algebraic[40], 2.00000)))+(states[3]*states[25])/(constants[124]*constants[122]*(power(constants[306], 2.00000)))+(states[3]*states[23])/(constants[93]*constants[365]*(power(algebraic[41], 2.00000)))+(states[3]*states[18])/(constants[88]*constants[357]*(power(algebraic[35], 2.00000)))+(states[3]*states[26])/(constants[368]*constants[366]*(power(algebraic[42], 2.00000)))+(states[3]*states[28])/(constants[95]*constants[369]*(power(algebraic[43], 2.00000)))+(states[3]*states[29])/(constants[372]*constants[370]*(power(algebraic[44], 2.00000)))+(states[3]*states[31])/(constants[102]*constants[376]*(power(algebraic[48], 2.00000)))+(states[3]*states[30])/(constants[103]*constants[378]*(power(algebraic[49], 2.00000)))+(states[3]*states[3])/(constants[142]*constants[140]*(power(constants[284], 2.00000)))+(states[3]*states[2])/(constants[139]*constants[137]*(power(constants[283], 2.00000)))+(states[3]*states[56])/(constants[127]*constants[125]*(power(constants[279], 2.00000)))+(states[3]*constants[291])/(constants[130]*constants[128]*(power(constants[280], 2.00000)))+(states[3]*constants[310])/(constants[133]*constants[131]*(power(constants[281], 2.00000)))+(states[3]*states[19])/(constants[97]*constants[373]*(power(algebraic[45], 2.00000))))
algebraic[192] = -((states[2]*states[1])/(constants[135]*constants[134]*(power(constants[282], 2.00000)))+(states[2]*states[0])/(constants[311]*constants[296]*(power(algebraic[0], 2.00000)))+(states[2]*states[8])/(constants[336]*constants[332]*(power(algebraic[21], 2.00000)))+(states[2]*states[10])/(constants[347]*constants[344]*(power(algebraic[26], 2.00000)))+(states[2]*states[11])/(constants[311]*constants[296]*(power(algebraic[29], 2.00000)))+(states[2]*states[12])/(constants[336]*constants[332]*(power(algebraic[31], 2.00000)))+(states[2]*states[13])/(constants[353]*constants[351]*(power(algebraic[33], 2.00000)))+(states[2]*states[21])/(constants[111]*constants[110]*(power(constants[277], 2.00000)))+(states[2]*algebraic[113])/(constants[108]*constants[107]*(power(constants[276], 2.00000)))+(states[2]*states[27])/(constants[114]*constants[113]*(power(constants[304], 2.00000)))+(states[2]*algebraic[118])/(constants[117]*constants[116]*(power(constants[320], 2.00000)))+(states[2]*states[17])/(constants[374]*constants[98]*(power(algebraic[46], 2.00000)))+(states[2]*states[22])/(constants[359]*constants[89]*(power(algebraic[36], 2.00000)))+(states[2]*states[20])/(constants[120]*constants[119]*(power(constants[278], 2.00000)))+(states[2]*states[15])/(constants[361]*constants[360]*(power(algebraic[37], 2.00000)))+(states[2]*states[24])/(constants[364]*constants[363]*(power(algebraic[40], 2.00000)))+(states[2]*states[25])/(constants[123]*constants[122]*(power(constants[306], 2.00000)))+(states[2]*states[23])/(constants[92]*constants[365]*(power(algebraic[41], 2.00000)))+(states[2]*states[18])/(constants[87]*constants[357]*(power(algebraic[35], 2.00000)))+(states[2]*states[26])/(constants[367]*constants[366]*(power(algebraic[42], 2.00000)))+(states[2]*states[28])/(constants[94]*constants[369]*(power(algebraic[43], 2.00000)))+(states[2]*states[29])/(constants[371]*constants[370]*(power(algebraic[44], 2.00000)))+(states[2]*states[31])/(constants[377]*constants[376]*(power(algebraic[48], 2.00000)))+(states[2]*states[30])/(constants[379]*constants[378]*(power(algebraic[49], 2.00000)))+(states[2]*states[3])/(constants[141]*constants[140]*(power(constants[284], 2.00000)))+(states[2]*states[2])/(constants[138]*constants[137]*(power(constants[283], 2.00000)))+(states[2]*states[56])/(constants[126]*constants[125]*(power(constants[279], 2.00000)))+(states[2]*constants[291])/(constants[129]*constants[128]*(power(constants[280], 2.00000)))+(states[2]*constants[310])/(constants[132]*constants[131]*(power(constants[281], 2.00000)))+(states[2]*states[19])/(constants[96]*constants[373]*(power(algebraic[45], 2.00000))))
algebraic[206] = (((algebraic[203]*algebraic[196]*algebraic[193]+algebraic[205]*algebraic[189]*algebraic[192]+algebraic[204]*algebraic[190]*algebraic[195])-algebraic[204]*algebraic[205]*algebraic[203])-algebraic[190]*algebraic[196]*algebraic[192])-algebraic[189]*algebraic[193]*algebraic[195]
algebraic[145] = ((algebraic[114]+algebraic[115]+algebraic[116]+algebraic[117])-algebraic[139])/constants[3]
algebraic[146] = -algebraic[145]
algebraic[147] = ((algebraic[119]+algebraic[139])-algebraic[122])/constants[3]
algebraic[148] = -algebraic[147]
algebraic[183] = ((algebraic[81]+algebraic[143])-algebraic[178])/constants[3]
algebraic[184] = ((-algebraic[81]-algebraic[143])+algebraic[178])/constants[3]
algebraic[123] = (algebraic[66]-algebraic[81])/constants[3]
algebraic[124] = (-algebraic[66]+algebraic[81])/constants[3]
algebraic[185] = (((-algebraic[66]-algebraic[143])+algebraic[179])-algebraic[176])/constants[3]
algebraic[159] = (-algebraic[114]+algebraic[155])/constants[3]
algebraic[125] = (algebraic[114]-algebraic[56])/constants[3]
algebraic[186] = (-algebraic[59]+algebraic[175]+algebraic[56])/constants[3]
algebraic[126] = (algebraic[59]-algebraic[115])/constants[3]
algebraic[127] = (((algebraic[115]-algebraic[116])-algebraic[86])+algebraic[87])/constants[3]
algebraic[128] = (algebraic[116]-algebraic[66])/constants[3]
algebraic[129] = ((-algebraic[114]-algebraic[116])+algebraic[66]+algebraic[56])/constants[3]
algebraic[187] = ((algebraic[66]-algebraic[119])+algebraic[177])/constants[3]
algebraic[130] = (algebraic[119]-algebraic[73])/constants[3]
algebraic[188] = (((((algebraic[73]-algebraic[117])+algebraic[176])-algebraic[87])-algebraic[175])-algebraic[177])/constants[3]
algebraic[131] = (-algebraic[56]+algebraic[117]+algebraic[86])/constants[3]
algebraic[160] = ((algebraic[86]+algebraic[156])-algebraic[158])/constants[3]
algebraic[161] = (-algebraic[86]+algebraic[158])/constants[3]
algebraic[181] = constants[50]*(algebraic[162]*states[1]-states[3]*algebraic[132])
algebraic[198] = (((((((((-algebraic[114]+2.00000*algebraic[56])-algebraic[116])+algebraic[66]+algebraic[117]+algebraic[155]+algebraic[156]+algebraic[175])-algebraic[158])-5.00000*algebraic[139])-2.00000*algebraic[122])-4.00000*algebraic[141])+(constants[7]-1.00000)*algebraic[143]+2.00000*algebraic[179]+algebraic[144])-algebraic[181])/constants[3]
algebraic[200] = -((states[1]*constants[292])/(constants[134]*constants[282])+(states[1]*algebraic[183])/(constants[296]*algebraic[0])+(states[1]*algebraic[184])/(constants[332]*algebraic[21])+(states[1]*constants[262])/(constants[344]*algebraic[26])+(states[1]*algebraic[123])/(constants[296]*algebraic[29])+(states[1]*algebraic[124])/(constants[332]*algebraic[31])+(states[1]*algebraic[185])/(constants[351]*algebraic[33])+(states[1]*algebraic[145])/(constants[110]*constants[277])+(states[1]*algebraic[146])/(constants[107]*constants[276])+(states[1]*algebraic[147])/(constants[113]*constants[304])+(states[1]*algebraic[148])/(constants[116]*constants[320])+(states[1]*algebraic[159])/(constants[98]*algebraic[46])+(states[1]*algebraic[131])/(constants[89]*algebraic[36])+(states[1]*algebraic[125])/(constants[119]*constants[278])+(states[1]*algebraic[186])/(constants[360]*algebraic[37])+(states[1]*algebraic[126])/(constants[363]*algebraic[40])+(states[1]*algebraic[127])/(constants[122]*constants[306])+(states[1]*algebraic[128])/(constants[365]*algebraic[41])+(states[1]*algebraic[129])/(constants[357]*algebraic[35])+(states[1]*algebraic[187])/(constants[366]*algebraic[42])+(states[1]*algebraic[130])/(constants[369]*algebraic[43])+(states[1]*algebraic[188])/(constants[370]*algebraic[44])+(states[1]*algebraic[160])/(constants[376]*algebraic[48])+(states[1]*algebraic[161])/(constants[378]*algebraic[49])+(states[1]*constants[294])/(constants[140]*constants[284])+(states[1]*constants[293])/(constants[137]*constants[283])+(states[1]*constants[263])/(constants[125]*constants[279])+(states[1]*constants[270])/(constants[128]*constants[280])+(states[1]*constants[269])/(constants[131]*constants[281])+(states[1]*constants[264])/(constants[373]*algebraic[45]))+algebraic[198]
algebraic[201] = -((states[2]*constants[292])/(constants[135]*constants[282])+(states[2]*algebraic[183])/(constants[311]*algebraic[0])+(states[2]*algebraic[184])/(constants[336]*algebraic[21])+(states[2]*constants[262])/(constants[347]*algebraic[26])+(states[2]*algebraic[123])/(constants[311]*algebraic[29])+(states[2]*algebraic[124])/(constants[336]*algebraic[31])+(states[2]*algebraic[185])/(constants[353]*algebraic[33])+(states[2]*algebraic[145])/(constants[111]*constants[277])+(states[2]*algebraic[146])/(constants[108]*constants[276])+(states[2]*algebraic[147])/(constants[114]*constants[304])+(states[2]*algebraic[148])/(constants[117]*constants[320])+(states[2]*algebraic[159])/(constants[374]*algebraic[46])+(states[2]*algebraic[131])/(constants[359]*algebraic[36])+(states[2]*algebraic[125])/(constants[120]*constants[278])+(states[2]*algebraic[186])/(constants[361]*algebraic[37])+(states[2]*algebraic[126])/(constants[364]*algebraic[40])+(states[2]*algebraic[127])/(constants[123]*constants[306])+(states[2]*algebraic[128])/(constants[92]*algebraic[41])+(states[2]*algebraic[129])/(constants[87]*algebraic[35])+(states[2]*algebraic[187])/(constants[367]*algebraic[42])+(states[2]*algebraic[130])/(constants[94]*algebraic[43])+(states[2]*algebraic[188])/(constants[371]*algebraic[44])+(states[2]*algebraic[160])/(constants[377]*algebraic[48])+(states[2]*algebraic[161])/(constants[379]*algebraic[49])+(states[2]*constants[294])/(constants[141]*constants[284])+(states[2]*constants[293])/(constants[138]*constants[283])+(states[2]*constants[263])/(constants[126]*constants[279])+(states[2]*constants[270])/(constants[129]*constants[280])+(states[2]*constants[269])/(constants[132]*constants[281])+(states[2]*constants[264])/(constants[96]*algebraic[45]))+constants[295]
algebraic[199] = algebraic[181]/constants[3]
algebraic[202] = -((states[3]*constants[292])/(constants[136]*constants[282])+(states[3]*algebraic[183])/(constants[325]*algebraic[0])+(states[3]*algebraic[184])/(constants[340]*algebraic[21])+(states[3]*constants[262])/(constants[349]*algebraic[26])+(states[3]*algebraic[123])/(constants[325]*algebraic[29])+(states[3]*algebraic[124])/(constants[340]*algebraic[31])+(states[3]*algebraic[185])/(constants[355]*algebraic[33])+(states[3]*algebraic[145])/(constants[112]*constants[277])+(states[3]*algebraic[146])/(constants[109]*constants[276])+(states[3]*algebraic[147])/(constants[115]*constants[304])+(states[3]*algebraic[148])/(constants[118]*constants[320])+(states[3]*algebraic[159])/(constants[99]*algebraic[46])+(states[3]*algebraic[131])/(constants[90]*algebraic[36])+(states[3]*algebraic[125])/(constants[121]*constants[278])+(states[3]*algebraic[186])/(constants[362]*algebraic[37])+(states[3]*algebraic[126])/(constants[91]*algebraic[40])+(states[3]*algebraic[127])/(constants[124]*constants[306])+(states[3]*algebraic[128])/(constants[93]*algebraic[41])+(states[3]*algebraic[129])/(constants[88]*algebraic[35])+(states[3]*algebraic[187])/(constants[368]*algebraic[42])+(states[3]*algebraic[130])/(constants[95]*algebraic[43])+(states[3]*algebraic[188])/(constants[372]*algebraic[44])+(states[3]*algebraic[160])/(constants[102]*algebraic[48])+(states[3]*algebraic[161])/(constants[103]*algebraic[49])+(states[3]*constants[294])/(constants[142]*constants[284])+(states[3]*constants[293])/(constants[139]*constants[283])+(states[3]*constants[263])/(constants[127]*constants[279])+(states[3]*constants[270])/(constants[130]*constants[280])+(states[3]*constants[269])/(constants[133]*constants[281])+(states[3]*constants[264])/(constants[97]*algebraic[45]))+algebraic[199]
algebraic[1] = 1.00000+states[5]/constants[351]+states[6]/constants[353]+states[7]/constants[355]
algebraic[2] = 1.00000+states[5]/constants[357]
algebraic[3] = 1.00000+states[6]/constants[359]
algebraic[4] = 1.00000+states[5]/constants[360]+states[6]/constants[361]+states[7]/constants[362]
algebraic[5] = 1.00000+states[5]/constants[363]+states[6]/constants[364]
algebraic[6] = 1.00000+states[5]/constants[365]
algebraic[7] = 1.00000+states[5]/constants[366]+states[6]/constants[367]+states[7]/constants[368]
algebraic[8] = 1.00000+states[5]/constants[369]
algebraic[9] = 1.00000+states[5]/constants[370]+states[6]/constants[371]+states[7]/constants[372]
algebraic[10] = 1.00000+states[5]/constants[373]
algebraic[11] = 1.00000+states[6]/constants[374]
algebraic[12] = 1.00000+states[1]/constants[375]
algebraic[13] = 1.00000+states[5]/constants[376]+states[6]/constants[377]
algebraic[14] = 1.00000+states[5]/constants[378]+states[6]/constants[379]
algebraic[19] = algebraic[17]
algebraic[24] = algebraic[22]
algebraic[133] = 1.00000+algebraic[132]/constants[357]
algebraic[134] = 1.00000+algebraic[132]/constants[365]
algebraic[135] = 1.00000+algebraic[132]/constants[369]
algebraic[136] = 1.00000+algebraic[132]/constants[373]
algebraic[137] = 1.00000+algebraic[132]/constants[375]
algebraic[150] = 1.00000+algebraic[149]/constants[359]
algebraic[151] = 1.00000+algebraic[132]/constants[363]+algebraic[149]/constants[364]
algebraic[164] = algebraic[163]
algebraic[166] = algebraic[165]
algebraic[171] = (states[16]*(1.00000+algebraic[132]/constants[360]))/algebraic[170]
return algebraic
def custom_piecewise(cases):
"""Compute result of a piecewise function"""
return select(cases[0::2],cases[1::2])
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)