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