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 =179;
end
% There are a total of 17 entries in each of the rate and state variable arrays.
% There are a total of 230 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 (day)');
    LEGEND_ALGEBRAIC(:,11) = strpad('U_Pd_PdPs in component Q_Pd (mole_per_day)');
    LEGEND_STATES(:,1) = strpad('Q_Pd in component Q_Pd (mole)');
    LEGEND_ALGEBRAIC(:,1) = strpad('C_Pd in component Q_Pd (molar)');
    LEGEND_ALGEBRAIC(:,9) = strpad('U_Pd_PdEx in component Q_Pd (mole_per_day)');
    LEGEND_CONSTANTS(:,156) = strpad('P_Pd_InPd in component Q_Pd (mole_per_day)');
    LEGEND_ALGEBRAIC(:,10) = strpad('v_PdPs in component Q_Pd (mole_per_day)');
    LEGEND_CONSTANTS(:,1) = strpad('v_PdPs_star in component Q_Pd (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,175) = strpad('M_Pd_PdPs in component Q_Pd (molar)');
    LEGEND_CONSTANTS(:,2) = strpad('M_Pd_PdPs_star in component Q_Pd (molar)');
    LEGEND_CONSTANTS(:,3) = strpad('Y_Pd_InPd in component Q_Pd (molar)');
    LEGEND_CONSTANTS(:,4) = strpad('D_Pd in component Q_Pd (litre_per_day)');
    LEGEND_CONSTANTS(:,155) = strpad('k_PdEx in component Q_Pd (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,5) = strpad('T_Pd in component Q_Pd (day)');
    LEGEND_CONSTANTS(:,6) = strpad('T_Pd_star in component Q_Pd (day)');
    LEGEND_STATES(:,2) = strpad('Q_Ma in component Q_Ma (gram)');
    LEGEND_STATES(:,3) = strpad('Q_Mc in component Q_Mc (gram)');
    LEGEND_CONSTANTS(:,7) = strpad('V_Ru in component model_parameters (litre)');
    LEGEND_CONSTANTS(:,8) = strpad('k_SoEx in component model_parameters (first_order_rate_constant)');
    LEGEND_ALGEBRAIC(:,12) = strpad('C_Ps in component Q_Ps (molar)');
    LEGEND_ALGEBRAIC(:,49) = strpad('U_McPs_PsAm in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,37) = strpad('U_Ps_PsMa in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,50) = strpad('U_Ps_PsMc in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,36) = strpad('U_MaPs_PsAm in component Q_Ps (mole_per_day)');
    LEGEND_CONSTANTS(:,9) = strpad('M_Ha_McMa in component Q_Ps (molar)');
    LEGEND_STATES(:,4) = strpad('Q_Ps in component Q_Ps (mole)');
    LEGEND_CONSTANTS(:,177) = strpad('P_Ps_InPs in component Q_Ps (mole_per_day)');
    LEGEND_CONSTANTS(:,180) = strpad('P_Ps_SaPs in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,13) = strpad('P_Ps_PdPs in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,63) = strpad('P_Ps_MaMd in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,79) = strpad('P_Ps_McPs in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,14) = strpad('U_Ps_PsEx in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,77) = strpad('U_Mc_McPs in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,15) = strpad('v_Ma_PsAm in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,16) = strpad('v_Mc_PsAm in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,17) = strpad('v_PsMa in component Q_Ps (mole_per_day)');
    LEGEND_ALGEBRAIC(:,18) = strpad('v_PsMc in component Q_Ps (mole_per_day)');
    LEGEND_CONSTANTS(:,10) = strpad('v_PsAm_star in component Q_Ps (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,11) = strpad('v_PsMa_star in component Q_Ps (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,12) = strpad('v_PsMc_star in component Q_Ps (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,157) = strpad('k_PsEx in component Q_Ps (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,13) = strpad('Y_Ps_InPs in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,14) = strpad('Y_Ps_SaPs in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,15) = strpad('Y_Ps_PdPs in component Q_Ps (dimensionless)');
    LEGEND_CONSTANTS(:,16) = strpad('Y_Ps_MaMd in component Q_Ps (mole_per_gram)');
    LEGEND_CONSTANTS(:,17) = strpad('Y_Ps_McPs in component Q_Ps (dimensionless)');
    LEGEND_CONSTANTS(:,18) = strpad('D_Ps in component Q_Ps (litre_per_day)');
    LEGEND_CONSTANTS(:,19) = strpad('M_Ps_PsAm in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,20) = strpad('M_Ps_PsMa in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,21) = strpad('M_Ps_PsMc in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,22) = strpad('M_Ha_PsMa in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,23) = strpad('M_Hc_PsMc in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,24) = strpad('J_Ha_PsAm in component Q_Ps (molar)');
    LEGEND_CONSTANTS(:,25) = strpad('J_Hc_PsAm in component Q_Ps (molar)');
    LEGEND_ALGEBRAIC(:,76) = strpad('U_Mc_McEg in component Q_Mc (gram_per_day)');
    LEGEND_ALGEBRAIC(:,62) = strpad('U_Ma_MaMd in component Q_Ma (gram_per_day)');
    LEGEND_ALGEBRAIC(:,55) = strpad('U_Hc_PsMc in component Q_Hc (mole_per_day)');
    LEGEND_ALGEBRAIC(:,83) = strpad('U_Ha_McMa in component Q_Ha (mole_per_day)');
    LEGEND_ALGEBRAIC(:,35) = strpad('C_Ha in component Q_Ha (molar)');
    LEGEND_ALGEBRAIC(:,48) = strpad('C_Hc in component Q_Hc (molar)');
    LEGEND_CONSTANTS(:,26) = strpad('k_FlEx in component model_parameters (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,27) = strpad('J_Ha_McAm in component model_parameters (molar)');
    LEGEND_STATES(:,5) = strpad('Q_Pu in component Q_Pu (mole)');
    LEGEND_ALGEBRAIC(:,2) = strpad('C_Pu in component Q_Pu (molar)');
    LEGEND_CONSTANTS(:,158) = strpad('P_Pu_InPu in component Q_Pu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,3) = strpad('U_Pu_PuEx in component Q_Pu (mole_per_day)');
    LEGEND_CONSTANTS(:,28) = strpad('Y_Pu_InPu in component Q_Pu (molar)');
    LEGEND_CONSTANTS(:,29) = strpad('D_Pu in component Q_Pu (litre_per_day)');
    LEGEND_CONSTANTS(:,159) = strpad('k_PuEx in component Q_Pu (first_order_rate_constant)');
    LEGEND_ALGEBRAIC(:,19) = strpad('C_Am in component Q_Am (molar)');
    LEGEND_ALGEBRAIC(:,52) = strpad('U_Am_AmMc in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,39) = strpad('U_Am_AmMa in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,78) = strpad('U_Mc_McAm in component Q_Am (mole_per_day)');
    LEGEND_STATES(:,6) = strpad('Q_Am in component Q_Am (mole)');
    LEGEND_CONSTANTS(:,181) = strpad('P_Am_InAm in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,23) = strpad('P_Am_UeAm in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,38) = strpad('P_MaAm_PsAm in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,51) = strpad('P_McAm_PsAm in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,80) = strpad('P_Am_McAm in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,20) = strpad('U_Am_AmAb in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,21) = strpad('U_Am_AmEx in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,22) = strpad('v_UeAm in component Q_Am (mole_per_day)');
    LEGEND_CONSTANTS(:,182) = strpad('v_AmAb in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,24) = strpad('v_AmMa in component Q_Am (mole_per_day)');
    LEGEND_ALGEBRAIC(:,25) = strpad('v_AmMc in component Q_Am (mole_per_day)');
    LEGEND_CONSTANTS(:,30) = strpad('v_UeAm_star in component Q_Am (mole_per_litre_day)');
    LEGEND_CONSTANTS(:,31) = strpad('v_AmAb_star in component Q_Am (mole_per_litre_day)');
    LEGEND_CONSTANTS(:,32) = strpad('v_AmMa_star in component Q_Am (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,33) = strpad('v_AmMc_star in component Q_Am (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,160) = strpad('k_AmEx in component Q_Am (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,34) = strpad('Y_Am_InAm in component Q_Am (molar)');
    LEGEND_CONSTANTS(:,35) = strpad('Y_Am_UeAm in component Q_Am (dimensionless)');
    LEGEND_CONSTANTS(:,36) = strpad('Y_Am_PsAm in component Q_Am (dimensionless)');
    LEGEND_CONSTANTS(:,37) = strpad('Y_Am_McAm in component Q_Am (dimensionless)');
    LEGEND_CONSTANTS(:,38) = strpad('M_Am_AmMa in component Q_Am (molar)');
    LEGEND_CONSTANTS(:,39) = strpad('M_Am_AmMc in component Q_Am (molar)');
    LEGEND_CONSTANTS(:,40) = strpad('M_Am_AmAb in component Q_Am (molar)');
    LEGEND_CONSTANTS(:,41) = strpad('M_pH_AmAb in component Q_Am (dimensionless)');
    LEGEND_CONSTANTS(:,42) = strpad('phi_pH_AmAb in component Q_Am (dimensionless)');
    LEGEND_CONSTANTS(:,43) = strpad('M_Ha_AmMa in component Q_Am (molar)');
    LEGEND_CONSTANTS(:,44) = strpad('M_Hc_AmMc in component Q_Am (molar)');
    LEGEND_CONSTANTS(:,45) = strpad('J_Am_UeAm in component Q_Am (molar)');
    LEGEND_CONSTANTS(:,46) = strpad('D_Am in component Q_Am (litre_per_day)');
    LEGEND_CONSTANTS(:,47) = strpad('pH in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,48) = strpad('f_Po_Ma in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,49) = strpad('f_Ni_Fe in component model_parameters (dimensionless)');
    LEGEND_STATES(:,7) = strpad('Q_Li in component Q_Li (mole)');
    LEGEND_ALGEBRAIC(:,4) = strpad('C_Li in component Q_Li (molar)');
    LEGEND_CONSTANTS(:,183) = strpad('P_Li_InLi in component Q_Li (mole_per_day)');
    LEGEND_ALGEBRAIC(:,64) = strpad('P_Li_MaMd in component Q_Li (mole_per_day)');
    LEGEND_ALGEBRAIC(:,26) = strpad('U_Li_LiEx in component Q_Li (mole_per_day)');
    LEGEND_CONSTANTS(:,50) = strpad('Y_Li_InLi in component Q_Li (molar)');
    LEGEND_CONSTANTS(:,51) = strpad('Y_Li_MaMd in component Q_Li (mole_per_gram)');
    LEGEND_CONSTANTS(:,52) = strpad('D_Li in component Q_Li (litre_per_day)');
    LEGEND_CONSTANTS(:,161) = strpad('k_LiEx in component Q_Li (first_order_rate_constant)');
    LEGEND_ALGEBRAIC(:,30) = strpad('U_Sd_SdHa in component Q_Sd (mole_per_day)');
    LEGEND_STATES(:,8) = strpad('Q_Sd in component Q_Sd (mole)');
    LEGEND_ALGEBRAIC(:,27) = strpad('C_Sd in component Q_Sd (molar)');
    LEGEND_CONSTANTS(:,184) = strpad('P_Sd_InSd in component Q_Sd (mole_per_day)');
    LEGEND_ALGEBRAIC(:,28) = strpad('U_Sd_SdEx in component Q_Sd (mole_per_day)');
    LEGEND_CONSTANTS(:,53) = strpad('Y_Sd_InSd in component Q_Sd (molar)');
    LEGEND_CONSTANTS(:,54) = strpad('D_Sd in component Q_Sd (litre_per_day)');
    LEGEND_CONSTANTS(:,162) = strpad('k_SdEx in component Q_Sd (first_order_rate_constant)');
    LEGEND_ALGEBRAIC(:,29) = strpad('v_SdHa in component Q_Sd (mole_per_day)');
    LEGEND_CONSTANTS(:,55) = strpad('v_SdHa_star in component Q_Sd (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,185) = strpad('M_Sd_SdHa in component Q_Sd (molar)');
    LEGEND_CONSTANTS(:,56) = strpad('M_Sd_SdHa_star in component Q_Sd (molar)');
    LEGEND_CONSTANTS(:,57) = strpad('T_Sd in component Q_Sd (day)');
    LEGEND_CONSTANTS(:,58) = strpad('T_Sd_star in component Q_Sd (day)');
    LEGEND_ALGEBRAIC(:,34) = strpad('U_Fd_FdHc in component Q_Fd (mole_per_day)');
    LEGEND_STATES(:,9) = strpad('Q_Fd in component Q_Fd (mole)');
    LEGEND_ALGEBRAIC(:,31) = strpad('C_Fd in component Q_Fd (molar)');
    LEGEND_CONSTANTS(:,186) = strpad('P_Fd_InFd in component Q_Fd (mole_per_day)');
    LEGEND_ALGEBRAIC(:,32) = strpad('U_Fd_FdEx in component Q_Fd (mole_per_day)');
    LEGEND_CONSTANTS(:,59) = strpad('Y_Fd_InFd in component Q_Fd (molar)');
    LEGEND_CONSTANTS(:,60) = strpad('D_Fd in component Q_Fd (litre_per_day)');
    LEGEND_CONSTANTS(:,163) = strpad('k_FdEx in component Q_Fd (first_order_rate_constant)');
    LEGEND_ALGEBRAIC(:,33) = strpad('v_FdHc in component Q_Fd (mole_per_day)');
    LEGEND_CONSTANTS(:,61) = strpad('v_FdHc_star in component Q_Fd (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,187) = strpad('M_Fd_FdHc in component Q_Fd (molar)');
    LEGEND_CONSTANTS(:,62) = strpad('M_Fd_FdHc_star in component Q_Fd (molar)');
    LEGEND_CONSTANTS(:,63) = strpad('T_Fd in component Q_Fd (day)');
    LEGEND_CONSTANTS(:,64) = strpad('T_Fd_star in component Q_Fd (day)');
    LEGEND_CONSTANTS(:,65) = strpad('phi_pH_FdHc in component Q_Fd (dimensionless)');
    LEGEND_CONSTANTS(:,66) = strpad('M_pH_FdHc in component Q_Fd (dimensionless)');
    LEGEND_CONSTANTS(:,67) = strpad('pm in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,68) = strpad('f in component model_parameters (dimensionless)');
    LEGEND_STATES(:,10) = strpad('Q_Fu in component Q_Fu (mole)');
    LEGEND_ALGEBRAIC(:,5) = strpad('C_Fu in component Q_Fu (molar)');
    LEGEND_CONSTANTS(:,164) = strpad('P_Fu_InFu in component Q_Fu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,6) = strpad('U_Fu_FuEx in component Q_Fu (mole_per_day)');
    LEGEND_CONSTANTS(:,69) = strpad('Y_Fu_InFu in component Q_Fu (molar)');
    LEGEND_CONSTANTS(:,70) = strpad('D_Fu in component Q_Fu (litre_per_day)');
    LEGEND_CONSTANTS(:,165) = strpad('k_FuEx in component Q_Fu (first_order_rate_constant)');
    LEGEND_ALGEBRAIC(:,45) = strpad('U_Ha_HaAs in component Q_Ha (mole_per_day)');
    LEGEND_ALGEBRAIC(:,41) = strpad('U_Ha_AmMa in component Q_Ha (mole_per_day)');
    LEGEND_ALGEBRAIC(:,42) = strpad('U_Ha_PsMa in component Q_Ha (mole_per_day)');
    LEGEND_ALGEBRAIC(:,47) = strpad('U_Ha_HaVa in component Q_Ha (mole_per_day)');
    LEGEND_CONSTANTS(:,188) = strpad('P_Ha_WrHa in component Q_Ha (mole_per_day)');
    LEGEND_CONSTANTS(:,190) = strpad('P_Ha_LaHa in component Q_Ha (mole_per_day)');
    LEGEND_CONSTANTS(:,189) = strpad('P_Ha_SrHa in component Q_Ha (mole_per_day)');
    LEGEND_ALGEBRAIC(:,40) = strpad('P_Ha_SdHa in component Q_Ha (mole_per_day)');
    LEGEND_ALGEBRAIC(:,68) = strpad('P_Ha_MaMd in component Q_Ha (mole_per_day)');
    LEGEND_STATES(:,11) = strpad('Q_Ha in component Q_Ha (mole)');
    LEGEND_ALGEBRAIC(:,43) = strpad('U_Ha_HaEx in component Q_Ha (mole_per_day)');
    LEGEND_CONSTANTS(:,71) = strpad('R_Ha_AmMa in component Q_Ha (dimensionless)');
    LEGEND_CONSTANTS(:,72) = strpad('R_Ha_PsMa in component Q_Ha (dimensionless)');
    LEGEND_CONSTANTS(:,73) = strpad('R_Ha_McMa in component Q_Ha (mole_per_gram)');
    LEGEND_CONSTANTS(:,74) = strpad('M_Ha_HaAs in component Q_Ha (molar)');
    LEGEND_CONSTANTS(:,75) = strpad('M_Ha_HaVa in component Q_Ha (molar)');
    LEGEND_ALGEBRAIC(:,44) = strpad('v_HaAs in component Q_Ha (mole_per_day)');
    LEGEND_ALGEBRAIC(:,46) = strpad('v_HaVa in component Q_Ha (mole_per_day)');
    LEGEND_CONSTANTS(:,76) = strpad('v_HaAs_star in component Q_Ha (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,77) = strpad('v_HaVa_star in component Q_Ha (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,166) = strpad('k_HaEx in component Q_Ha (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,78) = strpad('Y_Ha_WrHa in component Q_Ha (molar)');
    LEGEND_CONSTANTS(:,79) = strpad('Y_Ha_SrHa in component Q_Ha (molar)');
    LEGEND_CONSTANTS(:,80) = strpad('Y_Ha_LaHa in component Q_Ha (molar)');
    LEGEND_CONSTANTS(:,81) = strpad('Y_Ha_SdHa in component Q_Ha (dimensionless)');
    LEGEND_CONSTANTS(:,82) = strpad('Y_Ha_MaMd in component Q_Ha (mole_per_gram)');
    LEGEND_CONSTANTS(:,83) = strpad('D_Wr in component Q_Ha (litre_per_day)');
    LEGEND_CONSTANTS(:,84) = strpad('D_Sr in component Q_Ha (litre_per_day)');
    LEGEND_CONSTANTS(:,85) = strpad('D_La in component Q_Ha (litre_per_day)');
    LEGEND_CONSTANTS(:,86) = strpad('J_Am_HaVa in component Q_Ha (molar)');
    LEGEND_CONSTANTS(:,87) = strpad('J_Ps_HaVa in component Q_Ha (molar)');
    LEGEND_ALGEBRAIC(:,67) = strpad('U_As_MaMd in component Q_As (gram_per_day)');
    LEGEND_ALGEBRAIC(:,81) = strpad('U_Mc_McMa in component Q_Ma (gram_per_day)');
    LEGEND_ALGEBRAIC(:,54) = strpad('U_Hc_AmMc in component Q_Hc (mole_per_day)');
    LEGEND_ALGEBRAIC(:,58) = strpad('U_Hc_HcVa in component Q_Hc (mole_per_day)');
    LEGEND_STATES(:,12) = strpad('Q_Hc in component Q_Hc (mole)');
    LEGEND_ALGEBRAIC(:,53) = strpad('P_Hc_FdHc in component Q_Hc (mole_per_day)');
    LEGEND_ALGEBRAIC(:,56) = strpad('U_Hc_HcEx in component Q_Hc (mole_per_day)');
    LEGEND_CONSTANTS(:,88) = strpad('R_Hc_AmMc in component Q_Hc (dimensionless)');
    LEGEND_CONSTANTS(:,89) = strpad('R_Hc_PsMc in component Q_Hc (dimensionless)');
    LEGEND_ALGEBRAIC(:,57) = strpad('v_HcVa in component Q_Hc (mole_per_day)');
    LEGEND_CONSTANTS(:,90) = strpad('v_HcVa_star in component Q_Hc (mole_per_gram_day)');
    LEGEND_CONSTANTS(:,167) = strpad('k_HcEx in component Q_Hc (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,91) = strpad('Y_Hc_FdHc in component Q_Hc (dimensionless)');
    LEGEND_CONSTANTS(:,92) = strpad('M_Hc_HcVa in component Q_Hc (molar)');
    LEGEND_CONSTANTS(:,93) = strpad('J_Am_HcVa in component Q_Hc (molar)');
    LEGEND_CONSTANTS(:,94) = strpad('J_Ps_HcVa in component Q_Hc (molar)');
    LEGEND_ALGEBRAIC(:,59) = strpad('C_Ma in component Q_Ma (gram_per_litre)');
    LEGEND_ALGEBRAIC(:,60) = strpad('P_Ma_AmMa in component Q_Ma (gram_per_day)');
    LEGEND_ALGEBRAIC(:,61) = strpad('P_Ma_PsMa in component Q_Ma (gram_per_day)');
    LEGEND_ALGEBRAIC(:,82) = strpad('P_Ma_McMa in component Q_Ma (gram_per_day)');
    LEGEND_ALGEBRAIC(:,65) = strpad('U_Ma_PoEx in component Q_Ma (gram_per_day)');
    LEGEND_ALGEBRAIC(:,7) = strpad('U_Ma_MaEx in component Q_Ma (gram_per_day)');
    LEGEND_CONSTANTS(:,178) = strpad('k_MaMd in component Q_Ma (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,95) = strpad('Y_Ma_AmMa in component Q_Ma (gram_per_mole)');
    LEGEND_CONSTANTS(:,96) = strpad('Y_Ma_PsMa in component Q_Ma (gram_per_mole)');
    LEGEND_CONSTANTS(:,97) = strpad('Y_Ma_McMa in component Q_Ma (dimensionless)');
    LEGEND_CONSTANTS(:,168) = strpad('k_MaEx in component model_parameters (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,176) = strpad('k_PoEx in component model_parameters (first_order_rate_constant)');
    LEGEND_STATES(:,13) = strpad('Q_As in component Q_As (gram)');
    LEGEND_ALGEBRAIC(:,8) = strpad('C_As in component Q_As (gram_per_litre)');
    LEGEND_ALGEBRAIC(:,66) = strpad('P_As_HaAs in component Q_As (gram_per_day)');
    LEGEND_ALGEBRAIC(:,71) = strpad('U_As_PoEx in component Q_As (gram_per_day)');
    LEGEND_ALGEBRAIC(:,69) = strpad('U_As_AsEx in component Q_As (gram_per_day)');
    LEGEND_CONSTANTS(:,169) = strpad('k_AsEx in component Q_As (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,179) = strpad('k_AsMd in component Q_As (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,98) = strpad('Y_As_HaAs in component Q_As (gram_per_mole)');
    LEGEND_ALGEBRAIC(:,70) = strpad('C_Mc in component Q_Mc (gram_per_litre)');
    LEGEND_ALGEBRAIC(:,72) = strpad('P_Mc_AmMc in component Q_Mc (gram_per_day)');
    LEGEND_ALGEBRAIC(:,73) = strpad('P_Mc_PsMc in component Q_Mc (gram_per_day)');
    LEGEND_ALGEBRAIC(:,74) = strpad('U_Mc_McEx in component Q_Mc (gram_per_day)');
    LEGEND_CONSTANTS(:,170) = strpad('k_McEx in component Q_Mc (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,99) = strpad('Y_Mc_AmMc in component Q_Mc (gram_per_mole)');
    LEGEND_CONSTANTS(:,100) = strpad('Y_Mc_PsMc in component Q_Mc (gram_per_mole)');
    LEGEND_ALGEBRAIC(:,75) = strpad('v_McEg in component Q_Mc (gram_per_day)');
    LEGEND_CONSTANTS(:,101) = strpad('v_McEg_star in component Q_Mc (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,102) = strpad('M_Mc_McEg in component Q_Mc (gram_per_litre)');
    LEGEND_STATES(:,14) = strpad('Q_Ac in component Q_Ac (mole)');
    LEGEND_ALGEBRAIC(:,84) = strpad('C_Ac in component Q_Ac (molar)');
    LEGEND_CONSTANTS(:,191) = strpad('P_Ac_InAc in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,160) = strpad('P_Ac_AmMa in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,120) = strpad('P_Ac_AmMc in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,168) = strpad('P_Ac_PsMa in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,132) = strpad('P_Ac_PsMc in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,176) = strpad('P_Ac_HaAs in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,152) = strpad('P_Ac_HaVa in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,88) = strpad('P_Ac_HcVa in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,92) = strpad('P_MaAc_PsAm in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,96) = strpad('P_McAc_PsAm in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,144) = strpad('P_Ac_McMa in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,100) = strpad('P_Ac_McAm in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,104) = strpad('U_Ac_AcAb in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,108) = strpad('U_Ac_AcEx in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,112) = strpad('U_Hf_AmMa in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,116) = strpad('U_Hf_AmMc in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,124) = strpad('U_Hf_PsMa in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,128) = strpad('U_Hf_PsMc in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,136) = strpad('U_Hf_HaAs in component Q_Ac (mole_per_day)');
    LEGEND_ALGEBRAIC(:,140) = strpad('U_Hf_McMa in component Q_Ac (mole_per_day)');
    LEGEND_CONSTANTS(:,103) = strpad('M_Ac_AcAb in component Q_Ac (molar)');
    LEGEND_CONSTANTS(:,200) = strpad('v_AcAb in component Q_Ac (mole_per_day)');
    LEGEND_CONSTANTS(:,104) = strpad('v_AcAb_star in component Q_Ac (mole_per_litre_day)');
    LEGEND_CONSTANTS(:,192) = strpad('Y_Ac_LaAc in component Q_Ac (dimensionless)');
    LEGEND_ALGEBRAIC(:,148) = strpad('Y_Ac_HaVa in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,194) = strpad('Y_Ac_PsVa in component Q_Ac (dimensionless)');
    LEGEND_ALGEBRAIC(:,156) = strpad('Y_Ac_AmMa in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,195) = strpad('Y_Ac_AmMc in component Q_Ac (dimensionless)');
    LEGEND_ALGEBRAIC(:,164) = strpad('Y_Ac_PsMa in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,196) = strpad('Y_Ac_PsMc in component Q_Ac (dimensionless)');
    LEGEND_ALGEBRAIC(:,172) = strpad('Y_Ac_HaAs in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,197) = strpad('Y_Ac_McMa in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,198) = strpad('Y_Ac_PsAm in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,199) = strpad('Y_Ac_McAm in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,105) = strpad('Y_Ac_WrAc in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,106) = strpad('Y_Ac_InAc in component Q_Ac (molar)');
    LEGEND_CONSTANTS(:,107) = strpad('Y_Ac_StAc in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,108) = strpad('Y_Ac_CeAc in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,109) = strpad('Y_Ac_HeAc in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,110) = strpad('Y_Ac_PsAc in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,193) = strpad('Y_Ac_HcVa in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,171) = strpad('k_AcEx in component Q_Ac (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,111) = strpad('D_Ac in component Q_Ac (litre_per_day)');
    LEGEND_CONSTANTS(:,112) = strpad('J_pH_AcAb in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,113) = strpad('phi_pH_AcAb in component Q_Ac (dimensionless)');
    LEGEND_CONSTANTS(:,114) = strpad('f_Hf_AmMa in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,115) = strpad('f_Hf_AmMc in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,116) = strpad('f_Hf_PsMa in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,117) = strpad('f_Hf_PsMc in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,118) = strpad('f_Hf_HaAs in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,119) = strpad('f_Hf_McMa in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,120) = strpad('f_Lc_Le in component model_parameters (dimensionless)');
    LEGEND_CONSTANTS(:,121) = strpad('f_Ce_Fd in component model_parameters (dimensionless)');
    LEGEND_STATES(:,15) = strpad('Q_Pr in component Q_Pr (mole)');
    LEGEND_ALGEBRAIC(:,85) = strpad('C_Pr in component Q_Pr (molar)');
    LEGEND_CONSTANTS(:,201) = strpad('P_Pr_InPr in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,161) = strpad('P_Pr_AmMa in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,121) = strpad('P_Pr_AmMc in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,169) = strpad('P_Pr_PsMa in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,133) = strpad('P_Pr_PsMc in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,177) = strpad('P_Pr_HaAs in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,153) = strpad('P_Pr_HaVa in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,89) = strpad('P_Pr_HcVa in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,93) = strpad('P_MaPr_PsAm in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,97) = strpad('P_McPr_PsAm in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,145) = strpad('P_Pr_McMa in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,101) = strpad('P_Pr_McAm in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,105) = strpad('U_Pr_PrAb in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,109) = strpad('U_Pr_PrEx in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,113) = strpad('U_Hf_AmMa in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,117) = strpad('U_Hf_AmMc in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,125) = strpad('U_Hf_PsMa in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,129) = strpad('U_Hf_PsMc in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,137) = strpad('U_Hf_HaAs in component Q_Pr (mole_per_day)');
    LEGEND_ALGEBRAIC(:,141) = strpad('U_Hf_McMa in component Q_Pr (mole_per_day)');
    LEGEND_CONSTANTS(:,122) = strpad('M_Pr_PrAb in component Q_Pr (molar)');
    LEGEND_CONSTANTS(:,210) = strpad('v_PrAb in component Q_Pr (mole_per_day)');
    LEGEND_CONSTANTS(:,123) = strpad('v_PrAb_star in component Q_Pr (mole_per_litre_day)');
    LEGEND_CONSTANTS(:,202) = strpad('Y_Pr_LaPr in component Q_Pr (dimensionless)');
    LEGEND_ALGEBRAIC(:,149) = strpad('Y_Pr_HaVa in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,204) = strpad('Y_Pr_PsVa in component Q_Pr (dimensionless)');
    LEGEND_ALGEBRAIC(:,157) = strpad('Y_Pr_AmMa in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,205) = strpad('Y_Pr_AmMc in component Q_Pr (dimensionless)');
    LEGEND_ALGEBRAIC(:,165) = strpad('Y_Pr_PsMa in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,206) = strpad('Y_Pr_PsMc in component Q_Pr (dimensionless)');
    LEGEND_ALGEBRAIC(:,173) = strpad('Y_Pr_HaAs in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,207) = strpad('Y_Pr_McMa in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,208) = strpad('Y_Pr_PsAm in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,209) = strpad('Y_Pr_McAm in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,124) = strpad('Y_Pr_WrPr in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,125) = strpad('Y_Pr_InPr in component Q_Pr (molar)');
    LEGEND_CONSTANTS(:,203) = strpad('Y_Pr_HcVa in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,126) = strpad('Y_Pr_StPr in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,127) = strpad('Y_Pr_CePr in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,128) = strpad('Y_Pr_HePr in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,129) = strpad('Y_Pr_PsPr in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,172) = strpad('k_PrEx in component Q_Pr (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,130) = strpad('D_Pr in component Q_Pr (litre_per_day)');
    LEGEND_CONSTANTS(:,131) = strpad('phi_pH_PrAb in component Q_Pr (dimensionless)');
    LEGEND_CONSTANTS(:,132) = strpad('J_pH_PrAb in component Q_Pr (dimensionless)');
    LEGEND_STATES(:,16) = strpad('Q_Bu in component Q_Bu (mole)');
    LEGEND_ALGEBRAIC(:,86) = strpad('C_Bu in component Q_Bu (molar)');
    LEGEND_CONSTANTS(:,211) = strpad('P_Bu_InBu in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,162) = strpad('P_Bu_AmMa in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,122) = strpad('P_Bu_AmMc in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,170) = strpad('P_Bu_PsMa in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,134) = strpad('P_Bu_PsMc in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,178) = strpad('P_Bu_HaAs in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,154) = strpad('P_Bu_HaVa in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,90) = strpad('P_Bu_HcVa in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,94) = strpad('P_MaBu_PsAm in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,98) = strpad('P_McBu_PsAm in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,146) = strpad('P_Bu_McMa in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,102) = strpad('P_Bu_McAm in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,106) = strpad('U_Bu_BuAb in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,110) = strpad('U_Bu_BuEx in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,114) = strpad('U_Hf_AmMa in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,118) = strpad('U_Hf_AmMc in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,126) = strpad('U_Hf_PsMa in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,130) = strpad('U_Hf_PsMc in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,138) = strpad('U_Hf_HaAs in component Q_Bu (mole_per_day)');
    LEGEND_ALGEBRAIC(:,142) = strpad('U_Hf_McMa in component Q_Bu (mole_per_day)');
    LEGEND_CONSTANTS(:,133) = strpad('M_Bu_BuAb in component Q_Bu (molar)');
    LEGEND_CONSTANTS(:,220) = strpad('v_BuAb in component Q_Bu (mole_per_day)');
    LEGEND_CONSTANTS(:,134) = strpad('v_BuAb_star in component Q_Bu (mole_per_litre_day)');
    LEGEND_CONSTANTS(:,212) = strpad('Y_Bu_LaBu in component Q_Bu (dimensionless)');
    LEGEND_ALGEBRAIC(:,150) = strpad('Y_Bu_HaVa in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,214) = strpad('Y_Bu_PsVa in component Q_Bu (dimensionless)');
    LEGEND_ALGEBRAIC(:,158) = strpad('Y_Bu_AmMa in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,215) = strpad('Y_Bu_AmMc in component Q_Bu (dimensionless)');
    LEGEND_ALGEBRAIC(:,166) = strpad('Y_Bu_PsMa in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,216) = strpad('Y_Bu_PsMc in component Q_Bu (dimensionless)');
    LEGEND_ALGEBRAIC(:,174) = strpad('Y_Bu_HaAs in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,217) = strpad('Y_Bu_McMa in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,218) = strpad('Y_Bu_PsAm in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,219) = strpad('Y_Bu_McAm in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,135) = strpad('Y_Bu_WrBu in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,136) = strpad('Y_Bu_InBu in component Q_Bu (molar)');
    LEGEND_CONSTANTS(:,213) = strpad('Y_Bu_HcVa in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,137) = strpad('Y_Bu_StBu in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,138) = strpad('Y_Bu_CeBu in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,139) = strpad('Y_Bu_HeBu in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,140) = strpad('Y_Bu_PsBu in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,173) = strpad('k_BuEx in component Q_Bu (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,141) = strpad('D_Bu in component Q_Bu (litre_per_day)');
    LEGEND_CONSTANTS(:,142) = strpad('phi_pH_BuAb in component Q_Bu (dimensionless)');
    LEGEND_CONSTANTS(:,143) = strpad('J_pH_BuAb in component Q_Bu (dimensionless)');
    LEGEND_STATES(:,17) = strpad('Q_Vl in component Q_Vl (mole)');
    LEGEND_ALGEBRAIC(:,87) = strpad('C_Vl in component Q_Vl (molar)');
    LEGEND_CONSTANTS(:,221) = strpad('P_Vl_InVl in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,163) = strpad('P_Vl_AmMa in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,123) = strpad('P_Vl_AmMc in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,171) = strpad('P_Vl_PsMa in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,135) = strpad('P_Vl_PsMc in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,179) = strpad('P_Vl_HaAs in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,155) = strpad('P_Vl_HaVa in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,91) = strpad('P_Vl_HcVa in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,95) = strpad('P_MaVl_PsAm in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,99) = strpad('P_McVl_PsAm in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,147) = strpad('P_Vl_McMa in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,103) = strpad('P_Vl_McAm in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,107) = strpad('U_Vl_VlAb in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,111) = strpad('U_Vl_VlEx in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,115) = strpad('U_Hf_AmMa in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,119) = strpad('U_Hf_AmMc in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,127) = strpad('U_Hf_PsMa in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,131) = strpad('U_Hf_PsMc in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,139) = strpad('U_Hf_HaAs in component Q_Vl (mole_per_day)');
    LEGEND_ALGEBRAIC(:,143) = strpad('U_Hf_McMa in component Q_Vl (mole_per_day)');
    LEGEND_CONSTANTS(:,144) = strpad('M_Vl_VlAb in component Q_Vl (molar)');
    LEGEND_CONSTANTS(:,230) = strpad('v_VlAb in component Q_Vl (mole_per_day)');
    LEGEND_CONSTANTS(:,145) = strpad('v_VlAb_star in component Q_Vl (mole_per_litre_day)');
    LEGEND_CONSTANTS(:,222) = strpad('Y_Vl_LaVl in component Q_Vl (dimensionless)');
    LEGEND_ALGEBRAIC(:,151) = strpad('Y_Vl_HaVa in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,224) = strpad('Y_Vl_PsVa in component Q_Vl (dimensionless)');
    LEGEND_ALGEBRAIC(:,159) = strpad('Y_Vl_AmMa in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,225) = strpad('Y_Vl_AmMc in component Q_Vl (dimensionless)');
    LEGEND_ALGEBRAIC(:,167) = strpad('Y_Vl_PsMa in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,226) = strpad('Y_Vl_PsMc in component Q_Vl (dimensionless)');
    LEGEND_ALGEBRAIC(:,175) = strpad('Y_Vl_HaAs in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,227) = strpad('Y_Vl_McMa in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,228) = strpad('Y_Vl_PsAm in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,229) = strpad('Y_Vl_McAm in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,146) = strpad('Y_Vl_WrVl in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,147) = strpad('Y_Vl_InVl in component Q_Vl (molar)');
    LEGEND_CONSTANTS(:,223) = strpad('Y_Vl_HcVa in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,148) = strpad('Y_Vl_StVl in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,149) = strpad('Y_Vl_CeVl in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,150) = strpad('Y_Vl_HeVl in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,151) = strpad('Y_Vl_PsVl in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,174) = strpad('k_VlEx in component Q_Vl (first_order_rate_constant)');
    LEGEND_CONSTANTS(:,152) = strpad('D_Vl in component Q_Vl (litre_per_day)');
    LEGEND_CONSTANTS(:,153) = strpad('phi_pH_VlAb in component Q_Vl (dimensionless)');
    LEGEND_CONSTANTS(:,154) = strpad('J_pH_VlAb in component Q_Vl (dimensionless)');
    LEGEND_RATES(:,1) = strpad('d/dt Q_Pd in component Q_Pd (mole)');
    LEGEND_RATES(:,4) = strpad('d/dt Q_Ps in component Q_Ps (mole)');
    LEGEND_RATES(:,5) = strpad('d/dt Q_Pu in component Q_Pu (mole)');
    LEGEND_RATES(:,6) = strpad('d/dt Q_Am in component Q_Am (mole)');
    LEGEND_RATES(:,7) = strpad('d/dt Q_Li in component Q_Li (mole)');
    LEGEND_RATES(:,8) = strpad('d/dt Q_Sd in component Q_Sd (mole)');
    LEGEND_RATES(:,9) = strpad('d/dt Q_Fd in component Q_Fd (mole)');
    LEGEND_RATES(:,10) = strpad('d/dt Q_Fu in component Q_Fu (mole)');
    LEGEND_RATES(:,11) = strpad('d/dt Q_Ha in component Q_Ha (mole)');
    LEGEND_RATES(:,12) = strpad('d/dt Q_Hc in component Q_Hc (mole)');
    LEGEND_RATES(:,2) = strpad('d/dt Q_Ma in component Q_Ma (gram)');
    LEGEND_RATES(:,13) = strpad('d/dt Q_As in component Q_As (gram)');
    LEGEND_RATES(:,3) = strpad('d/dt Q_Mc in component Q_Mc (gram)');
    LEGEND_RATES(:,14) = strpad('d/dt Q_Ac in component Q_Ac (mole)');
    LEGEND_RATES(:,15) = strpad('d/dt Q_Pr in component Q_Pr (mole)');
    LEGEND_RATES(:,16) = strpad('d/dt Q_Bu in component Q_Bu (mole)');
    LEGEND_RATES(:,17) = strpad('d/dt Q_Vl in component Q_Vl (mole)');
    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 = [];
    STATES(:,1) = 1.000;
    CONSTANTS(:,1) = 0.0576;
    CONSTANTS(:,2) = 0.264;
    CONSTANTS(:,3) = 0.0091;
    CONSTANTS(:,4) = 1.000;
    CONSTANTS(:,5) = 1.000;
    CONSTANTS(:,6) = 0.66;
    STATES(:,2) = 1.000;
    STATES(:,3) = 1.000;
    CONSTANTS(:,7) = 1.000;
    CONSTANTS(:,8) = 1.000;
    CONSTANTS(:,9) = 0.0248;
    STATES(:,4) = 1.000;
    CONSTANTS(:,10) = 0.0144;
    CONSTANTS(:,11) = 0.0576;
    CONSTANTS(:,12) = 0.0576;
    CONSTANTS(:,13) = 0.0091;
    CONSTANTS(:,14) = 0.0010;
    CONSTANTS(:,15) = 1.0;
    CONSTANTS(:,16) = 0.0067;
    CONSTANTS(:,17) = 0.0067;
    CONSTANTS(:,18) = 1.000;
    CONSTANTS(:,19) = 0.0289;
    CONSTANTS(:,20) = 0.0224;
    CONSTANTS(:,21) = 0.0224;
    CONSTANTS(:,22) = 0.0248;
    CONSTANTS(:,23) = 0.0248;
    CONSTANTS(:,24) = 0.0165;
    CONSTANTS(:,25) = 0.0165;
    CONSTANTS(:,26) = 1.000;
    CONSTANTS(:,27) = 0.0165;
    STATES(:,5) = 1.000;
    CONSTANTS(:,28) = 0.0091;
    CONSTANTS(:,29) = 1.000;
    STATES(:,6) = 1.000;
    CONSTANTS(:,30) = 0.00165;
    CONSTANTS(:,31) = 1.10;
    CONSTANTS(:,32) = 0.0528;
    CONSTANTS(:,33) = 0.0528;
    CONSTANTS(:,34) = 0.0588;
    CONSTANTS(:,35) = 2.0;
    CONSTANTS(:,36) = 1.257;
    CONSTANTS(:,37) = 0.0084;
    CONSTANTS(:,38) = 0.00135;
    CONSTANTS(:,39) = 0.00135;
    CONSTANTS(:,40) = 0.0132;
    CONSTANTS(:,41) = 7.5;
    CONSTANTS(:,42) = 7.85;
    CONSTANTS(:,43) = 0.0159;
    CONSTANTS(:,44) = 0.0159;
    CONSTANTS(:,45) = 0.00621;
    CONSTANTS(:,46) = 1.000;
    CONSTANTS(:,47) = 7.000;
    CONSTANTS(:,48) = 0.35;
    CONSTANTS(:,49) = 1.000;
    STATES(:,7) = 1.000;
    CONSTANTS(:,50) = 0.0015;
    CONSTANTS(:,51) = 0.00021;
    CONSTANTS(:,52) = 1.000;
    STATES(:,8) = 1.000;
    CONSTANTS(:,53) = 0.0062;
    CONSTANTS(:,54) = 1.000;
    CONSTANTS(:,55) = 0.2179;
    CONSTANTS(:,56) = 0.416;
    CONSTANTS(:,57) = 1.000;
    CONSTANTS(:,58) = 1.04;
    STATES(:,9) = 1.000;
    CONSTANTS(:,59) = 0.0062;
    CONSTANTS(:,60) = 1.000;
    CONSTANTS(:,61) = 0.1646;
    CONSTANTS(:,62) = 0.332;
    CONSTANTS(:,63) = 1.000;
    CONSTANTS(:,64) = 0.83;
    CONSTANTS(:,65) = 22.9;
    CONSTANTS(:,66) = 5.97;
    CONSTANTS(:,67) = 1.000;
    CONSTANTS(:,68) = 1.000;
    STATES(:,10) = 1.000;
    CONSTANTS(:,69) = 0.0062;
    CONSTANTS(:,70) = 1.000;
    STATES(:,11) = 1.000;
    CONSTANTS(:,71) = 1.793;
    CONSTANTS(:,72) = 1.291;
    CONSTANTS(:,73) = 0.0086;
    CONSTANTS(:,74) = 0.0268;
    CONSTANTS(:,75) = 0.055;
    CONSTANTS(:,76) = 0.053;
    CONSTANTS(:,77) = 0.1646;
    CONSTANTS(:,78) = 0.0062;
    CONSTANTS(:,79) = 0.0062;
    CONSTANTS(:,80) = 0.0025;
    CONSTANTS(:,81) = 1.0;
    CONSTANTS(:,82) = 0.0062;
    CONSTANTS(:,83) = 1.000;
    CONSTANTS(:,84) = 1.000;
    CONSTANTS(:,85) = 1.000;
    CONSTANTS(:,86) = 0.00861;
    CONSTANTS(:,87) = 0.01465;
    STATES(:,12) = 1.000;
    CONSTANTS(:,88) = 1.793;
    CONSTANTS(:,89) = 1.291;
    CONSTANTS(:,90) = 0.1646;
    CONSTANTS(:,91) = 1.0;
    CONSTANTS(:,92) = 0.055;
    CONSTANTS(:,93) = 0.00861;
    CONSTANTS(:,94) = 0.01465;
    CONSTANTS(:,95) = 118.91;
    CONSTANTS(:,96) = 149.48;
    CONSTANTS(:,97) = 149.48;
    STATES(:,13) = 1.000;
    CONSTANTS(:,98) = 112.5;
    CONSTANTS(:,99) = 118.91;
    CONSTANTS(:,100) = 149.48;
    CONSTANTS(:,101) = 15.439;
    CONSTANTS(:,102) = 34.694;
    STATES(:,14) = 1.000;
    CONSTANTS(:,103) = 0.338;
    CONSTANTS(:,104) = 7.86;
    CONSTANTS(:,105) = 1.000;
    CONSTANTS(:,106) = 0.0167;
    CONSTANTS(:,107) = 1.000;
    CONSTANTS(:,108) = 1.000;
    CONSTANTS(:,109) = 1.000;
    CONSTANTS(:,110) = 1.000;
    CONSTANTS(:,111) = 1.000;
    CONSTANTS(:,112) = 6.45;
    CONSTANTS(:,113) = 6.48;
    CONSTANTS(:,114) = 0.526;
    CONSTANTS(:,115) = 0.526;
    CONSTANTS(:,116) = 0.711;
    CONSTANTS(:,117) = 0.711;
    CONSTANTS(:,118) = 0.306;
    CONSTANTS(:,119) = 0.711;
    CONSTANTS(:,120) = 2.250;
    CONSTANTS(:,121) = 1.000;
    STATES(:,15) = 1.000;
    CONSTANTS(:,122) = 0.338;
    CONSTANTS(:,123) = 7.86;
    CONSTANTS(:,124) = 1.000;
    CONSTANTS(:,125) = 0.0135;
    CONSTANTS(:,126) = 1.000;
    CONSTANTS(:,127) = 1.000;
    CONSTANTS(:,128) = 1.000;
    CONSTANTS(:,129) = 1.000;
    CONSTANTS(:,130) = 1.000;
    CONSTANTS(:,131) = 6.48;
    CONSTANTS(:,132) = 6.45;
    STATES(:,16) = 1.000;
    CONSTANTS(:,133) = 0.338;
    CONSTANTS(:,134) = 7.86;
    CONSTANTS(:,135) = 1.000;
    CONSTANTS(:,136) = 0.0114;
    CONSTANTS(:,137) = 1.000;
    CONSTANTS(:,138) = 1.000;
    CONSTANTS(:,139) = 1.000;
    CONSTANTS(:,140) = 1.000;
    CONSTANTS(:,141) = 1.000;
    CONSTANTS(:,142) = 6.48;
    CONSTANTS(:,143) = 6.45;
    STATES(:,17) = 1.000;
    CONSTANTS(:,144) = 0.338;
    CONSTANTS(:,145) = 7.86;
    CONSTANTS(:,146) = 1.000;
    CONSTANTS(:,147) = 0.0098;
    CONSTANTS(:,148) = 1.000;
    CONSTANTS(:,149) = 1.000;
    CONSTANTS(:,150) = 1.000;
    CONSTANTS(:,151) = 1.000;
    CONSTANTS(:,152) = 1.000;
    CONSTANTS(:,153) = 6.48;
    CONSTANTS(:,154) = 6.45;
    CONSTANTS(:,155) = CONSTANTS(:,8);
    CONSTANTS(:,156) =  CONSTANTS(:,3).*CONSTANTS(:,4);
    CONSTANTS(:,157) = CONSTANTS(:,26);
    CONSTANTS(:,158) =  CONSTANTS(:,28).*CONSTANTS(:,29);
    CONSTANTS(:,159) = CONSTANTS(:,8);
    CONSTANTS(:,160) = CONSTANTS(:,26);
    CONSTANTS(:,161) = CONSTANTS(:,26);
    CONSTANTS(:,162) = CONSTANTS(:,8);
    CONSTANTS(:,163) = CONSTANTS(:,8);
    CONSTANTS(:,164) =  CONSTANTS(:,69).*CONSTANTS(:,70);
    CONSTANTS(:,165) = CONSTANTS(:,8);
    CONSTANTS(:,166) = CONSTANTS(:,26);
    CONSTANTS(:,167) = CONSTANTS(:,26);
    CONSTANTS(:,168) = CONSTANTS(:,26);
    CONSTANTS(:,169) = CONSTANTS(:,26);
    CONSTANTS(:,170) = CONSTANTS(:,8);
    CONSTANTS(:,171) = CONSTANTS(:,26);
    CONSTANTS(:,172) = CONSTANTS(:,26);
    CONSTANTS(:,173) = CONSTANTS(:,26);
    CONSTANTS(:,174) = CONSTANTS(:,26);
    CONSTANTS(:,175) =  CONSTANTS(:,2).*(CONSTANTS(:,5)./CONSTANTS(:,6));
    CONSTANTS(:,176) = CONSTANTS(:,8)./2.00000;
    CONSTANTS(:,177) =  CONSTANTS(:,13).*CONSTANTS(:,18);
    CONSTANTS(:,178) = CONSTANTS(:,168) - CONSTANTS(:,176);
    CONSTANTS(:,179) = CONSTANTS(:,168) - CONSTANTS(:,176);
    CONSTANTS(:,180) =  CONSTANTS(:,14).*CONSTANTS(:,18);
    CONSTANTS(:,181) =  CONSTANTS(:,34).*CONSTANTS(:,46);
    CONSTANTS(:,182) =  CONSTANTS(:,31).*power(1.00000, 0.250000).*power(CONSTANTS(:,7), 0.750000).*(1.00000+power(CONSTANTS(:,41)./CONSTANTS(:,47), CONSTANTS(:,42)));
    CONSTANTS(:,183) =  CONSTANTS(:,50).*CONSTANTS(:,52);
    CONSTANTS(:,184) =  CONSTANTS(:,53).*CONSTANTS(:,54);
    CONSTANTS(:,185) =  CONSTANTS(:,56).*(CONSTANTS(:,57)./CONSTANTS(:,58));
    CONSTANTS(:,186) =  CONSTANTS(:,59).*CONSTANTS(:,60);
    CONSTANTS(:,187) =  CONSTANTS(:,62).*(CONSTANTS(:,63)./CONSTANTS(:,64));
    CONSTANTS(:,188) =  CONSTANTS(:,78).*CONSTANTS(:,83);
    CONSTANTS(:,189) =  CONSTANTS(:,79).*CONSTANTS(:,84);
    CONSTANTS(:,190) =  CONSTANTS(:,80).*CONSTANTS(:,85);
    CONSTANTS(:,191) =  CONSTANTS(:,106).*CONSTANTS(:,111);
    CONSTANTS(:,192) =  CONSTANTS(:,120).*CONSTANTS(:,105);
    CONSTANTS(:,193) =  CONSTANTS(:,121).*CONSTANTS(:,108)+ (1.00000 - CONSTANTS(:,121)).*CONSTANTS(:,109);
    CONSTANTS(:,194) = CONSTANTS(:,110);
    CONSTANTS(:,195) = CONSTANTS(:,193);
    CONSTANTS(:,196) = CONSTANTS(:,193);
    CONSTANTS(:,197) = CONSTANTS(:,193);
    CONSTANTS(:,198) = CONSTANTS(:,194);
    CONSTANTS(:,199) = CONSTANTS(:,194);
    CONSTANTS(:,200) =  CONSTANTS(:,104).*power(1.00000, 0.250000).*(power(CONSTANTS(:,7), 0.750000)./(1.00000+power(CONSTANTS(:,47)./CONSTANTS(:,112), CONSTANTS(:,113))));
    CONSTANTS(:,201) =  CONSTANTS(:,125).*CONSTANTS(:,130);
    CONSTANTS(:,202) =  CONSTANTS(:,120).*CONSTANTS(:,124);
    CONSTANTS(:,203) =  CONSTANTS(:,121).*CONSTANTS(:,127)+ (1.00000 - CONSTANTS(:,121)).*CONSTANTS(:,128);
    CONSTANTS(:,204) = CONSTANTS(:,129);
    CONSTANTS(:,205) = CONSTANTS(:,203);
    CONSTANTS(:,206) = CONSTANTS(:,203);
    CONSTANTS(:,207) = CONSTANTS(:,203);
    CONSTANTS(:,208) = CONSTANTS(:,204);
    CONSTANTS(:,209) = CONSTANTS(:,204);
    CONSTANTS(:,210) =  CONSTANTS(:,123).*power(1.00000, 0.250000).*(power(CONSTANTS(:,7), 0.750000)./(1.00000+power(CONSTANTS(:,47)./CONSTANTS(:,132), CONSTANTS(:,131))));
    CONSTANTS(:,211) =  CONSTANTS(:,136).*CONSTANTS(:,141);
    CONSTANTS(:,212) =  CONSTANTS(:,120).*CONSTANTS(:,135);
    CONSTANTS(:,213) =  CONSTANTS(:,121).*CONSTANTS(:,138)+ (1.00000 - CONSTANTS(:,121)).*CONSTANTS(:,139);
    CONSTANTS(:,214) = CONSTANTS(:,140);
    CONSTANTS(:,215) = CONSTANTS(:,213);
    CONSTANTS(:,216) = CONSTANTS(:,213);
    CONSTANTS(:,217) = CONSTANTS(:,213);
    CONSTANTS(:,218) = CONSTANTS(:,214);
    CONSTANTS(:,219) = CONSTANTS(:,214);
    CONSTANTS(:,220) =  CONSTANTS(:,134).*power(1.00000, 0.250000).*(power(CONSTANTS(:,7), 0.750000)./(1.00000+power(CONSTANTS(:,47)./CONSTANTS(:,143), CONSTANTS(:,142))));
    CONSTANTS(:,221) =  CONSTANTS(:,147).*CONSTANTS(:,152);
    CONSTANTS(:,222) =  CONSTANTS(:,120).*CONSTANTS(:,146);
    CONSTANTS(:,223) =  CONSTANTS(:,121).*CONSTANTS(:,149)+ (1.00000 - CONSTANTS(:,121)).*CONSTANTS(:,150);
    CONSTANTS(:,224) = CONSTANTS(:,151);
    CONSTANTS(:,225) = CONSTANTS(:,223);
    CONSTANTS(:,226) = CONSTANTS(:,223);
    CONSTANTS(:,227) = CONSTANTS(:,223);
    CONSTANTS(:,228) = CONSTANTS(:,224);
    CONSTANTS(:,229) = CONSTANTS(:,224);
    CONSTANTS(:,230) =  CONSTANTS(:,145).*power(1.00000, 0.250000).*(power(CONSTANTS(:,7), 0.750000)./(1.00000+power(CONSTANTS(:,47)./CONSTANTS(:,154), CONSTANTS(:,153))));
    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
    ALGEBRAIC(:,3) =  CONSTANTS(:,159).*STATES(:,5);
    RATES(:,5) = CONSTANTS(:,158) - ALGEBRAIC(:,3);
    ALGEBRAIC(:,6) =  CONSTANTS(:,165).*STATES(:,10);
    RATES(:,10) = CONSTANTS(:,164) - ALGEBRAIC(:,6);
    ALGEBRAIC(:,1) = STATES(:,1)./CONSTANTS(:,7);
    ALGEBRAIC(:,10) =  CONSTANTS(:,1).*(STATES(:,2)+STATES(:,3));
    ALGEBRAIC(:,11) = ALGEBRAIC(:,10)./(1.00000+CONSTANTS(:,175)./ALGEBRAIC(:,1));
    ALGEBRAIC(:,9) =  CONSTANTS(:,155).*STATES(:,1);
    RATES(:,1) = CONSTANTS(:,156) - (ALGEBRAIC(:,11)+ALGEBRAIC(:,9));
    ALGEBRAIC(:,27) = STATES(:,8)./CONSTANTS(:,7);
    ALGEBRAIC(:,29) =  CONSTANTS(:,55).*STATES(:,2);
    ALGEBRAIC(:,30) = ALGEBRAIC(:,29)./(1.00000+CONSTANTS(:,185)./ALGEBRAIC(:,27));
    ALGEBRAIC(:,28) =  CONSTANTS(:,162).*STATES(:,8);
    RATES(:,8) = CONSTANTS(:,184) - (ALGEBRAIC(:,30)+ALGEBRAIC(:,28));
    ALGEBRAIC(:,31) = STATES(:,9)./CONSTANTS(:,7);
    ALGEBRAIC(:,33) =  CONSTANTS(:,61).*STATES(:,3).*((1.00000 - ( VOI.*CONSTANTS(:,68))./1.00000)+(( VOI.*CONSTANTS(:,68))./1.00000)./(1.00000+power(CONSTANTS(:,66)./CONSTANTS(:,67), CONSTANTS(:,65))));
    ALGEBRAIC(:,34) = ALGEBRAIC(:,33)./(1.00000+CONSTANTS(:,187)./ALGEBRAIC(:,31));
    ALGEBRAIC(:,32) =  CONSTANTS(:,163).*STATES(:,9);
    RATES(:,9) = CONSTANTS(:,186) - (ALGEBRAIC(:,34)+ALGEBRAIC(:,32));
    ALGEBRAIC(:,12) = STATES(:,4)./CONSTANTS(:,7);
    ALGEBRAIC(:,18) =  CONSTANTS(:,12).*STATES(:,3);
    ALGEBRAIC(:,48) = STATES(:,12)./CONSTANTS(:,7);
    ALGEBRAIC(:,50) = ALGEBRAIC(:,18)./(1.00000+CONSTANTS(:,21)./ALGEBRAIC(:,12)+CONSTANTS(:,23)./ALGEBRAIC(:,48));
    ALGEBRAIC(:,55) =  CONSTANTS(:,89).*ALGEBRAIC(:,50);
    ALGEBRAIC(:,19) = STATES(:,6)./CONSTANTS(:,7);
    ALGEBRAIC(:,25) =  CONSTANTS(:,33).*STATES(:,3);
    ALGEBRAIC(:,52) = ALGEBRAIC(:,25)./(1.00000+CONSTANTS(:,39)./ALGEBRAIC(:,19)+CONSTANTS(:,44)./ALGEBRAIC(:,48));
    ALGEBRAIC(:,54) =  CONSTANTS(:,88).*ALGEBRAIC(:,52);
    ALGEBRAIC(:,57) =  CONSTANTS(:,90).*STATES(:,3);
    ALGEBRAIC(:,58) = ALGEBRAIC(:,57)./( (1.00000+CONSTANTS(:,92)./ALGEBRAIC(:,48)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,93)+ALGEBRAIC(:,12)./CONSTANTS(:,94)));
    ALGEBRAIC(:,53) =  CONSTANTS(:,91).*ALGEBRAIC(:,34);
    ALGEBRAIC(:,56) =  CONSTANTS(:,167).*STATES(:,12);
    RATES(:,12) = ALGEBRAIC(:,53) - (ALGEBRAIC(:,54)+ALGEBRAIC(:,55)+ALGEBRAIC(:,58)+ALGEBRAIC(:,56));
    ALGEBRAIC(:,62) =  CONSTANTS(:,178).*STATES(:,2).*CONSTANTS(:,48);
    ALGEBRAIC(:,64) =  CONSTANTS(:,51).*ALGEBRAIC(:,62);
    ALGEBRAIC(:,26) =  CONSTANTS(:,161).*STATES(:,7);
    RATES(:,7) = (CONSTANTS(:,183)+ALGEBRAIC(:,64)) - ALGEBRAIC(:,26);
    ALGEBRAIC(:,67) =  CONSTANTS(:,179).*STATES(:,13).*CONSTANTS(:,48);
    ALGEBRAIC(:,35) = STATES(:,11)./CONSTANTS(:,7);
    ALGEBRAIC(:,44) =  CONSTANTS(:,76).*STATES(:,2);
    ALGEBRAIC(:,45) = ALGEBRAIC(:,44)./(1.00000+CONSTANTS(:,74)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,66) =  CONSTANTS(:,98).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,71) =  CONSTANTS(:,176).*STATES(:,13).*CONSTANTS(:,48);
    ALGEBRAIC(:,69) =  CONSTANTS(:,169).*STATES(:,13).*(1.00000 - CONSTANTS(:,48));
    RATES(:,13) = ALGEBRAIC(:,66) - (ALGEBRAIC(:,67)+ALGEBRAIC(:,69)+ALGEBRAIC(:,71));
    ALGEBRAIC(:,59) = STATES(:,2)./CONSTANTS(:,7);
    ALGEBRAIC(:,70) = STATES(:,3)./CONSTANTS(:,7);
    ALGEBRAIC(:,75) = ( CONSTANTS(:,101).*STATES(:,2).*CONSTANTS(:,48).*ALGEBRAIC(:,70))./(ALGEBRAIC(:,70)+ ALGEBRAIC(:,59).*(1.00000 - CONSTANTS(:,48)));
    ALGEBRAIC(:,76) = ALGEBRAIC(:,75)./(1.00000+CONSTANTS(:,102)./(ALGEBRAIC(:,70)+ ALGEBRAIC(:,59).*(1.00000 - CONSTANTS(:,48))));
    ALGEBRAIC(:,72) =  CONSTANTS(:,99).*ALGEBRAIC(:,52);
    ALGEBRAIC(:,73) =  CONSTANTS(:,100).*ALGEBRAIC(:,50);
    ALGEBRAIC(:,74) =  CONSTANTS(:,170).*STATES(:,3);
    RATES(:,3) = (ALGEBRAIC(:,72)+ALGEBRAIC(:,73)) - (ALGEBRAIC(:,76)+ALGEBRAIC(:,74));
    ALGEBRAIC(:,16) =  CONSTANTS(:,10).*STATES(:,3);
    ALGEBRAIC(:,49) = ALGEBRAIC(:,16)./(1.00000+CONSTANTS(:,19)./ALGEBRAIC(:,12)+ALGEBRAIC(:,48)./CONSTANTS(:,25));
    ALGEBRAIC(:,17) =  CONSTANTS(:,11).*STATES(:,2);
    ALGEBRAIC(:,37) = ALGEBRAIC(:,17)./(1.00000+CONSTANTS(:,20)./ALGEBRAIC(:,12)+CONSTANTS(:,22)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,15) =  CONSTANTS(:,10).*STATES(:,2);
    ALGEBRAIC(:,36) = ALGEBRAIC(:,15)./(1.00000+CONSTANTS(:,19)./ALGEBRAIC(:,12)+ALGEBRAIC(:,35)./CONSTANTS(:,24));
    ALGEBRAIC(:,13) =  CONSTANTS(:,15).*ALGEBRAIC(:,11);
    ALGEBRAIC(:,63) =  CONSTANTS(:,16).*ALGEBRAIC(:,62);
    ALGEBRAIC(:,77) =  ALGEBRAIC(:,76).*(1.00000 - (1.00000./(1.00000+ALGEBRAIC(:,35)./CONSTANTS(:,27))+1.00000./(1.00000+CONSTANTS(:,9)./ALGEBRAIC(:,35))));
    ALGEBRAIC(:,79) =  CONSTANTS(:,17).*ALGEBRAIC(:,77);
    ALGEBRAIC(:,14) =  CONSTANTS(:,157).*STATES(:,4);
    RATES(:,4) = (CONSTANTS(:,177)+ALGEBRAIC(:,13)+ALGEBRAIC(:,63)+CONSTANTS(:,180)+ALGEBRAIC(:,79)) - (ALGEBRAIC(:,36)+ALGEBRAIC(:,49)+ALGEBRAIC(:,37)+ALGEBRAIC(:,50)+ALGEBRAIC(:,14));
    ALGEBRAIC(:,24) =  CONSTANTS(:,32).*STATES(:,2).*(1.00000 - CONSTANTS(:,48));
    ALGEBRAIC(:,39) = ALGEBRAIC(:,24)./(1.00000+CONSTANTS(:,38)./ALGEBRAIC(:,19)+CONSTANTS(:,43)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,22) =  CONSTANTS(:,30).*CONSTANTS(:,7).*(CONSTANTS(:,49)./(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,45)));
    ALGEBRAIC(:,23) =  CONSTANTS(:,35).*ALGEBRAIC(:,22);
    ALGEBRAIC(:,38) =  CONSTANTS(:,36).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,51) =  CONSTANTS(:,36).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,78) = ( 1.00000.*ALGEBRAIC(:,76))./(1.00000+ALGEBRAIC(:,35)./CONSTANTS(:,27));
    ALGEBRAIC(:,80) =  CONSTANTS(:,37).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,20) = CONSTANTS(:,182)./(1.00000+CONSTANTS(:,40)./ALGEBRAIC(:,19));
    ALGEBRAIC(:,21) =  CONSTANTS(:,160).*STATES(:,6);
    RATES(:,6) = (CONSTANTS(:,181)+ALGEBRAIC(:,38)+ALGEBRAIC(:,51)+ALGEBRAIC(:,23)+ALGEBRAIC(:,80)) - (ALGEBRAIC(:,20)+ALGEBRAIC(:,39)+ALGEBRAIC(:,52)+ALGEBRAIC(:,21));
    ALGEBRAIC(:,81) = ALGEBRAIC(:,76)./(1.00000+CONSTANTS(:,9)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,83) =  CONSTANTS(:,73).*ALGEBRAIC(:,81);
    ALGEBRAIC(:,41) =  CONSTANTS(:,71).*ALGEBRAIC(:,39);
    ALGEBRAIC(:,42) =  CONSTANTS(:,72).*ALGEBRAIC(:,37);
    ALGEBRAIC(:,46) =  CONSTANTS(:,77).*STATES(:,2);
    ALGEBRAIC(:,47) =  ALGEBRAIC(:,46).*((1.00000 - CONSTANTS(:,48))./( (1.00000+CONSTANTS(:,75)./ALGEBRAIC(:,35)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,86)).*(ALGEBRAIC(:,12)./CONSTANTS(:,87))))+( ALGEBRAIC(:,46).*CONSTANTS(:,48))./(1.00000+CONSTANTS(:,75)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,40) =  CONSTANTS(:,81).*ALGEBRAIC(:,30);
    ALGEBRAIC(:,68) =  CONSTANTS(:,82).*ALGEBRAIC(:,67);
    ALGEBRAIC(:,43) =  CONSTANTS(:,166).*STATES(:,11);
    RATES(:,11) = (CONSTANTS(:,188)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)+CONSTANTS(:,190)) - (ALGEBRAIC(:,41)+ALGEBRAIC(:,42)+ALGEBRAIC(:,83)+ALGEBRAIC(:,45)+ALGEBRAIC(:,47)+ALGEBRAIC(:,43));
    ALGEBRAIC(:,60) =  CONSTANTS(:,95).*ALGEBRAIC(:,39);
    ALGEBRAIC(:,61) =  CONSTANTS(:,96).*ALGEBRAIC(:,37);
    ALGEBRAIC(:,82) =  CONSTANTS(:,97).*ALGEBRAIC(:,81);
    ALGEBRAIC(:,65) =  CONSTANTS(:,176).*STATES(:,2).*CONSTANTS(:,48);
    RATES(:,2) = (ALGEBRAIC(:,60)+ALGEBRAIC(:,61)+ALGEBRAIC(:,82)) - (ALGEBRAIC(:,62)+ALGEBRAIC(:,81)+ALGEBRAIC(:,65));
    ALGEBRAIC(:,112) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,148) = ( CONSTANTS(:,105).*CONSTANTS(:,188)+ CONSTANTS(:,192).*CONSTANTS(:,190)+ CONSTANTS(:,107).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,156) = ALGEBRAIC(:,148);
    ALGEBRAIC(:,160) =  ALGEBRAIC(:,156).*ALGEBRAIC(:,112);
    ALGEBRAIC(:,116) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,120) =  CONSTANTS(:,195).*ALGEBRAIC(:,116);
    ALGEBRAIC(:,124) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,164) = ALGEBRAIC(:,148);
    ALGEBRAIC(:,168) =  ALGEBRAIC(:,164).*ALGEBRAIC(:,124);
    ALGEBRAIC(:,128) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,132) =  CONSTANTS(:,196).*ALGEBRAIC(:,128);
    ALGEBRAIC(:,136) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,172) = ALGEBRAIC(:,148);
    ALGEBRAIC(:,176) =  ALGEBRAIC(:,172).*ALGEBRAIC(:,136);
    ALGEBRAIC(:,152) =  ALGEBRAIC(:,148).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,88) =  CONSTANTS(:,193).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,92) =  CONSTANTS(:,198).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,96) =  CONSTANTS(:,198).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,140) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,144) =  CONSTANTS(:,197).*ALGEBRAIC(:,140);
    ALGEBRAIC(:,100) =  CONSTANTS(:,199).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,84) = STATES(:,14)./CONSTANTS(:,7);
    ALGEBRAIC(:,104) = CONSTANTS(:,200)./(1.00000+CONSTANTS(:,103)./ALGEBRAIC(:,84));
    ALGEBRAIC(:,108) =  CONSTANTS(:,171).*STATES(:,14);
    RATES(:,14) = (CONSTANTS(:,191)+ALGEBRAIC(:,160)+ALGEBRAIC(:,120)+ALGEBRAIC(:,168)+ALGEBRAIC(:,132)+ALGEBRAIC(:,176)+ALGEBRAIC(:,152)+ALGEBRAIC(:,88)+ALGEBRAIC(:,92)+ALGEBRAIC(:,96)+ALGEBRAIC(:,144)+ALGEBRAIC(:,100)) - (ALGEBRAIC(:,104)+ALGEBRAIC(:,108));
    ALGEBRAIC(:,113) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,149) = ( CONSTANTS(:,124).*CONSTANTS(:,188)+ CONSTANTS(:,202).*CONSTANTS(:,190)+ CONSTANTS(:,126).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,157) = ALGEBRAIC(:,149);
    ALGEBRAIC(:,161) =  ALGEBRAIC(:,157).*ALGEBRAIC(:,113);
    ALGEBRAIC(:,117) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,121) =  CONSTANTS(:,205).*ALGEBRAIC(:,117);
    ALGEBRAIC(:,125) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,165) = ALGEBRAIC(:,149);
    ALGEBRAIC(:,169) =  ALGEBRAIC(:,165).*ALGEBRAIC(:,125);
    ALGEBRAIC(:,129) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,133) =  CONSTANTS(:,206).*ALGEBRAIC(:,129);
    ALGEBRAIC(:,137) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,173) = ALGEBRAIC(:,149);
    ALGEBRAIC(:,177) =  ALGEBRAIC(:,173).*ALGEBRAIC(:,137);
    ALGEBRAIC(:,153) =  ALGEBRAIC(:,149).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,89) =  CONSTANTS(:,203).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,93) =  CONSTANTS(:,208).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,97) =  CONSTANTS(:,208).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,141) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,145) =  CONSTANTS(:,207).*ALGEBRAIC(:,141);
    ALGEBRAIC(:,101) =  CONSTANTS(:,209).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,85) = STATES(:,15)./CONSTANTS(:,7);
    ALGEBRAIC(:,105) = CONSTANTS(:,210)./(1.00000+CONSTANTS(:,122)./ALGEBRAIC(:,85));
    ALGEBRAIC(:,109) =  CONSTANTS(:,172).*STATES(:,15);
    RATES(:,15) = (CONSTANTS(:,201)+ALGEBRAIC(:,161)+ALGEBRAIC(:,121)+ALGEBRAIC(:,169)+ALGEBRAIC(:,133)+ALGEBRAIC(:,177)+ALGEBRAIC(:,153)+ALGEBRAIC(:,89)+ALGEBRAIC(:,93)+ALGEBRAIC(:,97)+ALGEBRAIC(:,145)+ALGEBRAIC(:,101)) - (ALGEBRAIC(:,105)+ALGEBRAIC(:,109));
    ALGEBRAIC(:,114) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,150) = ( CONSTANTS(:,135).*CONSTANTS(:,188)+ CONSTANTS(:,212).*CONSTANTS(:,190)+ CONSTANTS(:,137).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,158) = ALGEBRAIC(:,150);
    ALGEBRAIC(:,162) =  ALGEBRAIC(:,158).*ALGEBRAIC(:,114);
    ALGEBRAIC(:,118) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,122) =  CONSTANTS(:,215).*ALGEBRAIC(:,118);
    ALGEBRAIC(:,126) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,166) = ALGEBRAIC(:,150);
    ALGEBRAIC(:,170) =  ALGEBRAIC(:,166).*ALGEBRAIC(:,126);
    ALGEBRAIC(:,130) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,134) =  CONSTANTS(:,216).*ALGEBRAIC(:,130);
    ALGEBRAIC(:,138) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,174) = ALGEBRAIC(:,150);
    ALGEBRAIC(:,178) =  ALGEBRAIC(:,174).*ALGEBRAIC(:,138);
    ALGEBRAIC(:,154) =  ALGEBRAIC(:,150).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,90) =  CONSTANTS(:,213).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,94) =  CONSTANTS(:,218).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,98) =  CONSTANTS(:,218).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,142) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,146) =  CONSTANTS(:,217).*ALGEBRAIC(:,142);
    ALGEBRAIC(:,102) =  CONSTANTS(:,219).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,86) = STATES(:,16)./CONSTANTS(:,7);
    ALGEBRAIC(:,106) = CONSTANTS(:,220)./(1.00000+CONSTANTS(:,133)./ALGEBRAIC(:,86));
    ALGEBRAIC(:,110) =  CONSTANTS(:,173).*STATES(:,16);
    RATES(:,16) = (CONSTANTS(:,211)+ALGEBRAIC(:,162)+ALGEBRAIC(:,122)+ALGEBRAIC(:,170)+ALGEBRAIC(:,134)+ALGEBRAIC(:,178)+ALGEBRAIC(:,154)+ALGEBRAIC(:,90)+ALGEBRAIC(:,94)+ALGEBRAIC(:,98)+ALGEBRAIC(:,146)+ALGEBRAIC(:,102)) - (ALGEBRAIC(:,106)+ALGEBRAIC(:,110));
    ALGEBRAIC(:,115) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,151) = ( CONSTANTS(:,146).*CONSTANTS(:,188)+ CONSTANTS(:,222).*CONSTANTS(:,190)+ CONSTANTS(:,148).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,159) = ALGEBRAIC(:,151);
    ALGEBRAIC(:,163) =  ALGEBRAIC(:,159).*ALGEBRAIC(:,115);
    ALGEBRAIC(:,119) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,123) =  CONSTANTS(:,225).*ALGEBRAIC(:,119);
    ALGEBRAIC(:,127) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,167) = ALGEBRAIC(:,151);
    ALGEBRAIC(:,171) =  ALGEBRAIC(:,167).*ALGEBRAIC(:,127);
    ALGEBRAIC(:,131) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,135) =  CONSTANTS(:,226).*ALGEBRAIC(:,131);
    ALGEBRAIC(:,139) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,175) = ALGEBRAIC(:,151);
    ALGEBRAIC(:,179) =  ALGEBRAIC(:,175).*ALGEBRAIC(:,139);
    ALGEBRAIC(:,155) =  ALGEBRAIC(:,151).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,91) =  CONSTANTS(:,223).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,95) =  CONSTANTS(:,228).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,99) =  CONSTANTS(:,228).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,143) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,147) =  CONSTANTS(:,227).*ALGEBRAIC(:,143);
    ALGEBRAIC(:,103) =  CONSTANTS(:,229).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,87) = STATES(:,17)./CONSTANTS(:,7);
    ALGEBRAIC(:,107) = CONSTANTS(:,230)./(1.00000+CONSTANTS(:,144)./ALGEBRAIC(:,87));
    ALGEBRAIC(:,111) =  CONSTANTS(:,174).*STATES(:,17);
    RATES(:,17) = (CONSTANTS(:,221)+ALGEBRAIC(:,163)+ALGEBRAIC(:,123)+ALGEBRAIC(:,171)+ALGEBRAIC(:,135)+ALGEBRAIC(:,179)+ALGEBRAIC(:,155)+ALGEBRAIC(:,91)+ALGEBRAIC(:,95)+ALGEBRAIC(:,99)+ALGEBRAIC(:,147)+ALGEBRAIC(:,103)) - (ALGEBRAIC(:,107)+ALGEBRAIC(:,111));
   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(:,3) =  CONSTANTS(:,159).*STATES(:,5);
    ALGEBRAIC(:,6) =  CONSTANTS(:,165).*STATES(:,10);
    ALGEBRAIC(:,1) = STATES(:,1)./CONSTANTS(:,7);
    ALGEBRAIC(:,10) =  CONSTANTS(:,1).*(STATES(:,2)+STATES(:,3));
    ALGEBRAIC(:,11) = ALGEBRAIC(:,10)./(1.00000+CONSTANTS(:,175)./ALGEBRAIC(:,1));
    ALGEBRAIC(:,9) =  CONSTANTS(:,155).*STATES(:,1);
    ALGEBRAIC(:,27) = STATES(:,8)./CONSTANTS(:,7);
    ALGEBRAIC(:,29) =  CONSTANTS(:,55).*STATES(:,2);
    ALGEBRAIC(:,30) = ALGEBRAIC(:,29)./(1.00000+CONSTANTS(:,185)./ALGEBRAIC(:,27));
    ALGEBRAIC(:,28) =  CONSTANTS(:,162).*STATES(:,8);
    ALGEBRAIC(:,31) = STATES(:,9)./CONSTANTS(:,7);
    ALGEBRAIC(:,33) =  CONSTANTS(:,61).*STATES(:,3).*((1.00000 - ( VOI.*CONSTANTS(:,68))./1.00000)+(( VOI.*CONSTANTS(:,68))./1.00000)./(1.00000+power(CONSTANTS(:,66)./CONSTANTS(:,67), CONSTANTS(:,65))));
    ALGEBRAIC(:,34) = ALGEBRAIC(:,33)./(1.00000+CONSTANTS(:,187)./ALGEBRAIC(:,31));
    ALGEBRAIC(:,32) =  CONSTANTS(:,163).*STATES(:,9);
    ALGEBRAIC(:,12) = STATES(:,4)./CONSTANTS(:,7);
    ALGEBRAIC(:,18) =  CONSTANTS(:,12).*STATES(:,3);
    ALGEBRAIC(:,48) = STATES(:,12)./CONSTANTS(:,7);
    ALGEBRAIC(:,50) = ALGEBRAIC(:,18)./(1.00000+CONSTANTS(:,21)./ALGEBRAIC(:,12)+CONSTANTS(:,23)./ALGEBRAIC(:,48));
    ALGEBRAIC(:,55) =  CONSTANTS(:,89).*ALGEBRAIC(:,50);
    ALGEBRAIC(:,19) = STATES(:,6)./CONSTANTS(:,7);
    ALGEBRAIC(:,25) =  CONSTANTS(:,33).*STATES(:,3);
    ALGEBRAIC(:,52) = ALGEBRAIC(:,25)./(1.00000+CONSTANTS(:,39)./ALGEBRAIC(:,19)+CONSTANTS(:,44)./ALGEBRAIC(:,48));
    ALGEBRAIC(:,54) =  CONSTANTS(:,88).*ALGEBRAIC(:,52);
    ALGEBRAIC(:,57) =  CONSTANTS(:,90).*STATES(:,3);
    ALGEBRAIC(:,58) = ALGEBRAIC(:,57)./( (1.00000+CONSTANTS(:,92)./ALGEBRAIC(:,48)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,93)+ALGEBRAIC(:,12)./CONSTANTS(:,94)));
    ALGEBRAIC(:,53) =  CONSTANTS(:,91).*ALGEBRAIC(:,34);
    ALGEBRAIC(:,56) =  CONSTANTS(:,167).*STATES(:,12);
    ALGEBRAIC(:,62) =  CONSTANTS(:,178).*STATES(:,2).*CONSTANTS(:,48);
    ALGEBRAIC(:,64) =  CONSTANTS(:,51).*ALGEBRAIC(:,62);
    ALGEBRAIC(:,26) =  CONSTANTS(:,161).*STATES(:,7);
    ALGEBRAIC(:,67) =  CONSTANTS(:,179).*STATES(:,13).*CONSTANTS(:,48);
    ALGEBRAIC(:,35) = STATES(:,11)./CONSTANTS(:,7);
    ALGEBRAIC(:,44) =  CONSTANTS(:,76).*STATES(:,2);
    ALGEBRAIC(:,45) = ALGEBRAIC(:,44)./(1.00000+CONSTANTS(:,74)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,66) =  CONSTANTS(:,98).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,71) =  CONSTANTS(:,176).*STATES(:,13).*CONSTANTS(:,48);
    ALGEBRAIC(:,69) =  CONSTANTS(:,169).*STATES(:,13).*(1.00000 - CONSTANTS(:,48));
    ALGEBRAIC(:,59) = STATES(:,2)./CONSTANTS(:,7);
    ALGEBRAIC(:,70) = STATES(:,3)./CONSTANTS(:,7);
    ALGEBRAIC(:,75) = ( CONSTANTS(:,101).*STATES(:,2).*CONSTANTS(:,48).*ALGEBRAIC(:,70))./(ALGEBRAIC(:,70)+ ALGEBRAIC(:,59).*(1.00000 - CONSTANTS(:,48)));
    ALGEBRAIC(:,76) = ALGEBRAIC(:,75)./(1.00000+CONSTANTS(:,102)./(ALGEBRAIC(:,70)+ ALGEBRAIC(:,59).*(1.00000 - CONSTANTS(:,48))));
    ALGEBRAIC(:,72) =  CONSTANTS(:,99).*ALGEBRAIC(:,52);
    ALGEBRAIC(:,73) =  CONSTANTS(:,100).*ALGEBRAIC(:,50);
    ALGEBRAIC(:,74) =  CONSTANTS(:,170).*STATES(:,3);
    ALGEBRAIC(:,16) =  CONSTANTS(:,10).*STATES(:,3);
    ALGEBRAIC(:,49) = ALGEBRAIC(:,16)./(1.00000+CONSTANTS(:,19)./ALGEBRAIC(:,12)+ALGEBRAIC(:,48)./CONSTANTS(:,25));
    ALGEBRAIC(:,17) =  CONSTANTS(:,11).*STATES(:,2);
    ALGEBRAIC(:,37) = ALGEBRAIC(:,17)./(1.00000+CONSTANTS(:,20)./ALGEBRAIC(:,12)+CONSTANTS(:,22)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,15) =  CONSTANTS(:,10).*STATES(:,2);
    ALGEBRAIC(:,36) = ALGEBRAIC(:,15)./(1.00000+CONSTANTS(:,19)./ALGEBRAIC(:,12)+ALGEBRAIC(:,35)./CONSTANTS(:,24));
    ALGEBRAIC(:,13) =  CONSTANTS(:,15).*ALGEBRAIC(:,11);
    ALGEBRAIC(:,63) =  CONSTANTS(:,16).*ALGEBRAIC(:,62);
    ALGEBRAIC(:,77) =  ALGEBRAIC(:,76).*(1.00000 - (1.00000./(1.00000+ALGEBRAIC(:,35)./CONSTANTS(:,27))+1.00000./(1.00000+CONSTANTS(:,9)./ALGEBRAIC(:,35))));
    ALGEBRAIC(:,79) =  CONSTANTS(:,17).*ALGEBRAIC(:,77);
    ALGEBRAIC(:,14) =  CONSTANTS(:,157).*STATES(:,4);
    ALGEBRAIC(:,24) =  CONSTANTS(:,32).*STATES(:,2).*(1.00000 - CONSTANTS(:,48));
    ALGEBRAIC(:,39) = ALGEBRAIC(:,24)./(1.00000+CONSTANTS(:,38)./ALGEBRAIC(:,19)+CONSTANTS(:,43)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,22) =  CONSTANTS(:,30).*CONSTANTS(:,7).*(CONSTANTS(:,49)./(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,45)));
    ALGEBRAIC(:,23) =  CONSTANTS(:,35).*ALGEBRAIC(:,22);
    ALGEBRAIC(:,38) =  CONSTANTS(:,36).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,51) =  CONSTANTS(:,36).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,78) = ( 1.00000.*ALGEBRAIC(:,76))./(1.00000+ALGEBRAIC(:,35)./CONSTANTS(:,27));
    ALGEBRAIC(:,80) =  CONSTANTS(:,37).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,20) = CONSTANTS(:,182)./(1.00000+CONSTANTS(:,40)./ALGEBRAIC(:,19));
    ALGEBRAIC(:,21) =  CONSTANTS(:,160).*STATES(:,6);
    ALGEBRAIC(:,81) = ALGEBRAIC(:,76)./(1.00000+CONSTANTS(:,9)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,83) =  CONSTANTS(:,73).*ALGEBRAIC(:,81);
    ALGEBRAIC(:,41) =  CONSTANTS(:,71).*ALGEBRAIC(:,39);
    ALGEBRAIC(:,42) =  CONSTANTS(:,72).*ALGEBRAIC(:,37);
    ALGEBRAIC(:,46) =  CONSTANTS(:,77).*STATES(:,2);
    ALGEBRAIC(:,47) =  ALGEBRAIC(:,46).*((1.00000 - CONSTANTS(:,48))./( (1.00000+CONSTANTS(:,75)./ALGEBRAIC(:,35)).*(1.00000+ALGEBRAIC(:,19)./CONSTANTS(:,86)).*(ALGEBRAIC(:,12)./CONSTANTS(:,87))))+( ALGEBRAIC(:,46).*CONSTANTS(:,48))./(1.00000+CONSTANTS(:,75)./ALGEBRAIC(:,35));
    ALGEBRAIC(:,40) =  CONSTANTS(:,81).*ALGEBRAIC(:,30);
    ALGEBRAIC(:,68) =  CONSTANTS(:,82).*ALGEBRAIC(:,67);
    ALGEBRAIC(:,43) =  CONSTANTS(:,166).*STATES(:,11);
    ALGEBRAIC(:,60) =  CONSTANTS(:,95).*ALGEBRAIC(:,39);
    ALGEBRAIC(:,61) =  CONSTANTS(:,96).*ALGEBRAIC(:,37);
    ALGEBRAIC(:,82) =  CONSTANTS(:,97).*ALGEBRAIC(:,81);
    ALGEBRAIC(:,65) =  CONSTANTS(:,176).*STATES(:,2).*CONSTANTS(:,48);
    ALGEBRAIC(:,112) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,148) = ( CONSTANTS(:,105).*CONSTANTS(:,188)+ CONSTANTS(:,192).*CONSTANTS(:,190)+ CONSTANTS(:,107).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,156) = ALGEBRAIC(:,148);
    ALGEBRAIC(:,160) =  ALGEBRAIC(:,156).*ALGEBRAIC(:,112);
    ALGEBRAIC(:,116) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,120) =  CONSTANTS(:,195).*ALGEBRAIC(:,116);
    ALGEBRAIC(:,124) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,164) = ALGEBRAIC(:,148);
    ALGEBRAIC(:,168) =  ALGEBRAIC(:,164).*ALGEBRAIC(:,124);
    ALGEBRAIC(:,128) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,132) =  CONSTANTS(:,196).*ALGEBRAIC(:,128);
    ALGEBRAIC(:,136) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,172) = ALGEBRAIC(:,148);
    ALGEBRAIC(:,176) =  ALGEBRAIC(:,172).*ALGEBRAIC(:,136);
    ALGEBRAIC(:,152) =  ALGEBRAIC(:,148).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,88) =  CONSTANTS(:,193).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,92) =  CONSTANTS(:,198).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,96) =  CONSTANTS(:,198).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,140) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,144) =  CONSTANTS(:,197).*ALGEBRAIC(:,140);
    ALGEBRAIC(:,100) =  CONSTANTS(:,199).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,84) = STATES(:,14)./CONSTANTS(:,7);
    ALGEBRAIC(:,104) = CONSTANTS(:,200)./(1.00000+CONSTANTS(:,103)./ALGEBRAIC(:,84));
    ALGEBRAIC(:,108) =  CONSTANTS(:,171).*STATES(:,14);
    ALGEBRAIC(:,113) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,149) = ( CONSTANTS(:,124).*CONSTANTS(:,188)+ CONSTANTS(:,202).*CONSTANTS(:,190)+ CONSTANTS(:,126).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,157) = ALGEBRAIC(:,149);
    ALGEBRAIC(:,161) =  ALGEBRAIC(:,157).*ALGEBRAIC(:,113);
    ALGEBRAIC(:,117) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,121) =  CONSTANTS(:,205).*ALGEBRAIC(:,117);
    ALGEBRAIC(:,125) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,165) = ALGEBRAIC(:,149);
    ALGEBRAIC(:,169) =  ALGEBRAIC(:,165).*ALGEBRAIC(:,125);
    ALGEBRAIC(:,129) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,133) =  CONSTANTS(:,206).*ALGEBRAIC(:,129);
    ALGEBRAIC(:,137) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,173) = ALGEBRAIC(:,149);
    ALGEBRAIC(:,177) =  ALGEBRAIC(:,173).*ALGEBRAIC(:,137);
    ALGEBRAIC(:,153) =  ALGEBRAIC(:,149).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,89) =  CONSTANTS(:,203).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,93) =  CONSTANTS(:,208).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,97) =  CONSTANTS(:,208).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,141) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,145) =  CONSTANTS(:,207).*ALGEBRAIC(:,141);
    ALGEBRAIC(:,101) =  CONSTANTS(:,209).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,85) = STATES(:,15)./CONSTANTS(:,7);
    ALGEBRAIC(:,105) = CONSTANTS(:,210)./(1.00000+CONSTANTS(:,122)./ALGEBRAIC(:,85));
    ALGEBRAIC(:,109) =  CONSTANTS(:,172).*STATES(:,15);
    ALGEBRAIC(:,114) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,150) = ( CONSTANTS(:,135).*CONSTANTS(:,188)+ CONSTANTS(:,212).*CONSTANTS(:,190)+ CONSTANTS(:,137).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,158) = ALGEBRAIC(:,150);
    ALGEBRAIC(:,162) =  ALGEBRAIC(:,158).*ALGEBRAIC(:,114);
    ALGEBRAIC(:,118) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,122) =  CONSTANTS(:,215).*ALGEBRAIC(:,118);
    ALGEBRAIC(:,126) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,166) = ALGEBRAIC(:,150);
    ALGEBRAIC(:,170) =  ALGEBRAIC(:,166).*ALGEBRAIC(:,126);
    ALGEBRAIC(:,130) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,134) =  CONSTANTS(:,216).*ALGEBRAIC(:,130);
    ALGEBRAIC(:,138) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,174) = ALGEBRAIC(:,150);
    ALGEBRAIC(:,178) =  ALGEBRAIC(:,174).*ALGEBRAIC(:,138);
    ALGEBRAIC(:,154) =  ALGEBRAIC(:,150).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,90) =  CONSTANTS(:,213).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,94) =  CONSTANTS(:,218).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,98) =  CONSTANTS(:,218).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,142) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,146) =  CONSTANTS(:,217).*ALGEBRAIC(:,142);
    ALGEBRAIC(:,102) =  CONSTANTS(:,219).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,86) = STATES(:,16)./CONSTANTS(:,7);
    ALGEBRAIC(:,106) = CONSTANTS(:,220)./(1.00000+CONSTANTS(:,133)./ALGEBRAIC(:,86));
    ALGEBRAIC(:,110) =  CONSTANTS(:,173).*STATES(:,16);
    ALGEBRAIC(:,115) =  CONSTANTS(:,114).*ALGEBRAIC(:,41);
    ALGEBRAIC(:,151) = ( CONSTANTS(:,146).*CONSTANTS(:,188)+ CONSTANTS(:,222).*CONSTANTS(:,190)+ CONSTANTS(:,148).*(CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68)))./(CONSTANTS(:,188)+CONSTANTS(:,190)+CONSTANTS(:,189)+ALGEBRAIC(:,40)+ALGEBRAIC(:,68));
    ALGEBRAIC(:,159) = ALGEBRAIC(:,151);
    ALGEBRAIC(:,163) =  ALGEBRAIC(:,159).*ALGEBRAIC(:,115);
    ALGEBRAIC(:,119) =  CONSTANTS(:,115).*ALGEBRAIC(:,54);
    ALGEBRAIC(:,123) =  CONSTANTS(:,225).*ALGEBRAIC(:,119);
    ALGEBRAIC(:,127) =  CONSTANTS(:,116).*ALGEBRAIC(:,42);
    ALGEBRAIC(:,167) = ALGEBRAIC(:,151);
    ALGEBRAIC(:,171) =  ALGEBRAIC(:,167).*ALGEBRAIC(:,127);
    ALGEBRAIC(:,131) =  CONSTANTS(:,117).*ALGEBRAIC(:,55);
    ALGEBRAIC(:,135) =  CONSTANTS(:,226).*ALGEBRAIC(:,131);
    ALGEBRAIC(:,139) =  CONSTANTS(:,118).*ALGEBRAIC(:,45);
    ALGEBRAIC(:,175) = ALGEBRAIC(:,151);
    ALGEBRAIC(:,179) =  ALGEBRAIC(:,175).*ALGEBRAIC(:,139);
    ALGEBRAIC(:,155) =  ALGEBRAIC(:,151).*ALGEBRAIC(:,47);
    ALGEBRAIC(:,91) =  CONSTANTS(:,223).*ALGEBRAIC(:,58);
    ALGEBRAIC(:,95) =  CONSTANTS(:,228).*ALGEBRAIC(:,36);
    ALGEBRAIC(:,99) =  CONSTANTS(:,228).*ALGEBRAIC(:,49);
    ALGEBRAIC(:,143) =  CONSTANTS(:,119).*ALGEBRAIC(:,83);
    ALGEBRAIC(:,147) =  CONSTANTS(:,227).*ALGEBRAIC(:,143);
    ALGEBRAIC(:,103) =  CONSTANTS(:,229).*ALGEBRAIC(:,78);
    ALGEBRAIC(:,87) = STATES(:,17)./CONSTANTS(:,7);
    ALGEBRAIC(:,107) = CONSTANTS(:,230)./(1.00000+CONSTANTS(:,144)./ALGEBRAIC(:,87));
    ALGEBRAIC(:,111) =  CONSTANTS(:,174).*STATES(:,17);
    ALGEBRAIC(:,2) = STATES(:,5)./CONSTANTS(:,7);
    ALGEBRAIC(:,4) = STATES(:,7)./CONSTANTS(:,7);
    ALGEBRAIC(:,5) = STATES(:,10)./CONSTANTS(:,7);
    ALGEBRAIC(:,7) =  CONSTANTS(:,168).*STATES(:,2).*(1.00000 - CONSTANTS(:,48));
    ALGEBRAIC(:,8) = STATES(:,13)./CONSTANTS(:,7);
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