Generated Code

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

The raw code is available.

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