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 =207;
end
% There are a total of 60 entries in each of the rate and state variable arrays.
% There are a total of 388 entries in the constant variable array.
%

function [VOI, STATES, ALGEBRAIC, CONSTANTS] = solveModel()
    % Create ALGEBRAIC of correct size
    global algebraicVariableCount;  algebraicVariableCount = getAlgebraicVariableCount();
    % Initialise constants and state variables
    [INIT_STATES, CONSTANTS] = initConsts;

    % Set timespan to solve over 
    tspan = [0, 10];

    % Set numerical accuracy options for ODE solver
    options = odeset('RelTol', 1e-06, 'AbsTol', 1e-06, 'MaxStep', 1);

    % Solve model with ODE solver
    [VOI, STATES] = ode15s(@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS), tspan, INIT_STATES, options);

    % Compute algebraic variables
    [RATES, ALGEBRAIC] = computeRates(VOI, STATES, CONSTANTS);
    ALGEBRAIC = computeAlgebraic(ALGEBRAIC, CONSTANTS, STATES, VOI);

    % Plot state variables against variable of integration
    [LEGEND_STATES, LEGEND_ALGEBRAIC, LEGEND_VOI, LEGEND_CONSTANTS] = createLegends();
    figure();
    plot(VOI, STATES);
    xlabel(LEGEND_VOI);
    l = legend(LEGEND_STATES);
    set(l,'Interpreter','none');
end

