# Size of variable arrays: sizeAlgebraic = 187 sizeStates = 60 sizeConstants = 327 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] = "X_pdh in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[3] = "X_cits in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[4] = "X_acon in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[5] = "X_isod in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[6] = "X_akgd in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[7] = "K_ir_akgd in component adjustable_parameters (molar)" legend_constants[8] = "X_scoas in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[9] = "X_sdh in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[10] = "X_fuma in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[11] = "X_mdh in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[12] = "X_ndk in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[13] = "X_got in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[14] = "X_PYRH in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_constants[15] = "X_GLUH in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_constants[16] = "X_CITMAL in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_constants[17] = "X_AKGMAL in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[18] = "X_SUCMAL in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_constants[19] = "X_MALPI in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_constants[20] = "X_ASPGLU in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[21] = "X_C1 in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_constants[22] = "X_C3 in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_half_molar)" legend_constants[23] = "k_PI_1 in component adjustable_parameters (molar)" legend_constants[24] = "k_PI_2 in component adjustable_parameters (molar)" legend_constants[25] = "X_C4 in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar)" legend_constants[26] = "X_F1 in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar)" legend_constants[27] = "X_ANT in component adjustable_parameters (mole_per_second_per_l_mito_volume)" legend_constants[28] = "X_PiHt in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar)" legend_constants[29] = "k_PiHt in component adjustable_parameters (molar)" legend_constants[30] = "X_KH in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_constants[31] = "X_Hle in component adjustable_parameters (mole_per_second_per_l_mito_volume_per_molar_per_mv)" legend_constants[245] = "X_HK in component adjustable_parameters (mole_per_second_per_l_cyto_volume)" legend_constants[32] = "J_AtC in component adjustable_parameters (mole_per_second_per_l_cyto_volume)" legend_constants[249] = "RT in component fixed_parameters (kilojoule_per_mole)" legend_constants[33] = "F in component fixed_parameters (kilojoule_per_mole_per_millivolt)" legend_constants[295] = "V_cyto in component fixed_parameters (cyto_per_cell)" legend_constants[291] = "V_mito in component fixed_parameters (mito_per_cell)" legend_constants[34] = "W_x in component fixed_parameters (l_water_per_l_mito)" legend_constants[35] = "W_i in component fixed_parameters (l_water_per_l_mito)" legend_constants[246] = "W_c in component fixed_parameters (l_water_per_l_cyto)" legend_constants[36] = "gamma in component fixed_parameters (per_micrometer)" legend_constants[37] = "rho_m in component fixed_parameters (mg_per_mitochondria)" legend_constants[38] = "n_A in component fixed_parameters (dimensionless)" legend_constants[39] = "p_PI in component fixed_parameters (micrometer_per_second)" legend_constants[40] = "p_A in component fixed_parameters (micrometer_per_second)" legend_constants[41] = "p_TI in component fixed_parameters (micrometer_per_second)" legend_constants[42] = "k_mADP in component fixed_parameters (molar)" legend_constants[43] = "theta in component fixed_parameters (dimensionless)" legend_constants[44] = "k_O2 in component fixed_parameters (molar)" legend_constants[45] = "NAD_tot in component fixed_parameters (molar)" legend_constants[46] = "Q_tot in component fixed_parameters (molar)" legend_constants[47] = "cytC_tot in component fixed_parameters (molar)" legend_constants[48] = "A_tot in component fixed_parameters (molar)" legend_constants[49] = "CR_tot in component fixed_parameters (molar)" legend_constants[50] = "CO2_tot_x in component fixed_parameters (molar)" legend_constants[51] = "FAD_tot in component fixed_parameters (molar)" legend_constants[52] = "C_IM in component fixed_parameters (mole_per_l_mito_volume_per_mv)" legend_algebraic[157] = "J_C1 in component J_C1 (mole_per_second_per_l_mito_volume)" legend_algebraic[3] = "k_eqC1 in component k_eqC1 (dimensionless)" legend_states[0] = "H_x in component dH_x_dt (molar)" legend_states[1] = "H_i in component dH_i_dt (molar)" legend_states[2] = "deltaPsi in component ddeltaPsi_dt (millivolt)" legend_algebraic[143] = "NAD_x in component dNAD_x_dt (molar)" legend_states[3] = "NADH_x in component dNADH_x_dt (molar)" legend_states[4] = "QH2_x in component dQH2_x_dt (molar)" legend_algebraic[156] = "COQ_x in component dCOQ_x_dt (molar)" legend_algebraic[0] = "k_eqC1_nought in component k_eqC1 (dimensionless)" legend_constants[53] = "deltaG_C1 in component k_eqC1 (kilojoule_per_mole)" legend_algebraic[164] = "J_C3 in component J_C3 (mole_per_second_per_l_mito_volume)" legend_algebraic[5] = "k_eqC3 in component k_eqC3 (dimensionless)" legend_states[5] = "PI_x in component dPI_x_dt (molar)" legend_states[6] = "Cred_i in component dCred_i_dt (molar)" legend_algebraic[163] = "Cox_i in component dCox_i_dt (molar)" legend_algebraic[4] = "k_eqC3_nought in component k_eqC3 (dimensionless)" legend_constants[54] = "deltaG_C3 in component k_eqC3 (kilojoule_per_mole)" legend_algebraic[165] = "J_C4 in component J_C4 (mole_per_second_per_l_mito_volume)" legend_algebraic[7] = "k_eqC4 in component k_eqC4 (dimensionless)" legend_states[7] = "O2_x in component dO2_x_dt (molar)" legend_algebraic[6] = "k_eqC4_nought in component k_eqC4 (dimensionless)" legend_constants[55] = "deltaG_C4 in component k_eqC4 (kilojoule_per_mole)" legend_algebraic[78] = "J_F1 in component J_F1 (mole_per_second_per_l_mito_volume)" legend_algebraic[77] = "k_eqF1 in component k_eqF1 (per_molar)" legend_algebraic[36] = "P_ATP_x in component P_ATP_x (dimensionless)" legend_algebraic[41] = "P_ADP_x in component P_ADP_x (dimensionless)" legend_algebraic[76] = "P_PI_x in component P_PI_x (dimensionless)" legend_states[8] = "ATP_x in component dATP_x_dt (molar)" legend_states[9] = "ADP_x in component dADP_x_dt (molar)" legend_algebraic[8] = "k_eqF1_nought in component k_eqF1 (dimensionless)" legend_constants[56] = "deltaG_F1 in component k_eqF1 (kilojoule_per_mole)" legend_algebraic[144] = "J_pdh in component J_pdh (mole_per_second_per_l_mito_volume)" legend_states[10] = "PYR_x in component dPYR_x_dt (molar)" legend_states[11] = "COASH_x in component dCOASH_x_dt (molar)" legend_states[12] = "CO2_tot_x in component dCO2_tot_x_dt (molar)" legend_states[13] = "ACCOA_x in component dACCOA_x_dt (molar)" legend_algebraic[117] = "P_CO2_tot_x in component P_CO2_tot_x (dimensionless)" legend_algebraic[120] = "P_PYR_x in component P_PYR_x (dimensionless)" legend_algebraic[83] = "P_COASH_x in component P_COASH_x (dimensionless)" legend_algebraic[121] = "K_eq_pdh in component J_pdh (dimensionless)" legend_constants[57] = "K0_eq_pdh in component J_pdh (molar)" legend_constants[58] = "K_mA in component J_pdh (molar)" legend_constants[59] = "K_mB in component J_pdh (molar)" legend_constants[60] = "K_mC in component J_pdh (molar)" legend_algebraic[9] = "alpha_i1 in component J_pdh (dimensionless)" legend_algebraic[10] = "alpha_i2 in component J_pdh (dimensionless)" legend_constants[61] = "K_iACCOA in component J_pdh (molar)" legend_constants[62] = "K_iNADH in component J_pdh (molar)" legend_constants[63] = "minCon in component J_pdh (molar)" legend_algebraic[89] = "J_cits in component J_cits (mole_per_second_per_l_mito_volume)" legend_states[14] = "OAA_x in component dOAA_x_dt (molar)" legend_states[15] = "CIT_x in component dCIT_x_dt (molar)" legend_states[16] = "AMP_x in component dAMP_x_dt (molar)" legend_states[17] = "SCOA_x in component dSCOA_x_dt (molar)" legend_algebraic[38] = "fATP_x in component P_ATP_x (molar)" legend_algebraic[43] = "fADP_x in component P_ADP_x (molar)" legend_algebraic[52] = "fAMP_x in component P_AMP_x (molar)" legend_algebraic[87] = "uCIT_x in component P_CIT_x (molar)" legend_algebraic[85] = "P_CIT_x in component P_CIT_x (dimensionless)" legend_algebraic[84] = "P_OAA_x in component P_OAA_x (dimensionless)" legend_algebraic[86] = "K_eq_cits in component J_cits (dimensionless)" legend_constants[64] = "K0_eq_cits in component J_cits (molar_squared)" legend_constants[65] = "K_mA in component J_cits (molar)" legend_constants[66] = "K_mB in component J_cits (molar)" legend_constants[67] = "K_ia in component J_cits (molar)" legend_constants[68] = "K_iCIT in component J_cits (molar)" legend_constants[69] = "K_iATP in component J_cits (molar)" legend_constants[70] = "K_iADP in component J_cits (molar)" legend_constants[71] = "K_iAMP in component J_cits (molar)" legend_constants[72] = "K_iCOASH in component J_cits (molar)" legend_constants[73] = "K_iSCOA in component J_cits (molar)" legend_algebraic[88] = "alpha_i1 in component J_cits (dimensionless)" legend_algebraic[54] = "alpha_i2 in component J_cits (dimensionless)" legend_algebraic[97] = "J_acon in component J_acon (mole_per_second_per_l_mito_volume)" legend_states[18] = "ICIT_x in component dICIT_x_dt (molar)" legend_algebraic[94] = "P_ICIT_x in component P_ICIT_x (dimensionless)" legend_algebraic[96] = "V_mr in component J_acon (mole_per_second_per_l_mito_volume)" legend_algebraic[95] = "K_eq_acon in component J_acon (dimensionless)" legend_constants[74] = "K0_eq_acon in component J_acon (dimensionless)" legend_constants[75] = "K_mA in component J_acon (molar)" legend_constants[76] = "K_mP in component J_acon (molar)" legend_algebraic[145] = "J_isod in component J_isod (mole_per_second_per_l_mito_volume)" legend_states[19] = "AKG_x in component dAKG_x_dt (molar)" legend_algebraic[118] = "K_eq_isod in component J_isod (molar)" legend_constants[77] = "K0_eq_isod in component J_isod (molar_cubed)" legend_constants[78] = "K_mA in component J_isod (molar)" legend_constants[79] = "K_mB in component J_isod (molar)" legend_constants[80] = "K_ib in component J_isod (molar)" legend_constants[81] = "K_iq in component J_isod (molar)" legend_algebraic[46] = "alpha_i in component J_isod (dimensionless)" legend_constants[82] = "n_H in component J_isod (dimensionless)" legend_constants[83] = "K_aATP in component J_isod (molar)" legend_constants[84] = "K_aADP in component J_isod (molar)" legend_constants[85] = "minCon in component J_isod (molar)" legend_algebraic[146] = "J_akgd in component J_akgd (mole_per_second_per_l_mito_volume)" legend_algebraic[98] = "P_SCOA_x in component P_SCOA_x (dimensionless)" legend_algebraic[119] = "K_eq_akgd in component J_akgd (dimensionless)" legend_constants[86] = "K0_eq_akgd in component J_akgd (molar)" legend_constants[87] = "K_mA in component J_akgd (molar)" legend_constants[88] = "K_mB in component J_akgd (molar)" legend_constants[89] = "K_mC in component J_akgd (molar)" legend_constants[90] = "K_iq in component J_akgd (molar)" legend_constants[91] = "K_ir in component J_akgd (molar)" legend_constants[92] = "K_ir2 in component J_akgd (molar)" legend_constants[93] = "K_iATP in component J_akgd (molar)" legend_constants[94] = "K_aADP in component J_akgd (molar)" legend_algebraic[48] = "alpha_i in component J_akgd (dimensionless)" legend_constants[95] = "minCon in component J_akgd (molar)" legend_algebraic[102] = "J_scoas in component J_scoas (mole_per_second_per_l_mito_volume)" legend_states[20] = "GDP_x in component dGDP_x_dt (molar)" legend_states[21] = "SUC_x in component dSUC_x_dt (molar)" legend_states[22] = "GTP_x in component dGTP_x_dt (molar)" legend_algebraic[99] = "P_SUC_x in component P_SUC_x (dimensionless)" legend_algebraic[66] = "P_GTP_x in component P_GTP_x (dimensionless)" legend_algebraic[72] = "P_GDP_x in component P_GDP_x (dimensionless)" legend_algebraic[100] = "K_eq_scoas in component J_scoas (dimensionless)" legend_constants[96] = "K0_eq_scoas in component J_scoas (molar)" legend_algebraic[101] = "V_mr in component J_scoas (mole_per_second_per_l_mito_volume)" legend_constants[97] = "K_iGDP in component J_scoas (molar)" legend_constants[98] = "K_iSCOA in component J_scoas (molar)" legend_constants[99] = "K_iPI in component J_scoas (molar)" legend_constants[100] = "K_iCOASH in component J_scoas (molar)" legend_constants[101] = "K_iSUC in component J_scoas (molar)" legend_constants[102] = "K_iGTP in component J_scoas (molar)" legend_constants[103] = "K_mGDP in component J_scoas (molar)" legend_constants[104] = "K_mSCOA in component J_scoas (molar)" legend_constants[105] = "K_mPI in component J_scoas (molar)" legend_constants[106] = "K_mCOASH in component J_scoas (molar)" legend_constants[107] = "K_mSUC in component J_scoas (molar)" legend_constants[108] = "K_mGTP in component J_scoas (molar)" legend_algebraic[158] = "J_sdh in component J_sdh (mole_per_second_per_l_mito_volume)" legend_states[23] = "FUM_x in component dFUM_x_dt (molar)" legend_algebraic[104] = "P_FUM_x in component P_FUM_x (dimensionless)" legend_algebraic[105] = "K_eq_sdh in component J_sdh (dimensionless)" legend_constants[109] = "K0_eq_sdh in component J_sdh (dimensionless)" legend_algebraic[106] = "V_mr in component J_sdh (mole_per_second_per_l_mito_volume)" legend_constants[110] = "K_ia in component J_sdh (molar)" legend_constants[111] = "K_iq in component J_sdh (molar)" legend_constants[112] = "K_mSUC in component J_sdh (molar)" legend_constants[113] = "K_mCOQ in component J_sdh (molar)" legend_constants[114] = "K_mQH2 in component J_sdh (molar)" legend_constants[115] = "K_mFUM in component J_sdh (molar)" legend_algebraic[11] = "alpha_i in component J_sdh (dimensionless)" legend_constants[116] = "K_iOAA in component J_sdh (molar)" legend_constants[117] = "K_aSUC in component J_sdh (molar)" legend_constants[118] = "K_aFUM in component J_sdh (molar)" legend_algebraic[110] = "J_fum in component J_fum (mole_per_second_per_l_mito_volume)" legend_states[24] = "MAL_x in component dMAL_x_dt (molar)" legend_algebraic[68] = "fGTP_x in component P_GTP_x (molar)" legend_algebraic[74] = "fGDP_x in component P_GDP_x (molar)" legend_algebraic[107] = "P_MAL_x in component P_MAL_x (dimensionless)" legend_algebraic[108] = "K_eq_fum in component J_fum (dimensionless)" legend_constants[119] = "K0_eq_fum in component J_fum (dimensionless)" legend_algebraic[109] = "V_mr in component J_fum (mole_per_second_per_l_mito_volume)" legend_constants[120] = "K_mFUM in component J_fum (molar)" legend_constants[121] = "K_mMAL in component J_fum (molar)" legend_algebraic[75] = "alpha_i in component J_fum (dimensionless)" legend_constants[122] = "K_iCIT in component J_fum (molar)" legend_constants[123] = "K_iATP in component J_fum (molar)" legend_constants[124] = "K_iADP in component J_fum (molar)" legend_constants[125] = "K_iGTP in component J_fum (molar)" legend_constants[126] = "K_iGDP in component J_fum (molar)" legend_algebraic[147] = "J_mdh in component J_mdh (mole_per_second_per_l_mito_volume)" legend_algebraic[111] = "K_eq_mdh in component J_mdh (dimensionless)" legend_constants[127] = "K0_eq_mdh in component J_mdh (molar)" legend_algebraic[112] = "V_mr in component J_mdh (mole_per_second_per_l_mito_volume)" legend_constants[128] = "K_ia in component J_mdh (molar)" legend_constants[129] = "K_ib in component J_mdh (molar)" legend_constants[130] = "K_ip in component J_mdh (molar)" legend_constants[131] = "K_iq in component J_mdh (molar)" legend_constants[132] = "K_mNAD in component J_mdh (molar)" legend_constants[133] = "K_mMAL in component J_mdh (molar)" legend_constants[134] = "K_mOAA in component J_mdh (molar)" legend_constants[135] = "K_mNADH in component J_mdh (molar)" legend_algebraic[55] = "alpha_i in component J_mdh (dimensionless)" legend_constants[136] = "K_iATP in component J_mdh (molar)" legend_constants[137] = "K_iADP in component J_mdh (molar)" legend_constants[138] = "K_iAMP in component J_mdh (molar)" legend_algebraic[57] = "J_ndk_f in component J_ndk (mole_per_second_per_l_mito_volume)" legend_algebraic[58] = "J_ndk in component J_ndk (mole_per_second_per_l_mito_volume)" legend_constants[139] = "K_eq_ndk in component J_ndk (dimensionless)" legend_constants[247] = "V_mr in component J_ndk (mole_per_second_per_l_mito_volume)" legend_constants[140] = "K_ia in component J_ndk (molar)" legend_constants[141] = "K_ib in component J_ndk (molar)" legend_constants[142] = "K_ip in component J_ndk (molar)" legend_constants[143] = "K_iq in component J_ndk (molar)" legend_constants[144] = "K_mA in component J_ndk (molar)" legend_constants[145] = "K_mB in component J_ndk (molar)" legend_constants[146] = "K_mP in component J_ndk (molar)" legend_constants[147] = "K_mQ in component J_ndk (molar)" legend_algebraic[56] = "alpha_i in component J_ndk (dimensionless)" legend_constants[148] = "K_iAMP in component J_ndk (molar)" legend_constants[149] = "minCon in component J_ndk (molar)" legend_algebraic[130] = "J_got_f in component J_got (mole_per_second_per_l_mito_volume)" legend_algebraic[131] = "J_got in component J_got (mole_per_second_per_l_mito_volume)" legend_states[25] = "ASP_x in component dASP_x_dt (molar)" legend_states[26] = "GLU_x in component dGLU_x_dt (molar)" legend_algebraic[127] = "P_ASP_x in component P_ASP_x (dimensionless)" legend_algebraic[124] = "P_GLU_x in component P_GLU_x (dimensionless)" legend_algebraic[128] = "K_eq_got in component J_got (dimensionless)" legend_constants[150] = "K0_eq_got in component J_got (dimensionless)" legend_algebraic[129] = "V_mr in component J_got (mole_per_second_per_l_mito_volume)" legend_constants[151] = "K_ia in component J_got (molar)" legend_constants[152] = "K_ib in component J_got (molar)" legend_constants[153] = "K_ip in component J_got (molar)" legend_constants[154] = "K_iq in component J_got (molar)" legend_constants[155] = "K_mA in component J_got (molar)" legend_constants[156] = "K_mB in component J_got (molar)" legend_constants[157] = "K_mP in component J_got (molar)" legend_constants[158] = "K_mQ in component J_got (molar)" legend_algebraic[12] = "alpha_i in component J_got (dimensionless)" legend_constants[159] = "K_iAKG in component J_got (molar)" legend_constants[160] = "minCon in component J_got (molar)" legend_algebraic[123] = "J_PYRH in component J_PYRH (mole_per_second_per_l_mito_volume)" legend_states[27] = "PYR_i in component dPYR_i_dt (molar)" legend_algebraic[122] = "P_PYR_i in component P_PYR_i (dimensionless)" legend_algebraic[126] = "J_GLUH in component J_GLUH (mole_per_second_per_l_mito_volume)" legend_states[28] = "GLU_i in component dGLU_i_dt (molar)" legend_algebraic[125] = "P_GLU_i in component P_GLU_i (dimensionless)" legend_algebraic[114] = "J_CITMAL in component J_CITMAL (mole_per_second_per_l_mito_volume)" legend_algebraic[92] = "hCIT_i in component P_CIT_i (molar)" legend_states[29] = "MAL_i in component dMAL_i_dt (molar)" legend_algebraic[113] = "P_MAL_i in component P_MAL_i (dimensionless)" legend_algebraic[90] = "hCIT_x in component P_CIT_x (molar)" legend_algebraic[13] = "J_AKGMAL in component J_AKGMAL (mole_per_second_per_l_mito_volume)" legend_constants[161] = "K_mAKGi in component J_AKGMAL (molar)" legend_constants[162] = "K_mAKGx in component J_AKGMAL (molar)" legend_constants[163] = "K_mMALi in component J_AKGMAL (molar)" legend_constants[164] = "K_mMALx in component J_AKGMAL (molar)" legend_states[30] = "AKG_i in component dAKG_i_dt (molar)" legend_algebraic[115] = "J_MALPI in component J_MALPI (mole_per_second_per_l_mito_volume)" legend_states[31] = "PI_i in component dPI_i_dt (molar)" legend_algebraic[80] = "P_PI_i in component P_PI_i (dimensionless)" legend_algebraic[134] = "J_ASPGLU in component J_ASPGLU (mole_per_second_per_l_mito_volume)" legend_algebraic[133] = "K_eqASPGLU in component J_ASPGLU (dimensionless)" legend_constants[248] = "K_hASPGLU in component J_ASPGLU (molar)" legend_constants[165] = "K_iASPi in component J_ASPGLU (molar)" legend_constants[166] = "K_iASPx in component J_ASPGLU (molar)" legend_constants[167] = "K_iGLUi in component J_ASPGLU (molar)" legend_constants[168] = "K_iGLUx in component J_ASPGLU (molar)" legend_states[32] = "ASP_i in component dASP_i_dt (molar)" legend_constants[169] = "m in component J_ASPGLU (dimensionless)" legend_algebraic[132] = "P_ASP_i in component P_ASP_i (dimensionless)" legend_algebraic[116] = "J_SUCMAL in component J_SUCMAL (mole_per_second_per_l_mito_volume)" legend_states[33] = "SUC_i in component dSUC_i_dt (molar)" legend_algebraic[103] = "P_SUC_i in component P_SUC_i (dimensionless)" legend_algebraic[14] = "J_ATPt in component J_ATPt (mole_per_second_per_l_mito_volume)" legend_states[34] = "ATP_c in component dATP_c_dt (molar)" legend_states[35] = "ATP_i in component dATP_i_dt (molar)" legend_algebraic[16] = "J_ADPt in component J_ADPt (mole_per_second_per_l_mito_volume)" legend_states[36] = "ADP_c in component dADP_c_dt (molar)" legend_states[37] = "ADP_i in component dADP_i_dt (molar)" legend_algebraic[17] = "J_AMPt in component J_AMPt (mole_per_second_per_l_mito_volume)" legend_states[38] = "AMP_c in component dAMP_c_dt (molar)" legend_states[39] = "AMP_i in component dAMP_i_dt (molar)" legend_algebraic[18] = "J_PIt in component J_PIt (mole_per_second_per_l_mito_volume)" legend_states[40] = "PI_c in component dPI_c_dt (molar)" legend_algebraic[20] = "J_PYRt in component J_PYRt (mole_per_second_per_l_mito_volume)" legend_states[41] = "PYR_c in component dPYR_c_dt (molar)" legend_algebraic[22] = "J_CITt in component J_CITt (mole_per_second_per_l_mito_volume)" legend_states[42] = "CIT_c in component dCIT_c_dt (molar)" legend_states[43] = "CIT_i in component dCIT_i_dt (molar)" legend_algebraic[1] = "J_ICITt in component J_ICITt (mole_per_second_per_l_mito_volume)" legend_states[44] = "ICIT_c in component dICIT_c_dt (molar)" legend_states[45] = "ICIT_i in component dICIT_i_dt (molar)" legend_algebraic[24] = "J_MALt in component J_MALt (mole_per_second_per_l_mito_volume)" legend_states[46] = "MAL_c in component dMAL_c_dt (molar)" legend_algebraic[15] = "J_AKGt in component J_AKGt (mole_per_second_per_l_mito_volume)" legend_states[47] = "AKG_c in component dAKG_c_dt (molar)" legend_algebraic[26] = "J_SUCt in component J_SUCt (mole_per_second_per_l_mito_volume)" legend_states[48] = "SUC_c in component dSUC_c_dt (molar)" legend_algebraic[28] = "J_GLUt in component J_GLUt (mole_per_second_per_l_mito_volume)" legend_states[49] = "GLU_c in component dGLU_c_dt (molar)" legend_algebraic[30] = "J_ASPt in component J_ASPt (mole_per_second_per_l_mito_volume)" legend_states[50] = "ASP_c in component dASP_c_dt (molar)" legend_algebraic[2] = "J_FUMt in component J_FUMt (mole_per_second_per_l_mito_volume)" legend_constants[170] = "p_TI in component J_FUMt (micrometer_per_second)" legend_states[51] = "FUM_c in component dFUM_c_dt (molar)" legend_states[52] = "FUM_i in component dFUM_i_dt (molar)" legend_algebraic[63] = "J_ANT in component J_ANT (mole_per_second_per_l_mito_volume)" legend_algebraic[60] = "P_ADP_i in component P_ADP_i (dimensionless)" legend_algebraic[59] = "P_ATP_i in component P_ATP_i (dimensionless)" legend_constants[171] = "minCond in component J_ANT (molar)" legend_algebraic[82] = "J_PIHt in component J_PIHt (mole_per_second_per_l_mito_volume)" legend_algebraic[81] = "hPI_i in component P_PI_i (molar)" legend_algebraic[79] = "hPI_x in component P_PI_x (molar)" legend_algebraic[32] = "J_Hle in component J_Hle (mole_per_second_per_l_mito_volume)" legend_constants[172] = "minCond in component J_Hle (millivolt)" legend_algebraic[34] = "J_KH in component J_KH (mole_per_second_per_l_mito_volume)" legend_states[53] = "K_i in component dK_i_dt (molar)" legend_states[54] = "K_x in component dK_x_dt (molar)" legend_algebraic[73] = "J_AKi in component J_AKi (mole_per_second_per_l_mito_volume)" legend_constants[173] = "X_AK in component J_AKi (mole_per_second_per_l_mito_volume_per_molar_per_molar)" legend_algebraic[69] = "K_AKi in component J_AKi (dimensionless)" legend_constants[174] = "K0_AKi in component J_AKi (dimensionless)" legend_algebraic[67] = "P_AMP_i in component P_AMP_i (dimensionless)" legend_algebraic[47] = "J_AKc in component J_AKc (mole_per_second_per_l_cyto_volume)" legend_constants[175] = "X_AK in component J_AKc (mole_per_second_per_l_cyto_volume_per_molar_per_molar)" legend_algebraic[44] = "K_AKc in component J_AKc (dimensionless)" legend_constants[176] = "K0_AKc in component J_AKc (dimensionless)" legend_algebraic[27] = "P_ATP_c in component P_ATP_c (dimensionless)" legend_algebraic[31] = "P_ADP_c in component P_ADP_c (dimensionless)" legend_algebraic[42] = "P_AMP_c in component P_AMP_c (dimensionless)" legend_algebraic[51] = "J_CKc in component J_CKc (mole_per_second_per_l_cyto_volume)" legend_constants[177] = "X_CK in component J_CKc (mole_per_second_per_l_cyto_volume_per_molar_per_molar)" legend_constants[293] = "K_CK in component J_CKc (per_molar)" legend_algebraic[19] = "H_c in component dH_c_dt (molar)" legend_states[55] = "PCr_c in component dPCr_c_dt (molar)" legend_algebraic[49] = "Cr_c in component dCr_c_dt (molar)" legend_algebraic[39] = "J_HK in component J_HK (mole_per_second_per_l_cyto_volume)" legend_algebraic[29] = "mATP_c in component P_ATP_c (molar)" legend_algebraic[37] = "mADP_c in component P_ADP_c (molar)" legend_states[56] = "GLC_c in component dGLC_c_dt (molar)" legend_states[57] = "G6P_c in component dG6P_c_dt (molar)" legend_constants[178] = "K_ia in component J_HK (molar)" legend_constants[179] = "K_ib in component J_HK (molar)" legend_constants[180] = "K_ip in component J_HK (molar)" legend_constants[181] = "K_iq in component J_HK (molar)" legend_constants[182] = "K_mA in component J_HK (molar)" legend_constants[183] = "K_mB in component J_HK (molar)" legend_constants[184] = "K_mP in component J_HK (molar)" legend_constants[185] = "K_mQ in component J_HK (molar)" legend_constants[186] = "Ki_G6P in component J_HK (molar)" legend_algebraic[21] = "P_G6P_c in component P_G6P_c (dimensionless)" legend_constants[300] = "K_Mg in component P_ATP_x (molar)" legend_constants[306] = "K_Mg in component P_ADP_x (molar)" legend_constants[292] = "Keq_HK in component J_HK (molar)" legend_algebraic[33] = "Kapp_HK in component J_HK (dimensionless)" legend_algebraic[35] = "Kapp_HK_m in component J_HK (dimensionless)" legend_algebraic[40] = "ATP_xref in component P_ATP_x (molar)" legend_states[58] = "Mg_x in component dMg_x_dt (molar)" legend_constants[296] = "K_H in component P_ATP_x (molar)" legend_constants[302] = "K_K in component P_ATP_x (molar)" legend_algebraic[45] = "ADP_xref in component P_ADP_x (molar)" legend_constants[304] = "K_H in component P_ADP_x (molar)" legend_constants[308] = "K_K in component P_ADP_x (molar)" legend_algebraic[50] = "P_AMP_x in component P_AMP_x (dimensionless)" legend_algebraic[53] = "AMP_xref in component P_AMP_x (molar)" legend_constants[310] = "K_H in component P_AMP_x (molar)" legend_constants[312] = "K_Mg in component P_AMP_x (molar)" legend_constants[314] = "K_K in component P_AMP_x (molar)" legend_algebraic[61] = "ATP_iref in component P_ATP_i (molar)" legend_algebraic[62] = "fATP_i in component P_ATP_i (molar)" legend_states[59] = "Mg_i in component dMg_i_dt (molar)" legend_algebraic[64] = "ADP_iref in component P_ADP_i (molar)" legend_algebraic[65] = "fADP_i in component P_ADP_i (molar)" legend_algebraic[70] = "AMP_iref in component P_AMP_i (molar)" legend_algebraic[71] = "fAMP_i in component P_AMP_i (molar)" legend_algebraic[23] = "Mg_c in component dMg_c_dt (molar)" legend_algebraic[25] = "K_c in component dK_c_dt (molar)" legend_constants[250] = "K_H in component P_PI_x (molar)" legend_constants[251] = "K_Mg in component P_PI_x (molar)" legend_constants[252] = "K_K in component P_PI_x (molar)" legend_constants[253] = "K_H in component P_COASH_x (molar)" legend_constants[187] = "K_Mg in component P_COASH_x (molar)" legend_constants[188] = "K_K in component P_COASH_x (molar)" legend_constants[254] = "K_Mg in component P_OAA_x (molar)" legend_constants[189] = "K_H in component P_OAA_x (molar)" legend_constants[190] = "K_K in component P_OAA_x (molar)" legend_constants[255] = "K_H in component P_CIT_x (molar)" legend_constants[256] = "K_Mg in component P_CIT_x (molar)" legend_constants[257] = "K_K in component P_CIT_x (molar)" legend_algebraic[91] = "P_CIT_i in component P_CIT_i (dimensionless)" legend_algebraic[93] = "uCIT_i in component P_CIT_i (molar)" legend_constants[258] = "K_H in component P_ICIT_x (molar)" legend_constants[259] = "K_Mg in component P_ICIT_x (molar)" legend_constants[191] = "K_K in component P_ICIT_x (molar)" legend_constants[260] = "K_H in component P_SCOA_x (molar)" legend_constants[192] = "K_Mg in component P_SCOA_x (molar)" legend_constants[193] = "K_K in component P_SCOA_x (molar)" legend_constants[261] = "K_H in component P_SUC_x (molar)" legend_constants[262] = "K_Mg in component P_SUC_x (molar)" legend_constants[263] = "K_K in component P_SUC_x (molar)" legend_constants[264] = "K_H in component P_FUM_x (molar)" legend_constants[194] = "K_Mg in component P_FUM_x (molar)" legend_constants[195] = "K_K in component P_FUM_x (molar)" legend_constants[265] = "K_H in component P_MAL_x (molar)" legend_constants[266] = "K_Mg in component P_MAL_x (molar)" legend_constants[267] = "K_K in component P_MAL_x (molar)" legend_constants[268] = "K_H in component P_CO2_tot_x (molar)" legend_constants[196] = "K_Mg in component P_CO2_tot_x (molar)" legend_constants[197] = "K_K in component P_CO2_tot_x (molar)" legend_constants[269] = "K_Mg in component P_PYR_x (molar)" legend_constants[198] = "K_H in component P_PYR_x (molar)" legend_constants[199] = "K_K in component P_PYR_x (molar)" legend_constants[270] = "K_H in component P_GLU_x (molar)" legend_constants[271] = "K_Mg in component P_GLU_x (molar)" legend_constants[200] = "K_K in component P_GLU_x (molar)" legend_constants[272] = "K_H in component P_ASP_x (molar)" legend_constants[273] = "K_Mg in component P_ASP_x (molar)" legend_constants[201] = "K_K in component P_ASP_x (molar)" legend_constants[274] = "P_NAD_x in component P_NAD_x (dimensionless)" legend_constants[202] = "K_H in component P_NAD_x (molar)" legend_constants[203] = "K_Mg in component P_NAD_x (molar)" legend_constants[204] = "K_K in component P_NAD_x (molar)" legend_constants[275] = "P_NADH_x in component P_NADH_x (dimensionless)" legend_constants[205] = "K_H in component P_NADH_x (molar)" legend_constants[206] = "K_Mg in component P_NADH_x (molar)" legend_constants[207] = "K_K in component P_NADH_x (molar)" legend_constants[276] = "P_QH2_x in component P_QH2_x (dimensionless)" legend_constants[208] = "K_H in component P_QH2_x (molar)" legend_constants[209] = "K_Mg in component P_QH2_x (molar)" legend_constants[210] = "K_K in component P_QH2_x (molar)" legend_constants[277] = "P_COQ_x in component P_COQ_x (dimensionless)" legend_constants[211] = "K_H in component P_COQ_x (molar)" legend_constants[212] = "K_Mg in component P_COQ_x (molar)" legend_constants[213] = "K_K in component P_COQ_x (molar)" legend_constants[278] = "P_ACCOA_x in component P_ACCOA_x (dimensionless)" legend_constants[214] = "K_H in component P_ACCOA_x (molar)" legend_constants[215] = "K_Mg in component P_ACCOA_x (molar)" legend_constants[216] = "K_K in component P_ACCOA_x (molar)" legend_constants[279] = "P_AKG_x in component P_AKG_x (dimensionless)" legend_constants[217] = "K_H in component P_AKG_x (molar)" legend_constants[218] = "K_Mg in component P_AKG_x (molar)" legend_constants[219] = "K_K in component P_AKG_x (molar)" legend_constants[280] = "P_O2_x in component P_O2_x (dimensionless)" legend_constants[220] = "K_H in component P_O2_x (molar)" legend_constants[221] = "K_Mg in component P_O2_x (molar)" legend_constants[222] = "K_K in component P_O2_x (molar)" legend_constants[281] = "P_FADH2_x in component P_FADH2_x (dimensionless)" legend_constants[223] = "K_H in component P_FADH2_x (molar)" legend_constants[224] = "K_Mg in component P_FADH2_x (molar)" legend_constants[225] = "K_K in component P_FADH2_x (molar)" legend_constants[282] = "P_FAD_x in component P_FAD_x (dimensionless)" legend_constants[226] = "K_H in component P_FAD_x (molar)" legend_constants[227] = "K_Mg in component P_FAD_x (molar)" legend_constants[228] = "K_K in component P_FAD_x (molar)" legend_constants[283] = "P_H_x in component P_H_x (dimensionless)" legend_constants[229] = "K_H in component P_H_x (molar)" legend_constants[230] = "K_Mg in component P_H_x (molar)" legend_constants[231] = "K_K in component P_H_x (molar)" legend_constants[284] = "P_Mg_x in component P_Mg_x (dimensionless)" legend_constants[232] = "K_H in component P_Mg_x (molar)" legend_constants[233] = "K_Mg in component P_Mg_x (molar)" legend_constants[234] = "K_K in component P_Mg_x (molar)" legend_constants[285] = "P_K_x in component P_K_x (dimensionless)" legend_constants[235] = "K_H in component P_K_x (molar)" legend_constants[236] = "K_Mg in component P_K_x (molar)" legend_constants[237] = "K_K in component P_K_x (molar)" legend_constants[316] = "K_H in component P_G6P_c (molar)" legend_algebraic[135] = "dATP_x_dt in component dATP_x_dt (molar_per_second)" legend_algebraic[136] = "dADP_x_dt in component dADP_x_dt (molar_per_second)" legend_constants[238] = "dAMP_x_dt in component dAMP_x_dt (molar_per_second)" legend_algebraic[137] = "dGTP_x_dt in component dGTP_x_dt (molar_per_second)" legend_algebraic[138] = "dGDP_x_dt in component dGDP_x_dt (molar_per_second)" legend_algebraic[139] = "dPI_x_dt in component dPI_x_dt (molar_per_second)" legend_algebraic[159] = "dNADH_x_dt in component dNADH_x_dt (molar_per_second)" legend_algebraic[166] = "dQH2_x_dt in component dQH2_x_dt (molar_per_second)" legend_algebraic[148] = "dPYR_x_dt in component dPYR_x_dt (molar_per_second)" legend_algebraic[149] = "dACCOA_x_dt in component dACCOA_x_dt (molar_per_second)" legend_algebraic[140] = "dCIT_x_dt in component dCIT_x_dt (molar_per_second)" legend_algebraic[150] = "dICIT_x_dt in component dICIT_x_dt (molar_per_second)" legend_algebraic[151] = "dAKG_x_dt in component dAKG_x_dt (molar_per_second)" legend_algebraic[152] = "dSCOA_x_dt in component dSCOA_x_dt (molar_per_second)" legend_algebraic[153] = "dCOASH_x_dt in component dCOASH_x_dt (molar_per_second)" legend_algebraic[160] = "dSUC_x_dt in component dSUC_x_dt (molar_per_second)" legend_algebraic[161] = "dFUM_x_dt in component dFUM_x_dt (molar_per_second)" legend_algebraic[154] = "dMAL_x_dt in component dMAL_x_dt (molar_per_second)" legend_algebraic[155] = "dOAA_x_dt in component dOAA_x_dt (molar_per_second)" legend_algebraic[141] = "dGLU_x_dt in component dGLU_x_dt (molar_per_second)" legend_algebraic[142] = "dASP_x_dt in component dASP_x_dt (molar_per_second)" legend_constants[239] = "dO2_x_dt in component dO2_x_dt (molar_per_second)" legend_constants[240] = "dCO2_tot_x_dt in component dCO2_tot_x_dt (molar_per_second)" legend_algebraic[162] = "dNAD_x_dt in component dNAD_x_dt (molar_per_second)" legend_algebraic[167] = "dCOQ_x_dt in component dCOQ_x_dt (molar_per_second)" legend_constants[286] = "FAD_x in component dFAD_x_dt (molar)" legend_constants[241] = "dFAD_x_dt in component dFAD_x_dt (molar_per_second)" legend_constants[242] = "dFADH2_x_dt in component dFAD_x_dt (molar_per_second)" legend_constants[294] = "FADH2_x in component dFAD_x_dt (molar)" legend_constants[287] = "dH_x_dt in component dH_x_dt (molar_per_second)" legend_algebraic[186] = "D in component BINDING_IONS (dimensionless)" legend_algebraic[180] = "phi_H in component BINDING_IONS (molar_per_second)" legend_algebraic[181] = "phi_Mg in component BINDING_IONS (molar_per_second)" legend_algebraic[182] = "phi_K in component BINDING_IONS (molar_per_second)" legend_algebraic[168] = "dHB_dMg in component BINDING_IONS (dimensionless)" legend_algebraic[169] = "dHB_dK in component BINDING_IONS (dimensionless)" legend_algebraic[172] = "dMgB_dK in component BINDING_IONS (dimensionless)" legend_algebraic[175] = "dKB_dMg in component BINDING_IONS (dimensionless)" legend_algebraic[185] = "alpha_K in component BINDING_IONS (dimensionless)" legend_algebraic[184] = "alpha_Mg in component BINDING_IONS (dimensionless)" legend_constants[288] = "dMg_x_dt in component dMg_x_dt (molar_per_second)" legend_algebraic[183] = "alpha_H in component BINDING_IONS (dimensionless)" legend_algebraic[174] = "dKB_dH in component BINDING_IONS (dimensionless)" legend_algebraic[171] = "dMgB_dH in component BINDING_IONS (dimensionless)" legend_constants[289] = "dK_x_dt in component dK_x_dt (molar_per_second)" legend_constants[297] = "Rm_cyto in component dATP_c_dt (mito_per_cyto)" legend_constants[301] = "Rm_cyto in component dADP_c_dt (mito_per_cyto)" legend_constants[303] = "Rm_cyto in component dPI_c_dt (mito_per_cyto)" legend_constants[305] = "Rm_cyto in component dPYR_c_dt (mito_per_cyto)" legend_constants[307] = "Rm_cyto in component dCIT_c_dt (mito_per_cyto)" legend_constants[298] = "Rm_cyto in component dICIT_c_dt (mito_per_cyto)" legend_constants[309] = "Rm_cyto in component dAKG_c_dt (mito_per_cyto)" legend_constants[311] = "Rm_cyto in component dSUC_c_dt (mito_per_cyto)" legend_constants[299] = "Rm_cyto in component dFUM_c_dt (mito_per_cyto)" legend_constants[313] = "Rm_cyto in component dMAL_c_dt (mito_per_cyto)" legend_constants[315] = "Rm_cyto in component dGLU_c_dt (mito_per_cyto)" legend_constants[317] = "Rm_cyto in component dASP_c_dt (mito_per_cyto)" legend_constants[318] = "multiplier in component dPCr_c_dt (dimensionless)" legend_constants[319] = "multiplier in component dGLC_c_dt (dimensionless)" legend_constants[320] = "multiplier in component dG6P_c_dt (dimensionless)" legend_algebraic[170] = "dHB_dH in component BINDING_IONS (dimensionless)" legend_algebraic[173] = "dMgB_dMg in component BINDING_IONS (dimensionless)" legend_algebraic[176] = "dKB_dK in component BINDING_IONS (dimensionless)" legend_algebraic[177] = "nK_JK in component BINDING_IONS (molar_per_second)" legend_algebraic[178] = "J_H_t in component BINDING_IONS (molar_per_second)" legend_constants[290] = "J_Mg_t in component BINDING_IONS (molar_per_second)" legend_algebraic[179] = "J_K_t in component BINDING_IONS (molar_per_second)" legend_constants[243] = "B_x in component BINDING_IONS (molar)" legend_constants[244] = "K_BX in component BINDING_IONS (molar)" legend_rates[2] = "d/dt deltaPsi in component ddeltaPsi_dt (millivolt)" legend_rates[8] = "d/dt ATP_x in component dATP_x_dt (molar)" legend_rates[9] = "d/dt ADP_x in component dADP_x_dt (molar)" legend_rates[16] = "d/dt AMP_x in component dAMP_x_dt (molar)" legend_rates[22] = "d/dt GTP_x in component dGTP_x_dt (molar)" legend_rates[20] = "d/dt GDP_x in component dGDP_x_dt (molar)" legend_rates[5] = "d/dt PI_x in component dPI_x_dt (molar)" legend_rates[3] = "d/dt NADH_x in component dNADH_x_dt (molar)" legend_rates[4] = "d/dt QH2_x in component dQH2_x_dt (molar)" legend_rates[10] = "d/dt PYR_x in component dPYR_x_dt (molar)" legend_rates[13] = "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[18] = "d/dt ICIT_x in component dICIT_x_dt (molar)" legend_rates[19] = "d/dt AKG_x in component dAKG_x_dt (molar)" legend_rates[17] = "d/dt SCOA_x in component dSCOA_x_dt (molar)" legend_rates[11] = "d/dt COASH_x in component dCOASH_x_dt (molar)" legend_rates[21] = "d/dt SUC_x in component dSUC_x_dt (molar)" legend_rates[23] = "d/dt FUM_x in component dFUM_x_dt (molar)" legend_rates[24] = "d/dt MAL_x in component dMAL_x_dt (molar)" legend_rates[14] = "d/dt OAA_x in component dOAA_x_dt (molar)" legend_rates[26] = "d/dt GLU_x in component dGLU_x_dt (molar)" legend_rates[25] = "d/dt ASP_x in component dASP_x_dt (molar)" legend_rates[7] = "d/dt O2_x in component dO2_x_dt (molar)" legend_rates[12] = "d/dt CO2_tot_x in component dCO2_tot_x_dt (molar)" legend_rates[0] = "d/dt H_x in component dH_x_dt (molar)" legend_rates[58] = "d/dt Mg_x in component dMg_x_dt (molar)" legend_rates[54] = "d/dt K_x in component dK_x_dt (molar)" legend_rates[1] = "d/dt H_i in component dH_i_dt (molar)" legend_rates[59] = "d/dt Mg_i in component dMg_i_dt (molar)" legend_rates[53] = "d/dt K_i in component dK_i_dt (molar)" legend_rates[6] = "d/dt Cred_i in component dCred_i_dt (molar)" legend_rates[35] = "d/dt ATP_i in component dATP_i_dt (molar)" legend_rates[37] = "d/dt ADP_i in component dADP_i_dt (molar)" legend_rates[39] = "d/dt AMP_i in component dAMP_i_dt (molar)" legend_rates[31] = "d/dt PI_i in component dPI_i_dt (molar)" legend_rates[27] = "d/dt PYR_i in component dPYR_i_dt (molar)" legend_rates[43] = "d/dt CIT_i in component dCIT_i_dt (molar)" legend_rates[45] = "d/dt ICIT_i in component dICIT_i_dt (molar)" legend_rates[30] = "d/dt AKG_i in component dAKG_i_dt (molar)" legend_rates[33] = "d/dt SUC_i in component dSUC_i_dt (molar)" legend_rates[52] = "d/dt FUM_i in component dFUM_i_dt (molar)" legend_rates[29] = "d/dt MAL_i in component dMAL_i_dt (molar)" legend_rates[28] = "d/dt GLU_i in component dGLU_i_dt (molar)" legend_rates[32] = "d/dt ASP_i in component dASP_i_dt (molar)" legend_rates[34] = "d/dt ATP_c in component dATP_c_dt (molar)" legend_rates[36] = "d/dt ADP_c in component dADP_c_dt (molar)" legend_rates[38] = "d/dt AMP_c in component dAMP_c_dt (molar)" legend_rates[40] = "d/dt PI_c in component dPI_c_dt (molar)" legend_rates[41] = "d/dt PYR_c in component dPYR_c_dt (molar)" legend_rates[42] = "d/dt CIT_c in component dCIT_c_dt (molar)" legend_rates[44] = "d/dt ICIT_c in component dICIT_c_dt (molar)" legend_rates[47] = "d/dt AKG_c in component dAKG_c_dt (molar)" legend_rates[48] = "d/dt SUC_c in component dSUC_c_dt (molar)" legend_rates[51] = "d/dt FUM_c in component dFUM_c_dt (molar)" legend_rates[46] = "d/dt MAL_c in component dMAL_c_dt (molar)" legend_rates[49] = "d/dt GLU_c in component dGLU_c_dt (molar)" legend_rates[50] = "d/dt ASP_c in component dASP_c_dt (molar)" legend_rates[55] = "d/dt PCr_c in component dPCr_c_dt (molar)" legend_rates[56] = "d/dt GLC_c in component dGLC_c_dt (molar)" legend_rates[57] = "d/dt G6P_c in component dG6P_c_dt (molar)" return (legend_states, legend_algebraic, legend_voi, legend_constants) def initConsts(): constants = [0.0] * sizeConstants; states = [0.0] * sizeStates; constants[0] = 2 constants[1] = 2 constants[2] = 1.22e-1 constants[3] = 11.6 constants[4] = 3.21e-2 constants[5] = 4.25e-1 constants[6] = 7.7e-2 constants[7] = 6.04e-7 constants[8] = 5.82e-1 constants[9] = 6.23e-2 constants[10] = 7.12e-3 constants[11] = 6.94e-2 constants[12] = 2.65e-2 constants[13] = 7.96 constants[14] = 4.12e8 constants[15] = 3.26e8 constants[16] = 7.31e1 constants[17] = 3.46e-1 constants[18] = 9.54e1 constants[19] = 1.58e1 constants[20] = 7.48e-5 constants[21] = 2.47e4 constants[22] = 6.65e-1 constants[23] = 2.81e-5 constants[24] = 3.14e-3 constants[25] = 9.93e-5 constants[26] = 5.95e3 constants[27] = 7.27e-3 constants[28] = 2.01e7 constants[29] = 1.01e-3 constants[30] = 5.65e6 constants[31] = 3.05e2 constants[32] = 0 constants[33] = 0.096484 constants[34] = 0.6514 constants[35] = 0.0724 constants[36] = 5.99 constants[37] = 2.725e5 constants[38] = 3 constants[39] = 327 constants[40] = 85 constants[41] = 85 constants[42] = 3.5e-6 constants[43] = 0.6 constants[44] = 1.2e-4 constants[45] = 2.97e-3 constants[46] = 1.35e-3 constants[47] = 2.7e-3 constants[48] = 10e-3 constants[49] = 42.7e-3 constants[50] = 21.4e-3 constants[51] = 1e-4 constants[52] = 6.75e-6 states[0] = 5.79e-8 states[1] = 1e-7 states[2] = -1.41e1 states[3] = 8.48e-21 states[4] = 1e-32 constants[53] = -109.7 states[5] = 2.99e-4 states[6] = 9.7e-8 constants[54] = 46.69 states[7] = 6.5e-5 constants[55] = -202.2 states[8] = 1.03e-11 states[9] = 1e-2 constants[56] = -4.51 states[10] = 1.21e-11 states[11] = 3.01e-3 states[12] = 2.14e-2 states[13] = 2.55e-12 constants[57] = 5.02e-4 constants[58] = 38.3e-6 constants[59] = 9.9e-6 constants[60] = 60.7e-6 constants[61] = 40.2e-6 constants[62] = 40e-6 constants[63] = 1e-32 states[14] = 8.92e-4 states[15] = 6.82e-5 states[16] = 1e-6 states[17] = 6.96e-10 constants[64] = 7.34e-8 constants[65] = 4e-6 constants[66] = 14e-6 constants[67] = 3.33e-6 constants[68] = 1600e-6 constants[69] = 900e-6 constants[70] = 1800e-6 constants[71] = 6000e-6 constants[72] = 67e-6 constants[73] = 140e-6 states[18] = 9.36e-7 constants[74] = 7.59e-2 constants[75] = 1161e-6 constants[76] = 434e-6 states[19] = 4.8e-11 constants[77] = 3.5e-16 constants[78] = 74e-6 constants[79] = 183e-6 constants[80] = 23.8e-6 constants[81] = 29e-6 constants[82] = 3 constants[83] = 91e-6 constants[84] = 50e-6 constants[85] = 1e-32 constants[86] = 6.93e-3 constants[87] = 80e-6 constants[88] = 55e-6 constants[89] = 21e-6 constants[90] = 6.9e-6 constants[91] = 6.036773936330217e-7 constants[92] = 1e3 constants[93] = 50e-6 constants[94] = 100e-6 constants[95] = 1e-32 states[20] = 5e-3 states[21] = 2.17e-7 states[22] = 2.72e-10 constants[96] = 9.54e-9 constants[97] = 5.5e-6 constants[98] = 100e-6 constants[99] = 2000e-6 constants[100] = 20e-6 constants[101] = 3000e-6 constants[102] = 11.1e-6 constants[103] = 16e-6 constants[104] = 55e-6 constants[105] = 660e-6 constants[106] = 20e-6 constants[107] = 880e-6 constants[108] = 11.1e-6 states[23] = 7.85e-9 constants[109] = 1.69 constants[110] = 120e-6 constants[111] = 1275e-6 constants[112] = 467e-6 constants[113] = 480e-6 constants[114] = 2.45e-6 constants[115] = 1200e-6 constants[116] = 1.5e-6 constants[117] = 450e-6 constants[118] = 375e-6 states[24] = 3.22e-8 constants[119] = 4.4 constants[120] = 44.7e-6 constants[121] = 197.7e-6 constants[122] = 3500e-6 constants[123] = 40e-6 constants[124] = 400e-6 constants[125] = 80e-6 constants[126] = 330e-6 constants[127] = 2.27e-12 constants[128] = 279e-6 constants[129] = 360e-6 constants[130] = 5.5e-6 constants[131] = 3.18e-6 constants[132] = 90.55e-6 constants[133] = 250e-6 constants[134] = 6.128e-6 constants[135] = 2.58e-6 constants[136] = 183.2e-6 constants[137] = 394.4e-6 constants[138] = 420e-6 constants[139] = 1 constants[140] = 170e-6 constants[141] = 143.6e-6 constants[142] = 146.6e-6 constants[143] = 156.5e-6 constants[144] = 111e-6 constants[145] = 100e-6 constants[146] = 260e-6 constants[147] = 278e-6 constants[148] = 650e-6 constants[149] = 1e-32 states[25] = 1.02e-4 states[26] = 9.93e-12 constants[150] = 1.77 constants[151] = 3480e-6 constants[152] = 710e-6 constants[153] = 50e-6 constants[154] = 8400e-6 constants[155] = 3900e-6 constants[156] = 430e-6 constants[157] = 88e-6 constants[158] = 8900e-6 constants[159] = 16.6e-3 constants[160] = 1e-32 states[27] = 9.18e-10 states[28] = 5.75e-12 states[29] = 1e-9 constants[161] = 0.3e-3 constants[162] = 0.17e-3 constants[163] = 1.4e-3 constants[164] = 0.7e-3 states[30] = 1e-9 states[31] = 1.16e-4 constants[165] = 0.028e-3 constants[166] = 2.8e-3 constants[167] = 0.18e-3 constants[168] = 1.6e-3 states[32] = 4.44e-9 constants[169] = 1.8 states[33] = 1e-9 states[34] = 1e-9 states[35] = 1e-9 states[36] = 1e-9 states[37] = 1e-9 states[38] = 1e-9 states[39] = 1e-9 states[40] = 1.16e-4 states[41] = 1e-9 states[42] = 1e-9 states[43] = 1e-9 states[44] = 1e-9 states[45] = 1e-9 states[46] = 1e-9 states[47] = 1e-9 states[48] = 1e-9 states[49] = 5.75e-12 states[50] = 4.44e-9 constants[170] = 0 states[51] = 1e-9 states[52] = 1e-9 constants[171] = 1e-9 constants[172] = 1e-9 states[53] = 150e-3 states[54] = 8.69e-2 constants[173] = 1e7 constants[174] = 0.13279 constants[175] = 1e7 constants[176] = 0.13279 constants[177] = 1e7 states[55] = 1e-9 states[56] = 1e-9 states[57] = 1e-9 constants[178] = 1e-3 constants[179] = 47e-6 constants[180] = 47e-6 constants[181] = 1e-3 constants[182] = 1e-3 constants[183] = 47e-6 constants[184] = 47e-6 constants[185] = 1e-3 constants[186] = 10e-6 states[58] = 4.64e-3 states[59] = 5e-3 constants[187] = 1 constants[188] = 1 constants[189] = 1 constants[190] = 1 constants[191] = 1 constants[192] = 1 constants[193] = 1 constants[194] = 1 constants[195] = 1 constants[196] = 1 constants[197] = 1 constants[198] = 1 constants[199] = 1 constants[200] = 1 constants[201] = 1 constants[202] = 1 constants[203] = 1 constants[204] = 1 constants[205] = 1 constants[206] = 1 constants[207] = 1 constants[208] = 1 constants[209] = 1 constants[210] = 1 constants[211] = 1 constants[212] = 1 constants[213] = 1 constants[214] = 1 constants[215] = 1 constants[216] = 1 constants[217] = 1 constants[218] = 1 constants[219] = 1 constants[220] = 1 constants[221] = 1 constants[222] = 1 constants[223] = 1 constants[224] = 1 constants[225] = 1 constants[226] = 1 constants[227] = 1 constants[228] = 1 constants[229] = 1 constants[230] = 1 constants[231] = 1 constants[232] = 1 constants[233] = 1 constants[234] = 1 constants[235] = 1 constants[236] = 1 constants[237] = 1 constants[238] = 0 constants[239] = 0 constants[240] = 0 constants[241] = 0 constants[242] = 0 constants[243] = 0.02 constants[244] = 1e-7 constants[245] = custom_piecewise([equal(constants[0] , 1.00000) & equal(constants[1] , 3.00000), 0.191000 , True, 0.00000]) constants[246] = 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.842508 , True, 80.0000]) constants[247] = (constants[12]*constants[147]*constants[142])/(constants[139]*constants[140]*constants[145]) constants[248] = 1.00000*(power(10.0000, -6.50000)) constants[249] = custom_piecewise([equal(constants[0] , 1.00000), 2.50370 , True, 2.57858]) constants[250] = 1.00000*(power(10.0000, -6.71000)) constants[251] = 1.00000*(power(10.0000, -1.69000)) constants[252] = 1.00000*(power(10.0000, 0.00740000)) constants[253] = 1.00000*(power(10.0000, -8.13000)) constants[254] = 1.00000*(power(10.0000, -0.00510000)) constants[255] = 1.00000*(power(10.0000, -5.63000)) constants[256] = 1.00000*(power(10.0000, -3.37000)) constants[257] = 1.00000*(power(10.0000, -0.339000)) constants[258] = 1.00000*(power(10.0000, -5.64000)) constants[259] = 1.00000*(power(10.0000, -2.46000)) constants[260] = 1.00000*(power(10.0000, -3.96000)) constants[261] = 1.00000*(power(10.0000, -5.13000)) constants[262] = 1.00000*(power(10.0000, -1.17000)) constants[263] = 1.00000*(power(10.0000, -0.503000)) constants[264] = 1.00000*(power(10.0000, -4.10000)) constants[265] = 1.00000*(power(10.0000, -4.75000)) constants[266] = 1.00000*(power(10.0000, -1.55000)) constants[267] = 1.00000*(power(10.0000, 0.170000)) constants[268] = 1.00000*(power(10.0000, -9.82000)) constants[269] = 1.00000*(power(10.0000, -1.02000)) constants[270] = 1.00000*(power(10.0000, -4.06000)) constants[271] = 1.00000*(power(10.0000, -1.82000)) constants[272] = 1.00000*(power(10.0000, -3.65000)) constants[273] = 1.00000*(power(10.0000, -2.32000)) constants[274] = 1.00000 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 constants[286] = constants[51]/2.00000 constants[287] = 0.00000 constants[288] = 0.00000 constants[289] = 0.00000 constants[290] = 0.00000 constants[321] = 0.00000 constants[322] = 0.00000 constants[323] = 0.00000 constants[324] = 0.00000 constants[325] = 0.00000 constants[326] = 0.00000 constants[291] = custom_piecewise([equal(constants[0] , 3.00000), 0.288200 , equal(constants[0] , 4.00000), 0.0560000 , True, 1.00000]) constants[292] = exp(-15.3800/constants[249])*1.00000 constants[293] = exp(50.7000/constants[249])/1.00000 constants[294] = constants[286] constants[295] = custom_piecewise([equal(constants[0] , 3.00000), 0.680100 , equal(constants[0] , 4.00000), 0.894000 , True, 1.00000]) constants[296] = 1.00000*(power(10.0000, -6.59000)) constants[297] = 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[291]/constants[295]]) constants[298] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) & constants[1] != 0.00000 | equal(constants[0] , 2.00000) & constants[1] != 0.00000 | equal(constants[0] , 6.00000) & constants[1] != 0.00000, 1.00000 , True, constants[291]/constants[295]]) constants[299] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) & constants[1] != 0.00000 | equal(constants[0] , 2.00000) & constants[1] != 0.00000 | equal(constants[0] , 6.00000) & constants[1] != 0.00000, 1.00000 , True, constants[291]/constants[295]]) constants[300] = 1.00000*(power(10.0000, -3.82000)) constants[301] = 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[291]/constants[295]]) constants[302] = 1.00000*(power(10.0000, -1.87000)) constants[303] = 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[291]/constants[295]]) constants[304] = 1.00000*(power(10.0000, -6.42000)) constants[305] = 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) & constants[1] != 0.00000 | equal(constants[0] , 2.00000) & constants[1] != 0.00000 | equal(constants[0] , 6.00000) & constants[1] != 0.00000, 1.00000 , True, constants[291]/constants[295]]) constants[306] = 1.00000*(power(10.0000, -2.79000)) constants[307] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) & constants[1] != 0.00000 | equal(constants[0] , 2.00000) & constants[1] != 0.00000 | equal(constants[0] , 6.00000) & constants[1] != 0.00000, 1.00000 , True, constants[291]/constants[295]]) constants[308] = 1.00000*(power(10.0000, -1.53000)) constants[309] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) & constants[1] != 0.00000 | equal(constants[0] , 2.00000) & constants[1] != 0.00000 | equal(constants[0] , 6.00000) & constants[1] != 0.00000, 1.00000 , True, constants[291]/constants[295]]) constants[310] = 1.00000*(power(10.0000, -6.22000)) constants[311] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) & constants[1] != 0.00000 | equal(constants[0] , 2.00000) & constants[1] != 0.00000 | equal(constants[0] , 6.00000) & constants[1] != 0.00000, 1.00000 , True, constants[291]/constants[295]]) constants[312] = 1.00000*(power(10.0000, -1.86000)) constants[313] = custom_piecewise([equal(constants[1] , 0.00000) | equal(constants[0] , 5.00000), 0.00000 , equal(constants[0] , 1.00000) & constants[1] != 0.00000 | equal(constants[0] , 2.00000) & constants[1] != 0.00000 | equal(constants[0] , 6.00000) & constants[1] != 0.00000, 1.00000 , True, constants[291]/constants[295]]) constants[314] = 1.00000*(power(10.0000, -1.05000)) constants[315] = 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[291]/constants[295]]) constants[316] = 1.00000*(power(10.0000, -5.91000)) constants[317] = 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[291]/constants[295]]) constants[318] = 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, 1.00000]) constants[319] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), 1.00000+1.00000/constants[246] , True, 0.00000]) constants[320] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), 1.00000+1.00000/constants[246] , True, 0.00000]) return (states, constants) def computeRates(voi, states, constants): rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic rates[16] = constants[321] rates[7] = constants[322] rates[12] = constants[323] rates[1] = constants[324] rates[59] = constants[325] rates[53] = constants[326] algebraic[1] = constants[36]*constants[41]*(states[44]-states[45]) rates[45] = algebraic[1]/constants[35] algebraic[2] = constants[36]*constants[170]*(states[51]-states[52]) rates[52] = algebraic[2]/constants[35] rates[44] = (-constants[298]*algebraic[1])/constants[246] rates[51] = (-constants[299]*algebraic[2])/constants[246] algebraic[13] = (constants[17]*(states[30]*states[24]-states[19]*states[29]))/(constants[161]*constants[164]*(2.00000+states[29]/constants[163]+states[24]/constants[164]+states[30]/constants[161]+states[19]/constants[162]+(states[29]*states[19])/(constants[163]*constants[162])+(states[24]*states[30])/(constants[164]*constants[161]))) algebraic[15] = constants[36]*constants[41]*(states[47]-states[30]) rates[30] = (-algebraic[13]+algebraic[15])/constants[35] rates[47] = (-constants[309]*algebraic[15])/constants[246] algebraic[18] = constants[36]*constants[39]*(states[40]-states[31]) rates[40] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), (-constants[303]*algebraic[18])/constants[246] , True, ((-constants[291]/constants[295])*algebraic[18]+constants[32])/constants[246]]) algebraic[20] = constants[36]*constants[41]*(states[41]-states[27]) rates[41] = (-constants[305]*algebraic[20])/constants[246] algebraic[22] = constants[36]*constants[41]*(states[42]-states[43]) rates[42] = (-constants[307]*algebraic[22])/constants[246] algebraic[24] = constants[36]*constants[41]*(states[46]-states[29]) rates[46] = (-constants[313]*algebraic[24])/constants[246] algebraic[26] = constants[36]*constants[41]*(states[48]-states[33]) rates[48] = (-constants[311]*algebraic[26])/constants[246] algebraic[28] = constants[36]*constants[41]*(states[49]-states[28]) rates[49] = (-constants[315]*algebraic[28])/constants[246] algebraic[30] = constants[36]*constants[41]*(states[50]-states[32]) rates[50] = (-constants[317]*algebraic[30])/constants[246] algebraic[19] = states[1] algebraic[23] = states[59] algebraic[25] = states[53] algebraic[27] = 1.00000+algebraic[19]/constants[296]+algebraic[23]/constants[300]+algebraic[25]/constants[302] algebraic[29] = ((states[34]*algebraic[23])/constants[300])/algebraic[27] algebraic[31] = 1.00000+algebraic[19]/constants[304]+algebraic[23]/constants[306]+algebraic[25]/constants[308] algebraic[37] = ((states[36]*algebraic[23])/constants[306])/algebraic[31] algebraic[21] = 1.00000+algebraic[19]/constants[316] algebraic[33] = (constants[292]*algebraic[21]*algebraic[31])/(algebraic[19]*algebraic[27]) algebraic[35] = (algebraic[33]*constants[300]*algebraic[27])/(constants[306]*algebraic[31]) algebraic[39] = ((constants[245]/(constants[179]*constants[182]))*(algebraic[29]*states[56]-(states[57]*algebraic[37])/algebraic[35]))/(1.00000+algebraic[29]/constants[178]+states[56]/constants[179]+(algebraic[29]*states[56])/(constants[179]*constants[182])+states[57]/constants[180]+algebraic[37]/constants[181]+(states[57]*algebraic[37])/(constants[181]*constants[184])+(states[57]*states[56])/(constants[186]*constants[179])) rates[56] = (-constants[319]*algebraic[39])/constants[246] rates[57] = (constants[320]*algebraic[39])/constants[246] algebraic[17] = constants[36]*constants[40]*(states[38]-states[39]) algebraic[42] = 1.00000+algebraic[19]/constants[310]+algebraic[23]/constants[312]+algebraic[25]/constants[314] algebraic[44] = (constants[176]*algebraic[27]*algebraic[42])/(algebraic[31]*algebraic[31]) algebraic[47] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[175]*(algebraic[44]*states[36]*states[36]-states[38]*states[34])]) rates[38] = 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[291]/constants[295])*algebraic[17]+algebraic[47])/constants[246]]) algebraic[14] = constants[36]*constants[40]*(states[34]-states[35]) algebraic[49] = constants[49]-states[55] algebraic[51] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[177]*(((constants[293]*states[36])/algebraic[31])*states[55]*algebraic[19]-(states[34]/algebraic[27])*algebraic[49])]) rates[34] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), (-(1.00000+1.00000/constants[246])*algebraic[39]-constants[297]*algebraic[14])/constants[246] , True, ((-constants[297]*algebraic[14]+algebraic[51]+algebraic[47])-constants[32])/constants[246]]) algebraic[16] = constants[36]*constants[40]*(states[36]-states[37]) rates[36] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000) | equal(constants[0] , 6.00000), ((1.00000+1.00000/constants[246])*algebraic[39]-constants[301]*algebraic[16])/constants[246] , True, (((-constants[301]*algebraic[16]-algebraic[51])-algebraic[47])+constants[32])/constants[246]]) rates[55] = (-constants[318]*algebraic[51])/constants[246] algebraic[36] = 1.00000+states[0]/constants[296]+states[58]/constants[300]+states[54]/constants[302] algebraic[41] = 1.00000+states[0]/constants[304]+states[58]/constants[306]+states[54]/constants[308] algebraic[60] = 1.00000+states[1]/constants[304]+states[59]/constants[306]+states[53]/constants[308] algebraic[59] = 1.00000+states[1]/constants[296]+states[59]/constants[300]+states[53]/constants[302] algebraic[63] = custom_piecewise([greater(states[37] , constants[171]) & greater(states[8] , constants[171]), (constants[27]/(1.00000+(constants[42]/states[37])/algebraic[60]))*((states[37]/algebraic[60])/(states[37]/algebraic[60]+(states[35]/algebraic[59])*exp(((constants[43]-1.00000)*constants[33]*states[2])/constants[249]))+(states[9]/algebraic[41])/(states[9]/algebraic[41]+(states[8]/algebraic[36])*exp((constants[43]*constants[33]*states[2])/constants[249]))) , True, 0.00000]) algebraic[67] = 1.00000+states[1]/constants[310]+states[59]/constants[312]+states[53]/constants[314] algebraic[69] = (constants[174]*algebraic[59]*algebraic[67])/(algebraic[60]*algebraic[60]) algebraic[73] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[173]*(algebraic[69]*states[37]*states[37]-states[39]*states[35])]) rates[35] = (algebraic[14]+algebraic[63]+algebraic[73])/constants[35] rates[37] = ((algebraic[16]-algebraic[63])-2.00000*algebraic[73])/constants[35] rates[39] = (algebraic[17]+algebraic[73])/constants[35] algebraic[76] = 1.00000+states[0]/constants[250]+states[58]/constants[251]+states[54]/constants[252] algebraic[8] = exp(-(constants[56]-constants[38]*constants[33]*states[2])/constants[249]) algebraic[77] = (algebraic[8]*(power(states[1], constants[38]))*algebraic[36])/((power(states[0], constants[38]-1.00000))*algebraic[41]*algebraic[76]*1.00000) algebraic[78] = constants[26]*(algebraic[77]*states[9]*states[5]-states[8]) algebraic[50] = 1.00000+states[0]/constants[310]+states[58]/constants[312]+states[54]/constants[314] algebraic[52] = (states[16]*(1.00000+states[0]/constants[310]))/algebraic[50] algebraic[56] = 1.00000+algebraic[52]/constants[148] algebraic[57] = custom_piecewise([greater(states[22] , constants[149]) & greater(states[9] , constants[149]), (((constants[12]*constants[247])/algebraic[56])*states[22]*states[9])/(constants[247]*constants[145]*states[22]+constants[247]*constants[144]*states[9]+(constants[12]*constants[147]*states[20])/constants[139]+(constants[12]*constants[146]*states[8])/constants[139]+constants[247]*states[22]*states[9]+(constants[12]*constants[147]*states[22]*states[20])/(constants[139]*constants[140])+(constants[12]*states[20]*states[8])/constants[139]+(constants[247]*constants[144]*states[9]*states[8])/constants[143]) , True, 0.00000]) algebraic[58] = custom_piecewise([greater(states[20] , constants[149]) & greater(states[8] , constants[149]), algebraic[57]-((((constants[12]*constants[247])/algebraic[56])*states[20]*states[8])/constants[139])/(constants[247]*constants[145]*states[22]+constants[247]*constants[144]*states[9]+(constants[12]*constants[147]*states[20])/constants[139]+(constants[12]*constants[146]*states[8])/constants[139]+constants[247]*states[22]*states[9]+(constants[12]*constants[147]*states[22]*states[20])/(constants[139]*constants[140])+(constants[12]*states[20]*states[8])/constants[139]+(constants[247]*constants[144]*states[9]*states[8])/constants[143]) , True, algebraic[57]]) rates[8] = ((algebraic[58]+algebraic[78])-algebraic[63])/constants[34] rates[9] = ((-algebraic[58]-algebraic[78])+algebraic[63])/constants[34] algebraic[83] = 1.00000+states[0]/constants[253] algebraic[98] = 1.00000+states[0]/constants[260] algebraic[99] = 1.00000+states[0]/constants[261]+states[58]/constants[262]+states[54]/constants[263] algebraic[66] = algebraic[36] algebraic[72] = algebraic[41] algebraic[100] = (constants[96]*algebraic[83]*algebraic[99]*algebraic[66])/(states[0]*algebraic[72]*algebraic[98]*algebraic[76]) algebraic[101] = (constants[8]*constants[106]*constants[101]*constants[102])/(algebraic[100]*constants[97]*constants[98]*constants[105]) algebraic[102] = (constants[8]*algebraic[101]*(states[20]*states[17]*states[5]-(states[11]*states[21]*states[22])/algebraic[100]))/(algebraic[101]*constants[97]*constants[98]*constants[106]+algebraic[101]*constants[98]*constants[106]*states[20]+algebraic[101]*constants[97]*constants[104]*states[5]+algebraic[101]*constants[106]*states[20]*states[17]+algebraic[101]*constants[104]*states[20]*states[5]+algebraic[101]*constants[103]*states[17]*states[5]+algebraic[101]*states[20]*states[17]*states[5]+(constants[8]*constants[102]*constants[107]*states[11])/algebraic[100]+(constants[8]*constants[101]*constants[106]*states[22])/algebraic[100]+(constants[8]*constants[108]*states[11]*states[21])/algebraic[100]+(constants[8]*constants[107]*states[11]*states[22])/algebraic[100]+(constants[8]*constants[106]*states[21]*states[22])/algebraic[100]+(constants[8]*states[11]*states[21]*states[22])/algebraic[100]+(constants[8]*constants[107]*constants[102]*states[20]*states[11])/(constants[97]*algebraic[100])+(algebraic[101]*constants[97]*constants[104]*states[5]*states[22])/constants[102]+(constants[8]*constants[107]*constants[102]*states[20]*states[17]*states[21])/(constants[97]*constants[98]*algebraic[100])+(algebraic[101]*constants[103]*states[17]*states[5]*states[22])/constants[102]+(constants[8]*constants[108]*states[20]*states[11]*states[21])/(constants[97]*algebraic[100])+(algebraic[101]*constants[97]*constants[104]*states[5]*states[21]*states[22])/(constants[101]*constants[102])+(constants[8]*constants[102]*constants[107]*states[20]*states[17]*states[5]*states[11])/(constants[97]*constants[98]*constants[99]*algebraic[100])+(constants[8]*constants[100]*constants[108]*states[20]*states[17]*states[5]*states[21])/(constants[97]*constants[98]*constants[100]*algebraic[100])+(constants[8]*constants[108]*states[20]*states[17]*states[11]*states[21])/(constants[97]*constants[98]*algebraic[100])+(constants[8]*constants[103]*states[17]*states[5]*states[21]*states[22])/(constants[101]*constants[102])+(algebraic[101]*constants[103]*constants[100]*states[17]*states[11]*states[21]*states[22])/(constants[100]*constants[101]*constants[102])+(algebraic[101]*constants[97]*constants[104]*states[5]*states[11]*states[21]*states[22])/(constants[100]*constants[101]*constants[102])+(constants[8]*constants[108]*states[20]*states[17]*states[5]*states[11]*states[21])/(constants[97]*constants[98]*constants[100]*algebraic[100])+(constants[8]*constants[103]*states[17]*states[5]*states[11]*states[21]*states[22])/(constants[100]*constants[101]*constants[102])) rates[22] = (algebraic[102]-algebraic[58])/constants[34] rates[20] = (-algebraic[102]+algebraic[58])/constants[34] algebraic[85] = 1.00000+states[0]/constants[255]+states[58]/constants[256]+states[54]/constants[257] algebraic[84] = 1.00000+states[58]/constants[254] algebraic[86] = (constants[64]*algebraic[85])/(states[0]*states[0]*algebraic[84]) algebraic[87] = (states[15]*(1.00000+states[0]/constants[255]))/algebraic[85] algebraic[88] = 1.00000+algebraic[87]/constants[68] algebraic[38] = (states[8]*(1.00000+states[0]/constants[296]))/algebraic[36] algebraic[43] = (states[9]*(1.00000+states[0]/constants[304]))/algebraic[41] algebraic[54] = 1.00000+algebraic[38]/constants[69]+algebraic[43]/constants[70]+algebraic[52]/constants[71]+states[11]/constants[72]+states[17]/constants[73] algebraic[89] = (constants[3]*(states[13]*states[14]-(states[11]*states[15])/algebraic[86]))/(constants[67]*constants[66]*algebraic[88]+constants[65]*algebraic[88]*states[13]+constants[66]*algebraic[54]*states[14]+states[14]*states[13]) algebraic[94] = 1.00000+states[0]/constants[258]+states[58]/constants[259] algebraic[95] = (constants[74]*algebraic[94])/algebraic[85] algebraic[96] = (constants[4]*constants[76])/(constants[75]*algebraic[95]) algebraic[97] = (constants[4]*algebraic[96]*(states[15]-states[18]/algebraic[95]))/(constants[75]*algebraic[96]+algebraic[96]*states[15]+(constants[4]*states[18])/algebraic[95]) algebraic[107] = 1.00000+states[0]/constants[265]+states[58]/constants[266]+states[54]/constants[267] algebraic[91] = 1.00000+states[1]/constants[255]+states[59]/constants[256]+states[53]/constants[257] algebraic[92] = ((states[43]*states[1])/constants[255])/algebraic[91] algebraic[113] = 1.00000+states[1]/constants[265]+states[59]/constants[266]+states[53]/constants[267] algebraic[90] = ((states[15]*states[0])/constants[255])/algebraic[85] algebraic[114] = constants[16]*((algebraic[92]*states[24])/algebraic[107]-(algebraic[90]*states[29])/algebraic[113]) rates[15] = (-algebraic[97]+algebraic[114]+algebraic[89])/constants[34] rates[43] = (-algebraic[114]+algebraic[22])/constants[35] algebraic[80] = 1.00000+states[1]/constants[250]+states[59]/constants[251]+states[53]/constants[252] algebraic[115] = constants[19]*(((states[29]/algebraic[113])*states[5])/algebraic[76]-((states[24]/algebraic[107])*states[31])/algebraic[80]) algebraic[81] = ((states[31]*states[1])/constants[250])/algebraic[80] algebraic[79] = ((states[5]*states[0])/constants[250])/algebraic[76] algebraic[82] = (constants[28]*(algebraic[81]*states[1]-algebraic[79]*states[0]))/(constants[29]*(1.00000+algebraic[81]/constants[29])*(1.00000+algebraic[79]/constants[29])) rates[5] = (((-algebraic[102]-algebraic[78])+algebraic[82])-algebraic[115])/constants[34] rates[31] = (-algebraic[82]+algebraic[18]+algebraic[115])/constants[35] algebraic[103] = 1.00000+states[1]/constants[261]+states[59]/constants[262]+states[53]/constants[263] algebraic[116] = constants[18]*(((states[33]/algebraic[103])*states[24])/algebraic[107]-((states[21]/algebraic[99])*states[29])/algebraic[113]) rates[33] = (algebraic[26]-algebraic[116])/constants[35] rates[29] = (-algebraic[115]+algebraic[24]+algebraic[13]+algebraic[114]+algebraic[116])/constants[35] algebraic[120] = 1.00000+states[58]/constants[269] algebraic[122] = 1.00000+states[59]/constants[269] algebraic[123] = constants[14]*((states[27]/algebraic[122])*states[1]-(states[10]/algebraic[120])*states[0]) rates[27] = (-algebraic[123]+algebraic[20])/constants[35] algebraic[127] = 1.00000+states[0]/constants[272]+states[58]/constants[273] algebraic[124] = 1.00000+states[0]/constants[270]+states[58]/constants[271] algebraic[128] = (constants[150]*algebraic[124]*algebraic[84])/algebraic[127] algebraic[129] = (constants[13]*constants[158]*constants[153])/(algebraic[128]*constants[151]*constants[156]) algebraic[12] = 1.00000+states[19]/constants[159] algebraic[130] = custom_piecewise([greater(states[25] , constants[160]) & greater(states[19] , constants[160]), (constants[13]*algebraic[129]*states[25]*states[19])/(algebraic[129]*constants[156]*states[25]+algebraic[129]*constants[155]*algebraic[12]*states[19]+(constants[13]*constants[158]*algebraic[12]*states[14])/algebraic[128]+(constants[13]*constants[157]*states[26])/algebraic[128]+algebraic[129]*states[25]*states[19]+(constants[13]*constants[158]*states[25]*states[14])/(algebraic[128]*constants[151])+(constants[13]*states[14]*states[26])/algebraic[128]+(algebraic[129]*constants[155]*states[19]*states[26])/constants[154]) , True, 0.00000]) algebraic[131] = custom_piecewise([greater(states[14] , constants[160]) & greater(states[26] , constants[160]), algebraic[130]-((constants[13]*algebraic[129]*states[14]*states[26])/algebraic[128])/(algebraic[129]*constants[156]*states[25]+algebraic[129]*constants[155]*algebraic[12]*states[19]+(constants[13]*constants[158]*algebraic[12]*states[14])/algebraic[128]+(constants[13]*constants[157]*states[26])/algebraic[128]+algebraic[129]*states[25]*states[19]+(constants[13]*constants[158]*states[25]*states[14])/(algebraic[128]*constants[151])+(constants[13]*states[14]*states[26])/algebraic[128]+(algebraic[129]*constants[155]*states[19]*states[26])/constants[154]) , True, algebraic[130]]) algebraic[125] = 1.00000+states[1]/constants[270]+states[59]/constants[271] algebraic[126] = constants[15]*((states[28]/algebraic[125])*states[1]-(states[26]/algebraic[124])*states[0]) algebraic[132] = 1.00000+states[1]/constants[272]+states[59]/constants[273] algebraic[133] = (exp((-constants[33]*states[2])/constants[249])*algebraic[127]*algebraic[125])/(algebraic[132]*algebraic[124]) algebraic[134] = (constants[20]*(algebraic[133]*states[32]*states[26]*states[0]-states[25]*states[28]*states[1]))/(algebraic[133]*constants[165]*constants[168]*constants[248]*(2.00000*constants[169]+(constants[169]*states[32])/constants[165]+(states[32]*states[26]*states[0])/(constants[165]*constants[168]*constants[248])+(constants[169]*states[25]*states[1])/(constants[166]*constants[248])+(states[25]*states[28]*states[1])/(constants[166]*constants[167]*constants[248])+(constants[169]*states[25])/constants[166]+(constants[169]*states[32]*states[0])/(constants[165]*constants[248])+(constants[169]*states[0])/constants[248]+(constants[169]*states[28]*states[1])/(constants[167]*constants[248])+(constants[169]*states[1])/constants[248]+(constants[169]*states[26]*states[0])/(constants[168]*constants[248]))) rates[26] = ((algebraic[131]+algebraic[126])-algebraic[134])/constants[34] rates[25] = (-algebraic[131]+algebraic[134])/constants[34] rates[28] = (-algebraic[126]+algebraic[134]+algebraic[28])/constants[35] rates[32] = (-algebraic[134]+algebraic[30])/constants[35] algebraic[143] = custom_piecewise([greater(states[3] , constants[45]), 0.00000 , True, constants[45]-states[3]]) algebraic[117] = 1.00000+states[0]/constants[268] algebraic[121] = (constants[57]*algebraic[117])/(states[0]*algebraic[120]*algebraic[83]) algebraic[9] = 1.00000+states[13]/constants[61] algebraic[10] = 1.00000+states[3]/constants[62] algebraic[144] = custom_piecewise([greater(states[10] , constants[63]) & greater(states[11] , constants[63]) & greater(algebraic[143] , constants[63]), (constants[2]*(states[10]*states[11]*algebraic[143]-(states[12]*states[13]*states[3])/algebraic[121]))/(constants[60]*algebraic[10]*states[10]*states[11]+constants[59]*algebraic[9]*states[10]*algebraic[143]+constants[58]*states[11]*algebraic[143]+states[10]*states[11]*algebraic[143]) , True, 0.00000]) rates[10] = (-algebraic[144]+algebraic[123])/constants[34] rates[13] = (algebraic[144]-algebraic[89])/constants[34] algebraic[118] = (constants[77]*algebraic[117])/(states[0]*states[0]*algebraic[94]) algebraic[46] = 1.00000+(constants[84]/algebraic[43])*(1.00000+algebraic[38]/constants[83]) algebraic[145] = custom_piecewise([greater(algebraic[143] , constants[85]) & greater(states[18] , constants[85]), (constants[5]*(1.00000-(states[12]*states[19]*states[3])/(algebraic[118]*algebraic[143]*states[18])))/(1.00000+(power(constants[79]/states[18], constants[82]))*algebraic[46]+(constants[78]/algebraic[143])*(1.00000+(power(constants[80]/states[18], constants[82]))*algebraic[46]+(states[3]*algebraic[46])/constants[81])) , True, 0.00000]) rates[18] = (algebraic[97]-algebraic[145])/constants[34] algebraic[119] = (constants[86]*algebraic[117]*algebraic[98])/(states[0]*algebraic[83]) algebraic[48] = 1.00000+(constants[94]/algebraic[43])*(1.00000+algebraic[38]/constants[93]) algebraic[146] = custom_piecewise([greater(states[19] , constants[95]) & greater(states[11] , constants[95]) & greater(algebraic[143] , constants[95]), ((constants[6]*(1.00000-(states[12]*states[17]*states[3])/(algebraic[119]*states[19]*states[11]*algebraic[143])))/(1.00000+(constants[87]/states[19])*algebraic[48]+(constants[88]/states[11])*(1.00000+states[17]/constants[90])+(constants[89]/algebraic[143])*(1.00000+states[3]/constants[91])))/(1.00000+states[3]/constants[92]) , True, 0.00000]) rates[19] = (((algebraic[145]-algebraic[146])-algebraic[131])+algebraic[13])/constants[34] rates[17] = (algebraic[146]-algebraic[102])/constants[34] rates[11] = ((-algebraic[144]-algebraic[146])+algebraic[102]+algebraic[89])/constants[34] algebraic[104] = 1.00000+states[0]/constants[264] algebraic[108] = (constants[119]*algebraic[107])/algebraic[104] algebraic[109] = (constants[10]*constants[121])/(algebraic[108]*constants[120]) algebraic[68] = algebraic[38] algebraic[74] = algebraic[43] algebraic[75] = 1.00000+states[15]/constants[122]+algebraic[38]/constants[123]+algebraic[43]/constants[124]+algebraic[68]/constants[125]+algebraic[74]/constants[126] algebraic[110] = (algebraic[109]*constants[10]*(states[23]-states[24]/algebraic[108]))/(constants[120]*algebraic[109]*algebraic[75]+algebraic[109]*states[23]+(constants[10]*states[24])/algebraic[108]) algebraic[111] = (constants[127]*algebraic[84])/(states[0]*algebraic[107]) algebraic[112] = (constants[11]*constants[131]*constants[134])/(algebraic[111]*constants[128]*constants[133]) algebraic[55] = 1.00000+algebraic[38]/constants[136]+algebraic[43]/constants[137]+algebraic[52]/constants[138] algebraic[147] = (constants[11]*algebraic[112]*(algebraic[143]*states[24]-(states[14]*states[3])/algebraic[111]))/(algebraic[112]*constants[128]*constants[133]*algebraic[55]+algebraic[112]*constants[133]*algebraic[143]+algebraic[112]*constants[132]*algebraic[55]*states[24]+(constants[11]*constants[135]*algebraic[55]*states[14])/algebraic[111]+(constants[11]*constants[134]*states[3])/algebraic[111]+algebraic[112]*algebraic[143]*states[24]+(constants[11]*constants[135]*algebraic[143]*states[14])/(algebraic[111]*constants[128])+(constants[11]*states[14]*states[3])/algebraic[111]+(algebraic[112]*constants[132]*states[24]*states[3])/constants[131]+(algebraic[112]*algebraic[143]*states[24]*states[14])/constants[130]+(constants[11]*states[24]*states[14]*states[3])/(constants[129]*algebraic[111])) rates[24] = (((((algebraic[110]-algebraic[147])+algebraic[115])-algebraic[13])-algebraic[114])-algebraic[116])/constants[34] rates[14] = (-algebraic[89]+algebraic[147]+algebraic[131])/constants[34] algebraic[0] = exp(-(constants[53]+4.00000*constants[33]*states[2])/constants[249]) algebraic[3] = ((algebraic[0]*(power(states[0], 5.00000)))/(power(states[1], 4.00000)))/1.00000 algebraic[156] = constants[46]-states[4] algebraic[157] = constants[21]*(algebraic[3]*states[3]*algebraic[156]-algebraic[143]*states[4]) rates[3] = ((algebraic[144]+algebraic[145]+algebraic[146]+algebraic[147])-algebraic[157])/constants[34] algebraic[105] = (constants[109]*algebraic[104])/algebraic[99] algebraic[106] = (constants[9]*constants[114]*constants[111])/(algebraic[105]*constants[110]*constants[113]) algebraic[11] = (1.00000+states[14]/constants[116]+states[21]/constants[117]+states[23]/constants[118])/(1.00000+states[21]/constants[117]+states[23]/constants[118]) algebraic[158] = (constants[9]*algebraic[106]*(states[21]*algebraic[156]-(states[4]*states[23])/algebraic[105]))/(algebraic[106]*constants[110]*constants[113]*algebraic[11]+algebraic[106]*constants[113]*states[21]+algebraic[106]*constants[112]*algebraic[11]*algebraic[156]+(constants[9]*constants[115]*algebraic[11]*states[4])/algebraic[105]+(constants[9]*constants[114]*states[23])/algebraic[105]+algebraic[106]*states[21]*algebraic[156]+(constants[9]*constants[115]*states[21]*states[4])/(algebraic[105]*constants[110])+(algebraic[106]*constants[112]*algebraic[156]*states[23])/constants[111]+(constants[9]*states[4]*states[23])/algebraic[105]) rates[21] = ((algebraic[102]-algebraic[158])+algebraic[116])/constants[34] rates[23] = (algebraic[158]-algebraic[110])/constants[34] algebraic[4] = exp(-(constants[54]+2.00000*constants[33]*states[2])/constants[249]) algebraic[5] = (1.00000*1.00000*algebraic[4]*(power(states[0], 2.00000)))/(power(states[1], 4.00000)) algebraic[163] = constants[47]-states[6] algebraic[164] = ((constants[22]*(1.00000+states[5]/constants[23]))/(1.00000+states[5]/constants[24]))*((power(algebraic[5], 0.500000))*algebraic[163]*(power(states[4], 0.500000))-states[6]*(power(algebraic[156], 0.500000))) rates[4] = ((algebraic[158]+algebraic[157])-algebraic[164])/constants[34] algebraic[6] = exp(-(constants[55]+4.00000*constants[33]*states[2])/constants[249]) algebraic[7] = ((algebraic[6]*(power(states[0], 4.00000)))/(power(states[1], 2.00000)))/(1.00000*1.00000) algebraic[165] = ((((constants[25]*states[7])/(states[7]+constants[44]))*exp((constants[33]*states[2])/constants[249])*states[6])/constants[47])*((power(algebraic[7], 1.00000/2.00000))*states[6]*(power(states[7], 1.00000/4.00000))-algebraic[163]) algebraic[32] = custom_piecewise([greater(fabs(states[2]) , constants[172]), (constants[31]*states[2]*(states[1]*exp((constants[33]*states[2])/constants[249])-states[0]))/(exp((constants[33]*states[2])/constants[249])-1.00000) , True, (constants[31]*constants[249]*(states[1]-states[0]))/constants[33]]) rates[2] = (((((4.00000*algebraic[157]+2.00000*algebraic[164]+4.00000*algebraic[165])-constants[38]*algebraic[78])-algebraic[63])-algebraic[32])+algebraic[134])/constants[52] rates[6] = (2.00000*algebraic[164]-2.00000*algebraic[165])/constants[35] algebraic[168] = -((states[0]*states[0])/(constants[229]*constants[230]*(power(constants[283], 2.00000)))+(states[0]*states[8])/(constants[296]*constants[300]*(power(algebraic[36], 2.00000)))+(states[0]*states[9])/(constants[304]*constants[306]*(power(algebraic[41], 2.00000)))+(states[0]*states[16])/(constants[310]*constants[312]*(power(algebraic[50], 2.00000)))+(states[0]*states[22])/(constants[296]*constants[300]*(power(algebraic[66], 2.00000)))+(states[0]*states[20])/(constants[304]*constants[306]*(power(algebraic[72], 2.00000)))+(states[0]*states[5])/(constants[250]*constants[251]*(power(algebraic[76], 2.00000)))+(states[0]*states[3])/(constants[205]*constants[206]*(power(constants[275], 2.00000)))+(states[0]*algebraic[143])/(constants[202]*constants[203]*(power(constants[274], 2.00000)))+(states[0]*states[4])/(constants[208]*constants[209]*(power(constants[276], 2.00000)))+(states[0]*algebraic[156])/(constants[211]*constants[212]*(power(constants[277], 2.00000)))+(states[0]*states[10])/(constants[198]*constants[269]*(power(algebraic[120], 2.00000)))+(states[0]*states[14])/(constants[189]*constants[254]*(power(algebraic[84], 2.00000)))+(states[0]*states[13])/(constants[214]*constants[215]*(power(constants[278], 2.00000)))+(states[0]*states[15])/(constants[255]*constants[256]*(power(algebraic[85], 2.00000)))+(states[0]*states[18])/(constants[258]*constants[259]*(power(algebraic[94], 2.00000)))+(states[0]*states[19])/(constants[217]*constants[218]*(power(constants[279], 2.00000)))+(states[0]*states[17])/(constants[260]*constants[192]*(power(algebraic[98], 2.00000)))+(states[0]*states[11])/(constants[253]*constants[187]*(power(algebraic[83], 2.00000)))+(states[0]*states[21])/(constants[261]*constants[262]*(power(algebraic[99], 2.00000)))+(states[0]*states[23])/(constants[264]*constants[194]*(power(algebraic[104], 2.00000)))+(states[0]*states[24])/(constants[265]*constants[266]*(power(algebraic[107], 2.00000)))+(states[0]*states[26])/(constants[270]*constants[271]*(power(algebraic[124], 2.00000)))+(states[0]*states[25])/(constants[272]*constants[273]*(power(algebraic[127], 2.00000)))+(states[0]*states[54])/(constants[235]*constants[236]*(power(constants[285], 2.00000)))+(states[0]*states[58])/(constants[232]*constants[233]*(power(constants[284], 2.00000)))+(states[0]*states[7])/(constants[220]*constants[221]*(power(constants[280], 2.00000)))+(states[0]*constants[294])/(constants[223]*constants[224]*(power(constants[281], 2.00000)))+(states[0]*constants[286])/(constants[226]*constants[227]*(power(constants[282], 2.00000)))+(states[0]*states[12])/(constants[268]*constants[196]*(power(algebraic[117], 2.00000)))) algebraic[169] = -((states[0]*states[0])/(constants[229]*constants[231]*(power(constants[283], 2.00000)))+(states[0]*states[8])/(constants[296]*constants[302]*(power(algebraic[36], 2.00000)))+(states[0]*states[9])/(constants[304]*constants[308]*(power(algebraic[41], 2.00000)))+(states[0]*states[16])/(constants[310]*constants[314]*(power(algebraic[50], 2.00000)))+(states[0]*states[22])/(constants[296]*constants[302]*(power(algebraic[66], 2.00000)))+(states[0]*states[20])/(constants[304]*constants[308]*(power(algebraic[72], 2.00000)))+(states[0]*states[5])/(constants[250]*constants[252]*(power(algebraic[76], 2.00000)))+(states[0]*states[3])/(constants[205]*constants[207]*(power(constants[275], 2.00000)))+(states[0]*algebraic[143])/(constants[202]*constants[204]*(power(constants[274], 2.00000)))+(states[0]*states[4])/(constants[208]*constants[210]*(power(constants[276], 2.00000)))+(states[0]*algebraic[156])/(constants[211]*constants[213]*(power(constants[277], 2.00000)))+(states[0]*states[10])/(constants[198]*constants[199]*(power(algebraic[120], 2.00000)))+(states[0]*states[14])/(constants[189]*constants[190]*(power(algebraic[84], 2.00000)))+(states[0]*states[13])/(constants[214]*constants[216]*(power(constants[278], 2.00000)))+(states[0]*states[15])/(constants[255]*constants[257]*(power(algebraic[85], 2.00000)))+(states[0]*states[18])/(constants[258]*constants[191]*(power(algebraic[94], 2.00000)))+(states[0]*states[19])/(constants[217]*constants[219]*(power(constants[279], 2.00000)))+(states[0]*states[17])/(constants[260]*constants[193]*(power(algebraic[98], 2.00000)))+(states[0]*states[11])/(constants[253]*constants[188]*(power(algebraic[83], 2.00000)))+(states[0]*states[21])/(constants[261]*constants[263]*(power(algebraic[99], 2.00000)))+(states[0]*states[23])/(constants[264]*constants[195]*(power(algebraic[104], 2.00000)))+(states[0]*states[24])/(constants[265]*constants[267]*(power(algebraic[107], 2.00000)))+(states[0]*states[26])/(constants[270]*constants[200]*(power(algebraic[124], 2.00000)))+(states[0]*states[25])/(constants[272]*constants[201]*(power(algebraic[127], 2.00000)))+(states[0]*states[54])/(constants[235]*constants[237]*(power(constants[285], 2.00000)))+(states[0]*states[58])/(constants[232]*constants[234]*(power(constants[284], 2.00000)))+(states[0]*states[7])/(constants[220]*constants[222]*(power(constants[280], 2.00000)))+(states[0]*constants[294])/(constants[223]*constants[225]*(power(constants[281], 2.00000)))+(states[0]*constants[286])/(constants[226]*constants[228]*(power(constants[282], 2.00000)))+(states[0]*states[12])/(constants[268]*constants[197]*(power(algebraic[117], 2.00000)))) algebraic[172] = -((states[58]*states[0])/(constants[230]*constants[231]*(power(constants[283], 2.00000)))+(states[58]*states[8])/(constants[300]*constants[302]*(power(algebraic[36], 2.00000)))+(states[58]*states[9])/(constants[306]*constants[308]*(power(algebraic[41], 2.00000)))+(states[58]*states[16])/(constants[312]*constants[314]*(power(algebraic[50], 2.00000)))+(states[58]*states[22])/(constants[300]*constants[302]*(power(algebraic[66], 2.00000)))+(states[58]*states[20])/(constants[306]*constants[308]*(power(algebraic[72], 2.00000)))+(states[58]*states[5])/(constants[251]*constants[252]*(power(algebraic[76], 2.00000)))+(states[58]*states[3])/(constants[206]*constants[207]*(power(constants[275], 2.00000)))+(states[58]*algebraic[143])/(constants[203]*constants[204]*(power(constants[274], 2.00000)))+(states[58]*states[4])/(constants[209]*constants[210]*(power(constants[276], 2.00000)))+(states[58]*algebraic[156])/(constants[212]*constants[213]*(power(constants[277], 2.00000)))+(states[58]*states[10])/(constants[269]*constants[199]*(power(algebraic[120], 2.00000)))+(states[58]*states[14])/(constants[254]*constants[190]*(power(algebraic[84], 2.00000)))+(states[58]*states[13])/(constants[215]*constants[216]*(power(constants[278], 2.00000)))+(states[58]*states[15])/(constants[256]*constants[257]*(power(algebraic[85], 2.00000)))+(states[58]*states[18])/(constants[259]*constants[191]*(power(algebraic[94], 2.00000)))+(states[58]*states[19])/(constants[218]*constants[219]*(power(constants[279], 2.00000)))+(states[58]*states[17])/(constants[192]*constants[193]*(power(algebraic[98], 2.00000)))+(states[58]*states[11])/(constants[187]*constants[188]*(power(algebraic[83], 2.00000)))+(states[58]*states[21])/(constants[262]*constants[263]*(power(algebraic[99], 2.00000)))+(states[58]*states[23])/(constants[194]*constants[195]*(power(algebraic[104], 2.00000)))+(states[58]*states[24])/(constants[266]*constants[267]*(power(algebraic[107], 2.00000)))+(states[58]*states[26])/(constants[271]*constants[200]*(power(algebraic[124], 2.00000)))+(states[58]*states[25])/(constants[273]*constants[201]*(power(algebraic[127], 2.00000)))+(states[58]*states[54])/(constants[236]*constants[237]*(power(constants[285], 2.00000)))+(states[58]*states[58])/(constants[233]*constants[234]*(power(constants[284], 2.00000)))+(states[58]*states[7])/(constants[221]*constants[222]*(power(constants[280], 2.00000)))+(states[58]*constants[294])/(constants[224]*constants[225]*(power(constants[281], 2.00000)))+(states[58]*constants[286])/(constants[227]*constants[228]*(power(constants[282], 2.00000)))+(states[58]*states[12])/(constants[196]*constants[197]*(power(algebraic[117], 2.00000)))) algebraic[175] = -((states[54]*states[0])/(constants[231]*constants[230]*(power(constants[283], 2.00000)))+(states[54]*states[8])/(constants[302]*constants[300]*(power(algebraic[36], 2.00000)))+(states[54]*states[9])/(constants[308]*constants[306]*(power(algebraic[41], 2.00000)))+(states[54]*states[16])/(constants[314]*constants[312]*(power(algebraic[50], 2.00000)))+(states[54]*states[22])/(constants[302]*constants[300]*(power(algebraic[66], 2.00000)))+(states[54]*states[20])/(constants[308]*constants[306]*(power(algebraic[72], 2.00000)))+(states[54]*states[5])/(constants[252]*constants[251]*(power(algebraic[76], 2.00000)))+(states[54]*states[3])/(constants[207]*constants[206]*(power(constants[275], 2.00000)))+(states[54]*algebraic[143])/(constants[204]*constants[203]*(power(constants[274], 2.00000)))+(states[54]*states[4])/(constants[210]*constants[209]*(power(constants[276], 2.00000)))+(states[54]*algebraic[156])/(constants[213]*constants[212]*(power(constants[277], 2.00000)))+(states[54]*states[10])/(constants[199]*constants[269]*(power(algebraic[120], 2.00000)))+(states[54]*states[14])/(constants[190]*constants[254]*(power(algebraic[84], 2.00000)))+(states[54]*states[13])/(constants[216]*constants[215]*(power(constants[278], 2.00000)))+(states[54]*states[15])/(constants[257]*constants[256]*(power(algebraic[85], 2.00000)))+(states[54]*states[18])/(constants[191]*constants[259]*(power(algebraic[94], 2.00000)))+(states[54]*states[19])/(constants[219]*constants[218]*(power(constants[279], 2.00000)))+(states[54]*states[17])/(constants[193]*constants[192]*(power(algebraic[98], 2.00000)))+(states[54]*states[11])/(constants[188]*constants[187]*(power(algebraic[83], 2.00000)))+(states[54]*states[21])/(constants[263]*constants[262]*(power(algebraic[99], 2.00000)))+(states[54]*states[23])/(constants[195]*constants[194]*(power(algebraic[104], 2.00000)))+(states[54]*states[24])/(constants[267]*constants[266]*(power(algebraic[107], 2.00000)))+(states[54]*states[26])/(constants[200]*constants[271]*(power(algebraic[124], 2.00000)))+(states[54]*states[25])/(constants[201]*constants[273]*(power(algebraic[127], 2.00000)))+(states[54]*states[54])/(constants[237]*constants[236]*(power(constants[285], 2.00000)))+(states[54]*states[58])/(constants[234]*constants[233]*(power(constants[284], 2.00000)))+(states[54]*states[7])/(constants[222]*constants[221]*(power(constants[280], 2.00000)))+(states[54]*constants[294])/(constants[225]*constants[224]*(power(constants[281], 2.00000)))+(states[54]*constants[286])/(constants[228]*constants[227]*(power(constants[282], 2.00000)))+(states[54]*states[12])/(constants[197]*constants[196]*(power(algebraic[117], 2.00000)))) algebraic[176] = ((1.00000+states[0]/constants[229]+states[58]/constants[230])*states[0])/(constants[231]*(power(constants[283], 2.00000)))+((1.00000+states[0]/constants[296]+states[58]/constants[300])*states[8])/(constants[302]*(power(algebraic[36], 2.00000)))+((1.00000+states[0]/constants[304]+states[58]/constants[306])*states[9])/(constants[308]*(power(algebraic[41], 2.00000)))+((1.00000+states[0]/constants[310]+states[58]/constants[312])*states[16])/(constants[314]*(power(algebraic[50], 2.00000)))+((1.00000+states[0]/constants[296]+states[58]/constants[300])*states[22])/(constants[302]*(power(algebraic[66], 2.00000)))+((1.00000+states[0]/constants[304]+states[58]/constants[306])*states[20])/(constants[308]*(power(algebraic[72], 2.00000)))+((1.00000+states[0]/constants[250]+states[58]/constants[251])*states[5])/(constants[252]*(power(algebraic[76], 2.00000)))+((1.00000+states[0]/constants[205]+states[58]/constants[206])*states[3])/(constants[207]*(power(constants[275], 2.00000)))+((1.00000+states[0]/constants[202]+states[58]/constants[203])*algebraic[143])/(constants[204]*(power(constants[274], 2.00000)))+((1.00000+states[0]/constants[208]+states[58]/constants[209])*states[4])/(constants[210]*(power(constants[276], 2.00000)))+((1.00000+states[0]/constants[211]+states[58]/constants[212])*algebraic[156])/(constants[213]*(power(constants[277], 2.00000)))+((1.00000+states[0]/constants[198]+states[58]/constants[269])*states[10])/(constants[199]*(power(algebraic[120], 2.00000)))+((1.00000+states[0]/constants[189]+states[58]/constants[254])*states[14])/(constants[190]*(power(algebraic[84], 2.00000)))+((1.00000+states[0]/constants[214]+states[58]/constants[215])*states[13])/(constants[216]*(power(constants[278], 2.00000)))+((1.00000+states[0]/constants[255]+states[58]/constants[256])*states[15])/(constants[257]*(power(algebraic[85], 2.00000)))+((1.00000+states[0]/constants[258]+states[58]/constants[259])*states[18])/(constants[191]*(power(algebraic[94], 2.00000)))+((1.00000+states[0]/constants[217]+states[58]/constants[218])*states[19])/(constants[219]*(power(constants[279], 2.00000)))+((1.00000+states[0]/constants[260]+states[58]/constants[192])*states[17])/(constants[193]*(power(algebraic[98], 2.00000)))+((1.00000+states[0]/constants[253]+states[58]/constants[187])*states[11])/(constants[188]*(power(algebraic[83], 2.00000)))+((1.00000+states[0]/constants[261]+states[58]/constants[262])*states[21])/(constants[263]*(power(algebraic[99], 2.00000)))+((1.00000+states[0]/constants[264]+states[58]/constants[194])*states[23])/(constants[195]*(power(algebraic[104], 2.00000)))+((1.00000+states[0]/constants[265]+states[58]/constants[266])*states[24])/(constants[267]*(power(algebraic[107], 2.00000)))+((1.00000+states[0]/constants[270]+states[58]/constants[271])*states[26])/(constants[200]*(power(algebraic[124], 2.00000)))+((1.00000+states[0]/constants[272]+states[58]/constants[273])*states[25])/(constants[201]*(power(algebraic[127], 2.00000)))+((1.00000+states[0]/constants[235]+states[58]/constants[236])*states[54])/(constants[237]*(power(constants[285], 2.00000)))+((1.00000+states[0]/constants[232]+states[58]/constants[233])*states[58])/(constants[234]*(power(constants[284], 2.00000)))+((1.00000+states[0]/constants[220]+states[58]/constants[221])*states[7])/(constants[222]*(power(constants[280], 2.00000)))+((1.00000+states[0]/constants[223]+states[58]/constants[224])*constants[294])/(constants[225]*(power(constants[281], 2.00000)))+((1.00000+states[0]/constants[226]+states[58]/constants[227])*constants[286])/(constants[228]*(power(constants[282], 2.00000)))+((1.00000+states[0]/constants[268]+states[58]/constants[196])*states[12])/(constants[197]*(power(algebraic[117], 2.00000))) algebraic[185] = 1.00000+algebraic[176] algebraic[173] = ((1.00000+states[0]/constants[229]+states[54]/constants[231])*states[0])/(constants[230]*(power(constants[283], 2.00000)))+((1.00000+states[0]/constants[296]+states[54]/constants[302])*states[8])/(constants[300]*(power(algebraic[36], 2.00000)))+((1.00000+states[0]/constants[304]+states[54]/constants[308])*states[9])/(constants[306]*(power(algebraic[41], 2.00000)))+((1.00000+states[0]/constants[310]+states[54]/constants[314])*states[16])/(constants[312]*(power(algebraic[50], 2.00000)))+((1.00000+states[0]/constants[296]+states[54]/constants[302])*states[22])/(constants[300]*(power(algebraic[66], 2.00000)))+((1.00000+states[0]/constants[304]+states[54]/constants[308])*states[20])/(constants[306]*(power(algebraic[72], 2.00000)))+((1.00000+states[0]/constants[250]+states[54]/constants[252])*states[5])/(constants[251]*(power(algebraic[76], 2.00000)))+((1.00000+states[0]/constants[205]+states[54]/constants[207])*states[3])/(constants[206]*(power(constants[275], 2.00000)))+((1.00000+states[0]/constants[202]+states[54]/constants[204])*algebraic[143])/(constants[203]*(power(constants[274], 2.00000)))+((1.00000+states[0]/constants[208]+states[54]/constants[210])*states[4])/(constants[209]*(power(constants[276], 2.00000)))+((1.00000+states[0]/constants[211]+states[54]/constants[213])*algebraic[156])/(constants[212]*(power(constants[277], 2.00000)))+((1.00000+states[0]/constants[198]+states[54]/constants[199])*states[10])/(constants[269]*(power(algebraic[120], 2.00000)))+((1.00000+states[0]/constants[189]+states[54]/constants[190])*states[14])/(constants[254]*(power(algebraic[84], 2.00000)))+((1.00000+states[0]/constants[214]+states[54]/constants[216])*states[13])/(constants[215]*(power(constants[278], 2.00000)))+((1.00000+states[0]/constants[255]+states[54]/constants[257])*states[15])/(constants[256]*(power(algebraic[85], 2.00000)))+((1.00000+states[0]/constants[258]+states[54]/constants[191])*states[18])/(constants[259]*(power(algebraic[94], 2.00000)))+((1.00000+states[0]/constants[217]+states[54]/constants[219])*states[19])/(constants[218]*(power(constants[279], 2.00000)))+((1.00000+states[0]/constants[260]+states[54]/constants[193])*states[17])/(constants[192]*(power(algebraic[98], 2.00000)))+((1.00000+states[0]/constants[253]+states[54]/constants[188])*states[11])/(constants[187]*(power(algebraic[83], 2.00000)))+((1.00000+states[0]/constants[261]+states[54]/constants[263])*states[21])/(constants[262]*(power(algebraic[99], 2.00000)))+((1.00000+states[0]/constants[264]+states[54]/constants[195])*states[23])/(constants[194]*(power(algebraic[104], 2.00000)))+((1.00000+states[0]/constants[265]+states[54]/constants[267])*states[24])/(constants[266]*(power(algebraic[107], 2.00000)))+((1.00000+states[0]/constants[270]+states[54]/constants[200])*states[26])/(constants[271]*(power(algebraic[124], 2.00000)))+((1.00000+states[0]/constants[272]+states[54]/constants[201])*states[25])/(constants[273]*(power(algebraic[127], 2.00000)))+((1.00000+states[0]/constants[235]+states[54]/constants[237])*states[54])/(constants[236]*(power(constants[285], 2.00000)))+((1.00000+states[0]/constants[232]+states[54]/constants[234])*states[58])/(constants[233]*(power(constants[284], 2.00000)))+((1.00000+states[0]/constants[220]+states[54]/constants[222])*states[7])/(constants[221]*(power(constants[280], 2.00000)))+((1.00000+states[0]/constants[223]+states[54]/constants[225])*constants[294])/(constants[224]*(power(constants[281], 2.00000)))+((1.00000+states[0]/constants[226]+states[54]/constants[228])*constants[286])/(constants[227]*(power(constants[282], 2.00000)))+((1.00000+states[0]/constants[268]+states[54]/constants[197])*states[12])/(constants[196]*(power(algebraic[117], 2.00000))) algebraic[184] = 1.00000+algebraic[173] algebraic[170] = ((1.00000+states[58]/constants[230]+states[54]/constants[231])*states[0])/(constants[229]*(power(constants[283], 2.00000)))+((1.00000+states[58]/constants[300]+states[54]/constants[302])*states[8])/(constants[296]*(power(algebraic[36], 2.00000)))+((1.00000+states[58]/constants[306]+states[54]/constants[308])*states[9])/(constants[304]*(power(algebraic[41], 2.00000)))+((1.00000+states[58]/constants[312]+states[54]/constants[314])*states[16])/(constants[310]*(power(algebraic[50], 2.00000)))+((1.00000+states[58]/constants[300]+states[54]/constants[302])*states[22])/(constants[296]*(power(algebraic[66], 2.00000)))+((1.00000+states[58]/constants[306]+states[54]/constants[308])*states[20])/(constants[304]*(power(algebraic[72], 2.00000)))+((1.00000+states[58]/constants[251]+states[54]/constants[252])*states[5])/(constants[250]*(power(algebraic[76], 2.00000)))+((1.00000+states[58]/constants[206]+states[54]/constants[207])*states[3])/(constants[205]*(power(constants[275], 2.00000)))+((1.00000+states[58]/constants[203]+states[54]/constants[204])*algebraic[143])/(constants[202]*(power(constants[274], 2.00000)))+((1.00000+states[58]/constants[209]+states[54]/constants[210])*states[4])/(constants[208]*(power(constants[276], 2.00000)))+((1.00000+states[58]/constants[212]+states[54]/constants[213])*algebraic[156])/(constants[211]*(power(constants[277], 2.00000)))+((1.00000+states[58]/constants[269]+states[54]/constants[199])*states[10])/(constants[198]*(power(algebraic[120], 2.00000)))+((1.00000+states[58]/constants[254]+states[54]/constants[190])*states[14])/(constants[189]*(power(algebraic[84], 2.00000)))+((1.00000+states[58]/constants[215]+states[54]/constants[216])*states[13])/(constants[214]*(power(constants[278], 2.00000)))+((1.00000+states[58]/constants[256]+states[54]/constants[257])*states[15])/(constants[255]*(power(algebraic[85], 2.00000)))+((1.00000+states[58]/constants[259]+states[54]/constants[191])*states[18])/(constants[258]*(power(algebraic[94], 2.00000)))+((1.00000+states[58]/constants[218]+states[54]/constants[219])*states[19])/(constants[217]*(power(constants[279], 2.00000)))+((1.00000+states[58]/constants[192]+states[54]/constants[193])*states[17])/(constants[260]*(power(algebraic[98], 2.00000)))+((1.00000+states[58]/constants[187]+states[54]/constants[188])*states[11])/(constants[253]*(power(algebraic[83], 2.00000)))+((1.00000+states[58]/constants[262]+states[54]/constants[263])*states[21])/(constants[261]*(power(algebraic[99], 2.00000)))+((1.00000+states[58]/constants[194]+states[54]/constants[195])*states[23])/(constants[264]*(power(algebraic[104], 2.00000)))+((1.00000+states[58]/constants[266]+states[54]/constants[267])*states[24])/(constants[265]*(power(algebraic[107], 2.00000)))+((1.00000+states[58]/constants[271]+states[54]/constants[200])*states[26])/(constants[270]*(power(algebraic[124], 2.00000)))+((1.00000+states[58]/constants[273]+states[54]/constants[201])*states[25])/(constants[272]*(power(algebraic[127], 2.00000)))+((1.00000+states[58]/constants[236]+states[54]/constants[237])*states[54])/(constants[235]*(power(constants[285], 2.00000)))+((1.00000+states[58]/constants[233]+states[54]/constants[234])*states[58])/(constants[232]*(power(constants[284], 2.00000)))+((1.00000+states[58]/constants[221]+states[54]/constants[222])*states[7])/(constants[220]*(power(constants[280], 2.00000)))+((1.00000+states[58]/constants[224]+states[54]/constants[225])*constants[294])/(constants[223]*(power(constants[281], 2.00000)))+((1.00000+states[58]/constants[227]+states[54]/constants[228])*constants[286])/(constants[226]*(power(constants[282], 2.00000)))+((1.00000+states[58]/constants[196]+states[54]/constants[197])*states[12])/(constants[268]*(power(algebraic[117], 2.00000))) algebraic[183] = 1.00000+algebraic[170]+constants[243]/(constants[244]*(power(1.00000+states[0]/constants[244], 2.00000))) algebraic[174] = -((states[54]*states[0])/(constants[231]*constants[229]*(power(constants[283], 2.00000)))+(states[54]*states[8])/(constants[302]*constants[296]*(power(algebraic[36], 2.00000)))+(states[54]*states[9])/(constants[308]*constants[304]*(power(algebraic[41], 2.00000)))+(states[54]*states[16])/(constants[314]*constants[310]*(power(algebraic[50], 2.00000)))+(states[54]*states[22])/(constants[302]*constants[296]*(power(algebraic[66], 2.00000)))+(states[54]*states[20])/(constants[308]*constants[304]*(power(algebraic[72], 2.00000)))+(states[54]*states[5])/(constants[252]*constants[250]*(power(algebraic[76], 2.00000)))+(states[54]*states[3])/(constants[207]*constants[205]*(power(constants[275], 2.00000)))+(states[54]*algebraic[143])/(constants[204]*constants[202]*(power(constants[274], 2.00000)))+(states[54]*states[4])/(constants[210]*constants[208]*(power(constants[276], 2.00000)))+(states[54]*algebraic[156])/(constants[213]*constants[211]*(power(constants[277], 2.00000)))+(states[54]*states[10])/(constants[199]*constants[198]*(power(algebraic[120], 2.00000)))+(states[54]*states[14])/(constants[190]*constants[189]*(power(algebraic[84], 2.00000)))+(states[54]*states[13])/(constants[216]*constants[214]*(power(constants[278], 2.00000)))+(states[54]*states[15])/(constants[257]*constants[255]*(power(algebraic[85], 2.00000)))+(states[54]*states[18])/(constants[191]*constants[258]*(power(algebraic[94], 2.00000)))+(states[54]*states[19])/(constants[219]*constants[217]*(power(constants[279], 2.00000)))+(states[54]*states[17])/(constants[193]*constants[260]*(power(algebraic[98], 2.00000)))+(states[54]*states[11])/(constants[188]*constants[253]*(power(algebraic[83], 2.00000)))+(states[54]*states[21])/(constants[263]*constants[261]*(power(algebraic[99], 2.00000)))+(states[54]*states[23])/(constants[195]*constants[264]*(power(algebraic[104], 2.00000)))+(states[54]*states[24])/(constants[267]*constants[265]*(power(algebraic[107], 2.00000)))+(states[54]*states[26])/(constants[200]*constants[270]*(power(algebraic[124], 2.00000)))+(states[54]*states[25])/(constants[201]*constants[272]*(power(algebraic[127], 2.00000)))+(states[54]*states[54])/(constants[237]*constants[235]*(power(constants[285], 2.00000)))+(states[54]*states[58])/(constants[234]*constants[232]*(power(constants[284], 2.00000)))+(states[54]*states[7])/(constants[222]*constants[220]*(power(constants[280], 2.00000)))+(states[54]*constants[294])/(constants[225]*constants[223]*(power(constants[281], 2.00000)))+(states[54]*constants[286])/(constants[228]*constants[226]*(power(constants[282], 2.00000)))+(states[54]*states[12])/(constants[197]*constants[268]*(power(algebraic[117], 2.00000)))) algebraic[171] = -((states[58]*states[0])/(constants[230]*constants[229]*(power(constants[283], 2.00000)))+(states[58]*states[8])/(constants[300]*constants[296]*(power(algebraic[36], 2.00000)))+(states[58]*states[9])/(constants[306]*constants[304]*(power(algebraic[41], 2.00000)))+(states[58]*states[16])/(constants[312]*constants[310]*(power(algebraic[50], 2.00000)))+(states[58]*states[22])/(constants[300]*constants[296]*(power(algebraic[66], 2.00000)))+(states[58]*states[20])/(constants[306]*constants[304]*(power(algebraic[72], 2.00000)))+(states[58]*states[5])/(constants[251]*constants[250]*(power(algebraic[76], 2.00000)))+(states[58]*states[3])/(constants[206]*constants[205]*(power(constants[275], 2.00000)))+(states[58]*algebraic[143])/(constants[203]*constants[202]*(power(constants[274], 2.00000)))+(states[58]*states[4])/(constants[209]*constants[208]*(power(constants[276], 2.00000)))+(states[58]*algebraic[156])/(constants[212]*constants[211]*(power(constants[277], 2.00000)))+(states[58]*states[10])/(constants[269]*constants[198]*(power(algebraic[120], 2.00000)))+(states[58]*states[14])/(constants[254]*constants[189]*(power(algebraic[84], 2.00000)))+(states[58]*states[13])/(constants[215]*constants[214]*(power(constants[278], 2.00000)))+(states[58]*states[15])/(constants[256]*constants[255]*(power(algebraic[85], 2.00000)))+(states[58]*states[18])/(constants[259]*constants[258]*(power(algebraic[94], 2.00000)))+(states[58]*states[19])/(constants[218]*constants[217]*(power(constants[279], 2.00000)))+(states[58]*states[17])/(constants[192]*constants[260]*(power(algebraic[98], 2.00000)))+(states[58]*states[11])/(constants[187]*constants[253]*(power(algebraic[83], 2.00000)))+(states[58]*states[21])/(constants[262]*constants[261]*(power(algebraic[99], 2.00000)))+(states[58]*states[23])/(constants[194]*constants[264]*(power(algebraic[104], 2.00000)))+(states[58]*states[24])/(constants[266]*constants[265]*(power(algebraic[107], 2.00000)))+(states[58]*states[26])/(constants[271]*constants[270]*(power(algebraic[124], 2.00000)))+(states[58]*states[25])/(constants[273]*constants[272]*(power(algebraic[127], 2.00000)))+(states[58]*states[54])/(constants[236]*constants[235]*(power(constants[285], 2.00000)))+(states[58]*states[58])/(constants[233]*constants[232]*(power(constants[284], 2.00000)))+(states[58]*states[7])/(constants[221]*constants[220]*(power(constants[280], 2.00000)))+(states[58]*constants[294])/(constants[224]*constants[223]*(power(constants[281], 2.00000)))+(states[58]*constants[286])/(constants[227]*constants[226]*(power(constants[282], 2.00000)))+(states[58]*states[12])/(constants[196]*constants[268]*(power(algebraic[117], 2.00000)))) algebraic[186] = (((algebraic[183]*algebraic[175]*algebraic[172]+algebraic[185]*algebraic[168]*algebraic[171]+algebraic[184]*algebraic[169]*algebraic[174])-algebraic[184]*algebraic[185]*algebraic[183])-algebraic[169]*algebraic[175]*algebraic[171])-algebraic[168]*algebraic[172]*algebraic[174] algebraic[135] = ((algebraic[58]+algebraic[78])-algebraic[63])/constants[34] algebraic[136] = ((-algebraic[58]-algebraic[78])+algebraic[63])/constants[34] algebraic[137] = (algebraic[102]-algebraic[58])/constants[34] algebraic[138] = (-algebraic[102]+algebraic[58])/constants[34] algebraic[139] = (((-algebraic[102]-algebraic[78])+algebraic[82])-algebraic[115])/constants[34] algebraic[159] = ((algebraic[144]+algebraic[145]+algebraic[146]+algebraic[147])-algebraic[157])/constants[34] algebraic[166] = ((algebraic[158]+algebraic[157])-algebraic[164])/constants[34] algebraic[148] = (-algebraic[144]+algebraic[123])/constants[34] algebraic[149] = (algebraic[144]-algebraic[89])/constants[34] algebraic[140] = (-algebraic[97]+algebraic[114]+algebraic[89])/constants[34] algebraic[150] = (algebraic[97]-algebraic[145])/constants[34] algebraic[151] = (((algebraic[145]-algebraic[146])-algebraic[131])+algebraic[13])/constants[34] algebraic[152] = (algebraic[146]-algebraic[102])/constants[34] algebraic[153] = ((-algebraic[144]-algebraic[146])+algebraic[102]+algebraic[89])/constants[34] algebraic[160] = ((algebraic[102]-algebraic[158])+algebraic[116])/constants[34] algebraic[161] = (algebraic[158]-algebraic[110])/constants[34] algebraic[154] = (((((algebraic[110]-algebraic[147])+algebraic[115])-algebraic[13])-algebraic[114])-algebraic[116])/constants[34] algebraic[155] = (-algebraic[89]+algebraic[147]+algebraic[131])/constants[34] algebraic[141] = ((algebraic[131]+algebraic[126])-algebraic[134])/constants[34] algebraic[142] = (-algebraic[131]+algebraic[134])/constants[34] algebraic[162] = -algebraic[159] algebraic[167] = -algebraic[166] algebraic[177] = (((-algebraic[144]+2.00000*algebraic[89])-algebraic[146])+algebraic[102]+algebraic[147])/constants[34] algebraic[34] = constants[30]*(states[53]*states[0]-states[54]*states[1]) algebraic[178] = (((((((algebraic[123]+algebraic[126]+algebraic[114])-algebraic[134])-5.00000*algebraic[157])-2.00000*algebraic[164])-4.00000*algebraic[165])+(constants[38]-1.00000)*algebraic[78]+2.00000*algebraic[82]+algebraic[32])-algebraic[34])/constants[34] algebraic[180] = (states[0]*constants[287])/(constants[229]*constants[283])+(states[0]*algebraic[135])/(constants[296]*algebraic[36])+(states[0]*algebraic[136])/(constants[304]*algebraic[41])+(states[0]*constants[238])/(constants[310]*algebraic[50])+(states[0]*algebraic[137])/(constants[296]*algebraic[66])+(states[0]*algebraic[138])/(constants[304]*algebraic[72])+(states[0]*algebraic[139])/(constants[250]*algebraic[76])+(states[0]*algebraic[159])/(constants[205]*constants[275])+(states[0]*algebraic[162])/(constants[202]*constants[274])+(states[0]*algebraic[166])/(constants[208]*constants[276])+(states[0]*algebraic[167])/(constants[211]*constants[277])+(states[0]*algebraic[148])/(constants[198]*algebraic[120])+(states[0]*algebraic[155])/(constants[189]*algebraic[84])+(states[0]*algebraic[149])/(constants[214]*constants[278])+(states[0]*algebraic[140])/(constants[255]*algebraic[85])+(states[0]*algebraic[150])/(constants[258]*algebraic[94])+(states[0]*algebraic[151])/(constants[217]*constants[279])+(states[0]*algebraic[152])/(constants[260]*algebraic[98])+(states[0]*algebraic[153])/(constants[253]*algebraic[83])+(states[0]*algebraic[160])/(constants[261]*algebraic[99])+(states[0]*algebraic[161])/(constants[264]*algebraic[104])+(states[0]*algebraic[154])/(constants[265]*algebraic[107])+(states[0]*algebraic[141])/(constants[270]*algebraic[124])+(states[0]*algebraic[142])/(constants[272]*algebraic[127])+(states[0]*constants[289])/(constants[235]*constants[285])+(states[0]*constants[288])/(constants[232]*constants[284])+(states[0]*constants[239])/(constants[220]*constants[280])+(states[0]*constants[242])/(constants[223]*constants[281])+(states[0]*constants[241])/(constants[226]*constants[282])+(states[0]*constants[240])/(constants[268]*algebraic[117])+algebraic[177]+algebraic[178] algebraic[181] = (states[58]*constants[287])/(constants[230]*constants[283])+(states[58]*algebraic[135])/(constants[300]*algebraic[36])+(states[58]*algebraic[136])/(constants[306]*algebraic[41])+(states[58]*constants[238])/(constants[312]*algebraic[50])+(states[58]*algebraic[137])/(constants[300]*algebraic[66])+(states[58]*algebraic[138])/(constants[306]*algebraic[72])+(states[58]*algebraic[139])/(constants[251]*algebraic[76])+(states[58]*algebraic[159])/(constants[206]*constants[275])+(states[58]*algebraic[162])/(constants[203]*constants[274])+(states[58]*algebraic[166])/(constants[209]*constants[276])+(states[58]*algebraic[167])/(constants[212]*constants[277])+(states[58]*algebraic[148])/(constants[269]*algebraic[120])+(states[58]*algebraic[155])/(constants[254]*algebraic[84])+(states[58]*algebraic[149])/(constants[215]*constants[278])+(states[58]*algebraic[140])/(constants[256]*algebraic[85])+(states[58]*algebraic[150])/(constants[259]*algebraic[94])+(states[58]*algebraic[151])/(constants[218]*constants[279])+(states[58]*algebraic[152])/(constants[192]*algebraic[98])+(states[58]*algebraic[153])/(constants[187]*algebraic[83])+(states[58]*algebraic[160])/(constants[262]*algebraic[99])+(states[58]*algebraic[161])/(constants[194]*algebraic[104])+(states[58]*algebraic[154])/(constants[266]*algebraic[107])+(states[58]*algebraic[141])/(constants[271]*algebraic[124])+(states[58]*algebraic[142])/(constants[273]*algebraic[127])+(states[58]*constants[289])/(constants[236]*constants[285])+(states[58]*constants[288])/(constants[233]*constants[284])+(states[58]*constants[239])/(constants[221]*constants[280])+(states[58]*constants[242])/(constants[224]*constants[281])+(states[58]*constants[241])/(constants[227]*constants[282])+(states[58]*constants[240])/(constants[196]*algebraic[117])+constants[290] algebraic[179] = algebraic[34]/constants[34] algebraic[182] = (states[54]*constants[287])/(constants[231]*constants[283])+(states[54]*algebraic[135])/(constants[302]*algebraic[36])+(states[54]*algebraic[136])/(constants[308]*algebraic[41])+(states[54]*constants[238])/(constants[314]*algebraic[50])+(states[54]*algebraic[137])/(constants[302]*algebraic[66])+(states[54]*algebraic[138])/(constants[308]*algebraic[72])+(states[54]*algebraic[139])/(constants[252]*algebraic[76])+(states[54]*algebraic[159])/(constants[207]*constants[275])+(states[54]*algebraic[162])/(constants[204]*constants[274])+(states[54]*algebraic[166])/(constants[210]*constants[276])+(states[54]*algebraic[167])/(constants[213]*constants[277])+(states[54]*algebraic[148])/(constants[199]*algebraic[120])+(states[54]*algebraic[155])/(constants[190]*algebraic[84])+(states[54]*algebraic[149])/(constants[216]*constants[278])+(states[54]*algebraic[140])/(constants[257]*algebraic[85])+(states[54]*algebraic[150])/(constants[191]*algebraic[94])+(states[54]*algebraic[151])/(constants[219]*constants[279])+(states[54]*algebraic[152])/(constants[193]*algebraic[98])+(states[54]*algebraic[153])/(constants[188]*algebraic[83])+(states[54]*algebraic[160])/(constants[263]*algebraic[99])+(states[54]*algebraic[161])/(constants[195]*algebraic[104])+(states[54]*algebraic[154])/(constants[267]*algebraic[107])+(states[54]*algebraic[141])/(constants[200]*algebraic[124])+(states[54]*algebraic[142])/(constants[201]*algebraic[127])+(states[54]*constants[289])/(constants[237]*constants[285])+(states[54]*constants[288])/(constants[234]*constants[284])+(states[54]*constants[239])/(constants[222]*constants[280])+(states[54]*constants[242])/(constants[225]*constants[281])+(states[54]*constants[241])/(constants[228]*constants[282])+(states[54]*constants[240])/(constants[197]*algebraic[117])+algebraic[179] rates[0] = ((algebraic[175]*algebraic[172]-algebraic[184]*algebraic[185])*algebraic[180]+(algebraic[185]*algebraic[168]-algebraic[169]*algebraic[175])*algebraic[181]+(algebraic[184]*algebraic[169]-algebraic[168]*algebraic[172])*algebraic[182])/algebraic[186] rates[58] = ((algebraic[185]*algebraic[171]-algebraic[174]*algebraic[172])*algebraic[180]+(algebraic[174]*algebraic[169]-algebraic[183]*algebraic[185])*algebraic[181]+(algebraic[183]*algebraic[172]-algebraic[169]*algebraic[171])*algebraic[182])/algebraic[186] rates[54] = ((algebraic[184]*algebraic[174]-algebraic[175]*algebraic[171])*algebraic[180]+(algebraic[183]*algebraic[175]-algebraic[174]*algebraic[168])*algebraic[181]+(algebraic[171]*algebraic[168]-algebraic[183]*algebraic[184])*algebraic[182])/algebraic[186] return(rates) def computeAlgebraic(constants, states, voi): algebraic = array([[0.0] * len(voi)] * sizeAlgebraic) states = array(states) voi = array(voi) algebraic[1] = constants[36]*constants[41]*(states[44]-states[45]) algebraic[2] = constants[36]*constants[170]*(states[51]-states[52]) algebraic[13] = (constants[17]*(states[30]*states[24]-states[19]*states[29]))/(constants[161]*constants[164]*(2.00000+states[29]/constants[163]+states[24]/constants[164]+states[30]/constants[161]+states[19]/constants[162]+(states[29]*states[19])/(constants[163]*constants[162])+(states[24]*states[30])/(constants[164]*constants[161]))) algebraic[15] = constants[36]*constants[41]*(states[47]-states[30]) algebraic[18] = constants[36]*constants[39]*(states[40]-states[31]) algebraic[20] = constants[36]*constants[41]*(states[41]-states[27]) algebraic[22] = constants[36]*constants[41]*(states[42]-states[43]) algebraic[24] = constants[36]*constants[41]*(states[46]-states[29]) algebraic[26] = constants[36]*constants[41]*(states[48]-states[33]) algebraic[28] = constants[36]*constants[41]*(states[49]-states[28]) algebraic[30] = constants[36]*constants[41]*(states[50]-states[32]) algebraic[19] = states[1] algebraic[23] = states[59] algebraic[25] = states[53] algebraic[27] = 1.00000+algebraic[19]/constants[296]+algebraic[23]/constants[300]+algebraic[25]/constants[302] algebraic[29] = ((states[34]*algebraic[23])/constants[300])/algebraic[27] algebraic[31] = 1.00000+algebraic[19]/constants[304]+algebraic[23]/constants[306]+algebraic[25]/constants[308] algebraic[37] = ((states[36]*algebraic[23])/constants[306])/algebraic[31] algebraic[21] = 1.00000+algebraic[19]/constants[316] algebraic[33] = (constants[292]*algebraic[21]*algebraic[31])/(algebraic[19]*algebraic[27]) algebraic[35] = (algebraic[33]*constants[300]*algebraic[27])/(constants[306]*algebraic[31]) algebraic[39] = ((constants[245]/(constants[179]*constants[182]))*(algebraic[29]*states[56]-(states[57]*algebraic[37])/algebraic[35]))/(1.00000+algebraic[29]/constants[178]+states[56]/constants[179]+(algebraic[29]*states[56])/(constants[179]*constants[182])+states[57]/constants[180]+algebraic[37]/constants[181]+(states[57]*algebraic[37])/(constants[181]*constants[184])+(states[57]*states[56])/(constants[186]*constants[179])) algebraic[17] = constants[36]*constants[40]*(states[38]-states[39]) algebraic[42] = 1.00000+algebraic[19]/constants[310]+algebraic[23]/constants[312]+algebraic[25]/constants[314] algebraic[44] = (constants[176]*algebraic[27]*algebraic[42])/(algebraic[31]*algebraic[31]) algebraic[47] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[175]*(algebraic[44]*states[36]*states[36]-states[38]*states[34])]) algebraic[14] = constants[36]*constants[40]*(states[34]-states[35]) algebraic[49] = constants[49]-states[55] algebraic[51] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[177]*(((constants[293]*states[36])/algebraic[31])*states[55]*algebraic[19]-(states[34]/algebraic[27])*algebraic[49])]) algebraic[16] = constants[36]*constants[40]*(states[36]-states[37]) algebraic[36] = 1.00000+states[0]/constants[296]+states[58]/constants[300]+states[54]/constants[302] algebraic[41] = 1.00000+states[0]/constants[304]+states[58]/constants[306]+states[54]/constants[308] algebraic[60] = 1.00000+states[1]/constants[304]+states[59]/constants[306]+states[53]/constants[308] algebraic[59] = 1.00000+states[1]/constants[296]+states[59]/constants[300]+states[53]/constants[302] algebraic[63] = custom_piecewise([greater(states[37] , constants[171]) & greater(states[8] , constants[171]), (constants[27]/(1.00000+(constants[42]/states[37])/algebraic[60]))*((states[37]/algebraic[60])/(states[37]/algebraic[60]+(states[35]/algebraic[59])*exp(((constants[43]-1.00000)*constants[33]*states[2])/constants[249]))+(states[9]/algebraic[41])/(states[9]/algebraic[41]+(states[8]/algebraic[36])*exp((constants[43]*constants[33]*states[2])/constants[249]))) , True, 0.00000]) algebraic[67] = 1.00000+states[1]/constants[310]+states[59]/constants[312]+states[53]/constants[314] algebraic[69] = (constants[174]*algebraic[59]*algebraic[67])/(algebraic[60]*algebraic[60]) algebraic[73] = custom_piecewise([equal(constants[0] , 1.00000) | equal(constants[0] , 2.00000) | equal(constants[0] , 5.00000), 0.00000 , True, constants[173]*(algebraic[69]*states[37]*states[37]-states[39]*states[35])]) algebraic[76] = 1.00000+states[0]/constants[250]+states[58]/constants[251]+states[54]/constants[252] algebraic[8] = exp(-(constants[56]-constants[38]*constants[33]*states[2])/constants[249]) algebraic[77] = (algebraic[8]*(power(states[1], constants[38]))*algebraic[36])/((power(states[0], constants[38]-1.00000))*algebraic[41]*algebraic[76]*1.00000) algebraic[78] = constants[26]*(algebraic[77]*states[9]*states[5]-states[8]) algebraic[50] = 1.00000+states[0]/constants[310]+states[58]/constants[312]+states[54]/constants[314] algebraic[52] = (states[16]*(1.00000+states[0]/constants[310]))/algebraic[50] algebraic[56] = 1.00000+algebraic[52]/constants[148] algebraic[57] = custom_piecewise([greater(states[22] , constants[149]) & greater(states[9] , constants[149]), (((constants[12]*constants[247])/algebraic[56])*states[22]*states[9])/(constants[247]*constants[145]*states[22]+constants[247]*constants[144]*states[9]+(constants[12]*constants[147]*states[20])/constants[139]+(constants[12]*constants[146]*states[8])/constants[139]+constants[247]*states[22]*states[9]+(constants[12]*constants[147]*states[22]*states[20])/(constants[139]*constants[140])+(constants[12]*states[20]*states[8])/constants[139]+(constants[247]*constants[144]*states[9]*states[8])/constants[143]) , True, 0.00000]) algebraic[58] = custom_piecewise([greater(states[20] , constants[149]) & greater(states[8] , constants[149]), algebraic[57]-((((constants[12]*constants[247])/algebraic[56])*states[20]*states[8])/constants[139])/(constants[247]*constants[145]*states[22]+constants[247]*constants[144]*states[9]+(constants[12]*constants[147]*states[20])/constants[139]+(constants[12]*constants[146]*states[8])/constants[139]+constants[247]*states[22]*states[9]+(constants[12]*constants[147]*states[22]*states[20])/(constants[139]*constants[140])+(constants[12]*states[20]*states[8])/constants[139]+(constants[247]*constants[144]*states[9]*states[8])/constants[143]) , True, algebraic[57]]) algebraic[83] = 1.00000+states[0]/constants[253] algebraic[98] = 1.00000+states[0]/constants[260] algebraic[99] = 1.00000+states[0]/constants[261]+states[58]/constants[262]+states[54]/constants[263] algebraic[66] = algebraic[36] algebraic[72] = algebraic[41] algebraic[100] = (constants[96]*algebraic[83]*algebraic[99]*algebraic[66])/(states[0]*algebraic[72]*algebraic[98]*algebraic[76]) algebraic[101] = (constants[8]*constants[106]*constants[101]*constants[102])/(algebraic[100]*constants[97]*constants[98]*constants[105]) algebraic[102] = (constants[8]*algebraic[101]*(states[20]*states[17]*states[5]-(states[11]*states[21]*states[22])/algebraic[100]))/(algebraic[101]*constants[97]*constants[98]*constants[106]+algebraic[101]*constants[98]*constants[106]*states[20]+algebraic[101]*constants[97]*constants[104]*states[5]+algebraic[101]*constants[106]*states[20]*states[17]+algebraic[101]*constants[104]*states[20]*states[5]+algebraic[101]*constants[103]*states[17]*states[5]+algebraic[101]*states[20]*states[17]*states[5]+(constants[8]*constants[102]*constants[107]*states[11])/algebraic[100]+(constants[8]*constants[101]*constants[106]*states[22])/algebraic[100]+(constants[8]*constants[108]*states[11]*states[21])/algebraic[100]+(constants[8]*constants[107]*states[11]*states[22])/algebraic[100]+(constants[8]*constants[106]*states[21]*states[22])/algebraic[100]+(constants[8]*states[11]*states[21]*states[22])/algebraic[100]+(constants[8]*constants[107]*constants[102]*states[20]*states[11])/(constants[97]*algebraic[100])+(algebraic[101]*constants[97]*constants[104]*states[5]*states[22])/constants[102]+(constants[8]*constants[107]*constants[102]*states[20]*states[17]*states[21])/(constants[97]*constants[98]*algebraic[100])+(algebraic[101]*constants[103]*states[17]*states[5]*states[22])/constants[102]+(constants[8]*constants[108]*states[20]*states[11]*states[21])/(constants[97]*algebraic[100])+(algebraic[101]*constants[97]*constants[104]*states[5]*states[21]*states[22])/(constants[101]*constants[102])+(constants[8]*constants[102]*constants[107]*states[20]*states[17]*states[5]*states[11])/(constants[97]*constants[98]*constants[99]*algebraic[100])+(constants[8]*constants[100]*constants[108]*states[20]*states[17]*states[5]*states[21])/(constants[97]*constants[98]*constants[100]*algebraic[100])+(constants[8]*constants[108]*states[20]*states[17]*states[11]*states[21])/(constants[97]*constants[98]*algebraic[100])+(constants[8]*constants[103]*states[17]*states[5]*states[21]*states[22])/(constants[101]*constants[102])+(algebraic[101]*constants[103]*constants[100]*states[17]*states[11]*states[21]*states[22])/(constants[100]*constants[101]*constants[102])+(algebraic[101]*constants[97]*constants[104]*states[5]*states[11]*states[21]*states[22])/(constants[100]*constants[101]*constants[102])+(constants[8]*constants[108]*states[20]*states[17]*states[5]*states[11]*states[21])/(constants[97]*constants[98]*constants[100]*algebraic[100])+(constants[8]*constants[103]*states[17]*states[5]*states[11]*states[21]*states[22])/(constants[100]*constants[101]*constants[102])) algebraic[85] = 1.00000+states[0]/constants[255]+states[58]/constants[256]+states[54]/constants[257] algebraic[84] = 1.00000+states[58]/constants[254] algebraic[86] = (constants[64]*algebraic[85])/(states[0]*states[0]*algebraic[84]) algebraic[87] = (states[15]*(1.00000+states[0]/constants[255]))/algebraic[85] algebraic[88] = 1.00000+algebraic[87]/constants[68] algebraic[38] = (states[8]*(1.00000+states[0]/constants[296]))/algebraic[36] algebraic[43] = (states[9]*(1.00000+states[0]/constants[304]))/algebraic[41] algebraic[54] = 1.00000+algebraic[38]/constants[69]+algebraic[43]/constants[70]+algebraic[52]/constants[71]+states[11]/constants[72]+states[17]/constants[73] algebraic[89] = (constants[3]*(states[13]*states[14]-(states[11]*states[15])/algebraic[86]))/(constants[67]*constants[66]*algebraic[88]+constants[65]*algebraic[88]*states[13]+constants[66]*algebraic[54]*states[14]+states[14]*states[13]) algebraic[94] = 1.00000+states[0]/constants[258]+states[58]/constants[259] algebraic[95] = (constants[74]*algebraic[94])/algebraic[85] algebraic[96] = (constants[4]*constants[76])/(constants[75]*algebraic[95]) algebraic[97] = (constants[4]*algebraic[96]*(states[15]-states[18]/algebraic[95]))/(constants[75]*algebraic[96]+algebraic[96]*states[15]+(constants[4]*states[18])/algebraic[95]) algebraic[107] = 1.00000+states[0]/constants[265]+states[58]/constants[266]+states[54]/constants[267] algebraic[91] = 1.00000+states[1]/constants[255]+states[59]/constants[256]+states[53]/constants[257] algebraic[92] = ((states[43]*states[1])/constants[255])/algebraic[91] algebraic[113] = 1.00000+states[1]/constants[265]+states[59]/constants[266]+states[53]/constants[267] algebraic[90] = ((states[15]*states[0])/constants[255])/algebraic[85] algebraic[114] = constants[16]*((algebraic[92]*states[24])/algebraic[107]-(algebraic[90]*states[29])/algebraic[113]) algebraic[80] = 1.00000+states[1]/constants[250]+states[59]/constants[251]+states[53]/constants[252] algebraic[115] = constants[19]*(((states[29]/algebraic[113])*states[5])/algebraic[76]-((states[24]/algebraic[107])*states[31])/algebraic[80]) algebraic[81] = ((states[31]*states[1])/constants[250])/algebraic[80] algebraic[79] = ((states[5]*states[0])/constants[250])/algebraic[76] algebraic[82] = (constants[28]*(algebraic[81]*states[1]-algebraic[79]*states[0]))/(constants[29]*(1.00000+algebraic[81]/constants[29])*(1.00000+algebraic[79]/constants[29])) algebraic[103] = 1.00000+states[1]/constants[261]+states[59]/constants[262]+states[53]/constants[263] algebraic[116] = constants[18]*(((states[33]/algebraic[103])*states[24])/algebraic[107]-((states[21]/algebraic[99])*states[29])/algebraic[113]) algebraic[120] = 1.00000+states[58]/constants[269] algebraic[122] = 1.00000+states[59]/constants[269] algebraic[123] = constants[14]*((states[27]/algebraic[122])*states[1]-(states[10]/algebraic[120])*states[0]) algebraic[127] = 1.00000+states[0]/constants[272]+states[58]/constants[273] algebraic[124] = 1.00000+states[0]/constants[270]+states[58]/constants[271] algebraic[128] = (constants[150]*algebraic[124]*algebraic[84])/algebraic[127] algebraic[129] = (constants[13]*constants[158]*constants[153])/(algebraic[128]*constants[151]*constants[156]) algebraic[12] = 1.00000+states[19]/constants[159] algebraic[130] = custom_piecewise([greater(states[25] , constants[160]) & greater(states[19] , constants[160]), (constants[13]*algebraic[129]*states[25]*states[19])/(algebraic[129]*constants[156]*states[25]+algebraic[129]*constants[155]*algebraic[12]*states[19]+(constants[13]*constants[158]*algebraic[12]*states[14])/algebraic[128]+(constants[13]*constants[157]*states[26])/algebraic[128]+algebraic[129]*states[25]*states[19]+(constants[13]*constants[158]*states[25]*states[14])/(algebraic[128]*constants[151])+(constants[13]*states[14]*states[26])/algebraic[128]+(algebraic[129]*constants[155]*states[19]*states[26])/constants[154]) , True, 0.00000]) algebraic[131] = custom_piecewise([greater(states[14] , constants[160]) & greater(states[26] , constants[160]), algebraic[130]-((constants[13]*algebraic[129]*states[14]*states[26])/algebraic[128])/(algebraic[129]*constants[156]*states[25]+algebraic[129]*constants[155]*algebraic[12]*states[19]+(constants[13]*constants[158]*algebraic[12]*states[14])/algebraic[128]+(constants[13]*constants[157]*states[26])/algebraic[128]+algebraic[129]*states[25]*states[19]+(constants[13]*constants[158]*states[25]*states[14])/(algebraic[128]*constants[151])+(constants[13]*states[14]*states[26])/algebraic[128]+(algebraic[129]*constants[155]*states[19]*states[26])/constants[154]) , True, algebraic[130]]) algebraic[125] = 1.00000+states[1]/constants[270]+states[59]/constants[271] algebraic[126] = constants[15]*((states[28]/algebraic[125])*states[1]-(states[26]/algebraic[124])*states[0]) algebraic[132] = 1.00000+states[1]/constants[272]+states[59]/constants[273] algebraic[133] = (exp((-constants[33]*states[2])/constants[249])*algebraic[127]*algebraic[125])/(algebraic[132]*algebraic[124]) algebraic[134] = (constants[20]*(algebraic[133]*states[32]*states[26]*states[0]-states[25]*states[28]*states[1]))/(algebraic[133]*constants[165]*constants[168]*constants[248]*(2.00000*constants[169]+(constants[169]*states[32])/constants[165]+(states[32]*states[26]*states[0])/(constants[165]*constants[168]*constants[248])+(constants[169]*states[25]*states[1])/(constants[166]*constants[248])+(states[25]*states[28]*states[1])/(constants[166]*constants[167]*constants[248])+(constants[169]*states[25])/constants[166]+(constants[169]*states[32]*states[0])/(constants[165]*constants[248])+(constants[169]*states[0])/constants[248]+(constants[169]*states[28]*states[1])/(constants[167]*constants[248])+(constants[169]*states[1])/constants[248]+(constants[169]*states[26]*states[0])/(constants[168]*constants[248]))) algebraic[143] = custom_piecewise([greater(states[3] , constants[45]), 0.00000 , True, constants[45]-states[3]]) algebraic[117] = 1.00000+states[0]/constants[268] algebraic[121] = (constants[57]*algebraic[117])/(states[0]*algebraic[120]*algebraic[83]) algebraic[9] = 1.00000+states[13]/constants[61] algebraic[10] = 1.00000+states[3]/constants[62] algebraic[144] = custom_piecewise([greater(states[10] , constants[63]) & greater(states[11] , constants[63]) & greater(algebraic[143] , constants[63]), (constants[2]*(states[10]*states[11]*algebraic[143]-(states[12]*states[13]*states[3])/algebraic[121]))/(constants[60]*algebraic[10]*states[10]*states[11]+constants[59]*algebraic[9]*states[10]*algebraic[143]+constants[58]*states[11]*algebraic[143]+states[10]*states[11]*algebraic[143]) , True, 0.00000]) algebraic[118] = (constants[77]*algebraic[117])/(states[0]*states[0]*algebraic[94]) algebraic[46] = 1.00000+(constants[84]/algebraic[43])*(1.00000+algebraic[38]/constants[83]) algebraic[145] = custom_piecewise([greater(algebraic[143] , constants[85]) & greater(states[18] , constants[85]), (constants[5]*(1.00000-(states[12]*states[19]*states[3])/(algebraic[118]*algebraic[143]*states[18])))/(1.00000+(power(constants[79]/states[18], constants[82]))*algebraic[46]+(constants[78]/algebraic[143])*(1.00000+(power(constants[80]/states[18], constants[82]))*algebraic[46]+(states[3]*algebraic[46])/constants[81])) , True, 0.00000]) algebraic[119] = (constants[86]*algebraic[117]*algebraic[98])/(states[0]*algebraic[83]) algebraic[48] = 1.00000+(constants[94]/algebraic[43])*(1.00000+algebraic[38]/constants[93]) algebraic[146] = custom_piecewise([greater(states[19] , constants[95]) & greater(states[11] , constants[95]) & greater(algebraic[143] , constants[95]), ((constants[6]*(1.00000-(states[12]*states[17]*states[3])/(algebraic[119]*states[19]*states[11]*algebraic[143])))/(1.00000+(constants[87]/states[19])*algebraic[48]+(constants[88]/states[11])*(1.00000+states[17]/constants[90])+(constants[89]/algebraic[143])*(1.00000+states[3]/constants[91])))/(1.00000+states[3]/constants[92]) , True, 0.00000]) algebraic[104] = 1.00000+states[0]/constants[264] algebraic[108] = (constants[119]*algebraic[107])/algebraic[104] algebraic[109] = (constants[10]*constants[121])/(algebraic[108]*constants[120]) algebraic[68] = algebraic[38] algebraic[74] = algebraic[43] algebraic[75] = 1.00000+states[15]/constants[122]+algebraic[38]/constants[123]+algebraic[43]/constants[124]+algebraic[68]/constants[125]+algebraic[74]/constants[126] algebraic[110] = (algebraic[109]*constants[10]*(states[23]-states[24]/algebraic[108]))/(constants[120]*algebraic[109]*algebraic[75]+algebraic[109]*states[23]+(constants[10]*states[24])/algebraic[108]) algebraic[111] = (constants[127]*algebraic[84])/(states[0]*algebraic[107]) algebraic[112] = (constants[11]*constants[131]*constants[134])/(algebraic[111]*constants[128]*constants[133]) algebraic[55] = 1.00000+algebraic[38]/constants[136]+algebraic[43]/constants[137]+algebraic[52]/constants[138] algebraic[147] = (constants[11]*algebraic[112]*(algebraic[143]*states[24]-(states[14]*states[3])/algebraic[111]))/(algebraic[112]*constants[128]*constants[133]*algebraic[55]+algebraic[112]*constants[133]*algebraic[143]+algebraic[112]*constants[132]*algebraic[55]*states[24]+(constants[11]*constants[135]*algebraic[55]*states[14])/algebraic[111]+(constants[11]*constants[134]*states[3])/algebraic[111]+algebraic[112]*algebraic[143]*states[24]+(constants[11]*constants[135]*algebraic[143]*states[14])/(algebraic[111]*constants[128])+(constants[11]*states[14]*states[3])/algebraic[111]+(algebraic[112]*constants[132]*states[24]*states[3])/constants[131]+(algebraic[112]*algebraic[143]*states[24]*states[14])/constants[130]+(constants[11]*states[24]*states[14]*states[3])/(constants[129]*algebraic[111])) algebraic[0] = exp(-(constants[53]+4.00000*constants[33]*states[2])/constants[249]) algebraic[3] = ((algebraic[0]*(power(states[0], 5.00000)))/(power(states[1], 4.00000)))/1.00000 algebraic[156] = constants[46]-states[4] algebraic[157] = constants[21]*(algebraic[3]*states[3]*algebraic[156]-algebraic[143]*states[4]) algebraic[105] = (constants[109]*algebraic[104])/algebraic[99] algebraic[106] = (constants[9]*constants[114]*constants[111])/(algebraic[105]*constants[110]*constants[113]) algebraic[11] = (1.00000+states[14]/constants[116]+states[21]/constants[117]+states[23]/constants[118])/(1.00000+states[21]/constants[117]+states[23]/constants[118]) algebraic[158] = (constants[9]*algebraic[106]*(states[21]*algebraic[156]-(states[4]*states[23])/algebraic[105]))/(algebraic[106]*constants[110]*constants[113]*algebraic[11]+algebraic[106]*constants[113]*states[21]+algebraic[106]*constants[112]*algebraic[11]*algebraic[156]+(constants[9]*constants[115]*algebraic[11]*states[4])/algebraic[105]+(constants[9]*constants[114]*states[23])/algebraic[105]+algebraic[106]*states[21]*algebraic[156]+(constants[9]*constants[115]*states[21]*states[4])/(algebraic[105]*constants[110])+(algebraic[106]*constants[112]*algebraic[156]*states[23])/constants[111]+(constants[9]*states[4]*states[23])/algebraic[105]) algebraic[4] = exp(-(constants[54]+2.00000*constants[33]*states[2])/constants[249]) algebraic[5] = (1.00000*1.00000*algebraic[4]*(power(states[0], 2.00000)))/(power(states[1], 4.00000)) algebraic[163] = constants[47]-states[6] algebraic[164] = ((constants[22]*(1.00000+states[5]/constants[23]))/(1.00000+states[5]/constants[24]))*((power(algebraic[5], 0.500000))*algebraic[163]*(power(states[4], 0.500000))-states[6]*(power(algebraic[156], 0.500000))) algebraic[6] = exp(-(constants[55]+4.00000*constants[33]*states[2])/constants[249]) algebraic[7] = ((algebraic[6]*(power(states[0], 4.00000)))/(power(states[1], 2.00000)))/(1.00000*1.00000) algebraic[165] = ((((constants[25]*states[7])/(states[7]+constants[44]))*exp((constants[33]*states[2])/constants[249])*states[6])/constants[47])*((power(algebraic[7], 1.00000/2.00000))*states[6]*(power(states[7], 1.00000/4.00000))-algebraic[163]) algebraic[32] = custom_piecewise([greater(fabs(states[2]) , constants[172]), (constants[31]*states[2]*(states[1]*exp((constants[33]*states[2])/constants[249])-states[0]))/(exp((constants[33]*states[2])/constants[249])-1.00000) , True, (constants[31]*constants[249]*(states[1]-states[0]))/constants[33]]) algebraic[168] = -((states[0]*states[0])/(constants[229]*constants[230]*(power(constants[283], 2.00000)))+(states[0]*states[8])/(constants[296]*constants[300]*(power(algebraic[36], 2.00000)))+(states[0]*states[9])/(constants[304]*constants[306]*(power(algebraic[41], 2.00000)))+(states[0]*states[16])/(constants[310]*constants[312]*(power(algebraic[50], 2.00000)))+(states[0]*states[22])/(constants[296]*constants[300]*(power(algebraic[66], 2.00000)))+(states[0]*states[20])/(constants[304]*constants[306]*(power(algebraic[72], 2.00000)))+(states[0]*states[5])/(constants[250]*constants[251]*(power(algebraic[76], 2.00000)))+(states[0]*states[3])/(constants[205]*constants[206]*(power(constants[275], 2.00000)))+(states[0]*algebraic[143])/(constants[202]*constants[203]*(power(constants[274], 2.00000)))+(states[0]*states[4])/(constants[208]*constants[209]*(power(constants[276], 2.00000)))+(states[0]*algebraic[156])/(constants[211]*constants[212]*(power(constants[277], 2.00000)))+(states[0]*states[10])/(constants[198]*constants[269]*(power(algebraic[120], 2.00000)))+(states[0]*states[14])/(constants[189]*constants[254]*(power(algebraic[84], 2.00000)))+(states[0]*states[13])/(constants[214]*constants[215]*(power(constants[278], 2.00000)))+(states[0]*states[15])/(constants[255]*constants[256]*(power(algebraic[85], 2.00000)))+(states[0]*states[18])/(constants[258]*constants[259]*(power(algebraic[94], 2.00000)))+(states[0]*states[19])/(constants[217]*constants[218]*(power(constants[279], 2.00000)))+(states[0]*states[17])/(constants[260]*constants[192]*(power(algebraic[98], 2.00000)))+(states[0]*states[11])/(constants[253]*constants[187]*(power(algebraic[83], 2.00000)))+(states[0]*states[21])/(constants[261]*constants[262]*(power(algebraic[99], 2.00000)))+(states[0]*states[23])/(constants[264]*constants[194]*(power(algebraic[104], 2.00000)))+(states[0]*states[24])/(constants[265]*constants[266]*(power(algebraic[107], 2.00000)))+(states[0]*states[26])/(constants[270]*constants[271]*(power(algebraic[124], 2.00000)))+(states[0]*states[25])/(constants[272]*constants[273]*(power(algebraic[127], 2.00000)))+(states[0]*states[54])/(constants[235]*constants[236]*(power(constants[285], 2.00000)))+(states[0]*states[58])/(constants[232]*constants[233]*(power(constants[284], 2.00000)))+(states[0]*states[7])/(constants[220]*constants[221]*(power(constants[280], 2.00000)))+(states[0]*constants[294])/(constants[223]*constants[224]*(power(constants[281], 2.00000)))+(states[0]*constants[286])/(constants[226]*constants[227]*(power(constants[282], 2.00000)))+(states[0]*states[12])/(constants[268]*constants[196]*(power(algebraic[117], 2.00000)))) algebraic[169] = -((states[0]*states[0])/(constants[229]*constants[231]*(power(constants[283], 2.00000)))+(states[0]*states[8])/(constants[296]*constants[302]*(power(algebraic[36], 2.00000)))+(states[0]*states[9])/(constants[304]*constants[308]*(power(algebraic[41], 2.00000)))+(states[0]*states[16])/(constants[310]*constants[314]*(power(algebraic[50], 2.00000)))+(states[0]*states[22])/(constants[296]*constants[302]*(power(algebraic[66], 2.00000)))+(states[0]*states[20])/(constants[304]*constants[308]*(power(algebraic[72], 2.00000)))+(states[0]*states[5])/(constants[250]*constants[252]*(power(algebraic[76], 2.00000)))+(states[0]*states[3])/(constants[205]*constants[207]*(power(constants[275], 2.00000)))+(states[0]*algebraic[143])/(constants[202]*constants[204]*(power(constants[274], 2.00000)))+(states[0]*states[4])/(constants[208]*constants[210]*(power(constants[276], 2.00000)))+(states[0]*algebraic[156])/(constants[211]*constants[213]*(power(constants[277], 2.00000)))+(states[0]*states[10])/(constants[198]*constants[199]*(power(algebraic[120], 2.00000)))+(states[0]*states[14])/(constants[189]*constants[190]*(power(algebraic[84], 2.00000)))+(states[0]*states[13])/(constants[214]*constants[216]*(power(constants[278], 2.00000)))+(states[0]*states[15])/(constants[255]*constants[257]*(power(algebraic[85], 2.00000)))+(states[0]*states[18])/(constants[258]*constants[191]*(power(algebraic[94], 2.00000)))+(states[0]*states[19])/(constants[217]*constants[219]*(power(constants[279], 2.00000)))+(states[0]*states[17])/(constants[260]*constants[193]*(power(algebraic[98], 2.00000)))+(states[0]*states[11])/(constants[253]*constants[188]*(power(algebraic[83], 2.00000)))+(states[0]*states[21])/(constants[261]*constants[263]*(power(algebraic[99], 2.00000)))+(states[0]*states[23])/(constants[264]*constants[195]*(power(algebraic[104], 2.00000)))+(states[0]*states[24])/(constants[265]*constants[267]*(power(algebraic[107], 2.00000)))+(states[0]*states[26])/(constants[270]*constants[200]*(power(algebraic[124], 2.00000)))+(states[0]*states[25])/(constants[272]*constants[201]*(power(algebraic[127], 2.00000)))+(states[0]*states[54])/(constants[235]*constants[237]*(power(constants[285], 2.00000)))+(states[0]*states[58])/(constants[232]*constants[234]*(power(constants[284], 2.00000)))+(states[0]*states[7])/(constants[220]*constants[222]*(power(constants[280], 2.00000)))+(states[0]*constants[294])/(constants[223]*constants[225]*(power(constants[281], 2.00000)))+(states[0]*constants[286])/(constants[226]*constants[228]*(power(constants[282], 2.00000)))+(states[0]*states[12])/(constants[268]*constants[197]*(power(algebraic[117], 2.00000)))) algebraic[172] = -((states[58]*states[0])/(constants[230]*constants[231]*(power(constants[283], 2.00000)))+(states[58]*states[8])/(constants[300]*constants[302]*(power(algebraic[36], 2.00000)))+(states[58]*states[9])/(constants[306]*constants[308]*(power(algebraic[41], 2.00000)))+(states[58]*states[16])/(constants[312]*constants[314]*(power(algebraic[50], 2.00000)))+(states[58]*states[22])/(constants[300]*constants[302]*(power(algebraic[66], 2.00000)))+(states[58]*states[20])/(constants[306]*constants[308]*(power(algebraic[72], 2.00000)))+(states[58]*states[5])/(constants[251]*constants[252]*(power(algebraic[76], 2.00000)))+(states[58]*states[3])/(constants[206]*constants[207]*(power(constants[275], 2.00000)))+(states[58]*algebraic[143])/(constants[203]*constants[204]*(power(constants[274], 2.00000)))+(states[58]*states[4])/(constants[209]*constants[210]*(power(constants[276], 2.00000)))+(states[58]*algebraic[156])/(constants[212]*constants[213]*(power(constants[277], 2.00000)))+(states[58]*states[10])/(constants[269]*constants[199]*(power(algebraic[120], 2.00000)))+(states[58]*states[14])/(constants[254]*constants[190]*(power(algebraic[84], 2.00000)))+(states[58]*states[13])/(constants[215]*constants[216]*(power(constants[278], 2.00000)))+(states[58]*states[15])/(constants[256]*constants[257]*(power(algebraic[85], 2.00000)))+(states[58]*states[18])/(constants[259]*constants[191]*(power(algebraic[94], 2.00000)))+(states[58]*states[19])/(constants[218]*constants[219]*(power(constants[279], 2.00000)))+(states[58]*states[17])/(constants[192]*constants[193]*(power(algebraic[98], 2.00000)))+(states[58]*states[11])/(constants[187]*constants[188]*(power(algebraic[83], 2.00000)))+(states[58]*states[21])/(constants[262]*constants[263]*(power(algebraic[99], 2.00000)))+(states[58]*states[23])/(constants[194]*constants[195]*(power(algebraic[104], 2.00000)))+(states[58]*states[24])/(constants[266]*constants[267]*(power(algebraic[107], 2.00000)))+(states[58]*states[26])/(constants[271]*constants[200]*(power(algebraic[124], 2.00000)))+(states[58]*states[25])/(constants[273]*constants[201]*(power(algebraic[127], 2.00000)))+(states[58]*states[54])/(constants[236]*constants[237]*(power(constants[285], 2.00000)))+(states[58]*states[58])/(constants[233]*constants[234]*(power(constants[284], 2.00000)))+(states[58]*states[7])/(constants[221]*constants[222]*(power(constants[280], 2.00000)))+(states[58]*constants[294])/(constants[224]*constants[225]*(power(constants[281], 2.00000)))+(states[58]*constants[286])/(constants[227]*constants[228]*(power(constants[282], 2.00000)))+(states[58]*states[12])/(constants[196]*constants[197]*(power(algebraic[117], 2.00000)))) algebraic[175] = -((states[54]*states[0])/(constants[231]*constants[230]*(power(constants[283], 2.00000)))+(states[54]*states[8])/(constants[302]*constants[300]*(power(algebraic[36], 2.00000)))+(states[54]*states[9])/(constants[308]*constants[306]*(power(algebraic[41], 2.00000)))+(states[54]*states[16])/(constants[314]*constants[312]*(power(algebraic[50], 2.00000)))+(states[54]*states[22])/(constants[302]*constants[300]*(power(algebraic[66], 2.00000)))+(states[54]*states[20])/(constants[308]*constants[306]*(power(algebraic[72], 2.00000)))+(states[54]*states[5])/(constants[252]*constants[251]*(power(algebraic[76], 2.00000)))+(states[54]*states[3])/(constants[207]*constants[206]*(power(constants[275], 2.00000)))+(states[54]*algebraic[143])/(constants[204]*constants[203]*(power(constants[274], 2.00000)))+(states[54]*states[4])/(constants[210]*constants[209]*(power(constants[276], 2.00000)))+(states[54]*algebraic[156])/(constants[213]*constants[212]*(power(constants[277], 2.00000)))+(states[54]*states[10])/(constants[199]*constants[269]*(power(algebraic[120], 2.00000)))+(states[54]*states[14])/(constants[190]*constants[254]*(power(algebraic[84], 2.00000)))+(states[54]*states[13])/(constants[216]*constants[215]*(power(constants[278], 2.00000)))+(states[54]*states[15])/(constants[257]*constants[256]*(power(algebraic[85], 2.00000)))+(states[54]*states[18])/(constants[191]*constants[259]*(power(algebraic[94], 2.00000)))+(states[54]*states[19])/(constants[219]*constants[218]*(power(constants[279], 2.00000)))+(states[54]*states[17])/(constants[193]*constants[192]*(power(algebraic[98], 2.00000)))+(states[54]*states[11])/(constants[188]*constants[187]*(power(algebraic[83], 2.00000)))+(states[54]*states[21])/(constants[263]*constants[262]*(power(algebraic[99], 2.00000)))+(states[54]*states[23])/(constants[195]*constants[194]*(power(algebraic[104], 2.00000)))+(states[54]*states[24])/(constants[267]*constants[266]*(power(algebraic[107], 2.00000)))+(states[54]*states[26])/(constants[200]*constants[271]*(power(algebraic[124], 2.00000)))+(states[54]*states[25])/(constants[201]*constants[273]*(power(algebraic[127], 2.00000)))+(states[54]*states[54])/(constants[237]*constants[236]*(power(constants[285], 2.00000)))+(states[54]*states[58])/(constants[234]*constants[233]*(power(constants[284], 2.00000)))+(states[54]*states[7])/(constants[222]*constants[221]*(power(constants[280], 2.00000)))+(states[54]*constants[294])/(constants[225]*constants[224]*(power(constants[281], 2.00000)))+(states[54]*constants[286])/(constants[228]*constants[227]*(power(constants[282], 2.00000)))+(states[54]*states[12])/(constants[197]*constants[196]*(power(algebraic[117], 2.00000)))) algebraic[176] = ((1.00000+states[0]/constants[229]+states[58]/constants[230])*states[0])/(constants[231]*(power(constants[283], 2.00000)))+((1.00000+states[0]/constants[296]+states[58]/constants[300])*states[8])/(constants[302]*(power(algebraic[36], 2.00000)))+((1.00000+states[0]/constants[304]+states[58]/constants[306])*states[9])/(constants[308]*(power(algebraic[41], 2.00000)))+((1.00000+states[0]/constants[310]+states[58]/constants[312])*states[16])/(constants[314]*(power(algebraic[50], 2.00000)))+((1.00000+states[0]/constants[296]+states[58]/constants[300])*states[22])/(constants[302]*(power(algebraic[66], 2.00000)))+((1.00000+states[0]/constants[304]+states[58]/constants[306])*states[20])/(constants[308]*(power(algebraic[72], 2.00000)))+((1.00000+states[0]/constants[250]+states[58]/constants[251])*states[5])/(constants[252]*(power(algebraic[76], 2.00000)))+((1.00000+states[0]/constants[205]+states[58]/constants[206])*states[3])/(constants[207]*(power(constants[275], 2.00000)))+((1.00000+states[0]/constants[202]+states[58]/constants[203])*algebraic[143])/(constants[204]*(power(constants[274], 2.00000)))+((1.00000+states[0]/constants[208]+states[58]/constants[209])*states[4])/(constants[210]*(power(constants[276], 2.00000)))+((1.00000+states[0]/constants[211]+states[58]/constants[212])*algebraic[156])/(constants[213]*(power(constants[277], 2.00000)))+((1.00000+states[0]/constants[198]+states[58]/constants[269])*states[10])/(constants[199]*(power(algebraic[120], 2.00000)))+((1.00000+states[0]/constants[189]+states[58]/constants[254])*states[14])/(constants[190]*(power(algebraic[84], 2.00000)))+((1.00000+states[0]/constants[214]+states[58]/constants[215])*states[13])/(constants[216]*(power(constants[278], 2.00000)))+((1.00000+states[0]/constants[255]+states[58]/constants[256])*states[15])/(constants[257]*(power(algebraic[85], 2.00000)))+((1.00000+states[0]/constants[258]+states[58]/constants[259])*states[18])/(constants[191]*(power(algebraic[94], 2.00000)))+((1.00000+states[0]/constants[217]+states[58]/constants[218])*states[19])/(constants[219]*(power(constants[279], 2.00000)))+((1.00000+states[0]/constants[260]+states[58]/constants[192])*states[17])/(constants[193]*(power(algebraic[98], 2.00000)))+((1.00000+states[0]/constants[253]+states[58]/constants[187])*states[11])/(constants[188]*(power(algebraic[83], 2.00000)))+((1.00000+states[0]/constants[261]+states[58]/constants[262])*states[21])/(constants[263]*(power(algebraic[99], 2.00000)))+((1.00000+states[0]/constants[264]+states[58]/constants[194])*states[23])/(constants[195]*(power(algebraic[104], 2.00000)))+((1.00000+states[0]/constants[265]+states[58]/constants[266])*states[24])/(constants[267]*(power(algebraic[107], 2.00000)))+((1.00000+states[0]/constants[270]+states[58]/constants[271])*states[26])/(constants[200]*(power(algebraic[124], 2.00000)))+((1.00000+states[0]/constants[272]+states[58]/constants[273])*states[25])/(constants[201]*(power(algebraic[127], 2.00000)))+((1.00000+states[0]/constants[235]+states[58]/constants[236])*states[54])/(constants[237]*(power(constants[285], 2.00000)))+((1.00000+states[0]/constants[232]+states[58]/constants[233])*states[58])/(constants[234]*(power(constants[284], 2.00000)))+((1.00000+states[0]/constants[220]+states[58]/constants[221])*states[7])/(constants[222]*(power(constants[280], 2.00000)))+((1.00000+states[0]/constants[223]+states[58]/constants[224])*constants[294])/(constants[225]*(power(constants[281], 2.00000)))+((1.00000+states[0]/constants[226]+states[58]/constants[227])*constants[286])/(constants[228]*(power(constants[282], 2.00000)))+((1.00000+states[0]/constants[268]+states[58]/constants[196])*states[12])/(constants[197]*(power(algebraic[117], 2.00000))) algebraic[185] = 1.00000+algebraic[176] algebraic[173] = ((1.00000+states[0]/constants[229]+states[54]/constants[231])*states[0])/(constants[230]*(power(constants[283], 2.00000)))+((1.00000+states[0]/constants[296]+states[54]/constants[302])*states[8])/(constants[300]*(power(algebraic[36], 2.00000)))+((1.00000+states[0]/constants[304]+states[54]/constants[308])*states[9])/(constants[306]*(power(algebraic[41], 2.00000)))+((1.00000+states[0]/constants[310]+states[54]/constants[314])*states[16])/(constants[312]*(power(algebraic[50], 2.00000)))+((1.00000+states[0]/constants[296]+states[54]/constants[302])*states[22])/(constants[300]*(power(algebraic[66], 2.00000)))+((1.00000+states[0]/constants[304]+states[54]/constants[308])*states[20])/(constants[306]*(power(algebraic[72], 2.00000)))+((1.00000+states[0]/constants[250]+states[54]/constants[252])*states[5])/(constants[251]*(power(algebraic[76], 2.00000)))+((1.00000+states[0]/constants[205]+states[54]/constants[207])*states[3])/(constants[206]*(power(constants[275], 2.00000)))+((1.00000+states[0]/constants[202]+states[54]/constants[204])*algebraic[143])/(constants[203]*(power(constants[274], 2.00000)))+((1.00000+states[0]/constants[208]+states[54]/constants[210])*states[4])/(constants[209]*(power(constants[276], 2.00000)))+((1.00000+states[0]/constants[211]+states[54]/constants[213])*algebraic[156])/(constants[212]*(power(constants[277], 2.00000)))+((1.00000+states[0]/constants[198]+states[54]/constants[199])*states[10])/(constants[269]*(power(algebraic[120], 2.00000)))+((1.00000+states[0]/constants[189]+states[54]/constants[190])*states[14])/(constants[254]*(power(algebraic[84], 2.00000)))+((1.00000+states[0]/constants[214]+states[54]/constants[216])*states[13])/(constants[215]*(power(constants[278], 2.00000)))+((1.00000+states[0]/constants[255]+states[54]/constants[257])*states[15])/(constants[256]*(power(algebraic[85], 2.00000)))+((1.00000+states[0]/constants[258]+states[54]/constants[191])*states[18])/(constants[259]*(power(algebraic[94], 2.00000)))+((1.00000+states[0]/constants[217]+states[54]/constants[219])*states[19])/(constants[218]*(power(constants[279], 2.00000)))+((1.00000+states[0]/constants[260]+states[54]/constants[193])*states[17])/(constants[192]*(power(algebraic[98], 2.00000)))+((1.00000+states[0]/constants[253]+states[54]/constants[188])*states[11])/(constants[187]*(power(algebraic[83], 2.00000)))+((1.00000+states[0]/constants[261]+states[54]/constants[263])*states[21])/(constants[262]*(power(algebraic[99], 2.00000)))+((1.00000+states[0]/constants[264]+states[54]/constants[195])*states[23])/(constants[194]*(power(algebraic[104], 2.00000)))+((1.00000+states[0]/constants[265]+states[54]/constants[267])*states[24])/(constants[266]*(power(algebraic[107], 2.00000)))+((1.00000+states[0]/constants[270]+states[54]/constants[200])*states[26])/(constants[271]*(power(algebraic[124], 2.00000)))+((1.00000+states[0]/constants[272]+states[54]/constants[201])*states[25])/(constants[273]*(power(algebraic[127], 2.00000)))+((1.00000+states[0]/constants[235]+states[54]/constants[237])*states[54])/(constants[236]*(power(constants[285], 2.00000)))+((1.00000+states[0]/constants[232]+states[54]/constants[234])*states[58])/(constants[233]*(power(constants[284], 2.00000)))+((1.00000+states[0]/constants[220]+states[54]/constants[222])*states[7])/(constants[221]*(power(constants[280], 2.00000)))+((1.00000+states[0]/constants[223]+states[54]/constants[225])*constants[294])/(constants[224]*(power(constants[281], 2.00000)))+((1.00000+states[0]/constants[226]+states[54]/constants[228])*constants[286])/(constants[227]*(power(constants[282], 2.00000)))+((1.00000+states[0]/constants[268]+states[54]/constants[197])*states[12])/(constants[196]*(power(algebraic[117], 2.00000))) algebraic[184] = 1.00000+algebraic[173] algebraic[170] = ((1.00000+states[58]/constants[230]+states[54]/constants[231])*states[0])/(constants[229]*(power(constants[283], 2.00000)))+((1.00000+states[58]/constants[300]+states[54]/constants[302])*states[8])/(constants[296]*(power(algebraic[36], 2.00000)))+((1.00000+states[58]/constants[306]+states[54]/constants[308])*states[9])/(constants[304]*(power(algebraic[41], 2.00000)))+((1.00000+states[58]/constants[312]+states[54]/constants[314])*states[16])/(constants[310]*(power(algebraic[50], 2.00000)))+((1.00000+states[58]/constants[300]+states[54]/constants[302])*states[22])/(constants[296]*(power(algebraic[66], 2.00000)))+((1.00000+states[58]/constants[306]+states[54]/constants[308])*states[20])/(constants[304]*(power(algebraic[72], 2.00000)))+((1.00000+states[58]/constants[251]+states[54]/constants[252])*states[5])/(constants[250]*(power(algebraic[76], 2.00000)))+((1.00000+states[58]/constants[206]+states[54]/constants[207])*states[3])/(constants[205]*(power(constants[275], 2.00000)))+((1.00000+states[58]/constants[203]+states[54]/constants[204])*algebraic[143])/(constants[202]*(power(constants[274], 2.00000)))+((1.00000+states[58]/constants[209]+states[54]/constants[210])*states[4])/(constants[208]*(power(constants[276], 2.00000)))+((1.00000+states[58]/constants[212]+states[54]/constants[213])*algebraic[156])/(constants[211]*(power(constants[277], 2.00000)))+((1.00000+states[58]/constants[269]+states[54]/constants[199])*states[10])/(constants[198]*(power(algebraic[120], 2.00000)))+((1.00000+states[58]/constants[254]+states[54]/constants[190])*states[14])/(constants[189]*(power(algebraic[84], 2.00000)))+((1.00000+states[58]/constants[215]+states[54]/constants[216])*states[13])/(constants[214]*(power(constants[278], 2.00000)))+((1.00000+states[58]/constants[256]+states[54]/constants[257])*states[15])/(constants[255]*(power(algebraic[85], 2.00000)))+((1.00000+states[58]/constants[259]+states[54]/constants[191])*states[18])/(constants[258]*(power(algebraic[94], 2.00000)))+((1.00000+states[58]/constants[218]+states[54]/constants[219])*states[19])/(constants[217]*(power(constants[279], 2.00000)))+((1.00000+states[58]/constants[192]+states[54]/constants[193])*states[17])/(constants[260]*(power(algebraic[98], 2.00000)))+((1.00000+states[58]/constants[187]+states[54]/constants[188])*states[11])/(constants[253]*(power(algebraic[83], 2.00000)))+((1.00000+states[58]/constants[262]+states[54]/constants[263])*states[21])/(constants[261]*(power(algebraic[99], 2.00000)))+((1.00000+states[58]/constants[194]+states[54]/constants[195])*states[23])/(constants[264]*(power(algebraic[104], 2.00000)))+((1.00000+states[58]/constants[266]+states[54]/constants[267])*states[24])/(constants[265]*(power(algebraic[107], 2.00000)))+((1.00000+states[58]/constants[271]+states[54]/constants[200])*states[26])/(constants[270]*(power(algebraic[124], 2.00000)))+((1.00000+states[58]/constants[273]+states[54]/constants[201])*states[25])/(constants[272]*(power(algebraic[127], 2.00000)))+((1.00000+states[58]/constants[236]+states[54]/constants[237])*states[54])/(constants[235]*(power(constants[285], 2.00000)))+((1.00000+states[58]/constants[233]+states[54]/constants[234])*states[58])/(constants[232]*(power(constants[284], 2.00000)))+((1.00000+states[58]/constants[221]+states[54]/constants[222])*states[7])/(constants[220]*(power(constants[280], 2.00000)))+((1.00000+states[58]/constants[224]+states[54]/constants[225])*constants[294])/(constants[223]*(power(constants[281], 2.00000)))+((1.00000+states[58]/constants[227]+states[54]/constants[228])*constants[286])/(constants[226]*(power(constants[282], 2.00000)))+((1.00000+states[58]/constants[196]+states[54]/constants[197])*states[12])/(constants[268]*(power(algebraic[117], 2.00000))) algebraic[183] = 1.00000+algebraic[170]+constants[243]/(constants[244]*(power(1.00000+states[0]/constants[244], 2.00000))) algebraic[174] = -((states[54]*states[0])/(constants[231]*constants[229]*(power(constants[283], 2.00000)))+(states[54]*states[8])/(constants[302]*constants[296]*(power(algebraic[36], 2.00000)))+(states[54]*states[9])/(constants[308]*constants[304]*(power(algebraic[41], 2.00000)))+(states[54]*states[16])/(constants[314]*constants[310]*(power(algebraic[50], 2.00000)))+(states[54]*states[22])/(constants[302]*constants[296]*(power(algebraic[66], 2.00000)))+(states[54]*states[20])/(constants[308]*constants[304]*(power(algebraic[72], 2.00000)))+(states[54]*states[5])/(constants[252]*constants[250]*(power(algebraic[76], 2.00000)))+(states[54]*states[3])/(constants[207]*constants[205]*(power(constants[275], 2.00000)))+(states[54]*algebraic[143])/(constants[204]*constants[202]*(power(constants[274], 2.00000)))+(states[54]*states[4])/(constants[210]*constants[208]*(power(constants[276], 2.00000)))+(states[54]*algebraic[156])/(constants[213]*constants[211]*(power(constants[277], 2.00000)))+(states[54]*states[10])/(constants[199]*constants[198]*(power(algebraic[120], 2.00000)))+(states[54]*states[14])/(constants[190]*constants[189]*(power(algebraic[84], 2.00000)))+(states[54]*states[13])/(constants[216]*constants[214]*(power(constants[278], 2.00000)))+(states[54]*states[15])/(constants[257]*constants[255]*(power(algebraic[85], 2.00000)))+(states[54]*states[18])/(constants[191]*constants[258]*(power(algebraic[94], 2.00000)))+(states[54]*states[19])/(constants[219]*constants[217]*(power(constants[279], 2.00000)))+(states[54]*states[17])/(constants[193]*constants[260]*(power(algebraic[98], 2.00000)))+(states[54]*states[11])/(constants[188]*constants[253]*(power(algebraic[83], 2.00000)))+(states[54]*states[21])/(constants[263]*constants[261]*(power(algebraic[99], 2.00000)))+(states[54]*states[23])/(constants[195]*constants[264]*(power(algebraic[104], 2.00000)))+(states[54]*states[24])/(constants[267]*constants[265]*(power(algebraic[107], 2.00000)))+(states[54]*states[26])/(constants[200]*constants[270]*(power(algebraic[124], 2.00000)))+(states[54]*states[25])/(constants[201]*constants[272]*(power(algebraic[127], 2.00000)))+(states[54]*states[54])/(constants[237]*constants[235]*(power(constants[285], 2.00000)))+(states[54]*states[58])/(constants[234]*constants[232]*(power(constants[284], 2.00000)))+(states[54]*states[7])/(constants[222]*constants[220]*(power(constants[280], 2.00000)))+(states[54]*constants[294])/(constants[225]*constants[223]*(power(constants[281], 2.00000)))+(states[54]*constants[286])/(constants[228]*constants[226]*(power(constants[282], 2.00000)))+(states[54]*states[12])/(constants[197]*constants[268]*(power(algebraic[117], 2.00000)))) algebraic[171] = -((states[58]*states[0])/(constants[230]*constants[229]*(power(constants[283], 2.00000)))+(states[58]*states[8])/(constants[300]*constants[296]*(power(algebraic[36], 2.00000)))+(states[58]*states[9])/(constants[306]*constants[304]*(power(algebraic[41], 2.00000)))+(states[58]*states[16])/(constants[312]*constants[310]*(power(algebraic[50], 2.00000)))+(states[58]*states[22])/(constants[300]*constants[296]*(power(algebraic[66], 2.00000)))+(states[58]*states[20])/(constants[306]*constants[304]*(power(algebraic[72], 2.00000)))+(states[58]*states[5])/(constants[251]*constants[250]*(power(algebraic[76], 2.00000)))+(states[58]*states[3])/(constants[206]*constants[205]*(power(constants[275], 2.00000)))+(states[58]*algebraic[143])/(constants[203]*constants[202]*(power(constants[274], 2.00000)))+(states[58]*states[4])/(constants[209]*constants[208]*(power(constants[276], 2.00000)))+(states[58]*algebraic[156])/(constants[212]*constants[211]*(power(constants[277], 2.00000)))+(states[58]*states[10])/(constants[269]*constants[198]*(power(algebraic[120], 2.00000)))+(states[58]*states[14])/(constants[254]*constants[189]*(power(algebraic[84], 2.00000)))+(states[58]*states[13])/(constants[215]*constants[214]*(power(constants[278], 2.00000)))+(states[58]*states[15])/(constants[256]*constants[255]*(power(algebraic[85], 2.00000)))+(states[58]*states[18])/(constants[259]*constants[258]*(power(algebraic[94], 2.00000)))+(states[58]*states[19])/(constants[218]*constants[217]*(power(constants[279], 2.00000)))+(states[58]*states[17])/(constants[192]*constants[260]*(power(algebraic[98], 2.00000)))+(states[58]*states[11])/(constants[187]*constants[253]*(power(algebraic[83], 2.00000)))+(states[58]*states[21])/(constants[262]*constants[261]*(power(algebraic[99], 2.00000)))+(states[58]*states[23])/(constants[194]*constants[264]*(power(algebraic[104], 2.00000)))+(states[58]*states[24])/(constants[266]*constants[265]*(power(algebraic[107], 2.00000)))+(states[58]*states[26])/(constants[271]*constants[270]*(power(algebraic[124], 2.00000)))+(states[58]*states[25])/(constants[273]*constants[272]*(power(algebraic[127], 2.00000)))+(states[58]*states[54])/(constants[236]*constants[235]*(power(constants[285], 2.00000)))+(states[58]*states[58])/(constants[233]*constants[232]*(power(constants[284], 2.00000)))+(states[58]*states[7])/(constants[221]*constants[220]*(power(constants[280], 2.00000)))+(states[58]*constants[294])/(constants[224]*constants[223]*(power(constants[281], 2.00000)))+(states[58]*constants[286])/(constants[227]*constants[226]*(power(constants[282], 2.00000)))+(states[58]*states[12])/(constants[196]*constants[268]*(power(algebraic[117], 2.00000)))) algebraic[186] = (((algebraic[183]*algebraic[175]*algebraic[172]+algebraic[185]*algebraic[168]*algebraic[171]+algebraic[184]*algebraic[169]*algebraic[174])-algebraic[184]*algebraic[185]*algebraic[183])-algebraic[169]*algebraic[175]*algebraic[171])-algebraic[168]*algebraic[172]*algebraic[174] algebraic[135] = ((algebraic[58]+algebraic[78])-algebraic[63])/constants[34] algebraic[136] = ((-algebraic[58]-algebraic[78])+algebraic[63])/constants[34] algebraic[137] = (algebraic[102]-algebraic[58])/constants[34] algebraic[138] = (-algebraic[102]+algebraic[58])/constants[34] algebraic[139] = (((-algebraic[102]-algebraic[78])+algebraic[82])-algebraic[115])/constants[34] algebraic[159] = ((algebraic[144]+algebraic[145]+algebraic[146]+algebraic[147])-algebraic[157])/constants[34] algebraic[166] = ((algebraic[158]+algebraic[157])-algebraic[164])/constants[34] algebraic[148] = (-algebraic[144]+algebraic[123])/constants[34] algebraic[149] = (algebraic[144]-algebraic[89])/constants[34] algebraic[140] = (-algebraic[97]+algebraic[114]+algebraic[89])/constants[34] algebraic[150] = (algebraic[97]-algebraic[145])/constants[34] algebraic[151] = (((algebraic[145]-algebraic[146])-algebraic[131])+algebraic[13])/constants[34] algebraic[152] = (algebraic[146]-algebraic[102])/constants[34] algebraic[153] = ((-algebraic[144]-algebraic[146])+algebraic[102]+algebraic[89])/constants[34] algebraic[160] = ((algebraic[102]-algebraic[158])+algebraic[116])/constants[34] algebraic[161] = (algebraic[158]-algebraic[110])/constants[34] algebraic[154] = (((((algebraic[110]-algebraic[147])+algebraic[115])-algebraic[13])-algebraic[114])-algebraic[116])/constants[34] algebraic[155] = (-algebraic[89]+algebraic[147]+algebraic[131])/constants[34] algebraic[141] = ((algebraic[131]+algebraic[126])-algebraic[134])/constants[34] algebraic[142] = (-algebraic[131]+algebraic[134])/constants[34] algebraic[162] = -algebraic[159] algebraic[167] = -algebraic[166] algebraic[177] = (((-algebraic[144]+2.00000*algebraic[89])-algebraic[146])+algebraic[102]+algebraic[147])/constants[34] algebraic[34] = constants[30]*(states[53]*states[0]-states[54]*states[1]) algebraic[178] = (((((((algebraic[123]+algebraic[126]+algebraic[114])-algebraic[134])-5.00000*algebraic[157])-2.00000*algebraic[164])-4.00000*algebraic[165])+(constants[38]-1.00000)*algebraic[78]+2.00000*algebraic[82]+algebraic[32])-algebraic[34])/constants[34] algebraic[180] = (states[0]*constants[287])/(constants[229]*constants[283])+(states[0]*algebraic[135])/(constants[296]*algebraic[36])+(states[0]*algebraic[136])/(constants[304]*algebraic[41])+(states[0]*constants[238])/(constants[310]*algebraic[50])+(states[0]*algebraic[137])/(constants[296]*algebraic[66])+(states[0]*algebraic[138])/(constants[304]*algebraic[72])+(states[0]*algebraic[139])/(constants[250]*algebraic[76])+(states[0]*algebraic[159])/(constants[205]*constants[275])+(states[0]*algebraic[162])/(constants[202]*constants[274])+(states[0]*algebraic[166])/(constants[208]*constants[276])+(states[0]*algebraic[167])/(constants[211]*constants[277])+(states[0]*algebraic[148])/(constants[198]*algebraic[120])+(states[0]*algebraic[155])/(constants[189]*algebraic[84])+(states[0]*algebraic[149])/(constants[214]*constants[278])+(states[0]*algebraic[140])/(constants[255]*algebraic[85])+(states[0]*algebraic[150])/(constants[258]*algebraic[94])+(states[0]*algebraic[151])/(constants[217]*constants[279])+(states[0]*algebraic[152])/(constants[260]*algebraic[98])+(states[0]*algebraic[153])/(constants[253]*algebraic[83])+(states[0]*algebraic[160])/(constants[261]*algebraic[99])+(states[0]*algebraic[161])/(constants[264]*algebraic[104])+(states[0]*algebraic[154])/(constants[265]*algebraic[107])+(states[0]*algebraic[141])/(constants[270]*algebraic[124])+(states[0]*algebraic[142])/(constants[272]*algebraic[127])+(states[0]*constants[289])/(constants[235]*constants[285])+(states[0]*constants[288])/(constants[232]*constants[284])+(states[0]*constants[239])/(constants[220]*constants[280])+(states[0]*constants[242])/(constants[223]*constants[281])+(states[0]*constants[241])/(constants[226]*constants[282])+(states[0]*constants[240])/(constants[268]*algebraic[117])+algebraic[177]+algebraic[178] algebraic[181] = (states[58]*constants[287])/(constants[230]*constants[283])+(states[58]*algebraic[135])/(constants[300]*algebraic[36])+(states[58]*algebraic[136])/(constants[306]*algebraic[41])+(states[58]*constants[238])/(constants[312]*algebraic[50])+(states[58]*algebraic[137])/(constants[300]*algebraic[66])+(states[58]*algebraic[138])/(constants[306]*algebraic[72])+(states[58]*algebraic[139])/(constants[251]*algebraic[76])+(states[58]*algebraic[159])/(constants[206]*constants[275])+(states[58]*algebraic[162])/(constants[203]*constants[274])+(states[58]*algebraic[166])/(constants[209]*constants[276])+(states[58]*algebraic[167])/(constants[212]*constants[277])+(states[58]*algebraic[148])/(constants[269]*algebraic[120])+(states[58]*algebraic[155])/(constants[254]*algebraic[84])+(states[58]*algebraic[149])/(constants[215]*constants[278])+(states[58]*algebraic[140])/(constants[256]*algebraic[85])+(states[58]*algebraic[150])/(constants[259]*algebraic[94])+(states[58]*algebraic[151])/(constants[218]*constants[279])+(states[58]*algebraic[152])/(constants[192]*algebraic[98])+(states[58]*algebraic[153])/(constants[187]*algebraic[83])+(states[58]*algebraic[160])/(constants[262]*algebraic[99])+(states[58]*algebraic[161])/(constants[194]*algebraic[104])+(states[58]*algebraic[154])/(constants[266]*algebraic[107])+(states[58]*algebraic[141])/(constants[271]*algebraic[124])+(states[58]*algebraic[142])/(constants[273]*algebraic[127])+(states[58]*constants[289])/(constants[236]*constants[285])+(states[58]*constants[288])/(constants[233]*constants[284])+(states[58]*constants[239])/(constants[221]*constants[280])+(states[58]*constants[242])/(constants[224]*constants[281])+(states[58]*constants[241])/(constants[227]*constants[282])+(states[58]*constants[240])/(constants[196]*algebraic[117])+constants[290] algebraic[179] = algebraic[34]/constants[34] algebraic[182] = (states[54]*constants[287])/(constants[231]*constants[283])+(states[54]*algebraic[135])/(constants[302]*algebraic[36])+(states[54]*algebraic[136])/(constants[308]*algebraic[41])+(states[54]*constants[238])/(constants[314]*algebraic[50])+(states[54]*algebraic[137])/(constants[302]*algebraic[66])+(states[54]*algebraic[138])/(constants[308]*algebraic[72])+(states[54]*algebraic[139])/(constants[252]*algebraic[76])+(states[54]*algebraic[159])/(constants[207]*constants[275])+(states[54]*algebraic[162])/(constants[204]*constants[274])+(states[54]*algebraic[166])/(constants[210]*constants[276])+(states[54]*algebraic[167])/(constants[213]*constants[277])+(states[54]*algebraic[148])/(constants[199]*algebraic[120])+(states[54]*algebraic[155])/(constants[190]*algebraic[84])+(states[54]*algebraic[149])/(constants[216]*constants[278])+(states[54]*algebraic[140])/(constants[257]*algebraic[85])+(states[54]*algebraic[150])/(constants[191]*algebraic[94])+(states[54]*algebraic[151])/(constants[219]*constants[279])+(states[54]*algebraic[152])/(constants[193]*algebraic[98])+(states[54]*algebraic[153])/(constants[188]*algebraic[83])+(states[54]*algebraic[160])/(constants[263]*algebraic[99])+(states[54]*algebraic[161])/(constants[195]*algebraic[104])+(states[54]*algebraic[154])/(constants[267]*algebraic[107])+(states[54]*algebraic[141])/(constants[200]*algebraic[124])+(states[54]*algebraic[142])/(constants[201]*algebraic[127])+(states[54]*constants[289])/(constants[237]*constants[285])+(states[54]*constants[288])/(constants[234]*constants[284])+(states[54]*constants[239])/(constants[222]*constants[280])+(states[54]*constants[242])/(constants[225]*constants[281])+(states[54]*constants[241])/(constants[228]*constants[282])+(states[54]*constants[240])/(constants[197]*algebraic[117])+algebraic[179] algebraic[40] = states[8]/algebraic[36] algebraic[45] = states[9]/algebraic[41] algebraic[53] = states[16]/algebraic[50] algebraic[61] = states[35]/algebraic[59] algebraic[62] = (states[35]*(1.00000+states[1]/constants[296]))/algebraic[59] algebraic[64] = states[37]/algebraic[60] algebraic[65] = (states[37]*(1.00000+states[1]/constants[304]))/algebraic[60] algebraic[70] = states[39]/algebraic[67] algebraic[71] = (states[39]*(1.00000+states[1]/constants[310]))/algebraic[67] algebraic[93] = (states[43]*(1.00000+states[1]/constants[255]))/algebraic[91] 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)