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