function [LEGEND_STATES, LEGEND_ALGEBRAIC, LEGEND_VOI, LEGEND_CONSTANTS] = createLegends()
    LEGEND_STATES = ''; LEGEND_ALGEBRAIC = ''; LEGEND_VOI = ''; LEGEND_CONSTANTS = '';
    LEGEND_VOI = strpad('time in component Environment (second)');
    LEGEND_CONSTANTS(:,1) = strpad('ExpType in component Environment (dimensionless)');
    LEGEND_CONSTANTS(:,2) = strpad('StateType in component Environment (dimensionless)');
    LEGEND_CONSTANTS(:,3) = strpad('F in component parameters (kilojoule_per_mole_per_millivolt)');
    LEGEND_CONSTANTS(:,324) = strpad('V_cyto in component parameters (cyto_per_cell)');
    LEGEND_CONSTANTS(:,308) = strpad('V_mito in component parameters (mito_per_cell)');
    LEGEND_CONSTANTS(:,4) = strpad('W_x in component parameters (l_water_per_l_mito)');
    LEGEND_CONSTANTS(:,5) = strpad('W_i in component parameters (l_water_per_l_mito)');
    LEGEND_CONSTANTS(:,287) = strpad('W_c in component parameters (l_water_per_l_cyto)');
    LEGEND_CONSTANTS(:,6) = strpad('gamma in component parameters (per_micrometer)');
    LEGEND_CONSTANTS(:,7) = strpad('rho_m in component parameters (mg_per_mitochondria)');
    LEGEND_CONSTANTS(:,8) = strpad('n_A in component parameters (dimensionless)');
    LEGEND_CONSTANTS(:,9) = strpad('p_PI in component parameters (micrometer_per_second)');
    LEGEND_CONSTANTS(:,10) = strpad('p_A in component parameters (micrometer_per_second)');
    LEGEND_CONSTANTS(:,11) = strpad('p_TI in component parameters (micrometer_per_second)');
    LEGEND_CONSTANTS(:,12) = strpad('k_mADP in component parameters (molar)');
    LEGEND_CONSTANTS(:,13) = strpad('theta in component parameters (dimensionless)');
    LEGEND_CONSTANTS(:,14) = strpad('k_O2 in component parameters (molar)');
    LEGEND_CONSTANTS(:,15) = strpad('NAD_tot in component parameters (molar)');
    LEGEND_CONSTANTS(:,16) = strpad('Q_tot in component parameters (molar)');
    LEGEND_CONSTANTS(:,17) = strpad('cytC_tot in component parameters (molar)');
    LEGEND_CONSTANTS(:,18) = strpad('A_tot in component parameters (molar)');
    LEGEND_CONSTANTS(:,19) = strpad('CR_tot in component parameters (molar)');
    LEGEND_CONSTANTS(:,20) = strpad('CO2_tot_x in component parameters (molar)');
    LEGEND_CONSTANTS(:,21) = strpad('FAD_tot in component parameters (molar)');
    LEGEND_CONSTANTS(:,22) = strpad('C_IM in component parameters (mole_per_l_mito_volume_per_mv)');
    LEGEND_CONSTANTS(:,23) = strpad('X_pdh in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,24) = strpad('X_cits in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,25) = strpad('X_acon in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,26) = strpad('X_isod in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,27) = strpad('X_akgd in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,28) = strpad('K_ir_akgd in component parameters (molar)');
    LEGEND_CONSTANTS(:,29) = strpad('X_scoas in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,30) = strpad('X_sdh in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,31) = strpad('X_fum in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,32) = strpad('X_mdh in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,33) = strpad('X_ndk in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,34) = strpad('X_got in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,35) = strpad('X_PYRH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,36) = strpad('X_GLUH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,37) = strpad('X_CITMAL in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,38) = strpad('X_AKGMAL in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,39) = strpad('X_SUCMAL in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,40) = strpad('X_MALPI in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,41) = strpad('X_ASPGLU in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,42) = strpad('X_C1 in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,43) = strpad('X_C3 in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_half_molar)');
    LEGEND_CONSTANTS(:,44) = strpad('k_PI_1 in component parameters (molar)');
    LEGEND_CONSTANTS(:,45) = strpad('k_PI_2 in component parameters (molar)');
    LEGEND_CONSTANTS(:,46) = strpad('X_C4 in component parameters (mole_per_second_per_l_mito_volume_per_molar)');
    LEGEND_CONSTANTS(:,47) = strpad('X_F1 in component parameters (mole_per_second_per_l_mito_volume_per_molar)');
    LEGEND_CONSTANTS(:,48) = strpad('X_ANT in component parameters (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,49) = strpad('X_PiHt in component parameters (mole_per_second_per_l_mito_volume_per_molar)');
    LEGEND_CONSTANTS(:,50) = strpad('k_PiHt in component parameters (molar)');
    LEGEND_CONSTANTS(:,51) = strpad('X_KH in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,52) = strpad('X_Hle in component parameters (mole_per_second_per_l_mito_volume_per_molar_per_mv)');
    LEGEND_CONSTANTS(:,274) = strpad('X_HK in component parameters (mole_per_second_per_l_cyto_volume)');
    LEGEND_CONSTANTS(:,53) = strpad('J_AtC in component parameters (mole_per_second_per_l_cyto_volume)');
    LEGEND_CONSTANTS(:,275) = strpad('RT in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,54) = strpad('dGf1_H2O in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,55) = strpad('dGf1_O2 in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,56) = strpad('dGf1_NADH in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,57) = strpad('dGf1_NAD in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,58) = strpad('dGf1_QH2 in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,59) = strpad('dGf1_COQ in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,60) = strpad('dGf1_ATP in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,61) = strpad('dGf1_ADP in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,62) = strpad('dGf1_AMP in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,63) = strpad('dGf1_GTP in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,64) = strpad('dGf1_GDP in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,65) = strpad('dGf1_Cred in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,66) = strpad('dGf1_Cox in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,67) = strpad('dGf1_Pi in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,68) = strpad('dGf1_PCr in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,69) = strpad('dGf1_Cr in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,70) = strpad('dGf1_FADH2 in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,71) = strpad('dGf1_FAD in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,72) = strpad('dGf1_COASH in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,73) = strpad('dGf1_ACCOA in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,74) = strpad('dGf1_OAA in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,75) = strpad('dGf1_CIT in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,76) = strpad('dGf1_ICIT in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,77) = strpad('dGf1_AKG in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,78) = strpad('dGf1_SCOA in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,79) = strpad('dGf1_SUC in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,80) = strpad('dGf1_FUM in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,81) = strpad('dGf1_MAL in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,82) = strpad('dGf1_ASP in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,83) = strpad('dGf1_GLU in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,84) = strpad('dGf1_CO2_tot in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,85) = strpad('dGf1_PYR in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,86) = strpad('dGf1_GLC in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,87) = strpad('dGf1_G6P in component thermodynamic_data (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,1) = strpad('P_ATP_x in component P_ATP (dimensionless)');
    LEGEND_ALGEBRAIC(:,164) = strpad('P_ATP_i in component P_ATP (dimensionless)');
    LEGEND_ALGEBRAIC(:,18) = strpad('P_ATP_c in component P_ATP (dimensionless)');
    LEGEND_ALGEBRAIC(:,19) = strpad('uATP_x in component P_ATP (molar)');
    LEGEND_ALGEBRAIC(:,21) = strpad('mATP_c in component P_ATP (molar)');
    LEGEND_CONSTANTS(:,297) = strpad('K_H in component P_ATP (molar)');
    LEGEND_CONSTANTS(:,312) = strpad('K_Mg in component P_ATP (molar)');
    LEGEND_CONSTANTS(:,326) = strpad('K_K in component P_ATP (molar)');
    LEGEND_STATES(:,1) = strpad('ATP_x in component dATP_x_dt (molar)');
    LEGEND_STATES(:,2) = strpad('H_x in component dH_x_dt (molar)');
    LEGEND_STATES(:,3) = strpad('Mg_x in component dMg_x_dt (molar)');
    LEGEND_STATES(:,4) = strpad('K_x in component dK_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,133) = strpad('H_i in component dH_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,150) = strpad('Mg_i in component dMg_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,163) = strpad('K_i in component dK_i_dt (molar)');
    LEGEND_STATES(:,5) = strpad('ATP_c in component dATP_c_dt (molar)');
    LEGEND_STATES(:,6) = strpad('H_c in component dH_c_dt (molar)');
    LEGEND_STATES(:,7) = strpad('Mg_c in component dMg_c_dt (molar)');
    LEGEND_STATES(:,8) = strpad('K_c in component dK_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,22) = strpad('P_ADP_x in component P_ADP (dimensionless)');
    LEGEND_ALGEBRAIC(:,166) = strpad('P_ADP_i in component P_ADP (dimensionless)');
    LEGEND_ALGEBRAIC(:,23) = strpad('P_ADP_c in component P_ADP (dimensionless)');
    LEGEND_ALGEBRAIC(:,24) = strpad('uADP_x in component P_ADP (molar)');
    LEGEND_ALGEBRAIC(:,26) = strpad('mADP_c in component P_ADP (molar)');
    LEGEND_CONSTANTS(:,333) = strpad('K_H in component P_ADP (molar)');
    LEGEND_CONSTANTS(:,337) = strpad('K_Mg in component P_ADP (molar)');
    LEGEND_CONSTANTS(:,341) = strpad('K_K in component P_ADP (molar)');
    LEGEND_STATES(:,9) = strpad('ADP_x in component dADP_x_dt (molar)');
    LEGEND_STATES(:,10) = strpad('ADP_c in component dADP_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,27) = strpad('P_AMP_x in component P_AMP (dimensionless)');
    LEGEND_ALGEBRAIC(:,168) = strpad('P_AMP_i in component P_AMP (dimensionless)');
    LEGEND_ALGEBRAIC(:,28) = strpad('P_AMP_c in component P_AMP (dimensionless)');
    LEGEND_ALGEBRAIC(:,29) = strpad('uAMP_x in component P_AMP (molar)');
    LEGEND_CONSTANTS(:,345) = strpad('K_H in component P_AMP (molar)');
    LEGEND_CONSTANTS(:,348) = strpad('K_Mg in component P_AMP (molar)');
    LEGEND_CONSTANTS(:,350) = strpad('K_K in component P_AMP (molar)');
    LEGEND_STATES(:,11) = strpad('AMP_x in component dAMP_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,30) = strpad('P_GTP_x in component P_GTP (dimensionless)');
    LEGEND_ALGEBRAIC(:,165) = strpad('P_GTP_i in component P_GTP (dimensionless)');
    LEGEND_ALGEBRAIC(:,20) = strpad('P_GTP_c in component P_GTP (dimensionless)');
    LEGEND_ALGEBRAIC(:,31) = strpad('uGTP_x in component P_GTP (molar)');
    LEGEND_CONSTANTS(:,314) = strpad('K_H in component P_GTP (molar)');
    LEGEND_CONSTANTS(:,327) = strpad('K_Mg in component P_GTP (molar)');
    LEGEND_CONSTANTS(:,334) = strpad('K_K in component P_GTP (molar)');
    LEGEND_STATES(:,12) = strpad('GTP_x in component dGTP_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,32) = strpad('P_GDP_x in component P_GDP (dimensionless)');
    LEGEND_ALGEBRAIC(:,167) = strpad('P_GDP_i in component P_GDP (dimensionless)');
    LEGEND_ALGEBRAIC(:,25) = strpad('P_GDP_c in component P_GDP (dimensionless)');
    LEGEND_ALGEBRAIC(:,33) = strpad('uGDP_x in component P_GDP (molar)');
    LEGEND_CONSTANTS(:,338) = strpad('K_H in component P_GDP (molar)');
    LEGEND_CONSTANTS(:,342) = strpad('K_Mg in component P_GDP (molar)');
    LEGEND_CONSTANTS(:,346) = strpad('K_K in component P_GDP (molar)');
    LEGEND_STATES(:,13) = strpad('GDP_x in component dGDP_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,34) = strpad('P_PI_x in component P_PI (dimensionless)');
    LEGEND_ALGEBRAIC(:,169) = strpad('P_PI_i in component P_PI (dimensionless)');
    LEGEND_ALGEBRAIC(:,2) = strpad('P_PI_c in component P_PI (dimensionless)');
    LEGEND_ALGEBRAIC(:,35) = strpad('hPI_x in component P_PI (molar)');
    LEGEND_ALGEBRAIC(:,170) = strpad('hPI_i in component P_PI (molar)');
    LEGEND_STATES(:,14) = strpad('PI_x in component dPI_x_dt (molar)');
    LEGEND_STATES(:,15) = strpad('PI_i in component dPI_i_dt (molar)');
    LEGEND_CONSTANTS(:,352) = strpad('K_H in component P_PI (molar)');
    LEGEND_CONSTANTS(:,354) = strpad('K_Mg in component P_PI (molar)');
    LEGEND_CONSTANTS(:,356) = strpad('K_K in component P_PI (molar)');
    LEGEND_ALGEBRAIC(:,36) = strpad('P_COASH_x in component P_COASH (dimensionless)');
    LEGEND_ALGEBRAIC(:,134) = strpad('P_COASH_i in component P_COASH (dimensionless)');
    LEGEND_ALGEBRAIC(:,3) = strpad('P_COASH_c in component P_COASH (dimensionless)');
    LEGEND_CONSTANTS(:,358) = strpad('K_H in component P_COASH (molar)');
    LEGEND_CONSTANTS(:,88) = strpad('K_Mg in component P_COASH (molar)');
    LEGEND_CONSTANTS(:,89) = strpad('K_K in component P_COASH (molar)');
    LEGEND_ALGEBRAIC(:,37) = strpad('P_OAA_x in component P_OAA (dimensionless)');
    LEGEND_ALGEBRAIC(:,151) = strpad('P_OAA_i in component P_OAA (dimensionless)');
    LEGEND_ALGEBRAIC(:,4) = strpad('P_OAA_c in component P_OAA (dimensionless)');
    LEGEND_CONSTANTS(:,360) = strpad('K_Mg in component P_OAA (molar)');
    LEGEND_CONSTANTS(:,90) = strpad('K_H in component P_OAA (molar)');
    LEGEND_CONSTANTS(:,91) = strpad('K_K in component P_OAA (molar)');
    LEGEND_ALGEBRAIC(:,38) = strpad('P_CIT_x in component P_CIT (dimensionless)');
    LEGEND_ALGEBRAIC(:,171) = strpad('P_CIT_i in component P_CIT (dimensionless)');
    LEGEND_ALGEBRAIC(:,5) = strpad('P_CIT_c in component P_CIT (dimensionless)');
    LEGEND_ALGEBRAIC(:,39) = strpad('uCIT_x in component P_CIT (molar)');
    LEGEND_ALGEBRAIC(:,172) = strpad('uCIT_i in component P_CIT (molar)');
    LEGEND_ALGEBRAIC(:,40) = strpad('hCIT_x in component P_CIT (molar)');
    LEGEND_ALGEBRAIC(:,173) = strpad('hCIT_i in component P_CIT (molar)');
    LEGEND_STATES(:,16) = strpad('CIT_x in component dCIT_x_dt (molar)');
    LEGEND_STATES(:,17) = strpad('CIT_i in component dCIT_i_dt (molar)');
    LEGEND_CONSTANTS(:,361) = strpad('K_H in component P_CIT (molar)');
    LEGEND_CONSTANTS(:,362) = strpad('K_Mg in component P_CIT (molar)');
    LEGEND_CONSTANTS(:,363) = strpad('K_K in component P_CIT (molar)');
    LEGEND_ALGEBRAIC(:,41) = strpad('P_ICIT_x in component P_ICIT (dimensionless)');
    LEGEND_ALGEBRAIC(:,152) = strpad('P_ICIT_i in component P_ICIT (dimensionless)');
    LEGEND_ALGEBRAIC(:,6) = strpad('P_ICIT_c in component P_ICIT (dimensionless)');
    LEGEND_CONSTANTS(:,364) = strpad('K_H in component P_ICIT (molar)');
    LEGEND_CONSTANTS(:,365) = strpad('K_Mg in component P_ICIT (molar)');
    LEGEND_CONSTANTS(:,92) = strpad('K_K in component P_ICIT (molar)');
    LEGEND_ALGEBRAIC(:,42) = strpad('P_SCOA_x in component P_SCOA (dimensionless)');
    LEGEND_ALGEBRAIC(:,135) = strpad('P_SCOA_i in component P_SCOA (dimensionless)');
    LEGEND_ALGEBRAIC(:,7) = strpad('P_SCOA_c in component P_SCOA (dimensionless)');
    LEGEND_CONSTANTS(:,366) = strpad('K_H in component P_SCOA (molar)');
    LEGEND_CONSTANTS(:,93) = strpad('K_Mg in component P_SCOA (molar)');
    LEGEND_CONSTANTS(:,94) = strpad('K_K in component P_SCOA (molar)');
    LEGEND_ALGEBRAIC(:,43) = strpad('P_SUC_x in component P_SUC (dimensionless)');
    LEGEND_ALGEBRAIC(:,174) = strpad('P_SUC_i in component P_SUC (dimensionless)');
    LEGEND_ALGEBRAIC(:,8) = strpad('P_SUC_c in component P_SUC (dimensionless)');
    LEGEND_CONSTANTS(:,367) = strpad('K_H in component P_SUC (molar)');
    LEGEND_CONSTANTS(:,368) = strpad('K_Mg in component P_SUC (molar)');
    LEGEND_CONSTANTS(:,369) = strpad('K_K in component P_SUC (molar)');
    LEGEND_ALGEBRAIC(:,44) = strpad('P_FUM_x in component P_FUM (dimensionless)');
    LEGEND_ALGEBRAIC(:,136) = strpad('P_FUM_i in component P_FUM (dimensionless)');
    LEGEND_ALGEBRAIC(:,9) = strpad('P_FUM_c in component P_FUM (dimensionless)');
    LEGEND_CONSTANTS(:,370) = strpad('K_H in component P_FUM (molar)');
    LEGEND_CONSTANTS(:,95) = strpad('K_Mg in component P_FUM (molar)');
    LEGEND_CONSTANTS(:,96) = strpad('K_K in component P_FUM (molar)');
    LEGEND_ALGEBRAIC(:,45) = strpad('P_MAL_x in component P_MAL (dimensionless)');
    LEGEND_ALGEBRAIC(:,175) = strpad('P_MAL_i in component P_MAL (dimensionless)');
    LEGEND_ALGEBRAIC(:,10) = strpad('P_MAL_c in component P_MAL (dimensionless)');
    LEGEND_CONSTANTS(:,371) = strpad('K_H in component P_MAL (molar)');
    LEGEND_CONSTANTS(:,372) = strpad('K_Mg in component P_MAL (molar)');
    LEGEND_CONSTANTS(:,373) = strpad('K_K in component P_MAL (molar)');
    LEGEND_ALGEBRAIC(:,46) = strpad('P_CO2_tot_x in component P_CO2_tot (dimensionless)');
    LEGEND_ALGEBRAIC(:,137) = strpad('P_CO2_tot_i in component P_CO2_tot (dimensionless)');
    LEGEND_ALGEBRAIC(:,11) = strpad('P_CO2_tot_c in component P_CO2_tot (dimensionless)');
    LEGEND_CONSTANTS(:,374) = strpad('K_H in component P_CO2_tot (molar)');
    LEGEND_CONSTANTS(:,97) = strpad('K_Mg in component P_CO2_tot (molar)');
    LEGEND_CONSTANTS(:,98) = strpad('K_K in component P_CO2_tot (molar)');
    LEGEND_ALGEBRAIC(:,47) = strpad('P_PYR_x in component P_PYR (dimensionless)');
    LEGEND_ALGEBRAIC(:,153) = strpad('P_PYR_i in component P_PYR (dimensionless)');
    LEGEND_ALGEBRAIC(:,12) = strpad('P_PYR_c in component P_PYR (dimensionless)');
    LEGEND_CONSTANTS(:,375) = strpad('K_Mg in component P_PYR (molar)');
    LEGEND_CONSTANTS(:,99) = strpad('K_H in component P_PYR (molar)');
    LEGEND_CONSTANTS(:,100) = strpad('K_K in component P_PYR (molar)');
    LEGEND_ALGEBRAIC(:,13) = strpad('P_G6P_x in component P_G6P (dimensionless)');
    LEGEND_ALGEBRAIC(:,138) = strpad('P_G6P_i in component P_G6P (dimensionless)');
    LEGEND_ALGEBRAIC(:,48) = strpad('P_G6P_c in component P_G6P (dimensionless)');
    LEGEND_CONSTANTS(:,376) = strpad('K_H in component P_G6P (molar)');
    LEGEND_CONSTANTS(:,101) = strpad('K_Mg in component P_G6P (molar)');
    LEGEND_CONSTANTS(:,102) = strpad('K_K in component P_G6P (molar)');
    LEGEND_ALGEBRAIC(:,49) = strpad('P_GLU_x in component P_GLU (dimensionless)');
    LEGEND_ALGEBRAIC(:,154) = strpad('P_GLU_i in component P_GLU (dimensionless)');
    LEGEND_ALGEBRAIC(:,14) = strpad('P_GLU_c in component P_GLU (dimensionless)');
    LEGEND_CONSTANTS(:,377) = strpad('K_H in component P_GLU (molar)');
    LEGEND_CONSTANTS(:,378) = strpad('K_Mg in component P_GLU (molar)');
    LEGEND_CONSTANTS(:,103) = strpad('K_K in component P_GLU (molar)');
    LEGEND_ALGEBRAIC(:,50) = strpad('P_ASP_x in component P_ASP (dimensionless)');
    LEGEND_ALGEBRAIC(:,155) = strpad('P_ASP_i in component P_ASP (dimensionless)');
    LEGEND_ALGEBRAIC(:,15) = strpad('P_ASP_c in component P_ASP (dimensionless)');
    LEGEND_CONSTANTS(:,379) = strpad('K_H in component P_ASP (molar)');
    LEGEND_CONSTANTS(:,380) = strpad('K_Mg in component P_ASP (molar)');
    LEGEND_CONSTANTS(:,104) = strpad('K_K in component P_ASP (molar)');
    LEGEND_CONSTANTS(:,276) = strpad('P_GLC_c in component P_GLC (dimensionless)');
    LEGEND_CONSTANTS(:,105) = strpad('K_H in component P_GLC (molar)');
    LEGEND_CONSTANTS(:,106) = strpad('K_Mg in component P_GLC (molar)');
    LEGEND_CONSTANTS(:,107) = strpad('K_K in component P_GLC (molar)');
    LEGEND_CONSTANTS(:,277) = strpad('P_NAD_x in component P_NAD (dimensionless)');
    LEGEND_CONSTANTS(:,108) = strpad('K_H in component P_NAD (molar)');
    LEGEND_CONSTANTS(:,109) = strpad('K_Mg in component P_NAD (molar)');
    LEGEND_CONSTANTS(:,110) = strpad('K_K in component P_NAD (molar)');
    LEGEND_CONSTANTS(:,278) = strpad('P_NADH_x in component P_NADH (dimensionless)');
    LEGEND_CONSTANTS(:,111) = strpad('K_H in component P_NADH (molar)');
    LEGEND_CONSTANTS(:,112) = strpad('K_Mg in component P_NADH (molar)');
    LEGEND_CONSTANTS(:,113) = strpad('K_K in component P_NADH (molar)');
    LEGEND_CONSTANTS(:,305) = strpad('P_QH2_x in component P_QH2 (dimensionless)');
    LEGEND_CONSTANTS(:,114) = strpad('K_H in component P_QH2 (molar)');
    LEGEND_CONSTANTS(:,115) = strpad('K_Mg in component P_QH2 (molar)');
    LEGEND_CONSTANTS(:,116) = strpad('K_K in component P_QH2 (molar)');
    LEGEND_CONSTANTS(:,321) = strpad('P_COQ_x in component P_COQ (dimensionless)');
    LEGEND_CONSTANTS(:,117) = strpad('K_H in component P_COQ (molar)');
    LEGEND_CONSTANTS(:,118) = strpad('K_Mg in component P_COQ (molar)');
    LEGEND_CONSTANTS(:,119) = strpad('K_K in component P_COQ (molar)');
    LEGEND_CONSTANTS(:,279) = strpad('P_ACCOA_x in component P_ACCOA (dimensionless)');
    LEGEND_CONSTANTS(:,120) = strpad('K_H in component P_ACCOA (molar)');
    LEGEND_CONSTANTS(:,121) = strpad('K_Mg in component P_ACCOA (molar)');
    LEGEND_CONSTANTS(:,122) = strpad('K_K in component P_ACCOA (molar)');
    LEGEND_CONSTANTS(:,307) = strpad('P_AKG_x in component P_AKG (dimensionless)');
    LEGEND_CONSTANTS(:,123) = strpad('K_H in component P_AKG (molar)');
    LEGEND_CONSTANTS(:,124) = strpad('K_Mg in component P_AKG (molar)');
    LEGEND_CONSTANTS(:,125) = strpad('K_K in component P_AKG (molar)');
    LEGEND_CONSTANTS(:,280) = strpad('P_O2_x in component P_O2 (dimensionless)');
    LEGEND_CONSTANTS(:,126) = strpad('K_H in component P_O2 (molar)');
    LEGEND_CONSTANTS(:,127) = strpad('K_Mg in component P_O2 (molar)');
    LEGEND_CONSTANTS(:,128) = strpad('K_K in component P_O2 (molar)');
    LEGEND_CONSTANTS(:,281) = strpad('P_FADH2_x in component P_FADH2 (dimensionless)');
    LEGEND_CONSTANTS(:,129) = strpad('K_H in component P_FADH2 (molar)');
    LEGEND_CONSTANTS(:,130) = strpad('K_Mg in component P_FADH2 (molar)');
    LEGEND_CONSTANTS(:,131) = strpad('K_K in component P_FADH2 (molar)');
    LEGEND_CONSTANTS(:,282) = strpad('P_FAD_x in component P_FAD (dimensionless)');
    LEGEND_CONSTANTS(:,132) = strpad('K_H in component P_FAD (molar)');
    LEGEND_CONSTANTS(:,133) = strpad('K_Mg in component P_FAD (molar)');
    LEGEND_CONSTANTS(:,134) = strpad('K_K in component P_FAD (molar)');
    LEGEND_CONSTANTS(:,283) = strpad('P_H_x in component P_H (dimensionless)');
    LEGEND_CONSTANTS(:,135) = strpad('K_H in component P_H (molar)');
    LEGEND_CONSTANTS(:,136) = strpad('K_Mg in component P_H (molar)');
    LEGEND_CONSTANTS(:,137) = strpad('K_K in component P_H (molar)');
    LEGEND_CONSTANTS(:,284) = strpad('P_Mg_x in component P_Mg (dimensionless)');
    LEGEND_CONSTANTS(:,138) = strpad('K_H in component P_Mg (molar)');
    LEGEND_CONSTANTS(:,139) = strpad('K_Mg in component P_Mg (molar)');
    LEGEND_CONSTANTS(:,140) = strpad('K_K in component P_Mg (molar)');
    LEGEND_CONSTANTS(:,285) = strpad('P_K_x in component P_K (dimensionless)');
    LEGEND_CONSTANTS(:,141) = strpad('K_H in component P_K (molar)');
    LEGEND_CONSTANTS(:,142) = strpad('K_Mg in component P_K (molar)');
    LEGEND_CONSTANTS(:,143) = strpad('K_K in component P_K (molar)');
    LEGEND_ALGEBRAIC(:,115) = strpad('J_pdh in component J_pdh (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,18) = strpad('PYR_x in component dPYR_x_dt (molar)');
    LEGEND_STATES(:,19) = strpad('COASH_x in component dCOASH_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,114) = strpad('NAD_x in component dNAD_x_dt (molar)');
    LEGEND_STATES(:,20) = strpad('CO2_tot_x in component dCO2_tot_x_dt (molar)');
    LEGEND_STATES(:,21) = strpad('ACCOA_x in component dACCOA_x_dt (molar)');
    LEGEND_STATES(:,22) = strpad('NADH_x in component dNADH_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,51) = strpad('K_eq_pdh in component J_pdh (dimensionless)');
    LEGEND_CONSTANTS(:,315) = strpad('K0_eq_pdh in component J_pdh (molar)');
    LEGEND_CONSTANTS(:,299) = strpad('dGr0_pdh in component J_pdh (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,144) = strpad('K_mA in component J_pdh (molar)');
    LEGEND_CONSTANTS(:,145) = strpad('K_mB in component J_pdh (molar)');
    LEGEND_CONSTANTS(:,146) = strpad('K_mC in component J_pdh (molar)');
    LEGEND_ALGEBRAIC(:,52) = strpad('alpha_i1 in component J_pdh (dimensionless)');
    LEGEND_ALGEBRAIC(:,53) = strpad('alpha_i2 in component J_pdh (dimensionless)');
    LEGEND_CONSTANTS(:,147) = strpad('K_iACCOA in component J_pdh (molar)');
    LEGEND_CONSTANTS(:,148) = strpad('K_iNADH in component J_pdh (molar)');
    LEGEND_CONSTANTS(:,149) = strpad('minCon in component J_pdh (molar)');
    LEGEND_ALGEBRAIC(:,57) = strpad('J_cits in component J_cits (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,23) = strpad('OAA_x in component dOAA_x_dt (molar)');
    LEGEND_STATES(:,24) = strpad('SCOA_x in component dSCOA_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,54) = strpad('K_eq_cits in component J_cits (dimensionless)');
    LEGEND_CONSTANTS(:,316) = strpad('K0_eq_cits in component J_cits (molar_squared)');
    LEGEND_CONSTANTS(:,300) = strpad('dGr0_cits in component J_cits (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,150) = strpad('K_mA in component J_cits (molar)');
    LEGEND_CONSTANTS(:,151) = strpad('K_mB in component J_cits (molar)');
    LEGEND_CONSTANTS(:,152) = strpad('K_ia in component J_cits (molar)');
    LEGEND_CONSTANTS(:,153) = strpad('K_iCIT in component J_cits (molar)');
    LEGEND_CONSTANTS(:,154) = strpad('K_iATP in component J_cits (molar)');
    LEGEND_CONSTANTS(:,155) = strpad('K_iADP in component J_cits (molar)');
    LEGEND_CONSTANTS(:,156) = strpad('K_iAMP in component J_cits (molar)');
    LEGEND_CONSTANTS(:,157) = strpad('K_iCOASH in component J_cits (molar)');
    LEGEND_CONSTANTS(:,158) = strpad('K_iSCOA in component J_cits (molar)');
    LEGEND_ALGEBRAIC(:,55) = strpad('alpha_i1 in component J_cits (dimensionless)');
    LEGEND_ALGEBRAIC(:,56) = strpad('alpha_i2 in component J_cits (dimensionless)');
    LEGEND_ALGEBRAIC(:,60) = strpad('J_acon in component J_acon (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,25) = strpad('ICIT_x in component dICIT_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,59) = strpad('V_mr in component J_acon (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,58) = strpad('K_eq_acon in component J_acon (dimensionless)');
    LEGEND_CONSTANTS(:,317) = strpad('K0_eq_acon in component J_acon (dimensionless)');
    LEGEND_CONSTANTS(:,301) = strpad('dGr0_acon in component J_acon (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,159) = strpad('K_mA in component J_acon (molar)');
    LEGEND_CONSTANTS(:,160) = strpad('K_mP in component J_acon (molar)');
    LEGEND_ALGEBRAIC(:,116) = strpad('J_isod in component J_isod (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,26) = strpad('AKG_x in component dAKG_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,61) = strpad('K_eq_isod in component J_isod (molar)');
    LEGEND_CONSTANTS(:,318) = strpad('K0_eq_isod in component J_isod (molar_cubed)');
    LEGEND_CONSTANTS(:,302) = strpad('dGr0_isod in component J_isod (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,161) = strpad('K_mA in component J_isod (molar)');
    LEGEND_CONSTANTS(:,162) = strpad('K_mB in component J_isod (molar)');
    LEGEND_CONSTANTS(:,163) = strpad('K_ib in component J_isod (molar)');
    LEGEND_CONSTANTS(:,164) = strpad('K_iq in component J_isod (molar)');
    LEGEND_ALGEBRAIC(:,62) = strpad('alpha_i in component J_isod (dimensionless)');
    LEGEND_CONSTANTS(:,165) = strpad('n_H in component J_isod (dimensionless)');
    LEGEND_CONSTANTS(:,166) = strpad('K_iATP in component J_isod (molar)');
    LEGEND_CONSTANTS(:,167) = strpad('K_aADP in component J_isod (molar)');
    LEGEND_CONSTANTS(:,168) = strpad('minCon in component J_isod (molar)');
    LEGEND_ALGEBRAIC(:,117) = strpad('J_akgd in component J_akgd (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,63) = strpad('K_eq_akgd in component J_akgd (dimensionless)');
    LEGEND_CONSTANTS(:,319) = strpad('K0_eq_akgd in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,303) = strpad('dGr0_akgd in component J_akgd (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,169) = strpad('K_mA in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,170) = strpad('K_mB in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,171) = strpad('K_mC in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,172) = strpad('K_iq in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,173) = strpad('K_ir1 in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,174) = strpad('K_ir2 in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,175) = strpad('K_iATP in component J_akgd (molar)');
    LEGEND_CONSTANTS(:,176) = strpad('K_aADP in component J_akgd (molar)');
    LEGEND_ALGEBRAIC(:,64) = strpad('alpha_i in component J_akgd (dimensionless)');
    LEGEND_CONSTANTS(:,177) = strpad('minCon in component J_akgd (molar)');
    LEGEND_ALGEBRAIC(:,67) = strpad('J_scoas in component J_scoas (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,27) = strpad('SUC_x in component dSUC_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,65) = strpad('K_eq_scoas in component J_scoas (dimensionless)');
    LEGEND_CONSTANTS(:,320) = strpad('K0_eq_scoas in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,304) = strpad('dGr0_scoas in component J_scoas (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,66) = strpad('V_mr in component J_scoas (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,178) = strpad('K_ia in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,179) = strpad('K_ib in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,180) = strpad('K_ic in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,181) = strpad('K_ip in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,182) = strpad('K_iq in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,183) = strpad('K_ir in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,184) = strpad('K_mA in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,185) = strpad('K_mB in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,186) = strpad('K_mC in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,187) = strpad('K_mP in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,188) = strpad('K_mQ in component J_scoas (molar)');
    LEGEND_CONSTANTS(:,189) = strpad('K_mR in component J_scoas (molar)');
    LEGEND_ALGEBRAIC(:,120) = strpad('J_sdh in component J_sdh (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,119) = strpad('COQ_x in component dCOQ_x_dt (molar)');
    LEGEND_STATES(:,28) = strpad('QH2_x in component dQH2_x_dt (molar)');
    LEGEND_STATES(:,29) = strpad('FUM_x in component dFUM_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,68) = strpad('K_eq_sdh in component J_sdh (dimensionless)');
    LEGEND_CONSTANTS(:,335) = strpad('K0_eq_sdh in component J_sdh (dimensionless)');
    LEGEND_CONSTANTS(:,328) = strpad('dGr0_sdh in component J_sdh (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,69) = strpad('V_mr in component J_sdh (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,190) = strpad('K_ia in component J_sdh (molar)');
    LEGEND_CONSTANTS(:,191) = strpad('K_iq in component J_sdh (molar)');
    LEGEND_CONSTANTS(:,192) = strpad('K_mSUC in component J_sdh (molar)');
    LEGEND_CONSTANTS(:,193) = strpad('K_mCOQ in component J_sdh (molar)');
    LEGEND_CONSTANTS(:,194) = strpad('K_mQH2 in component J_sdh (molar)');
    LEGEND_CONSTANTS(:,195) = strpad('K_mFUM in component J_sdh (molar)');
    LEGEND_ALGEBRAIC(:,70) = strpad('alpha_i in component J_sdh (dimensionless)');
    LEGEND_CONSTANTS(:,196) = strpad('K_iOAA in component J_sdh (molar)');
    LEGEND_CONSTANTS(:,197) = strpad('K_aSUC in component J_sdh (molar)');
    LEGEND_CONSTANTS(:,198) = strpad('K_aFUM in component J_sdh (molar)');
    LEGEND_ALGEBRAIC(:,74) = strpad('J_fum in component J_fum (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,30) = strpad('MAL_x in component dMAL_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,71) = strpad('K_eq_fum in component J_fum (dimensionless)');
    LEGEND_CONSTANTS(:,343) = strpad('K0_eq_fum in component J_fum (dimensionless)');
    LEGEND_CONSTANTS(:,339) = strpad('dGr0_fum in component J_fum (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,72) = strpad('V_mr in component J_fum (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,199) = strpad('K_mFUM in component J_fum (molar)');
    LEGEND_CONSTANTS(:,200) = strpad('K_mMAL in component J_fum (molar)');
    LEGEND_ALGEBRAIC(:,73) = strpad('alpha_i in component J_fum (dimensionless)');
    LEGEND_CONSTANTS(:,201) = strpad('K_iCIT in component J_fum (molar)');
    LEGEND_CONSTANTS(:,202) = strpad('K_iATP in component J_fum (molar)');
    LEGEND_CONSTANTS(:,203) = strpad('K_iADP in component J_fum (molar)');
    LEGEND_CONSTANTS(:,204) = strpad('K_iGTP in component J_fum (molar)');
    LEGEND_CONSTANTS(:,205) = strpad('K_iGDP in component J_fum (molar)');
    LEGEND_ALGEBRAIC(:,118) = strpad('J_mdh in component J_mdh (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,75) = strpad('K_eq_mdh in component J_mdh (dimensionless)');
    LEGEND_CONSTANTS(:,322) = strpad('K0_eq_mdh in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,306) = strpad('dGr0_mdh in component J_mdh (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,76) = strpad('V_mr in component J_mdh (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,206) = strpad('K_ia in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,207) = strpad('K_ib in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,208) = strpad('K_ip in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,209) = strpad('K_iq in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,210) = strpad('K_mNAD in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,211) = strpad('K_mMAL in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,212) = strpad('K_mOAA in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,213) = strpad('K_mNADH in component J_mdh (molar)');
    LEGEND_ALGEBRAIC(:,77) = strpad('alpha_i in component J_mdh (dimensionless)');
    LEGEND_CONSTANTS(:,214) = strpad('K_iATP in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,215) = strpad('K_iADP in component J_mdh (molar)');
    LEGEND_CONSTANTS(:,216) = strpad('K_iAMP in component J_mdh (molar)');
    LEGEND_ALGEBRAIC(:,81) = strpad('J_ndk_f in component J_ndk (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,82) = strpad('J_ndk in component J_ndk (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,78) = strpad('K_eq_ndk in component J_ndk (dimensionless)');
    LEGEND_CONSTANTS(:,313) = strpad('K0_eq_ndk in component J_ndk (dimensionless)');
    LEGEND_CONSTANTS(:,298) = strpad('dGr0_ndk in component J_ndk (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,79) = strpad('V_mr in component J_ndk (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,217) = strpad('K_ia in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,218) = strpad('K_ib in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,219) = strpad('K_ip in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,220) = strpad('K_iq in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,221) = strpad('K_mA in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,222) = strpad('K_mB in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,223) = strpad('K_mP in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,224) = strpad('K_mQ in component J_ndk (molar)');
    LEGEND_ALGEBRAIC(:,80) = strpad('alpha_i in component J_ndk (dimensionless)');
    LEGEND_CONSTANTS(:,225) = strpad('K_iAMP in component J_ndk (molar)');
    LEGEND_CONSTANTS(:,226) = strpad('minCon in component J_ndk (molar)');
    LEGEND_ALGEBRAIC(:,86) = strpad('J_got_f in component J_got (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,87) = strpad('J_got in component J_got (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,31) = strpad('ASP_x in component dASP_x_dt (molar)');
    LEGEND_STATES(:,32) = strpad('GLU_x in component dGLU_x_dt (molar)');
    LEGEND_ALGEBRAIC(:,83) = strpad('K_eq_got in component J_got (dimensionless)');
    LEGEND_CONSTANTS(:,329) = strpad('K0_eq_got in component J_got (dimensionless)');
    LEGEND_CONSTANTS(:,323) = strpad('dGr0_got in component J_got (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,84) = strpad('V_mr in component J_got (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,227) = strpad('K_ia in component J_got (molar)');
    LEGEND_CONSTANTS(:,228) = strpad('K_ib in component J_got (molar)');
    LEGEND_CONSTANTS(:,229) = strpad('K_ip in component J_got (molar)');
    LEGEND_CONSTANTS(:,230) = strpad('K_iq in component J_got (molar)');
    LEGEND_CONSTANTS(:,231) = strpad('K_mA in component J_got (molar)');
    LEGEND_CONSTANTS(:,232) = strpad('K_mB in component J_got (molar)');
    LEGEND_CONSTANTS(:,233) = strpad('K_mP in component J_got (molar)');
    LEGEND_CONSTANTS(:,234) = strpad('K_mQ in component J_got (molar)');
    LEGEND_ALGEBRAIC(:,85) = strpad('alpha_i in component J_got (dimensionless)');
    LEGEND_CONSTANTS(:,235) = strpad('K_iAKG in component J_got (molar)');
    LEGEND_CONSTANTS(:,236) = strpad('minCon in component J_got (molar)');
    LEGEND_ALGEBRAIC(:,156) = strpad('J_PYRH in component J_PYRH (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,33) = strpad('PYR_i in component dPYR_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,157) = strpad('J_GLUH in component J_GLUH (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,34) = strpad('GLU_i in component dGLU_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,176) = strpad('J_CITMAL in component J_CITMAL (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,35) = strpad('MAL_i in component dMAL_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,88) = strpad('J_AKGMAL in component J_AKGMAL (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,237) = strpad('K_mAKGi in component J_AKGMAL (molar)');
    LEGEND_CONSTANTS(:,238) = strpad('K_mAKGx in component J_AKGMAL (molar)');
    LEGEND_CONSTANTS(:,239) = strpad('K_mMALi in component J_AKGMAL (molar)');
    LEGEND_CONSTANTS(:,240) = strpad('K_mMALx in component J_AKGMAL (molar)');
    LEGEND_STATES(:,36) = strpad('AKG_i in component dAKG_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,177) = strpad('J_MALPI in component J_MALPI (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,159) = strpad('J_ASPGLU in component J_ASPGLU (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,158) = strpad('K_eqASPGLU in component J_ASPGLU (dimensionless)');
    LEGEND_CONSTANTS(:,286) = strpad('K_hASPGLU in component J_ASPGLU (molar)');
    LEGEND_CONSTANTS(:,241) = strpad('K_iASPi in component J_ASPGLU (molar)');
    LEGEND_CONSTANTS(:,242) = strpad('K_iASPx in component J_ASPGLU (molar)');
    LEGEND_CONSTANTS(:,243) = strpad('K_iGLUi in component J_ASPGLU (molar)');
    LEGEND_CONSTANTS(:,244) = strpad('K_iGLUx in component J_ASPGLU (molar)');
    LEGEND_STATES(:,37) = strpad('ASP_i in component dASP_i_dt (molar)');
    LEGEND_CONSTANTS(:,245) = strpad('m in component J_ASPGLU (dimensionless)');
    LEGEND_CONSTANTS(:,246) = strpad('pKa_ASPGLU in component J_ASPGLU (dimensionless)');
    LEGEND_STATES(:,38) = strpad('dPsi in component ddPsi_dt (millivolt)');
    LEGEND_ALGEBRAIC(:,178) = strpad('J_SUCMAL in component J_SUCMAL (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,39) = strpad('SUC_i in component dSUC_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,89) = strpad('J_PYRt in component J_PYRt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,40) = strpad('PYR_c in component dPYR_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,91) = strpad('J_CITt in component J_CITt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,41) = strpad('CIT_c in component dCIT_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,92) = strpad('J_MALt in component J_MALt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,42) = strpad('MAL_c in component dMAL_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,90) = strpad('J_AKGt in component J_AKGt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,43) = strpad('AKG_c in component dAKG_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,93) = strpad('J_SUCt in component J_SUCt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,44) = strpad('SUC_c in component dSUC_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,94) = strpad('J_GLUt in component J_GLUt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,45) = strpad('GLU_c in component dGLU_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,95) = strpad('J_ASPt in component J_ASPt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,46) = strpad('ASP_c in component dASP_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,16) = strpad('J_FUMt in component J_FUMt (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,247) = strpad('zero in component J_FUMt (micrometer_per_second)');
    LEGEND_STATES(:,47) = strpad('FUM_c in component dFUM_c_dt (molar)');
    LEGEND_STATES(:,48) = strpad('FUM_i in component dFUM_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,17) = strpad('J_ICITt in component J_ICITt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,49) = strpad('ICIT_c in component dICIT_c_dt (molar)');
    LEGEND_STATES(:,50) = strpad('ICIT_i in component dICIT_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,96) = strpad('J_ADPt in component J_ADPt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,51) = strpad('ADP_i in component dADP_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,97) = strpad('J_ATPt in component J_ATPt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,52) = strpad('ATP_i in component dATP_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,98) = strpad('J_AMPt in component J_AMPt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,53) = strpad('AMP_c in component dAMP_c_dt (molar)');
    LEGEND_STATES(:,54) = strpad('AMP_i in component dAMP_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,99) = strpad('J_PIt in component J_PIt (mole_per_second_per_l_mito_volume)');
    LEGEND_STATES(:,55) = strpad('PI_c in component dPI_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,140) = strpad('J_C1 in component J_C1 (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,288) = strpad('dGr0_C1 in component J_C1 (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,101) = strpad('K_eq_C1 in component J_C1 (per_molar)');
    LEGEND_ALGEBRAIC(:,139) = strpad('K_app_C1 in component J_C1 (dimensionless)');
    LEGEND_ALGEBRAIC(:,123) = strpad('J_C3 in component J_C3 (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,289) = strpad('dGr0_C3 in component J_C3 (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,103) = strpad('K_eq_C3 in component J_C3 (molar_squared)');
    LEGEND_ALGEBRAIC(:,105) = strpad('K_app_C3 in component J_C3 (dimensionless)');
    LEGEND_ALGEBRAIC(:,107) = strpad('QH2_x in component J_C3 (molar)');
    LEGEND_ALGEBRAIC(:,121) = strpad('COQ_x in component J_C3 (molar)');
    LEGEND_STATES(:,56) = strpad('Cred_i in component dCred_i_dt (molar)');
    LEGEND_ALGEBRAIC(:,122) = strpad('Cox_i in component dCox_i_dt (molar)');
    LEGEND_CONSTANTS(:,248) = strpad('minCon in component J_C3 (molar)');
    LEGEND_ALGEBRAIC(:,142) = strpad('J_C4 in component J_C4 (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,290) = strpad('dGr0_C4 in component J_C4 (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,109) = strpad('K_eq_C4 in component J_C4 (dimensionless)');
    LEGEND_ALGEBRAIC(:,141) = strpad('K_app_C4 in component J_C4 (dimensionless)');
    LEGEND_ALGEBRAIC(:,111) = strpad('O2_x in component J_C4 (molar)');
    LEGEND_STATES(:,57) = strpad('O2_x in component dO2_x_dt (molar)');
    LEGEND_CONSTANTS(:,249) = strpad('minCon in component J_C4 (molar)');
    LEGEND_ALGEBRAIC(:,144) = strpad('J_F1 in component J_F1 (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,291) = strpad('dGr0_F1 in component J_F1 (kilojoule_per_mole)');
    LEGEND_ALGEBRAIC(:,113) = strpad('K_eq_F1 in component J_F1 (per_molar_squared)');
    LEGEND_ALGEBRAIC(:,143) = strpad('K_app_F1 in component J_F1 (per_molar)');
    LEGEND_ALGEBRAIC(:,179) = strpad('J_ANT in component J_ANT (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,250) = strpad('minCond in component J_ANT (molar)');
    LEGEND_ALGEBRAIC(:,180) = strpad('J_PIHt in component J_PIHt (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,145) = strpad('J_Hle in component J_Hle (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,251) = strpad('minCond in component J_Hle (millivolt)');
    LEGEND_ALGEBRAIC(:,182) = strpad('J_KH in component J_KH (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,104) = strpad('J_HK in component J_HK (mole_per_second_per_l_cyto_volume)');
    LEGEND_CONSTANTS(:,309) = strpad('dGr0_HK in component J_HK (kilojoule_per_mole)');
    LEGEND_STATES(:,58) = strpad('GLC_c in component dGLC_c_dt (molar)');
    LEGEND_STATES(:,59) = strpad('G6P_c in component dG6P_c_dt (molar)');
    LEGEND_CONSTANTS(:,252) = strpad('K_ia in component J_HK (molar)');
    LEGEND_CONSTANTS(:,253) = strpad('K_ib in component J_HK (molar)');
    LEGEND_CONSTANTS(:,254) = strpad('K_ip in component J_HK (molar)');
    LEGEND_CONSTANTS(:,255) = strpad('K_iq in component J_HK (molar)');
    LEGEND_CONSTANTS(:,256) = strpad('K_mA in component J_HK (molar)');
    LEGEND_CONSTANTS(:,257) = strpad('K_mB in component J_HK (molar)');
    LEGEND_CONSTANTS(:,258) = strpad('K_mP in component J_HK (molar)');
    LEGEND_CONSTANTS(:,259) = strpad('K_mQ in component J_HK (molar)');
    LEGEND_CONSTANTS(:,260) = strpad('Ki_G6P in component J_HK (molar)');
    LEGEND_CONSTANTS(:,325) = strpad('K_eq_HK in component J_HK (molar)');
    LEGEND_ALGEBRAIC(:,100) = strpad('Kapp_HK in component J_HK (dimensionless)');
    LEGEND_ALGEBRAIC(:,102) = strpad('Kapp_HK_m in component J_HK (dimensionless)');
    LEGEND_ALGEBRAIC(:,108) = strpad('J_CKc in component J_CKc (mole_per_second_per_l_cyto_volume)');
    LEGEND_CONSTANTS(:,261) = strpad('X_CK in component J_CKc (mole_per_second_per_l_cyto_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,310) = strpad('K_CK in component J_CKc (per_molar)');
    LEGEND_STATES(:,60) = strpad('PCr_c in component dPCr_c_dt (molar)');
    LEGEND_ALGEBRAIC(:,106) = strpad('Cr_c in component J_CKc (molar)');
    LEGEND_ALGEBRAIC(:,183) = strpad('J_AKi in component J_AK (mole_per_second_per_l_mito_volume)');
    LEGEND_ALGEBRAIC(:,112) = strpad('J_AKc in component J_AK (mole_per_second_per_l_mito_volume)');
    LEGEND_CONSTANTS(:,262) = strpad('X_AK in component J_AK (mole_per_second_per_l_mito_volume_per_molar_per_molar)');
    LEGEND_CONSTANTS(:,381) = strpad('dGr0_AK in component J_AK (kilojoule_per_mole)');
    LEGEND_CONSTANTS(:,382) = strpad('K_eq_AK in component J_AK (dimensionless)');
    LEGEND_ALGEBRAIC(:,181) = strpad('Kapp_AKi in component J_AK (dimensionless)');
    LEGEND_ALGEBRAIC(:,110) = strpad('Kapp_AKc in component J_AK (dimensionless)');
    LEGEND_ALGEBRAIC(:,147) = strpad('dNAD_x_dt in component dNAD_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,146) = strpad('dNADH_x_dt in component dNADH_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,149) = strpad('dCOQ_x_dt in component dCOQ_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,148) = strpad('dQH2_x_dt in component dQH2_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,184) = strpad('dATP_x_dt in component dATP_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,185) = strpad('dADP_x_dt in component dADP_x_dt (molar_per_second)');
    LEGEND_CONSTANTS(:,263) = strpad('dAMP_x_dt in component dAMP_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,124) = strpad('dGTP_x_dt in component dGTP_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,125) = strpad('dGDP_x_dt in component dGDP_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,186) = strpad('dPI_x_dt in component dPI_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,160) = strpad('dPYR_x_dt in component dPYR_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,126) = strpad('dACCOA_x_dt in component dACCOA_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,187) = strpad('dCIT_x_dt in component dCIT_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,127) = strpad('dICIT_x_dt in component dICIT_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,128) = strpad('dAKG_x_dt in component dAKG_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,129) = strpad('dSCOA_x_dt in component dSCOA_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,130) = strpad('dCOASH_x_dt in component dCOASH_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,188) = strpad('dSUC_x_dt in component dSUC_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,131) = strpad('dFUM_x_dt in component dFUM_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,189) = strpad('dMAL_x_dt in component dMAL_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,132) = strpad('dOAA_x_dt in component dOAA_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,161) = strpad('dGLU_x_dt in component dGLU_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,162) = strpad('dASP_x_dt in component dASP_x_dt (molar_per_second)');
    LEGEND_CONSTANTS(:,264) = strpad('dO2_x_dt in component dO2_x_dt (molar_per_second)');
    LEGEND_CONSTANTS(:,265) = strpad('dCO2_tot_x_dt in component dCO2_tot_x_dt (molar_per_second)');
    LEGEND_CONSTANTS(:,332) = strpad('Rm_cyto in component dPYR_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,336) = strpad('Rm_cyto in component dCIT_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,331) = strpad('Rm_cyto in component dICIT_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,340) = strpad('Rm_cyto in component dAKG_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,344) = strpad('Rm_cyto in component dSUC_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,330) = strpad('Rm_cyto in component dFUM_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,347) = strpad('Rm_cyto in component dMAL_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,349) = strpad('Rm_cyto in component dGLU_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,351) = strpad('Rm_cyto in component dASP_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,353) = strpad('Rm_cyto in component dPI_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,266) = strpad('conversion_factor in component dPI_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,355) = strpad('Rm_cyto in component dATP_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,267) = strpad('conversion_factor in component dATP_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,357) = strpad('Rm_cyto in component dADP_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,268) = strpad('conversion_factor in component dADP_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,359) = strpad('Rm_cyto in component dAMP_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,269) = strpad('conversion_factor in component dAMP_c_dt (mito_per_cyto)');
    LEGEND_CONSTANTS(:,311) = strpad('FAD_x in component dFAD_x_dt (molar)');
    LEGEND_CONSTANTS(:,270) = strpad('dFAD_x_dt in component dFAD_x_dt (molar_per_second)');
    LEGEND_CONSTANTS(:,271) = strpad('dFADH2_x_dt in component dFAD_x_dt (molar_per_second)');
    LEGEND_CONSTANTS(:,292) = strpad('FADH2_x in component dFAD_x_dt (molar)');
    LEGEND_CONSTANTS(:,293) = strpad('dH_x_dt in component dH_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,207) = strpad('D in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,201) = strpad('phi_H in component BINDING_IONS (molar_per_second)');
    LEGEND_ALGEBRAIC(:,202) = strpad('phi_Mg in component BINDING_IONS (molar_per_second)');
    LEGEND_ALGEBRAIC(:,203) = strpad('phi_K in component BINDING_IONS (molar_per_second)');
    LEGEND_ALGEBRAIC(:,190) = strpad('dHB_dMg in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,191) = strpad('dHB_dK in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,194) = strpad('dMgB_dK in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,197) = strpad('dKB_dMg in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,206) = strpad('alpha_K in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,205) = strpad('alpha_Mg in component BINDING_IONS (dimensionless)');
    LEGEND_CONSTANTS(:,294) = strpad('dMg_x_dt in component dMg_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,204) = strpad('alpha_H in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,196) = strpad('dKB_dH in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,193) = strpad('dMgB_dH in component BINDING_IONS (dimensionless)');
    LEGEND_CONSTANTS(:,295) = strpad('dK_x_dt in component dK_x_dt (molar_per_second)');
    LEGEND_ALGEBRAIC(:,192) = strpad('dHB_dH in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,195) = strpad('dMgB_dMg in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,198) = strpad('dKB_dK in component BINDING_IONS (dimensionless)');
    LEGEND_ALGEBRAIC(:,199) = strpad('J_H_t in component BINDING_IONS (molar_per_second)');
    LEGEND_CONSTANTS(:,296) = strpad('J_Mg_t in component BINDING_IONS (molar_per_second)');
    LEGEND_ALGEBRAIC(:,200) = strpad('J_K_t in component BINDING_IONS (molar_per_second)');
    LEGEND_CONSTANTS(:,272) = strpad('B_x in component BINDING_IONS (molar)');
    LEGEND_CONSTANTS(:,273) = strpad('K_BX in component BINDING_IONS (molar)');
    LEGEND_RATES(:,38) = strpad('d/dt dPsi in component ddPsi_dt (millivolt)');
    LEGEND_RATES(:,1) = strpad('d/dt ATP_x in component dATP_x_dt (molar)');
    LEGEND_RATES(:,9) = strpad('d/dt ADP_x in component dADP_x_dt (molar)');
    LEGEND_RATES(:,11) = strpad('d/dt AMP_x in component dAMP_x_dt (molar)');
    LEGEND_RATES(:,12) = strpad('d/dt GTP_x in component dGTP_x_dt (molar)');
    LEGEND_RATES(:,13) = strpad('d/dt GDP_x in component dGDP_x_dt (molar)');
    LEGEND_RATES(:,14) = strpad('d/dt PI_x in component dPI_x_dt (molar)');
    LEGEND_RATES(:,22) = strpad('d/dt NADH_x in component dNADH_x_dt (molar)');
    LEGEND_RATES(:,28) = strpad('d/dt QH2_x in component dQH2_x_dt (molar)');
    LEGEND_RATES(:,18) = strpad('d/dt PYR_x in component dPYR_x_dt (molar)');
    LEGEND_RATES(:,21) = strpad('d/dt ACCOA_x in component dACCOA_x_dt (molar)');
    LEGEND_RATES(:,16) = strpad('d/dt CIT_x in component dCIT_x_dt (molar)');
    LEGEND_RATES(:,25) = strpad('d/dt ICIT_x in component dICIT_x_dt (molar)');
    LEGEND_RATES(:,26) = strpad('d/dt AKG_x in component dAKG_x_dt (molar)');
    LEGEND_RATES(:,24) = strpad('d/dt SCOA_x in component dSCOA_x_dt (molar)');
    LEGEND_RATES(:,19) = strpad('d/dt COASH_x in component dCOASH_x_dt (molar)');
    LEGEND_RATES(:,27) = strpad('d/dt SUC_x in component dSUC_x_dt (molar)');
    LEGEND_RATES(:,29) = strpad('d/dt FUM_x in component dFUM_x_dt (molar)');
    LEGEND_RATES(:,30) = strpad('d/dt MAL_x in component dMAL_x_dt (molar)');
    LEGEND_RATES(:,23) = strpad('d/dt OAA_x in component dOAA_x_dt (molar)');
    LEGEND_RATES(:,32) = strpad('d/dt GLU_x in component dGLU_x_dt (molar)');
    LEGEND_RATES(:,31) = strpad('d/dt ASP_x in component dASP_x_dt (molar)');
    LEGEND_RATES(:,57) = strpad('d/dt O2_x in component dO2_x_dt (molar)');
    LEGEND_RATES(:,20) = strpad('d/dt CO2_tot_x in component dCO2_tot_x_dt (molar)');
    LEGEND_RATES(:,56) = strpad('d/dt Cred_i in component dCred_i_dt (molar)');
    LEGEND_RATES(:,52) = strpad('d/dt ATP_i in component dATP_i_dt (molar)');
    LEGEND_RATES(:,51) = strpad('d/dt ADP_i in component dADP_i_dt (molar)');
    LEGEND_RATES(:,54) = strpad('d/dt AMP_i in component dAMP_i_dt (molar)');
    LEGEND_RATES(:,15) = strpad('d/dt PI_i in component dPI_i_dt (molar)');
    LEGEND_RATES(:,33) = strpad('d/dt PYR_i in component dPYR_i_dt (molar)');
    LEGEND_RATES(:,17) = strpad('d/dt CIT_i in component dCIT_i_dt (molar)');
    LEGEND_RATES(:,50) = strpad('d/dt ICIT_i in component dICIT_i_dt (molar)');
    LEGEND_RATES(:,36) = strpad('d/dt AKG_i in component dAKG_i_dt (molar)');
    LEGEND_RATES(:,39) = strpad('d/dt SUC_i in component dSUC_i_dt (molar)');
    LEGEND_RATES(:,48) = strpad('d/dt FUM_i in component dFUM_i_dt (molar)');
    LEGEND_RATES(:,35) = strpad('d/dt MAL_i in component dMAL_i_dt (molar)');
    LEGEND_RATES(:,34) = strpad('d/dt GLU_i in component dGLU_i_dt (molar)');
    LEGEND_RATES(:,37) = strpad('d/dt ASP_i in component dASP_i_dt (molar)');
    LEGEND_RATES(:,40) = strpad('d/dt PYR_c in component dPYR_c_dt (molar)');
    LEGEND_RATES(:,41) = strpad('d/dt CIT_c in component dCIT_c_dt (molar)');
    LEGEND_RATES(:,49) = strpad('d/dt ICIT_c in component dICIT_c_dt (molar)');
    LEGEND_RATES(:,43) = strpad('d/dt AKG_c in component dAKG_c_dt (molar)');
    LEGEND_RATES(:,44) = strpad('d/dt SUC_c in component dSUC_c_dt (molar)');
    LEGEND_RATES(:,47) = strpad('d/dt FUM_c in component dFUM_c_dt (molar)');
    LEGEND_RATES(:,42) = strpad('d/dt MAL_c in component dMAL_c_dt (molar)');
    LEGEND_RATES(:,45) = strpad('d/dt GLU_c in component dGLU_c_dt (molar)');
    LEGEND_RATES(:,46) = strpad('d/dt ASP_c in component dASP_c_dt (molar)');
    LEGEND_RATES(:,6) = strpad('d/dt H_c in component dH_c_dt (molar)');
    LEGEND_RATES(:,7) = strpad('d/dt Mg_c in component dMg_c_dt (molar)');
    LEGEND_RATES(:,8) = strpad('d/dt K_c in component dK_c_dt (molar)');
    LEGEND_RATES(:,55) = strpad('d/dt PI_c in component dPI_c_dt (molar)');
    LEGEND_RATES(:,5) = strpad('d/dt ATP_c in component dATP_c_dt (molar)');
    LEGEND_RATES(:,10) = strpad('d/dt ADP_c in component dADP_c_dt (molar)');
    LEGEND_RATES(:,58) = strpad('d/dt GLC_c in component dGLC_c_dt (molar)');
    LEGEND_RATES(:,59) = strpad('d/dt G6P_c in component dG6P_c_dt (molar)');
    LEGEND_RATES(:,60) = strpad('d/dt PCr_c in component dPCr_c_dt (molar)');
    LEGEND_RATES(:,53) = strpad('d/dt AMP_c in component dAMP_c_dt (molar)');
    LEGEND_RATES(:,2) = strpad('d/dt H_x in component dH_x_dt (molar)');
    LEGEND_RATES(:,3) = strpad('d/dt Mg_x in component dMg_x_dt (molar)');
    LEGEND_RATES(:,4) = strpad('d/dt K_x in component dK_x_dt (molar)');
    LEGEND_STATES  = LEGEND_STATES';
    LEGEND_ALGEBRAIC = LEGEND_ALGEBRAIC';
    LEGEND_RATES = LEGEND_RATES';
    LEGEND_CONSTANTS = LEGEND_CONSTANTS';
end

function [STATES, CONSTANTS] = initConsts()
    VOI = 0; CONSTANTS = []; STATES = []; ALGEBRAIC = [];
    CONSTANTS(:,1) = 1;
    CONSTANTS(:,2) = 2;
    CONSTANTS(:,3) = 0.096484;
    CONSTANTS(:,4) = 0.651384;
    CONSTANTS(:,5) = 0.072376;
    CONSTANTS(:,6) = 5.99;
    CONSTANTS(:,7) = 2.725e5;
    CONSTANTS(:,8) = 3;
    CONSTANTS(:,9) = 327;
    CONSTANTS(:,10) = 85;
    CONSTANTS(:,11) = 85;
    CONSTANTS(:,12) = 3.5e-6;
    CONSTANTS(:,13) = 0.6;
    CONSTANTS(:,14) = 1.2e-4;
    CONSTANTS(:,15) = 2.97e-3;
    CONSTANTS(:,16) = 1.35e-3;
    CONSTANTS(:,17) = 2.7e-3;
    CONSTANTS(:,18) = 10e-3;
    CONSTANTS(:,19) = 42.7e-3;
    CONSTANTS(:,20) = 21.4e-3;
    CONSTANTS(:,21) = 1e-4;
    CONSTANTS(:,22) = 6.75e-6;
    CONSTANTS(:,23) = 0.205093643103884;
    CONSTANTS(:,24) = 9.82530750761898;
    CONSTANTS(:,25) = 0.027665411180519;
    CONSTANTS(:,26) = 0.492047553306503;
    CONSTANTS(:,27) = 0.087630003425524;
    CONSTANTS(:,28) = 6.04e-7;
    CONSTANTS(:,29) = 0.448520452008943;
    CONSTANTS(:,30) = 0.08578867843174;
    CONSTANTS(:,31) = 0.00707728025683;
    CONSTANTS(:,32) = 0.077060004692409;
    CONSTANTS(:,33) = 0.025543515147654;
    CONSTANTS(:,34) = 4.551391837332316;
    CONSTANTS(:,35) = 2.996779273769102e8;
    CONSTANTS(:,36) = 2.725114215177993e8;
    CONSTANTS(:,37) = 80.432110322866379;
    CONSTANTS(:,38) = 0.282091450119444;
    CONSTANTS(:,39) = 96.995618247389871;
    CONSTANTS(:,40) = 16.062441848715917;
    CONSTANTS(:,41) = 5.865621634754698e-5;
    CONSTANTS(:,42) = 3.236481869934019e4;
    CONSTANTS(:,43) = 0.790813746051711;
    CONSTANTS(:,44) = 3.248913595933763e-5;
    CONSTANTS(:,45) = 0.003474542589709;
    CONSTANTS(:,46) = 1.076088893810943e-4;
    CONSTANTS(:,47) = 8.120393248498723e3;
    CONSTANTS(:,48) = 0.007522124149343;
    CONSTANTS(:,49) = 3.335568500871044e7;
    CONSTANTS(:,50) = 0.00161251709168;
    CONSTANTS(:,51) = 4.758023126099895e6;
    CONSTANTS(:,52) = 3.02092849850504e2;
    CONSTANTS(:,53) = 0;
    CONSTANTS(:,54) = -235.74;
    CONSTANTS(:,55) = 16.4;
    CONSTANTS(:,56) = 39.31;
    CONSTANTS(:,57) = 18.1;
    CONSTANTS(:,58) = -23.3;
    CONSTANTS(:,59) = 65.17;
    CONSTANTS(:,60) = -2771;
    CONSTANTS(:,61) = -1903.96;
    CONSTANTS(:,62) = -1034.66;
    CONSTANTS(:,63) = -2771;
    CONSTANTS(:,64) = -1903.96;
    CONSTANTS(:,65) = -27.41;
    CONSTANTS(:,66) = -6.52;
    CONSTANTS(:,67) = -1098.27;
    CONSTANTS(:,68) = 0;
    CONSTANTS(:,69) = -252.68;
    CONSTANTS(:,70) = -67.6;
    CONSTANTS(:,71) = 19.55;
    CONSTANTS(:,72) = -0.72;
    CONSTANTS(:,73) = -178.19;
    CONSTANTS(:,74) = -794.74;
    CONSTANTS(:,75) = -1165.59;
    CONSTANTS(:,76) = -1158.94;
    CONSTANTS(:,77) = -793.41;
    CONSTANTS(:,78) = -507.55;
    CONSTANTS(:,79) = -690.44;
    CONSTANTS(:,80) = -603.32;
    CONSTANTS(:,81) = -842.66;
    CONSTANTS(:,82) = -692.26;
    CONSTANTS(:,83) = -692.4;
    CONSTANTS(:,84) = -530.71;
    CONSTANTS(:,85) = -470.82;
    CONSTANTS(:,86) = -907.21;
    CONSTANTS(:,87) = -1758.87;
    STATES(:,1) = 1.025e-11;
    STATES(:,2) = 5.7919e-8;
    STATES(:,3) = 0.0046388;
    STATES(:,4) = 0.086878;
    STATES(:,5) = 1e-9;
    STATES(:,6) = 1e-7;
    STATES(:,7) = 0.005;
    STATES(:,8) = 0.15;
    STATES(:,9) = 0.01;
    STATES(:,10) = 1e-9;
    STATES(:,11) = 1e-6;
    STATES(:,12) = 2.7157e-10;
    STATES(:,13) = 0.005;
    STATES(:,14) = 2.9918e-4;
    STATES(:,15) = 1.1574e-4;
    CONSTANTS(:,88) = 1e20;
    CONSTANTS(:,89) = 1e20;
    CONSTANTS(:,90) = 1e20;
    CONSTANTS(:,91) = 1e20;
    STATES(:,16) = 6.815100000000001e-5;
    STATES(:,17) = 1e-9;
    CONSTANTS(:,92) = 1e20;
    CONSTANTS(:,93) = 1e20;
    CONSTANTS(:,94) = 1e20;
    CONSTANTS(:,95) = 1e20;
    CONSTANTS(:,96) = 1e20;
    CONSTANTS(:,97) = 1e20;
    CONSTANTS(:,98) = 1e20;
    CONSTANTS(:,99) = 1e20;
    CONSTANTS(:,100) = 1e20;
    CONSTANTS(:,101) = 1e20;
    CONSTANTS(:,102) = 1e20;
    CONSTANTS(:,103) = 1e20;
    CONSTANTS(:,104) = 1e20;
    CONSTANTS(:,105) = 1e20;
    CONSTANTS(:,106) = 1e20;
    CONSTANTS(:,107) = 1e20;
    CONSTANTS(:,108) = 1e20;
    CONSTANTS(:,109) = 1e20;
    CONSTANTS(:,110) = 1e20;
    CONSTANTS(:,111) = 1e20;
    CONSTANTS(:,112) = 1e20;
    CONSTANTS(:,113) = 1e20;
    CONSTANTS(:,114) = 1e20;
    CONSTANTS(:,115) = 1e20;
    CONSTANTS(:,116) = 1e20;
    CONSTANTS(:,117) = 1e20;
    CONSTANTS(:,118) = 1e20;
    CONSTANTS(:,119) = 1e20;
    CONSTANTS(:,120) = 1e20;
    CONSTANTS(:,121) = 1e20;
    CONSTANTS(:,122) = 1e20;
    CONSTANTS(:,123) = 1e20;
    CONSTANTS(:,124) = 1e20;
    CONSTANTS(:,125) = 1e20;
    CONSTANTS(:,126) = 1e20;
    CONSTANTS(:,127) = 1e20;
    CONSTANTS(:,128) = 1e20;
    CONSTANTS(:,129) = 1e20;
    CONSTANTS(:,130) = 1e20;
    CONSTANTS(:,131) = 1e20;
    CONSTANTS(:,132) = 1e20;
    CONSTANTS(:,133) = 1e20;
    CONSTANTS(:,134) = 1e20;
    CONSTANTS(:,135) = 1e20;
    CONSTANTS(:,136) = 1e20;
    CONSTANTS(:,137) = 1e20;
    CONSTANTS(:,138) = 1e20;
    CONSTANTS(:,139) = 1e20;
    CONSTANTS(:,140) = 1e20;
    CONSTANTS(:,141) = 1e20;
    CONSTANTS(:,142) = 1e20;
    CONSTANTS(:,143) = 1e20;
    STATES(:,18) = 1.2076e-11;
    STATES(:,19) = 0.0030107;
    STATES(:,20) = 0.0214;
    STATES(:,21) = 2.5507e-12;
    STATES(:,22) = 8.484200000000001e-21;
    CONSTANTS(:,144) = 38.3e-6;
    CONSTANTS(:,145) = 9.9e-6;
    CONSTANTS(:,146) = 60.7e-6;
    CONSTANTS(:,147) = 40.2e-6;
    CONSTANTS(:,148) = 40e-6;
    CONSTANTS(:,149) = 1e-32;
    STATES(:,23) = 8.9236e-4;
    STATES(:,24) = 6.9627e-10;
    CONSTANTS(:,150) = 4e-6;
    CONSTANTS(:,151) = 14e-6;
    CONSTANTS(:,152) = 3.33e-6;
    CONSTANTS(:,153) = 1600e-6;
    CONSTANTS(:,154) = 900e-6;
    CONSTANTS(:,155) = 1800e-6;
    CONSTANTS(:,156) = 6000e-6;
    CONSTANTS(:,157) = 67e-6;
    CONSTANTS(:,158) = 140e-6;
    STATES(:,25) = 9.3558e-7;
    CONSTANTS(:,159) = 1161e-6;
    CONSTANTS(:,160) = 434e-6;
    STATES(:,26) = 4.7979e-11;
    CONSTANTS(:,161) = 74e-6;
    CONSTANTS(:,162) = 183e-6;
    CONSTANTS(:,163) = 23.8e-6;
    CONSTANTS(:,164) = 29e-6;
    CONSTANTS(:,165) = 3;
    CONSTANTS(:,166) = 91e-6;
    CONSTANTS(:,167) = 50e-6;
    CONSTANTS(:,168) = 1e-32;
    CONSTANTS(:,169) = 80e-6;
    CONSTANTS(:,170) = 55e-6;
    CONSTANTS(:,171) = 21e-6;
    CONSTANTS(:,172) = 6.9e-6;
    CONSTANTS(:,173) = 6.010527470838273e-7;
    CONSTANTS(:,174) = 1e3;
    CONSTANTS(:,175) = 50e-6;
    CONSTANTS(:,176) = 100e-6;
    CONSTANTS(:,177) = 1e-32;
    STATES(:,27) = 2.1681e-7;
    CONSTANTS(:,178) = 5.5e-6;
    CONSTANTS(:,179) = 100e-6;
    CONSTANTS(:,180) = 2000e-6;
    CONSTANTS(:,181) = 20e-6;
    CONSTANTS(:,182) = 3000e-6;
    CONSTANTS(:,183) = 11.1e-6;
    CONSTANTS(:,184) = 16e-6;
    CONSTANTS(:,185) = 55e-6;
    CONSTANTS(:,186) = 660e-6;
    CONSTANTS(:,187) = 20e-6;
    CONSTANTS(:,188) = 880e-6;
    CONSTANTS(:,189) = 11.1e-6;
    STATES(:,28) = 0;
    STATES(:,29) = 7.8476e-9;
    CONSTANTS(:,190) = 120e-6;
    CONSTANTS(:,191) = 1275e-6;
    CONSTANTS(:,192) = 467e-6;
    CONSTANTS(:,193) = 480e-6;
    CONSTANTS(:,194) = 2.45e-6;
    CONSTANTS(:,195) = 1200e-6;
    CONSTANTS(:,196) = 1.5e-6;
    CONSTANTS(:,197) = 450e-6;
    CONSTANTS(:,198) = 375e-6;
    STATES(:,30) = 3.2191e-8;
    CONSTANTS(:,199) = 44.7e-6;
    CONSTANTS(:,200) = 197.7e-6;
    CONSTANTS(:,201) = 3500e-6;
    CONSTANTS(:,202) = 40e-6;
    CONSTANTS(:,203) = 400e-6;
    CONSTANTS(:,204) = 80e-6;
    CONSTANTS(:,205) = 330e-6;
    CONSTANTS(:,206) = 279e-6;
    CONSTANTS(:,207) = 360e-6;
    CONSTANTS(:,208) = 5.5e-6;
    CONSTANTS(:,209) = 3.18e-6;
    CONSTANTS(:,210) = 90.55e-6;
    CONSTANTS(:,211) = 250e-6;
    CONSTANTS(:,212) = 6.128e-6;
    CONSTANTS(:,213) = 2.58e-6;
    CONSTANTS(:,214) = 183.2e-6;
    CONSTANTS(:,215) = 394.4e-6;
    CONSTANTS(:,216) = 420e-6;
    CONSTANTS(:,217) = 170e-6;
    CONSTANTS(:,218) = 143.6e-6;
    CONSTANTS(:,219) = 146.6e-6;
    CONSTANTS(:,220) = 156.5e-6;
    CONSTANTS(:,221) = 111e-6;
    CONSTANTS(:,222) = 100e-6;
    CONSTANTS(:,223) = 260e-6;
    CONSTANTS(:,224) = 278e-6;
    CONSTANTS(:,225) = 650e-6;
    CONSTANTS(:,226) = 1e-32;
    STATES(:,31) = 1.0217e-4;
    STATES(:,32) = 9.9311e-12;
    CONSTANTS(:,227) = 3480e-6;
    CONSTANTS(:,228) = 710e-6;
    CONSTANTS(:,229) = 50e-6;
    CONSTANTS(:,230) = 8400e-6;
    CONSTANTS(:,231) = 3900e-6;
    CONSTANTS(:,232) = 430e-6;
    CONSTANTS(:,233) = 88e-6;
    CONSTANTS(:,234) = 8900e-6;
    CONSTANTS(:,235) = 16.6e-3;
    CONSTANTS(:,236) = 1e-32;
    STATES(:,33) = 9.1804e-10;
    STATES(:,34) = 5.7539e-12;
    STATES(:,35) = 1e-9;
    CONSTANTS(:,237) = 0.3e-3;
    CONSTANTS(:,238) = 0.17e-3;
    CONSTANTS(:,239) = 1.4e-3;
    CONSTANTS(:,240) = 0.7e-3;
    STATES(:,36) = 9.999500000000001e-10;
    CONSTANTS(:,241) = 28e-6;
    CONSTANTS(:,242) = 2.8e-3;
    CONSTANTS(:,243) = 180e-6;
    CONSTANTS(:,244) = 1.6e-3;
    STATES(:,37) = 4.4404e-9;
    CONSTANTS(:,245) = 1.8;
    CONSTANTS(:,246) = 6.5;
    STATES(:,38) = -14.087;
    STATES(:,39) = 1e-9;
    STATES(:,40) = 1e-9;
    STATES(:,41) = 1e-9;
    STATES(:,42) = 1e-9;
    STATES(:,43) = 1e-9;
    STATES(:,44) = 1e-9;
    STATES(:,45) = 5.754e-12;
    STATES(:,46) = 4.4404e-9;
    CONSTANTS(:,247) = 0;
    STATES(:,47) = 1e-9;
    STATES(:,48) = 1e-9;
    STATES(:,49) = 1e-9;
    STATES(:,50) = 1e-9;
    STATES(:,51) = 1e-9;
    STATES(:,52) = 1e-9;
    STATES(:,53) = 1e-9;
    STATES(:,54) = 1e-9;
    STATES(:,55) = 1.1574e-4;
    STATES(:,56) = 9.6996e-8;
    CONSTANTS(:,248) = 1e-32;
    STATES(:,57) = 6.499999999999999e-5;
    CONSTANTS(:,249) = 1e-32;
    CONSTANTS(:,250) = 1e-9;
    CONSTANTS(:,251) = 1e-9;
    STATES(:,58) = 1e-9;
    STATES(:,59) = 1e-9;
    CONSTANTS(:,252) = 1e-3;
    CONSTANTS(:,253) = 47e-6;
    CONSTANTS(:,254) = 47e-6;
    CONSTANTS(:,255) = 1e-3;
    CONSTANTS(:,256) = 1e-3;
    CONSTANTS(:,257) = 47e-6;
    CONSTANTS(:,258) = 47e-6;
    CONSTANTS(:,259) = 1e-3;
    CONSTANTS(:,260) = 10e-6;
    CONSTANTS(:,261) = 1e7;
    STATES(:,60) = 1e-9;
    CONSTANTS(:,262) = 1e7;
    CONSTANTS(:,263) = 0;
    CONSTANTS(:,264) = 0;
    CONSTANTS(:,265) = 0;
    CONSTANTS(:,266) = 1;
    CONSTANTS(:,267) = 1;
    CONSTANTS(:,268) = 1;
    CONSTANTS(:,269) = 1;
    CONSTANTS(:,270) = 0;
    CONSTANTS(:,271) = 0;
    CONSTANTS(:,272) = 0.02;
    CONSTANTS(:,273) = 1e-7;
    CONSTANTS(:,274) = piecewise({CONSTANTS(:,1)==1.00000&CONSTANTS(:,2)==3.00000, 0.107108 }, 0.00000);
    CONSTANTS(:,275) = piecewise({CONSTANTS(:,1)==1.00000, ( 8.31400.*(28.0000+273.150))./1000.00 }, ( 8.31400.*(37.0000+273.150))./1000.00);
    CONSTANTS(:,276) = 1.00000;
    CONSTANTS(:,277) = 1.00000;
    CONSTANTS(:,278) = 1.00000;
    CONSTANTS(:,279) = 1.00000;
    CONSTANTS(:,280) = 1.00000;
    CONSTANTS(:,281) = 1.00000;
    CONSTANTS(:,282) = 1.00000;
    CONSTANTS(:,283) = 1.00000;
    CONSTANTS(:,284) = 1.00000;
    CONSTANTS(:,285) = 1.00000;
    CONSTANTS(:,286) =  1.00000.*power(10.0000,  - CONSTANTS(:,246));
    CONSTANTS(:,287) = piecewise({CONSTANTS(:,1)==1.00000, 80.0000 , CONSTANTS(:,1)==2.00000, 800.000 , CONSTANTS(:,1)==3.00000||CONSTANTS(:,1)==4.00000,  0.807000.*1.04400 }, 80.0000);
    CONSTANTS(:,288) = ((CONSTANTS(:,57)+CONSTANTS(:,58)) - CONSTANTS(:,56)) - CONSTANTS(:,59);
    CONSTANTS(:,289) = ((CONSTANTS(:,59)+ 2.00000.*CONSTANTS(:,65)) - CONSTANTS(:,58)) -  2.00000.*CONSTANTS(:,66);
    CONSTANTS(:,290) = (( 2.00000.*CONSTANTS(:,66)+CONSTANTS(:,54)) -  2.00000.*CONSTANTS(:,65)) -  0.500000.*CONSTANTS(:,55);
    CONSTANTS(:,291) = ((CONSTANTS(:,60)+CONSTANTS(:,54)) - CONSTANTS(:,61)) - CONSTANTS(:,67);
    CONSTANTS(:,292) = CONSTANTS(:,21)./2.00000;
    CONSTANTS(:,293) = 0.00000;
    CONSTANTS(:,294) = 0.00000;
    CONSTANTS(:,295) = 0.00000;
    CONSTANTS(:,296) = 0.00000;
    CONSTANTS(:,382) = 0.00000;
    CONSTANTS(:,383) = 0.00000;
    CONSTANTS(:,384) = 0.00000;
    CONSTANTS(:,385) = 0.00000;
    CONSTANTS(:,386) = 0.00000;
    CONSTANTS(:,387) = 0.00000;
    CONSTANTS(:,297) =  1.00000.*power(10.0000,  - 6.59000);
    CONSTANTS(:,298) = ((CONSTANTS(:,64)+CONSTANTS(:,60)) - CONSTANTS(:,63)) - CONSTANTS(:,61);
    CONSTANTS(:,299) = ((((CONSTANTS(:,84)+CONSTANTS(:,73)+CONSTANTS(:,56)) - CONSTANTS(:,85)) - CONSTANTS(:,72)) - CONSTANTS(:,57)) - CONSTANTS(:,54);
    CONSTANTS(:,300) = (((CONSTANTS(:,72)+CONSTANTS(:,75)) - CONSTANTS(:,73)) - CONSTANTS(:,74)) - CONSTANTS(:,54);
    CONSTANTS(:,301) = CONSTANTS(:,76) - CONSTANTS(:,75);
    CONSTANTS(:,302) = (((CONSTANTS(:,77)+CONSTANTS(:,56)+CONSTANTS(:,84)) - CONSTANTS(:,76)) - CONSTANTS(:,57)) - CONSTANTS(:,54);
    CONSTANTS(:,303) = ((((CONSTANTS(:,84)+CONSTANTS(:,78)+CONSTANTS(:,56)) - CONSTANTS(:,77)) - CONSTANTS(:,72)) - CONSTANTS(:,57)) - CONSTANTS(:,54);
    CONSTANTS(:,304) = (((CONSTANTS(:,72)+CONSTANTS(:,79)+CONSTANTS(:,63)) - CONSTANTS(:,64)) - CONSTANTS(:,78)) - CONSTANTS(:,67);
    CONSTANTS(:,305) = 1.00000;
    CONSTANTS(:,306) = ((CONSTANTS(:,74)+CONSTANTS(:,56)) - CONSTANTS(:,57)) - CONSTANTS(:,81);
    CONSTANTS(:,307) = 1.00000;
    CONSTANTS(:,308) = piecewise({CONSTANTS(:,1)==3.00000, 0.288200 , CONSTANTS(:,1)==4.00000, 0.0560000 }, 0.00000);
    CONSTANTS(:,309) = ((CONSTANTS(:,87)+CONSTANTS(:,61)) - CONSTANTS(:,86)) - CONSTANTS(:,60);
    CONSTANTS(:,310) = exp(50.7800./CONSTANTS(:,275))./1.00000;
    CONSTANTS(:,311) = CONSTANTS(:,21) - CONSTANTS(:,292);
    CONSTANTS(:,312) =  1.00000.*power(10.0000,  - 3.82000);
    CONSTANTS(:,313) = exp( - CONSTANTS(:,298)./CONSTANTS(:,275));
    CONSTANTS(:,314) = CONSTANTS(:,297);
    CONSTANTS(:,315) =  1.00000.*exp( - CONSTANTS(:,299)./CONSTANTS(:,275));
    CONSTANTS(:,316) =  1.00000.*exp( - CONSTANTS(:,300)./CONSTANTS(:,275));
    CONSTANTS(:,317) = exp( - CONSTANTS(:,301)./CONSTANTS(:,275));
    CONSTANTS(:,318) =  1.00000.*exp( - CONSTANTS(:,302)./CONSTANTS(:,275));
    CONSTANTS(:,319) =  1.00000.*exp( - CONSTANTS(:,303)./CONSTANTS(:,275));
    CONSTANTS(:,320) =  1.00000.*exp( - CONSTANTS(:,304)./CONSTANTS(:,275));
    CONSTANTS(:,321) = 1.00000;
    CONSTANTS(:,322) =  1.00000.*exp( - CONSTANTS(:,306)./CONSTANTS(:,275));
    CONSTANTS(:,323) = ((CONSTANTS(:,74)+CONSTANTS(:,83)) - CONSTANTS(:,82)) - CONSTANTS(:,77);
    CONSTANTS(:,324) = piecewise({CONSTANTS(:,1)==3.00000, 0.680100 , CONSTANTS(:,1)==4.00000, 0.894000 }, 1.00000);
    CONSTANTS(:,325) =  1.00000.*exp( - CONSTANTS(:,309)./CONSTANTS(:,275));
    CONSTANTS(:,326) =  1.00000.*power(10.0000,  - 1.01300);
    CONSTANTS(:,327) = CONSTANTS(:,312);
    CONSTANTS(:,328) = ((CONSTANTS(:,58)+CONSTANTS(:,80)) - CONSTANTS(:,79)) - CONSTANTS(:,59);
    CONSTANTS(:,329) = exp( - CONSTANTS(:,323)./CONSTANTS(:,275));
    CONSTANTS(:,330) = piecewise({CONSTANTS(:,2)==0.00000||CONSTANTS(:,1)==5.00000, 0.00000 , CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,331) = piecewise({CONSTANTS(:,2)==0.00000||CONSTANTS(:,1)==5.00000, 0.00000 , CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,332) = piecewise({CONSTANTS(:,2)==0.00000||CONSTANTS(:,1)==3.00000||CONSTANTS(:,1)==4.00000||CONSTANTS(:,1)==5.00000, 0.00000 , CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,333) =  1.00000.*power(10.0000,  - 6.42000);
    CONSTANTS(:,334) = CONSTANTS(:,326);
    CONSTANTS(:,335) = exp( - CONSTANTS(:,328)./CONSTANTS(:,275));
    CONSTANTS(:,336) = piecewise({CONSTANTS(:,2)==0.00000||CONSTANTS(:,1)==5.00000, 0.00000 , CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,337) =  1.00000.*power(10.0000,  - 2.79000);
    CONSTANTS(:,338) = CONSTANTS(:,333);
    CONSTANTS(:,339) = (CONSTANTS(:,81) - CONSTANTS(:,80)) - CONSTANTS(:,54);
    CONSTANTS(:,340) = piecewise({CONSTANTS(:,2)==0.00000||CONSTANTS(:,1)==5.00000, 0.00000 , CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,341) =  1.00000.*power(10.0000,  - 0.882000);
    CONSTANTS(:,342) = CONSTANTS(:,337);
    CONSTANTS(:,343) = exp( - CONSTANTS(:,339)./CONSTANTS(:,275));
    CONSTANTS(:,344) = piecewise({CONSTANTS(:,2)==0.00000||CONSTANTS(:,1)==5.00000, 0.00000 , CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,345) =  1.00000.*power(10.0000,  - 6.22000);
    CONSTANTS(:,346) = CONSTANTS(:,341);
    CONSTANTS(:,347) = piecewise({CONSTANTS(:,2)==0.00000||CONSTANTS(:,1)==5.00000, 0.00000 , CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,348) =  1.00000.*power(10.0000,  - 1.86000);
    CONSTANTS(:,349) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,350) =  1.00000.*power(10.0000,  - 0.621500);
    CONSTANTS(:,351) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, 1.00000 }, CONSTANTS(:,308)./CONSTANTS(:,324));
    CONSTANTS(:,352) =  1.00000.*power(10.0000,  - 6.71000);
    CONSTANTS(:,353) = CONSTANTS(:,308)./CONSTANTS(:,324);
    CONSTANTS(:,354) =  1.00000.*power(10.0000,  - 1.69000);
    CONSTANTS(:,355) = CONSTANTS(:,308)./CONSTANTS(:,324);
    CONSTANTS(:,356) =  1.00000.*power(10.0000, 0.00740000);
    CONSTANTS(:,357) = CONSTANTS(:,308)./CONSTANTS(:,324);
    CONSTANTS(:,358) =  1.00000.*power(10.0000,  - 8.13000);
    CONSTANTS(:,359) = CONSTANTS(:,308)./CONSTANTS(:,324);
    CONSTANTS(:,360) =  1.00000.*power(10.0000,  - 0.862900);
    CONSTANTS(:,361) =  1.00000.*power(10.0000,  - 5.63000);
    CONSTANTS(:,362) =  1.00000.*power(10.0000,  - 3.37000);
    CONSTANTS(:,363) =  1.00000.*power(10.0000,  - 0.339000);
    CONSTANTS(:,364) =  1.00000.*power(10.0000,  - 5.64000);
    CONSTANTS(:,365) =  1.00000.*power(10.0000,  - 2.46000);
    CONSTANTS(:,366) =  1.00000.*power(10.0000,  - 3.96000);
    CONSTANTS(:,367) =  1.00000.*power(10.0000,  - 5.13000);
    CONSTANTS(:,368) =  1.00000.*power(10.0000,  - 1.17000);
    CONSTANTS(:,369) =  1.00000.*power(10.0000,  - 0.352500);
    CONSTANTS(:,370) =  1.00000.*power(10.0000,  - 4.10000);
    CONSTANTS(:,371) =  1.00000.*power(10.0000,  - 4.75000);
    CONSTANTS(:,372) =  1.00000.*power(10.0000,  - 1.55000);
    CONSTANTS(:,373) =  1.00000.*power(10.0000, 0.107000);
    CONSTANTS(:,374) =  1.00000.*power(10.0000,  - 9.82000);
    CONSTANTS(:,375) =  1.00000.*power(10.0000,  - 1.02000);
    CONSTANTS(:,376) =  1.00000.*power(10.0000,  - 5.91000);
    CONSTANTS(:,377) =  1.00000.*power(10.0000,  - 4.06000);
    CONSTANTS(:,378) =  1.00000.*power(10.0000,  - 1.82000);
    CONSTANTS(:,379) =  1.00000.*power(10.0000,  - 3.65000);
    CONSTANTS(:,380) =  1.00000.*power(10.0000,  - 2.32000);
    CONSTANTS(:,381) = (CONSTANTS(:,60)+CONSTANTS(:,62)) -  2.00000.*CONSTANTS(:,61);
    CONSTANTS(:,382) = exp( - CONSTANTS(:,381)./CONSTANTS(:,275));
    if (isempty(STATES)), warning('Initial values for states not set');, end
end

function [RATES, ALGEBRAIC] = computeRates(VOI, STATES, CONSTANTS)
    global algebraicVariableCount;
    statesSize = size(STATES);
    statesColumnCount = statesSize(2);
    if ( statesColumnCount == 1)
        STATES = STATES';
        ALGEBRAIC = zeros(1, algebraicVariableCount);
        utilOnes = 1;
    else
        statesRowCount = statesSize(1);
        ALGEBRAIC = zeros(statesRowCount, algebraicVariableCount);
        RATES = zeros(statesRowCount, statesColumnCount);
        utilOnes = ones(statesRowCount, 1);
    end
    RATES(:,11) = CONSTANTS(:,382);
    RATES(:,57) = CONSTANTS(:,383);
    RATES(:,20) = CONSTANTS(:,384);
    RATES(:,6) = CONSTANTS(:,385);
    RATES(:,7) = CONSTANTS(:,386);
    RATES(:,8) = CONSTANTS(:,387);
    ALGEBRAIC(:,17) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,49) - STATES(:,50));
    RATES(:,50) = ALGEBRAIC(:,17)./CONSTANTS(:,5);
    ALGEBRAIC(:,16) =  CONSTANTS(:,6).*CONSTANTS(:,247).*(STATES(:,47) - STATES(:,48));
    RATES(:,48) = ALGEBRAIC(:,16)./CONSTANTS(:,5);
    RATES(:,49) = (  - CONSTANTS(:,331).*ALGEBRAIC(:,17))./CONSTANTS(:,287);
    RATES(:,47) = (  - CONSTANTS(:,330).*ALGEBRAIC(:,16))./CONSTANTS(:,287);
    ALGEBRAIC(:,1) = 1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,3)./CONSTANTS(:,312)+STATES(:,4)./CONSTANTS(:,326);
    ALGEBRAIC(:,30) = ALGEBRAIC(:,1);
    ALGEBRAIC(:,22) = 1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,3)./CONSTANTS(:,337)+STATES(:,4)./CONSTANTS(:,341);
    ALGEBRAIC(:,32) = ALGEBRAIC(:,22);
    ALGEBRAIC(:,34) = 1.00000+STATES(:,2)./CONSTANTS(:,352)+STATES(:,3)./CONSTANTS(:,354)+STATES(:,4)./CONSTANTS(:,356);
    ALGEBRAIC(:,36) = 1.00000+STATES(:,2)./CONSTANTS(:,358);
    ALGEBRAIC(:,42) = 1.00000+STATES(:,2)./CONSTANTS(:,366);
    ALGEBRAIC(:,43) = 1.00000+STATES(:,2)./CONSTANTS(:,367)+STATES(:,3)./CONSTANTS(:,368)+STATES(:,4)./CONSTANTS(:,369);
    ALGEBRAIC(:,65) = ( CONSTANTS(:,320).*ALGEBRAIC(:,36).*ALGEBRAIC(:,43).*ALGEBRAIC(:,30))./( STATES(:,2).*ALGEBRAIC(:,32).*ALGEBRAIC(:,42).*ALGEBRAIC(:,34));
    ALGEBRAIC(:,66) = ( CONSTANTS(:,29).*CONSTANTS(:,187).*CONSTANTS(:,182).*CONSTANTS(:,183))./( ALGEBRAIC(:,65).*CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,186));
    ALGEBRAIC(:,67) = ( CONSTANTS(:,29).*ALGEBRAIC(:,66).*( STATES(:,13).*STATES(:,24).*STATES(:,14) - ( STATES(:,19).*STATES(:,27).*STATES(:,12))./ALGEBRAIC(:,65)))./( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,186)+ ALGEBRAIC(:,66).*CONSTANTS(:,179).*CONSTANTS(:,186).*STATES(:,13)+ ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14)+ ALGEBRAIC(:,66).*CONSTANTS(:,186).*STATES(:,13).*STATES(:,24)+ ALGEBRAIC(:,66).*CONSTANTS(:,185).*STATES(:,13).*STATES(:,14)+ ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14)+ ALGEBRAIC(:,66).*STATES(:,13).*STATES(:,24).*STATES(:,14)+( CONSTANTS(:,29).*CONSTANTS(:,183).*CONSTANTS(:,188).*STATES(:,19))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,182).*CONSTANTS(:,187).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,19).*STATES(:,27))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,188).*STATES(:,19).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,187).*STATES(:,27).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*STATES(:,19).*STATES(:,27).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,188).*CONSTANTS(:,183).*STATES(:,13).*STATES(:,19))./( CONSTANTS(:,178).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14).*STATES(:,12))./CONSTANTS(:,183)+( CONSTANTS(:,29).*CONSTANTS(:,188).*CONSTANTS(:,183).*STATES(:,13).*STATES(:,24).*STATES(:,19))./( CONSTANTS(:,178).*CONSTANTS(:,179).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14).*STATES(:,12))./CONSTANTS(:,183)+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,19).*STATES(:,27))./( CONSTANTS(:,178).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,182).*CONSTANTS(:,183))+( CONSTANTS(:,29).*CONSTANTS(:,183).*CONSTANTS(:,188).*STATES(:,13).*STATES(:,24).*STATES(:,14).*STATES(:,19))./( CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,180).*ALGEBRAIC(:,65))+( CONSTANTS(:,29).*CONSTANTS(:,181).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,24).*STATES(:,14).*STATES(:,27))./( CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,180).*ALGEBRAIC(:,65))+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,24).*STATES(:,19).*STATES(:,27))./( CONSTANTS(:,178).*CONSTANTS(:,179).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,182).*CONSTANTS(:,183))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*CONSTANTS(:,180).*STATES(:,24).*STATES(:,19).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,181).*CONSTANTS(:,182).*CONSTANTS(:,183))+( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14).*STATES(:,19).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,181).*CONSTANTS(:,182).*CONSTANTS(:,183))+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,24).*STATES(:,14).*STATES(:,19).*STATES(:,27))./( CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,180).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14).*STATES(:,19).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,181).*CONSTANTS(:,182).*CONSTANTS(:,183)));
    ALGEBRAIC(:,78) = ( CONSTANTS(:,313).*ALGEBRAIC(:,32).*ALGEBRAIC(:,1))./( ALGEBRAIC(:,30).*ALGEBRAIC(:,22));
    ALGEBRAIC(:,79) = ( CONSTANTS(:,33).*CONSTANTS(:,224).*CONSTANTS(:,219))./( ALGEBRAIC(:,78).*CONSTANTS(:,217).*CONSTANTS(:,222));
    ALGEBRAIC(:,27) = 1.00000+STATES(:,2)./CONSTANTS(:,345)+STATES(:,3)./CONSTANTS(:,348)+STATES(:,4)./CONSTANTS(:,350);
    ALGEBRAIC(:,29) = ( STATES(:,11).*(1.00000+STATES(:,2)./CONSTANTS(:,345)))./ALGEBRAIC(:,27);
    ALGEBRAIC(:,80) = 1.00000+ALGEBRAIC(:,29)./CONSTANTS(:,225);
    ALGEBRAIC(:,81) = piecewise({STATES(:,12)>CONSTANTS(:,226)&STATES(:,9)>CONSTANTS(:,226), ( (( CONSTANTS(:,33).*ALGEBRAIC(:,79))./ALGEBRAIC(:,80)).*STATES(:,12).*STATES(:,9))./( ALGEBRAIC(:,79).*CONSTANTS(:,222).*STATES(:,12)+ ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,13))./ALGEBRAIC(:,78)+( CONSTANTS(:,33).*CONSTANTS(:,223).*STATES(:,1))./ALGEBRAIC(:,78)+ ALGEBRAIC(:,79).*STATES(:,12).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,12).*STATES(:,13))./( ALGEBRAIC(:,78).*CONSTANTS(:,217))+( CONSTANTS(:,33).*STATES(:,13).*STATES(:,1))./ALGEBRAIC(:,78)+( ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9).*STATES(:,1))./CONSTANTS(:,220)) }, 0.00000);
    ALGEBRAIC(:,82) = piecewise({STATES(:,13)>CONSTANTS(:,226)&STATES(:,1)>CONSTANTS(:,226), ALGEBRAIC(:,81) - (( (( CONSTANTS(:,33).*ALGEBRAIC(:,79))./ALGEBRAIC(:,80)).*STATES(:,13).*STATES(:,1))./ALGEBRAIC(:,78))./( ALGEBRAIC(:,79).*CONSTANTS(:,222).*STATES(:,12)+ ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,13))./ALGEBRAIC(:,78)+( CONSTANTS(:,33).*CONSTANTS(:,223).*STATES(:,1))./ALGEBRAIC(:,78)+ ALGEBRAIC(:,79).*STATES(:,12).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,12).*STATES(:,13))./( ALGEBRAIC(:,78).*CONSTANTS(:,217))+( CONSTANTS(:,33).*STATES(:,13).*STATES(:,1))./ALGEBRAIC(:,78)+( ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9).*STATES(:,1))./CONSTANTS(:,220)) }, ALGEBRAIC(:,81));
    RATES(:,12) = (ALGEBRAIC(:,67) - ALGEBRAIC(:,82))./CONSTANTS(:,4);
    RATES(:,13) = ( - ALGEBRAIC(:,67)+ALGEBRAIC(:,82))./CONSTANTS(:,4);
    ALGEBRAIC(:,88) = ( CONSTANTS(:,38).*( STATES(:,36).*STATES(:,30) -  STATES(:,26).*STATES(:,35)))./( CONSTANTS(:,237).*CONSTANTS(:,240).*(2.00000+STATES(:,35)./CONSTANTS(:,239)+STATES(:,30)./CONSTANTS(:,240)+STATES(:,36)./CONSTANTS(:,237)+STATES(:,26)./CONSTANTS(:,238)+( STATES(:,35).*STATES(:,26))./( CONSTANTS(:,239).*CONSTANTS(:,238))+( STATES(:,30).*STATES(:,36))./( CONSTANTS(:,240).*CONSTANTS(:,237))));
    ALGEBRAIC(:,90) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,43) - STATES(:,36));
    RATES(:,36) = ( - ALGEBRAIC(:,88)+ALGEBRAIC(:,90))./CONSTANTS(:,5);
    ALGEBRAIC(:,89) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,40) - STATES(:,33));
    RATES(:,40) = (  - CONSTANTS(:,332).*ALGEBRAIC(:,89))./CONSTANTS(:,287);
    RATES(:,43) = (  - CONSTANTS(:,340).*ALGEBRAIC(:,90))./CONSTANTS(:,287);
    ALGEBRAIC(:,91) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,41) - STATES(:,17));
    RATES(:,41) = (  - CONSTANTS(:,336).*ALGEBRAIC(:,91))./CONSTANTS(:,287);
    ALGEBRAIC(:,92) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,42) - STATES(:,35));
    RATES(:,42) = (  - CONSTANTS(:,347).*ALGEBRAIC(:,92))./CONSTANTS(:,287);
    ALGEBRAIC(:,93) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,44) - STATES(:,39));
    RATES(:,44) = (  - CONSTANTS(:,344).*ALGEBRAIC(:,93))./CONSTANTS(:,287);
    ALGEBRAIC(:,94) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,45) - STATES(:,34));
    RATES(:,45) = (  - CONSTANTS(:,349).*ALGEBRAIC(:,94))./CONSTANTS(:,287);
    ALGEBRAIC(:,95) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,46) - STATES(:,37));
    RATES(:,46) = (  - CONSTANTS(:,351).*ALGEBRAIC(:,95))./CONSTANTS(:,287);
    ALGEBRAIC(:,99) =  CONSTANTS(:,6).*CONSTANTS(:,9).*(STATES(:,55) - STATES(:,15));
    RATES(:,55) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, (  - ALGEBRAIC(:,99).*CONSTANTS(:,266))./CONSTANTS(:,287) }, ( CONSTANTS(:,353).*ALGEBRAIC(:,99)+CONSTANTS(:,53))./CONSTANTS(:,287));
    ALGEBRAIC(:,18) = 1.00000+STATES(:,6)./CONSTANTS(:,297)+STATES(:,7)./CONSTANTS(:,312)+STATES(:,8)./CONSTANTS(:,326);
    ALGEBRAIC(:,21) = (( STATES(:,5).*STATES(:,7))./CONSTANTS(:,312))./ALGEBRAIC(:,18);
    ALGEBRAIC(:,23) = 1.00000+STATES(:,6)./CONSTANTS(:,333)+STATES(:,7)./CONSTANTS(:,337)+STATES(:,8)./CONSTANTS(:,341);
    ALGEBRAIC(:,26) = (( STATES(:,10).*STATES(:,7))./CONSTANTS(:,337))./ALGEBRAIC(:,23);
    ALGEBRAIC(:,48) = 1.00000+STATES(:,6)./CONSTANTS(:,376);
    ALGEBRAIC(:,100) = ( CONSTANTS(:,325).*ALGEBRAIC(:,48).*ALGEBRAIC(:,23))./( STATES(:,6).*ALGEBRAIC(:,18).*CONSTANTS(:,276));
    ALGEBRAIC(:,102) = ( ALGEBRAIC(:,100).*CONSTANTS(:,312).*ALGEBRAIC(:,18))./( CONSTANTS(:,337).*ALGEBRAIC(:,23));
    ALGEBRAIC(:,104) = ( (CONSTANTS(:,274)./( CONSTANTS(:,253).*CONSTANTS(:,256))).*( ALGEBRAIC(:,21).*STATES(:,58) - ( STATES(:,59).*ALGEBRAIC(:,26))./ALGEBRAIC(:,102)))./(1.00000+ALGEBRAIC(:,21)./CONSTANTS(:,252)+STATES(:,58)./CONSTANTS(:,253)+( ALGEBRAIC(:,21).*STATES(:,58))./( CONSTANTS(:,253).*CONSTANTS(:,256))+STATES(:,59)./CONSTANTS(:,254)+ALGEBRAIC(:,26)./CONSTANTS(:,255)+( STATES(:,59).*ALGEBRAIC(:,26))./( CONSTANTS(:,255).*CONSTANTS(:,258))+( STATES(:,59).*STATES(:,58))./( CONSTANTS(:,260).*CONSTANTS(:,253)));
    RATES(:,58) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, (  - (1.00000+CONSTANTS(:,287)./1.00000).*ALGEBRAIC(:,104))./CONSTANTS(:,287) }, 0.00000);
    RATES(:,59) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, ( (1.00000+CONSTANTS(:,287)./1.00000).*ALGEBRAIC(:,104))./CONSTANTS(:,287) }, 0.00000);
    ALGEBRAIC(:,106) = CONSTANTS(:,19) - STATES(:,60);
    ALGEBRAIC(:,108) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000, 0.00000 },  CONSTANTS(:,261).*( (( CONSTANTS(:,310).*STATES(:,10))./ALGEBRAIC(:,23)).*STATES(:,60).*STATES(:,6) -  (STATES(:,5)./ALGEBRAIC(:,18)).*ALGEBRAIC(:,106)));
    RATES(:,60) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, 0.00000 },  - ALGEBRAIC(:,108)./CONSTANTS(:,287));
    ALGEBRAIC(:,97) =  CONSTANTS(:,6).*CONSTANTS(:,10).*(STATES(:,5) - STATES(:,52));
    ALGEBRAIC(:,28) = 1.00000+STATES(:,6)./CONSTANTS(:,345)+STATES(:,7)./CONSTANTS(:,348)+STATES(:,8)./CONSTANTS(:,350);
    ALGEBRAIC(:,110) = ( CONSTANTS(:,382).*ALGEBRAIC(:,18).*ALGEBRAIC(:,28))./( ALGEBRAIC(:,23).*ALGEBRAIC(:,23));
    ALGEBRAIC(:,112) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000, 0.00000 },  CONSTANTS(:,262).*( ALGEBRAIC(:,110).*STATES(:,10).*STATES(:,10) -  STATES(:,53).*STATES(:,5)));
    RATES(:,5) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, (  - (1.00000+CONSTANTS(:,287)./1.00000).*ALGEBRAIC(:,104) -  CONSTANTS(:,267).*ALGEBRAIC(:,97))./CONSTANTS(:,287) }, ((  - CONSTANTS(:,355).*ALGEBRAIC(:,97)+ALGEBRAIC(:,108)+ CONSTANTS(:,267).*ALGEBRAIC(:,112)) - CONSTANTS(:,53))./CONSTANTS(:,287));
    ALGEBRAIC(:,96) =  CONSTANTS(:,6).*CONSTANTS(:,10).*(STATES(:,10) - STATES(:,51));
    RATES(:,10) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, ( (1.00000+CONSTANTS(:,287)./1.00000).*ALGEBRAIC(:,104) -  CONSTANTS(:,268).*ALGEBRAIC(:,96))./CONSTANTS(:,287) }, (((  - CONSTANTS(:,357).*ALGEBRAIC(:,96) - ALGEBRAIC(:,108)) -  2.00000.*CONSTANTS(:,268).*ALGEBRAIC(:,112))+CONSTANTS(:,53))./CONSTANTS(:,287));
    ALGEBRAIC(:,98) =  CONSTANTS(:,6).*CONSTANTS(:,10).*(STATES(:,53) - STATES(:,54));
    RATES(:,53) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000||CONSTANTS(:,1)==6.00000, 0.00000 }, (  - CONSTANTS(:,359).*ALGEBRAIC(:,98)+ CONSTANTS(:,269).*ALGEBRAIC(:,112))./CONSTANTS(:,287));
    ALGEBRAIC(:,114) = CONSTANTS(:,15) - STATES(:,22);
    ALGEBRAIC(:,46) = 1.00000+STATES(:,2)./CONSTANTS(:,374);
    ALGEBRAIC(:,47) = 1.00000+STATES(:,3)./CONSTANTS(:,375);
    ALGEBRAIC(:,51) = ( CONSTANTS(:,315).*ALGEBRAIC(:,46).*CONSTANTS(:,279).*CONSTANTS(:,278))./( STATES(:,2).*ALGEBRAIC(:,47).*ALGEBRAIC(:,36).*CONSTANTS(:,277));
    ALGEBRAIC(:,52) = 1.00000+STATES(:,21)./CONSTANTS(:,147);
    ALGEBRAIC(:,53) = 1.00000+STATES(:,22)./CONSTANTS(:,148);
    ALGEBRAIC(:,115) = piecewise({STATES(:,18)>CONSTANTS(:,149)&STATES(:,19)>CONSTANTS(:,149)&ALGEBRAIC(:,114)>CONSTANTS(:,149), ( CONSTANTS(:,23).*( STATES(:,18).*STATES(:,19).*ALGEBRAIC(:,114) - ( STATES(:,20).*STATES(:,21).*STATES(:,22))./ALGEBRAIC(:,51)))./( CONSTANTS(:,146).*ALGEBRAIC(:,53).*STATES(:,18).*STATES(:,19)+ CONSTANTS(:,145).*ALGEBRAIC(:,52).*STATES(:,18).*ALGEBRAIC(:,114)+ CONSTANTS(:,144).*STATES(:,19).*ALGEBRAIC(:,114)+ STATES(:,18).*STATES(:,19).*ALGEBRAIC(:,114)) }, 0.00000);
    ALGEBRAIC(:,37) = 1.00000+STATES(:,3)./CONSTANTS(:,360);
    ALGEBRAIC(:,38) = 1.00000+STATES(:,2)./CONSTANTS(:,361)+STATES(:,3)./CONSTANTS(:,362)+STATES(:,4)./CONSTANTS(:,363);
    ALGEBRAIC(:,54) = ( CONSTANTS(:,316).*ALGEBRAIC(:,38).*ALGEBRAIC(:,36))./( power(STATES(:,2), 2.00000).*CONSTANTS(:,279).*ALGEBRAIC(:,37));
    ALGEBRAIC(:,39) = ( STATES(:,16).*(1.00000+STATES(:,2)./CONSTANTS(:,361)))./ALGEBRAIC(:,38);
    ALGEBRAIC(:,55) = 1.00000+ALGEBRAIC(:,39)./CONSTANTS(:,153);
    ALGEBRAIC(:,19) = ( STATES(:,1).*(1.00000+STATES(:,2)./CONSTANTS(:,297)))./ALGEBRAIC(:,1);
    ALGEBRAIC(:,24) = ( STATES(:,9).*(1.00000+STATES(:,2)./CONSTANTS(:,333)))./ALGEBRAIC(:,22);
    ALGEBRAIC(:,56) = 1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,154)+ALGEBRAIC(:,24)./CONSTANTS(:,155)+ALGEBRAIC(:,29)./CONSTANTS(:,156)+STATES(:,19)./CONSTANTS(:,157)+STATES(:,24)./CONSTANTS(:,158);
    ALGEBRAIC(:,57) = ( CONSTANTS(:,24).*( STATES(:,21).*STATES(:,23) - ( STATES(:,19).*STATES(:,16))./ALGEBRAIC(:,54)))./( CONSTANTS(:,152).*CONSTANTS(:,151).*ALGEBRAIC(:,55)+ CONSTANTS(:,150).*ALGEBRAIC(:,55).*STATES(:,21)+ CONSTANTS(:,151).*ALGEBRAIC(:,56).*STATES(:,23)+ STATES(:,23).*STATES(:,21));
    RATES(:,21) = (ALGEBRAIC(:,115) - ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,41) = 1.00000+STATES(:,2)./CONSTANTS(:,364)+STATES(:,3)./CONSTANTS(:,365);
    ALGEBRAIC(:,58) = ( CONSTANTS(:,317).*ALGEBRAIC(:,41))./ALGEBRAIC(:,38);
    ALGEBRAIC(:,59) = ( CONSTANTS(:,25).*CONSTANTS(:,160))./( CONSTANTS(:,159).*ALGEBRAIC(:,58));
    ALGEBRAIC(:,60) = ( CONSTANTS(:,25).*ALGEBRAIC(:,59).*(STATES(:,16) - STATES(:,25)./ALGEBRAIC(:,58)))./( CONSTANTS(:,159).*ALGEBRAIC(:,59)+ ALGEBRAIC(:,59).*STATES(:,16)+( CONSTANTS(:,25).*STATES(:,25))./ALGEBRAIC(:,58));
    ALGEBRAIC(:,61) = ( CONSTANTS(:,318).*ALGEBRAIC(:,46).*CONSTANTS(:,307).*CONSTANTS(:,278))./( power(STATES(:,2), 2.00000).*CONSTANTS(:,277).*ALGEBRAIC(:,41));
    ALGEBRAIC(:,62) = 1.00000+ (CONSTANTS(:,167)./ALGEBRAIC(:,24)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,166));
    ALGEBRAIC(:,116) = piecewise({ALGEBRAIC(:,114)>CONSTANTS(:,168)&STATES(:,25)>CONSTANTS(:,168), ( CONSTANTS(:,26).*(1.00000 - ( STATES(:,20).*STATES(:,26).*STATES(:,22))./( ALGEBRAIC(:,61).*ALGEBRAIC(:,114).*STATES(:,25))))./(1.00000+ power(CONSTANTS(:,162)./STATES(:,25), CONSTANTS(:,165)).*ALGEBRAIC(:,62)+ (CONSTANTS(:,161)./ALGEBRAIC(:,114)).*(1.00000+ power(CONSTANTS(:,163)./STATES(:,25), CONSTANTS(:,165)).*ALGEBRAIC(:,62)+( STATES(:,22).*ALGEBRAIC(:,62))./CONSTANTS(:,164))) }, 0.00000);
    RATES(:,25) = (ALGEBRAIC(:,60) - ALGEBRAIC(:,116))./CONSTANTS(:,4);
    ALGEBRAIC(:,63) = ( CONSTANTS(:,319).*ALGEBRAIC(:,46).*ALGEBRAIC(:,42).*CONSTANTS(:,278))./( STATES(:,2).*ALGEBRAIC(:,36).*CONSTANTS(:,307).*CONSTANTS(:,277));
    ALGEBRAIC(:,64) = 1.00000+ (CONSTANTS(:,176)./ALGEBRAIC(:,24)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,175));
    ALGEBRAIC(:,117) = piecewise({STATES(:,26)>CONSTANTS(:,177)&STATES(:,19)>CONSTANTS(:,177)&ALGEBRAIC(:,114)>CONSTANTS(:,177), (( CONSTANTS(:,27).*(1.00000 - ( STATES(:,20).*STATES(:,24).*STATES(:,22))./( ALGEBRAIC(:,63).*STATES(:,26).*STATES(:,19).*ALGEBRAIC(:,114))))./(1.00000+ (CONSTANTS(:,169)./STATES(:,26)).*ALGEBRAIC(:,64)+ (CONSTANTS(:,170)./STATES(:,19)).*(1.00000+STATES(:,24)./CONSTANTS(:,172))+ (CONSTANTS(:,171)./ALGEBRAIC(:,114)).*(1.00000+STATES(:,22)./CONSTANTS(:,173))))./(1.00000+STATES(:,22)./CONSTANTS(:,174)) }, 0.00000);
    ALGEBRAIC(:,49) = 1.00000+STATES(:,2)./CONSTANTS(:,377)+STATES(:,3)./CONSTANTS(:,378);
    ALGEBRAIC(:,50) = 1.00000+STATES(:,2)./CONSTANTS(:,379)+STATES(:,3)./CONSTANTS(:,380);
    ALGEBRAIC(:,83) = ( CONSTANTS(:,329).*ALGEBRAIC(:,49).*ALGEBRAIC(:,37))./( ALGEBRAIC(:,50).*CONSTANTS(:,307));
    ALGEBRAIC(:,84) = ( CONSTANTS(:,34).*CONSTANTS(:,234).*CONSTANTS(:,229))./( ALGEBRAIC(:,83).*CONSTANTS(:,227).*CONSTANTS(:,232));
    ALGEBRAIC(:,85) = 1.00000+STATES(:,26)./CONSTANTS(:,235);
    ALGEBRAIC(:,86) = piecewise({STATES(:,31)>CONSTANTS(:,236)&STATES(:,26)>CONSTANTS(:,236), ( CONSTANTS(:,34).*ALGEBRAIC(:,84).*STATES(:,31).*STATES(:,26))./( ALGEBRAIC(:,84).*CONSTANTS(:,232).*STATES(:,31)+ ALGEBRAIC(:,84).*CONSTANTS(:,231).*ALGEBRAIC(:,85).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*ALGEBRAIC(:,85).*STATES(:,23))./ALGEBRAIC(:,83)+( CONSTANTS(:,34).*CONSTANTS(:,233).*STATES(:,32))./ALGEBRAIC(:,83)+ ALGEBRAIC(:,84).*STATES(:,31).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*STATES(:,31).*STATES(:,23))./( ALGEBRAIC(:,83).*CONSTANTS(:,227))+( CONSTANTS(:,34).*STATES(:,23).*STATES(:,32))./ALGEBRAIC(:,83)+( ALGEBRAIC(:,84).*CONSTANTS(:,231).*STATES(:,26).*STATES(:,32))./CONSTANTS(:,230)) }, 0.00000);
    ALGEBRAIC(:,87) = piecewise({STATES(:,23)>CONSTANTS(:,236)&STATES(:,32)>CONSTANTS(:,236), ALGEBRAIC(:,86) - (( CONSTANTS(:,34).*ALGEBRAIC(:,84).*STATES(:,23).*STATES(:,32))./ALGEBRAIC(:,83))./( ALGEBRAIC(:,84).*CONSTANTS(:,232).*STATES(:,31)+ ALGEBRAIC(:,84).*CONSTANTS(:,231).*ALGEBRAIC(:,85).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*ALGEBRAIC(:,85).*STATES(:,23))./ALGEBRAIC(:,83)+( CONSTANTS(:,34).*CONSTANTS(:,233).*STATES(:,32))./ALGEBRAIC(:,83)+ ALGEBRAIC(:,84).*STATES(:,31).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*STATES(:,31).*STATES(:,23))./( ALGEBRAIC(:,83).*CONSTANTS(:,227))+( CONSTANTS(:,34).*STATES(:,23).*STATES(:,32))./ALGEBRAIC(:,83)+( ALGEBRAIC(:,84).*CONSTANTS(:,231).*STATES(:,26).*STATES(:,32))./CONSTANTS(:,230)) }, ALGEBRAIC(:,86));
    RATES(:,26) = (((ALGEBRAIC(:,116) - ALGEBRAIC(:,117)) - ALGEBRAIC(:,87))+ALGEBRAIC(:,88))./CONSTANTS(:,4);
    RATES(:,24) = (ALGEBRAIC(:,117) - ALGEBRAIC(:,67))./CONSTANTS(:,4);
    RATES(:,19) = (( - ALGEBRAIC(:,115) - ALGEBRAIC(:,117))+ALGEBRAIC(:,67)+ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,45) = 1.00000+STATES(:,2)./CONSTANTS(:,371)+STATES(:,3)./CONSTANTS(:,372)+STATES(:,4)./CONSTANTS(:,373);
    ALGEBRAIC(:,75) = ( CONSTANTS(:,322).*ALGEBRAIC(:,37).*CONSTANTS(:,278))./( STATES(:,2).*ALGEBRAIC(:,45).*CONSTANTS(:,277));
    ALGEBRAIC(:,76) = ( CONSTANTS(:,32).*CONSTANTS(:,209).*CONSTANTS(:,212))./( ALGEBRAIC(:,75).*CONSTANTS(:,206).*CONSTANTS(:,211));
    ALGEBRAIC(:,77) = 1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,214)+ALGEBRAIC(:,24)./CONSTANTS(:,215)+ALGEBRAIC(:,29)./CONSTANTS(:,216);
    ALGEBRAIC(:,118) = ( CONSTANTS(:,32).*ALGEBRAIC(:,76).*( ALGEBRAIC(:,114).*STATES(:,30) - ( STATES(:,23).*STATES(:,22))./ALGEBRAIC(:,75)))./( ALGEBRAIC(:,76).*CONSTANTS(:,206).*CONSTANTS(:,211).*ALGEBRAIC(:,77)+ ALGEBRAIC(:,76).*CONSTANTS(:,211).*ALGEBRAIC(:,114)+ ALGEBRAIC(:,76).*CONSTANTS(:,210).*ALGEBRAIC(:,77).*STATES(:,30)+( CONSTANTS(:,32).*CONSTANTS(:,213).*ALGEBRAIC(:,77).*STATES(:,23))./ALGEBRAIC(:,75)+( CONSTANTS(:,32).*CONSTANTS(:,212).*STATES(:,22))./ALGEBRAIC(:,75)+ ALGEBRAIC(:,76).*ALGEBRAIC(:,114).*STATES(:,30)+( CONSTANTS(:,32).*CONSTANTS(:,213).*ALGEBRAIC(:,114).*STATES(:,23))./( ALGEBRAIC(:,75).*CONSTANTS(:,206))+( CONSTANTS(:,32).*STATES(:,23).*STATES(:,22))./ALGEBRAIC(:,75)+( ALGEBRAIC(:,76).*CONSTANTS(:,210).*STATES(:,30).*STATES(:,22))./CONSTANTS(:,209)+( ALGEBRAIC(:,76).*ALGEBRAIC(:,114).*STATES(:,30).*STATES(:,23))./CONSTANTS(:,208)+( CONSTANTS(:,32).*STATES(:,30).*STATES(:,23).*STATES(:,22))./( CONSTANTS(:,207).*ALGEBRAIC(:,75)));
    RATES(:,23) = ( - ALGEBRAIC(:,57)+ALGEBRAIC(:,118)+ALGEBRAIC(:,87))./CONSTANTS(:,4);
    ALGEBRAIC(:,119) = CONSTANTS(:,16) - STATES(:,28);
    ALGEBRAIC(:,44) = 1.00000+STATES(:,2)./CONSTANTS(:,370);
    ALGEBRAIC(:,68) = ( CONSTANTS(:,335).*ALGEBRAIC(:,44).*CONSTANTS(:,305))./( ALGEBRAIC(:,43).*CONSTANTS(:,321));
    ALGEBRAIC(:,69) = ( CONSTANTS(:,30).*CONSTANTS(:,194).*CONSTANTS(:,191))./( ALGEBRAIC(:,68).*CONSTANTS(:,190).*CONSTANTS(:,193));
    ALGEBRAIC(:,70) = (1.00000+STATES(:,23)./CONSTANTS(:,196)+STATES(:,27)./CONSTANTS(:,197)+STATES(:,29)./CONSTANTS(:,198))./(1.00000+STATES(:,27)./CONSTANTS(:,197)+STATES(:,29)./CONSTANTS(:,198));
    ALGEBRAIC(:,120) = ( CONSTANTS(:,30).*ALGEBRAIC(:,69).*( STATES(:,27).*ALGEBRAIC(:,119) - ( STATES(:,28).*STATES(:,29))./ALGEBRAIC(:,68)))./( ALGEBRAIC(:,69).*CONSTANTS(:,190).*CONSTANTS(:,193).*ALGEBRAIC(:,70)+ ALGEBRAIC(:,69).*CONSTANTS(:,193).*STATES(:,27)+ ALGEBRAIC(:,69).*CONSTANTS(:,192).*ALGEBRAIC(:,70).*ALGEBRAIC(:,119)+( CONSTANTS(:,30).*CONSTANTS(:,195).*ALGEBRAIC(:,70).*STATES(:,28))./ALGEBRAIC(:,68)+( CONSTANTS(:,30).*CONSTANTS(:,194).*STATES(:,29))./ALGEBRAIC(:,68)+ ALGEBRAIC(:,69).*STATES(:,27).*ALGEBRAIC(:,119)+( CONSTANTS(:,30).*CONSTANTS(:,195).*STATES(:,27).*STATES(:,28))./( ALGEBRAIC(:,68).*CONSTANTS(:,190))+( ALGEBRAIC(:,69).*CONSTANTS(:,192).*ALGEBRAIC(:,119).*STATES(:,29))./CONSTANTS(:,191)+( CONSTANTS(:,30).*STATES(:,28).*STATES(:,29))./ALGEBRAIC(:,68));
    ALGEBRAIC(:,71) = ( CONSTANTS(:,343).*ALGEBRAIC(:,45))./ALGEBRAIC(:,44);
    ALGEBRAIC(:,72) = ( CONSTANTS(:,31).*CONSTANTS(:,200))./( ALGEBRAIC(:,71).*CONSTANTS(:,199));
    ALGEBRAIC(:,31) = ( STATES(:,12).*(1.00000+STATES(:,2)./CONSTANTS(:,314)))./ALGEBRAIC(:,30);
    ALGEBRAIC(:,33) = ( STATES(:,13).*(1.00000+STATES(:,2)./CONSTANTS(:,338)))./ALGEBRAIC(:,32);
    ALGEBRAIC(:,73) = 1.00000+STATES(:,16)./CONSTANTS(:,201)+ALGEBRAIC(:,19)./CONSTANTS(:,202)+ALGEBRAIC(:,24)./CONSTANTS(:,203)+ALGEBRAIC(:,31)./CONSTANTS(:,204)+ALGEBRAIC(:,33)./CONSTANTS(:,205);
    ALGEBRAIC(:,74) = ( ALGEBRAIC(:,72).*CONSTANTS(:,31).*(STATES(:,29) - STATES(:,30)./ALGEBRAIC(:,71)))./( CONSTANTS(:,199).*ALGEBRAIC(:,72).*ALGEBRAIC(:,73)+ ALGEBRAIC(:,72).*STATES(:,29)+( CONSTANTS(:,31).*STATES(:,30))./ALGEBRAIC(:,71));
    RATES(:,29) = (ALGEBRAIC(:,120) - ALGEBRAIC(:,74))./CONSTANTS(:,4);
    ALGEBRAIC(:,133) = STATES(:,6);
    ALGEBRAIC(:,101) =  1.00000.*exp( - (CONSTANTS(:,288)+ 4.00000.*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275));
    ALGEBRAIC(:,139) = ( ALGEBRAIC(:,101).*power(STATES(:,2), 5.00000))./power(ALGEBRAIC(:,133), 4.00000);
    ALGEBRAIC(:,140) =  CONSTANTS(:,42).*( ALGEBRAIC(:,139).*STATES(:,22).*ALGEBRAIC(:,119) -  ALGEBRAIC(:,114).*STATES(:,28));
    RATES(:,22) = ((ALGEBRAIC(:,115)+ALGEBRAIC(:,116)+ALGEBRAIC(:,117)+ALGEBRAIC(:,118)) - ALGEBRAIC(:,140))./CONSTANTS(:,4);
    ALGEBRAIC(:,103) =  1.00000.*exp( - (CONSTANTS(:,289)+ 2.00000.*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275));
    ALGEBRAIC(:,105) = ( ALGEBRAIC(:,103).*power(STATES(:,2), 2.00000))./power(STATES(:,2), 4.00000);
    ALGEBRAIC(:,107) = piecewise({STATES(:,28)>CONSTANTS(:,248), STATES(:,28) }, CONSTANTS(:,248));
    ALGEBRAIC(:,121) = piecewise({ALGEBRAIC(:,119)>CONSTANTS(:,248), ALGEBRAIC(:,119) }, CONSTANTS(:,248));
    ALGEBRAIC(:,122) = CONSTANTS(:,17) - STATES(:,56);
    ALGEBRAIC(:,123) =  (( CONSTANTS(:,43).*(1.00000+STATES(:,14)./CONSTANTS(:,44)))./(1.00000+STATES(:,14)./CONSTANTS(:,45))).*( power(ALGEBRAIC(:,105), 1.0 ./ 2).*ALGEBRAIC(:,122).*power(ALGEBRAIC(:,107), 1.0 ./ 2) -  STATES(:,56).*power(ALGEBRAIC(:,121), 1.0 ./ 2));
    RATES(:,28) = ((ALGEBRAIC(:,120)+ALGEBRAIC(:,140)) - ALGEBRAIC(:,123))./CONSTANTS(:,4);
    ALGEBRAIC(:,109) = exp( - (CONSTANTS(:,290)+ 4.00000.*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275));
    ALGEBRAIC(:,141) = ( ALGEBRAIC(:,109).*power(STATES(:,2), 4.00000))./( power(ALGEBRAIC(:,133), 2.00000).*power(1.00000, 2.00000));
    ALGEBRAIC(:,111) = piecewise({STATES(:,57)>CONSTANTS(:,249), STATES(:,57) }, CONSTANTS(:,249));
    ALGEBRAIC(:,142) =  CONSTANTS(:,46).*( (ALGEBRAIC(:,111)./(ALGEBRAIC(:,111)+CONSTANTS(:,14))).*exp(( CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)).*(STATES(:,56)./CONSTANTS(:,17))).*( power(ALGEBRAIC(:,141), 1.00000./2.00000).*STATES(:,56).*power(ALGEBRAIC(:,111)./1.00000, 1.00000./4.00000) - ALGEBRAIC(:,122));
    RATES(:,56) = ( 2.00000.*ALGEBRAIC(:,123) -  2.00000.*ALGEBRAIC(:,142))./CONSTANTS(:,5);
    ALGEBRAIC(:,150) = STATES(:,7);
    ALGEBRAIC(:,153) = 1.00000+ALGEBRAIC(:,150)./CONSTANTS(:,375);
    ALGEBRAIC(:,156) =  CONSTANTS(:,35).*( (STATES(:,33)./ALGEBRAIC(:,153)).*ALGEBRAIC(:,133) -  (STATES(:,18)./ALGEBRAIC(:,47)).*STATES(:,2));
    RATES(:,18) = ( - ALGEBRAIC(:,115)+ALGEBRAIC(:,156))./CONSTANTS(:,4);
    RATES(:,33) = ( - ALGEBRAIC(:,156)+ALGEBRAIC(:,89))./CONSTANTS(:,5);
    ALGEBRAIC(:,154) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,377)+ALGEBRAIC(:,150)./CONSTANTS(:,378);
    ALGEBRAIC(:,157) =  CONSTANTS(:,36).*( (STATES(:,34)./ALGEBRAIC(:,154)).*ALGEBRAIC(:,133) -  (STATES(:,32)./ALGEBRAIC(:,49)).*STATES(:,2));
    ALGEBRAIC(:,155) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,379)+ALGEBRAIC(:,150)./CONSTANTS(:,380);
    ALGEBRAIC(:,158) = ( exp((  - CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)).*ALGEBRAIC(:,50).*ALGEBRAIC(:,154))./( ALGEBRAIC(:,155).*ALGEBRAIC(:,49));
    ALGEBRAIC(:,159) = ( CONSTANTS(:,41).*( ALGEBRAIC(:,158).*STATES(:,37).*STATES(:,32).*STATES(:,2) -  STATES(:,31).*STATES(:,34).*ALGEBRAIC(:,133)))./( ALGEBRAIC(:,158).*CONSTANTS(:,241).*CONSTANTS(:,244).*CONSTANTS(:,286).*( 2.00000.*CONSTANTS(:,245)+( CONSTANTS(:,245).*STATES(:,37))./CONSTANTS(:,241)+( STATES(:,37).*STATES(:,32).*STATES(:,2))./( CONSTANTS(:,241).*CONSTANTS(:,244).*CONSTANTS(:,286))+( CONSTANTS(:,245).*STATES(:,31).*ALGEBRAIC(:,133))./( CONSTANTS(:,242).*CONSTANTS(:,286))+( STATES(:,31).*STATES(:,34).*ALGEBRAIC(:,133))./( CONSTANTS(:,242).*CONSTANTS(:,243).*CONSTANTS(:,286))+( CONSTANTS(:,245).*STATES(:,31))./CONSTANTS(:,242)+( CONSTANTS(:,245).*STATES(:,37).*STATES(:,2))./( CONSTANTS(:,241).*CONSTANTS(:,286))+( CONSTANTS(:,245).*STATES(:,2))./CONSTANTS(:,286)+( CONSTANTS(:,245).*STATES(:,34).*ALGEBRAIC(:,133))./( CONSTANTS(:,243).*CONSTANTS(:,286))+( CONSTANTS(:,245).*ALGEBRAIC(:,133))./CONSTANTS(:,286)+( CONSTANTS(:,245).*STATES(:,32).*STATES(:,2))./( CONSTANTS(:,244).*CONSTANTS(:,286))));
    RATES(:,32) = ((ALGEBRAIC(:,87)+ALGEBRAIC(:,157)) - ALGEBRAIC(:,159))./CONSTANTS(:,4);
    RATES(:,31) = ( - ALGEBRAIC(:,87)+ALGEBRAIC(:,159))./CONSTANTS(:,4);
    RATES(:,34) = ( - ALGEBRAIC(:,157)+ALGEBRAIC(:,159)+ALGEBRAIC(:,94))./CONSTANTS(:,5);
    RATES(:,37) = ( - ALGEBRAIC(:,159)+ALGEBRAIC(:,95))./CONSTANTS(:,5);
    ALGEBRAIC(:,40) = (( STATES(:,16).*STATES(:,2))./CONSTANTS(:,361))./ALGEBRAIC(:,38);
    ALGEBRAIC(:,163) = STATES(:,8);
    ALGEBRAIC(:,171) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,361)+ALGEBRAIC(:,150)./CONSTANTS(:,362)+ALGEBRAIC(:,163)./CONSTANTS(:,363);
    ALGEBRAIC(:,173) = (( STATES(:,17).*ALGEBRAIC(:,133))./CONSTANTS(:,361))./ALGEBRAIC(:,171);
    ALGEBRAIC(:,175) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,371)+ALGEBRAIC(:,150)./CONSTANTS(:,372)+ALGEBRAIC(:,163)./CONSTANTS(:,373);
    ALGEBRAIC(:,176) =  CONSTANTS(:,37).*(( ALGEBRAIC(:,173).*STATES(:,30))./ALGEBRAIC(:,45) - ( ALGEBRAIC(:,40).*STATES(:,35))./ALGEBRAIC(:,175));
    RATES(:,16) = ( - ALGEBRAIC(:,60)+ALGEBRAIC(:,176)+ALGEBRAIC(:,57))./CONSTANTS(:,4);
    RATES(:,17) = ( - ALGEBRAIC(:,176)+ALGEBRAIC(:,91))./CONSTANTS(:,5);
    ALGEBRAIC(:,174) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,367)+ALGEBRAIC(:,150)./CONSTANTS(:,368)+ALGEBRAIC(:,163)./CONSTANTS(:,369);
    ALGEBRAIC(:,178) =  CONSTANTS(:,39).*(( (STATES(:,39)./ALGEBRAIC(:,174)).*STATES(:,30))./ALGEBRAIC(:,45) - ( (STATES(:,27)./ALGEBRAIC(:,43)).*STATES(:,35))./ALGEBRAIC(:,175));
    RATES(:,27) = ((ALGEBRAIC(:,67) - ALGEBRAIC(:,120))+ALGEBRAIC(:,178))./CONSTANTS(:,4);
    ALGEBRAIC(:,169) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,352)+ALGEBRAIC(:,150)./CONSTANTS(:,354)+ALGEBRAIC(:,163)./CONSTANTS(:,356);
    ALGEBRAIC(:,177) =  CONSTANTS(:,40).*(( (STATES(:,35)./ALGEBRAIC(:,175)).*STATES(:,14))./ALGEBRAIC(:,34) - ( (STATES(:,30)./ALGEBRAIC(:,45)).*STATES(:,15))./ALGEBRAIC(:,169));
    RATES(:,30) = (((((ALGEBRAIC(:,74) - ALGEBRAIC(:,118))+ALGEBRAIC(:,177)) - ALGEBRAIC(:,88)) - ALGEBRAIC(:,176)) - ALGEBRAIC(:,178))./CONSTANTS(:,4);
    RATES(:,39) = (ALGEBRAIC(:,93) - ALGEBRAIC(:,178))./CONSTANTS(:,5);
    RATES(:,35) = ( - ALGEBRAIC(:,177)+ALGEBRAIC(:,92)+ALGEBRAIC(:,88)+ALGEBRAIC(:,176)+ALGEBRAIC(:,178))./CONSTANTS(:,5);
    ALGEBRAIC(:,113) = exp( - (CONSTANTS(:,291) -  CONSTANTS(:,8).*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275))./1.00000;
    ALGEBRAIC(:,143) = ( (( ALGEBRAIC(:,113).*power(ALGEBRAIC(:,133), CONSTANTS(:,8)))./power(STATES(:,2), CONSTANTS(:,8) - 1.00000)).*ALGEBRAIC(:,1))./( ALGEBRAIC(:,22).*ALGEBRAIC(:,34));
    ALGEBRAIC(:,144) =  CONSTANTS(:,47).*( ALGEBRAIC(:,143).*STATES(:,9).*STATES(:,14) - STATES(:,1));
    ALGEBRAIC(:,164) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,297)+ALGEBRAIC(:,150)./CONSTANTS(:,312)+ALGEBRAIC(:,163)./CONSTANTS(:,326);
    ALGEBRAIC(:,166) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,333)+ALGEBRAIC(:,150)./CONSTANTS(:,337)+ALGEBRAIC(:,163)./CONSTANTS(:,341);
    ALGEBRAIC(:,179) = piecewise({STATES(:,51)>CONSTANTS(:,250)&STATES(:,1)>CONSTANTS(:,250),  CONSTANTS(:,48).*((STATES(:,51)./ALGEBRAIC(:,166))./(STATES(:,51)./ALGEBRAIC(:,166)+ (STATES(:,52)./ALGEBRAIC(:,164)).*exp(( CONSTANTS(:,13).*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275))) - (STATES(:,9)./ALGEBRAIC(:,22))./(STATES(:,9)./ALGEBRAIC(:,22)+ (STATES(:,1)./ALGEBRAIC(:,1)).*exp(( CONSTANTS(:,13).*(CONSTANTS(:,3) - 1.00000).*STATES(:,38))./CONSTANTS(:,275)))).*((STATES(:,51)./ALGEBRAIC(:,166))./(STATES(:,51)./ALGEBRAIC(:,166)+CONSTANTS(:,12))) }, 0.00000);
    ALGEBRAIC(:,145) = piecewise({abs(STATES(:,38))>CONSTANTS(:,251), ( CONSTANTS(:,52).*STATES(:,38).*( ALGEBRAIC(:,133).*exp(( CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)) - STATES(:,2)))./(exp(( CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)) - 1.00000) }, ( CONSTANTS(:,52).*CONSTANTS(:,275).*(ALGEBRAIC(:,133) - STATES(:,2)))./CONSTANTS(:,3));
    RATES(:,38) = ((((( 4.00000.*ALGEBRAIC(:,140)+ 2.00000.*ALGEBRAIC(:,123)+ 4.00000.*ALGEBRAIC(:,142)) -  CONSTANTS(:,8).*ALGEBRAIC(:,144)) - ALGEBRAIC(:,179)) - ALGEBRAIC(:,145))+ALGEBRAIC(:,159))./CONSTANTS(:,22);
    RATES(:,1) = ((ALGEBRAIC(:,82)+ALGEBRAIC(:,144)) - ALGEBRAIC(:,179))./CONSTANTS(:,4);
    RATES(:,9) = (( - ALGEBRAIC(:,82) - ALGEBRAIC(:,144))+ALGEBRAIC(:,179))./CONSTANTS(:,4);
    ALGEBRAIC(:,35) = (( STATES(:,14).*STATES(:,2))./CONSTANTS(:,352))./ALGEBRAIC(:,34);
    ALGEBRAIC(:,170) = (( STATES(:,15).*ALGEBRAIC(:,133))./CONSTANTS(:,352))./ALGEBRAIC(:,169);
    ALGEBRAIC(:,180) = ( CONSTANTS(:,49).*( ALGEBRAIC(:,170).*ALGEBRAIC(:,133) -  ALGEBRAIC(:,35).*STATES(:,2)))./( CONSTANTS(:,50).*(1.00000+ALGEBRAIC(:,170)./CONSTANTS(:,50)).*(1.00000+ALGEBRAIC(:,35)./CONSTANTS(:,50)));
    RATES(:,14) = ((( - ALGEBRAIC(:,67) - ALGEBRAIC(:,144))+ALGEBRAIC(:,180)) - ALGEBRAIC(:,177))./CONSTANTS(:,4);
    RATES(:,15) = ( - ALGEBRAIC(:,180)+ALGEBRAIC(:,99)+ALGEBRAIC(:,177))./CONSTANTS(:,5);
    ALGEBRAIC(:,168) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,345)+ALGEBRAIC(:,150)./CONSTANTS(:,348)+ALGEBRAIC(:,163)./CONSTANTS(:,350);
    ALGEBRAIC(:,181) = ( CONSTANTS(:,382).*ALGEBRAIC(:,164).*ALGEBRAIC(:,168))./( ALGEBRAIC(:,166).*ALGEBRAIC(:,166));
    ALGEBRAIC(:,183) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000, 0.00000 },  CONSTANTS(:,262).*( ALGEBRAIC(:,181).*STATES(:,51).*STATES(:,51) -  STATES(:,54).*STATES(:,52)));
    RATES(:,52) = (ALGEBRAIC(:,97)+ALGEBRAIC(:,179)+ALGEBRAIC(:,183))./CONSTANTS(:,5);
    RATES(:,51) = ((ALGEBRAIC(:,96) - ALGEBRAIC(:,179)) -  2.00000.*ALGEBRAIC(:,183))./CONSTANTS(:,5);
    RATES(:,54) = (ALGEBRAIC(:,98)+ALGEBRAIC(:,183))./CONSTANTS(:,5);
    ALGEBRAIC(:,190) =  - (( STATES(:,2).*STATES(:,2))./( CONSTANTS(:,135).*CONSTANTS(:,136).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,2).*STATES(:,1))./( CONSTANTS(:,297).*CONSTANTS(:,312).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,2).*STATES(:,9))./( CONSTANTS(:,333).*CONSTANTS(:,337).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,2).*STATES(:,11))./( CONSTANTS(:,345).*CONSTANTS(:,348).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,2).*STATES(:,12))./( CONSTANTS(:,297).*CONSTANTS(:,312).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,2).*STATES(:,13))./( CONSTANTS(:,333).*CONSTANTS(:,337).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,2).*STATES(:,14))./( CONSTANTS(:,352).*CONSTANTS(:,354).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,2).*STATES(:,22))./( CONSTANTS(:,111).*CONSTANTS(:,112).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,114))./( CONSTANTS(:,108).*CONSTANTS(:,109).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,2).*STATES(:,28))./( CONSTANTS(:,114).*CONSTANTS(:,115).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,119))./( CONSTANTS(:,117).*CONSTANTS(:,118).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,2).*STATES(:,18))./( CONSTANTS(:,99).*CONSTANTS(:,375).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,2).*STATES(:,23))./( CONSTANTS(:,90).*CONSTANTS(:,360).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,2).*STATES(:,21))./( CONSTANTS(:,120).*CONSTANTS(:,121).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,2).*STATES(:,16))./( CONSTANTS(:,361).*CONSTANTS(:,362).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,2).*STATES(:,25))./( CONSTANTS(:,364).*CONSTANTS(:,365).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,2).*STATES(:,26))./( CONSTANTS(:,123).*CONSTANTS(:,124).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,2).*STATES(:,24))./( CONSTANTS(:,366).*CONSTANTS(:,93).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,2).*STATES(:,19))./( CONSTANTS(:,358).*CONSTANTS(:,88).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,2).*STATES(:,27))./( CONSTANTS(:,367).*CONSTANTS(:,368).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,2).*STATES(:,29))./( CONSTANTS(:,370).*CONSTANTS(:,95).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,2).*STATES(:,30))./( CONSTANTS(:,371).*CONSTANTS(:,372).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,2).*STATES(:,32))./( CONSTANTS(:,377).*CONSTANTS(:,378).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,2).*STATES(:,31))./( CONSTANTS(:,379).*CONSTANTS(:,380).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,2).*STATES(:,4))./( CONSTANTS(:,141).*CONSTANTS(:,142).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,2).*STATES(:,3))./( CONSTANTS(:,138).*CONSTANTS(:,139).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,2).*STATES(:,57))./( CONSTANTS(:,126).*CONSTANTS(:,127).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,2).*CONSTANTS(:,292))./( CONSTANTS(:,129).*CONSTANTS(:,130).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,2).*CONSTANTS(:,311))./( CONSTANTS(:,132).*CONSTANTS(:,133).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,2).*STATES(:,20))./( CONSTANTS(:,374).*CONSTANTS(:,97).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,191) =  - (( STATES(:,2).*STATES(:,2))./( CONSTANTS(:,135).*CONSTANTS(:,137).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,2).*STATES(:,1))./( CONSTANTS(:,297).*CONSTANTS(:,326).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,2).*STATES(:,9))./( CONSTANTS(:,333).*CONSTANTS(:,341).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,2).*STATES(:,11))./( CONSTANTS(:,345).*CONSTANTS(:,350).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,2).*STATES(:,12))./( CONSTANTS(:,297).*CONSTANTS(:,326).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,2).*STATES(:,13))./( CONSTANTS(:,333).*CONSTANTS(:,341).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,2).*STATES(:,14))./( CONSTANTS(:,352).*CONSTANTS(:,356).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,2).*STATES(:,22))./( CONSTANTS(:,111).*CONSTANTS(:,113).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,114))./( CONSTANTS(:,108).*CONSTANTS(:,110).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,2).*STATES(:,28))./( CONSTANTS(:,114).*CONSTANTS(:,116).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,119))./( CONSTANTS(:,117).*CONSTANTS(:,119).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,2).*STATES(:,18))./( CONSTANTS(:,99).*CONSTANTS(:,100).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,2).*STATES(:,23))./( CONSTANTS(:,90).*CONSTANTS(:,91).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,2).*STATES(:,21))./( CONSTANTS(:,120).*CONSTANTS(:,122).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,2).*STATES(:,16))./( CONSTANTS(:,361).*CONSTANTS(:,363).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,2).*STATES(:,25))./( CONSTANTS(:,364).*CONSTANTS(:,92).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,2).*STATES(:,26))./( CONSTANTS(:,123).*CONSTANTS(:,125).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,2).*STATES(:,24))./( CONSTANTS(:,366).*CONSTANTS(:,94).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,2).*STATES(:,19))./( CONSTANTS(:,358).*CONSTANTS(:,89).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,2).*STATES(:,27))./( CONSTANTS(:,367).*CONSTANTS(:,369).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,2).*STATES(:,29))./( CONSTANTS(:,370).*CONSTANTS(:,96).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,2).*STATES(:,30))./( CONSTANTS(:,371).*CONSTANTS(:,373).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,2).*STATES(:,32))./( CONSTANTS(:,377).*CONSTANTS(:,103).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,2).*STATES(:,31))./( CONSTANTS(:,379).*CONSTANTS(:,104).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,2).*STATES(:,4))./( CONSTANTS(:,141).*CONSTANTS(:,143).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,2).*STATES(:,3))./( CONSTANTS(:,138).*CONSTANTS(:,140).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,2).*STATES(:,57))./( CONSTANTS(:,126).*CONSTANTS(:,128).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,2).*CONSTANTS(:,292))./( CONSTANTS(:,129).*CONSTANTS(:,131).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,2).*CONSTANTS(:,311))./( CONSTANTS(:,132).*CONSTANTS(:,134).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,2).*STATES(:,20))./( CONSTANTS(:,374).*CONSTANTS(:,98).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,194) =  - (( STATES(:,3).*STATES(:,2))./( CONSTANTS(:,136).*CONSTANTS(:,137).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,3).*STATES(:,1))./( CONSTANTS(:,312).*CONSTANTS(:,326).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,3).*STATES(:,9))./( CONSTANTS(:,337).*CONSTANTS(:,341).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,3).*STATES(:,11))./( CONSTANTS(:,348).*CONSTANTS(:,350).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,3).*STATES(:,12))./( CONSTANTS(:,312).*CONSTANTS(:,326).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,3).*STATES(:,13))./( CONSTANTS(:,337).*CONSTANTS(:,341).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,3).*STATES(:,14))./( CONSTANTS(:,354).*CONSTANTS(:,356).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,3).*STATES(:,22))./( CONSTANTS(:,112).*CONSTANTS(:,113).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,114))./( CONSTANTS(:,109).*CONSTANTS(:,110).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,3).*STATES(:,28))./( CONSTANTS(:,115).*CONSTANTS(:,116).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,119))./( CONSTANTS(:,118).*CONSTANTS(:,119).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,3).*STATES(:,18))./( CONSTANTS(:,375).*CONSTANTS(:,100).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,3).*STATES(:,23))./( CONSTANTS(:,360).*CONSTANTS(:,91).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,3).*STATES(:,21))./( CONSTANTS(:,121).*CONSTANTS(:,122).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,3).*STATES(:,16))./( CONSTANTS(:,362).*CONSTANTS(:,363).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,3).*STATES(:,25))./( CONSTANTS(:,365).*CONSTANTS(:,92).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,3).*STATES(:,26))./( CONSTANTS(:,124).*CONSTANTS(:,125).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,3).*STATES(:,24))./( CONSTANTS(:,93).*CONSTANTS(:,94).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,3).*STATES(:,19))./( CONSTANTS(:,88).*CONSTANTS(:,89).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,3).*STATES(:,27))./( CONSTANTS(:,368).*CONSTANTS(:,369).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,3).*STATES(:,29))./( CONSTANTS(:,95).*CONSTANTS(:,96).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,3).*STATES(:,30))./( CONSTANTS(:,372).*CONSTANTS(:,373).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,3).*STATES(:,32))./( CONSTANTS(:,378).*CONSTANTS(:,103).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,3).*STATES(:,31))./( CONSTANTS(:,380).*CONSTANTS(:,104).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,3).*STATES(:,4))./( CONSTANTS(:,142).*CONSTANTS(:,143).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,3).*STATES(:,3))./( CONSTANTS(:,139).*CONSTANTS(:,140).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,3).*STATES(:,57))./( CONSTANTS(:,127).*CONSTANTS(:,128).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,3).*CONSTANTS(:,292))./( CONSTANTS(:,130).*CONSTANTS(:,131).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,3).*CONSTANTS(:,311))./( CONSTANTS(:,133).*CONSTANTS(:,134).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,3).*STATES(:,20))./( CONSTANTS(:,97).*CONSTANTS(:,98).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,197) =  - (( STATES(:,4).*STATES(:,2))./( CONSTANTS(:,137).*CONSTANTS(:,136).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,4).*STATES(:,1))./( CONSTANTS(:,326).*CONSTANTS(:,312).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,4).*STATES(:,9))./( CONSTANTS(:,341).*CONSTANTS(:,337).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,4).*STATES(:,11))./( CONSTANTS(:,350).*CONSTANTS(:,348).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,4).*STATES(:,12))./( CONSTANTS(:,326).*CONSTANTS(:,312).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,4).*STATES(:,13))./( CONSTANTS(:,341).*CONSTANTS(:,337).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,4).*STATES(:,14))./( CONSTANTS(:,356).*CONSTANTS(:,354).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,4).*STATES(:,22))./( CONSTANTS(:,113).*CONSTANTS(:,112).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,114))./( CONSTANTS(:,110).*CONSTANTS(:,109).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,4).*STATES(:,28))./( CONSTANTS(:,116).*CONSTANTS(:,115).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,119))./( CONSTANTS(:,119).*CONSTANTS(:,118).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,4).*STATES(:,18))./( CONSTANTS(:,100).*CONSTANTS(:,375).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,4).*STATES(:,23))./( CONSTANTS(:,91).*CONSTANTS(:,360).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,4).*STATES(:,21))./( CONSTANTS(:,122).*CONSTANTS(:,121).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,4).*STATES(:,16))./( CONSTANTS(:,363).*CONSTANTS(:,362).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,4).*STATES(:,25))./( CONSTANTS(:,92).*CONSTANTS(:,365).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,4).*STATES(:,26))./( CONSTANTS(:,125).*CONSTANTS(:,124).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,4).*STATES(:,24))./( CONSTANTS(:,94).*CONSTANTS(:,93).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,4).*STATES(:,19))./( CONSTANTS(:,89).*CONSTANTS(:,88).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,4).*STATES(:,27))./( CONSTANTS(:,369).*CONSTANTS(:,368).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,4).*STATES(:,29))./( CONSTANTS(:,96).*CONSTANTS(:,95).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,4).*STATES(:,30))./( CONSTANTS(:,373).*CONSTANTS(:,372).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,4).*STATES(:,32))./( CONSTANTS(:,103).*CONSTANTS(:,378).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,4).*STATES(:,31))./( CONSTANTS(:,104).*CONSTANTS(:,380).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,4).*STATES(:,4))./( CONSTANTS(:,143).*CONSTANTS(:,142).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,4).*STATES(:,3))./( CONSTANTS(:,140).*CONSTANTS(:,139).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,4).*STATES(:,57))./( CONSTANTS(:,128).*CONSTANTS(:,127).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,4).*CONSTANTS(:,292))./( CONSTANTS(:,131).*CONSTANTS(:,130).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,4).*CONSTANTS(:,311))./( CONSTANTS(:,134).*CONSTANTS(:,133).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,4).*STATES(:,20))./( CONSTANTS(:,98).*CONSTANTS(:,97).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,198) = ( (1.00000+STATES(:,2)./CONSTANTS(:,135)+STATES(:,3)./CONSTANTS(:,136)).*STATES(:,2))./( CONSTANTS(:,137).*power(CONSTANTS(:,283), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,3)./CONSTANTS(:,312)).*STATES(:,1))./( CONSTANTS(:,326).*power(ALGEBRAIC(:,1), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,3)./CONSTANTS(:,337)).*STATES(:,9))./( CONSTANTS(:,341).*power(ALGEBRAIC(:,22), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,345)+STATES(:,3)./CONSTANTS(:,348)).*STATES(:,11))./( CONSTANTS(:,350).*power(ALGEBRAIC(:,27), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,3)./CONSTANTS(:,312)).*STATES(:,12))./( CONSTANTS(:,326).*power(ALGEBRAIC(:,30), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,3)./CONSTANTS(:,337)).*STATES(:,13))./( CONSTANTS(:,341).*power(ALGEBRAIC(:,32), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,352)+STATES(:,3)./CONSTANTS(:,354)).*STATES(:,14))./( CONSTANTS(:,356).*power(ALGEBRAIC(:,34), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,111)+STATES(:,3)./CONSTANTS(:,112)).*STATES(:,22))./( CONSTANTS(:,113).*power(CONSTANTS(:,278), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,108)+STATES(:,3)./CONSTANTS(:,109)).*ALGEBRAIC(:,114))./( CONSTANTS(:,110).*power(CONSTANTS(:,277), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,114)+STATES(:,3)./CONSTANTS(:,115)).*STATES(:,28))./( CONSTANTS(:,116).*power(CONSTANTS(:,305), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,117)+STATES(:,3)./CONSTANTS(:,118)).*ALGEBRAIC(:,119))./( CONSTANTS(:,119).*power(CONSTANTS(:,321), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,99)+STATES(:,3)./CONSTANTS(:,375)).*STATES(:,18))./( CONSTANTS(:,100).*power(ALGEBRAIC(:,47), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,90)+STATES(:,3)./CONSTANTS(:,360)).*STATES(:,23))./( CONSTANTS(:,91).*power(ALGEBRAIC(:,37), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,120)+STATES(:,3)./CONSTANTS(:,121)).*STATES(:,21))./( CONSTANTS(:,122).*power(CONSTANTS(:,279), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,361)+STATES(:,3)./CONSTANTS(:,362)).*STATES(:,16))./( CONSTANTS(:,363).*power(ALGEBRAIC(:,38), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,364)+STATES(:,3)./CONSTANTS(:,365)).*STATES(:,25))./( CONSTANTS(:,92).*power(ALGEBRAIC(:,41), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,123)+STATES(:,3)./CONSTANTS(:,124)).*STATES(:,26))./( CONSTANTS(:,125).*power(CONSTANTS(:,307), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,366)+STATES(:,3)./CONSTANTS(:,93)).*STATES(:,24))./( CONSTANTS(:,94).*power(ALGEBRAIC(:,42), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,358)+STATES(:,3)./CONSTANTS(:,88)).*STATES(:,19))./( CONSTANTS(:,89).*power(ALGEBRAIC(:,36), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,367)+STATES(:,3)./CONSTANTS(:,368)).*STATES(:,27))./( CONSTANTS(:,369).*power(ALGEBRAIC(:,43), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,370)+STATES(:,3)./CONSTANTS(:,95)).*STATES(:,29))./( CONSTANTS(:,96).*power(ALGEBRAIC(:,44), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,371)+STATES(:,3)./CONSTANTS(:,372)).*STATES(:,30))./( CONSTANTS(:,373).*power(ALGEBRAIC(:,45), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,377)+STATES(:,3)./CONSTANTS(:,378)).*STATES(:,32))./( CONSTANTS(:,103).*power(ALGEBRAIC(:,49), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,379)+STATES(:,3)./CONSTANTS(:,380)).*STATES(:,31))./( CONSTANTS(:,104).*power(ALGEBRAIC(:,50), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,141)+STATES(:,3)./CONSTANTS(:,142)).*STATES(:,4))./( CONSTANTS(:,143).*power(CONSTANTS(:,285), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,138)+STATES(:,3)./CONSTANTS(:,139)).*STATES(:,3))./( CONSTANTS(:,140).*power(CONSTANTS(:,284), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,126)+STATES(:,3)./CONSTANTS(:,127)).*STATES(:,57))./( CONSTANTS(:,128).*power(CONSTANTS(:,280), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,129)+STATES(:,3)./CONSTANTS(:,130)).*CONSTANTS(:,292))./( CONSTANTS(:,131).*power(CONSTANTS(:,281), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,132)+STATES(:,3)./CONSTANTS(:,133)).*CONSTANTS(:,311))./( CONSTANTS(:,134).*power(CONSTANTS(:,282), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,374)+STATES(:,3)./CONSTANTS(:,97)).*STATES(:,20))./( CONSTANTS(:,98).*power(ALGEBRAIC(:,46), 2.00000));
    ALGEBRAIC(:,206) = 1.00000+ALGEBRAIC(:,198);
    ALGEBRAIC(:,195) = ( (1.00000+STATES(:,2)./CONSTANTS(:,135)+STATES(:,4)./CONSTANTS(:,137)).*STATES(:,2))./( CONSTANTS(:,136).*power(CONSTANTS(:,283), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,1))./( CONSTANTS(:,312).*power(ALGEBRAIC(:,1), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,9))./( CONSTANTS(:,337).*power(ALGEBRAIC(:,22), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,345)+STATES(:,4)./CONSTANTS(:,350)).*STATES(:,11))./( CONSTANTS(:,348).*power(ALGEBRAIC(:,27), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,12))./( CONSTANTS(:,312).*power(ALGEBRAIC(:,30), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,13))./( CONSTANTS(:,337).*power(ALGEBRAIC(:,32), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,352)+STATES(:,4)./CONSTANTS(:,356)).*STATES(:,14))./( CONSTANTS(:,354).*power(ALGEBRAIC(:,34), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,111)+STATES(:,4)./CONSTANTS(:,113)).*STATES(:,22))./( CONSTANTS(:,112).*power(CONSTANTS(:,278), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,108)+STATES(:,4)./CONSTANTS(:,110)).*ALGEBRAIC(:,114))./( CONSTANTS(:,109).*power(CONSTANTS(:,277), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,114)+STATES(:,4)./CONSTANTS(:,116)).*STATES(:,28))./( CONSTANTS(:,115).*power(CONSTANTS(:,305), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,117)+STATES(:,4)./CONSTANTS(:,119)).*ALGEBRAIC(:,119))./( CONSTANTS(:,118).*power(CONSTANTS(:,321), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,99)+STATES(:,4)./CONSTANTS(:,100)).*STATES(:,18))./( CONSTANTS(:,375).*power(ALGEBRAIC(:,47), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,90)+STATES(:,4)./CONSTANTS(:,91)).*STATES(:,23))./( CONSTANTS(:,360).*power(ALGEBRAIC(:,37), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,120)+STATES(:,4)./CONSTANTS(:,122)).*STATES(:,21))./( CONSTANTS(:,121).*power(CONSTANTS(:,279), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,361)+STATES(:,4)./CONSTANTS(:,363)).*STATES(:,16))./( CONSTANTS(:,362).*power(ALGEBRAIC(:,38), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,364)+STATES(:,4)./CONSTANTS(:,92)).*STATES(:,25))./( CONSTANTS(:,365).*power(ALGEBRAIC(:,41), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,123)+STATES(:,4)./CONSTANTS(:,125)).*STATES(:,26))./( CONSTANTS(:,124).*power(CONSTANTS(:,307), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,366)+STATES(:,4)./CONSTANTS(:,94)).*STATES(:,24))./( CONSTANTS(:,93).*power(ALGEBRAIC(:,42), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,358)+STATES(:,4)./CONSTANTS(:,89)).*STATES(:,19))./( CONSTANTS(:,88).*power(ALGEBRAIC(:,36), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,367)+STATES(:,4)./CONSTANTS(:,369)).*STATES(:,27))./( CONSTANTS(:,368).*power(ALGEBRAIC(:,43), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,370)+STATES(:,4)./CONSTANTS(:,96)).*STATES(:,29))./( CONSTANTS(:,95).*power(ALGEBRAIC(:,44), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,371)+STATES(:,4)./CONSTANTS(:,373)).*STATES(:,30))./( CONSTANTS(:,372).*power(ALGEBRAIC(:,45), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,377)+STATES(:,4)./CONSTANTS(:,103)).*STATES(:,32))./( CONSTANTS(:,378).*power(ALGEBRAIC(:,49), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,379)+STATES(:,4)./CONSTANTS(:,104)).*STATES(:,31))./( CONSTANTS(:,380).*power(ALGEBRAIC(:,50), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,141)+STATES(:,4)./CONSTANTS(:,143)).*STATES(:,4))./( CONSTANTS(:,142).*power(CONSTANTS(:,285), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,138)+STATES(:,4)./CONSTANTS(:,140)).*STATES(:,3))./( CONSTANTS(:,139).*power(CONSTANTS(:,284), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,126)+STATES(:,4)./CONSTANTS(:,128)).*STATES(:,57))./( CONSTANTS(:,127).*power(CONSTANTS(:,280), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,129)+STATES(:,4)./CONSTANTS(:,131)).*CONSTANTS(:,292))./( CONSTANTS(:,130).*power(CONSTANTS(:,281), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,132)+STATES(:,4)./CONSTANTS(:,134)).*CONSTANTS(:,311))./( CONSTANTS(:,133).*power(CONSTANTS(:,282), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,374)+STATES(:,4)./CONSTANTS(:,98)).*STATES(:,20))./( CONSTANTS(:,97).*power(ALGEBRAIC(:,46), 2.00000));
    ALGEBRAIC(:,205) = 1.00000+ALGEBRAIC(:,195);
    ALGEBRAIC(:,192) = ( (1.00000+STATES(:,3)./CONSTANTS(:,136)+STATES(:,4)./CONSTANTS(:,137)).*STATES(:,2))./( CONSTANTS(:,135).*power(CONSTANTS(:,283), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,312)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,1))./( CONSTANTS(:,297).*power(ALGEBRAIC(:,1), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,337)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,9))./( CONSTANTS(:,333).*power(ALGEBRAIC(:,22), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,348)+STATES(:,4)./CONSTANTS(:,350)).*STATES(:,11))./( CONSTANTS(:,345).*power(ALGEBRAIC(:,27), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,312)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,12))./( CONSTANTS(:,297).*power(ALGEBRAIC(:,30), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,337)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,13))./( CONSTANTS(:,333).*power(ALGEBRAIC(:,32), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,354)+STATES(:,4)./CONSTANTS(:,356)).*STATES(:,14))./( CONSTANTS(:,352).*power(ALGEBRAIC(:,34), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,112)+STATES(:,4)./CONSTANTS(:,113)).*STATES(:,22))./( CONSTANTS(:,111).*power(CONSTANTS(:,278), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,109)+STATES(:,4)./CONSTANTS(:,110)).*ALGEBRAIC(:,114))./( CONSTANTS(:,108).*power(CONSTANTS(:,277), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,115)+STATES(:,4)./CONSTANTS(:,116)).*STATES(:,28))./( CONSTANTS(:,114).*power(CONSTANTS(:,305), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,118)+STATES(:,4)./CONSTANTS(:,119)).*ALGEBRAIC(:,119))./( CONSTANTS(:,117).*power(CONSTANTS(:,321), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,375)+STATES(:,4)./CONSTANTS(:,100)).*STATES(:,18))./( CONSTANTS(:,99).*power(ALGEBRAIC(:,47), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,360)+STATES(:,4)./CONSTANTS(:,91)).*STATES(:,23))./( CONSTANTS(:,90).*power(ALGEBRAIC(:,37), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,121)+STATES(:,4)./CONSTANTS(:,122)).*STATES(:,21))./( CONSTANTS(:,120).*power(CONSTANTS(:,279), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,362)+STATES(:,4)./CONSTANTS(:,363)).*STATES(:,16))./( CONSTANTS(:,361).*power(ALGEBRAIC(:,38), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,365)+STATES(:,4)./CONSTANTS(:,92)).*STATES(:,25))./( CONSTANTS(:,364).*power(ALGEBRAIC(:,41), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,124)+STATES(:,4)./CONSTANTS(:,125)).*STATES(:,26))./( CONSTANTS(:,123).*power(CONSTANTS(:,307), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,93)+STATES(:,4)./CONSTANTS(:,94)).*STATES(:,24))./( CONSTANTS(:,366).*power(ALGEBRAIC(:,42), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,88)+STATES(:,4)./CONSTANTS(:,89)).*STATES(:,19))./( CONSTANTS(:,358).*power(ALGEBRAIC(:,36), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,368)+STATES(:,4)./CONSTANTS(:,369)).*STATES(:,27))./( CONSTANTS(:,367).*power(ALGEBRAIC(:,43), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,95)+STATES(:,4)./CONSTANTS(:,96)).*STATES(:,29))./( CONSTANTS(:,370).*power(ALGEBRAIC(:,44), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,372)+STATES(:,4)./CONSTANTS(:,373)).*STATES(:,30))./( CONSTANTS(:,371).*power(ALGEBRAIC(:,45), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,378)+STATES(:,4)./CONSTANTS(:,103)).*STATES(:,32))./( CONSTANTS(:,377).*power(ALGEBRAIC(:,49), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,380)+STATES(:,4)./CONSTANTS(:,104)).*STATES(:,31))./( CONSTANTS(:,379).*power(ALGEBRAIC(:,50), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,142)+STATES(:,4)./CONSTANTS(:,143)).*STATES(:,4))./( CONSTANTS(:,141).*power(CONSTANTS(:,285), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,139)+STATES(:,4)./CONSTANTS(:,140)).*STATES(:,3))./( CONSTANTS(:,138).*power(CONSTANTS(:,284), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,127)+STATES(:,4)./CONSTANTS(:,128)).*STATES(:,57))./( CONSTANTS(:,126).*power(CONSTANTS(:,280), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,130)+STATES(:,4)./CONSTANTS(:,131)).*CONSTANTS(:,292))./( CONSTANTS(:,129).*power(CONSTANTS(:,281), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,133)+STATES(:,4)./CONSTANTS(:,134)).*CONSTANTS(:,311))./( CONSTANTS(:,132).*power(CONSTANTS(:,282), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,97)+STATES(:,4)./CONSTANTS(:,98)).*STATES(:,20))./( CONSTANTS(:,374).*power(ALGEBRAIC(:,46), 2.00000));
    ALGEBRAIC(:,204) = 1.00000+ALGEBRAIC(:,192)+CONSTANTS(:,272)./( CONSTANTS(:,273).*power(1.00000+STATES(:,2)./CONSTANTS(:,273), 2.00000));
    ALGEBRAIC(:,196) =  - (( STATES(:,4).*STATES(:,2))./( CONSTANTS(:,137).*CONSTANTS(:,135).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,4).*STATES(:,1))./( CONSTANTS(:,326).*CONSTANTS(:,297).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,4).*STATES(:,9))./( CONSTANTS(:,341).*CONSTANTS(:,333).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,4).*STATES(:,11))./( CONSTANTS(:,350).*CONSTANTS(:,345).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,4).*STATES(:,12))./( CONSTANTS(:,326).*CONSTANTS(:,297).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,4).*STATES(:,13))./( CONSTANTS(:,341).*CONSTANTS(:,333).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,4).*STATES(:,14))./( CONSTANTS(:,356).*CONSTANTS(:,352).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,4).*STATES(:,22))./( CONSTANTS(:,113).*CONSTANTS(:,111).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,114))./( CONSTANTS(:,110).*CONSTANTS(:,108).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,4).*STATES(:,28))./( CONSTANTS(:,116).*CONSTANTS(:,114).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,119))./( CONSTANTS(:,119).*CONSTANTS(:,117).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,4).*STATES(:,18))./( CONSTANTS(:,100).*CONSTANTS(:,99).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,4).*STATES(:,23))./( CONSTANTS(:,91).*CONSTANTS(:,90).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,4).*STATES(:,21))./( CONSTANTS(:,122).*CONSTANTS(:,120).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,4).*STATES(:,16))./( CONSTANTS(:,363).*CONSTANTS(:,361).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,4).*STATES(:,25))./( CONSTANTS(:,92).*CONSTANTS(:,364).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,4).*STATES(:,26))./( CONSTANTS(:,125).*CONSTANTS(:,123).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,4).*STATES(:,24))./( CONSTANTS(:,94).*CONSTANTS(:,366).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,4).*STATES(:,19))./( CONSTANTS(:,89).*CONSTANTS(:,358).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,4).*STATES(:,27))./( CONSTANTS(:,369).*CONSTANTS(:,367).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,4).*STATES(:,29))./( CONSTANTS(:,96).*CONSTANTS(:,370).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,4).*STATES(:,30))./( CONSTANTS(:,373).*CONSTANTS(:,371).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,4).*STATES(:,32))./( CONSTANTS(:,103).*CONSTANTS(:,377).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,4).*STATES(:,31))./( CONSTANTS(:,104).*CONSTANTS(:,379).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,4).*STATES(:,4))./( CONSTANTS(:,143).*CONSTANTS(:,141).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,4).*STATES(:,3))./( CONSTANTS(:,140).*CONSTANTS(:,138).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,4).*STATES(:,57))./( CONSTANTS(:,128).*CONSTANTS(:,126).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,4).*CONSTANTS(:,292))./( CONSTANTS(:,131).*CONSTANTS(:,129).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,4).*CONSTANTS(:,311))./( CONSTANTS(:,134).*CONSTANTS(:,132).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,4).*STATES(:,20))./( CONSTANTS(:,98).*CONSTANTS(:,374).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,193) =  - (( STATES(:,3).*STATES(:,2))./( CONSTANTS(:,136).*CONSTANTS(:,135).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,3).*STATES(:,1))./( CONSTANTS(:,312).*CONSTANTS(:,297).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,3).*STATES(:,9))./( CONSTANTS(:,337).*CONSTANTS(:,333).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,3).*STATES(:,11))./( CONSTANTS(:,348).*CONSTANTS(:,345).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,3).*STATES(:,12))./( CONSTANTS(:,312).*CONSTANTS(:,297).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,3).*STATES(:,13))./( CONSTANTS(:,337).*CONSTANTS(:,333).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,3).*STATES(:,14))./( CONSTANTS(:,354).*CONSTANTS(:,352).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,3).*STATES(:,22))./( CONSTANTS(:,112).*CONSTANTS(:,111).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,114))./( CONSTANTS(:,109).*CONSTANTS(:,108).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,3).*STATES(:,28))./( CONSTANTS(:,115).*CONSTANTS(:,114).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,119))./( CONSTANTS(:,118).*CONSTANTS(:,117).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,3).*STATES(:,18))./( CONSTANTS(:,375).*CONSTANTS(:,99).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,3).*STATES(:,23))./( CONSTANTS(:,360).*CONSTANTS(:,90).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,3).*STATES(:,21))./( CONSTANTS(:,121).*CONSTANTS(:,120).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,3).*STATES(:,16))./( CONSTANTS(:,362).*CONSTANTS(:,361).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,3).*STATES(:,25))./( CONSTANTS(:,365).*CONSTANTS(:,364).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,3).*STATES(:,26))./( CONSTANTS(:,124).*CONSTANTS(:,123).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,3).*STATES(:,24))./( CONSTANTS(:,93).*CONSTANTS(:,366).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,3).*STATES(:,19))./( CONSTANTS(:,88).*CONSTANTS(:,358).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,3).*STATES(:,27))./( CONSTANTS(:,368).*CONSTANTS(:,367).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,3).*STATES(:,29))./( CONSTANTS(:,95).*CONSTANTS(:,370).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,3).*STATES(:,30))./( CONSTANTS(:,372).*CONSTANTS(:,371).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,3).*STATES(:,32))./( CONSTANTS(:,378).*CONSTANTS(:,377).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,3).*STATES(:,31))./( CONSTANTS(:,380).*CONSTANTS(:,379).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,3).*STATES(:,4))./( CONSTANTS(:,142).*CONSTANTS(:,141).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,3).*STATES(:,3))./( CONSTANTS(:,139).*CONSTANTS(:,138).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,3).*STATES(:,57))./( CONSTANTS(:,127).*CONSTANTS(:,126).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,3).*CONSTANTS(:,292))./( CONSTANTS(:,130).*CONSTANTS(:,129).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,3).*CONSTANTS(:,311))./( CONSTANTS(:,133).*CONSTANTS(:,132).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,3).*STATES(:,20))./( CONSTANTS(:,97).*CONSTANTS(:,374).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,207) = ((( ALGEBRAIC(:,204).*ALGEBRAIC(:,197).*ALGEBRAIC(:,194)+ ALGEBRAIC(:,206).*ALGEBRAIC(:,190).*ALGEBRAIC(:,193)+ ALGEBRAIC(:,205).*ALGEBRAIC(:,191).*ALGEBRAIC(:,196)) -  ALGEBRAIC(:,205).*ALGEBRAIC(:,206).*ALGEBRAIC(:,204)) -  ALGEBRAIC(:,191).*ALGEBRAIC(:,197).*ALGEBRAIC(:,193)) -  ALGEBRAIC(:,190).*ALGEBRAIC(:,194).*ALGEBRAIC(:,196);
    ALGEBRAIC(:,146) = ((ALGEBRAIC(:,115)+ALGEBRAIC(:,116)+ALGEBRAIC(:,117)+ALGEBRAIC(:,118)) - ALGEBRAIC(:,140))./CONSTANTS(:,4);
    ALGEBRAIC(:,147) =  - ALGEBRAIC(:,146);
    ALGEBRAIC(:,148) = ((ALGEBRAIC(:,120)+ALGEBRAIC(:,140)) - ALGEBRAIC(:,123))./CONSTANTS(:,4);
    ALGEBRAIC(:,149) =  - ALGEBRAIC(:,148);
    ALGEBRAIC(:,184) = ((ALGEBRAIC(:,82)+ALGEBRAIC(:,144)) - ALGEBRAIC(:,179))./CONSTANTS(:,4);
    ALGEBRAIC(:,185) = (( - ALGEBRAIC(:,82) - ALGEBRAIC(:,144))+ALGEBRAIC(:,179))./CONSTANTS(:,4);
    ALGEBRAIC(:,124) = (ALGEBRAIC(:,67) - ALGEBRAIC(:,82))./CONSTANTS(:,4);
    ALGEBRAIC(:,125) = ( - ALGEBRAIC(:,67)+ALGEBRAIC(:,82))./CONSTANTS(:,4);
    ALGEBRAIC(:,186) = ((( - ALGEBRAIC(:,67) - ALGEBRAIC(:,144))+ALGEBRAIC(:,180)) - ALGEBRAIC(:,177))./CONSTANTS(:,4);
    ALGEBRAIC(:,160) = ( - ALGEBRAIC(:,115)+ALGEBRAIC(:,156))./CONSTANTS(:,4);
    ALGEBRAIC(:,126) = (ALGEBRAIC(:,115) - ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,187) = ( - ALGEBRAIC(:,60)+ALGEBRAIC(:,176)+ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,127) = (ALGEBRAIC(:,60) - ALGEBRAIC(:,116))./CONSTANTS(:,4);
    ALGEBRAIC(:,128) = (((ALGEBRAIC(:,116) - ALGEBRAIC(:,117)) - ALGEBRAIC(:,87))+ALGEBRAIC(:,88))./CONSTANTS(:,4);
    ALGEBRAIC(:,129) = (ALGEBRAIC(:,117) - ALGEBRAIC(:,67))./CONSTANTS(:,4);
    ALGEBRAIC(:,130) = (( - ALGEBRAIC(:,115) - ALGEBRAIC(:,117))+ALGEBRAIC(:,67)+ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,188) = ((ALGEBRAIC(:,67) - ALGEBRAIC(:,120))+ALGEBRAIC(:,178))./CONSTANTS(:,4);
    ALGEBRAIC(:,131) = (ALGEBRAIC(:,120) - ALGEBRAIC(:,74))./CONSTANTS(:,4);
    ALGEBRAIC(:,189) = (((((ALGEBRAIC(:,74) - ALGEBRAIC(:,118))+ALGEBRAIC(:,177)) - ALGEBRAIC(:,88)) - ALGEBRAIC(:,176)) - ALGEBRAIC(:,178))./CONSTANTS(:,4);
    ALGEBRAIC(:,132) = ( - ALGEBRAIC(:,57)+ALGEBRAIC(:,118)+ALGEBRAIC(:,87))./CONSTANTS(:,4);
    ALGEBRAIC(:,161) = ((ALGEBRAIC(:,87)+ALGEBRAIC(:,157)) - ALGEBRAIC(:,159))./CONSTANTS(:,4);
    ALGEBRAIC(:,162) = ( - ALGEBRAIC(:,87)+ALGEBRAIC(:,159))./CONSTANTS(:,4);
    ALGEBRAIC(:,182) =  CONSTANTS(:,51).*( ALGEBRAIC(:,163).*STATES(:,2) -  STATES(:,4).*ALGEBRAIC(:,133));
    ALGEBRAIC(:,199) = ((((((((( - ALGEBRAIC(:,115)+ 2.00000.*ALGEBRAIC(:,57)) - ALGEBRAIC(:,117))+ALGEBRAIC(:,67)+ALGEBRAIC(:,118)+ALGEBRAIC(:,156)+ALGEBRAIC(:,157)+ALGEBRAIC(:,176)) - ALGEBRAIC(:,159)) -  5.00000.*ALGEBRAIC(:,140)) -  2.00000.*ALGEBRAIC(:,123)) -  4.00000.*ALGEBRAIC(:,142))+ (CONSTANTS(:,8) - 1.00000).*ALGEBRAIC(:,144)+ 2.00000.*ALGEBRAIC(:,180)+ALGEBRAIC(:,145)) - ALGEBRAIC(:,182))./CONSTANTS(:,4);
    ALGEBRAIC(:,201) =  - (( STATES(:,2).*CONSTANTS(:,293))./( CONSTANTS(:,135).*CONSTANTS(:,283))+( STATES(:,2).*ALGEBRAIC(:,184))./( CONSTANTS(:,297).*ALGEBRAIC(:,1))+( STATES(:,2).*ALGEBRAIC(:,185))./( CONSTANTS(:,333).*ALGEBRAIC(:,22))+( STATES(:,2).*CONSTANTS(:,263))./( CONSTANTS(:,345).*ALGEBRAIC(:,27))+( STATES(:,2).*ALGEBRAIC(:,124))./( CONSTANTS(:,297).*ALGEBRAIC(:,30))+( STATES(:,2).*ALGEBRAIC(:,125))./( CONSTANTS(:,333).*ALGEBRAIC(:,32))+( STATES(:,2).*ALGEBRAIC(:,186))./( CONSTANTS(:,352).*ALGEBRAIC(:,34))+( STATES(:,2).*ALGEBRAIC(:,146))./( CONSTANTS(:,111).*CONSTANTS(:,278))+( STATES(:,2).*ALGEBRAIC(:,147))./( CONSTANTS(:,108).*CONSTANTS(:,277))+( STATES(:,2).*ALGEBRAIC(:,148))./( CONSTANTS(:,114).*CONSTANTS(:,305))+( STATES(:,2).*ALGEBRAIC(:,149))./( CONSTANTS(:,117).*CONSTANTS(:,321))+( STATES(:,2).*ALGEBRAIC(:,160))./( CONSTANTS(:,99).*ALGEBRAIC(:,47))+( STATES(:,2).*ALGEBRAIC(:,132))./( CONSTANTS(:,90).*ALGEBRAIC(:,37))+( STATES(:,2).*ALGEBRAIC(:,126))./( CONSTANTS(:,120).*CONSTANTS(:,279))+( STATES(:,2).*ALGEBRAIC(:,187))./( CONSTANTS(:,361).*ALGEBRAIC(:,38))+( STATES(:,2).*ALGEBRAIC(:,127))./( CONSTANTS(:,364).*ALGEBRAIC(:,41))+( STATES(:,2).*ALGEBRAIC(:,128))./( CONSTANTS(:,123).*CONSTANTS(:,307))+( STATES(:,2).*ALGEBRAIC(:,129))./( CONSTANTS(:,366).*ALGEBRAIC(:,42))+( STATES(:,2).*ALGEBRAIC(:,130))./( CONSTANTS(:,358).*ALGEBRAIC(:,36))+( STATES(:,2).*ALGEBRAIC(:,188))./( CONSTANTS(:,367).*ALGEBRAIC(:,43))+( STATES(:,2).*ALGEBRAIC(:,131))./( CONSTANTS(:,370).*ALGEBRAIC(:,44))+( STATES(:,2).*ALGEBRAIC(:,189))./( CONSTANTS(:,371).*ALGEBRAIC(:,45))+( STATES(:,2).*ALGEBRAIC(:,161))./( CONSTANTS(:,377).*ALGEBRAIC(:,49))+( STATES(:,2).*ALGEBRAIC(:,162))./( CONSTANTS(:,379).*ALGEBRAIC(:,50))+( STATES(:,2).*CONSTANTS(:,295))./( CONSTANTS(:,141).*CONSTANTS(:,285))+( STATES(:,2).*CONSTANTS(:,294))./( CONSTANTS(:,138).*CONSTANTS(:,284))+( STATES(:,2).*CONSTANTS(:,264))./( CONSTANTS(:,126).*CONSTANTS(:,280))+( STATES(:,2).*CONSTANTS(:,271))./( CONSTANTS(:,129).*CONSTANTS(:,281))+( STATES(:,2).*CONSTANTS(:,270))./( CONSTANTS(:,132).*CONSTANTS(:,282))+( STATES(:,2).*CONSTANTS(:,265))./( CONSTANTS(:,374).*ALGEBRAIC(:,46)))+ALGEBRAIC(:,199);
    ALGEBRAIC(:,202) =  - (( STATES(:,3).*CONSTANTS(:,293))./( CONSTANTS(:,136).*CONSTANTS(:,283))+( STATES(:,3).*ALGEBRAIC(:,184))./( CONSTANTS(:,312).*ALGEBRAIC(:,1))+( STATES(:,3).*ALGEBRAIC(:,185))./( CONSTANTS(:,337).*ALGEBRAIC(:,22))+( STATES(:,3).*CONSTANTS(:,263))./( CONSTANTS(:,348).*ALGEBRAIC(:,27))+( STATES(:,3).*ALGEBRAIC(:,124))./( CONSTANTS(:,312).*ALGEBRAIC(:,30))+( STATES(:,3).*ALGEBRAIC(:,125))./( CONSTANTS(:,337).*ALGEBRAIC(:,32))+( STATES(:,3).*ALGEBRAIC(:,186))./( CONSTANTS(:,354).*ALGEBRAIC(:,34))+( STATES(:,3).*ALGEBRAIC(:,146))./( CONSTANTS(:,112).*CONSTANTS(:,278))+( STATES(:,3).*ALGEBRAIC(:,147))./( CONSTANTS(:,109).*CONSTANTS(:,277))+( STATES(:,3).*ALGEBRAIC(:,148))./( CONSTANTS(:,115).*CONSTANTS(:,305))+( STATES(:,3).*ALGEBRAIC(:,149))./( CONSTANTS(:,118).*CONSTANTS(:,321))+( STATES(:,3).*ALGEBRAIC(:,160))./( CONSTANTS(:,375).*ALGEBRAIC(:,47))+( STATES(:,3).*ALGEBRAIC(:,132))./( CONSTANTS(:,360).*ALGEBRAIC(:,37))+( STATES(:,3).*ALGEBRAIC(:,126))./( CONSTANTS(:,121).*CONSTANTS(:,279))+( STATES(:,3).*ALGEBRAIC(:,187))./( CONSTANTS(:,362).*ALGEBRAIC(:,38))+( STATES(:,3).*ALGEBRAIC(:,127))./( CONSTANTS(:,365).*ALGEBRAIC(:,41))+( STATES(:,3).*ALGEBRAIC(:,128))./( CONSTANTS(:,124).*CONSTANTS(:,307))+( STATES(:,3).*ALGEBRAIC(:,129))./( CONSTANTS(:,93).*ALGEBRAIC(:,42))+( STATES(:,3).*ALGEBRAIC(:,130))./( CONSTANTS(:,88).*ALGEBRAIC(:,36))+( STATES(:,3).*ALGEBRAIC(:,188))./( CONSTANTS(:,368).*ALGEBRAIC(:,43))+( STATES(:,3).*ALGEBRAIC(:,131))./( CONSTANTS(:,95).*ALGEBRAIC(:,44))+( STATES(:,3).*ALGEBRAIC(:,189))./( CONSTANTS(:,372).*ALGEBRAIC(:,45))+( STATES(:,3).*ALGEBRAIC(:,161))./( CONSTANTS(:,378).*ALGEBRAIC(:,49))+( STATES(:,3).*ALGEBRAIC(:,162))./( CONSTANTS(:,380).*ALGEBRAIC(:,50))+( STATES(:,3).*CONSTANTS(:,295))./( CONSTANTS(:,142).*CONSTANTS(:,285))+( STATES(:,3).*CONSTANTS(:,294))./( CONSTANTS(:,139).*CONSTANTS(:,284))+( STATES(:,3).*CONSTANTS(:,264))./( CONSTANTS(:,127).*CONSTANTS(:,280))+( STATES(:,3).*CONSTANTS(:,271))./( CONSTANTS(:,130).*CONSTANTS(:,281))+( STATES(:,3).*CONSTANTS(:,270))./( CONSTANTS(:,133).*CONSTANTS(:,282))+( STATES(:,3).*CONSTANTS(:,265))./( CONSTANTS(:,97).*ALGEBRAIC(:,46)))+CONSTANTS(:,296);
    ALGEBRAIC(:,200) = ALGEBRAIC(:,182)./CONSTANTS(:,4);
    ALGEBRAIC(:,203) =  - (( STATES(:,4).*CONSTANTS(:,293))./( CONSTANTS(:,137).*CONSTANTS(:,283))+( STATES(:,4).*ALGEBRAIC(:,184))./( CONSTANTS(:,326).*ALGEBRAIC(:,1))+( STATES(:,4).*ALGEBRAIC(:,185))./( CONSTANTS(:,341).*ALGEBRAIC(:,22))+( STATES(:,4).*CONSTANTS(:,263))./( CONSTANTS(:,350).*ALGEBRAIC(:,27))+( STATES(:,4).*ALGEBRAIC(:,124))./( CONSTANTS(:,326).*ALGEBRAIC(:,30))+( STATES(:,4).*ALGEBRAIC(:,125))./( CONSTANTS(:,341).*ALGEBRAIC(:,32))+( STATES(:,4).*ALGEBRAIC(:,186))./( CONSTANTS(:,356).*ALGEBRAIC(:,34))+( STATES(:,4).*ALGEBRAIC(:,146))./( CONSTANTS(:,113).*CONSTANTS(:,278))+( STATES(:,4).*ALGEBRAIC(:,147))./( CONSTANTS(:,110).*CONSTANTS(:,277))+( STATES(:,4).*ALGEBRAIC(:,148))./( CONSTANTS(:,116).*CONSTANTS(:,305))+( STATES(:,4).*ALGEBRAIC(:,149))./( CONSTANTS(:,119).*CONSTANTS(:,321))+( STATES(:,4).*ALGEBRAIC(:,160))./( CONSTANTS(:,100).*ALGEBRAIC(:,47))+( STATES(:,4).*ALGEBRAIC(:,132))./( CONSTANTS(:,91).*ALGEBRAIC(:,37))+( STATES(:,4).*ALGEBRAIC(:,126))./( CONSTANTS(:,122).*CONSTANTS(:,279))+( STATES(:,4).*ALGEBRAIC(:,187))./( CONSTANTS(:,363).*ALGEBRAIC(:,38))+( STATES(:,4).*ALGEBRAIC(:,127))./( CONSTANTS(:,92).*ALGEBRAIC(:,41))+( STATES(:,4).*ALGEBRAIC(:,128))./( CONSTANTS(:,125).*CONSTANTS(:,307))+( STATES(:,4).*ALGEBRAIC(:,129))./( CONSTANTS(:,94).*ALGEBRAIC(:,42))+( STATES(:,4).*ALGEBRAIC(:,130))./( CONSTANTS(:,89).*ALGEBRAIC(:,36))+( STATES(:,4).*ALGEBRAIC(:,188))./( CONSTANTS(:,369).*ALGEBRAIC(:,43))+( STATES(:,4).*ALGEBRAIC(:,131))./( CONSTANTS(:,96).*ALGEBRAIC(:,44))+( STATES(:,4).*ALGEBRAIC(:,189))./( CONSTANTS(:,373).*ALGEBRAIC(:,45))+( STATES(:,4).*ALGEBRAIC(:,161))./( CONSTANTS(:,103).*ALGEBRAIC(:,49))+( STATES(:,4).*ALGEBRAIC(:,162))./( CONSTANTS(:,104).*ALGEBRAIC(:,50))+( STATES(:,4).*CONSTANTS(:,295))./( CONSTANTS(:,143).*CONSTANTS(:,285))+( STATES(:,4).*CONSTANTS(:,294))./( CONSTANTS(:,140).*CONSTANTS(:,284))+( STATES(:,4).*CONSTANTS(:,264))./( CONSTANTS(:,128).*CONSTANTS(:,280))+( STATES(:,4).*CONSTANTS(:,271))./( CONSTANTS(:,131).*CONSTANTS(:,281))+( STATES(:,4).*CONSTANTS(:,270))./( CONSTANTS(:,134).*CONSTANTS(:,282))+( STATES(:,4).*CONSTANTS(:,265))./( CONSTANTS(:,98).*ALGEBRAIC(:,46)))+ALGEBRAIC(:,200);
    RATES(:,2) = ( ( ALGEBRAIC(:,197).*ALGEBRAIC(:,194) -  ALGEBRAIC(:,205).*ALGEBRAIC(:,206)).*ALGEBRAIC(:,201)+ ( ALGEBRAIC(:,206).*ALGEBRAIC(:,190) -  ALGEBRAIC(:,191).*ALGEBRAIC(:,197)).*ALGEBRAIC(:,202)+ ( ALGEBRAIC(:,205).*ALGEBRAIC(:,191) -  ALGEBRAIC(:,190).*ALGEBRAIC(:,194)).*ALGEBRAIC(:,203))./ALGEBRAIC(:,207);
    RATES(:,3) = ( ( ALGEBRAIC(:,206).*ALGEBRAIC(:,193) -  ALGEBRAIC(:,196).*ALGEBRAIC(:,194)).*ALGEBRAIC(:,201)+ ( ALGEBRAIC(:,196).*ALGEBRAIC(:,191) -  ALGEBRAIC(:,204).*ALGEBRAIC(:,206)).*ALGEBRAIC(:,202)+ ( ALGEBRAIC(:,204).*ALGEBRAIC(:,194) -  ALGEBRAIC(:,191).*ALGEBRAIC(:,193)).*ALGEBRAIC(:,203))./ALGEBRAIC(:,207);
    RATES(:,4) = ( ( ALGEBRAIC(:,205).*ALGEBRAIC(:,196) -  ALGEBRAIC(:,197).*ALGEBRAIC(:,193)).*ALGEBRAIC(:,201)+ ( ALGEBRAIC(:,204).*ALGEBRAIC(:,197) -  ALGEBRAIC(:,196).*ALGEBRAIC(:,190)).*ALGEBRAIC(:,202)+ ( ALGEBRAIC(:,193).*ALGEBRAIC(:,190) -  ALGEBRAIC(:,204).*ALGEBRAIC(:,205)).*ALGEBRAIC(:,203))./ALGEBRAIC(:,207);
   RATES = RATES';
end

% Calculate algebraic variables
function ALGEBRAIC = computeAlgebraic(ALGEBRAIC, CONSTANTS, STATES, VOI)
    statesSize = size(STATES);
    statesColumnCount = statesSize(2);
    if ( statesColumnCount == 1)
        STATES = STATES';
        utilOnes = 1;
    else
        statesRowCount = statesSize(1);
        utilOnes = ones(statesRowCount, 1);
    end
    ALGEBRAIC(:,17) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,49) - STATES(:,50));
    ALGEBRAIC(:,16) =  CONSTANTS(:,6).*CONSTANTS(:,247).*(STATES(:,47) - STATES(:,48));
    ALGEBRAIC(:,1) = 1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,3)./CONSTANTS(:,312)+STATES(:,4)./CONSTANTS(:,326);
    ALGEBRAIC(:,30) = ALGEBRAIC(:,1);
    ALGEBRAIC(:,22) = 1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,3)./CONSTANTS(:,337)+STATES(:,4)./CONSTANTS(:,341);
    ALGEBRAIC(:,32) = ALGEBRAIC(:,22);
    ALGEBRAIC(:,34) = 1.00000+STATES(:,2)./CONSTANTS(:,352)+STATES(:,3)./CONSTANTS(:,354)+STATES(:,4)./CONSTANTS(:,356);
    ALGEBRAIC(:,36) = 1.00000+STATES(:,2)./CONSTANTS(:,358);
    ALGEBRAIC(:,42) = 1.00000+STATES(:,2)./CONSTANTS(:,366);
    ALGEBRAIC(:,43) = 1.00000+STATES(:,2)./CONSTANTS(:,367)+STATES(:,3)./CONSTANTS(:,368)+STATES(:,4)./CONSTANTS(:,369);
    ALGEBRAIC(:,65) = ( CONSTANTS(:,320).*ALGEBRAIC(:,36).*ALGEBRAIC(:,43).*ALGEBRAIC(:,30))./( STATES(:,2).*ALGEBRAIC(:,32).*ALGEBRAIC(:,42).*ALGEBRAIC(:,34));
    ALGEBRAIC(:,66) = ( CONSTANTS(:,29).*CONSTANTS(:,187).*CONSTANTS(:,182).*CONSTANTS(:,183))./( ALGEBRAIC(:,65).*CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,186));
    ALGEBRAIC(:,67) = ( CONSTANTS(:,29).*ALGEBRAIC(:,66).*( STATES(:,13).*STATES(:,24).*STATES(:,14) - ( STATES(:,19).*STATES(:,27).*STATES(:,12))./ALGEBRAIC(:,65)))./( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,186)+ ALGEBRAIC(:,66).*CONSTANTS(:,179).*CONSTANTS(:,186).*STATES(:,13)+ ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14)+ ALGEBRAIC(:,66).*CONSTANTS(:,186).*STATES(:,13).*STATES(:,24)+ ALGEBRAIC(:,66).*CONSTANTS(:,185).*STATES(:,13).*STATES(:,14)+ ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14)+ ALGEBRAIC(:,66).*STATES(:,13).*STATES(:,24).*STATES(:,14)+( CONSTANTS(:,29).*CONSTANTS(:,183).*CONSTANTS(:,188).*STATES(:,19))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,182).*CONSTANTS(:,187).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,19).*STATES(:,27))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,188).*STATES(:,19).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,187).*STATES(:,27).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*STATES(:,19).*STATES(:,27).*STATES(:,12))./ALGEBRAIC(:,65)+( CONSTANTS(:,29).*CONSTANTS(:,188).*CONSTANTS(:,183).*STATES(:,13).*STATES(:,19))./( CONSTANTS(:,178).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14).*STATES(:,12))./CONSTANTS(:,183)+( CONSTANTS(:,29).*CONSTANTS(:,188).*CONSTANTS(:,183).*STATES(:,13).*STATES(:,24).*STATES(:,19))./( CONSTANTS(:,178).*CONSTANTS(:,179).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14).*STATES(:,12))./CONSTANTS(:,183)+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,19).*STATES(:,27))./( CONSTANTS(:,178).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,182).*CONSTANTS(:,183))+( CONSTANTS(:,29).*CONSTANTS(:,183).*CONSTANTS(:,188).*STATES(:,13).*STATES(:,24).*STATES(:,14).*STATES(:,19))./( CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,180).*ALGEBRAIC(:,65))+( CONSTANTS(:,29).*CONSTANTS(:,181).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,24).*STATES(:,14).*STATES(:,27))./( CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,180).*ALGEBRAIC(:,65))+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,24).*STATES(:,19).*STATES(:,27))./( CONSTANTS(:,178).*CONSTANTS(:,179).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,182).*CONSTANTS(:,183))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*CONSTANTS(:,180).*STATES(:,24).*STATES(:,19).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,181).*CONSTANTS(:,182).*CONSTANTS(:,183))+( ALGEBRAIC(:,66).*CONSTANTS(:,178).*CONSTANTS(:,185).*STATES(:,14).*STATES(:,19).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,181).*CONSTANTS(:,182).*CONSTANTS(:,183))+( CONSTANTS(:,29).*CONSTANTS(:,189).*STATES(:,13).*STATES(:,24).*STATES(:,14).*STATES(:,19).*STATES(:,27))./( CONSTANTS(:,178).*CONSTANTS(:,179).*CONSTANTS(:,180).*ALGEBRAIC(:,65))+( ALGEBRAIC(:,66).*CONSTANTS(:,184).*STATES(:,24).*STATES(:,14).*STATES(:,19).*STATES(:,27).*STATES(:,12))./( CONSTANTS(:,181).*CONSTANTS(:,182).*CONSTANTS(:,183)));
    ALGEBRAIC(:,78) = ( CONSTANTS(:,313).*ALGEBRAIC(:,32).*ALGEBRAIC(:,1))./( ALGEBRAIC(:,30).*ALGEBRAIC(:,22));
    ALGEBRAIC(:,79) = ( CONSTANTS(:,33).*CONSTANTS(:,224).*CONSTANTS(:,219))./( ALGEBRAIC(:,78).*CONSTANTS(:,217).*CONSTANTS(:,222));
    ALGEBRAIC(:,27) = 1.00000+STATES(:,2)./CONSTANTS(:,345)+STATES(:,3)./CONSTANTS(:,348)+STATES(:,4)./CONSTANTS(:,350);
    ALGEBRAIC(:,29) = ( STATES(:,11).*(1.00000+STATES(:,2)./CONSTANTS(:,345)))./ALGEBRAIC(:,27);
    ALGEBRAIC(:,80) = 1.00000+ALGEBRAIC(:,29)./CONSTANTS(:,225);
    ALGEBRAIC(:,81) = piecewise({STATES(:,12)>CONSTANTS(:,226)&STATES(:,9)>CONSTANTS(:,226), ( (( CONSTANTS(:,33).*ALGEBRAIC(:,79))./ALGEBRAIC(:,80)).*STATES(:,12).*STATES(:,9))./( ALGEBRAIC(:,79).*CONSTANTS(:,222).*STATES(:,12)+ ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,13))./ALGEBRAIC(:,78)+( CONSTANTS(:,33).*CONSTANTS(:,223).*STATES(:,1))./ALGEBRAIC(:,78)+ ALGEBRAIC(:,79).*STATES(:,12).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,12).*STATES(:,13))./( ALGEBRAIC(:,78).*CONSTANTS(:,217))+( CONSTANTS(:,33).*STATES(:,13).*STATES(:,1))./ALGEBRAIC(:,78)+( ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9).*STATES(:,1))./CONSTANTS(:,220)) }, 0.00000);
    ALGEBRAIC(:,82) = piecewise({STATES(:,13)>CONSTANTS(:,226)&STATES(:,1)>CONSTANTS(:,226), ALGEBRAIC(:,81) - (( (( CONSTANTS(:,33).*ALGEBRAIC(:,79))./ALGEBRAIC(:,80)).*STATES(:,13).*STATES(:,1))./ALGEBRAIC(:,78))./( ALGEBRAIC(:,79).*CONSTANTS(:,222).*STATES(:,12)+ ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,13))./ALGEBRAIC(:,78)+( CONSTANTS(:,33).*CONSTANTS(:,223).*STATES(:,1))./ALGEBRAIC(:,78)+ ALGEBRAIC(:,79).*STATES(:,12).*STATES(:,9)+( CONSTANTS(:,33).*CONSTANTS(:,224).*STATES(:,12).*STATES(:,13))./( ALGEBRAIC(:,78).*CONSTANTS(:,217))+( CONSTANTS(:,33).*STATES(:,13).*STATES(:,1))./ALGEBRAIC(:,78)+( ALGEBRAIC(:,79).*CONSTANTS(:,221).*STATES(:,9).*STATES(:,1))./CONSTANTS(:,220)) }, ALGEBRAIC(:,81));
    ALGEBRAIC(:,88) = ( CONSTANTS(:,38).*( STATES(:,36).*STATES(:,30) -  STATES(:,26).*STATES(:,35)))./( CONSTANTS(:,237).*CONSTANTS(:,240).*(2.00000+STATES(:,35)./CONSTANTS(:,239)+STATES(:,30)./CONSTANTS(:,240)+STATES(:,36)./CONSTANTS(:,237)+STATES(:,26)./CONSTANTS(:,238)+( STATES(:,35).*STATES(:,26))./( CONSTANTS(:,239).*CONSTANTS(:,238))+( STATES(:,30).*STATES(:,36))./( CONSTANTS(:,240).*CONSTANTS(:,237))));
    ALGEBRAIC(:,90) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,43) - STATES(:,36));
    ALGEBRAIC(:,89) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,40) - STATES(:,33));
    ALGEBRAIC(:,91) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,41) - STATES(:,17));
    ALGEBRAIC(:,92) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,42) - STATES(:,35));
    ALGEBRAIC(:,93) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,44) - STATES(:,39));
    ALGEBRAIC(:,94) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,45) - STATES(:,34));
    ALGEBRAIC(:,95) =  CONSTANTS(:,6).*CONSTANTS(:,11).*(STATES(:,46) - STATES(:,37));
    ALGEBRAIC(:,99) =  CONSTANTS(:,6).*CONSTANTS(:,9).*(STATES(:,55) - STATES(:,15));
    ALGEBRAIC(:,18) = 1.00000+STATES(:,6)./CONSTANTS(:,297)+STATES(:,7)./CONSTANTS(:,312)+STATES(:,8)./CONSTANTS(:,326);
    ALGEBRAIC(:,21) = (( STATES(:,5).*STATES(:,7))./CONSTANTS(:,312))./ALGEBRAIC(:,18);
    ALGEBRAIC(:,23) = 1.00000+STATES(:,6)./CONSTANTS(:,333)+STATES(:,7)./CONSTANTS(:,337)+STATES(:,8)./CONSTANTS(:,341);
    ALGEBRAIC(:,26) = (( STATES(:,10).*STATES(:,7))./CONSTANTS(:,337))./ALGEBRAIC(:,23);
    ALGEBRAIC(:,48) = 1.00000+STATES(:,6)./CONSTANTS(:,376);
    ALGEBRAIC(:,100) = ( CONSTANTS(:,325).*ALGEBRAIC(:,48).*ALGEBRAIC(:,23))./( STATES(:,6).*ALGEBRAIC(:,18).*CONSTANTS(:,276));
    ALGEBRAIC(:,102) = ( ALGEBRAIC(:,100).*CONSTANTS(:,312).*ALGEBRAIC(:,18))./( CONSTANTS(:,337).*ALGEBRAIC(:,23));
    ALGEBRAIC(:,104) = ( (CONSTANTS(:,274)./( CONSTANTS(:,253).*CONSTANTS(:,256))).*( ALGEBRAIC(:,21).*STATES(:,58) - ( STATES(:,59).*ALGEBRAIC(:,26))./ALGEBRAIC(:,102)))./(1.00000+ALGEBRAIC(:,21)./CONSTANTS(:,252)+STATES(:,58)./CONSTANTS(:,253)+( ALGEBRAIC(:,21).*STATES(:,58))./( CONSTANTS(:,253).*CONSTANTS(:,256))+STATES(:,59)./CONSTANTS(:,254)+ALGEBRAIC(:,26)./CONSTANTS(:,255)+( STATES(:,59).*ALGEBRAIC(:,26))./( CONSTANTS(:,255).*CONSTANTS(:,258))+( STATES(:,59).*STATES(:,58))./( CONSTANTS(:,260).*CONSTANTS(:,253)));
    ALGEBRAIC(:,106) = CONSTANTS(:,19) - STATES(:,60);
    ALGEBRAIC(:,108) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000, 0.00000 },  CONSTANTS(:,261).*( (( CONSTANTS(:,310).*STATES(:,10))./ALGEBRAIC(:,23)).*STATES(:,60).*STATES(:,6) -  (STATES(:,5)./ALGEBRAIC(:,18)).*ALGEBRAIC(:,106)));
    ALGEBRAIC(:,97) =  CONSTANTS(:,6).*CONSTANTS(:,10).*(STATES(:,5) - STATES(:,52));
    ALGEBRAIC(:,28) = 1.00000+STATES(:,6)./CONSTANTS(:,345)+STATES(:,7)./CONSTANTS(:,348)+STATES(:,8)./CONSTANTS(:,350);
    ALGEBRAIC(:,110) = ( CONSTANTS(:,382).*ALGEBRAIC(:,18).*ALGEBRAIC(:,28))./( ALGEBRAIC(:,23).*ALGEBRAIC(:,23));
    ALGEBRAIC(:,112) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000, 0.00000 },  CONSTANTS(:,262).*( ALGEBRAIC(:,110).*STATES(:,10).*STATES(:,10) -  STATES(:,53).*STATES(:,5)));
    ALGEBRAIC(:,96) =  CONSTANTS(:,6).*CONSTANTS(:,10).*(STATES(:,10) - STATES(:,51));
    ALGEBRAIC(:,98) =  CONSTANTS(:,6).*CONSTANTS(:,10).*(STATES(:,53) - STATES(:,54));
    ALGEBRAIC(:,114) = CONSTANTS(:,15) - STATES(:,22);
    ALGEBRAIC(:,46) = 1.00000+STATES(:,2)./CONSTANTS(:,374);
    ALGEBRAIC(:,47) = 1.00000+STATES(:,3)./CONSTANTS(:,375);
    ALGEBRAIC(:,51) = ( CONSTANTS(:,315).*ALGEBRAIC(:,46).*CONSTANTS(:,279).*CONSTANTS(:,278))./( STATES(:,2).*ALGEBRAIC(:,47).*ALGEBRAIC(:,36).*CONSTANTS(:,277));
    ALGEBRAIC(:,52) = 1.00000+STATES(:,21)./CONSTANTS(:,147);
    ALGEBRAIC(:,53) = 1.00000+STATES(:,22)./CONSTANTS(:,148);
    ALGEBRAIC(:,115) = piecewise({STATES(:,18)>CONSTANTS(:,149)&STATES(:,19)>CONSTANTS(:,149)&ALGEBRAIC(:,114)>CONSTANTS(:,149), ( CONSTANTS(:,23).*( STATES(:,18).*STATES(:,19).*ALGEBRAIC(:,114) - ( STATES(:,20).*STATES(:,21).*STATES(:,22))./ALGEBRAIC(:,51)))./( CONSTANTS(:,146).*ALGEBRAIC(:,53).*STATES(:,18).*STATES(:,19)+ CONSTANTS(:,145).*ALGEBRAIC(:,52).*STATES(:,18).*ALGEBRAIC(:,114)+ CONSTANTS(:,144).*STATES(:,19).*ALGEBRAIC(:,114)+ STATES(:,18).*STATES(:,19).*ALGEBRAIC(:,114)) }, 0.00000);
    ALGEBRAIC(:,37) = 1.00000+STATES(:,3)./CONSTANTS(:,360);
    ALGEBRAIC(:,38) = 1.00000+STATES(:,2)./CONSTANTS(:,361)+STATES(:,3)./CONSTANTS(:,362)+STATES(:,4)./CONSTANTS(:,363);
    ALGEBRAIC(:,54) = ( CONSTANTS(:,316).*ALGEBRAIC(:,38).*ALGEBRAIC(:,36))./( power(STATES(:,2), 2.00000).*CONSTANTS(:,279).*ALGEBRAIC(:,37));
    ALGEBRAIC(:,39) = ( STATES(:,16).*(1.00000+STATES(:,2)./CONSTANTS(:,361)))./ALGEBRAIC(:,38);
    ALGEBRAIC(:,55) = 1.00000+ALGEBRAIC(:,39)./CONSTANTS(:,153);
    ALGEBRAIC(:,19) = ( STATES(:,1).*(1.00000+STATES(:,2)./CONSTANTS(:,297)))./ALGEBRAIC(:,1);
    ALGEBRAIC(:,24) = ( STATES(:,9).*(1.00000+STATES(:,2)./CONSTANTS(:,333)))./ALGEBRAIC(:,22);
    ALGEBRAIC(:,56) = 1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,154)+ALGEBRAIC(:,24)./CONSTANTS(:,155)+ALGEBRAIC(:,29)./CONSTANTS(:,156)+STATES(:,19)./CONSTANTS(:,157)+STATES(:,24)./CONSTANTS(:,158);
    ALGEBRAIC(:,57) = ( CONSTANTS(:,24).*( STATES(:,21).*STATES(:,23) - ( STATES(:,19).*STATES(:,16))./ALGEBRAIC(:,54)))./( CONSTANTS(:,152).*CONSTANTS(:,151).*ALGEBRAIC(:,55)+ CONSTANTS(:,150).*ALGEBRAIC(:,55).*STATES(:,21)+ CONSTANTS(:,151).*ALGEBRAIC(:,56).*STATES(:,23)+ STATES(:,23).*STATES(:,21));
    ALGEBRAIC(:,41) = 1.00000+STATES(:,2)./CONSTANTS(:,364)+STATES(:,3)./CONSTANTS(:,365);
    ALGEBRAIC(:,58) = ( CONSTANTS(:,317).*ALGEBRAIC(:,41))./ALGEBRAIC(:,38);
    ALGEBRAIC(:,59) = ( CONSTANTS(:,25).*CONSTANTS(:,160))./( CONSTANTS(:,159).*ALGEBRAIC(:,58));
    ALGEBRAIC(:,60) = ( CONSTANTS(:,25).*ALGEBRAIC(:,59).*(STATES(:,16) - STATES(:,25)./ALGEBRAIC(:,58)))./( CONSTANTS(:,159).*ALGEBRAIC(:,59)+ ALGEBRAIC(:,59).*STATES(:,16)+( CONSTANTS(:,25).*STATES(:,25))./ALGEBRAIC(:,58));
    ALGEBRAIC(:,61) = ( CONSTANTS(:,318).*ALGEBRAIC(:,46).*CONSTANTS(:,307).*CONSTANTS(:,278))./( power(STATES(:,2), 2.00000).*CONSTANTS(:,277).*ALGEBRAIC(:,41));
    ALGEBRAIC(:,62) = 1.00000+ (CONSTANTS(:,167)./ALGEBRAIC(:,24)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,166));
    ALGEBRAIC(:,116) = piecewise({ALGEBRAIC(:,114)>CONSTANTS(:,168)&STATES(:,25)>CONSTANTS(:,168), ( CONSTANTS(:,26).*(1.00000 - ( STATES(:,20).*STATES(:,26).*STATES(:,22))./( ALGEBRAIC(:,61).*ALGEBRAIC(:,114).*STATES(:,25))))./(1.00000+ power(CONSTANTS(:,162)./STATES(:,25), CONSTANTS(:,165)).*ALGEBRAIC(:,62)+ (CONSTANTS(:,161)./ALGEBRAIC(:,114)).*(1.00000+ power(CONSTANTS(:,163)./STATES(:,25), CONSTANTS(:,165)).*ALGEBRAIC(:,62)+( STATES(:,22).*ALGEBRAIC(:,62))./CONSTANTS(:,164))) }, 0.00000);
    ALGEBRAIC(:,63) = ( CONSTANTS(:,319).*ALGEBRAIC(:,46).*ALGEBRAIC(:,42).*CONSTANTS(:,278))./( STATES(:,2).*ALGEBRAIC(:,36).*CONSTANTS(:,307).*CONSTANTS(:,277));
    ALGEBRAIC(:,64) = 1.00000+ (CONSTANTS(:,176)./ALGEBRAIC(:,24)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,175));
    ALGEBRAIC(:,117) = piecewise({STATES(:,26)>CONSTANTS(:,177)&STATES(:,19)>CONSTANTS(:,177)&ALGEBRAIC(:,114)>CONSTANTS(:,177), (( CONSTANTS(:,27).*(1.00000 - ( STATES(:,20).*STATES(:,24).*STATES(:,22))./( ALGEBRAIC(:,63).*STATES(:,26).*STATES(:,19).*ALGEBRAIC(:,114))))./(1.00000+ (CONSTANTS(:,169)./STATES(:,26)).*ALGEBRAIC(:,64)+ (CONSTANTS(:,170)./STATES(:,19)).*(1.00000+STATES(:,24)./CONSTANTS(:,172))+ (CONSTANTS(:,171)./ALGEBRAIC(:,114)).*(1.00000+STATES(:,22)./CONSTANTS(:,173))))./(1.00000+STATES(:,22)./CONSTANTS(:,174)) }, 0.00000);
    ALGEBRAIC(:,49) = 1.00000+STATES(:,2)./CONSTANTS(:,377)+STATES(:,3)./CONSTANTS(:,378);
    ALGEBRAIC(:,50) = 1.00000+STATES(:,2)./CONSTANTS(:,379)+STATES(:,3)./CONSTANTS(:,380);
    ALGEBRAIC(:,83) = ( CONSTANTS(:,329).*ALGEBRAIC(:,49).*ALGEBRAIC(:,37))./( ALGEBRAIC(:,50).*CONSTANTS(:,307));
    ALGEBRAIC(:,84) = ( CONSTANTS(:,34).*CONSTANTS(:,234).*CONSTANTS(:,229))./( ALGEBRAIC(:,83).*CONSTANTS(:,227).*CONSTANTS(:,232));
    ALGEBRAIC(:,85) = 1.00000+STATES(:,26)./CONSTANTS(:,235);
    ALGEBRAIC(:,86) = piecewise({STATES(:,31)>CONSTANTS(:,236)&STATES(:,26)>CONSTANTS(:,236), ( CONSTANTS(:,34).*ALGEBRAIC(:,84).*STATES(:,31).*STATES(:,26))./( ALGEBRAIC(:,84).*CONSTANTS(:,232).*STATES(:,31)+ ALGEBRAIC(:,84).*CONSTANTS(:,231).*ALGEBRAIC(:,85).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*ALGEBRAIC(:,85).*STATES(:,23))./ALGEBRAIC(:,83)+( CONSTANTS(:,34).*CONSTANTS(:,233).*STATES(:,32))./ALGEBRAIC(:,83)+ ALGEBRAIC(:,84).*STATES(:,31).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*STATES(:,31).*STATES(:,23))./( ALGEBRAIC(:,83).*CONSTANTS(:,227))+( CONSTANTS(:,34).*STATES(:,23).*STATES(:,32))./ALGEBRAIC(:,83)+( ALGEBRAIC(:,84).*CONSTANTS(:,231).*STATES(:,26).*STATES(:,32))./CONSTANTS(:,230)) }, 0.00000);
    ALGEBRAIC(:,87) = piecewise({STATES(:,23)>CONSTANTS(:,236)&STATES(:,32)>CONSTANTS(:,236), ALGEBRAIC(:,86) - (( CONSTANTS(:,34).*ALGEBRAIC(:,84).*STATES(:,23).*STATES(:,32))./ALGEBRAIC(:,83))./( ALGEBRAIC(:,84).*CONSTANTS(:,232).*STATES(:,31)+ ALGEBRAIC(:,84).*CONSTANTS(:,231).*ALGEBRAIC(:,85).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*ALGEBRAIC(:,85).*STATES(:,23))./ALGEBRAIC(:,83)+( CONSTANTS(:,34).*CONSTANTS(:,233).*STATES(:,32))./ALGEBRAIC(:,83)+ ALGEBRAIC(:,84).*STATES(:,31).*STATES(:,26)+( CONSTANTS(:,34).*CONSTANTS(:,234).*STATES(:,31).*STATES(:,23))./( ALGEBRAIC(:,83).*CONSTANTS(:,227))+( CONSTANTS(:,34).*STATES(:,23).*STATES(:,32))./ALGEBRAIC(:,83)+( ALGEBRAIC(:,84).*CONSTANTS(:,231).*STATES(:,26).*STATES(:,32))./CONSTANTS(:,230)) }, ALGEBRAIC(:,86));
    ALGEBRAIC(:,45) = 1.00000+STATES(:,2)./CONSTANTS(:,371)+STATES(:,3)./CONSTANTS(:,372)+STATES(:,4)./CONSTANTS(:,373);
    ALGEBRAIC(:,75) = ( CONSTANTS(:,322).*ALGEBRAIC(:,37).*CONSTANTS(:,278))./( STATES(:,2).*ALGEBRAIC(:,45).*CONSTANTS(:,277));
    ALGEBRAIC(:,76) = ( CONSTANTS(:,32).*CONSTANTS(:,209).*CONSTANTS(:,212))./( ALGEBRAIC(:,75).*CONSTANTS(:,206).*CONSTANTS(:,211));
    ALGEBRAIC(:,77) = 1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,214)+ALGEBRAIC(:,24)./CONSTANTS(:,215)+ALGEBRAIC(:,29)./CONSTANTS(:,216);
    ALGEBRAIC(:,118) = ( CONSTANTS(:,32).*ALGEBRAIC(:,76).*( ALGEBRAIC(:,114).*STATES(:,30) - ( STATES(:,23).*STATES(:,22))./ALGEBRAIC(:,75)))./( ALGEBRAIC(:,76).*CONSTANTS(:,206).*CONSTANTS(:,211).*ALGEBRAIC(:,77)+ ALGEBRAIC(:,76).*CONSTANTS(:,211).*ALGEBRAIC(:,114)+ ALGEBRAIC(:,76).*CONSTANTS(:,210).*ALGEBRAIC(:,77).*STATES(:,30)+( CONSTANTS(:,32).*CONSTANTS(:,213).*ALGEBRAIC(:,77).*STATES(:,23))./ALGEBRAIC(:,75)+( CONSTANTS(:,32).*CONSTANTS(:,212).*STATES(:,22))./ALGEBRAIC(:,75)+ ALGEBRAIC(:,76).*ALGEBRAIC(:,114).*STATES(:,30)+( CONSTANTS(:,32).*CONSTANTS(:,213).*ALGEBRAIC(:,114).*STATES(:,23))./( ALGEBRAIC(:,75).*CONSTANTS(:,206))+( CONSTANTS(:,32).*STATES(:,23).*STATES(:,22))./ALGEBRAIC(:,75)+( ALGEBRAIC(:,76).*CONSTANTS(:,210).*STATES(:,30).*STATES(:,22))./CONSTANTS(:,209)+( ALGEBRAIC(:,76).*ALGEBRAIC(:,114).*STATES(:,30).*STATES(:,23))./CONSTANTS(:,208)+( CONSTANTS(:,32).*STATES(:,30).*STATES(:,23).*STATES(:,22))./( CONSTANTS(:,207).*ALGEBRAIC(:,75)));
    ALGEBRAIC(:,119) = CONSTANTS(:,16) - STATES(:,28);
    ALGEBRAIC(:,44) = 1.00000+STATES(:,2)./CONSTANTS(:,370);
    ALGEBRAIC(:,68) = ( CONSTANTS(:,335).*ALGEBRAIC(:,44).*CONSTANTS(:,305))./( ALGEBRAIC(:,43).*CONSTANTS(:,321));
    ALGEBRAIC(:,69) = ( CONSTANTS(:,30).*CONSTANTS(:,194).*CONSTANTS(:,191))./( ALGEBRAIC(:,68).*CONSTANTS(:,190).*CONSTANTS(:,193));
    ALGEBRAIC(:,70) = (1.00000+STATES(:,23)./CONSTANTS(:,196)+STATES(:,27)./CONSTANTS(:,197)+STATES(:,29)./CONSTANTS(:,198))./(1.00000+STATES(:,27)./CONSTANTS(:,197)+STATES(:,29)./CONSTANTS(:,198));
    ALGEBRAIC(:,120) = ( CONSTANTS(:,30).*ALGEBRAIC(:,69).*( STATES(:,27).*ALGEBRAIC(:,119) - ( STATES(:,28).*STATES(:,29))./ALGEBRAIC(:,68)))./( ALGEBRAIC(:,69).*CONSTANTS(:,190).*CONSTANTS(:,193).*ALGEBRAIC(:,70)+ ALGEBRAIC(:,69).*CONSTANTS(:,193).*STATES(:,27)+ ALGEBRAIC(:,69).*CONSTANTS(:,192).*ALGEBRAIC(:,70).*ALGEBRAIC(:,119)+( CONSTANTS(:,30).*CONSTANTS(:,195).*ALGEBRAIC(:,70).*STATES(:,28))./ALGEBRAIC(:,68)+( CONSTANTS(:,30).*CONSTANTS(:,194).*STATES(:,29))./ALGEBRAIC(:,68)+ ALGEBRAIC(:,69).*STATES(:,27).*ALGEBRAIC(:,119)+( CONSTANTS(:,30).*CONSTANTS(:,195).*STATES(:,27).*STATES(:,28))./( ALGEBRAIC(:,68).*CONSTANTS(:,190))+( ALGEBRAIC(:,69).*CONSTANTS(:,192).*ALGEBRAIC(:,119).*STATES(:,29))./CONSTANTS(:,191)+( CONSTANTS(:,30).*STATES(:,28).*STATES(:,29))./ALGEBRAIC(:,68));
    ALGEBRAIC(:,71) = ( CONSTANTS(:,343).*ALGEBRAIC(:,45))./ALGEBRAIC(:,44);
    ALGEBRAIC(:,72) = ( CONSTANTS(:,31).*CONSTANTS(:,200))./( ALGEBRAIC(:,71).*CONSTANTS(:,199));
    ALGEBRAIC(:,31) = ( STATES(:,12).*(1.00000+STATES(:,2)./CONSTANTS(:,314)))./ALGEBRAIC(:,30);
    ALGEBRAIC(:,33) = ( STATES(:,13).*(1.00000+STATES(:,2)./CONSTANTS(:,338)))./ALGEBRAIC(:,32);
    ALGEBRAIC(:,73) = 1.00000+STATES(:,16)./CONSTANTS(:,201)+ALGEBRAIC(:,19)./CONSTANTS(:,202)+ALGEBRAIC(:,24)./CONSTANTS(:,203)+ALGEBRAIC(:,31)./CONSTANTS(:,204)+ALGEBRAIC(:,33)./CONSTANTS(:,205);
    ALGEBRAIC(:,74) = ( ALGEBRAIC(:,72).*CONSTANTS(:,31).*(STATES(:,29) - STATES(:,30)./ALGEBRAIC(:,71)))./( CONSTANTS(:,199).*ALGEBRAIC(:,72).*ALGEBRAIC(:,73)+ ALGEBRAIC(:,72).*STATES(:,29)+( CONSTANTS(:,31).*STATES(:,30))./ALGEBRAIC(:,71));
    ALGEBRAIC(:,133) = STATES(:,6);
    ALGEBRAIC(:,101) =  1.00000.*exp( - (CONSTANTS(:,288)+ 4.00000.*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275));
    ALGEBRAIC(:,139) = ( ALGEBRAIC(:,101).*power(STATES(:,2), 5.00000))./power(ALGEBRAIC(:,133), 4.00000);
    ALGEBRAIC(:,140) =  CONSTANTS(:,42).*( ALGEBRAIC(:,139).*STATES(:,22).*ALGEBRAIC(:,119) -  ALGEBRAIC(:,114).*STATES(:,28));
    ALGEBRAIC(:,103) =  1.00000.*exp( - (CONSTANTS(:,289)+ 2.00000.*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275));
    ALGEBRAIC(:,105) = ( ALGEBRAIC(:,103).*power(STATES(:,2), 2.00000))./power(STATES(:,2), 4.00000);
    ALGEBRAIC(:,107) = piecewise({STATES(:,28)>CONSTANTS(:,248), STATES(:,28) }, CONSTANTS(:,248));
    ALGEBRAIC(:,121) = piecewise({ALGEBRAIC(:,119)>CONSTANTS(:,248), ALGEBRAIC(:,119) }, CONSTANTS(:,248));
    ALGEBRAIC(:,122) = CONSTANTS(:,17) - STATES(:,56);
    ALGEBRAIC(:,123) =  (( CONSTANTS(:,43).*(1.00000+STATES(:,14)./CONSTANTS(:,44)))./(1.00000+STATES(:,14)./CONSTANTS(:,45))).*( power(ALGEBRAIC(:,105), 1.0 ./ 2).*ALGEBRAIC(:,122).*power(ALGEBRAIC(:,107), 1.0 ./ 2) -  STATES(:,56).*power(ALGEBRAIC(:,121), 1.0 ./ 2));
    ALGEBRAIC(:,109) = exp( - (CONSTANTS(:,290)+ 4.00000.*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275));
    ALGEBRAIC(:,141) = ( ALGEBRAIC(:,109).*power(STATES(:,2), 4.00000))./( power(ALGEBRAIC(:,133), 2.00000).*power(1.00000, 2.00000));
    ALGEBRAIC(:,111) = piecewise({STATES(:,57)>CONSTANTS(:,249), STATES(:,57) }, CONSTANTS(:,249));
    ALGEBRAIC(:,142) =  CONSTANTS(:,46).*( (ALGEBRAIC(:,111)./(ALGEBRAIC(:,111)+CONSTANTS(:,14))).*exp(( CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)).*(STATES(:,56)./CONSTANTS(:,17))).*( power(ALGEBRAIC(:,141), 1.00000./2.00000).*STATES(:,56).*power(ALGEBRAIC(:,111)./1.00000, 1.00000./4.00000) - ALGEBRAIC(:,122));
    ALGEBRAIC(:,150) = STATES(:,7);
    ALGEBRAIC(:,153) = 1.00000+ALGEBRAIC(:,150)./CONSTANTS(:,375);
    ALGEBRAIC(:,156) =  CONSTANTS(:,35).*( (STATES(:,33)./ALGEBRAIC(:,153)).*ALGEBRAIC(:,133) -  (STATES(:,18)./ALGEBRAIC(:,47)).*STATES(:,2));
    ALGEBRAIC(:,154) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,377)+ALGEBRAIC(:,150)./CONSTANTS(:,378);
    ALGEBRAIC(:,157) =  CONSTANTS(:,36).*( (STATES(:,34)./ALGEBRAIC(:,154)).*ALGEBRAIC(:,133) -  (STATES(:,32)./ALGEBRAIC(:,49)).*STATES(:,2));
    ALGEBRAIC(:,155) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,379)+ALGEBRAIC(:,150)./CONSTANTS(:,380);
    ALGEBRAIC(:,158) = ( exp((  - CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)).*ALGEBRAIC(:,50).*ALGEBRAIC(:,154))./( ALGEBRAIC(:,155).*ALGEBRAIC(:,49));
    ALGEBRAIC(:,159) = ( CONSTANTS(:,41).*( ALGEBRAIC(:,158).*STATES(:,37).*STATES(:,32).*STATES(:,2) -  STATES(:,31).*STATES(:,34).*ALGEBRAIC(:,133)))./( ALGEBRAIC(:,158).*CONSTANTS(:,241).*CONSTANTS(:,244).*CONSTANTS(:,286).*( 2.00000.*CONSTANTS(:,245)+( CONSTANTS(:,245).*STATES(:,37))./CONSTANTS(:,241)+( STATES(:,37).*STATES(:,32).*STATES(:,2))./( CONSTANTS(:,241).*CONSTANTS(:,244).*CONSTANTS(:,286))+( CONSTANTS(:,245).*STATES(:,31).*ALGEBRAIC(:,133))./( CONSTANTS(:,242).*CONSTANTS(:,286))+( STATES(:,31).*STATES(:,34).*ALGEBRAIC(:,133))./( CONSTANTS(:,242).*CONSTANTS(:,243).*CONSTANTS(:,286))+( CONSTANTS(:,245).*STATES(:,31))./CONSTANTS(:,242)+( CONSTANTS(:,245).*STATES(:,37).*STATES(:,2))./( CONSTANTS(:,241).*CONSTANTS(:,286))+( CONSTANTS(:,245).*STATES(:,2))./CONSTANTS(:,286)+( CONSTANTS(:,245).*STATES(:,34).*ALGEBRAIC(:,133))./( CONSTANTS(:,243).*CONSTANTS(:,286))+( CONSTANTS(:,245).*ALGEBRAIC(:,133))./CONSTANTS(:,286)+( CONSTANTS(:,245).*STATES(:,32).*STATES(:,2))./( CONSTANTS(:,244).*CONSTANTS(:,286))));
    ALGEBRAIC(:,40) = (( STATES(:,16).*STATES(:,2))./CONSTANTS(:,361))./ALGEBRAIC(:,38);
    ALGEBRAIC(:,163) = STATES(:,8);
    ALGEBRAIC(:,171) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,361)+ALGEBRAIC(:,150)./CONSTANTS(:,362)+ALGEBRAIC(:,163)./CONSTANTS(:,363);
    ALGEBRAIC(:,173) = (( STATES(:,17).*ALGEBRAIC(:,133))./CONSTANTS(:,361))./ALGEBRAIC(:,171);
    ALGEBRAIC(:,175) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,371)+ALGEBRAIC(:,150)./CONSTANTS(:,372)+ALGEBRAIC(:,163)./CONSTANTS(:,373);
    ALGEBRAIC(:,176) =  CONSTANTS(:,37).*(( ALGEBRAIC(:,173).*STATES(:,30))./ALGEBRAIC(:,45) - ( ALGEBRAIC(:,40).*STATES(:,35))./ALGEBRAIC(:,175));
    ALGEBRAIC(:,174) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,367)+ALGEBRAIC(:,150)./CONSTANTS(:,368)+ALGEBRAIC(:,163)./CONSTANTS(:,369);
    ALGEBRAIC(:,178) =  CONSTANTS(:,39).*(( (STATES(:,39)./ALGEBRAIC(:,174)).*STATES(:,30))./ALGEBRAIC(:,45) - ( (STATES(:,27)./ALGEBRAIC(:,43)).*STATES(:,35))./ALGEBRAIC(:,175));
    ALGEBRAIC(:,169) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,352)+ALGEBRAIC(:,150)./CONSTANTS(:,354)+ALGEBRAIC(:,163)./CONSTANTS(:,356);
    ALGEBRAIC(:,177) =  CONSTANTS(:,40).*(( (STATES(:,35)./ALGEBRAIC(:,175)).*STATES(:,14))./ALGEBRAIC(:,34) - ( (STATES(:,30)./ALGEBRAIC(:,45)).*STATES(:,15))./ALGEBRAIC(:,169));
    ALGEBRAIC(:,113) = exp( - (CONSTANTS(:,291) -  CONSTANTS(:,8).*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275))./1.00000;
    ALGEBRAIC(:,143) = ( (( ALGEBRAIC(:,113).*power(ALGEBRAIC(:,133), CONSTANTS(:,8)))./power(STATES(:,2), CONSTANTS(:,8) - 1.00000)).*ALGEBRAIC(:,1))./( ALGEBRAIC(:,22).*ALGEBRAIC(:,34));
    ALGEBRAIC(:,144) =  CONSTANTS(:,47).*( ALGEBRAIC(:,143).*STATES(:,9).*STATES(:,14) - STATES(:,1));
    ALGEBRAIC(:,164) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,297)+ALGEBRAIC(:,150)./CONSTANTS(:,312)+ALGEBRAIC(:,163)./CONSTANTS(:,326);
    ALGEBRAIC(:,166) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,333)+ALGEBRAIC(:,150)./CONSTANTS(:,337)+ALGEBRAIC(:,163)./CONSTANTS(:,341);
    ALGEBRAIC(:,179) = piecewise({STATES(:,51)>CONSTANTS(:,250)&STATES(:,1)>CONSTANTS(:,250),  CONSTANTS(:,48).*((STATES(:,51)./ALGEBRAIC(:,166))./(STATES(:,51)./ALGEBRAIC(:,166)+ (STATES(:,52)./ALGEBRAIC(:,164)).*exp(( CONSTANTS(:,13).*CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275))) - (STATES(:,9)./ALGEBRAIC(:,22))./(STATES(:,9)./ALGEBRAIC(:,22)+ (STATES(:,1)./ALGEBRAIC(:,1)).*exp(( CONSTANTS(:,13).*(CONSTANTS(:,3) - 1.00000).*STATES(:,38))./CONSTANTS(:,275)))).*((STATES(:,51)./ALGEBRAIC(:,166))./(STATES(:,51)./ALGEBRAIC(:,166)+CONSTANTS(:,12))) }, 0.00000);
    ALGEBRAIC(:,145) = piecewise({abs(STATES(:,38))>CONSTANTS(:,251), ( CONSTANTS(:,52).*STATES(:,38).*( ALGEBRAIC(:,133).*exp(( CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)) - STATES(:,2)))./(exp(( CONSTANTS(:,3).*STATES(:,38))./CONSTANTS(:,275)) - 1.00000) }, ( CONSTANTS(:,52).*CONSTANTS(:,275).*(ALGEBRAIC(:,133) - STATES(:,2)))./CONSTANTS(:,3));
    ALGEBRAIC(:,35) = (( STATES(:,14).*STATES(:,2))./CONSTANTS(:,352))./ALGEBRAIC(:,34);
    ALGEBRAIC(:,170) = (( STATES(:,15).*ALGEBRAIC(:,133))./CONSTANTS(:,352))./ALGEBRAIC(:,169);
    ALGEBRAIC(:,180) = ( CONSTANTS(:,49).*( ALGEBRAIC(:,170).*ALGEBRAIC(:,133) -  ALGEBRAIC(:,35).*STATES(:,2)))./( CONSTANTS(:,50).*(1.00000+ALGEBRAIC(:,170)./CONSTANTS(:,50)).*(1.00000+ALGEBRAIC(:,35)./CONSTANTS(:,50)));
    ALGEBRAIC(:,168) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,345)+ALGEBRAIC(:,150)./CONSTANTS(:,348)+ALGEBRAIC(:,163)./CONSTANTS(:,350);
    ALGEBRAIC(:,181) = ( CONSTANTS(:,382).*ALGEBRAIC(:,164).*ALGEBRAIC(:,168))./( ALGEBRAIC(:,166).*ALGEBRAIC(:,166));
    ALGEBRAIC(:,183) = piecewise({CONSTANTS(:,1)==1.00000||CONSTANTS(:,1)==2.00000||CONSTANTS(:,1)==5.00000, 0.00000 },  CONSTANTS(:,262).*( ALGEBRAIC(:,181).*STATES(:,51).*STATES(:,51) -  STATES(:,54).*STATES(:,52)));
    ALGEBRAIC(:,190) =  - (( STATES(:,2).*STATES(:,2))./( CONSTANTS(:,135).*CONSTANTS(:,136).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,2).*STATES(:,1))./( CONSTANTS(:,297).*CONSTANTS(:,312).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,2).*STATES(:,9))./( CONSTANTS(:,333).*CONSTANTS(:,337).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,2).*STATES(:,11))./( CONSTANTS(:,345).*CONSTANTS(:,348).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,2).*STATES(:,12))./( CONSTANTS(:,297).*CONSTANTS(:,312).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,2).*STATES(:,13))./( CONSTANTS(:,333).*CONSTANTS(:,337).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,2).*STATES(:,14))./( CONSTANTS(:,352).*CONSTANTS(:,354).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,2).*STATES(:,22))./( CONSTANTS(:,111).*CONSTANTS(:,112).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,114))./( CONSTANTS(:,108).*CONSTANTS(:,109).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,2).*STATES(:,28))./( CONSTANTS(:,114).*CONSTANTS(:,115).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,119))./( CONSTANTS(:,117).*CONSTANTS(:,118).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,2).*STATES(:,18))./( CONSTANTS(:,99).*CONSTANTS(:,375).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,2).*STATES(:,23))./( CONSTANTS(:,90).*CONSTANTS(:,360).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,2).*STATES(:,21))./( CONSTANTS(:,120).*CONSTANTS(:,121).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,2).*STATES(:,16))./( CONSTANTS(:,361).*CONSTANTS(:,362).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,2).*STATES(:,25))./( CONSTANTS(:,364).*CONSTANTS(:,365).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,2).*STATES(:,26))./( CONSTANTS(:,123).*CONSTANTS(:,124).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,2).*STATES(:,24))./( CONSTANTS(:,366).*CONSTANTS(:,93).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,2).*STATES(:,19))./( CONSTANTS(:,358).*CONSTANTS(:,88).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,2).*STATES(:,27))./( CONSTANTS(:,367).*CONSTANTS(:,368).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,2).*STATES(:,29))./( CONSTANTS(:,370).*CONSTANTS(:,95).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,2).*STATES(:,30))./( CONSTANTS(:,371).*CONSTANTS(:,372).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,2).*STATES(:,32))./( CONSTANTS(:,377).*CONSTANTS(:,378).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,2).*STATES(:,31))./( CONSTANTS(:,379).*CONSTANTS(:,380).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,2).*STATES(:,4))./( CONSTANTS(:,141).*CONSTANTS(:,142).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,2).*STATES(:,3))./( CONSTANTS(:,138).*CONSTANTS(:,139).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,2).*STATES(:,57))./( CONSTANTS(:,126).*CONSTANTS(:,127).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,2).*CONSTANTS(:,292))./( CONSTANTS(:,129).*CONSTANTS(:,130).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,2).*CONSTANTS(:,311))./( CONSTANTS(:,132).*CONSTANTS(:,133).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,2).*STATES(:,20))./( CONSTANTS(:,374).*CONSTANTS(:,97).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,191) =  - (( STATES(:,2).*STATES(:,2))./( CONSTANTS(:,135).*CONSTANTS(:,137).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,2).*STATES(:,1))./( CONSTANTS(:,297).*CONSTANTS(:,326).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,2).*STATES(:,9))./( CONSTANTS(:,333).*CONSTANTS(:,341).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,2).*STATES(:,11))./( CONSTANTS(:,345).*CONSTANTS(:,350).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,2).*STATES(:,12))./( CONSTANTS(:,297).*CONSTANTS(:,326).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,2).*STATES(:,13))./( CONSTANTS(:,333).*CONSTANTS(:,341).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,2).*STATES(:,14))./( CONSTANTS(:,352).*CONSTANTS(:,356).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,2).*STATES(:,22))./( CONSTANTS(:,111).*CONSTANTS(:,113).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,114))./( CONSTANTS(:,108).*CONSTANTS(:,110).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,2).*STATES(:,28))./( CONSTANTS(:,114).*CONSTANTS(:,116).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,2).*ALGEBRAIC(:,119))./( CONSTANTS(:,117).*CONSTANTS(:,119).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,2).*STATES(:,18))./( CONSTANTS(:,99).*CONSTANTS(:,100).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,2).*STATES(:,23))./( CONSTANTS(:,90).*CONSTANTS(:,91).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,2).*STATES(:,21))./( CONSTANTS(:,120).*CONSTANTS(:,122).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,2).*STATES(:,16))./( CONSTANTS(:,361).*CONSTANTS(:,363).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,2).*STATES(:,25))./( CONSTANTS(:,364).*CONSTANTS(:,92).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,2).*STATES(:,26))./( CONSTANTS(:,123).*CONSTANTS(:,125).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,2).*STATES(:,24))./( CONSTANTS(:,366).*CONSTANTS(:,94).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,2).*STATES(:,19))./( CONSTANTS(:,358).*CONSTANTS(:,89).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,2).*STATES(:,27))./( CONSTANTS(:,367).*CONSTANTS(:,369).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,2).*STATES(:,29))./( CONSTANTS(:,370).*CONSTANTS(:,96).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,2).*STATES(:,30))./( CONSTANTS(:,371).*CONSTANTS(:,373).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,2).*STATES(:,32))./( CONSTANTS(:,377).*CONSTANTS(:,103).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,2).*STATES(:,31))./( CONSTANTS(:,379).*CONSTANTS(:,104).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,2).*STATES(:,4))./( CONSTANTS(:,141).*CONSTANTS(:,143).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,2).*STATES(:,3))./( CONSTANTS(:,138).*CONSTANTS(:,140).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,2).*STATES(:,57))./( CONSTANTS(:,126).*CONSTANTS(:,128).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,2).*CONSTANTS(:,292))./( CONSTANTS(:,129).*CONSTANTS(:,131).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,2).*CONSTANTS(:,311))./( CONSTANTS(:,132).*CONSTANTS(:,134).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,2).*STATES(:,20))./( CONSTANTS(:,374).*CONSTANTS(:,98).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,194) =  - (( STATES(:,3).*STATES(:,2))./( CONSTANTS(:,136).*CONSTANTS(:,137).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,3).*STATES(:,1))./( CONSTANTS(:,312).*CONSTANTS(:,326).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,3).*STATES(:,9))./( CONSTANTS(:,337).*CONSTANTS(:,341).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,3).*STATES(:,11))./( CONSTANTS(:,348).*CONSTANTS(:,350).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,3).*STATES(:,12))./( CONSTANTS(:,312).*CONSTANTS(:,326).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,3).*STATES(:,13))./( CONSTANTS(:,337).*CONSTANTS(:,341).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,3).*STATES(:,14))./( CONSTANTS(:,354).*CONSTANTS(:,356).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,3).*STATES(:,22))./( CONSTANTS(:,112).*CONSTANTS(:,113).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,114))./( CONSTANTS(:,109).*CONSTANTS(:,110).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,3).*STATES(:,28))./( CONSTANTS(:,115).*CONSTANTS(:,116).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,119))./( CONSTANTS(:,118).*CONSTANTS(:,119).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,3).*STATES(:,18))./( CONSTANTS(:,375).*CONSTANTS(:,100).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,3).*STATES(:,23))./( CONSTANTS(:,360).*CONSTANTS(:,91).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,3).*STATES(:,21))./( CONSTANTS(:,121).*CONSTANTS(:,122).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,3).*STATES(:,16))./( CONSTANTS(:,362).*CONSTANTS(:,363).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,3).*STATES(:,25))./( CONSTANTS(:,365).*CONSTANTS(:,92).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,3).*STATES(:,26))./( CONSTANTS(:,124).*CONSTANTS(:,125).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,3).*STATES(:,24))./( CONSTANTS(:,93).*CONSTANTS(:,94).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,3).*STATES(:,19))./( CONSTANTS(:,88).*CONSTANTS(:,89).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,3).*STATES(:,27))./( CONSTANTS(:,368).*CONSTANTS(:,369).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,3).*STATES(:,29))./( CONSTANTS(:,95).*CONSTANTS(:,96).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,3).*STATES(:,30))./( CONSTANTS(:,372).*CONSTANTS(:,373).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,3).*STATES(:,32))./( CONSTANTS(:,378).*CONSTANTS(:,103).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,3).*STATES(:,31))./( CONSTANTS(:,380).*CONSTANTS(:,104).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,3).*STATES(:,4))./( CONSTANTS(:,142).*CONSTANTS(:,143).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,3).*STATES(:,3))./( CONSTANTS(:,139).*CONSTANTS(:,140).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,3).*STATES(:,57))./( CONSTANTS(:,127).*CONSTANTS(:,128).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,3).*CONSTANTS(:,292))./( CONSTANTS(:,130).*CONSTANTS(:,131).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,3).*CONSTANTS(:,311))./( CONSTANTS(:,133).*CONSTANTS(:,134).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,3).*STATES(:,20))./( CONSTANTS(:,97).*CONSTANTS(:,98).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,197) =  - (( STATES(:,4).*STATES(:,2))./( CONSTANTS(:,137).*CONSTANTS(:,136).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,4).*STATES(:,1))./( CONSTANTS(:,326).*CONSTANTS(:,312).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,4).*STATES(:,9))./( CONSTANTS(:,341).*CONSTANTS(:,337).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,4).*STATES(:,11))./( CONSTANTS(:,350).*CONSTANTS(:,348).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,4).*STATES(:,12))./( CONSTANTS(:,326).*CONSTANTS(:,312).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,4).*STATES(:,13))./( CONSTANTS(:,341).*CONSTANTS(:,337).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,4).*STATES(:,14))./( CONSTANTS(:,356).*CONSTANTS(:,354).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,4).*STATES(:,22))./( CONSTANTS(:,113).*CONSTANTS(:,112).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,114))./( CONSTANTS(:,110).*CONSTANTS(:,109).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,4).*STATES(:,28))./( CONSTANTS(:,116).*CONSTANTS(:,115).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,119))./( CONSTANTS(:,119).*CONSTANTS(:,118).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,4).*STATES(:,18))./( CONSTANTS(:,100).*CONSTANTS(:,375).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,4).*STATES(:,23))./( CONSTANTS(:,91).*CONSTANTS(:,360).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,4).*STATES(:,21))./( CONSTANTS(:,122).*CONSTANTS(:,121).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,4).*STATES(:,16))./( CONSTANTS(:,363).*CONSTANTS(:,362).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,4).*STATES(:,25))./( CONSTANTS(:,92).*CONSTANTS(:,365).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,4).*STATES(:,26))./( CONSTANTS(:,125).*CONSTANTS(:,124).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,4).*STATES(:,24))./( CONSTANTS(:,94).*CONSTANTS(:,93).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,4).*STATES(:,19))./( CONSTANTS(:,89).*CONSTANTS(:,88).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,4).*STATES(:,27))./( CONSTANTS(:,369).*CONSTANTS(:,368).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,4).*STATES(:,29))./( CONSTANTS(:,96).*CONSTANTS(:,95).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,4).*STATES(:,30))./( CONSTANTS(:,373).*CONSTANTS(:,372).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,4).*STATES(:,32))./( CONSTANTS(:,103).*CONSTANTS(:,378).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,4).*STATES(:,31))./( CONSTANTS(:,104).*CONSTANTS(:,380).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,4).*STATES(:,4))./( CONSTANTS(:,143).*CONSTANTS(:,142).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,4).*STATES(:,3))./( CONSTANTS(:,140).*CONSTANTS(:,139).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,4).*STATES(:,57))./( CONSTANTS(:,128).*CONSTANTS(:,127).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,4).*CONSTANTS(:,292))./( CONSTANTS(:,131).*CONSTANTS(:,130).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,4).*CONSTANTS(:,311))./( CONSTANTS(:,134).*CONSTANTS(:,133).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,4).*STATES(:,20))./( CONSTANTS(:,98).*CONSTANTS(:,97).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,198) = ( (1.00000+STATES(:,2)./CONSTANTS(:,135)+STATES(:,3)./CONSTANTS(:,136)).*STATES(:,2))./( CONSTANTS(:,137).*power(CONSTANTS(:,283), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,3)./CONSTANTS(:,312)).*STATES(:,1))./( CONSTANTS(:,326).*power(ALGEBRAIC(:,1), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,3)./CONSTANTS(:,337)).*STATES(:,9))./( CONSTANTS(:,341).*power(ALGEBRAIC(:,22), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,345)+STATES(:,3)./CONSTANTS(:,348)).*STATES(:,11))./( CONSTANTS(:,350).*power(ALGEBRAIC(:,27), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,3)./CONSTANTS(:,312)).*STATES(:,12))./( CONSTANTS(:,326).*power(ALGEBRAIC(:,30), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,3)./CONSTANTS(:,337)).*STATES(:,13))./( CONSTANTS(:,341).*power(ALGEBRAIC(:,32), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,352)+STATES(:,3)./CONSTANTS(:,354)).*STATES(:,14))./( CONSTANTS(:,356).*power(ALGEBRAIC(:,34), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,111)+STATES(:,3)./CONSTANTS(:,112)).*STATES(:,22))./( CONSTANTS(:,113).*power(CONSTANTS(:,278), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,108)+STATES(:,3)./CONSTANTS(:,109)).*ALGEBRAIC(:,114))./( CONSTANTS(:,110).*power(CONSTANTS(:,277), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,114)+STATES(:,3)./CONSTANTS(:,115)).*STATES(:,28))./( CONSTANTS(:,116).*power(CONSTANTS(:,305), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,117)+STATES(:,3)./CONSTANTS(:,118)).*ALGEBRAIC(:,119))./( CONSTANTS(:,119).*power(CONSTANTS(:,321), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,99)+STATES(:,3)./CONSTANTS(:,375)).*STATES(:,18))./( CONSTANTS(:,100).*power(ALGEBRAIC(:,47), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,90)+STATES(:,3)./CONSTANTS(:,360)).*STATES(:,23))./( CONSTANTS(:,91).*power(ALGEBRAIC(:,37), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,120)+STATES(:,3)./CONSTANTS(:,121)).*STATES(:,21))./( CONSTANTS(:,122).*power(CONSTANTS(:,279), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,361)+STATES(:,3)./CONSTANTS(:,362)).*STATES(:,16))./( CONSTANTS(:,363).*power(ALGEBRAIC(:,38), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,364)+STATES(:,3)./CONSTANTS(:,365)).*STATES(:,25))./( CONSTANTS(:,92).*power(ALGEBRAIC(:,41), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,123)+STATES(:,3)./CONSTANTS(:,124)).*STATES(:,26))./( CONSTANTS(:,125).*power(CONSTANTS(:,307), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,366)+STATES(:,3)./CONSTANTS(:,93)).*STATES(:,24))./( CONSTANTS(:,94).*power(ALGEBRAIC(:,42), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,358)+STATES(:,3)./CONSTANTS(:,88)).*STATES(:,19))./( CONSTANTS(:,89).*power(ALGEBRAIC(:,36), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,367)+STATES(:,3)./CONSTANTS(:,368)).*STATES(:,27))./( CONSTANTS(:,369).*power(ALGEBRAIC(:,43), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,370)+STATES(:,3)./CONSTANTS(:,95)).*STATES(:,29))./( CONSTANTS(:,96).*power(ALGEBRAIC(:,44), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,371)+STATES(:,3)./CONSTANTS(:,372)).*STATES(:,30))./( CONSTANTS(:,373).*power(ALGEBRAIC(:,45), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,377)+STATES(:,3)./CONSTANTS(:,378)).*STATES(:,32))./( CONSTANTS(:,103).*power(ALGEBRAIC(:,49), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,379)+STATES(:,3)./CONSTANTS(:,380)).*STATES(:,31))./( CONSTANTS(:,104).*power(ALGEBRAIC(:,50), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,141)+STATES(:,3)./CONSTANTS(:,142)).*STATES(:,4))./( CONSTANTS(:,143).*power(CONSTANTS(:,285), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,138)+STATES(:,3)./CONSTANTS(:,139)).*STATES(:,3))./( CONSTANTS(:,140).*power(CONSTANTS(:,284), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,126)+STATES(:,3)./CONSTANTS(:,127)).*STATES(:,57))./( CONSTANTS(:,128).*power(CONSTANTS(:,280), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,129)+STATES(:,3)./CONSTANTS(:,130)).*CONSTANTS(:,292))./( CONSTANTS(:,131).*power(CONSTANTS(:,281), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,132)+STATES(:,3)./CONSTANTS(:,133)).*CONSTANTS(:,311))./( CONSTANTS(:,134).*power(CONSTANTS(:,282), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,374)+STATES(:,3)./CONSTANTS(:,97)).*STATES(:,20))./( CONSTANTS(:,98).*power(ALGEBRAIC(:,46), 2.00000));
    ALGEBRAIC(:,206) = 1.00000+ALGEBRAIC(:,198);
    ALGEBRAIC(:,195) = ( (1.00000+STATES(:,2)./CONSTANTS(:,135)+STATES(:,4)./CONSTANTS(:,137)).*STATES(:,2))./( CONSTANTS(:,136).*power(CONSTANTS(:,283), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,1))./( CONSTANTS(:,312).*power(ALGEBRAIC(:,1), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,9))./( CONSTANTS(:,337).*power(ALGEBRAIC(:,22), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,345)+STATES(:,4)./CONSTANTS(:,350)).*STATES(:,11))./( CONSTANTS(:,348).*power(ALGEBRAIC(:,27), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,297)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,12))./( CONSTANTS(:,312).*power(ALGEBRAIC(:,30), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,333)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,13))./( CONSTANTS(:,337).*power(ALGEBRAIC(:,32), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,352)+STATES(:,4)./CONSTANTS(:,356)).*STATES(:,14))./( CONSTANTS(:,354).*power(ALGEBRAIC(:,34), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,111)+STATES(:,4)./CONSTANTS(:,113)).*STATES(:,22))./( CONSTANTS(:,112).*power(CONSTANTS(:,278), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,108)+STATES(:,4)./CONSTANTS(:,110)).*ALGEBRAIC(:,114))./( CONSTANTS(:,109).*power(CONSTANTS(:,277), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,114)+STATES(:,4)./CONSTANTS(:,116)).*STATES(:,28))./( CONSTANTS(:,115).*power(CONSTANTS(:,305), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,117)+STATES(:,4)./CONSTANTS(:,119)).*ALGEBRAIC(:,119))./( CONSTANTS(:,118).*power(CONSTANTS(:,321), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,99)+STATES(:,4)./CONSTANTS(:,100)).*STATES(:,18))./( CONSTANTS(:,375).*power(ALGEBRAIC(:,47), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,90)+STATES(:,4)./CONSTANTS(:,91)).*STATES(:,23))./( CONSTANTS(:,360).*power(ALGEBRAIC(:,37), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,120)+STATES(:,4)./CONSTANTS(:,122)).*STATES(:,21))./( CONSTANTS(:,121).*power(CONSTANTS(:,279), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,361)+STATES(:,4)./CONSTANTS(:,363)).*STATES(:,16))./( CONSTANTS(:,362).*power(ALGEBRAIC(:,38), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,364)+STATES(:,4)./CONSTANTS(:,92)).*STATES(:,25))./( CONSTANTS(:,365).*power(ALGEBRAIC(:,41), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,123)+STATES(:,4)./CONSTANTS(:,125)).*STATES(:,26))./( CONSTANTS(:,124).*power(CONSTANTS(:,307), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,366)+STATES(:,4)./CONSTANTS(:,94)).*STATES(:,24))./( CONSTANTS(:,93).*power(ALGEBRAIC(:,42), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,358)+STATES(:,4)./CONSTANTS(:,89)).*STATES(:,19))./( CONSTANTS(:,88).*power(ALGEBRAIC(:,36), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,367)+STATES(:,4)./CONSTANTS(:,369)).*STATES(:,27))./( CONSTANTS(:,368).*power(ALGEBRAIC(:,43), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,370)+STATES(:,4)./CONSTANTS(:,96)).*STATES(:,29))./( CONSTANTS(:,95).*power(ALGEBRAIC(:,44), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,371)+STATES(:,4)./CONSTANTS(:,373)).*STATES(:,30))./( CONSTANTS(:,372).*power(ALGEBRAIC(:,45), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,377)+STATES(:,4)./CONSTANTS(:,103)).*STATES(:,32))./( CONSTANTS(:,378).*power(ALGEBRAIC(:,49), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,379)+STATES(:,4)./CONSTANTS(:,104)).*STATES(:,31))./( CONSTANTS(:,380).*power(ALGEBRAIC(:,50), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,141)+STATES(:,4)./CONSTANTS(:,143)).*STATES(:,4))./( CONSTANTS(:,142).*power(CONSTANTS(:,285), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,138)+STATES(:,4)./CONSTANTS(:,140)).*STATES(:,3))./( CONSTANTS(:,139).*power(CONSTANTS(:,284), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,126)+STATES(:,4)./CONSTANTS(:,128)).*STATES(:,57))./( CONSTANTS(:,127).*power(CONSTANTS(:,280), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,129)+STATES(:,4)./CONSTANTS(:,131)).*CONSTANTS(:,292))./( CONSTANTS(:,130).*power(CONSTANTS(:,281), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,132)+STATES(:,4)./CONSTANTS(:,134)).*CONSTANTS(:,311))./( CONSTANTS(:,133).*power(CONSTANTS(:,282), 2.00000))+( (1.00000+STATES(:,2)./CONSTANTS(:,374)+STATES(:,4)./CONSTANTS(:,98)).*STATES(:,20))./( CONSTANTS(:,97).*power(ALGEBRAIC(:,46), 2.00000));
    ALGEBRAIC(:,205) = 1.00000+ALGEBRAIC(:,195);
    ALGEBRAIC(:,192) = ( (1.00000+STATES(:,3)./CONSTANTS(:,136)+STATES(:,4)./CONSTANTS(:,137)).*STATES(:,2))./( CONSTANTS(:,135).*power(CONSTANTS(:,283), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,312)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,1))./( CONSTANTS(:,297).*power(ALGEBRAIC(:,1), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,337)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,9))./( CONSTANTS(:,333).*power(ALGEBRAIC(:,22), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,348)+STATES(:,4)./CONSTANTS(:,350)).*STATES(:,11))./( CONSTANTS(:,345).*power(ALGEBRAIC(:,27), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,312)+STATES(:,4)./CONSTANTS(:,326)).*STATES(:,12))./( CONSTANTS(:,297).*power(ALGEBRAIC(:,30), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,337)+STATES(:,4)./CONSTANTS(:,341)).*STATES(:,13))./( CONSTANTS(:,333).*power(ALGEBRAIC(:,32), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,354)+STATES(:,4)./CONSTANTS(:,356)).*STATES(:,14))./( CONSTANTS(:,352).*power(ALGEBRAIC(:,34), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,112)+STATES(:,4)./CONSTANTS(:,113)).*STATES(:,22))./( CONSTANTS(:,111).*power(CONSTANTS(:,278), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,109)+STATES(:,4)./CONSTANTS(:,110)).*ALGEBRAIC(:,114))./( CONSTANTS(:,108).*power(CONSTANTS(:,277), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,115)+STATES(:,4)./CONSTANTS(:,116)).*STATES(:,28))./( CONSTANTS(:,114).*power(CONSTANTS(:,305), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,118)+STATES(:,4)./CONSTANTS(:,119)).*ALGEBRAIC(:,119))./( CONSTANTS(:,117).*power(CONSTANTS(:,321), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,375)+STATES(:,4)./CONSTANTS(:,100)).*STATES(:,18))./( CONSTANTS(:,99).*power(ALGEBRAIC(:,47), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,360)+STATES(:,4)./CONSTANTS(:,91)).*STATES(:,23))./( CONSTANTS(:,90).*power(ALGEBRAIC(:,37), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,121)+STATES(:,4)./CONSTANTS(:,122)).*STATES(:,21))./( CONSTANTS(:,120).*power(CONSTANTS(:,279), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,362)+STATES(:,4)./CONSTANTS(:,363)).*STATES(:,16))./( CONSTANTS(:,361).*power(ALGEBRAIC(:,38), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,365)+STATES(:,4)./CONSTANTS(:,92)).*STATES(:,25))./( CONSTANTS(:,364).*power(ALGEBRAIC(:,41), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,124)+STATES(:,4)./CONSTANTS(:,125)).*STATES(:,26))./( CONSTANTS(:,123).*power(CONSTANTS(:,307), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,93)+STATES(:,4)./CONSTANTS(:,94)).*STATES(:,24))./( CONSTANTS(:,366).*power(ALGEBRAIC(:,42), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,88)+STATES(:,4)./CONSTANTS(:,89)).*STATES(:,19))./( CONSTANTS(:,358).*power(ALGEBRAIC(:,36), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,368)+STATES(:,4)./CONSTANTS(:,369)).*STATES(:,27))./( CONSTANTS(:,367).*power(ALGEBRAIC(:,43), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,95)+STATES(:,4)./CONSTANTS(:,96)).*STATES(:,29))./( CONSTANTS(:,370).*power(ALGEBRAIC(:,44), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,372)+STATES(:,4)./CONSTANTS(:,373)).*STATES(:,30))./( CONSTANTS(:,371).*power(ALGEBRAIC(:,45), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,378)+STATES(:,4)./CONSTANTS(:,103)).*STATES(:,32))./( CONSTANTS(:,377).*power(ALGEBRAIC(:,49), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,380)+STATES(:,4)./CONSTANTS(:,104)).*STATES(:,31))./( CONSTANTS(:,379).*power(ALGEBRAIC(:,50), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,142)+STATES(:,4)./CONSTANTS(:,143)).*STATES(:,4))./( CONSTANTS(:,141).*power(CONSTANTS(:,285), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,139)+STATES(:,4)./CONSTANTS(:,140)).*STATES(:,3))./( CONSTANTS(:,138).*power(CONSTANTS(:,284), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,127)+STATES(:,4)./CONSTANTS(:,128)).*STATES(:,57))./( CONSTANTS(:,126).*power(CONSTANTS(:,280), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,130)+STATES(:,4)./CONSTANTS(:,131)).*CONSTANTS(:,292))./( CONSTANTS(:,129).*power(CONSTANTS(:,281), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,133)+STATES(:,4)./CONSTANTS(:,134)).*CONSTANTS(:,311))./( CONSTANTS(:,132).*power(CONSTANTS(:,282), 2.00000))+( (1.00000+STATES(:,3)./CONSTANTS(:,97)+STATES(:,4)./CONSTANTS(:,98)).*STATES(:,20))./( CONSTANTS(:,374).*power(ALGEBRAIC(:,46), 2.00000));
    ALGEBRAIC(:,204) = 1.00000+ALGEBRAIC(:,192)+CONSTANTS(:,272)./( CONSTANTS(:,273).*power(1.00000+STATES(:,2)./CONSTANTS(:,273), 2.00000));
    ALGEBRAIC(:,196) =  - (( STATES(:,4).*STATES(:,2))./( CONSTANTS(:,137).*CONSTANTS(:,135).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,4).*STATES(:,1))./( CONSTANTS(:,326).*CONSTANTS(:,297).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,4).*STATES(:,9))./( CONSTANTS(:,341).*CONSTANTS(:,333).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,4).*STATES(:,11))./( CONSTANTS(:,350).*CONSTANTS(:,345).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,4).*STATES(:,12))./( CONSTANTS(:,326).*CONSTANTS(:,297).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,4).*STATES(:,13))./( CONSTANTS(:,341).*CONSTANTS(:,333).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,4).*STATES(:,14))./( CONSTANTS(:,356).*CONSTANTS(:,352).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,4).*STATES(:,22))./( CONSTANTS(:,113).*CONSTANTS(:,111).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,114))./( CONSTANTS(:,110).*CONSTANTS(:,108).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,4).*STATES(:,28))./( CONSTANTS(:,116).*CONSTANTS(:,114).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,4).*ALGEBRAIC(:,119))./( CONSTANTS(:,119).*CONSTANTS(:,117).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,4).*STATES(:,18))./( CONSTANTS(:,100).*CONSTANTS(:,99).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,4).*STATES(:,23))./( CONSTANTS(:,91).*CONSTANTS(:,90).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,4).*STATES(:,21))./( CONSTANTS(:,122).*CONSTANTS(:,120).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,4).*STATES(:,16))./( CONSTANTS(:,363).*CONSTANTS(:,361).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,4).*STATES(:,25))./( CONSTANTS(:,92).*CONSTANTS(:,364).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,4).*STATES(:,26))./( CONSTANTS(:,125).*CONSTANTS(:,123).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,4).*STATES(:,24))./( CONSTANTS(:,94).*CONSTANTS(:,366).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,4).*STATES(:,19))./( CONSTANTS(:,89).*CONSTANTS(:,358).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,4).*STATES(:,27))./( CONSTANTS(:,369).*CONSTANTS(:,367).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,4).*STATES(:,29))./( CONSTANTS(:,96).*CONSTANTS(:,370).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,4).*STATES(:,30))./( CONSTANTS(:,373).*CONSTANTS(:,371).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,4).*STATES(:,32))./( CONSTANTS(:,103).*CONSTANTS(:,377).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,4).*STATES(:,31))./( CONSTANTS(:,104).*CONSTANTS(:,379).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,4).*STATES(:,4))./( CONSTANTS(:,143).*CONSTANTS(:,141).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,4).*STATES(:,3))./( CONSTANTS(:,140).*CONSTANTS(:,138).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,4).*STATES(:,57))./( CONSTANTS(:,128).*CONSTANTS(:,126).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,4).*CONSTANTS(:,292))./( CONSTANTS(:,131).*CONSTANTS(:,129).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,4).*CONSTANTS(:,311))./( CONSTANTS(:,134).*CONSTANTS(:,132).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,4).*STATES(:,20))./( CONSTANTS(:,98).*CONSTANTS(:,374).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,193) =  - (( STATES(:,3).*STATES(:,2))./( CONSTANTS(:,136).*CONSTANTS(:,135).*power(CONSTANTS(:,283), 2.00000))+( STATES(:,3).*STATES(:,1))./( CONSTANTS(:,312).*CONSTANTS(:,297).*power(ALGEBRAIC(:,1), 2.00000))+( STATES(:,3).*STATES(:,9))./( CONSTANTS(:,337).*CONSTANTS(:,333).*power(ALGEBRAIC(:,22), 2.00000))+( STATES(:,3).*STATES(:,11))./( CONSTANTS(:,348).*CONSTANTS(:,345).*power(ALGEBRAIC(:,27), 2.00000))+( STATES(:,3).*STATES(:,12))./( CONSTANTS(:,312).*CONSTANTS(:,297).*power(ALGEBRAIC(:,30), 2.00000))+( STATES(:,3).*STATES(:,13))./( CONSTANTS(:,337).*CONSTANTS(:,333).*power(ALGEBRAIC(:,32), 2.00000))+( STATES(:,3).*STATES(:,14))./( CONSTANTS(:,354).*CONSTANTS(:,352).*power(ALGEBRAIC(:,34), 2.00000))+( STATES(:,3).*STATES(:,22))./( CONSTANTS(:,112).*CONSTANTS(:,111).*power(CONSTANTS(:,278), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,114))./( CONSTANTS(:,109).*CONSTANTS(:,108).*power(CONSTANTS(:,277), 2.00000))+( STATES(:,3).*STATES(:,28))./( CONSTANTS(:,115).*CONSTANTS(:,114).*power(CONSTANTS(:,305), 2.00000))+( STATES(:,3).*ALGEBRAIC(:,119))./( CONSTANTS(:,118).*CONSTANTS(:,117).*power(CONSTANTS(:,321), 2.00000))+( STATES(:,3).*STATES(:,18))./( CONSTANTS(:,375).*CONSTANTS(:,99).*power(ALGEBRAIC(:,47), 2.00000))+( STATES(:,3).*STATES(:,23))./( CONSTANTS(:,360).*CONSTANTS(:,90).*power(ALGEBRAIC(:,37), 2.00000))+( STATES(:,3).*STATES(:,21))./( CONSTANTS(:,121).*CONSTANTS(:,120).*power(CONSTANTS(:,279), 2.00000))+( STATES(:,3).*STATES(:,16))./( CONSTANTS(:,362).*CONSTANTS(:,361).*power(ALGEBRAIC(:,38), 2.00000))+( STATES(:,3).*STATES(:,25))./( CONSTANTS(:,365).*CONSTANTS(:,364).*power(ALGEBRAIC(:,41), 2.00000))+( STATES(:,3).*STATES(:,26))./( CONSTANTS(:,124).*CONSTANTS(:,123).*power(CONSTANTS(:,307), 2.00000))+( STATES(:,3).*STATES(:,24))./( CONSTANTS(:,93).*CONSTANTS(:,366).*power(ALGEBRAIC(:,42), 2.00000))+( STATES(:,3).*STATES(:,19))./( CONSTANTS(:,88).*CONSTANTS(:,358).*power(ALGEBRAIC(:,36), 2.00000))+( STATES(:,3).*STATES(:,27))./( CONSTANTS(:,368).*CONSTANTS(:,367).*power(ALGEBRAIC(:,43), 2.00000))+( STATES(:,3).*STATES(:,29))./( CONSTANTS(:,95).*CONSTANTS(:,370).*power(ALGEBRAIC(:,44), 2.00000))+( STATES(:,3).*STATES(:,30))./( CONSTANTS(:,372).*CONSTANTS(:,371).*power(ALGEBRAIC(:,45), 2.00000))+( STATES(:,3).*STATES(:,32))./( CONSTANTS(:,378).*CONSTANTS(:,377).*power(ALGEBRAIC(:,49), 2.00000))+( STATES(:,3).*STATES(:,31))./( CONSTANTS(:,380).*CONSTANTS(:,379).*power(ALGEBRAIC(:,50), 2.00000))+( STATES(:,3).*STATES(:,4))./( CONSTANTS(:,142).*CONSTANTS(:,141).*power(CONSTANTS(:,285), 2.00000))+( STATES(:,3).*STATES(:,3))./( CONSTANTS(:,139).*CONSTANTS(:,138).*power(CONSTANTS(:,284), 2.00000))+( STATES(:,3).*STATES(:,57))./( CONSTANTS(:,127).*CONSTANTS(:,126).*power(CONSTANTS(:,280), 2.00000))+( STATES(:,3).*CONSTANTS(:,292))./( CONSTANTS(:,130).*CONSTANTS(:,129).*power(CONSTANTS(:,281), 2.00000))+( STATES(:,3).*CONSTANTS(:,311))./( CONSTANTS(:,133).*CONSTANTS(:,132).*power(CONSTANTS(:,282), 2.00000))+( STATES(:,3).*STATES(:,20))./( CONSTANTS(:,97).*CONSTANTS(:,374).*power(ALGEBRAIC(:,46), 2.00000)));
    ALGEBRAIC(:,207) = ((( ALGEBRAIC(:,204).*ALGEBRAIC(:,197).*ALGEBRAIC(:,194)+ ALGEBRAIC(:,206).*ALGEBRAIC(:,190).*ALGEBRAIC(:,193)+ ALGEBRAIC(:,205).*ALGEBRAIC(:,191).*ALGEBRAIC(:,196)) -  ALGEBRAIC(:,205).*ALGEBRAIC(:,206).*ALGEBRAIC(:,204)) -  ALGEBRAIC(:,191).*ALGEBRAIC(:,197).*ALGEBRAIC(:,193)) -  ALGEBRAIC(:,190).*ALGEBRAIC(:,194).*ALGEBRAIC(:,196);
    ALGEBRAIC(:,146) = ((ALGEBRAIC(:,115)+ALGEBRAIC(:,116)+ALGEBRAIC(:,117)+ALGEBRAIC(:,118)) - ALGEBRAIC(:,140))./CONSTANTS(:,4);
    ALGEBRAIC(:,147) =  - ALGEBRAIC(:,146);
    ALGEBRAIC(:,148) = ((ALGEBRAIC(:,120)+ALGEBRAIC(:,140)) - ALGEBRAIC(:,123))./CONSTANTS(:,4);
    ALGEBRAIC(:,149) =  - ALGEBRAIC(:,148);
    ALGEBRAIC(:,184) = ((ALGEBRAIC(:,82)+ALGEBRAIC(:,144)) - ALGEBRAIC(:,179))./CONSTANTS(:,4);
    ALGEBRAIC(:,185) = (( - ALGEBRAIC(:,82) - ALGEBRAIC(:,144))+ALGEBRAIC(:,179))./CONSTANTS(:,4);
    ALGEBRAIC(:,124) = (ALGEBRAIC(:,67) - ALGEBRAIC(:,82))./CONSTANTS(:,4);
    ALGEBRAIC(:,125) = ( - ALGEBRAIC(:,67)+ALGEBRAIC(:,82))./CONSTANTS(:,4);
    ALGEBRAIC(:,186) = ((( - ALGEBRAIC(:,67) - ALGEBRAIC(:,144))+ALGEBRAIC(:,180)) - ALGEBRAIC(:,177))./CONSTANTS(:,4);
    ALGEBRAIC(:,160) = ( - ALGEBRAIC(:,115)+ALGEBRAIC(:,156))./CONSTANTS(:,4);
    ALGEBRAIC(:,126) = (ALGEBRAIC(:,115) - ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,187) = ( - ALGEBRAIC(:,60)+ALGEBRAIC(:,176)+ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,127) = (ALGEBRAIC(:,60) - ALGEBRAIC(:,116))./CONSTANTS(:,4);
    ALGEBRAIC(:,128) = (((ALGEBRAIC(:,116) - ALGEBRAIC(:,117)) - ALGEBRAIC(:,87))+ALGEBRAIC(:,88))./CONSTANTS(:,4);
    ALGEBRAIC(:,129) = (ALGEBRAIC(:,117) - ALGEBRAIC(:,67))./CONSTANTS(:,4);
    ALGEBRAIC(:,130) = (( - ALGEBRAIC(:,115) - ALGEBRAIC(:,117))+ALGEBRAIC(:,67)+ALGEBRAIC(:,57))./CONSTANTS(:,4);
    ALGEBRAIC(:,188) = ((ALGEBRAIC(:,67) - ALGEBRAIC(:,120))+ALGEBRAIC(:,178))./CONSTANTS(:,4);
    ALGEBRAIC(:,131) = (ALGEBRAIC(:,120) - ALGEBRAIC(:,74))./CONSTANTS(:,4);
    ALGEBRAIC(:,189) = (((((ALGEBRAIC(:,74) - ALGEBRAIC(:,118))+ALGEBRAIC(:,177)) - ALGEBRAIC(:,88)) - ALGEBRAIC(:,176)) - ALGEBRAIC(:,178))./CONSTANTS(:,4);
    ALGEBRAIC(:,132) = ( - ALGEBRAIC(:,57)+ALGEBRAIC(:,118)+ALGEBRAIC(:,87))./CONSTANTS(:,4);
    ALGEBRAIC(:,161) = ((ALGEBRAIC(:,87)+ALGEBRAIC(:,157)) - ALGEBRAIC(:,159))./CONSTANTS(:,4);
    ALGEBRAIC(:,162) = ( - ALGEBRAIC(:,87)+ALGEBRAIC(:,159))./CONSTANTS(:,4);
    ALGEBRAIC(:,182) =  CONSTANTS(:,51).*( ALGEBRAIC(:,163).*STATES(:,2) -  STATES(:,4).*ALGEBRAIC(:,133));
    ALGEBRAIC(:,199) = ((((((((( - ALGEBRAIC(:,115)+ 2.00000.*ALGEBRAIC(:,57)) - ALGEBRAIC(:,117))+ALGEBRAIC(:,67)+ALGEBRAIC(:,118)+ALGEBRAIC(:,156)+ALGEBRAIC(:,157)+ALGEBRAIC(:,176)) - ALGEBRAIC(:,159)) -  5.00000.*ALGEBRAIC(:,140)) -  2.00000.*ALGEBRAIC(:,123)) -  4.00000.*ALGEBRAIC(:,142))+ (CONSTANTS(:,8) - 1.00000).*ALGEBRAIC(:,144)+ 2.00000.*ALGEBRAIC(:,180)+ALGEBRAIC(:,145)) - ALGEBRAIC(:,182))./CONSTANTS(:,4);
    ALGEBRAIC(:,201) =  - (( STATES(:,2).*CONSTANTS(:,293))./( CONSTANTS(:,135).*CONSTANTS(:,283))+( STATES(:,2).*ALGEBRAIC(:,184))./( CONSTANTS(:,297).*ALGEBRAIC(:,1))+( STATES(:,2).*ALGEBRAIC(:,185))./( CONSTANTS(:,333).*ALGEBRAIC(:,22))+( STATES(:,2).*CONSTANTS(:,263))./( CONSTANTS(:,345).*ALGEBRAIC(:,27))+( STATES(:,2).*ALGEBRAIC(:,124))./( CONSTANTS(:,297).*ALGEBRAIC(:,30))+( STATES(:,2).*ALGEBRAIC(:,125))./( CONSTANTS(:,333).*ALGEBRAIC(:,32))+( STATES(:,2).*ALGEBRAIC(:,186))./( CONSTANTS(:,352).*ALGEBRAIC(:,34))+( STATES(:,2).*ALGEBRAIC(:,146))./( CONSTANTS(:,111).*CONSTANTS(:,278))+( STATES(:,2).*ALGEBRAIC(:,147))./( CONSTANTS(:,108).*CONSTANTS(:,277))+( STATES(:,2).*ALGEBRAIC(:,148))./( CONSTANTS(:,114).*CONSTANTS(:,305))+( STATES(:,2).*ALGEBRAIC(:,149))./( CONSTANTS(:,117).*CONSTANTS(:,321))+( STATES(:,2).*ALGEBRAIC(:,160))./( CONSTANTS(:,99).*ALGEBRAIC(:,47))+( STATES(:,2).*ALGEBRAIC(:,132))./( CONSTANTS(:,90).*ALGEBRAIC(:,37))+( STATES(:,2).*ALGEBRAIC(:,126))./( CONSTANTS(:,120).*CONSTANTS(:,279))+( STATES(:,2).*ALGEBRAIC(:,187))./( CONSTANTS(:,361).*ALGEBRAIC(:,38))+( STATES(:,2).*ALGEBRAIC(:,127))./( CONSTANTS(:,364).*ALGEBRAIC(:,41))+( STATES(:,2).*ALGEBRAIC(:,128))./( CONSTANTS(:,123).*CONSTANTS(:,307))+( STATES(:,2).*ALGEBRAIC(:,129))./( CONSTANTS(:,366).*ALGEBRAIC(:,42))+( STATES(:,2).*ALGEBRAIC(:,130))./( CONSTANTS(:,358).*ALGEBRAIC(:,36))+( STATES(:,2).*ALGEBRAIC(:,188))./( CONSTANTS(:,367).*ALGEBRAIC(:,43))+( STATES(:,2).*ALGEBRAIC(:,131))./( CONSTANTS(:,370).*ALGEBRAIC(:,44))+( STATES(:,2).*ALGEBRAIC(:,189))./( CONSTANTS(:,371).*ALGEBRAIC(:,45))+( STATES(:,2).*ALGEBRAIC(:,161))./( CONSTANTS(:,377).*ALGEBRAIC(:,49))+( STATES(:,2).*ALGEBRAIC(:,162))./( CONSTANTS(:,379).*ALGEBRAIC(:,50))+( STATES(:,2).*CONSTANTS(:,295))./( CONSTANTS(:,141).*CONSTANTS(:,285))+( STATES(:,2).*CONSTANTS(:,294))./( CONSTANTS(:,138).*CONSTANTS(:,284))+( STATES(:,2).*CONSTANTS(:,264))./( CONSTANTS(:,126).*CONSTANTS(:,280))+( STATES(:,2).*CONSTANTS(:,271))./( CONSTANTS(:,129).*CONSTANTS(:,281))+( STATES(:,2).*CONSTANTS(:,270))./( CONSTANTS(:,132).*CONSTANTS(:,282))+( STATES(:,2).*CONSTANTS(:,265))./( CONSTANTS(:,374).*ALGEBRAIC(:,46)))+ALGEBRAIC(:,199);
    ALGEBRAIC(:,202) =  - (( STATES(:,3).*CONSTANTS(:,293))./( CONSTANTS(:,136).*CONSTANTS(:,283))+( STATES(:,3).*ALGEBRAIC(:,184))./( CONSTANTS(:,312).*ALGEBRAIC(:,1))+( STATES(:,3).*ALGEBRAIC(:,185))./( CONSTANTS(:,337).*ALGEBRAIC(:,22))+( STATES(:,3).*CONSTANTS(:,263))./( CONSTANTS(:,348).*ALGEBRAIC(:,27))+( STATES(:,3).*ALGEBRAIC(:,124))./( CONSTANTS(:,312).*ALGEBRAIC(:,30))+( STATES(:,3).*ALGEBRAIC(:,125))./( CONSTANTS(:,337).*ALGEBRAIC(:,32))+( STATES(:,3).*ALGEBRAIC(:,186))./( CONSTANTS(:,354).*ALGEBRAIC(:,34))+( STATES(:,3).*ALGEBRAIC(:,146))./( CONSTANTS(:,112).*CONSTANTS(:,278))+( STATES(:,3).*ALGEBRAIC(:,147))./( CONSTANTS(:,109).*CONSTANTS(:,277))+( STATES(:,3).*ALGEBRAIC(:,148))./( CONSTANTS(:,115).*CONSTANTS(:,305))+( STATES(:,3).*ALGEBRAIC(:,149))./( CONSTANTS(:,118).*CONSTANTS(:,321))+( STATES(:,3).*ALGEBRAIC(:,160))./( CONSTANTS(:,375).*ALGEBRAIC(:,47))+( STATES(:,3).*ALGEBRAIC(:,132))./( CONSTANTS(:,360).*ALGEBRAIC(:,37))+( STATES(:,3).*ALGEBRAIC(:,126))./( CONSTANTS(:,121).*CONSTANTS(:,279))+( STATES(:,3).*ALGEBRAIC(:,187))./( CONSTANTS(:,362).*ALGEBRAIC(:,38))+( STATES(:,3).*ALGEBRAIC(:,127))./( CONSTANTS(:,365).*ALGEBRAIC(:,41))+( STATES(:,3).*ALGEBRAIC(:,128))./( CONSTANTS(:,124).*CONSTANTS(:,307))+( STATES(:,3).*ALGEBRAIC(:,129))./( CONSTANTS(:,93).*ALGEBRAIC(:,42))+( STATES(:,3).*ALGEBRAIC(:,130))./( CONSTANTS(:,88).*ALGEBRAIC(:,36))+( STATES(:,3).*ALGEBRAIC(:,188))./( CONSTANTS(:,368).*ALGEBRAIC(:,43))+( STATES(:,3).*ALGEBRAIC(:,131))./( CONSTANTS(:,95).*ALGEBRAIC(:,44))+( STATES(:,3).*ALGEBRAIC(:,189))./( CONSTANTS(:,372).*ALGEBRAIC(:,45))+( STATES(:,3).*ALGEBRAIC(:,161))./( CONSTANTS(:,378).*ALGEBRAIC(:,49))+( STATES(:,3).*ALGEBRAIC(:,162))./( CONSTANTS(:,380).*ALGEBRAIC(:,50))+( STATES(:,3).*CONSTANTS(:,295))./( CONSTANTS(:,142).*CONSTANTS(:,285))+( STATES(:,3).*CONSTANTS(:,294))./( CONSTANTS(:,139).*CONSTANTS(:,284))+( STATES(:,3).*CONSTANTS(:,264))./( CONSTANTS(:,127).*CONSTANTS(:,280))+( STATES(:,3).*CONSTANTS(:,271))./( CONSTANTS(:,130).*CONSTANTS(:,281))+( STATES(:,3).*CONSTANTS(:,270))./( CONSTANTS(:,133).*CONSTANTS(:,282))+( STATES(:,3).*CONSTANTS(:,265))./( CONSTANTS(:,97).*ALGEBRAIC(:,46)))+CONSTANTS(:,296);
    ALGEBRAIC(:,200) = ALGEBRAIC(:,182)./CONSTANTS(:,4);
    ALGEBRAIC(:,203) =  - (( STATES(:,4).*CONSTANTS(:,293))./( CONSTANTS(:,137).*CONSTANTS(:,283))+( STATES(:,4).*ALGEBRAIC(:,184))./( CONSTANTS(:,326).*ALGEBRAIC(:,1))+( STATES(:,4).*ALGEBRAIC(:,185))./( CONSTANTS(:,341).*ALGEBRAIC(:,22))+( STATES(:,4).*CONSTANTS(:,263))./( CONSTANTS(:,350).*ALGEBRAIC(:,27))+( STATES(:,4).*ALGEBRAIC(:,124))./( CONSTANTS(:,326).*ALGEBRAIC(:,30))+( STATES(:,4).*ALGEBRAIC(:,125))./( CONSTANTS(:,341).*ALGEBRAIC(:,32))+( STATES(:,4).*ALGEBRAIC(:,186))./( CONSTANTS(:,356).*ALGEBRAIC(:,34))+( STATES(:,4).*ALGEBRAIC(:,146))./( CONSTANTS(:,113).*CONSTANTS(:,278))+( STATES(:,4).*ALGEBRAIC(:,147))./( CONSTANTS(:,110).*CONSTANTS(:,277))+( STATES(:,4).*ALGEBRAIC(:,148))./( CONSTANTS(:,116).*CONSTANTS(:,305))+( STATES(:,4).*ALGEBRAIC(:,149))./( CONSTANTS(:,119).*CONSTANTS(:,321))+( STATES(:,4).*ALGEBRAIC(:,160))./( CONSTANTS(:,100).*ALGEBRAIC(:,47))+( STATES(:,4).*ALGEBRAIC(:,132))./( CONSTANTS(:,91).*ALGEBRAIC(:,37))+( STATES(:,4).*ALGEBRAIC(:,126))./( CONSTANTS(:,122).*CONSTANTS(:,279))+( STATES(:,4).*ALGEBRAIC(:,187))./( CONSTANTS(:,363).*ALGEBRAIC(:,38))+( STATES(:,4).*ALGEBRAIC(:,127))./( CONSTANTS(:,92).*ALGEBRAIC(:,41))+( STATES(:,4).*ALGEBRAIC(:,128))./( CONSTANTS(:,125).*CONSTANTS(:,307))+( STATES(:,4).*ALGEBRAIC(:,129))./( CONSTANTS(:,94).*ALGEBRAIC(:,42))+( STATES(:,4).*ALGEBRAIC(:,130))./( CONSTANTS(:,89).*ALGEBRAIC(:,36))+( STATES(:,4).*ALGEBRAIC(:,188))./( CONSTANTS(:,369).*ALGEBRAIC(:,43))+( STATES(:,4).*ALGEBRAIC(:,131))./( CONSTANTS(:,96).*ALGEBRAIC(:,44))+( STATES(:,4).*ALGEBRAIC(:,189))./( CONSTANTS(:,373).*ALGEBRAIC(:,45))+( STATES(:,4).*ALGEBRAIC(:,161))./( CONSTANTS(:,103).*ALGEBRAIC(:,49))+( STATES(:,4).*ALGEBRAIC(:,162))./( CONSTANTS(:,104).*ALGEBRAIC(:,50))+( STATES(:,4).*CONSTANTS(:,295))./( CONSTANTS(:,143).*CONSTANTS(:,285))+( STATES(:,4).*CONSTANTS(:,294))./( CONSTANTS(:,140).*CONSTANTS(:,284))+( STATES(:,4).*CONSTANTS(:,264))./( CONSTANTS(:,128).*CONSTANTS(:,280))+( STATES(:,4).*CONSTANTS(:,271))./( CONSTANTS(:,131).*CONSTANTS(:,281))+( STATES(:,4).*CONSTANTS(:,270))./( CONSTANTS(:,134).*CONSTANTS(:,282))+( STATES(:,4).*CONSTANTS(:,265))./( CONSTANTS(:,98).*ALGEBRAIC(:,46)))+ALGEBRAIC(:,200);
    ALGEBRAIC(:,2) = 1.00000+STATES(:,6)./CONSTANTS(:,352)+STATES(:,7)./CONSTANTS(:,354)+STATES(:,8)./CONSTANTS(:,356);
    ALGEBRAIC(:,3) = 1.00000+STATES(:,6)./CONSTANTS(:,358);
    ALGEBRAIC(:,4) = 1.00000+STATES(:,7)./CONSTANTS(:,360);
    ALGEBRAIC(:,5) = 1.00000+STATES(:,6)./CONSTANTS(:,361)+STATES(:,7)./CONSTANTS(:,362)+STATES(:,8)./CONSTANTS(:,363);
    ALGEBRAIC(:,6) = 1.00000+STATES(:,6)./CONSTANTS(:,364)+STATES(:,7)./CONSTANTS(:,365);
    ALGEBRAIC(:,7) = 1.00000+STATES(:,6)./CONSTANTS(:,366);
    ALGEBRAIC(:,8) = 1.00000+STATES(:,6)./CONSTANTS(:,367)+STATES(:,7)./CONSTANTS(:,368)+STATES(:,8)./CONSTANTS(:,369);
    ALGEBRAIC(:,9) = 1.00000+STATES(:,6)./CONSTANTS(:,370);
    ALGEBRAIC(:,10) = 1.00000+STATES(:,6)./CONSTANTS(:,371)+STATES(:,7)./CONSTANTS(:,372)+STATES(:,8)./CONSTANTS(:,373);
    ALGEBRAIC(:,11) = 1.00000+STATES(:,6)./CONSTANTS(:,374);
    ALGEBRAIC(:,12) = 1.00000+STATES(:,7)./CONSTANTS(:,375);
    ALGEBRAIC(:,13) = 1.00000+STATES(:,2)./CONSTANTS(:,376);
    ALGEBRAIC(:,14) = 1.00000+STATES(:,6)./CONSTANTS(:,377)+STATES(:,7)./CONSTANTS(:,378);
    ALGEBRAIC(:,15) = 1.00000+STATES(:,6)./CONSTANTS(:,379)+STATES(:,7)./CONSTANTS(:,380);
    ALGEBRAIC(:,20) = ALGEBRAIC(:,18);
    ALGEBRAIC(:,25) = ALGEBRAIC(:,23);
    ALGEBRAIC(:,134) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,358);
    ALGEBRAIC(:,135) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,366);
    ALGEBRAIC(:,136) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,370);
    ALGEBRAIC(:,137) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,374);
    ALGEBRAIC(:,138) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,376);
    ALGEBRAIC(:,151) = 1.00000+ALGEBRAIC(:,150)./CONSTANTS(:,360);
    ALGEBRAIC(:,152) = 1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,364)+ALGEBRAIC(:,150)./CONSTANTS(:,365);
    ALGEBRAIC(:,165) = ALGEBRAIC(:,164);
    ALGEBRAIC(:,167) = ALGEBRAIC(:,166);
    ALGEBRAIC(:,172) = ( STATES(:,17).*(1.00000+ALGEBRAIC(:,133)./CONSTANTS(:,361)))./ALGEBRAIC(:,171);
end

% Compute result of a piecewise function
function x = piecewise(cases, default)
    set = [0];
    for i = 1:2:length(cases)
        if (length(cases{i+1}) == 1)
            x(cases{i} & ~set,:) = cases{i+1};
        else
            x(cases{i} & ~set,:) = cases{i+1}(cases{i} & ~set);
        end
        set = set | cases{i};
        if(set), break, end
    end
    if (length(default) == 1)
        x(~set,:) = default;
    else
        x(~set,:) = default(~set);
    end
end

% Pad out or shorten strings to a set length
function strout = strpad(strin)
    req_length = 160;
    insize = size(strin,2);
    if insize > req_length
        strout = strin(1:req_length);
    else
        strout = [strin, blanks(req_length - insize)];
    end
end