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 =312;
end
% There are a total of 25 entries in each of the rate and state variable arrays.
% There are a total of 342 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 (minute)');
LEGEND_ALGEBRAIC(:,1) = strpad('pH_cy in component environment (dimensionless)');
LEGEND_CONSTANTS(:,35) = strpad('addbuffer in component environment (dimensionless)');
LEGEND_CONSTANTS(:,36) = strpad('pHstat in component environment (dimensionless)');
LEGEND_CONSTANTS(:,3) = strpad('Par_97 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,4) = strpad('Par_98 in component parameters (dimensionless)');
LEGEND_STATES(:,1) = strpad('pH_calc in component differential_equations (dimensionless)');
LEGEND_CONSTANTS(:,1) = strpad('R in component global_parameters (kilojoule_per_kelvin_per_mole)');
LEGEND_CONSTANTS(:,2) = strpad('T1 in component global_parameters (kelvin)');
LEGEND_CONSTANTS(:,131) = strpad('T in component global_parameters (kelvin)');
LEGEND_CONSTANTS(:,126) = strpad('I in component global_parameters (molar)');
LEGEND_CONSTANTS(:,125) = strpad('Par_90 in component parameters (molar)');
LEGEND_CONSTANTS(:,130) = strpad('Par_94 in component parameters (kelvin)');
LEGEND_CONSTANTS(:,5) = strpad('Par_1 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,6) = strpad('Par_2 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,7) = strpad('Par_3 in component parameters (molar)');
LEGEND_CONSTANTS(:,37) = strpad('Par_4 in component parameters (molar)');
LEGEND_CONSTANTS(:,57) = strpad('Par_5 in component parameters (molar)');
LEGEND_CONSTANTS(:,8) = strpad('Par_6 in component parameters (molar)');
LEGEND_CONSTANTS(:,58) = strpad('Par_7 in component parameters (molar)');
LEGEND_CONSTANTS(:,9) = strpad('Par_8 in component parameters (molar)');
LEGEND_CONSTANTS(:,59) = strpad('Par_9 in component parameters (molar)');
LEGEND_CONSTANTS(:,60) = strpad('Par_10 in component parameters (molar)');
LEGEND_CONSTANTS(:,10) = strpad('Par_11 in component parameters (molar)');
LEGEND_CONSTANTS(:,61) = strpad('Par_12 in component parameters (molar)');
LEGEND_CONSTANTS(:,62) = strpad('Par_13 in component parameters (molar)');
LEGEND_CONSTANTS(:,11) = strpad('Par_14 in component parameters (molar)');
LEGEND_CONSTANTS(:,63) = strpad('Par_15 in component parameters (molar)');
LEGEND_CONSTANTS(:,12) = strpad('Par_16 in component parameters (molar)');
LEGEND_CONSTANTS(:,13) = strpad('Par_17 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,14) = strpad('Par_18 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,64) = strpad('Par_19 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,65) = strpad('Par_20 in component parameters (molar)');
LEGEND_CONSTANTS(:,66) = strpad('Par_21 in component parameters (molar)');
LEGEND_CONSTANTS(:,67) = strpad('Par_22 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,68) = strpad('Par_23 in component parameters (molar)');
LEGEND_CONSTANTS(:,69) = strpad('Par_24 in component parameters (molar)');
LEGEND_CONSTANTS(:,70) = strpad('Par_25 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,71) = strpad('Par_26 in component parameters (molar)');
LEGEND_CONSTANTS(:,72) = strpad('Par_27 in component parameters (molar)');
LEGEND_CONSTANTS(:,73) = strpad('Par_28 in component parameters (molar)');
LEGEND_CONSTANTS(:,74) = strpad('Par_29 in component parameters (molar)');
LEGEND_CONSTANTS(:,75) = strpad('Par_30 in component parameters (molar)');
LEGEND_CONSTANTS(:,15) = strpad('Par_31 in component parameters (molar)');
LEGEND_CONSTANTS(:,76) = strpad('Par_32 in component parameters (molar)');
LEGEND_CONSTANTS(:,16) = strpad('Par_33 in component parameters (molar)');
LEGEND_CONSTANTS(:,17) = strpad('Par_34 in component parameters (molar)');
LEGEND_CONSTANTS(:,18) = strpad('Par_35 in component parameters (molar)');
LEGEND_CONSTANTS(:,19) = strpad('Par_36 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,20) = strpad('Par_37 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,21) = strpad('Par_38 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,77) = strpad('Par_39 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,78) = strpad('Par_40 in component parameters (molar)');
LEGEND_CONSTANTS(:,79) = strpad('Par_41 in component parameters (molar)');
LEGEND_CONSTANTS(:,80) = strpad('Par_42 in component parameters (molar)');
LEGEND_CONSTANTS(:,81) = strpad('Par_43 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,82) = strpad('Par_44 in component parameters (molar)');
LEGEND_CONSTANTS(:,83) = strpad('Par_45 in component parameters (molar)');
LEGEND_CONSTANTS(:,84) = strpad('Par_46 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,85) = strpad('Par_47 in component parameters (molar)');
LEGEND_CONSTANTS(:,86) = strpad('Par_48 in component parameters (molar)');
LEGEND_CONSTANTS(:,87) = strpad('Par_49 in component parameters (molar)');
LEGEND_CONSTANTS(:,88) = strpad('Par_50 in component parameters (molar)');
LEGEND_CONSTANTS(:,89) = strpad('Par_51 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,90) = strpad('Par_52 in component parameters (molar)');
LEGEND_CONSTANTS(:,91) = strpad('Par_53 in component parameters (molar)');
LEGEND_CONSTANTS(:,92) = strpad('Par_54 in component parameters (molar)');
LEGEND_CONSTANTS(:,93) = strpad('Par_55 in component parameters (molar)');
LEGEND_CONSTANTS(:,94) = strpad('Par_56 in component parameters (molar)');
LEGEND_CONSTANTS(:,95) = strpad('Par_57 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,96) = strpad('Par_58 in component parameters (molar)');
LEGEND_CONSTANTS(:,97) = strpad('Par_59 in component parameters (molar)');
LEGEND_CONSTANTS(:,98) = strpad('Par_60 in component parameters (molar)');
LEGEND_CONSTANTS(:,99) = strpad('Par_61 in component parameters (molar)');
LEGEND_CONSTANTS(:,100) = strpad('Par_62 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,101) = strpad('Par_63 in component parameters (molar)');
LEGEND_CONSTANTS(:,102) = strpad('Par_64 in component parameters (molar)');
LEGEND_CONSTANTS(:,103) = strpad('Par_65 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,104) = strpad('Par_66 in component parameters (molar)');
LEGEND_CONSTANTS(:,105) = strpad('Par_67 in component parameters (molar)');
LEGEND_CONSTANTS(:,106) = strpad('Par_68 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,107) = strpad('Par_69 in component parameters (molar)');
LEGEND_CONSTANTS(:,108) = strpad('Par_70 in component parameters (molar)');
LEGEND_CONSTANTS(:,109) = strpad('Par_71 in component parameters (molar)');
LEGEND_CONSTANTS(:,110) = strpad('Par_72 in component parameters (molar)');
LEGEND_CONSTANTS(:,111) = strpad('Par_73 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,112) = strpad('Par_74 in component parameters (molar)');
LEGEND_CONSTANTS(:,113) = strpad('Par_75 in component parameters (molar)');
LEGEND_CONSTANTS(:,114) = strpad('Par_76 in component parameters (molar)');
LEGEND_CONSTANTS(:,115) = strpad('Par_77 in component parameters (molar)');
LEGEND_CONSTANTS(:,116) = strpad('Par_78 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,117) = strpad('Par_79 in component parameters (molar)');
LEGEND_CONSTANTS(:,118) = strpad('Par_80 in component parameters (molar)');
LEGEND_CONSTANTS(:,119) = strpad('Par_81 in component parameters (molar)');
LEGEND_CONSTANTS(:,22) = strpad('Par_82 in component parameters (molar)');
LEGEND_CONSTANTS(:,120) = strpad('Par_83 in component parameters (molar)');
LEGEND_CONSTANTS(:,121) = strpad('Par_84 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,122) = strpad('Par_85 in component parameters (molar)');
LEGEND_CONSTANTS(:,123) = strpad('Par_86 in component parameters (molar)');
LEGEND_CONSTANTS(:,124) = strpad('Par_87 in component parameters (molar)');
LEGEND_CONSTANTS(:,23) = strpad('Par_88 in component parameters (molar_per_minute)');
LEGEND_CONSTANTS(:,24) = strpad('Par_89 in component parameters (molar)');
LEGEND_CONSTANTS(:,127) = strpad('Par_91 in component parameters (molar)');
LEGEND_CONSTANTS(:,128) = strpad('Par_92 in component parameters (molar)');
LEGEND_CONSTANTS(:,129) = strpad('Par_93 in component parameters (molar)');
LEGEND_CONSTANTS(:,25) = strpad('Par_95 in component parameters (molar)');
LEGEND_CONSTANTS(:,132) = strpad('Par_96 in component parameters (molar)');
LEGEND_CONSTANTS(:,133) = strpad('Par_99 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,26) = strpad('Par_100 in component parameters (dimensionless)');
LEGEND_CONSTANTS(:,38) = strpad('mgT in component equilibrium_constants (molar)');
LEGEND_CONSTANTS(:,134) = strpad('k in component equilibrium_constants (molar)');
LEGEND_CONSTANTS(:,135) = strpad('c0 in component equilibrium_constants (molar)');
LEGEND_CONSTANTS(:,136) = strpad('RT2dadT in component correction_factors (kilojoule_half_liter_per_3_half_mole)');
LEGEND_CONSTANTS(:,137) = strpad('B in component correction_factors (per_half_molar)');
LEGEND_CONSTANTS(:,138) = strpad('Icorr in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,139) = strpad('I1 in component correction_factors (molar)');
LEGEND_CONSTANTS(:,140) = strpad('alphadebye in component correction_factors (per_half_molar)');
LEGEND_CONSTANTS(:,141) = strpad('IcorrpKa in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,142) = strpad('TcorrpKa in component correction_factors (mole_per_kilojoule)');
LEGEND_CONSTANTS(:,143) = strpad('RTalpha in component correction_factors (kilojoule_half_liter_per_3_half_mole)');
LEGEND_CONSTANTS(:,144) = strpad('IcorrdeltaGpof in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,145) = strpad('pKak_Pi in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,146) = strpad('deltaH1o_Pi in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,147) = strpad('deltaHmgo_Pi in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,148) = strpad('deltaH1_Pi in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,149) = strpad('deltaHmg_Pi in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,150) = strpad('pKa1_Pi in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,151) = strpad('pKamg_Pi in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,2) = strpad('P_Pi in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,3) = strpad('HPi2 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,4) = strpad('H2Pi1 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,5) = strpad('kPi in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,6) = strpad('mgPi in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,7) = strpad('Navg_Pi in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,8) = strpad('dNavgPidH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,9) = strpad('dNavgPidmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,10) = strpad('dmgPidmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,11) = strpad('dmgPidpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,152) = strpad('NH_HPi2 in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,153) = strpad('deltaGof_HPi2 in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,12) = strpad('deltaGpof_HPi2 in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,154) = strpad('deltaH1o_ATP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,155) = strpad('deltaHmgo_ATP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,156) = strpad('deltaHko_ATP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,157) = strpad('deltaH1_ATP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,158) = strpad('deltaHmg_ATP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,159) = strpad('deltaHk_ATP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,160) = strpad('pKa1_ATP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,161) = strpad('pKamg_ATP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,162) = strpad('pKak_ATP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,13) = strpad('P_ATP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,14) = strpad('ATP4 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,15) = strpad('HATP3 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,16) = strpad('mgATP2 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,17) = strpad('kATP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,18) = strpad('Navg_ATP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,19) = strpad('dNavgATPdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,20) = strpad('dNavgATPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,21) = strpad('dmgATP2dmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,22) = strpad('dmgATP2dpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,163) = strpad('NH_ATP4 in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,164) = strpad('deltaGof_ATP4 in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,23) = strpad('deltaGpof_ATP4 in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,165) = strpad('pKak_ADP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,166) = strpad('deltaH1o_ADP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,167) = strpad('deltaHmgo_ADP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,168) = strpad('deltaH1_ADP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,169) = strpad('deltaHmg_ADP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,170) = strpad('pKa1_ADP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,171) = strpad('pKamg_ADP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,24) = strpad('P_ADP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,25) = strpad('ADP3 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,26) = strpad('HADP2 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,27) = strpad('mgADP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,28) = strpad('kADP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,29) = strpad('Navg_ADP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,30) = strpad('dNavgADPdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,32) = strpad('dNavgADPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,33) = strpad('dmgADPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,34) = strpad('dmgADPdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,172) = strpad('NH_ADP3 in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,173) = strpad('deltaGof_ADP3 in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,35) = strpad('deltaGpof_ADP3 in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,174) = strpad('deltaH1o_AMP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,175) = strpad('deltaHmgo_AMP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,176) = strpad('deltaH1_AMP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,177) = strpad('deltaHmg_AMP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,178) = strpad('pKa1_AMP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,179) = strpad('pKamg_AMP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,36) = strpad('P_AMP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,37) = strpad('AMP2 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,38) = strpad('HAMP1 in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,39) = strpad('mgAMP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,40) = strpad('Navg_AMP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,41) = strpad('dNavgAMPdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,42) = strpad('dNavgAMPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,43) = strpad('dmgAMPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,44) = strpad('dmgAMPdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,180) = strpad('NH_AMP2 in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,181) = strpad('deltaGof_AMP2 in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,45) = strpad('deltaGpof_AMP2 in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,182) = strpad('pKak_PCR in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,183) = strpad('deltaH1o_PCR in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,184) = strpad('deltaHmgo_PCR in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,185) = strpad('deltaH1_PCR in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,186) = strpad('deltaHmg_PCR in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,187) = strpad('pKa1_PCR in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,188) = strpad('pKamg_PCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,46) = strpad('P_PCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,47) = strpad('HPCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,49) = strpad('H2PCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,48) = strpad('kPCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,50) = strpad('mgPCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,51) = strpad('Navg_PCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,52) = strpad('dNavgPCRdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,53) = strpad('dNavgPCRdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,54) = strpad('dmgPCRdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,55) = strpad('dmgPCRdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,189) = strpad('NH_HPCR in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,190) = strpad('pKa1_CR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,56) = strpad('P_CR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,57) = strpad('HCR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,58) = strpad('H2CR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,59) = strpad('Navg_CR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,60) = strpad('dNavgCRdH in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,27) = strpad('dNavgCRdmg in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,191) = strpad('NH_HCR in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,192) = strpad('deltaH1o_G1P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,193) = strpad('deltaHmgo_G1P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,194) = strpad('deltaH1_G1P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,195) = strpad('deltaHmg_G1P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,196) = strpad('pKa1_G1P in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,197) = strpad('pKamg_G1P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,61) = strpad('P_G1P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,62) = strpad('UG1P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,63) = strpad('HG1P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,64) = strpad('mgG1P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,65) = strpad('Navg_G1P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,66) = strpad('dNavgG1PdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,67) = strpad('dNavgG1Pdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,68) = strpad('dmgG1Pdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,69) = strpad('dmgG1PdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,198) = strpad('NH_UG1P in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,199) = strpad('deltaGof_UG1P in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,70) = strpad('deltaGpof_UG1P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,200) = strpad('pKa1_G6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,71) = strpad('P_G6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,72) = strpad('UG6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,73) = strpad('HG6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,74) = strpad('Navg_G6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,75) = strpad('dNavgG6PdH in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,28) = strpad('dNavgG6Pdmg in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,201) = strpad('NH_UG6P in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,202) = strpad('deltaGof_UG6P in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,76) = strpad('deltaGpof_UG6P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,203) = strpad('pKa1_F6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,77) = strpad('P_F6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,78) = strpad('UF6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,79) = strpad('HF6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,80) = strpad('Navg_F6P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,81) = strpad('dNavgF6PdH in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,29) = strpad('dNavgF6Pdmg in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,204) = strpad('NH_UF6P in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,205) = strpad('deltaGof_UF6P in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,82) = strpad('deltaGpof_UF6P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,206) = strpad('pKa1_FDP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,207) = strpad('pKa2_FDP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,208) = strpad('pKamg_FDP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,83) = strpad('P_FDP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,84) = strpad('UFDP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,85) = strpad('HFDP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,86) = strpad('H2FDP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,87) = strpad('mgFDP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,88) = strpad('Navg_FDP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,89) = strpad('dNavgFDPdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,90) = strpad('dNavgFDPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,91) = strpad('dmgFDPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,92) = strpad('dmgFDPdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,209) = strpad('NH_UFDP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,210) = strpad('deltaGof_UFDP in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,93) = strpad('deltaGpof_UFDP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,211) = strpad('pKa1_GAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,94) = strpad('P_GAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,95) = strpad('UGAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,96) = strpad('HGAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,97) = strpad('Navg_GAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,98) = strpad('dNavgGAPdH in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,30) = strpad('dNavgGAPdmg in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,212) = strpad('NH_UGAP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,213) = strpad('deltaGof_UGAP in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,99) = strpad('deltaGpof_UGAP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,214) = strpad('pKamg_G3P in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,215) = strpad('deltaH1o_G3P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,216) = strpad('deltaH1_G3P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,217) = strpad('pKa1_G3P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,100) = strpad('P_G3P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,101) = strpad('UG3P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,102) = strpad('HG3P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,103) = strpad('mgG3P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,104) = strpad('Navg_G3P in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,105) = strpad('dNavgG3PdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,106) = strpad('dNavgG3Pdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,107) = strpad('dmgG3Pdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,108) = strpad('dmgG3PdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,218) = strpad('NH_UG3P in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,219) = strpad('deltaGof_UG3P in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,109) = strpad('deltaGpof_UG3P in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,220) = strpad('pKa1_DHAP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,221) = strpad('pKamg_DHAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,110) = strpad('P_DHAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,111) = strpad('UDHAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,112) = strpad('HDHAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,113) = strpad('mgDHAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,114) = strpad('Navg_DHAP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,115) = strpad('dNavgDHAPdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,116) = strpad('dNavgDHAPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,117) = strpad('dmgDHAPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,118) = strpad('dmgDHAPdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,222) = strpad('NH_UDHAP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,223) = strpad('deltaGof_UDHAP in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,119) = strpad('deltaGpof_UDHAP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,224) = strpad('pKa1_13DPG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,120) = strpad('P_13DPG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,121) = strpad('U13DPG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,122) = strpad('H13DPG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,123) = strpad('Navg_13DPG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,124) = strpad('dNavg13DPGdH in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,31) = strpad('dNavg13DPGdmg in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,225) = strpad('NH_U13DPG in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,226) = strpad('deltaGof_U13DPG in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,125) = strpad('deltaGpof_U13DPG in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,227) = strpad('pKa1_3PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,126) = strpad('P_3PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,127) = strpad('U3PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,128) = strpad('H3PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,129) = strpad('Navg_3PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,130) = strpad('dNavg3PGdH in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,32) = strpad('dNavg3PGdmg in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,228) = strpad('NH_U3PG in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,229) = strpad('deltaGof_U3PG in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,131) = strpad('deltaGpof_U3PG in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,230) = strpad('pKa1_2PG in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,231) = strpad('pKamg_2PG in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,232) = strpad('pKak_2PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,132) = strpad('P_2PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,133) = strpad('U2PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,134) = strpad('H2PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,136) = strpad('mg2PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,135) = strpad('k2PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,137) = strpad('Navg_2PG in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,138) = strpad('dNavg2PGdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,139) = strpad('dNavg2PGdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,140) = strpad('dmg2PGdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,141) = strpad('dmg2PGdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,233) = strpad('NH_U2PG in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,234) = strpad('deltaGof_U2PG in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,142) = strpad('deltaGpof_U2PG in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,235) = strpad('pKa1_PEP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,236) = strpad('pKamg_PEP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,237) = strpad('pKak_PEP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,143) = strpad('P_PEP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,144) = strpad('UPEP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,145) = strpad('HPEP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,146) = strpad('kPEP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,147) = strpad('mgPEP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,148) = strpad('Navg_PEP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,149) = strpad('dNavgPEPdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,150) = strpad('dNavgPEPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,151) = strpad('dmgPEPdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,152) = strpad('dmgPEPdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,238) = strpad('NH_UPEP in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,239) = strpad('deltaGof_UPEP in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,153) = strpad('deltaGpof_UPEP in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,240) = strpad('pKa1_PYR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,154) = strpad('P_PYR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,155) = strpad('UPYR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,156) = strpad('HPYR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,157) = strpad('Navg_PYR in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,158) = strpad('dNavgPYRdH in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,33) = strpad('dNavgPYRdmg in component correction_factors (per_molar)');
LEGEND_CONSTANTS(:,241) = strpad('NH_UPYR in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,242) = strpad('deltaGof_UPYR in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,159) = strpad('deltaGpof_UPYR in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,243) = strpad('pKamg_LAC in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,244) = strpad('deltaH1o_LAC in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,245) = strpad('deltaH1_LAC in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,246) = strpad('pKa1_LAC in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,160) = strpad('P_LAC in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,161) = strpad('ULAC in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,162) = strpad('HLAC in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,163) = strpad('mgLAC in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,164) = strpad('Navg_LAC in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,165) = strpad('dNavgLACdH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,166) = strpad('dNavgLACdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,167) = strpad('dmgLACdmg in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,168) = strpad('dmgLACdpH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,247) = strpad('NH_ULAC in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,248) = strpad('deltaGof_ULAC in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,169) = strpad('deltaGpof_ULAC in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,249) = strpad('dNH_GLY in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,170) = strpad('deltaGpo_GLY in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,250) = strpad('NH_NAD in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,251) = strpad('deltaGof_NAD in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,171) = strpad('deltaGpof_NAD in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,252) = strpad('NH_NADH in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,253) = strpad('deltaGof_NADH in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,172) = strpad('deltaGpof_NADH in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,254) = strpad('NH_H2O in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,255) = strpad('deltaGof_H2O in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,173) = strpad('deltaGpof_H2O in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,256) = strpad('NH_H in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,257) = strpad('deltaGof_H in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,174) = strpad('deltaGpof_H in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,175) = strpad('deltaH_CK in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,258) = strpad('Kref_CK in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,259) = strpad('deltaHo_CKo in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,260) = strpad('deltaH1_CK in component correction_factors (kilojoule_per_mole)');
LEGEND_CONSTANTS(:,261) = strpad('Kref_CKI in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,262) = strpad('Kref_CKT in component correction_factors (dimensionless)');
LEGEND_CONSTANTS(:,263) = strpad('deltaGpo_CK in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,177) = strpad('Kapp_CK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,179) = strpad('deltaH_ADK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,180) = strpad('deltaGpo_ADK in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,181) = strpad('Kapp_ADK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,182) = strpad('deltaH_GP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,183) = strpad('deltaGpo_GP in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,184) = strpad('Kapp_GP in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,185) = strpad('deltaH_PGLM in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,186) = strpad('deltaGpo_PGLM in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,187) = strpad('Kapp_PGLM in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,188) = strpad('deltaH_PGI in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,189) = strpad('deltaGpo_PGI in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,190) = strpad('Kapp_PGI in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,191) = strpad('deltaH_PFK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,192) = strpad('deltaGpo_PFK in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,193) = strpad('Kapp_PFK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,194) = strpad('deltaH_ALD in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,195) = strpad('deltaGpo_ALD in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,196) = strpad('Kapp_ALD in component correction_factors (molar)');
LEGEND_ALGEBRAIC(:,197) = strpad('deltaH_TPI in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,198) = strpad('deltaGpo_TPI in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,199) = strpad('Kapp_TPI in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,200) = strpad('deltaH_GAPDH in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,201) = strpad('deltaGpo_GAPDH in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,202) = strpad('Kapp_GAPDH in component correction_factors (per_molar)');
LEGEND_ALGEBRAIC(:,203) = strpad('deltaH_G3PDH in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,204) = strpad('deltaGpo_G3PDH in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,205) = strpad('Kapp_G3PDH in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,206) = strpad('deltaH_PGK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,207) = strpad('deltaGpo_PGK in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,208) = strpad('Kapp_PGK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,209) = strpad('deltaH_PGM in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,210) = strpad('deltaGpo_PGM in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,211) = strpad('Kapp_PGM in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,212) = strpad('deltaH_ENOL in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,213) = strpad('deltaGpo_ENOL in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,214) = strpad('Kapp_ENOL in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,215) = strpad('deltaH_PK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,216) = strpad('deltaGpo_PK in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,217) = strpad('Kapp_PK in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,218) = strpad('deltaH_LDH in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,219) = strpad('deltaGpo_LDH in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,220) = strpad('Kapp_LDH in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,31) = strpad('deltaH_ATPase in component correction_factors (dimensionless)');
LEGEND_ALGEBRAIC(:,176) = strpad('deltaGpo_ATPase in component correction_factors (kilojoule_per_mole)');
LEGEND_ALGEBRAIC(:,178) = strpad('Kapp_ATPase in component correction_factors (dimensionless)');
LEGEND_STATES(:,2) = strpad('Mg in component differential_equations (molar)');
LEGEND_CONSTANTS(:,264) = strpad('Vfgly in component glycogen_phosphorylase (molar_per_minute)');
LEGEND_CONSTANTS(:,39) = strpad('expno in component glycogen_phosphorylase (dimensionless)');
LEGEND_ALGEBRAIC(:,221) = strpad('fracA in component glycogen_phosphorylase (dimensionless)');
LEGEND_CONSTANTS(:,40) = strpad('KgpA_glyf in component glycogen_phosphorylase (molar)');
LEGEND_CONSTANTS(:,265) = strpad('KgpA_pi in component glycogen_phosphorylase (molar)');
LEGEND_CONSTANTS(:,266) = strpad('KgpA_igly in component glycogen_phosphorylase (molar)');
LEGEND_CONSTANTS(:,41) = strpad('KgpA_ipi in component glycogen_phosphorylase (molar)');
LEGEND_CONSTANTS(:,267) = strpad('KgpA_glyb in component glycogen_phosphorylase (molar)');
LEGEND_CONSTANTS(:,42) = strpad('KgpA_g1p in component glycogen_phosphorylase (molar)');
LEGEND_CONSTANTS(:,268) = strpad('KgpA_ig1p in component glycogen_phosphorylase (molar)');
LEGEND_ALGEBRAIC(:,222) = strpad('Dglya in component glycogen_phosphorylase (dimensionless)');
LEGEND_ALGEBRAIC(:,223) = strpad('pa in component glycogen_phosphorylase (dimensionless)');
LEGEND_ALGEBRAIC(:,224) = strpad('VbglyA in component glycogen_phosphorylase (molar_per_minute)');
LEGEND_ALGEBRAIC(:,225) = strpad('glyAF in component glycogen_phosphorylase (per_minute)');
LEGEND_ALGEBRAIC(:,226) = strpad('glyAR in component glycogen_phosphorylase (per_minute)');
LEGEND_ALGEBRAIC(:,227) = strpad('flux_GPa in component glycogen_phosphorylase (molar_per_minute)');
LEGEND_STATES(:,3) = strpad('G1P in component differential_equations (molar)');
LEGEND_STATES(:,4) = strpad('Pi in component differential_equations (molar)');
LEGEND_STATES(:,5) = strpad('Gly in component differential_equations (molar)');
LEGEND_ALGEBRAIC(:,228) = strpad('fracB in component glycogen_phosphorylase_B (dimensionless)');
LEGEND_CONSTANTS(:,269) = strpad('KgpB_pi in component glycogen_phosphorylase_B (molar)');
LEGEND_CONSTANTS(:,43) = strpad('KgpB_ipi in component glycogen_phosphorylase_B (molar)');
LEGEND_CONSTANTS(:,270) = strpad('KgpB_iglyf in component glycogen_phosphorylase_B (molar)');
LEGEND_CONSTANTS(:,271) = strpad('KgpB_g1p in component glycogen_phosphorylase_B (molar)');
LEGEND_CONSTANTS(:,44) = strpad('KgpB_ig1p in component glycogen_phosphorylase_B (molar)');
LEGEND_CONSTANTS(:,272) = strpad('KgpB_iglyb in component glycogen_phosphorylase_B (molar)');
LEGEND_CONSTANTS(:,45) = strpad('Kgp_amp in component glycogen_phosphorylase_B (molar)');
LEGEND_CONSTANTS(:,46) = strpad('interactioncoeff in component glycogen_phosphorylase_B (dimensionless)');
LEGEND_CONSTANTS(:,47) = strpad('nH in component glycogen_phosphorylase_B (dimensionless)');
LEGEND_ALGEBRAIC(:,229) = strpad('M in component glycogen_phosphorylase_B (dimensionless)');
LEGEND_ALGEBRAIC(:,230) = strpad('Dglyb in component glycogen_phosphorylase_B (dimensionless)');
LEGEND_ALGEBRAIC(:,231) = strpad('pb in component glycogen_phosphorylase_B (dimensionless)');
LEGEND_ALGEBRAIC(:,232) = strpad('VbglyB in component glycogen_phosphorylase_B (molar_per_minute)');
LEGEND_ALGEBRAIC(:,233) = strpad('glyBF in component glycogen_phosphorylase_B (per_minute)');
LEGEND_ALGEBRAIC(:,234) = strpad('glyBR in component glycogen_phosphorylase_B (per_minute)');
LEGEND_ALGEBRAIC(:,235) = strpad('flux_GPb in component glycogen_phosphorylase_B (molar_per_minute)');
LEGEND_STATES(:,6) = strpad('AMP in component differential_equations (molar)');
LEGEND_CONSTANTS(:,273) = strpad('Vffpglm in component PGLM (molar_per_minute)');
LEGEND_CONSTANTS(:,274) = strpad('Kpglm_g1p in component PGLM (molar)');
LEGEND_CONSTANTS(:,275) = strpad('Kpglm_g6p in component PGLM (molar)');
LEGEND_ALGEBRAIC(:,236) = strpad('Vfpglm in component PGLM (molar_per_minute)');
LEGEND_ALGEBRAIC(:,238) = strpad('Vbpglm in component PGLM (molar_per_minute)');
LEGEND_ALGEBRAIC(:,239) = strpad('v_PGLM in component PGLM (molar_per_minute)');
LEGEND_STATES(:,7) = strpad('G6P in component differential_equations (molar)');
LEGEND_CONSTANTS(:,276) = strpad('Vbbpgi in component PGI (molar_per_minute)');
LEGEND_CONSTANTS(:,277) = strpad('Kpgi_g6p in component PGI (molar)');
LEGEND_CONSTANTS(:,278) = strpad('Kpgi_f6p in component PGI (molar)');
LEGEND_ALGEBRAIC(:,240) = strpad('Vbpgi in component PGI (molar_per_minute)');
LEGEND_ALGEBRAIC(:,241) = strpad('Vfpgi in component PGI (molar_per_minute)');
LEGEND_ALGEBRAIC(:,242) = strpad('v_PGI in component PGI (molar_per_minute)');
LEGEND_STATES(:,8) = strpad('F6P in component differential_equations (molar)');
LEGEND_CONSTANTS(:,279) = strpad('Vffpfk in component PFK (molar_per_minute)');
LEGEND_CONSTANTS(:,280) = strpad('Kpfk_f6p in component PFK (molar)');
LEGEND_CONSTANTS(:,281) = strpad('Kpfk_f6pT in component PFK (molar)');
LEGEND_CONSTANTS(:,282) = strpad('Kpfk_atp in component PFK (molar)');
LEGEND_CONSTANTS(:,284) = strpad('Kpfk_atpT in component PFK (molar)');
LEGEND_CONSTANTS(:,283) = strpad('Kpfk_fbp in component PFK (molar)');
LEGEND_CONSTANTS(:,48) = strpad('Kpfk_fbpT in component PFK (molar)');
LEGEND_CONSTANTS(:,285) = strpad('Kpfk_adp in component PFK (molar)');
LEGEND_CONSTANTS(:,49) = strpad('Kpfk_adpT in component PFK (molar)');
LEGEND_CONSTANTS(:,50) = strpad('Kpfki in component PFK (molar)');
LEGEND_CONSTANTS(:,51) = strpad('Kmpfk in component PFK (molar)');
LEGEND_CONSTANTS(:,52) = strpad('d in component PFK (dimensionless)');
LEGEND_CONSTANTS(:,53) = strpad('e_ in component PFK (dimensionless)');
LEGEND_CONSTANTS(:,54) = strpad('Lo in component PFK (dimensionless)');
LEGEND_ALGEBRAIC(:,243) = strpad('Vfpfk in component PFK (molar_per_minute)');
LEGEND_ALGEBRAIC(:,245) = strpad('Vbpfk in component PFK (molar_per_minute)');
LEGEND_ALGEBRAIC(:,246) = strpad('L in component PFK (dimensionless)');
LEGEND_CONSTANTS(:,286) = strpad('alpha in component PFK (dimensionless)');
LEGEND_ALGEBRAIC(:,247) = strpad('Delta in component PFK (dimensionless)');
LEGEND_ALGEBRAIC(:,248) = strpad('Deltap in component PFK (dimensionless)');
LEGEND_ALGEBRAIC(:,249) = strpad('v_PFK in component PFK (molar_per_minute)');
LEGEND_STATES(:,9) = strpad('FBP in component differential_equations (molar)');
LEGEND_STATES(:,10) = strpad('ADP in component differential_equations (molar)');
LEGEND_STATES(:,11) = strpad('ATP in component differential_equations (molar)');
LEGEND_CONSTANTS(:,287) = strpad('Vffald in component ALD (molar_per_minute)');
LEGEND_CONSTANTS(:,288) = strpad('Kald_fbp in component ALD (molar)');
LEGEND_CONSTANTS(:,289) = strpad('Kald_dhap in component ALD (molar)');
LEGEND_CONSTANTS(:,290) = strpad('Kald_gap in component ALD (molar)');
LEGEND_ALGEBRAIC(:,250) = strpad('Vfald in component ALD (molar_per_minute)');
LEGEND_ALGEBRAIC(:,252) = strpad('Vbald in component ALD (molar_per_minute)');
LEGEND_ALGEBRAIC(:,253) = strpad('v_ALD in component ALD (molar_per_minute)');
LEGEND_STATES(:,12) = strpad('DHAP in component differential_equations (molar)');
LEGEND_STATES(:,13) = strpad('GAP in component differential_equations (molar)');
LEGEND_CONSTANTS(:,291) = strpad('Vfftpi in component TPI (molar_per_minute)');
LEGEND_CONSTANTS(:,292) = strpad('Ktpi_gap in component TPI (molar)');
LEGEND_CONSTANTS(:,293) = strpad('Ktpi_dhap in component TPI (molar)');
LEGEND_CONSTANTS(:,294) = strpad('Vftpi in component TPI (molar_per_minute)');
LEGEND_ALGEBRAIC(:,254) = strpad('Vbtpi in component TPI (molar_per_minute)');
LEGEND_ALGEBRAIC(:,255) = strpad('v_TPI in component TPI (molar_per_minute)');
LEGEND_CONSTANTS(:,295) = strpad('Vbbg3pdh in component G3PDH (molar_per_minute)');
LEGEND_CONSTANTS(:,296) = strpad('Kg3pdh_g3p in component G3PDH (molar)');
LEGEND_CONSTANTS(:,297) = strpad('Kg3pdh_nad in component G3PDH (molar)');
LEGEND_CONSTANTS(:,298) = strpad('Kg3pdh_dhap in component G3PDH (molar)');
LEGEND_CONSTANTS(:,299) = strpad('Kg3pdh_nadh in component G3PDH (molar)');
LEGEND_ALGEBRAIC(:,256) = strpad('Dg3pdh in component G3PDH (dimensionless)');
LEGEND_CONSTANTS(:,300) = strpad('Vbg3pdh in component G3PDH (molar_per_minute)');
LEGEND_ALGEBRAIC(:,257) = strpad('Vfg3pdh in component G3PDH (molar_per_minute)');
LEGEND_ALGEBRAIC(:,258) = strpad('v_G3PDH in component G3PDH (molar_per_minute)');
LEGEND_STATES(:,14) = strpad('G3P in component differential_equations (molar)');
LEGEND_STATES(:,15) = strpad('NAD in component differential_equations (molar)');
LEGEND_STATES(:,16) = strpad('NADH in component differential_equations (molar)');
LEGEND_CONSTANTS(:,301) = strpad('Vffgad in component GAPDH (molar_per_minute)');
LEGEND_CONSTANTS(:,302) = strpad('Kgapdh_gap in component GAPDH (molar)');
LEGEND_CONSTANTS(:,303) = strpad('Kgapdh_nad in component GAPDH (molar)');
LEGEND_CONSTANTS(:,304) = strpad('Kgapdh_pi in component GAPDH (molar)');
LEGEND_CONSTANTS(:,305) = strpad('Kgapdh_bpg in component GAPDH (molar)');
LEGEND_CONSTANTS(:,306) = strpad('Kgapdh_nadh in component GAPDH (molar)');
LEGEND_ALGEBRAIC(:,259) = strpad('Dgap in component GAPDH (dimensionless)');
LEGEND_ALGEBRAIC(:,260) = strpad('Vfgad in component GAPDH (molar_per_minute)');
LEGEND_ALGEBRAIC(:,261) = strpad('Vbgad in component GAPDH (molar_per_minute)');
LEGEND_ALGEBRAIC(:,262) = strpad('v_GAPDH in component GAPDH (molar_per_minute)');
LEGEND_STATES(:,17) = strpad('BPG in component differential_equations (molar)');
LEGEND_CONSTANTS(:,307) = strpad('Vbbpgk in component PGK (molar_per_minute)');
LEGEND_CONSTANTS(:,308) = strpad('Kpgk_bpg in component PGK (molar)');
LEGEND_CONSTANTS(:,309) = strpad('Kpgk_adp in component PGK (molar)');
LEGEND_CONSTANTS(:,310) = strpad('Kpgk_3pg in component PGK (molar)');
LEGEND_CONSTANTS(:,311) = strpad('Kpgk_atp in component PGK (molar)');
LEGEND_CONSTANTS(:,312) = strpad('Vbpgk in component PGK (molar_per_minute)');
LEGEND_ALGEBRAIC(:,263) = strpad('Vfpgk in component PGK (molar_per_minute)');
LEGEND_ALGEBRAIC(:,265) = strpad('D_PGK in component PGK (dimensionless)');
LEGEND_ALGEBRAIC(:,266) = strpad('v_PGK in component PGK (molar_per_minute)');
LEGEND_STATES(:,18) = strpad('P3G in component differential_equations (molar)');
LEGEND_CONSTANTS(:,313) = strpad('Vffpgm in component PGM (molar_per_minute)');
LEGEND_CONSTANTS(:,314) = strpad('Kpgm_3pg in component PGM (molar)');
LEGEND_CONSTANTS(:,315) = strpad('Kpgm_2pg in component PGM (molar)');
LEGEND_ALGEBRAIC(:,267) = strpad('Vfpgm in component PGM (molar_per_minute)');
LEGEND_ALGEBRAIC(:,269) = strpad('Vbpgm in component PGM (molar_per_minute)');
LEGEND_ALGEBRAIC(:,270) = strpad('v_PGM in component PGM (molar_per_minute)');
LEGEND_STATES(:,19) = strpad('P2G in component differential_equations (molar)');
LEGEND_CONSTANTS(:,316) = strpad('Vffen in component ENOL (molar_per_minute)');
LEGEND_CONSTANTS(:,317) = strpad('Ken_2pg in component ENOL (molar)');
LEGEND_CONSTANTS(:,318) = strpad('Ken_pep in component ENOL (molar)');
LEGEND_CONSTANTS(:,319) = strpad('Vfen in component ENOL (molar_per_minute)');
LEGEND_ALGEBRAIC(:,271) = strpad('Vben in component ENOL (molar_per_minute)');
LEGEND_ALGEBRAIC(:,273) = strpad('v_ENOL in component ENOL (molar_per_minute)');
LEGEND_STATES(:,20) = strpad('PEP in component differential_equations (molar)');
LEGEND_CONSTANTS(:,320) = strpad('Vffpk in component PK (molar_per_minute)');
LEGEND_CONSTANTS(:,321) = strpad('Kpk_pep in component PK (molar)');
LEGEND_CONSTANTS(:,322) = strpad('Kpk_adp in component PK (molar)');
LEGEND_CONSTANTS(:,323) = strpad('Kpk_pyr in component PK (molar)');
LEGEND_CONSTANTS(:,324) = strpad('Kpk_atp in component PK (molar)');
LEGEND_ALGEBRAIC(:,274) = strpad('Vfpk in component PK (molar_per_minute)');
LEGEND_ALGEBRAIC(:,275) = strpad('Vbpk in component PK (molar_per_minute)');
LEGEND_ALGEBRAIC(:,276) = strpad('v_PK in component PK (molar_per_minute)');
LEGEND_STATES(:,21) = strpad('PYR in component differential_equations (molar)');
LEGEND_CONSTANTS(:,325) = strpad('Vffldh in component LDH (molar_per_minute)');
LEGEND_CONSTANTS(:,326) = strpad('Kldh_pyr in component LDH (molar)');
LEGEND_CONSTANTS(:,327) = strpad('Kldh_nadh in component LDH (molar)');
LEGEND_CONSTANTS(:,328) = strpad('Kldh_lac in component LDH (molar)');
LEGEND_CONSTANTS(:,329) = strpad('Kldh_nad in component LDH (molar)');
LEGEND_ALGEBRAIC(:,277) = strpad('Vfldh in component LDH (molar_per_minute)');
LEGEND_ALGEBRAIC(:,278) = strpad('Vbldh in component LDH (molar_per_minute)');
LEGEND_ALGEBRAIC(:,279) = strpad('v_LDH in component LDH (molar_per_minute)');
LEGEND_STATES(:,22) = strpad('LAC in component differential_equations (molar)');
LEGEND_ALGEBRAIC(:,280) = strpad('VmaxATPase in component ATPase (molar_per_minute)');
LEGEND_CONSTANTS(:,55) = strpad('Katp_ATPase in component ATPase (molar)');
LEGEND_ALGEBRAIC(:,285) = strpad('ATPase in component ATPase (molar_per_minute)');
LEGEND_CONSTANTS(:,330) = strpad('VforCK in component creatine_kinase (molar_per_minute)');
LEGEND_CONSTANTS(:,331) = strpad('Kck_pcr in component creatine_kinase (molar)');
LEGEND_CONSTANTS(:,332) = strpad('Kck_iatp in component creatine_kinase (molar)');
LEGEND_CONSTANTS(:,333) = strpad('Kck_iadp in component creatine_kinase (molar)');
LEGEND_CONSTANTS(:,56) = strpad('Kck_ipcr in component creatine_kinase (molar)');
LEGEND_CONSTANTS(:,334) = strpad('Kck_cr in component creatine_kinase (molar)');
LEGEND_ALGEBRAIC(:,286) = strpad('VrevCK in component creatine_kinase (molar_per_minute)');
LEGEND_ALGEBRAIC(:,287) = strpad('CK in component creatine_kinase (molar_per_minute)');
LEGEND_STATES(:,23) = strpad('Cr in component differential_equations (molar)');
LEGEND_STATES(:,24) = strpad('PCr in component differential_equations (molar)');
LEGEND_CONSTANTS(:,335) = strpad('Vfadk in component adenylate_kinase (molar_per_minute)');
LEGEND_CONSTANTS(:,336) = strpad('Kadk_amp in component adenylate_kinase (molar)');
LEGEND_CONSTANTS(:,337) = strpad('Kadk_atp in component adenylate_kinase (molar)');
LEGEND_CONSTANTS(:,338) = strpad('Kadk_adp in component adenylate_kinase (molar)');
LEGEND_ALGEBRAIC(:,288) = strpad('Vbadk in component adenylate_kinase (molar_per_minute)');
LEGEND_ALGEBRAIC(:,291) = strpad('ADK in component adenylate_kinase (molar_per_minute)');
LEGEND_CONSTANTS(:,339) = strpad('carnosine in component buffer_capacity (molar)');
LEGEND_CONSTANTS(:,340) = strpad('tris in component buffer_capacity (molar)');
LEGEND_CONSTANTS(:,341) = strpad('acetate in component buffer_capacity (molar)');
LEGEND_ALGEBRAIC(:,292) = strpad('bufcapfixed in component buffer_capacity (molar)');
LEGEND_ALGEBRAIC(:,293) = strpad('bufcapmetab in component buffer_capacity (molar)');
LEGEND_ALGEBRAIC(:,294) = strpad('protons_consumed in component buffer_capacity (molar_per_minute)');
LEGEND_ALGEBRAIC(:,289) = strpad('CKprtflux in component buffer_capacity (molar_per_minute)');
LEGEND_ALGEBRAIC(:,281) = strpad('glycprtflux in component buffer_capacity (molar_per_minute)');
LEGEND_ALGEBRAIC(:,295) = strpad('pHODEterm1 in component buffer_capacity (per_minute)');
LEGEND_ALGEBRAIC(:,296) = strpad('pHODEterm2 in component buffer_capacity (per_molar)');
LEGEND_ALGEBRAIC(:,297) = strpad('denom_mgODE in component buffer_capacity (dimensionless)');
LEGEND_ALGEBRAIC(:,298) = strpad('RHSterm1_mgODE in component buffer_capacity (molar)');
LEGEND_ALGEBRAIC(:,299) = strpad('denomMgpHODE in component buffer_capacity (dimensionless)');
LEGEND_ALGEBRAIC(:,312) = strpad('RHSterm2_mgODE in component differential_equations (molar_per_minute)');
LEGEND_CONSTANTS(:,34) = strpad('fixmg in component differential_equations (dimensionless)');
LEGEND_CONSTANTS(:,342) = strpad('fixpH in component differential_equations (dimensionless)');
LEGEND_STATES(:,25) = strpad('protonload in component differential_equations (molar)');
LEGEND_ALGEBRAIC(:,301) = strpad('dATPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,302) = strpad('dADPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,303) = strpad('dAMPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,307) = strpad('dDHAPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,306) = strpad('dFBPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,305) = strpad('dG1Pdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,308) = strpad('dG3Pdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,311) = strpad('dLACdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,309) = strpad('dP2Gdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,300) = strpad('dPCrdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,310) = strpad('dPEPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,304) = strpad('dPidt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,290) = strpad('dCrdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,282) = strpad('dNADdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,283) = strpad('dNADHdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,237) = strpad('dGlydt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,244) = strpad('dG6Pdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,251) = strpad('dF6Pdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,264) = strpad('dGAPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,268) = strpad('dBGPdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,272) = strpad('dP3Gdt in component differential_equations (molar_per_minute)');
LEGEND_ALGEBRAIC(:,284) = strpad('dPYRdt in component differential_equations (molar_per_minute)');
LEGEND_RATES(:,24) = strpad('d/dt PCr in component differential_equations (molar)');
LEGEND_RATES(:,23) = strpad('d/dt Cr in component differential_equations (molar)');
LEGEND_RATES(:,15) = strpad('d/dt NAD in component differential_equations (molar)');
LEGEND_RATES(:,16) = strpad('d/dt NADH in component differential_equations (molar)');
LEGEND_RATES(:,11) = strpad('d/dt ATP in component differential_equations (molar)');
LEGEND_RATES(:,10) = strpad('d/dt ADP in component differential_equations (molar)');
LEGEND_RATES(:,6) = strpad('d/dt AMP in component differential_equations (molar)');
LEGEND_RATES(:,4) = strpad('d/dt Pi in component differential_equations (molar)');
LEGEND_RATES(:,5) = strpad('d/dt Gly in component differential_equations (molar)');
LEGEND_RATES(:,3) = strpad('d/dt G1P in component differential_equations (molar)');
LEGEND_RATES(:,7) = strpad('d/dt G6P in component differential_equations (molar)');
LEGEND_RATES(:,8) = strpad('d/dt F6P in component differential_equations (molar)');
LEGEND_RATES(:,9) = strpad('d/dt FBP in component differential_equations (molar)');
LEGEND_RATES(:,12) = strpad('d/dt DHAP in component differential_equations (molar)');
LEGEND_RATES(:,14) = strpad('d/dt G3P in component differential_equations (molar)');
LEGEND_RATES(:,13) = strpad('d/dt GAP in component differential_equations (molar)');
LEGEND_RATES(:,17) = strpad('d/dt BPG in component differential_equations (molar)');
LEGEND_RATES(:,18) = strpad('d/dt P3G in component differential_equations (molar)');
LEGEND_RATES(:,19) = strpad('d/dt P2G in component differential_equations (molar)');
LEGEND_RATES(:,20) = strpad('d/dt PEP in component differential_equations (molar)');
LEGEND_RATES(:,21) = strpad('d/dt PYR in component differential_equations (molar)');
LEGEND_RATES(:,22) = strpad('d/dt LAC in component differential_equations (molar)');
LEGEND_RATES(:,2) = strpad('d/dt Mg in component differential_equations (molar)');
LEGEND_RATES(:,1) = strpad('d/dt pH_calc in component differential_equations (dimensionless)');
LEGEND_RATES(:,25) = strpad('d/dt protonload in component differential_equations (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 = [];
STATES(:,1) = 7.8;
CONSTANTS(:,1) = 8.314e-3;
CONSTANTS(:,2) = 298.15;
STATES(:,2) = 5.132658807e-4;
STATES(:,3) = 1e-9;
STATES(:,4) = 0.03;
STATES(:,5) = 0.04;
STATES(:,6) = 1e-9;
STATES(:,7) = 1e-9;
STATES(:,8) = 1e-9;
STATES(:,9) = 1e-9;
STATES(:,10) = 1e-9;
STATES(:,11) = 0.005;
STATES(:,12) = 1e-9;
STATES(:,13) = 1e-9;
STATES(:,14) = 1e-9;
STATES(:,15) = 0.0005;
STATES(:,16) = 1e-9;
STATES(:,17) = 1e-9;
STATES(:,18) = 1e-9;
STATES(:,19) = 1e-9;
STATES(:,20) = 1e-9;
STATES(:,21) = 1e-9;
STATES(:,22) = 1e-9;
STATES(:,23) = 0.029999999;
STATES(:,24) = 1e-9;
STATES(:,25) = 0;
CONSTANTS(:,3) = 1.00000;
CONSTANTS(:,4) = 7.40000;
CONSTANTS(:,5) = 0.0500000;
CONSTANTS(:,6) = 0.400000;
CONSTANTS(:,7) = 0.00170000;
CONSTANTS(:,8) = 0.00470000;
CONSTANTS(:,9) = 0.00270000;
CONSTANTS(:,10) = 0.00460000;
CONSTANTS(:,11) = 0.00740000;
CONSTANTS(:,12) = 0.00266055;
CONSTANTS(:,13) = 0.0200000;
CONSTANTS(:,14) = 1.75052;
CONSTANTS(:,15) = 0.00402000;
CONSTANTS(:,16) = 0.00270000;
CONSTANTS(:,17) = 0.000870000;
CONSTANTS(:,18) = 6.00000e-05;
CONSTANTS(:,19) = 0.0100000;
CONSTANTS(:,20) = 0.0100000;
CONSTANTS(:,21) = 13.0000;
CONSTANTS(:,22) = 0.00390000;
CONSTANTS(:,23) = 0.00000;
CONSTANTS(:,24) = 0.000100000;
CONSTANTS(:,25) = 0.00500000;
CONSTANTS(:,26) = 29.0000;
CONSTANTS(:,27) = 0.00000;
CONSTANTS(:,28) = 0.00000;
CONSTANTS(:,29) = 0.00000;
CONSTANTS(:,30) = 0.00000;
CONSTANTS(:,31) = 0.00000;
CONSTANTS(:,32) = 0.00000;
CONSTANTS(:,33) = 0.00000;
CONSTANTS(:,34) = 1.00000;
CONSTANTS(:,35) = CONSTANTS(:,3);
CONSTANTS(:,36) = CONSTANTS(:,4);
CONSTANTS(:,37) = 0.00400000;
CONSTANTS(:,38) = CONSTANTS(:,25);
CONSTANTS(:,39) = CONSTANTS(:,26);
CONSTANTS(:,40) = CONSTANTS(:,7);
CONSTANTS(:,41) = CONSTANTS(:,8);
CONSTANTS(:,42) = CONSTANTS(:,9);
CONSTANTS(:,43) = CONSTANTS(:,10);
CONSTANTS(:,44) = CONSTANTS(:,11);
CONSTANTS(:,45) = CONSTANTS(:,12);
CONSTANTS(:,46) = CONSTANTS(:,13);
CONSTANTS(:,47) = CONSTANTS(:,14);
CONSTANTS(:,48) = CONSTANTS(:,15);
CONSTANTS(:,49) = CONSTANTS(:,16);
CONSTANTS(:,50) = CONSTANTS(:,17);
CONSTANTS(:,51) = CONSTANTS(:,18);
CONSTANTS(:,52) = CONSTANTS(:,19);
CONSTANTS(:,53) = CONSTANTS(:,20);
CONSTANTS(:,54) = CONSTANTS(:,21);
CONSTANTS(:,55) = CONSTANTS(:,24);
CONSTANTS(:,56) = CONSTANTS(:,22);
CONSTANTS(:,57) = 0.00200000./1.50000;
CONSTANTS(:,58) = 0.000150000;
CONSTANTS(:,59) = 0.0101000;
CONSTANTS(:,60) = 0.000200000;
CONSTANTS(:,61) = 0.0150000;
CONSTANTS(:,62) = 0.00150000;
CONSTANTS(:,63) = 0.00440000;
CONSTANTS(:,64) = 0.480000;
CONSTANTS(:,65) = 6.30000e-05;
CONSTANTS(:,66) = 3.00000e-05;
CONSTANTS(:,67) = 0.880000;
CONSTANTS(:,68) = 0.000480000;
CONSTANTS(:,69) = 0.000119000;
CONSTANTS(:,70) = 0.0560000;
CONSTANTS(:,71) = 0.000180000;
CONSTANTS(:,72) = 0.0200000;
CONSTANTS(:,73) = 8.00000e-05;
CONSTANTS(:,74) = 0.000250000;
CONSTANTS(:,75) = 0.00402000;
CONSTANTS(:,76) = 0.00270000;
CONSTANTS(:,77) = 0.0106591;
CONSTANTS(:,78) = 5.00000e-05;
CONSTANTS(:,79) = 0.00200000;
CONSTANTS(:,80) = 0.00100000;
CONSTANTS(:,81) = 12.0000;
CONSTANTS(:,82) = 0.000320000;
CONSTANTS(:,83) = 0.000610000;
CONSTANTS(:,84) = 0.0825000;
CONSTANTS(:,85) = 0.000180000;
CONSTANTS(:,86) = 1.20000e-05;
CONSTANTS(:,87) = 0.000220000;
CONSTANTS(:,88) = 8.00000e-06;
CONSTANTS(:,89) = 1.26500;
CONSTANTS(:,90) = 2.50000e-06;
CONSTANTS(:,91) = 9.00000e-05;
CONSTANTS(:,92) = 0.000290000;
CONSTANTS(:,93) = 8.00000e-07;
CONSTANTS(:,94) = 3.30000e-06;
CONSTANTS(:,95) = 1.12000;
CONSTANTS(:,96) = 0.00200000;
CONSTANTS(:,97) = 8.00000e-06;
CONSTANTS(:,98) = 0.00120000;
CONSTANTS(:,99) = 0.000350000;
CONSTANTS(:,100) = 1.12000;
CONSTANTS(:,101) = 0.000200000;
CONSTANTS(:,102) = 1.40000e-05;
CONSTANTS(:,103) = 0.192000;
CONSTANTS(:,104) = 0.000100000;
CONSTANTS(:,105) = 0.000370000;
CONSTANTS(:,106) = 1.44000;
CONSTANTS(:,107) = 8.00000e-05;
CONSTANTS(:,108) = 0.000300000;
CONSTANTS(:,109) = 0.00705000;
CONSTANTS(:,110) = 0.00113000;
CONSTANTS(:,111) = 1.92000;
CONSTANTS(:,112) = 0.000335000;
CONSTANTS(:,113) = 2.00000e-06;
CONSTANTS(:,114) = 0.0170000;
CONSTANTS(:,115) = 0.000849000;
CONSTANTS(:,116) = 0.500000;
CONSTANTS(:,117) = 0.00111000;
CONSTANTS(:,118) = 0.00350000;
CONSTANTS(:,119) = 0.000135000;
CONSTANTS(:,120) = 0.00380000;
CONSTANTS(:,121) = 0.880000;
CONSTANTS(:,122) = 0.000320000;
CONSTANTS(:,123) = 0.000270000;
CONSTANTS(:,124) = 0.000350000;
CONSTANTS(:,125) = 0.100000;
CONSTANTS(:,126) = CONSTANTS(:,125);
CONSTANTS(:,127) = 0.0150000;
CONSTANTS(:,128) = 0.0250000;
CONSTANTS(:,129) = 0.0100000;
CONSTANTS(:,130) = 303.150;
CONSTANTS(:,131) = CONSTANTS(:,130);
CONSTANTS(:,132) = 0.0800000;
CONSTANTS(:,133) = 1.00000;
CONSTANTS(:,134) = CONSTANTS(:,132);
CONSTANTS(:,135) = 1.00000;
CONSTANTS(:,136) = 1.47750;
CONSTANTS(:,137) = 1.60000;
CONSTANTS(:,138) = ( CONSTANTS(:,136).*power(CONSTANTS(:,126), 1.0 ./ 2))./(1.00000+ CONSTANTS(:,137).*power(CONSTANTS(:,126), 1.0 ./ 2));
CONSTANTS(:,139) = 0.100000;
CONSTANTS(:,140) = 1.17582;
CONSTANTS(:,141) = ( 1.00000.*CONSTANTS(:,140).*(power(CONSTANTS(:,139), 1.0 ./ 2)./(1.00000+ CONSTANTS(:,137).*power(CONSTANTS(:,139), 1.0 ./ 2)) - power(CONSTANTS(:,126), 1.0 ./ 2)./(1.00000+ CONSTANTS(:,137).*power(CONSTANTS(:,126), 1.0 ./ 2))))./log(10.0000);
CONSTANTS(:,142) = (1.00000./CONSTANTS(:,131) - 1.00000./CONSTANTS(:,2))./( log(10.0000).*CONSTANTS(:,1));
CONSTANTS(:,143) = 2.91482;
CONSTANTS(:,144) = ( CONSTANTS(:,143).*power(CONSTANTS(:,126), 1.0 ./ 2))./(1.00000+ CONSTANTS(:,137).*power(CONSTANTS(:,126), 1.0 ./ 2));
CONSTANTS(:,145) = 0.500000;
CONSTANTS(:,146) = 3.00000;
CONSTANTS(:,147) = - 2.90000;
CONSTANTS(:,148) = CONSTANTS(:,146)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000));
CONSTANTS(:,149) = CONSTANTS(:,147)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(2.00000, 2.00000)) - power(0.00000, 2.00000));
CONSTANTS(:,150) = 6.75000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,148);
CONSTANTS(:,151) = 1.65000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,149);
CONSTANTS(:,152) = 1.00000;
CONSTANTS(:,153) = - 1096.10;
CONSTANTS(:,154) = - 5.00000;
CONSTANTS(:,155) = - 18.0000;
CONSTANTS(:,156) = - 1.00000;
CONSTANTS(:,157) = CONSTANTS(:,154)+ CONSTANTS(:,138).*((power(4.00000, 2.00000)+power(1.00000, 2.00000)) - power(3.00000, 2.00000));
CONSTANTS(:,158) = CONSTANTS(:,155)+ CONSTANTS(:,138).*((power(4.00000, 2.00000)+power(2.00000, 2.00000)) - power(2.00000, 2.00000));
CONSTANTS(:,159) = CONSTANTS(:,156)+ CONSTANTS(:,138).*((power(4.00000, 2.00000)+power(1.00000, 2.00000)) - power(3.00000, 2.00000));
CONSTANTS(:,160) = 6.48000+ (CONSTANTS(:,141)./1.00000).*((power(4.00000, 2.00000)+power(1.00000, 2.00000)) - power(3.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,157);
CONSTANTS(:,161) = 4.19000+ (CONSTANTS(:,141)./1.00000).*((power(4.00000, 2.00000)+power(2.00000, 2.00000)) - power(2.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,158);
CONSTANTS(:,162) = 1.17000+ (CONSTANTS(:,141)./1.00000).*((power(4.00000, 2.00000)+power(1.00000, 2.00000)) - power(3.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,159);
CONSTANTS(:,163) = 12.0000;
CONSTANTS(:,164) = - 2768.10;
CONSTANTS(:,165) = 1.00000;
CONSTANTS(:,166) = - 3.00000;
CONSTANTS(:,167) = - 15.0000;
CONSTANTS(:,168) = CONSTANTS(:,166)+ CONSTANTS(:,138).*((power(3.00000, 2.00000)+power(1.00000, 2.00000)) - power(2.00000, 2.00000));
CONSTANTS(:,169) = CONSTANTS(:,167)+ CONSTANTS(:,138).*((power(3.00000, 2.00000)+power(2.00000, 2.00000)) - power(1.00000, 2.00000));
CONSTANTS(:,170) = 6.38000+ (CONSTANTS(:,141)./1.00000).*((power(3.00000, 2.00000)+power(1.00000, 2.00000)) - power(2.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,168);
CONSTANTS(:,171) = 3.25000+ (CONSTANTS(:,141)./1.00000).*((power(3.00000, 2.00000)+power(2.00000, 2.00000)) - power(1.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,169);
CONSTANTS(:,172) = 12.0000;
CONSTANTS(:,173) = - 1906.13;
CONSTANTS(:,174) = - 3.00000;
CONSTANTS(:,175) = - 7.50000;
CONSTANTS(:,176) = CONSTANTS(:,174)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000));
CONSTANTS(:,177) = CONSTANTS(:,175)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(2.00000, 2.00000)) - power(0.00000, 2.00000));
CONSTANTS(:,178) = 6.29000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,176);
CONSTANTS(:,179) = 1.92000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(2.00000, 2.00000)) - power(0.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,177);
CONSTANTS(:,180) = 12.0000;
CONSTANTS(:,181) = - 1040.45;
CONSTANTS(:,182) = 0.310000;
CONSTANTS(:,183) = 2.66000;
CONSTANTS(:,184) = 8.19000;
CONSTANTS(:,185) = CONSTANTS(:,183)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000));
CONSTANTS(:,186) = CONSTANTS(:,184)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(2.00000, 2.00000)) - power(0.00000, 2.00000));
CONSTANTS(:,187) = 4.50000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,185);
CONSTANTS(:,188) = 1.60000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(2.00000, 2.00000)) - power(0.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,186);
CONSTANTS(:,189) = 8.00000;
CONSTANTS(:,190) = 2.30000;
CONSTANTS(:,191) = 9.00000;
CONSTANTS(:,192) = - 1.70000;
CONSTANTS(:,193) = - 12.0000;
CONSTANTS(:,194) = CONSTANTS(:,192)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000));
CONSTANTS(:,195) = CONSTANTS(:,193)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(2.00000, 2.00000)) - power(0.00000, 2.00000));
CONSTANTS(:,196) = 6.09000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,194);
CONSTANTS(:,197) = 2.48000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(2.00000, 2.00000)) - power(0.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,195);
CONSTANTS(:,198) = 11.0000;
CONSTANTS(:,199) = - 1756.87;
CONSTANTS(:,200) = 6.11000;
CONSTANTS(:,201) = 11.0000;
CONSTANTS(:,202) = - 1763.94;
CONSTANTS(:,203) = 5.89000;
CONSTANTS(:,204) = 11.0000;
CONSTANTS(:,205) = - 1760.80;
CONSTANTS(:,206) = 6.40000;
CONSTANTS(:,207) = 5.92000;
CONSTANTS(:,208) = 2.70000;
CONSTANTS(:,209) = 10.0000;
CONSTANTS(:,210) = - 2601.40;
CONSTANTS(:,211) = 6.45000;
CONSTANTS(:,212) = 5.00000;
CONSTANTS(:,213) = - 1288.60;
CONSTANTS(:,214) = 1.63000;
CONSTANTS(:,215) = - 3.10000;
CONSTANTS(:,216) = CONSTANTS(:,215)+ CONSTANTS(:,138).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000));
CONSTANTS(:,217) = 6.22000+ (CONSTANTS(:,141)./1.00000).*((power(2.00000, 2.00000)+power(1.00000, 2.00000)) - power(1.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,216);
CONSTANTS(:,218) = 7.00000;
CONSTANTS(:,219) = - 1339.25;
CONSTANTS(:,220) = 5.90000;
CONSTANTS(:,221) = 1.57000;
CONSTANTS(:,222) = 5.00000;
CONSTANTS(:,223) = - 1296.26;
CONSTANTS(:,224) = 7.50000;
CONSTANTS(:,225) = 4.00000;
CONSTANTS(:,226) = - 2356.14;
CONSTANTS(:,227) = 6.21000;
CONSTANTS(:,228) = 4.00000;
CONSTANTS(:,229) = - 1502.54;
CONSTANTS(:,230) = 7.00000;
CONSTANTS(:,231) = 2.45000;
CONSTANTS(:,232) = 1.18000;
CONSTANTS(:,233) = 4.00000;
CONSTANTS(:,234) = - 1496.38;
CONSTANTS(:,235) = 6.35000;
CONSTANTS(:,236) = 2.26000;
CONSTANTS(:,237) = 1.08000;
CONSTANTS(:,238) = 2.00000;
CONSTANTS(:,239) = - 1263.65;
CONSTANTS(:,240) = 2.49000;
CONSTANTS(:,241) = 3.00000;
CONSTANTS(:,242) = - 472.270;
CONSTANTS(:,243) = 0.980000;
CONSTANTS(:,244) = - 0.330000;
CONSTANTS(:,245) = CONSTANTS(:,244)+ CONSTANTS(:,138).*((power(1.00000, 2.00000)+power(1.00000, 2.00000)) - power(0.00000, 2.00000));
CONSTANTS(:,246) = 3.67000+ (CONSTANTS(:,141)./1.00000).*((power(1.00000, 2.00000)+power(1.00000, 2.00000)) - power(0.00000, 2.00000))+ CONSTANTS(:,142).*CONSTANTS(:,245);
CONSTANTS(:,247) = 5.00000;
CONSTANTS(:,248) = - 516.720;
CONSTANTS(:,249) = - 10.0000;
CONSTANTS(:,250) = 26.0000;
CONSTANTS(:,251) = 0.00000;
CONSTANTS(:,252) = 27.0000;
CONSTANTS(:,253) = 22.6500;
CONSTANTS(:,254) = 2.00000;
CONSTANTS(:,255) = - 237.190;
CONSTANTS(:,256) = 1.00000;
CONSTANTS(:,257) = 0.00000;
CONSTANTS(:,258) = 2.58000e+08;
CONSTANTS(:,259) = - 17.5500;
CONSTANTS(:,260) = CONSTANTS(:,259)+ CONSTANTS(:,138).*(((power(2.00000, 2.00000)+power(3.00000, 2.00000)+power(1.00000, 2.00000)) - power(4.00000, 2.00000)) - power(0.00000, 2.00000));
CONSTANTS(:,261) = exp(log(CONSTANTS(:,258))+( CONSTANTS(:,140).*power(CONSTANTS(:,126), 1.0 ./ 2).*(((power(2.00000, 2.00000)+power(3.00000, 2.00000)+power(1.00000, 2.00000)) - power(4.00000, 2.00000)) - power(0.00000, 2.00000)))./(1.00000+ CONSTANTS(:,137).*power(CONSTANTS(:,126), 1.0 ./ 2)));
CONSTANTS(:,262) = power(10.0000, arbitrary_log(CONSTANTS(:,261), 10) - CONSTANTS(:,142).*CONSTANTS(:,260));
CONSTANTS(:,263) = - CONSTANTS(:,1).*CONSTANTS(:,131).*log(CONSTANTS(:,262));
CONSTANTS(:,264) = CONSTANTS(:,5);
CONSTANTS(:,265) = CONSTANTS(:,37);
CONSTANTS(:,266) = CONSTANTS(:,57);
CONSTANTS(:,267) = CONSTANTS(:,58);
CONSTANTS(:,268) = CONSTANTS(:,59);
CONSTANTS(:,269) = CONSTANTS(:,60);
CONSTANTS(:,270) = CONSTANTS(:,61);
CONSTANTS(:,271) = CONSTANTS(:,62);
CONSTANTS(:,272) = CONSTANTS(:,63);
CONSTANTS(:,273) = CONSTANTS(:,64);
CONSTANTS(:,274) = CONSTANTS(:,65);
CONSTANTS(:,275) = CONSTANTS(:,66);
CONSTANTS(:,276) = CONSTANTS(:,67);
CONSTANTS(:,277) = CONSTANTS(:,68);
CONSTANTS(:,278) = CONSTANTS(:,69);
CONSTANTS(:,279) = CONSTANTS(:,70);
CONSTANTS(:,280) = CONSTANTS(:,71);
CONSTANTS(:,281) = CONSTANTS(:,72);
CONSTANTS(:,282) = CONSTANTS(:,73);
CONSTANTS(:,283) = CONSTANTS(:,75);
CONSTANTS(:,284) = CONSTANTS(:,74);
CONSTANTS(:,285) = CONSTANTS(:,76);
CONSTANTS(:,286) = ( CONSTANTS(:,280).*CONSTANTS(:,282))./( CONSTANTS(:,281).*CONSTANTS(:,284));
CONSTANTS(:,287) = CONSTANTS(:,77);
CONSTANTS(:,288) = CONSTANTS(:,78);
CONSTANTS(:,289) = CONSTANTS(:,79);
CONSTANTS(:,290) = CONSTANTS(:,80);
CONSTANTS(:,291) = CONSTANTS(:,81);
CONSTANTS(:,292) = CONSTANTS(:,82);
CONSTANTS(:,293) = CONSTANTS(:,83);
CONSTANTS(:,294) = CONSTANTS(:,291);
CONSTANTS(:,295) = CONSTANTS(:,84);
CONSTANTS(:,296) = CONSTANTS(:,85);
CONSTANTS(:,297) = CONSTANTS(:,86);
CONSTANTS(:,298) = CONSTANTS(:,87);
CONSTANTS(:,299) = CONSTANTS(:,88);
CONSTANTS(:,300) = CONSTANTS(:,295);
CONSTANTS(:,301) = CONSTANTS(:,89);
CONSTANTS(:,302) = CONSTANTS(:,90);
CONSTANTS(:,303) = CONSTANTS(:,91);
CONSTANTS(:,304) = CONSTANTS(:,92);
CONSTANTS(:,305) = CONSTANTS(:,93);
CONSTANTS(:,306) = CONSTANTS(:,94);
CONSTANTS(:,307) = CONSTANTS(:,95);
CONSTANTS(:,308) = CONSTANTS(:,96);
CONSTANTS(:,309) = CONSTANTS(:,97);
CONSTANTS(:,310) = CONSTANTS(:,98);
CONSTANTS(:,311) = CONSTANTS(:,99);
CONSTANTS(:,312) = CONSTANTS(:,307);
CONSTANTS(:,313) = CONSTANTS(:,100);
CONSTANTS(:,314) = CONSTANTS(:,101);
CONSTANTS(:,315) = CONSTANTS(:,102);
CONSTANTS(:,316) = CONSTANTS(:,103);
CONSTANTS(:,317) = CONSTANTS(:,104);
CONSTANTS(:,318) = CONSTANTS(:,105);
CONSTANTS(:,319) = CONSTANTS(:,316);
CONSTANTS(:,320) = CONSTANTS(:,106);
CONSTANTS(:,321) = CONSTANTS(:,107);
CONSTANTS(:,322) = CONSTANTS(:,108);
CONSTANTS(:,323) = CONSTANTS(:,109);
CONSTANTS(:,324) = CONSTANTS(:,110);
CONSTANTS(:,325) = CONSTANTS(:,111);
CONSTANTS(:,326) = CONSTANTS(:,112);
CONSTANTS(:,327) = CONSTANTS(:,113);
CONSTANTS(:,328) = CONSTANTS(:,114);
CONSTANTS(:,329) = CONSTANTS(:,115);
CONSTANTS(:,330) = CONSTANTS(:,116);
CONSTANTS(:,331) = CONSTANTS(:,117);
CONSTANTS(:,332) = CONSTANTS(:,118);
CONSTANTS(:,333) = CONSTANTS(:,119);
CONSTANTS(:,334) = CONSTANTS(:,120);
CONSTANTS(:,335) = CONSTANTS(:,121);
CONSTANTS(:,336) = CONSTANTS(:,122);
CONSTANTS(:,337) = CONSTANTS(:,123);
CONSTANTS(:,338) = CONSTANTS(:,124);
CONSTANTS(:,339) = CONSTANTS(:,128);
CONSTANTS(:,340) = CONSTANTS(:,127);
CONSTANTS(:,341) = CONSTANTS(:,129);
CONSTANTS(:,342) = CONSTANTS(:,133);
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(:,221) = piecewise({CONSTANTS(:,39) ~= 45.0000, CONSTANTS(:,6) , VOI<40.0000, 0.00100000 , VOI>=40.0000&VOI<80.0000, 0.00400000 , VOI>=80.0000&VOI<100.000, 0.0100000 , VOI>=100.000, 0.0400000 }, NaN);
ALGEBRAIC(:,222) = 1.00000+STATES(:,5)./CONSTANTS(:,40)+STATES(:,4)./CONSTANTS(:,265)+( STATES(:,5).*STATES(:,4))./( CONSTANTS(:,40).*CONSTANTS(:,41))+STATES(:,5)./CONSTANTS(:,267)+STATES(:,3)./CONSTANTS(:,42)+( STATES(:,5).*STATES(:,3))./( CONSTANTS(:,268).*CONSTANTS(:,267));
ALGEBRAIC(:,1) = piecewise({VOI<=1.00000||VOI>1.00000&CONSTANTS(:,35)==0.00000, STATES(:,1) }, CONSTANTS(:,36));
ALGEBRAIC(:,223) = 1.40400./(1.00000+power(10.0000, 5.94000 - ALGEBRAIC(:,1))+power(10.0000, ALGEBRAIC(:,1) - 7.29000));
ALGEBRAIC(:,225) = (( ALGEBRAIC(:,223).*CONSTANTS(:,264).*STATES(:,4))./( CONSTANTS(:,266).*CONSTANTS(:,265)))./ALGEBRAIC(:,222);
ALGEBRAIC(:,2) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,151))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,145));
ALGEBRAIC(:,61) = 1.00000+ power(10.0000, - ALGEBRAIC(:,1)).*power(10.0000, CONSTANTS(:,196))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,197));
ALGEBRAIC(:,12) = (CONSTANTS(:,153)+ CONSTANTS(:,152).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,152));
ALGEBRAIC(:,70) = (CONSTANTS(:,199)+ CONSTANTS(:,198).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,198));
ALGEBRAIC(:,170) = 655.700+ CONSTANTS(:,249).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1);
ALGEBRAIC(:,183) = (ALGEBRAIC(:,170)+ALGEBRAIC(:,70)) - ALGEBRAIC(:,12);
ALGEBRAIC(:,184) = ( exp( - ALGEBRAIC(:,183)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,61))./ALGEBRAIC(:,2);
ALGEBRAIC(:,224) = ( ALGEBRAIC(:,223).*CONSTANTS(:,264).*CONSTANTS(:,267).*CONSTANTS(:,268))./( CONSTANTS(:,266).*CONSTANTS(:,265).*ALGEBRAIC(:,184));
ALGEBRAIC(:,226) = (( ALGEBRAIC(:,224).*STATES(:,5))./( CONSTANTS(:,267).*CONSTANTS(:,268)))./ALGEBRAIC(:,222);
ALGEBRAIC(:,227) = ALGEBRAIC(:,221).*( STATES(:,5).*ALGEBRAIC(:,225) - STATES(:,3).*ALGEBRAIC(:,226));
ALGEBRAIC(:,228) = 1.00000 - ALGEBRAIC(:,221);
ALGEBRAIC(:,229) = (power(STATES(:,6)./CONSTANTS(:,45), CONSTANTS(:,47))./CONSTANTS(:,46))./(1.00000+power(STATES(:,6)./CONSTANTS(:,45), CONSTANTS(:,47))./CONSTANTS(:,46));
ALGEBRAIC(:,230) = 1.00000+STATES(:,5)./CONSTANTS(:,270)+STATES(:,4)./CONSTANTS(:,43)+STATES(:,5)./CONSTANTS(:,272)+STATES(:,3)./CONSTANTS(:,44)+( STATES(:,5).*STATES(:,4))./( CONSTANTS(:,270).*CONSTANTS(:,269))+( STATES(:,5).*STATES(:,3))./( CONSTANTS(:,271).*CONSTANTS(:,272));
ALGEBRAIC(:,231) = 1.75000./(1.00000+power(10.0000, 6.12000 - ALGEBRAIC(:,1))+power(10.0000, ALGEBRAIC(:,1) - 7.03000));
ALGEBRAIC(:,233) = (( ALGEBRAIC(:,231).*ALGEBRAIC(:,229).*CONSTANTS(:,264).*STATES(:,4))./( CONSTANTS(:,270).*CONSTANTS(:,269)))./ALGEBRAIC(:,230);
ALGEBRAIC(:,232) = ( ALGEBRAIC(:,231).*CONSTANTS(:,264).*CONSTANTS(:,271).*CONSTANTS(:,272))./( CONSTANTS(:,270).*CONSTANTS(:,269).*ALGEBRAIC(:,184));
ALGEBRAIC(:,234) = (( ALGEBRAIC(:,229).*ALGEBRAIC(:,232).*STATES(:,5))./( CONSTANTS(:,271).*CONSTANTS(:,272)))./ALGEBRAIC(:,230);
ALGEBRAIC(:,235) = ALGEBRAIC(:,228).*( STATES(:,5).*ALGEBRAIC(:,233) - STATES(:,3).*ALGEBRAIC(:,234));
ALGEBRAIC(:,237) = - (ALGEBRAIC(:,227)+ALGEBRAIC(:,235));
RATES(:,5) = ALGEBRAIC(:,237);
ALGEBRAIC(:,236) = ( CONSTANTS(:,273).*1.32900)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+6.64000)+power(10.0000, ALGEBRAIC(:,1) - 8.36000));
ALGEBRAIC(:,71) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,200));
ALGEBRAIC(:,76) = (CONSTANTS(:,202)+ CONSTANTS(:,201).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,201));
ALGEBRAIC(:,186) = ALGEBRAIC(:,76) - ALGEBRAIC(:,70);
ALGEBRAIC(:,187) = ( exp( - ALGEBRAIC(:,186)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,71))./ALGEBRAIC(:,61);
ALGEBRAIC(:,238) = ( ALGEBRAIC(:,236).*CONSTANTS(:,275))./( CONSTANTS(:,274).*ALGEBRAIC(:,187));
ALGEBRAIC(:,239) = (( ALGEBRAIC(:,236).*STATES(:,3))./CONSTANTS(:,274) - ( ALGEBRAIC(:,238).*STATES(:,7))./CONSTANTS(:,275))./(1.00000+STATES(:,3)./CONSTANTS(:,274)+STATES(:,7)./CONSTANTS(:,275));
ALGEBRAIC(:,240) = CONSTANTS(:,276)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+6.94000)+power(10.0000, ALGEBRAIC(:,1) - 9.35000));
ALGEBRAIC(:,77) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,203));
ALGEBRAIC(:,82) = (CONSTANTS(:,205)+ CONSTANTS(:,204).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,204));
ALGEBRAIC(:,189) = ALGEBRAIC(:,82) - ALGEBRAIC(:,76);
ALGEBRAIC(:,190) = ( exp( - ALGEBRAIC(:,189)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,77))./ALGEBRAIC(:,71);
ALGEBRAIC(:,241) = (( ALGEBRAIC(:,240).*CONSTANTS(:,277))./CONSTANTS(:,278)).*ALGEBRAIC(:,190);
ALGEBRAIC(:,242) = (( ALGEBRAIC(:,241).*STATES(:,7))./CONSTANTS(:,277) - ( ALGEBRAIC(:,240).*STATES(:,8))./CONSTANTS(:,278))./(1.00000+STATES(:,8)./CONSTANTS(:,278)+STATES(:,7)./CONSTANTS(:,277));
ALGEBRAIC(:,244) = ALGEBRAIC(:,239) - ALGEBRAIC(:,242);
RATES(:,7) = ALGEBRAIC(:,244);
ALGEBRAIC(:,243) = CONSTANTS(:,279)./(1.00000+power(ALGEBRAIC(:,1)./6.80000, - 30.0000));
ALGEBRAIC(:,13) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,161))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,162));
ALGEBRAIC(:,24) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,171))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,165));
ALGEBRAIC(:,83) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206))+power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,208));
ALGEBRAIC(:,23) = (CONSTANTS(:,164)+ CONSTANTS(:,163).*CONSTANTS(:,1).*CONSTANTS(:,131).*log(10.0000).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(16.0000 - CONSTANTS(:,163));
ALGEBRAIC(:,35) = (CONSTANTS(:,173)+ CONSTANTS(:,172).*CONSTANTS(:,1).*CONSTANTS(:,131).*log(10.0000).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,172));
ALGEBRAIC(:,93) = (CONSTANTS(:,210)+ CONSTANTS(:,209).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(16.0000 - CONSTANTS(:,209));
ALGEBRAIC(:,174) = CONSTANTS(:,257)+ CONSTANTS(:,256).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1);
ALGEBRAIC(:,192) = ((ALGEBRAIC(:,93)+ALGEBRAIC(:,35)+ALGEBRAIC(:,174)) - ALGEBRAIC(:,82)) - ALGEBRAIC(:,23);
ALGEBRAIC(:,193) = ( exp( - ALGEBRAIC(:,192)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,83).*ALGEBRAIC(:,24))./( ALGEBRAIC(:,77).*ALGEBRAIC(:,13).*power(10.0000, - ALGEBRAIC(:,1)));
ALGEBRAIC(:,245) = ( ALGEBRAIC(:,243).*CONSTANTS(:,283).*CONSTANTS(:,285))./( CONSTANTS(:,280).*CONSTANTS(:,282).*ALGEBRAIC(:,193));
ALGEBRAIC(:,246) = CONSTANTS(:,54).*power(( ((1.00000+STATES(:,11)./CONSTANTS(:,50))./(1.00000+( CONSTANTS(:,52).*STATES(:,11))./CONSTANTS(:,50))).*(1.00000+( CONSTANTS(:,53).*STATES(:,6))./CONSTANTS(:,51)))./(1.00000+STATES(:,6)./CONSTANTS(:,51)), 4.00000);
ALGEBRAIC(:,247) = (1.00000+STATES(:,8)./CONSTANTS(:,280)).*(1.00000+STATES(:,11)./CONSTANTS(:,282))+STATES(:,9)./CONSTANTS(:,283)+ (STATES(:,10)./CONSTANTS(:,285)).*(1.00000+STATES(:,9)./CONSTANTS(:,283));
ALGEBRAIC(:,248) = (1.00000+STATES(:,8)./CONSTANTS(:,281)).*(1.00000+STATES(:,11)./CONSTANTS(:,284))+STATES(:,9)./CONSTANTS(:,48)+ (STATES(:,10)./CONSTANTS(:,49)).*(1.00000+STATES(:,9)./CONSTANTS(:,48));
ALGEBRAIC(:,249) = ( ((( ALGEBRAIC(:,243).*STATES(:,8).*STATES(:,11))./( CONSTANTS(:,280).*CONSTANTS(:,282)))./ALGEBRAIC(:,247) - (( ALGEBRAIC(:,245).*STATES(:,10).*STATES(:,9))./( CONSTANTS(:,285).*CONSTANTS(:,283)))./ALGEBRAIC(:,247)).*(1.00000+ CONSTANTS(:,286).*ALGEBRAIC(:,246).*power(ALGEBRAIC(:,248)./ALGEBRAIC(:,247), 3.00000)))./(1.00000+ ALGEBRAIC(:,246).*power(ALGEBRAIC(:,248)./ALGEBRAIC(:,247), 4.00000));
ALGEBRAIC(:,251) = ALGEBRAIC(:,242) - ALGEBRAIC(:,249);
RATES(:,8) = ALGEBRAIC(:,251);
ALGEBRAIC(:,250) = ( CONSTANTS(:,287).*1.01300)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+5.32000)+power(10.0000, ALGEBRAIC(:,1) - 9.15000));
ALGEBRAIC(:,94) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,211));
ALGEBRAIC(:,110) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,221));
ALGEBRAIC(:,99) = (CONSTANTS(:,213)+ CONSTANTS(:,212).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,212));
ALGEBRAIC(:,119) = (CONSTANTS(:,223)+ CONSTANTS(:,222).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,222));
ALGEBRAIC(:,195) = (ALGEBRAIC(:,119)+ALGEBRAIC(:,99)) - ALGEBRAIC(:,93);
ALGEBRAIC(:,196) = ( 1.00000.*exp( - ALGEBRAIC(:,195)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,94).*ALGEBRAIC(:,83))./ALGEBRAIC(:,110);
ALGEBRAIC(:,252) = ( ALGEBRAIC(:,250).*CONSTANTS(:,290).*CONSTANTS(:,289))./( CONSTANTS(:,288).*ALGEBRAIC(:,196));
ALGEBRAIC(:,253) = (( ALGEBRAIC(:,250).*STATES(:,9))./CONSTANTS(:,288) - ( ALGEBRAIC(:,252).*STATES(:,13).*STATES(:,12))./( CONSTANTS(:,290).*CONSTANTS(:,289)))./(1.00000+STATES(:,9)./CONSTANTS(:,288)+STATES(:,13)./CONSTANTS(:,290)+STATES(:,12)./CONSTANTS(:,289));
ALGEBRAIC(:,198) = ALGEBRAIC(:,119) - ALGEBRAIC(:,99);
ALGEBRAIC(:,199) = ( exp( - ALGEBRAIC(:,198)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,110))./ALGEBRAIC(:,94);
ALGEBRAIC(:,254) = ( CONSTANTS(:,294).*CONSTANTS(:,293))./( CONSTANTS(:,292).*ALGEBRAIC(:,199));
ALGEBRAIC(:,255) = (( CONSTANTS(:,294).*STATES(:,13))./CONSTANTS(:,292) - ( ALGEBRAIC(:,254).*STATES(:,12))./CONSTANTS(:,293))./(1.00000+STATES(:,13)./CONSTANTS(:,292)+STATES(:,12)./CONSTANTS(:,293));
ALGEBRAIC(:,259) = 1.00000+STATES(:,4)./CONSTANTS(:,304)+STATES(:,13)./CONSTANTS(:,302)+STATES(:,15)./CONSTANTS(:,303)+( STATES(:,13).*STATES(:,15))./( CONSTANTS(:,302).*CONSTANTS(:,303))+( STATES(:,13).*STATES(:,15).*STATES(:,4))./( CONSTANTS(:,302).*CONSTANTS(:,303).*CONSTANTS(:,304))+STATES(:,17)./CONSTANTS(:,305)+STATES(:,16)./CONSTANTS(:,306)+( STATES(:,17).*STATES(:,16))./( CONSTANTS(:,306).*CONSTANTS(:,305));
ALGEBRAIC(:,260) = CONSTANTS(:,301).*0.000700000.*exp( ALGEBRAIC(:,1).*0.897900);
ALGEBRAIC(:,120) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,224));
ALGEBRAIC(:,125) = (CONSTANTS(:,226)+ CONSTANTS(:,225).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(16.0000 - CONSTANTS(:,225));
ALGEBRAIC(:,171) = (CONSTANTS(:,251)+ CONSTANTS(:,250).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(1.00000 - CONSTANTS(:,250));
ALGEBRAIC(:,172) = (CONSTANTS(:,253)+ CONSTANTS(:,252).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,252));
ALGEBRAIC(:,201) = (((ALGEBRAIC(:,125)+ALGEBRAIC(:,172)+ALGEBRAIC(:,174)) - ALGEBRAIC(:,12)) - ALGEBRAIC(:,99)) - ALGEBRAIC(:,171);
ALGEBRAIC(:,202) = ( exp( - ALGEBRAIC(:,201)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,120))./( ALGEBRAIC(:,2).*ALGEBRAIC(:,94).*power(10.0000, - ALGEBRAIC(:,1)).*1.00000);
ALGEBRAIC(:,261) = ( ALGEBRAIC(:,260).*CONSTANTS(:,305).*CONSTANTS(:,306))./( CONSTANTS(:,302).*CONSTANTS(:,304).*CONSTANTS(:,303).*ALGEBRAIC(:,202));
ALGEBRAIC(:,262) = (( ALGEBRAIC(:,260).*STATES(:,13).*STATES(:,15).*STATES(:,4))./( CONSTANTS(:,303).*CONSTANTS(:,302).*CONSTANTS(:,304)) - ( ALGEBRAIC(:,261).*STATES(:,17).*STATES(:,16))./( CONSTANTS(:,305).*CONSTANTS(:,306)))./ALGEBRAIC(:,259);
ALGEBRAIC(:,264) = (ALGEBRAIC(:,253) - ALGEBRAIC(:,255)) - ALGEBRAIC(:,262);
RATES(:,13) = ALGEBRAIC(:,264);
ALGEBRAIC(:,126) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+6.21000);
ALGEBRAIC(:,131) = (CONSTANTS(:,229)+ CONSTANTS(:,228).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,228));
ALGEBRAIC(:,207) = ((ALGEBRAIC(:,23)+ALGEBRAIC(:,131)) - ALGEBRAIC(:,125)) - ALGEBRAIC(:,35);
ALGEBRAIC(:,208) = ( exp( - ALGEBRAIC(:,207)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,13).*ALGEBRAIC(:,126))./( ALGEBRAIC(:,120).*ALGEBRAIC(:,24));
ALGEBRAIC(:,263) = (( CONSTANTS(:,312).*CONSTANTS(:,308).*CONSTANTS(:,309))./( CONSTANTS(:,310).*CONSTANTS(:,311))).*ALGEBRAIC(:,208);
ALGEBRAIC(:,265) = 1.00000+STATES(:,10)./CONSTANTS(:,309)+STATES(:,17)./CONSTANTS(:,308)+( STATES(:,17).*STATES(:,10))./( CONSTANTS(:,308).*CONSTANTS(:,309))+STATES(:,18)./CONSTANTS(:,310)+STATES(:,11)./CONSTANTS(:,311)+( STATES(:,18).*STATES(:,11))./( CONSTANTS(:,310).*CONSTANTS(:,311));
ALGEBRAIC(:,266) = (( ALGEBRAIC(:,263).*STATES(:,17).*STATES(:,10))./( CONSTANTS(:,309).*CONSTANTS(:,308)) - ( CONSTANTS(:,312).*STATES(:,11).*STATES(:,18))./( CONSTANTS(:,311).*CONSTANTS(:,310)))./ALGEBRAIC(:,265);
ALGEBRAIC(:,268) = ALGEBRAIC(:,262) - ALGEBRAIC(:,266);
RATES(:,17) = ALGEBRAIC(:,268);
ALGEBRAIC(:,267) = ( CONSTANTS(:,313).*0.989000)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+5.62000)+power(10.0000, ALGEBRAIC(:,1) - 8.74000));
ALGEBRAIC(:,132) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,231))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,232));
ALGEBRAIC(:,142) = (CONSTANTS(:,234)+ CONSTANTS(:,233).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,233));
ALGEBRAIC(:,210) = ALGEBRAIC(:,142) - ALGEBRAIC(:,131);
ALGEBRAIC(:,211) = ( exp( - ALGEBRAIC(:,210)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,132))./ALGEBRAIC(:,126);
ALGEBRAIC(:,269) = ( ALGEBRAIC(:,267).*CONSTANTS(:,315))./( CONSTANTS(:,314).*ALGEBRAIC(:,211));
ALGEBRAIC(:,270) = (( ALGEBRAIC(:,267).*STATES(:,18))./CONSTANTS(:,314) - ( ALGEBRAIC(:,269).*STATES(:,19))./CONSTANTS(:,315))./(1.00000+STATES(:,18)./CONSTANTS(:,314)+STATES(:,19)./CONSTANTS(:,315));
ALGEBRAIC(:,272) = ALGEBRAIC(:,266) - ALGEBRAIC(:,270);
RATES(:,18) = ALGEBRAIC(:,272);
ALGEBRAIC(:,256) = (1.00000+STATES(:,14)./CONSTANTS(:,296)+STATES(:,16)./CONSTANTS(:,299)).*(1.00000+STATES(:,12)./CONSTANTS(:,298)+STATES(:,15)./CONSTANTS(:,299));
ALGEBRAIC(:,100) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,214));
ALGEBRAIC(:,109) = (CONSTANTS(:,219)+ CONSTANTS(:,218).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,218));
ALGEBRAIC(:,204) = ((ALGEBRAIC(:,174)+ALGEBRAIC(:,172)+ALGEBRAIC(:,119)) - ALGEBRAIC(:,171)) - ALGEBRAIC(:,109);
ALGEBRAIC(:,205) = ( exp( - ALGEBRAIC(:,204)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,110))./( ALGEBRAIC(:,100).*power(10.0000, - ALGEBRAIC(:,1)));
ALGEBRAIC(:,257) = ( CONSTANTS(:,300).*CONSTANTS(:,296).*CONSTANTS(:,297).*ALGEBRAIC(:,205))./( CONSTANTS(:,298).*CONSTANTS(:,299));
ALGEBRAIC(:,258) = (( ALGEBRAIC(:,257).*STATES(:,14).*STATES(:,15))./( CONSTANTS(:,296).*CONSTANTS(:,297)) - ( CONSTANTS(:,300).*STATES(:,12).*STATES(:,16))./( CONSTANTS(:,298).*CONSTANTS(:,299)))./ALGEBRAIC(:,256);
ALGEBRAIC(:,277) = CONSTANTS(:,325).*( - 0.113400.*ALGEBRAIC(:,1)+1.60690);
ALGEBRAIC(:,154) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,240));
ALGEBRAIC(:,160) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,243));
ALGEBRAIC(:,159) = (CONSTANTS(:,242)+ CONSTANTS(:,241).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(1.00000 - CONSTANTS(:,241));
ALGEBRAIC(:,169) = (CONSTANTS(:,248)+ CONSTANTS(:,247).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(1.00000 - CONSTANTS(:,247));
ALGEBRAIC(:,219) = (((ALGEBRAIC(:,169)+ALGEBRAIC(:,171)) - ALGEBRAIC(:,159)) - ALGEBRAIC(:,172)) - ALGEBRAIC(:,174);
ALGEBRAIC(:,220) = ( exp( - ALGEBRAIC(:,219)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,160).*power(10.0000, - ALGEBRAIC(:,1)))./ALGEBRAIC(:,154);
ALGEBRAIC(:,278) = ( ALGEBRAIC(:,277).*CONSTANTS(:,328).*CONSTANTS(:,329))./( CONSTANTS(:,326).*CONSTANTS(:,327).*ALGEBRAIC(:,220));
ALGEBRAIC(:,279) = (( ALGEBRAIC(:,277).*STATES(:,21).*STATES(:,16))./( CONSTANTS(:,326).*CONSTANTS(:,327)) - ( ALGEBRAIC(:,278).*STATES(:,22).*STATES(:,15))./( CONSTANTS(:,328).*CONSTANTS(:,329)))./(1.00000+STATES(:,21)./CONSTANTS(:,326)+STATES(:,16)./CONSTANTS(:,327)+( STATES(:,21).*STATES(:,16))./( CONSTANTS(:,326).*CONSTANTS(:,327))+STATES(:,22)./CONSTANTS(:,328)+STATES(:,15)./CONSTANTS(:,329)+( STATES(:,22).*STATES(:,15))./( CONSTANTS(:,328).*CONSTANTS(:,329)));
ALGEBRAIC(:,282) = ( - ALGEBRAIC(:,262) - ALGEBRAIC(:,258))+ALGEBRAIC(:,279);
RATES(:,15) = ALGEBRAIC(:,282);
ALGEBRAIC(:,283) = (ALGEBRAIC(:,262)+ALGEBRAIC(:,258)) - ALGEBRAIC(:,279);
RATES(:,16) = ALGEBRAIC(:,283);
ALGEBRAIC(:,274) = ( CONSTANTS(:,320).*1.05000)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+5.58000)+power(10.0000, ALGEBRAIC(:,1) - 8.79000));
ALGEBRAIC(:,143) = 1.00000+power(10.0000, CONSTANTS(:,235) - ALGEBRAIC(:,1))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,236))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,237));
ALGEBRAIC(:,153) = (CONSTANTS(:,239)+ CONSTANTS(:,238).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,238));
ALGEBRAIC(:,216) = (((ALGEBRAIC(:,159)+ALGEBRAIC(:,23)) - ALGEBRAIC(:,174)) - ALGEBRAIC(:,153)) - ALGEBRAIC(:,35);
ALGEBRAIC(:,217) = ( exp( - ALGEBRAIC(:,216)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,154).*ALGEBRAIC(:,13).*power(10.0000, - ALGEBRAIC(:,1)))./( ALGEBRAIC(:,143).*ALGEBRAIC(:,24));
ALGEBRAIC(:,275) = ( ALGEBRAIC(:,274).*CONSTANTS(:,323).*CONSTANTS(:,324))./( CONSTANTS(:,321).*CONSTANTS(:,322).*ALGEBRAIC(:,217));
ALGEBRAIC(:,276) = (( ALGEBRAIC(:,274).*STATES(:,20).*STATES(:,10))./( CONSTANTS(:,321).*CONSTANTS(:,322)) - ( ALGEBRAIC(:,275).*STATES(:,21).*STATES(:,11))./( CONSTANTS(:,323).*CONSTANTS(:,324)))./(1.00000+STATES(:,20)./CONSTANTS(:,321)+STATES(:,10)./CONSTANTS(:,322)+( STATES(:,20).*STATES(:,10))./( CONSTANTS(:,321).*CONSTANTS(:,322))+STATES(:,11)./CONSTANTS(:,324)+STATES(:,21)./CONSTANTS(:,323)+( STATES(:,21).*STATES(:,11))./( CONSTANTS(:,323).*CONSTANTS(:,324)));
ALGEBRAIC(:,284) = ALGEBRAIC(:,276) - ALGEBRAIC(:,279);
RATES(:,21) = ALGEBRAIC(:,284);
ALGEBRAIC(:,46) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,188))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,182));
ALGEBRAIC(:,56) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,190));
ALGEBRAIC(:,177) = ( exp( - CONSTANTS(:,263)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*power(10.0000, - ALGEBRAIC(:,1)).*ALGEBRAIC(:,13).*ALGEBRAIC(:,56))./( ALGEBRAIC(:,46).*ALGEBRAIC(:,24));
ALGEBRAIC(:,286) = ( (CONSTANTS(:,330)./ALGEBRAIC(:,177)).*CONSTANTS(:,332).*CONSTANTS(:,334))./( CONSTANTS(:,333).*CONSTANTS(:,331));
ALGEBRAIC(:,287) = (( ALGEBRAIC(:,286).*STATES(:,11).*STATES(:,23))./( CONSTANTS(:,332).*CONSTANTS(:,334)) - ( CONSTANTS(:,330).*STATES(:,10).*STATES(:,24))./( CONSTANTS(:,333).*CONSTANTS(:,331)))./(1.00000+STATES(:,10)./CONSTANTS(:,333)+STATES(:,24)./CONSTANTS(:,56)+( STATES(:,24).*STATES(:,10))./( CONSTANTS(:,333).*CONSTANTS(:,331))+STATES(:,11)./CONSTANTS(:,332)+( STATES(:,23).*STATES(:,11))./( CONSTANTS(:,334).*CONSTANTS(:,332)));
ALGEBRAIC(:,290) = - ALGEBRAIC(:,287);
RATES(:,23) = ALGEBRAIC(:,290);
ALGEBRAIC(:,14) = 1.00000./ALGEBRAIC(:,13);
ALGEBRAIC(:,15) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160)).*ALGEBRAIC(:,14);
ALGEBRAIC(:,16) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,161)).*ALGEBRAIC(:,14);
ALGEBRAIC(:,17) = (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,162)).*ALGEBRAIC(:,14);
ALGEBRAIC(:,18) = 0.00000.*ALGEBRAIC(:,14)+ 1.00000.*ALGEBRAIC(:,15)+ 0.00000.*ALGEBRAIC(:,16)+ 0.00000.*ALGEBRAIC(:,17);
ALGEBRAIC(:,25) = 1.00000./ALGEBRAIC(:,24);
ALGEBRAIC(:,26) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170)).*ALGEBRAIC(:,25);
ALGEBRAIC(:,27) = (( ALGEBRAIC(:,25).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,171));
ALGEBRAIC(:,28) = (( ALGEBRAIC(:,25).*CONSTANTS(:,134))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,165));
ALGEBRAIC(:,29) = 0.00000.*ALGEBRAIC(:,25)+ 1.00000.*ALGEBRAIC(:,26)+ 0.00000.*ALGEBRAIC(:,27)+ 0.00000.*ALGEBRAIC(:,28);
ALGEBRAIC(:,47) = 1.00000./ALGEBRAIC(:,46);
ALGEBRAIC(:,49) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187)).*ALGEBRAIC(:,47);
ALGEBRAIC(:,51) = ALGEBRAIC(:,49);
ALGEBRAIC(:,57) = 1.00000./ALGEBRAIC(:,56);
ALGEBRAIC(:,58) = ALGEBRAIC(:,57).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,190));
ALGEBRAIC(:,59) = 0.00000.*ALGEBRAIC(:,57)+ 1.00000.*ALGEBRAIC(:,58);
ALGEBRAIC(:,175) = (((ALGEBRAIC(:,18)+ALGEBRAIC(:,59)) - ALGEBRAIC(:,51)) - ALGEBRAIC(:,29))+(((CONSTANTS(:,163)+CONSTANTS(:,191)) - CONSTANTS(:,189)) - CONSTANTS(:,172));
ALGEBRAIC(:,36) = 1.00000+power(10.0000, CONSTANTS(:,178) - ALGEBRAIC(:,1))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,179));
ALGEBRAIC(:,37) = 1.00000./ALGEBRAIC(:,36);
ALGEBRAIC(:,38) = ALGEBRAIC(:,37).*power(10.0000, CONSTANTS(:,178) - ALGEBRAIC(:,1));
ALGEBRAIC(:,39) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,179)).*ALGEBRAIC(:,37);
ALGEBRAIC(:,40) = 0.00000.*ALGEBRAIC(:,37)+ALGEBRAIC(:,38)+ 0.00000.*ALGEBRAIC(:,39);
ALGEBRAIC(:,179) = (( 2.00000.*ALGEBRAIC(:,29) - ALGEBRAIC(:,18)) - ALGEBRAIC(:,40))+(( 2.00000.*CONSTANTS(:,172) - CONSTANTS(:,163)) - CONSTANTS(:,180));
ALGEBRAIC(:,3) = 1.00000./ALGEBRAIC(:,2);
ALGEBRAIC(:,4) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150)).*ALGEBRAIC(:,3);
ALGEBRAIC(:,7) = 1.00000.*ALGEBRAIC(:,4);
ALGEBRAIC(:,62) = 1.00000./ALGEBRAIC(:,61);
ALGEBRAIC(:,63) = ALGEBRAIC(:,62).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,196));
ALGEBRAIC(:,65) = ALGEBRAIC(:,63);
ALGEBRAIC(:,182) = (ALGEBRAIC(:,65) - ALGEBRAIC(:,7))+(1.00000 - CONSTANTS(:,152));
ALGEBRAIC(:,72) = 1.00000./ALGEBRAIC(:,71);
ALGEBRAIC(:,73) = ALGEBRAIC(:,72).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,200));
ALGEBRAIC(:,74) = ALGEBRAIC(:,73);
ALGEBRAIC(:,185) = (ALGEBRAIC(:,74) - ALGEBRAIC(:,65))+(CONSTANTS(:,201) - CONSTANTS(:,198));
ALGEBRAIC(:,78) = 1.00000./ALGEBRAIC(:,77);
ALGEBRAIC(:,79) = ALGEBRAIC(:,78).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,203));
ALGEBRAIC(:,80) = ALGEBRAIC(:,79);
ALGEBRAIC(:,188) = (ALGEBRAIC(:,80) - ALGEBRAIC(:,74))+(CONSTANTS(:,204) - CONSTANTS(:,201));
ALGEBRAIC(:,84) = 1.00000./ALGEBRAIC(:,83);
ALGEBRAIC(:,85) = ALGEBRAIC(:,84).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206));
ALGEBRAIC(:,86) = ALGEBRAIC(:,84).*power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207));
ALGEBRAIC(:,88) = ALGEBRAIC(:,85)+ 2.00000.*ALGEBRAIC(:,86);
ALGEBRAIC(:,191) = (((((ALGEBRAIC(:,29)+ALGEBRAIC(:,88)) - ALGEBRAIC(:,80)) - ALGEBRAIC(:,18))+CONSTANTS(:,172)+CONSTANTS(:,209)) - CONSTANTS(:,204)) - CONSTANTS(:,163);
ALGEBRAIC(:,95) = 1.00000./ALGEBRAIC(:,94);
ALGEBRAIC(:,96) = ALGEBRAIC(:,95).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,211));
ALGEBRAIC(:,97) = ALGEBRAIC(:,96);
ALGEBRAIC(:,111) = 1.00000./ALGEBRAIC(:,110);
ALGEBRAIC(:,112) = ALGEBRAIC(:,111).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220));
ALGEBRAIC(:,114) = ALGEBRAIC(:,112);
ALGEBRAIC(:,194) = ((ALGEBRAIC(:,114)+ALGEBRAIC(:,97)) - ALGEBRAIC(:,88))+((CONSTANTS(:,222)+CONSTANTS(:,212)) - CONSTANTS(:,209));
ALGEBRAIC(:,197) = (ALGEBRAIC(:,114) - ALGEBRAIC(:,97))+(CONSTANTS(:,222) - CONSTANTS(:,212));
ALGEBRAIC(:,121) = 1.00000./ALGEBRAIC(:,120);
ALGEBRAIC(:,122) = ALGEBRAIC(:,121).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,224));
ALGEBRAIC(:,123) = ALGEBRAIC(:,122);
ALGEBRAIC(:,200) = ((ALGEBRAIC(:,123) - ALGEBRAIC(:,97)) - ALGEBRAIC(:,7))+((((CONSTANTS(:,225)+CONSTANTS(:,252)) - CONSTANTS(:,212)) - CONSTANTS(:,152)) - CONSTANTS(:,250));
ALGEBRAIC(:,101) = 1.00000./ALGEBRAIC(:,100);
ALGEBRAIC(:,102) = ALGEBRAIC(:,101).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217));
ALGEBRAIC(:,104) = ALGEBRAIC(:,102);
ALGEBRAIC(:,203) = (ALGEBRAIC(:,114) - ALGEBRAIC(:,104))+(((CONSTANTS(:,222)+CONSTANTS(:,252)) - CONSTANTS(:,250)) - CONSTANTS(:,218));
ALGEBRAIC(:,127) = 1.00000./ALGEBRAIC(:,126);
ALGEBRAIC(:,128) = ALGEBRAIC(:,127).*power(10.0000, - ALGEBRAIC(:,1)+6.21000);
ALGEBRAIC(:,129) = ALGEBRAIC(:,128);
ALGEBRAIC(:,206) = (((ALGEBRAIC(:,129)+ALGEBRAIC(:,18)) - ALGEBRAIC(:,123)) - ALGEBRAIC(:,29))+(((CONSTANTS(:,228)+CONSTANTS(:,163)) - CONSTANTS(:,225)) - CONSTANTS(:,172));
ALGEBRAIC(:,133) = 1.00000./ALGEBRAIC(:,132);
ALGEBRAIC(:,134) = ALGEBRAIC(:,133).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230));
ALGEBRAIC(:,137) = ALGEBRAIC(:,134);
ALGEBRAIC(:,209) = (ALGEBRAIC(:,137) - ALGEBRAIC(:,129))+(CONSTANTS(:,233) - CONSTANTS(:,228));
ALGEBRAIC(:,144) = 1.00000./ALGEBRAIC(:,143);
ALGEBRAIC(:,145) = ALGEBRAIC(:,144).*power(10.0000, CONSTANTS(:,235) - ALGEBRAIC(:,1));
ALGEBRAIC(:,148) = ALGEBRAIC(:,145);
ALGEBRAIC(:,212) = (ALGEBRAIC(:,148) - ALGEBRAIC(:,137))+((CONSTANTS(:,254)+CONSTANTS(:,238)) - CONSTANTS(:,233));
ALGEBRAIC(:,155) = 1.00000./ALGEBRAIC(:,154);
ALGEBRAIC(:,156) = ALGEBRAIC(:,155).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,240));
ALGEBRAIC(:,157) = ALGEBRAIC(:,156);
ALGEBRAIC(:,215) = (((ALGEBRAIC(:,157)+ALGEBRAIC(:,18)) - ALGEBRAIC(:,148)) - ALGEBRAIC(:,29))+(((CONSTANTS(:,241)+CONSTANTS(:,163)) - CONSTANTS(:,238)) - CONSTANTS(:,172));
ALGEBRAIC(:,161) = 1.00000./ALGEBRAIC(:,160);
ALGEBRAIC(:,162) = ALGEBRAIC(:,161).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246));
ALGEBRAIC(:,164) = ALGEBRAIC(:,162);
ALGEBRAIC(:,218) = (ALGEBRAIC(:,164) - ALGEBRAIC(:,157))+(((CONSTANTS(:,247)+CONSTANTS(:,250)) - CONSTANTS(:,241)) - CONSTANTS(:,252));
ALGEBRAIC(:,173) = (CONSTANTS(:,255)+ CONSTANTS(:,254).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(0.00000 - CONSTANTS(:,254));
ALGEBRAIC(:,213) = (ALGEBRAIC(:,173)+ALGEBRAIC(:,153)) - ALGEBRAIC(:,142);
ALGEBRAIC(:,214) = ( exp( - ALGEBRAIC(:,213)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,143))./ALGEBRAIC(:,132);
ALGEBRAIC(:,271) = ( CONSTANTS(:,319).*CONSTANTS(:,318))./( CONSTANTS(:,317).*ALGEBRAIC(:,214));
ALGEBRAIC(:,273) = (( CONSTANTS(:,319).*STATES(:,19))./CONSTANTS(:,317) - ( ALGEBRAIC(:,271).*STATES(:,20))./CONSTANTS(:,318))./(1.00000+STATES(:,20)./CONSTANTS(:,318)+STATES(:,19)./CONSTANTS(:,317));
ALGEBRAIC(:,45) = (CONSTANTS(:,181)+ CONSTANTS(:,180).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,180));
ALGEBRAIC(:,180) = ( 2.00000.*ALGEBRAIC(:,35) - ALGEBRAIC(:,23)) - ALGEBRAIC(:,45);
ALGEBRAIC(:,181) = ( exp( - ALGEBRAIC(:,180)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*power(ALGEBRAIC(:,24), 2.00000))./( ALGEBRAIC(:,13).*ALGEBRAIC(:,36));
ALGEBRAIC(:,288) = ( CONSTANTS(:,335).*power(CONSTANTS(:,338), 2.00000))./( CONSTANTS(:,336).*CONSTANTS(:,337).*ALGEBRAIC(:,181));
ALGEBRAIC(:,291) = (( CONSTANTS(:,335).*STATES(:,11).*STATES(:,6))./( CONSTANTS(:,337).*CONSTANTS(:,336)) - ALGEBRAIC(:,288).*power(STATES(:,10)./CONSTANTS(:,338), 2.00000))./(1.00000+STATES(:,11)./CONSTANTS(:,337)+STATES(:,6)./CONSTANTS(:,336)+( STATES(:,11).*STATES(:,6))./( CONSTANTS(:,337).*CONSTANTS(:,336))+( 2.00000.*STATES(:,10))./CONSTANTS(:,338)+power(STATES(:,10), 2.00000)./power(CONSTANTS(:,338), 2.00000));
ALGEBRAIC(:,294) = ALGEBRAIC(:,175).* - ALGEBRAIC(:,287)+ ALGEBRAIC(:,179).*ALGEBRAIC(:,291)+ ALGEBRAIC(:,185).*ALGEBRAIC(:,239)+ ALGEBRAIC(:,182).*(ALGEBRAIC(:,227)+ALGEBRAIC(:,235))+ ALGEBRAIC(:,188).*ALGEBRAIC(:,242)+ ALGEBRAIC(:,191).*ALGEBRAIC(:,249)+ ALGEBRAIC(:,194).*ALGEBRAIC(:,253)+ ALGEBRAIC(:,197).*ALGEBRAIC(:,255)+ ALGEBRAIC(:,200).*ALGEBRAIC(:,262)+ ALGEBRAIC(:,206).*ALGEBRAIC(:,266)+ ALGEBRAIC(:,209).*ALGEBRAIC(:,270)+ ALGEBRAIC(:,212).*ALGEBRAIC(:,273)+ ALGEBRAIC(:,215).*ALGEBRAIC(:,276)+ ALGEBRAIC(:,218).*ALGEBRAIC(:,279)+ ALGEBRAIC(:,203).*ALGEBRAIC(:,258);
RATES(:,25) = - ALGEBRAIC(:,294);
ALGEBRAIC(:,300) = ALGEBRAIC(:,287);
RATES(:,24) = ALGEBRAIC(:,300);
ALGEBRAIC(:,280) = piecewise({CONSTANTS(:,39) ~= 5.00000&CONSTANTS(:,39) ~= 45.0000, CONSTANTS(:,23) , CONSTANTS(:,39) ~= 5.00000&CONSTANTS(:,39)==45.0000&VOI>=110.000, CONSTANTS(:,23) , CONSTANTS(:,39)==5.00000&VOI>30.0000, CONSTANTS(:,23) }, 0.00000);
ALGEBRAIC(:,285) = ( ALGEBRAIC(:,280).*STATES(:,11))./(CONSTANTS(:,55)+STATES(:,11));
ALGEBRAIC(:,301) = ((( - ALGEBRAIC(:,287) - ALGEBRAIC(:,291)) - ALGEBRAIC(:,249))+ALGEBRAIC(:,266)+ALGEBRAIC(:,276)) - ALGEBRAIC(:,285);
RATES(:,11) = ALGEBRAIC(:,301);
ALGEBRAIC(:,302) = (((ALGEBRAIC(:,287)+ 2.00000.*ALGEBRAIC(:,291)+ALGEBRAIC(:,249)) - ALGEBRAIC(:,266)) - ALGEBRAIC(:,276))+ALGEBRAIC(:,285);
RATES(:,10) = ALGEBRAIC(:,302);
ALGEBRAIC(:,303) = - ALGEBRAIC(:,291);
RATES(:,6) = ALGEBRAIC(:,303);
ALGEBRAIC(:,304) = ( - (ALGEBRAIC(:,227)+ALGEBRAIC(:,235)) - ALGEBRAIC(:,262))+ALGEBRAIC(:,285);
RATES(:,4) = ALGEBRAIC(:,304);
ALGEBRAIC(:,305) = (ALGEBRAIC(:,227)+ALGEBRAIC(:,235)) - ALGEBRAIC(:,239);
RATES(:,3) = ALGEBRAIC(:,305);
ALGEBRAIC(:,306) = ALGEBRAIC(:,249) - ALGEBRAIC(:,253);
RATES(:,9) = ALGEBRAIC(:,306);
ALGEBRAIC(:,307) = ALGEBRAIC(:,253)+ALGEBRAIC(:,255)+ALGEBRAIC(:,258);
RATES(:,12) = ALGEBRAIC(:,307);
ALGEBRAIC(:,308) = - ALGEBRAIC(:,258);
RATES(:,14) = ALGEBRAIC(:,308);
ALGEBRAIC(:,309) = ALGEBRAIC(:,270) - ALGEBRAIC(:,273);
RATES(:,19) = ALGEBRAIC(:,309);
ALGEBRAIC(:,310) = ALGEBRAIC(:,273) - ALGEBRAIC(:,276);
RATES(:,20) = ALGEBRAIC(:,310);
ALGEBRAIC(:,311) = ALGEBRAIC(:,279);
RATES(:,22) = ALGEBRAIC(:,311);
ALGEBRAIC(:,292) = ( log(10.0000).*CONSTANTS(:,339).*power(10.0000, - ALGEBRAIC(:,1) - 6.87000))./power(power(10.0000, - ALGEBRAIC(:,1))+power(10.0000, - 6.87000), 2.00000)+( log(10.0000).*CONSTANTS(:,340).*power(10.0000, - ALGEBRAIC(:,1) - 8.30000))./power(power(10.0000, - ALGEBRAIC(:,1))+power(10.0000, - 8.30000), 2.00000)+( log(10.0000).*CONSTANTS(:,341).*power(10.0000, - ALGEBRAIC(:,1) - 4.80000))./power(power(10.0000, - ALGEBRAIC(:,1))+power(10.0000, - 4.80000), 2.00000);
ALGEBRAIC(:,8) = ( power(10.0000, CONSTANTS(:,150)).*(ALGEBRAIC(:,2) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,2), 2.00000));
ALGEBRAIC(:,19) = ( power(10.0000, CONSTANTS(:,160)).*(ALGEBRAIC(:,13) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,13), 2.00000));
ALGEBRAIC(:,30) = ( power(10.0000, CONSTANTS(:,170)).*(ALGEBRAIC(:,24) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,24), 2.00000));
ALGEBRAIC(:,41) = ( power(10.0000, CONSTANTS(:,178)).*(ALGEBRAIC(:,36) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,178))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,36), 2.00000));
ALGEBRAIC(:,52) = ( power(10.0000, CONSTANTS(:,187)).*(ALGEBRAIC(:,46) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,46), 2.00000));
ALGEBRAIC(:,60) = ( power(10.0000, CONSTANTS(:,190)).*(ALGEBRAIC(:,56) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,190))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,56), 2.00000));
ALGEBRAIC(:,66) = ( power(10.0000, CONSTANTS(:,196)).*(ALGEBRAIC(:,61) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,196))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,61), 2.00000));
ALGEBRAIC(:,75) = ( power(10.0000, CONSTANTS(:,200)).*(ALGEBRAIC(:,71) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,200))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,71), 2.00000));
ALGEBRAIC(:,81) = ( power(10.0000, CONSTANTS(:,203)).*(ALGEBRAIC(:,77) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,203))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,77), 2.00000));
ALGEBRAIC(:,89) = ( ALGEBRAIC(:,83).*(power(10.0000, CONSTANTS(:,206))+ 2.00000.*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))) - (power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206))+ 2.00000.*power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))).*(power(10.0000, CONSTANTS(:,206))+ 2.00000.*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,83), 2.00000));
ALGEBRAIC(:,98) = ( power(10.0000, CONSTANTS(:,211)).*(ALGEBRAIC(:,94) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,211))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,94), 2.00000));
ALGEBRAIC(:,105) = ( power(10.0000, CONSTANTS(:,217)).*(ALGEBRAIC(:,100) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,100), 2.00000));
ALGEBRAIC(:,115) = ( power(10.0000, CONSTANTS(:,220)).*(ALGEBRAIC(:,110) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,110), 2.00000));
ALGEBRAIC(:,124) = ( power(10.0000, CONSTANTS(:,224)).*(ALGEBRAIC(:,120) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,224))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,120), 2.00000));
ALGEBRAIC(:,130) = ( power(10.0000, CONSTANTS(:,227)).*(ALGEBRAIC(:,126) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,227))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,126), 2.00000));
ALGEBRAIC(:,138) = ( power(10.0000, CONSTANTS(:,230)).*(ALGEBRAIC(:,132) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,132), 2.00000));
ALGEBRAIC(:,149) = ( power(10.0000, CONSTANTS(:,235)).*(ALGEBRAIC(:,143) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,235))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,143), 2.00000));
ALGEBRAIC(:,158) = ( power(10.0000, CONSTANTS(:,240)).*(ALGEBRAIC(:,154) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,240))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,154), 2.00000));
ALGEBRAIC(:,165) = ( power(10.0000, CONSTANTS(:,246)).*(ALGEBRAIC(:,160) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,160), 2.00000));
ALGEBRAIC(:,293) = log(10.0000).*power(10.0000, - ALGEBRAIC(:,1)).*CONSTANTS(:,135).*(1.00000+ ALGEBRAIC(:,8).*STATES(:,4)+ ALGEBRAIC(:,19).*STATES(:,11)+ ALGEBRAIC(:,30).*STATES(:,10)+ ALGEBRAIC(:,41).*STATES(:,6)+ ALGEBRAIC(:,52).*STATES(:,24)+ ALGEBRAIC(:,60).*STATES(:,23)+ ALGEBRAIC(:,66).*STATES(:,3)+ ALGEBRAIC(:,75).*STATES(:,7)+ ALGEBRAIC(:,81).*STATES(:,8)+ ALGEBRAIC(:,89).*STATES(:,9)+ ALGEBRAIC(:,98).*STATES(:,13)+ ALGEBRAIC(:,115).*STATES(:,12)+ ALGEBRAIC(:,105).*STATES(:,14)+ ALGEBRAIC(:,124).*STATES(:,17)+ ALGEBRAIC(:,130).*STATES(:,18)+ ALGEBRAIC(:,138).*STATES(:,19)+ ALGEBRAIC(:,149).*STATES(:,20)+ ALGEBRAIC(:,158).*STATES(:,21)+ ALGEBRAIC(:,165).*STATES(:,22));
ALGEBRAIC(:,295) = ALGEBRAIC(:,294)./(ALGEBRAIC(:,292)+ALGEBRAIC(:,293));
ALGEBRAIC(:,9) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,150)+CONSTANTS(:,151)))./power(ALGEBRAIC(:,2), 2.00000);
ALGEBRAIC(:,20) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,160)+CONSTANTS(:,161)))./power(ALGEBRAIC(:,13), 2.00000);
ALGEBRAIC(:,32) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,170)+CONSTANTS(:,171)))./power(ALGEBRAIC(:,24), 2.00000);
ALGEBRAIC(:,42) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,178)+CONSTANTS(:,179)))./power(ALGEBRAIC(:,36), 2.00000);
ALGEBRAIC(:,53) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,187)+CONSTANTS(:,188)))./power(ALGEBRAIC(:,46), 2.00000);
ALGEBRAIC(:,67) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,196)+CONSTANTS(:,197)))./power(ALGEBRAIC(:,61), 2.00000);
ALGEBRAIC(:,90) = ( - (power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206))+ 2.00000.*power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))).*power(10.0000, CONSTANTS(:,208)))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,83), 2.00000));
ALGEBRAIC(:,106) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,217)+CONSTANTS(:,214)))./power(ALGEBRAIC(:,100), 2.00000);
ALGEBRAIC(:,116) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,220)+CONSTANTS(:,221)))./power(ALGEBRAIC(:,110), 2.00000);
ALGEBRAIC(:,139) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,230)+CONSTANTS(:,231)))./power(ALGEBRAIC(:,132), 2.00000);
ALGEBRAIC(:,150) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,235)+CONSTANTS(:,236)))./power(ALGEBRAIC(:,143), 2.00000);
ALGEBRAIC(:,166) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,246)+CONSTANTS(:,243)))./power(ALGEBRAIC(:,160), 2.00000);
ALGEBRAIC(:,296) = ( ALGEBRAIC(:,9).*STATES(:,4)+ ALGEBRAIC(:,20).*STATES(:,11)+ ALGEBRAIC(:,32).*STATES(:,10)+ ALGEBRAIC(:,42).*STATES(:,6)+ ALGEBRAIC(:,53).*STATES(:,24)+ CONSTANTS(:,27).*STATES(:,23)+ ALGEBRAIC(:,67).*STATES(:,3)+ CONSTANTS(:,28).*STATES(:,7)+ CONSTANTS(:,29).*STATES(:,8)+ ALGEBRAIC(:,90).*STATES(:,9)+ CONSTANTS(:,30).*STATES(:,13)+ ALGEBRAIC(:,116).*STATES(:,12)+ ALGEBRAIC(:,106).*STATES(:,14)+ CONSTANTS(:,31).*STATES(:,17)+ CONSTANTS(:,32).*STATES(:,18)+ ALGEBRAIC(:,139).*STATES(:,19)+ ALGEBRAIC(:,150).*STATES(:,20)+ CONSTANTS(:,33).*STATES(:,21)+ ALGEBRAIC(:,166).*STATES(:,22))./(ALGEBRAIC(:,292)+ALGEBRAIC(:,293));
ALGEBRAIC(:,11) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150)+CONSTANTS(:,151)).*log(10.0000))./power(ALGEBRAIC(:,2), 2.00000);
ALGEBRAIC(:,22) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160)+CONSTANTS(:,161)).*log(10.0000))./power(ALGEBRAIC(:,13), 2.00000);
ALGEBRAIC(:,34) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170)+CONSTANTS(:,171)).*log(10.0000))./power(ALGEBRAIC(:,24), 2.00000);
ALGEBRAIC(:,44) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,178)+CONSTANTS(:,179)).*log(10.0000))./power(ALGEBRAIC(:,36), 2.00000);
ALGEBRAIC(:,55) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187)+CONSTANTS(:,188)).*log(10.0000))./power(ALGEBRAIC(:,46), 2.00000);
ALGEBRAIC(:,69) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,196)+CONSTANTS(:,197)).*log(10.0000))./power(ALGEBRAIC(:,61), 2.00000);
ALGEBRAIC(:,92) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,208)).*( ALGEBRAIC(:,83).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206)).*log(10.0000) - power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207)).*2.00000.*log(10.0000)))./power(ALGEBRAIC(:,83), 2.00000);
ALGEBRAIC(:,108) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217)+CONSTANTS(:,214)).*log(10.0000))./power(ALGEBRAIC(:,100), 2.00000);
ALGEBRAIC(:,118) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220)+CONSTANTS(:,221)).*log(10.0000))./power(ALGEBRAIC(:,110), 2.00000);
ALGEBRAIC(:,141) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230)+CONSTANTS(:,231)).*log(10.0000))./power(ALGEBRAIC(:,132), 2.00000);
ALGEBRAIC(:,152) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,235)+CONSTANTS(:,236)).*log(10.0000))./power(ALGEBRAIC(:,143), 2.00000);
ALGEBRAIC(:,168) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246)+CONSTANTS(:,243)).*log(10.0000))./power(ALGEBRAIC(:,160), 2.00000);
ALGEBRAIC(:,10) = (( ALGEBRAIC(:,2).*power(10.0000, CONSTANTS(:,151)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,151)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,2), 2.00000);
ALGEBRAIC(:,21) = (( ALGEBRAIC(:,13).*power(10.0000, CONSTANTS(:,161)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,161)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,13), 2.00000);
ALGEBRAIC(:,33) = (( ALGEBRAIC(:,24).*power(10.0000, CONSTANTS(:,171)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,171)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,24), 2.00000);
ALGEBRAIC(:,43) = (( ALGEBRAIC(:,36).*power(10.0000, CONSTANTS(:,179)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,179)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,36), 2.00000);
ALGEBRAIC(:,54) = (( ALGEBRAIC(:,46).*power(10.0000, CONSTANTS(:,188)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,188)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,46), 2.00000);
ALGEBRAIC(:,68) = (( ALGEBRAIC(:,61).*power(10.0000, CONSTANTS(:,197)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,197)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,61), 2.00000);
ALGEBRAIC(:,91) = (( ALGEBRAIC(:,83).*power(10.0000, CONSTANTS(:,208)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,208)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,83), 2.00000);
ALGEBRAIC(:,107) = (( ALGEBRAIC(:,100).*power(10.0000, CONSTANTS(:,214)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,214)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,100), 2.00000);
ALGEBRAIC(:,117) = (( ALGEBRAIC(:,110).*power(10.0000, CONSTANTS(:,221)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,221)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,110), 2.00000);
ALGEBRAIC(:,140) = (( ALGEBRAIC(:,132).*power(10.0000, CONSTANTS(:,231)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,231)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,132), 2.00000);
ALGEBRAIC(:,151) = (( ALGEBRAIC(:,143).*power(10.0000, CONSTANTS(:,236)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,236)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,143), 2.00000);
ALGEBRAIC(:,167) = (( ALGEBRAIC(:,160).*power(10.0000, CONSTANTS(:,243)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,243)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,160), 2.00000);
ALGEBRAIC(:,297) = - 1.00000 - ( ALGEBRAIC(:,21).*STATES(:,11)+ ALGEBRAIC(:,33).*STATES(:,10)+ ALGEBRAIC(:,43).*STATES(:,6)+ ALGEBRAIC(:,10).*STATES(:,4)+ ALGEBRAIC(:,54).*STATES(:,24)+ ALGEBRAIC(:,68).*STATES(:,3)+ ALGEBRAIC(:,91).*STATES(:,9)+ ALGEBRAIC(:,107).*STATES(:,14)+ ALGEBRAIC(:,117).*STATES(:,12)+ ALGEBRAIC(:,151).*STATES(:,20)+ ALGEBRAIC(:,140).*STATES(:,19)+ ALGEBRAIC(:,167).*STATES(:,22));
ALGEBRAIC(:,298) = ( ALGEBRAIC(:,22).*STATES(:,11)+ ALGEBRAIC(:,34).*STATES(:,10)+ ALGEBRAIC(:,44).*STATES(:,6)+ ALGEBRAIC(:,11).*STATES(:,4)+ ALGEBRAIC(:,55).*STATES(:,24)+ ALGEBRAIC(:,69).*STATES(:,3)+ ALGEBRAIC(:,92).*STATES(:,9)+ ALGEBRAIC(:,108).*STATES(:,14)+ ALGEBRAIC(:,118).*STATES(:,12)+ ALGEBRAIC(:,152).*STATES(:,20)+ ALGEBRAIC(:,141).*STATES(:,19)+ ALGEBRAIC(:,168).*STATES(:,22))./ALGEBRAIC(:,297);
ALGEBRAIC(:,6) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,151)).*ALGEBRAIC(:,3);
ALGEBRAIC(:,50) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,188)).*ALGEBRAIC(:,47);
ALGEBRAIC(:,64) = (( ALGEBRAIC(:,62).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,197));
ALGEBRAIC(:,87) = (( ALGEBRAIC(:,84).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,208));
ALGEBRAIC(:,103) = (( ALGEBRAIC(:,101).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,214));
ALGEBRAIC(:,113) = (( ALGEBRAIC(:,111).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,221));
ALGEBRAIC(:,136) = (( ALGEBRAIC(:,133).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,231));
ALGEBRAIC(:,147) = (( ALGEBRAIC(:,144).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,236));
ALGEBRAIC(:,163) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,243)).*ALGEBRAIC(:,161);
ALGEBRAIC(:,312) = ( ALGEBRAIC(:,16).*ALGEBRAIC(:,301)+ ALGEBRAIC(:,27).*ALGEBRAIC(:,302)+ ALGEBRAIC(:,39).*ALGEBRAIC(:,303)+ ALGEBRAIC(:,6).*ALGEBRAIC(:,304)+ ALGEBRAIC(:,50).*ALGEBRAIC(:,300)+ ALGEBRAIC(:,64).*ALGEBRAIC(:,305)+ ALGEBRAIC(:,87).*ALGEBRAIC(:,306)+ ALGEBRAIC(:,103).*ALGEBRAIC(:,308)+ ALGEBRAIC(:,113).*ALGEBRAIC(:,307)+ ALGEBRAIC(:,147).*ALGEBRAIC(:,310)+ ALGEBRAIC(:,136).*ALGEBRAIC(:,309)+ ALGEBRAIC(:,163).*ALGEBRAIC(:,311))./ALGEBRAIC(:,297);
RATES(:,2) = piecewise({VOI<=1.00000, ( CONSTANTS(:,34).*(ALGEBRAIC(:,312)+ ALGEBRAIC(:,298).*ALGEBRAIC(:,295)))./(1.00000 - ALGEBRAIC(:,298).*ALGEBRAIC(:,296)) }, CONSTANTS(:,34).*ALGEBRAIC(:,312));
RATES(:,1) = ( CONSTANTS(:,342).*(ALGEBRAIC(:,295)+ ALGEBRAIC(:,312).*ALGEBRAIC(:,296)))./(1.00000 - ALGEBRAIC(:,298).*ALGEBRAIC(:,296));
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(:,221) = piecewise({CONSTANTS(:,39) ~= 45.0000, CONSTANTS(:,6) , VOI<40.0000, 0.00100000 , VOI>=40.0000&VOI<80.0000, 0.00400000 , VOI>=80.0000&VOI<100.000, 0.0100000 , VOI>=100.000, 0.0400000 }, NaN);
ALGEBRAIC(:,222) = 1.00000+STATES(:,5)./CONSTANTS(:,40)+STATES(:,4)./CONSTANTS(:,265)+( STATES(:,5).*STATES(:,4))./( CONSTANTS(:,40).*CONSTANTS(:,41))+STATES(:,5)./CONSTANTS(:,267)+STATES(:,3)./CONSTANTS(:,42)+( STATES(:,5).*STATES(:,3))./( CONSTANTS(:,268).*CONSTANTS(:,267));
ALGEBRAIC(:,1) = piecewise({VOI<=1.00000||VOI>1.00000&CONSTANTS(:,35)==0.00000, STATES(:,1) }, CONSTANTS(:,36));
ALGEBRAIC(:,223) = 1.40400./(1.00000+power(10.0000, 5.94000 - ALGEBRAIC(:,1))+power(10.0000, ALGEBRAIC(:,1) - 7.29000));
ALGEBRAIC(:,225) = (( ALGEBRAIC(:,223).*CONSTANTS(:,264).*STATES(:,4))./( CONSTANTS(:,266).*CONSTANTS(:,265)))./ALGEBRAIC(:,222);
ALGEBRAIC(:,2) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,151))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,145));
ALGEBRAIC(:,61) = 1.00000+ power(10.0000, - ALGEBRAIC(:,1)).*power(10.0000, CONSTANTS(:,196))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,197));
ALGEBRAIC(:,12) = (CONSTANTS(:,153)+ CONSTANTS(:,152).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,152));
ALGEBRAIC(:,70) = (CONSTANTS(:,199)+ CONSTANTS(:,198).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,198));
ALGEBRAIC(:,170) = 655.700+ CONSTANTS(:,249).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1);
ALGEBRAIC(:,183) = (ALGEBRAIC(:,170)+ALGEBRAIC(:,70)) - ALGEBRAIC(:,12);
ALGEBRAIC(:,184) = ( exp( - ALGEBRAIC(:,183)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,61))./ALGEBRAIC(:,2);
ALGEBRAIC(:,224) = ( ALGEBRAIC(:,223).*CONSTANTS(:,264).*CONSTANTS(:,267).*CONSTANTS(:,268))./( CONSTANTS(:,266).*CONSTANTS(:,265).*ALGEBRAIC(:,184));
ALGEBRAIC(:,226) = (( ALGEBRAIC(:,224).*STATES(:,5))./( CONSTANTS(:,267).*CONSTANTS(:,268)))./ALGEBRAIC(:,222);
ALGEBRAIC(:,227) = ALGEBRAIC(:,221).*( STATES(:,5).*ALGEBRAIC(:,225) - STATES(:,3).*ALGEBRAIC(:,226));
ALGEBRAIC(:,228) = 1.00000 - ALGEBRAIC(:,221);
ALGEBRAIC(:,229) = (power(STATES(:,6)./CONSTANTS(:,45), CONSTANTS(:,47))./CONSTANTS(:,46))./(1.00000+power(STATES(:,6)./CONSTANTS(:,45), CONSTANTS(:,47))./CONSTANTS(:,46));
ALGEBRAIC(:,230) = 1.00000+STATES(:,5)./CONSTANTS(:,270)+STATES(:,4)./CONSTANTS(:,43)+STATES(:,5)./CONSTANTS(:,272)+STATES(:,3)./CONSTANTS(:,44)+( STATES(:,5).*STATES(:,4))./( CONSTANTS(:,270).*CONSTANTS(:,269))+( STATES(:,5).*STATES(:,3))./( CONSTANTS(:,271).*CONSTANTS(:,272));
ALGEBRAIC(:,231) = 1.75000./(1.00000+power(10.0000, 6.12000 - ALGEBRAIC(:,1))+power(10.0000, ALGEBRAIC(:,1) - 7.03000));
ALGEBRAIC(:,233) = (( ALGEBRAIC(:,231).*ALGEBRAIC(:,229).*CONSTANTS(:,264).*STATES(:,4))./( CONSTANTS(:,270).*CONSTANTS(:,269)))./ALGEBRAIC(:,230);
ALGEBRAIC(:,232) = ( ALGEBRAIC(:,231).*CONSTANTS(:,264).*CONSTANTS(:,271).*CONSTANTS(:,272))./( CONSTANTS(:,270).*CONSTANTS(:,269).*ALGEBRAIC(:,184));
ALGEBRAIC(:,234) = (( ALGEBRAIC(:,229).*ALGEBRAIC(:,232).*STATES(:,5))./( CONSTANTS(:,271).*CONSTANTS(:,272)))./ALGEBRAIC(:,230);
ALGEBRAIC(:,235) = ALGEBRAIC(:,228).*( STATES(:,5).*ALGEBRAIC(:,233) - STATES(:,3).*ALGEBRAIC(:,234));
ALGEBRAIC(:,237) = - (ALGEBRAIC(:,227)+ALGEBRAIC(:,235));
ALGEBRAIC(:,236) = ( CONSTANTS(:,273).*1.32900)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+6.64000)+power(10.0000, ALGEBRAIC(:,1) - 8.36000));
ALGEBRAIC(:,71) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,200));
ALGEBRAIC(:,76) = (CONSTANTS(:,202)+ CONSTANTS(:,201).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,201));
ALGEBRAIC(:,186) = ALGEBRAIC(:,76) - ALGEBRAIC(:,70);
ALGEBRAIC(:,187) = ( exp( - ALGEBRAIC(:,186)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,71))./ALGEBRAIC(:,61);
ALGEBRAIC(:,238) = ( ALGEBRAIC(:,236).*CONSTANTS(:,275))./( CONSTANTS(:,274).*ALGEBRAIC(:,187));
ALGEBRAIC(:,239) = (( ALGEBRAIC(:,236).*STATES(:,3))./CONSTANTS(:,274) - ( ALGEBRAIC(:,238).*STATES(:,7))./CONSTANTS(:,275))./(1.00000+STATES(:,3)./CONSTANTS(:,274)+STATES(:,7)./CONSTANTS(:,275));
ALGEBRAIC(:,240) = CONSTANTS(:,276)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+6.94000)+power(10.0000, ALGEBRAIC(:,1) - 9.35000));
ALGEBRAIC(:,77) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,203));
ALGEBRAIC(:,82) = (CONSTANTS(:,205)+ CONSTANTS(:,204).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,204));
ALGEBRAIC(:,189) = ALGEBRAIC(:,82) - ALGEBRAIC(:,76);
ALGEBRAIC(:,190) = ( exp( - ALGEBRAIC(:,189)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,77))./ALGEBRAIC(:,71);
ALGEBRAIC(:,241) = (( ALGEBRAIC(:,240).*CONSTANTS(:,277))./CONSTANTS(:,278)).*ALGEBRAIC(:,190);
ALGEBRAIC(:,242) = (( ALGEBRAIC(:,241).*STATES(:,7))./CONSTANTS(:,277) - ( ALGEBRAIC(:,240).*STATES(:,8))./CONSTANTS(:,278))./(1.00000+STATES(:,8)./CONSTANTS(:,278)+STATES(:,7)./CONSTANTS(:,277));
ALGEBRAIC(:,244) = ALGEBRAIC(:,239) - ALGEBRAIC(:,242);
ALGEBRAIC(:,243) = CONSTANTS(:,279)./(1.00000+power(ALGEBRAIC(:,1)./6.80000, - 30.0000));
ALGEBRAIC(:,13) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,161))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,162));
ALGEBRAIC(:,24) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,171))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,165));
ALGEBRAIC(:,83) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206))+power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,208));
ALGEBRAIC(:,23) = (CONSTANTS(:,164)+ CONSTANTS(:,163).*CONSTANTS(:,1).*CONSTANTS(:,131).*log(10.0000).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(16.0000 - CONSTANTS(:,163));
ALGEBRAIC(:,35) = (CONSTANTS(:,173)+ CONSTANTS(:,172).*CONSTANTS(:,1).*CONSTANTS(:,131).*log(10.0000).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,172));
ALGEBRAIC(:,93) = (CONSTANTS(:,210)+ CONSTANTS(:,209).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(16.0000 - CONSTANTS(:,209));
ALGEBRAIC(:,174) = CONSTANTS(:,257)+ CONSTANTS(:,256).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1);
ALGEBRAIC(:,192) = ((ALGEBRAIC(:,93)+ALGEBRAIC(:,35)+ALGEBRAIC(:,174)) - ALGEBRAIC(:,82)) - ALGEBRAIC(:,23);
ALGEBRAIC(:,193) = ( exp( - ALGEBRAIC(:,192)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,83).*ALGEBRAIC(:,24))./( ALGEBRAIC(:,77).*ALGEBRAIC(:,13).*power(10.0000, - ALGEBRAIC(:,1)));
ALGEBRAIC(:,245) = ( ALGEBRAIC(:,243).*CONSTANTS(:,283).*CONSTANTS(:,285))./( CONSTANTS(:,280).*CONSTANTS(:,282).*ALGEBRAIC(:,193));
ALGEBRAIC(:,246) = CONSTANTS(:,54).*power(( ((1.00000+STATES(:,11)./CONSTANTS(:,50))./(1.00000+( CONSTANTS(:,52).*STATES(:,11))./CONSTANTS(:,50))).*(1.00000+( CONSTANTS(:,53).*STATES(:,6))./CONSTANTS(:,51)))./(1.00000+STATES(:,6)./CONSTANTS(:,51)), 4.00000);
ALGEBRAIC(:,247) = (1.00000+STATES(:,8)./CONSTANTS(:,280)).*(1.00000+STATES(:,11)./CONSTANTS(:,282))+STATES(:,9)./CONSTANTS(:,283)+ (STATES(:,10)./CONSTANTS(:,285)).*(1.00000+STATES(:,9)./CONSTANTS(:,283));
ALGEBRAIC(:,248) = (1.00000+STATES(:,8)./CONSTANTS(:,281)).*(1.00000+STATES(:,11)./CONSTANTS(:,284))+STATES(:,9)./CONSTANTS(:,48)+ (STATES(:,10)./CONSTANTS(:,49)).*(1.00000+STATES(:,9)./CONSTANTS(:,48));
ALGEBRAIC(:,249) = ( ((( ALGEBRAIC(:,243).*STATES(:,8).*STATES(:,11))./( CONSTANTS(:,280).*CONSTANTS(:,282)))./ALGEBRAIC(:,247) - (( ALGEBRAIC(:,245).*STATES(:,10).*STATES(:,9))./( CONSTANTS(:,285).*CONSTANTS(:,283)))./ALGEBRAIC(:,247)).*(1.00000+ CONSTANTS(:,286).*ALGEBRAIC(:,246).*power(ALGEBRAIC(:,248)./ALGEBRAIC(:,247), 3.00000)))./(1.00000+ ALGEBRAIC(:,246).*power(ALGEBRAIC(:,248)./ALGEBRAIC(:,247), 4.00000));
ALGEBRAIC(:,251) = ALGEBRAIC(:,242) - ALGEBRAIC(:,249);
ALGEBRAIC(:,250) = ( CONSTANTS(:,287).*1.01300)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+5.32000)+power(10.0000, ALGEBRAIC(:,1) - 9.15000));
ALGEBRAIC(:,94) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,211));
ALGEBRAIC(:,110) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,221));
ALGEBRAIC(:,99) = (CONSTANTS(:,213)+ CONSTANTS(:,212).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,212));
ALGEBRAIC(:,119) = (CONSTANTS(:,223)+ CONSTANTS(:,222).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,222));
ALGEBRAIC(:,195) = (ALGEBRAIC(:,119)+ALGEBRAIC(:,99)) - ALGEBRAIC(:,93);
ALGEBRAIC(:,196) = ( 1.00000.*exp( - ALGEBRAIC(:,195)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,94).*ALGEBRAIC(:,83))./ALGEBRAIC(:,110);
ALGEBRAIC(:,252) = ( ALGEBRAIC(:,250).*CONSTANTS(:,290).*CONSTANTS(:,289))./( CONSTANTS(:,288).*ALGEBRAIC(:,196));
ALGEBRAIC(:,253) = (( ALGEBRAIC(:,250).*STATES(:,9))./CONSTANTS(:,288) - ( ALGEBRAIC(:,252).*STATES(:,13).*STATES(:,12))./( CONSTANTS(:,290).*CONSTANTS(:,289)))./(1.00000+STATES(:,9)./CONSTANTS(:,288)+STATES(:,13)./CONSTANTS(:,290)+STATES(:,12)./CONSTANTS(:,289));
ALGEBRAIC(:,198) = ALGEBRAIC(:,119) - ALGEBRAIC(:,99);
ALGEBRAIC(:,199) = ( exp( - ALGEBRAIC(:,198)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,110))./ALGEBRAIC(:,94);
ALGEBRAIC(:,254) = ( CONSTANTS(:,294).*CONSTANTS(:,293))./( CONSTANTS(:,292).*ALGEBRAIC(:,199));
ALGEBRAIC(:,255) = (( CONSTANTS(:,294).*STATES(:,13))./CONSTANTS(:,292) - ( ALGEBRAIC(:,254).*STATES(:,12))./CONSTANTS(:,293))./(1.00000+STATES(:,13)./CONSTANTS(:,292)+STATES(:,12)./CONSTANTS(:,293));
ALGEBRAIC(:,259) = 1.00000+STATES(:,4)./CONSTANTS(:,304)+STATES(:,13)./CONSTANTS(:,302)+STATES(:,15)./CONSTANTS(:,303)+( STATES(:,13).*STATES(:,15))./( CONSTANTS(:,302).*CONSTANTS(:,303))+( STATES(:,13).*STATES(:,15).*STATES(:,4))./( CONSTANTS(:,302).*CONSTANTS(:,303).*CONSTANTS(:,304))+STATES(:,17)./CONSTANTS(:,305)+STATES(:,16)./CONSTANTS(:,306)+( STATES(:,17).*STATES(:,16))./( CONSTANTS(:,306).*CONSTANTS(:,305));
ALGEBRAIC(:,260) = CONSTANTS(:,301).*0.000700000.*exp( ALGEBRAIC(:,1).*0.897900);
ALGEBRAIC(:,120) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,224));
ALGEBRAIC(:,125) = (CONSTANTS(:,226)+ CONSTANTS(:,225).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(16.0000 - CONSTANTS(:,225));
ALGEBRAIC(:,171) = (CONSTANTS(:,251)+ CONSTANTS(:,250).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(1.00000 - CONSTANTS(:,250));
ALGEBRAIC(:,172) = (CONSTANTS(:,253)+ CONSTANTS(:,252).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,252));
ALGEBRAIC(:,201) = (((ALGEBRAIC(:,125)+ALGEBRAIC(:,172)+ALGEBRAIC(:,174)) - ALGEBRAIC(:,12)) - ALGEBRAIC(:,99)) - ALGEBRAIC(:,171);
ALGEBRAIC(:,202) = ( exp( - ALGEBRAIC(:,201)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,120))./( ALGEBRAIC(:,2).*ALGEBRAIC(:,94).*power(10.0000, - ALGEBRAIC(:,1)).*1.00000);
ALGEBRAIC(:,261) = ( ALGEBRAIC(:,260).*CONSTANTS(:,305).*CONSTANTS(:,306))./( CONSTANTS(:,302).*CONSTANTS(:,304).*CONSTANTS(:,303).*ALGEBRAIC(:,202));
ALGEBRAIC(:,262) = (( ALGEBRAIC(:,260).*STATES(:,13).*STATES(:,15).*STATES(:,4))./( CONSTANTS(:,303).*CONSTANTS(:,302).*CONSTANTS(:,304)) - ( ALGEBRAIC(:,261).*STATES(:,17).*STATES(:,16))./( CONSTANTS(:,305).*CONSTANTS(:,306)))./ALGEBRAIC(:,259);
ALGEBRAIC(:,264) = (ALGEBRAIC(:,253) - ALGEBRAIC(:,255)) - ALGEBRAIC(:,262);
ALGEBRAIC(:,126) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+6.21000);
ALGEBRAIC(:,131) = (CONSTANTS(:,229)+ CONSTANTS(:,228).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,228));
ALGEBRAIC(:,207) = ((ALGEBRAIC(:,23)+ALGEBRAIC(:,131)) - ALGEBRAIC(:,125)) - ALGEBRAIC(:,35);
ALGEBRAIC(:,208) = ( exp( - ALGEBRAIC(:,207)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,13).*ALGEBRAIC(:,126))./( ALGEBRAIC(:,120).*ALGEBRAIC(:,24));
ALGEBRAIC(:,263) = (( CONSTANTS(:,312).*CONSTANTS(:,308).*CONSTANTS(:,309))./( CONSTANTS(:,310).*CONSTANTS(:,311))).*ALGEBRAIC(:,208);
ALGEBRAIC(:,265) = 1.00000+STATES(:,10)./CONSTANTS(:,309)+STATES(:,17)./CONSTANTS(:,308)+( STATES(:,17).*STATES(:,10))./( CONSTANTS(:,308).*CONSTANTS(:,309))+STATES(:,18)./CONSTANTS(:,310)+STATES(:,11)./CONSTANTS(:,311)+( STATES(:,18).*STATES(:,11))./( CONSTANTS(:,310).*CONSTANTS(:,311));
ALGEBRAIC(:,266) = (( ALGEBRAIC(:,263).*STATES(:,17).*STATES(:,10))./( CONSTANTS(:,309).*CONSTANTS(:,308)) - ( CONSTANTS(:,312).*STATES(:,11).*STATES(:,18))./( CONSTANTS(:,311).*CONSTANTS(:,310)))./ALGEBRAIC(:,265);
ALGEBRAIC(:,268) = ALGEBRAIC(:,262) - ALGEBRAIC(:,266);
ALGEBRAIC(:,267) = ( CONSTANTS(:,313).*0.989000)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+5.62000)+power(10.0000, ALGEBRAIC(:,1) - 8.74000));
ALGEBRAIC(:,132) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,231))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,232));
ALGEBRAIC(:,142) = (CONSTANTS(:,234)+ CONSTANTS(:,233).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,233));
ALGEBRAIC(:,210) = ALGEBRAIC(:,142) - ALGEBRAIC(:,131);
ALGEBRAIC(:,211) = ( exp( - ALGEBRAIC(:,210)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,132))./ALGEBRAIC(:,126);
ALGEBRAIC(:,269) = ( ALGEBRAIC(:,267).*CONSTANTS(:,315))./( CONSTANTS(:,314).*ALGEBRAIC(:,211));
ALGEBRAIC(:,270) = (( ALGEBRAIC(:,267).*STATES(:,18))./CONSTANTS(:,314) - ( ALGEBRAIC(:,269).*STATES(:,19))./CONSTANTS(:,315))./(1.00000+STATES(:,18)./CONSTANTS(:,314)+STATES(:,19)./CONSTANTS(:,315));
ALGEBRAIC(:,272) = ALGEBRAIC(:,266) - ALGEBRAIC(:,270);
ALGEBRAIC(:,256) = (1.00000+STATES(:,14)./CONSTANTS(:,296)+STATES(:,16)./CONSTANTS(:,299)).*(1.00000+STATES(:,12)./CONSTANTS(:,298)+STATES(:,15)./CONSTANTS(:,299));
ALGEBRAIC(:,100) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,214));
ALGEBRAIC(:,109) = (CONSTANTS(:,219)+ CONSTANTS(:,218).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,218));
ALGEBRAIC(:,204) = ((ALGEBRAIC(:,174)+ALGEBRAIC(:,172)+ALGEBRAIC(:,119)) - ALGEBRAIC(:,171)) - ALGEBRAIC(:,109);
ALGEBRAIC(:,205) = ( exp( - ALGEBRAIC(:,204)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,110))./( ALGEBRAIC(:,100).*power(10.0000, - ALGEBRAIC(:,1)));
ALGEBRAIC(:,257) = ( CONSTANTS(:,300).*CONSTANTS(:,296).*CONSTANTS(:,297).*ALGEBRAIC(:,205))./( CONSTANTS(:,298).*CONSTANTS(:,299));
ALGEBRAIC(:,258) = (( ALGEBRAIC(:,257).*STATES(:,14).*STATES(:,15))./( CONSTANTS(:,296).*CONSTANTS(:,297)) - ( CONSTANTS(:,300).*STATES(:,12).*STATES(:,16))./( CONSTANTS(:,298).*CONSTANTS(:,299)))./ALGEBRAIC(:,256);
ALGEBRAIC(:,277) = CONSTANTS(:,325).*( - 0.113400.*ALGEBRAIC(:,1)+1.60690);
ALGEBRAIC(:,154) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,240));
ALGEBRAIC(:,160) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,243));
ALGEBRAIC(:,159) = (CONSTANTS(:,242)+ CONSTANTS(:,241).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(1.00000 - CONSTANTS(:,241));
ALGEBRAIC(:,169) = (CONSTANTS(:,248)+ CONSTANTS(:,247).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(1.00000 - CONSTANTS(:,247));
ALGEBRAIC(:,219) = (((ALGEBRAIC(:,169)+ALGEBRAIC(:,171)) - ALGEBRAIC(:,159)) - ALGEBRAIC(:,172)) - ALGEBRAIC(:,174);
ALGEBRAIC(:,220) = ( exp( - ALGEBRAIC(:,219)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,160).*power(10.0000, - ALGEBRAIC(:,1)))./ALGEBRAIC(:,154);
ALGEBRAIC(:,278) = ( ALGEBRAIC(:,277).*CONSTANTS(:,328).*CONSTANTS(:,329))./( CONSTANTS(:,326).*CONSTANTS(:,327).*ALGEBRAIC(:,220));
ALGEBRAIC(:,279) = (( ALGEBRAIC(:,277).*STATES(:,21).*STATES(:,16))./( CONSTANTS(:,326).*CONSTANTS(:,327)) - ( ALGEBRAIC(:,278).*STATES(:,22).*STATES(:,15))./( CONSTANTS(:,328).*CONSTANTS(:,329)))./(1.00000+STATES(:,21)./CONSTANTS(:,326)+STATES(:,16)./CONSTANTS(:,327)+( STATES(:,21).*STATES(:,16))./( CONSTANTS(:,326).*CONSTANTS(:,327))+STATES(:,22)./CONSTANTS(:,328)+STATES(:,15)./CONSTANTS(:,329)+( STATES(:,22).*STATES(:,15))./( CONSTANTS(:,328).*CONSTANTS(:,329)));
ALGEBRAIC(:,282) = ( - ALGEBRAIC(:,262) - ALGEBRAIC(:,258))+ALGEBRAIC(:,279);
ALGEBRAIC(:,283) = (ALGEBRAIC(:,262)+ALGEBRAIC(:,258)) - ALGEBRAIC(:,279);
ALGEBRAIC(:,274) = ( CONSTANTS(:,320).*1.05000)./(1.00000+power(10.0000, - ALGEBRAIC(:,1)+5.58000)+power(10.0000, ALGEBRAIC(:,1) - 8.79000));
ALGEBRAIC(:,143) = 1.00000+power(10.0000, CONSTANTS(:,235) - ALGEBRAIC(:,1))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,236))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,237));
ALGEBRAIC(:,153) = (CONSTANTS(:,239)+ CONSTANTS(:,238).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(9.00000 - CONSTANTS(:,238));
ALGEBRAIC(:,216) = (((ALGEBRAIC(:,159)+ALGEBRAIC(:,23)) - ALGEBRAIC(:,174)) - ALGEBRAIC(:,153)) - ALGEBRAIC(:,35);
ALGEBRAIC(:,217) = ( exp( - ALGEBRAIC(:,216)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,154).*ALGEBRAIC(:,13).*power(10.0000, - ALGEBRAIC(:,1)))./( ALGEBRAIC(:,143).*ALGEBRAIC(:,24));
ALGEBRAIC(:,275) = ( ALGEBRAIC(:,274).*CONSTANTS(:,323).*CONSTANTS(:,324))./( CONSTANTS(:,321).*CONSTANTS(:,322).*ALGEBRAIC(:,217));
ALGEBRAIC(:,276) = (( ALGEBRAIC(:,274).*STATES(:,20).*STATES(:,10))./( CONSTANTS(:,321).*CONSTANTS(:,322)) - ( ALGEBRAIC(:,275).*STATES(:,21).*STATES(:,11))./( CONSTANTS(:,323).*CONSTANTS(:,324)))./(1.00000+STATES(:,20)./CONSTANTS(:,321)+STATES(:,10)./CONSTANTS(:,322)+( STATES(:,20).*STATES(:,10))./( CONSTANTS(:,321).*CONSTANTS(:,322))+STATES(:,11)./CONSTANTS(:,324)+STATES(:,21)./CONSTANTS(:,323)+( STATES(:,21).*STATES(:,11))./( CONSTANTS(:,323).*CONSTANTS(:,324)));
ALGEBRAIC(:,284) = ALGEBRAIC(:,276) - ALGEBRAIC(:,279);
ALGEBRAIC(:,46) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,188))+ (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,182));
ALGEBRAIC(:,56) = 1.00000+power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,190));
ALGEBRAIC(:,177) = ( exp( - CONSTANTS(:,263)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*power(10.0000, - ALGEBRAIC(:,1)).*ALGEBRAIC(:,13).*ALGEBRAIC(:,56))./( ALGEBRAIC(:,46).*ALGEBRAIC(:,24));
ALGEBRAIC(:,286) = ( (CONSTANTS(:,330)./ALGEBRAIC(:,177)).*CONSTANTS(:,332).*CONSTANTS(:,334))./( CONSTANTS(:,333).*CONSTANTS(:,331));
ALGEBRAIC(:,287) = (( ALGEBRAIC(:,286).*STATES(:,11).*STATES(:,23))./( CONSTANTS(:,332).*CONSTANTS(:,334)) - ( CONSTANTS(:,330).*STATES(:,10).*STATES(:,24))./( CONSTANTS(:,333).*CONSTANTS(:,331)))./(1.00000+STATES(:,10)./CONSTANTS(:,333)+STATES(:,24)./CONSTANTS(:,56)+( STATES(:,24).*STATES(:,10))./( CONSTANTS(:,333).*CONSTANTS(:,331))+STATES(:,11)./CONSTANTS(:,332)+( STATES(:,23).*STATES(:,11))./( CONSTANTS(:,334).*CONSTANTS(:,332)));
ALGEBRAIC(:,290) = - ALGEBRAIC(:,287);
ALGEBRAIC(:,14) = 1.00000./ALGEBRAIC(:,13);
ALGEBRAIC(:,15) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160)).*ALGEBRAIC(:,14);
ALGEBRAIC(:,16) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,161)).*ALGEBRAIC(:,14);
ALGEBRAIC(:,17) = (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,162)).*ALGEBRAIC(:,14);
ALGEBRAIC(:,18) = 0.00000.*ALGEBRAIC(:,14)+ 1.00000.*ALGEBRAIC(:,15)+ 0.00000.*ALGEBRAIC(:,16)+ 0.00000.*ALGEBRAIC(:,17);
ALGEBRAIC(:,25) = 1.00000./ALGEBRAIC(:,24);
ALGEBRAIC(:,26) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170)).*ALGEBRAIC(:,25);
ALGEBRAIC(:,27) = (( ALGEBRAIC(:,25).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,171));
ALGEBRAIC(:,28) = (( ALGEBRAIC(:,25).*CONSTANTS(:,134))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,165));
ALGEBRAIC(:,29) = 0.00000.*ALGEBRAIC(:,25)+ 1.00000.*ALGEBRAIC(:,26)+ 0.00000.*ALGEBRAIC(:,27)+ 0.00000.*ALGEBRAIC(:,28);
ALGEBRAIC(:,47) = 1.00000./ALGEBRAIC(:,46);
ALGEBRAIC(:,49) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187)).*ALGEBRAIC(:,47);
ALGEBRAIC(:,51) = ALGEBRAIC(:,49);
ALGEBRAIC(:,57) = 1.00000./ALGEBRAIC(:,56);
ALGEBRAIC(:,58) = ALGEBRAIC(:,57).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,190));
ALGEBRAIC(:,59) = 0.00000.*ALGEBRAIC(:,57)+ 1.00000.*ALGEBRAIC(:,58);
ALGEBRAIC(:,175) = (((ALGEBRAIC(:,18)+ALGEBRAIC(:,59)) - ALGEBRAIC(:,51)) - ALGEBRAIC(:,29))+(((CONSTANTS(:,163)+CONSTANTS(:,191)) - CONSTANTS(:,189)) - CONSTANTS(:,172));
ALGEBRAIC(:,36) = 1.00000+power(10.0000, CONSTANTS(:,178) - ALGEBRAIC(:,1))+ (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,179));
ALGEBRAIC(:,37) = 1.00000./ALGEBRAIC(:,36);
ALGEBRAIC(:,38) = ALGEBRAIC(:,37).*power(10.0000, CONSTANTS(:,178) - ALGEBRAIC(:,1));
ALGEBRAIC(:,39) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,179)).*ALGEBRAIC(:,37);
ALGEBRAIC(:,40) = 0.00000.*ALGEBRAIC(:,37)+ALGEBRAIC(:,38)+ 0.00000.*ALGEBRAIC(:,39);
ALGEBRAIC(:,179) = (( 2.00000.*ALGEBRAIC(:,29) - ALGEBRAIC(:,18)) - ALGEBRAIC(:,40))+(( 2.00000.*CONSTANTS(:,172) - CONSTANTS(:,163)) - CONSTANTS(:,180));
ALGEBRAIC(:,3) = 1.00000./ALGEBRAIC(:,2);
ALGEBRAIC(:,4) = power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150)).*ALGEBRAIC(:,3);
ALGEBRAIC(:,7) = 1.00000.*ALGEBRAIC(:,4);
ALGEBRAIC(:,62) = 1.00000./ALGEBRAIC(:,61);
ALGEBRAIC(:,63) = ALGEBRAIC(:,62).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,196));
ALGEBRAIC(:,65) = ALGEBRAIC(:,63);
ALGEBRAIC(:,182) = (ALGEBRAIC(:,65) - ALGEBRAIC(:,7))+(1.00000 - CONSTANTS(:,152));
ALGEBRAIC(:,72) = 1.00000./ALGEBRAIC(:,71);
ALGEBRAIC(:,73) = ALGEBRAIC(:,72).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,200));
ALGEBRAIC(:,74) = ALGEBRAIC(:,73);
ALGEBRAIC(:,185) = (ALGEBRAIC(:,74) - ALGEBRAIC(:,65))+(CONSTANTS(:,201) - CONSTANTS(:,198));
ALGEBRAIC(:,78) = 1.00000./ALGEBRAIC(:,77);
ALGEBRAIC(:,79) = ALGEBRAIC(:,78).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,203));
ALGEBRAIC(:,80) = ALGEBRAIC(:,79);
ALGEBRAIC(:,188) = (ALGEBRAIC(:,80) - ALGEBRAIC(:,74))+(CONSTANTS(:,204) - CONSTANTS(:,201));
ALGEBRAIC(:,84) = 1.00000./ALGEBRAIC(:,83);
ALGEBRAIC(:,85) = ALGEBRAIC(:,84).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206));
ALGEBRAIC(:,86) = ALGEBRAIC(:,84).*power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207));
ALGEBRAIC(:,88) = ALGEBRAIC(:,85)+ 2.00000.*ALGEBRAIC(:,86);
ALGEBRAIC(:,191) = (((((ALGEBRAIC(:,29)+ALGEBRAIC(:,88)) - ALGEBRAIC(:,80)) - ALGEBRAIC(:,18))+CONSTANTS(:,172)+CONSTANTS(:,209)) - CONSTANTS(:,204)) - CONSTANTS(:,163);
ALGEBRAIC(:,95) = 1.00000./ALGEBRAIC(:,94);
ALGEBRAIC(:,96) = ALGEBRAIC(:,95).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,211));
ALGEBRAIC(:,97) = ALGEBRAIC(:,96);
ALGEBRAIC(:,111) = 1.00000./ALGEBRAIC(:,110);
ALGEBRAIC(:,112) = ALGEBRAIC(:,111).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220));
ALGEBRAIC(:,114) = ALGEBRAIC(:,112);
ALGEBRAIC(:,194) = ((ALGEBRAIC(:,114)+ALGEBRAIC(:,97)) - ALGEBRAIC(:,88))+((CONSTANTS(:,222)+CONSTANTS(:,212)) - CONSTANTS(:,209));
ALGEBRAIC(:,197) = (ALGEBRAIC(:,114) - ALGEBRAIC(:,97))+(CONSTANTS(:,222) - CONSTANTS(:,212));
ALGEBRAIC(:,121) = 1.00000./ALGEBRAIC(:,120);
ALGEBRAIC(:,122) = ALGEBRAIC(:,121).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,224));
ALGEBRAIC(:,123) = ALGEBRAIC(:,122);
ALGEBRAIC(:,200) = ((ALGEBRAIC(:,123) - ALGEBRAIC(:,97)) - ALGEBRAIC(:,7))+((((CONSTANTS(:,225)+CONSTANTS(:,252)) - CONSTANTS(:,212)) - CONSTANTS(:,152)) - CONSTANTS(:,250));
ALGEBRAIC(:,101) = 1.00000./ALGEBRAIC(:,100);
ALGEBRAIC(:,102) = ALGEBRAIC(:,101).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217));
ALGEBRAIC(:,104) = ALGEBRAIC(:,102);
ALGEBRAIC(:,203) = (ALGEBRAIC(:,114) - ALGEBRAIC(:,104))+(((CONSTANTS(:,222)+CONSTANTS(:,252)) - CONSTANTS(:,250)) - CONSTANTS(:,218));
ALGEBRAIC(:,127) = 1.00000./ALGEBRAIC(:,126);
ALGEBRAIC(:,128) = ALGEBRAIC(:,127).*power(10.0000, - ALGEBRAIC(:,1)+6.21000);
ALGEBRAIC(:,129) = ALGEBRAIC(:,128);
ALGEBRAIC(:,206) = (((ALGEBRAIC(:,129)+ALGEBRAIC(:,18)) - ALGEBRAIC(:,123)) - ALGEBRAIC(:,29))+(((CONSTANTS(:,228)+CONSTANTS(:,163)) - CONSTANTS(:,225)) - CONSTANTS(:,172));
ALGEBRAIC(:,133) = 1.00000./ALGEBRAIC(:,132);
ALGEBRAIC(:,134) = ALGEBRAIC(:,133).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230));
ALGEBRAIC(:,137) = ALGEBRAIC(:,134);
ALGEBRAIC(:,209) = (ALGEBRAIC(:,137) - ALGEBRAIC(:,129))+(CONSTANTS(:,233) - CONSTANTS(:,228));
ALGEBRAIC(:,144) = 1.00000./ALGEBRAIC(:,143);
ALGEBRAIC(:,145) = ALGEBRAIC(:,144).*power(10.0000, CONSTANTS(:,235) - ALGEBRAIC(:,1));
ALGEBRAIC(:,148) = ALGEBRAIC(:,145);
ALGEBRAIC(:,212) = (ALGEBRAIC(:,148) - ALGEBRAIC(:,137))+((CONSTANTS(:,254)+CONSTANTS(:,238)) - CONSTANTS(:,233));
ALGEBRAIC(:,155) = 1.00000./ALGEBRAIC(:,154);
ALGEBRAIC(:,156) = ALGEBRAIC(:,155).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,240));
ALGEBRAIC(:,157) = ALGEBRAIC(:,156);
ALGEBRAIC(:,215) = (((ALGEBRAIC(:,157)+ALGEBRAIC(:,18)) - ALGEBRAIC(:,148)) - ALGEBRAIC(:,29))+(((CONSTANTS(:,241)+CONSTANTS(:,163)) - CONSTANTS(:,238)) - CONSTANTS(:,172));
ALGEBRAIC(:,161) = 1.00000./ALGEBRAIC(:,160);
ALGEBRAIC(:,162) = ALGEBRAIC(:,161).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246));
ALGEBRAIC(:,164) = ALGEBRAIC(:,162);
ALGEBRAIC(:,218) = (ALGEBRAIC(:,164) - ALGEBRAIC(:,157))+(((CONSTANTS(:,247)+CONSTANTS(:,250)) - CONSTANTS(:,241)) - CONSTANTS(:,252));
ALGEBRAIC(:,173) = (CONSTANTS(:,255)+ CONSTANTS(:,254).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(0.00000 - CONSTANTS(:,254));
ALGEBRAIC(:,213) = (ALGEBRAIC(:,173)+ALGEBRAIC(:,153)) - ALGEBRAIC(:,142);
ALGEBRAIC(:,214) = ( exp( - ALGEBRAIC(:,213)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,143))./ALGEBRAIC(:,132);
ALGEBRAIC(:,271) = ( CONSTANTS(:,319).*CONSTANTS(:,318))./( CONSTANTS(:,317).*ALGEBRAIC(:,214));
ALGEBRAIC(:,273) = (( CONSTANTS(:,319).*STATES(:,19))./CONSTANTS(:,317) - ( ALGEBRAIC(:,271).*STATES(:,20))./CONSTANTS(:,318))./(1.00000+STATES(:,20)./CONSTANTS(:,318)+STATES(:,19)./CONSTANTS(:,317));
ALGEBRAIC(:,45) = (CONSTANTS(:,181)+ CONSTANTS(:,180).*log(10.0000).*CONSTANTS(:,1).*CONSTANTS(:,131).*ALGEBRAIC(:,1)) - CONSTANTS(:,144).*(4.00000 - CONSTANTS(:,180));
ALGEBRAIC(:,180) = ( 2.00000.*ALGEBRAIC(:,35) - ALGEBRAIC(:,23)) - ALGEBRAIC(:,45);
ALGEBRAIC(:,181) = ( exp( - ALGEBRAIC(:,180)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*power(ALGEBRAIC(:,24), 2.00000))./( ALGEBRAIC(:,13).*ALGEBRAIC(:,36));
ALGEBRAIC(:,288) = ( CONSTANTS(:,335).*power(CONSTANTS(:,338), 2.00000))./( CONSTANTS(:,336).*CONSTANTS(:,337).*ALGEBRAIC(:,181));
ALGEBRAIC(:,291) = (( CONSTANTS(:,335).*STATES(:,11).*STATES(:,6))./( CONSTANTS(:,337).*CONSTANTS(:,336)) - ALGEBRAIC(:,288).*power(STATES(:,10)./CONSTANTS(:,338), 2.00000))./(1.00000+STATES(:,11)./CONSTANTS(:,337)+STATES(:,6)./CONSTANTS(:,336)+( STATES(:,11).*STATES(:,6))./( CONSTANTS(:,337).*CONSTANTS(:,336))+( 2.00000.*STATES(:,10))./CONSTANTS(:,338)+power(STATES(:,10), 2.00000)./power(CONSTANTS(:,338), 2.00000));
ALGEBRAIC(:,294) = ALGEBRAIC(:,175).* - ALGEBRAIC(:,287)+ ALGEBRAIC(:,179).*ALGEBRAIC(:,291)+ ALGEBRAIC(:,185).*ALGEBRAIC(:,239)+ ALGEBRAIC(:,182).*(ALGEBRAIC(:,227)+ALGEBRAIC(:,235))+ ALGEBRAIC(:,188).*ALGEBRAIC(:,242)+ ALGEBRAIC(:,191).*ALGEBRAIC(:,249)+ ALGEBRAIC(:,194).*ALGEBRAIC(:,253)+ ALGEBRAIC(:,197).*ALGEBRAIC(:,255)+ ALGEBRAIC(:,200).*ALGEBRAIC(:,262)+ ALGEBRAIC(:,206).*ALGEBRAIC(:,266)+ ALGEBRAIC(:,209).*ALGEBRAIC(:,270)+ ALGEBRAIC(:,212).*ALGEBRAIC(:,273)+ ALGEBRAIC(:,215).*ALGEBRAIC(:,276)+ ALGEBRAIC(:,218).*ALGEBRAIC(:,279)+ ALGEBRAIC(:,203).*ALGEBRAIC(:,258);
ALGEBRAIC(:,300) = ALGEBRAIC(:,287);
ALGEBRAIC(:,280) = piecewise({CONSTANTS(:,39) ~= 5.00000&CONSTANTS(:,39) ~= 45.0000, CONSTANTS(:,23) , CONSTANTS(:,39) ~= 5.00000&CONSTANTS(:,39)==45.0000&VOI>=110.000, CONSTANTS(:,23) , CONSTANTS(:,39)==5.00000&VOI>30.0000, CONSTANTS(:,23) }, 0.00000);
ALGEBRAIC(:,285) = ( ALGEBRAIC(:,280).*STATES(:,11))./(CONSTANTS(:,55)+STATES(:,11));
ALGEBRAIC(:,301) = ((( - ALGEBRAIC(:,287) - ALGEBRAIC(:,291)) - ALGEBRAIC(:,249))+ALGEBRAIC(:,266)+ALGEBRAIC(:,276)) - ALGEBRAIC(:,285);
ALGEBRAIC(:,302) = (((ALGEBRAIC(:,287)+ 2.00000.*ALGEBRAIC(:,291)+ALGEBRAIC(:,249)) - ALGEBRAIC(:,266)) - ALGEBRAIC(:,276))+ALGEBRAIC(:,285);
ALGEBRAIC(:,303) = - ALGEBRAIC(:,291);
ALGEBRAIC(:,304) = ( - (ALGEBRAIC(:,227)+ALGEBRAIC(:,235)) - ALGEBRAIC(:,262))+ALGEBRAIC(:,285);
ALGEBRAIC(:,305) = (ALGEBRAIC(:,227)+ALGEBRAIC(:,235)) - ALGEBRAIC(:,239);
ALGEBRAIC(:,306) = ALGEBRAIC(:,249) - ALGEBRAIC(:,253);
ALGEBRAIC(:,307) = ALGEBRAIC(:,253)+ALGEBRAIC(:,255)+ALGEBRAIC(:,258);
ALGEBRAIC(:,308) = - ALGEBRAIC(:,258);
ALGEBRAIC(:,309) = ALGEBRAIC(:,270) - ALGEBRAIC(:,273);
ALGEBRAIC(:,310) = ALGEBRAIC(:,273) - ALGEBRAIC(:,276);
ALGEBRAIC(:,311) = ALGEBRAIC(:,279);
ALGEBRAIC(:,292) = ( log(10.0000).*CONSTANTS(:,339).*power(10.0000, - ALGEBRAIC(:,1) - 6.87000))./power(power(10.0000, - ALGEBRAIC(:,1))+power(10.0000, - 6.87000), 2.00000)+( log(10.0000).*CONSTANTS(:,340).*power(10.0000, - ALGEBRAIC(:,1) - 8.30000))./power(power(10.0000, - ALGEBRAIC(:,1))+power(10.0000, - 8.30000), 2.00000)+( log(10.0000).*CONSTANTS(:,341).*power(10.0000, - ALGEBRAIC(:,1) - 4.80000))./power(power(10.0000, - ALGEBRAIC(:,1))+power(10.0000, - 4.80000), 2.00000);
ALGEBRAIC(:,8) = ( power(10.0000, CONSTANTS(:,150)).*(ALGEBRAIC(:,2) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,2), 2.00000));
ALGEBRAIC(:,19) = ( power(10.0000, CONSTANTS(:,160)).*(ALGEBRAIC(:,13) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,13), 2.00000));
ALGEBRAIC(:,30) = ( power(10.0000, CONSTANTS(:,170)).*(ALGEBRAIC(:,24) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,24), 2.00000));
ALGEBRAIC(:,41) = ( power(10.0000, CONSTANTS(:,178)).*(ALGEBRAIC(:,36) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,178))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,36), 2.00000));
ALGEBRAIC(:,52) = ( power(10.0000, CONSTANTS(:,187)).*(ALGEBRAIC(:,46) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,46), 2.00000));
ALGEBRAIC(:,60) = ( power(10.0000, CONSTANTS(:,190)).*(ALGEBRAIC(:,56) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,190))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,56), 2.00000));
ALGEBRAIC(:,66) = ( power(10.0000, CONSTANTS(:,196)).*(ALGEBRAIC(:,61) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,196))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,61), 2.00000));
ALGEBRAIC(:,75) = ( power(10.0000, CONSTANTS(:,200)).*(ALGEBRAIC(:,71) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,200))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,71), 2.00000));
ALGEBRAIC(:,81) = ( power(10.0000, CONSTANTS(:,203)).*(ALGEBRAIC(:,77) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,203))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,77), 2.00000));
ALGEBRAIC(:,89) = ( ALGEBRAIC(:,83).*(power(10.0000, CONSTANTS(:,206))+ 2.00000.*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))) - (power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206))+ 2.00000.*power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))).*(power(10.0000, CONSTANTS(:,206))+ 2.00000.*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,83), 2.00000));
ALGEBRAIC(:,98) = ( power(10.0000, CONSTANTS(:,211)).*(ALGEBRAIC(:,94) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,211))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,94), 2.00000));
ALGEBRAIC(:,105) = ( power(10.0000, CONSTANTS(:,217)).*(ALGEBRAIC(:,100) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,100), 2.00000));
ALGEBRAIC(:,115) = ( power(10.0000, CONSTANTS(:,220)).*(ALGEBRAIC(:,110) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,110), 2.00000));
ALGEBRAIC(:,124) = ( power(10.0000, CONSTANTS(:,224)).*(ALGEBRAIC(:,120) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,224))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,120), 2.00000));
ALGEBRAIC(:,130) = ( power(10.0000, CONSTANTS(:,227)).*(ALGEBRAIC(:,126) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,227))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,126), 2.00000));
ALGEBRAIC(:,138) = ( power(10.0000, CONSTANTS(:,230)).*(ALGEBRAIC(:,132) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,132), 2.00000));
ALGEBRAIC(:,149) = ( power(10.0000, CONSTANTS(:,235)).*(ALGEBRAIC(:,143) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,235))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,143), 2.00000));
ALGEBRAIC(:,158) = ( power(10.0000, CONSTANTS(:,240)).*(ALGEBRAIC(:,154) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,240))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,154), 2.00000));
ALGEBRAIC(:,165) = ( power(10.0000, CONSTANTS(:,246)).*(ALGEBRAIC(:,160) - power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246))))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,160), 2.00000));
ALGEBRAIC(:,293) = log(10.0000).*power(10.0000, - ALGEBRAIC(:,1)).*CONSTANTS(:,135).*(1.00000+ ALGEBRAIC(:,8).*STATES(:,4)+ ALGEBRAIC(:,19).*STATES(:,11)+ ALGEBRAIC(:,30).*STATES(:,10)+ ALGEBRAIC(:,41).*STATES(:,6)+ ALGEBRAIC(:,52).*STATES(:,24)+ ALGEBRAIC(:,60).*STATES(:,23)+ ALGEBRAIC(:,66).*STATES(:,3)+ ALGEBRAIC(:,75).*STATES(:,7)+ ALGEBRAIC(:,81).*STATES(:,8)+ ALGEBRAIC(:,89).*STATES(:,9)+ ALGEBRAIC(:,98).*STATES(:,13)+ ALGEBRAIC(:,115).*STATES(:,12)+ ALGEBRAIC(:,105).*STATES(:,14)+ ALGEBRAIC(:,124).*STATES(:,17)+ ALGEBRAIC(:,130).*STATES(:,18)+ ALGEBRAIC(:,138).*STATES(:,19)+ ALGEBRAIC(:,149).*STATES(:,20)+ ALGEBRAIC(:,158).*STATES(:,21)+ ALGEBRAIC(:,165).*STATES(:,22));
ALGEBRAIC(:,295) = ALGEBRAIC(:,294)./(ALGEBRAIC(:,292)+ALGEBRAIC(:,293));
ALGEBRAIC(:,9) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,150)+CONSTANTS(:,151)))./power(ALGEBRAIC(:,2), 2.00000);
ALGEBRAIC(:,20) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,160)+CONSTANTS(:,161)))./power(ALGEBRAIC(:,13), 2.00000);
ALGEBRAIC(:,32) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,170)+CONSTANTS(:,171)))./power(ALGEBRAIC(:,24), 2.00000);
ALGEBRAIC(:,42) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,178)+CONSTANTS(:,179)))./power(ALGEBRAIC(:,36), 2.00000);
ALGEBRAIC(:,53) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,187)+CONSTANTS(:,188)))./power(ALGEBRAIC(:,46), 2.00000);
ALGEBRAIC(:,67) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,196)+CONSTANTS(:,197)))./power(ALGEBRAIC(:,61), 2.00000);
ALGEBRAIC(:,90) = ( - (power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206))+ 2.00000.*power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207))).*power(10.0000, CONSTANTS(:,208)))./( CONSTANTS(:,135).*power(ALGEBRAIC(:,83), 2.00000));
ALGEBRAIC(:,106) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,217)+CONSTANTS(:,214)))./power(ALGEBRAIC(:,100), 2.00000);
ALGEBRAIC(:,116) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,220)+CONSTANTS(:,221)))./power(ALGEBRAIC(:,110), 2.00000);
ALGEBRAIC(:,139) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,230)+CONSTANTS(:,231)))./power(ALGEBRAIC(:,132), 2.00000);
ALGEBRAIC(:,150) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,235)+CONSTANTS(:,236)))./power(ALGEBRAIC(:,143), 2.00000);
ALGEBRAIC(:,166) = ( ( - power(10.0000, - ALGEBRAIC(:,1))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,246)+CONSTANTS(:,243)))./power(ALGEBRAIC(:,160), 2.00000);
ALGEBRAIC(:,296) = ( ALGEBRAIC(:,9).*STATES(:,4)+ ALGEBRAIC(:,20).*STATES(:,11)+ ALGEBRAIC(:,32).*STATES(:,10)+ ALGEBRAIC(:,42).*STATES(:,6)+ ALGEBRAIC(:,53).*STATES(:,24)+ CONSTANTS(:,27).*STATES(:,23)+ ALGEBRAIC(:,67).*STATES(:,3)+ CONSTANTS(:,28).*STATES(:,7)+ CONSTANTS(:,29).*STATES(:,8)+ ALGEBRAIC(:,90).*STATES(:,9)+ CONSTANTS(:,30).*STATES(:,13)+ ALGEBRAIC(:,116).*STATES(:,12)+ ALGEBRAIC(:,106).*STATES(:,14)+ CONSTANTS(:,31).*STATES(:,17)+ CONSTANTS(:,32).*STATES(:,18)+ ALGEBRAIC(:,139).*STATES(:,19)+ ALGEBRAIC(:,150).*STATES(:,20)+ CONSTANTS(:,33).*STATES(:,21)+ ALGEBRAIC(:,166).*STATES(:,22))./(ALGEBRAIC(:,292)+ALGEBRAIC(:,293));
ALGEBRAIC(:,11) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,150)+CONSTANTS(:,151)).*log(10.0000))./power(ALGEBRAIC(:,2), 2.00000);
ALGEBRAIC(:,22) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,160)+CONSTANTS(:,161)).*log(10.0000))./power(ALGEBRAIC(:,13), 2.00000);
ALGEBRAIC(:,34) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,170)+CONSTANTS(:,171)).*log(10.0000))./power(ALGEBRAIC(:,24), 2.00000);
ALGEBRAIC(:,44) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,178)+CONSTANTS(:,179)).*log(10.0000))./power(ALGEBRAIC(:,36), 2.00000);
ALGEBRAIC(:,55) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,187)+CONSTANTS(:,188)).*log(10.0000))./power(ALGEBRAIC(:,46), 2.00000);
ALGEBRAIC(:,69) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,196)+CONSTANTS(:,197)).*log(10.0000))./power(ALGEBRAIC(:,61), 2.00000);
ALGEBRAIC(:,92) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,208)).*( ALGEBRAIC(:,83).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,206)).*log(10.0000) - power(10.0000, - 2.00000.*ALGEBRAIC(:,1)+CONSTANTS(:,206)+CONSTANTS(:,207)).*2.00000.*log(10.0000)))./power(ALGEBRAIC(:,83), 2.00000);
ALGEBRAIC(:,108) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,217)+CONSTANTS(:,214)).*log(10.0000))./power(ALGEBRAIC(:,100), 2.00000);
ALGEBRAIC(:,118) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,220)+CONSTANTS(:,221)).*log(10.0000))./power(ALGEBRAIC(:,110), 2.00000);
ALGEBRAIC(:,141) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,230)+CONSTANTS(:,231)).*log(10.0000))./power(ALGEBRAIC(:,132), 2.00000);
ALGEBRAIC(:,152) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,235)+CONSTANTS(:,236)).*log(10.0000))./power(ALGEBRAIC(:,143), 2.00000);
ALGEBRAIC(:,168) = ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, - ALGEBRAIC(:,1)+CONSTANTS(:,246)+CONSTANTS(:,243)).*log(10.0000))./power(ALGEBRAIC(:,160), 2.00000);
ALGEBRAIC(:,10) = (( ALGEBRAIC(:,2).*power(10.0000, CONSTANTS(:,151)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,151)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,2), 2.00000);
ALGEBRAIC(:,21) = (( ALGEBRAIC(:,13).*power(10.0000, CONSTANTS(:,161)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,161)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,13), 2.00000);
ALGEBRAIC(:,33) = (( ALGEBRAIC(:,24).*power(10.0000, CONSTANTS(:,171)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,171)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,24), 2.00000);
ALGEBRAIC(:,43) = (( ALGEBRAIC(:,36).*power(10.0000, CONSTANTS(:,179)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,179)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,36), 2.00000);
ALGEBRAIC(:,54) = (( ALGEBRAIC(:,46).*power(10.0000, CONSTANTS(:,188)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,188)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,46), 2.00000);
ALGEBRAIC(:,68) = (( ALGEBRAIC(:,61).*power(10.0000, CONSTANTS(:,197)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,197)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,61), 2.00000);
ALGEBRAIC(:,91) = (( ALGEBRAIC(:,83).*power(10.0000, CONSTANTS(:,208)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,208)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,83), 2.00000);
ALGEBRAIC(:,107) = (( ALGEBRAIC(:,100).*power(10.0000, CONSTANTS(:,214)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,214)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,100), 2.00000);
ALGEBRAIC(:,117) = (( ALGEBRAIC(:,110).*power(10.0000, CONSTANTS(:,221)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,221)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,110), 2.00000);
ALGEBRAIC(:,140) = (( ALGEBRAIC(:,132).*power(10.0000, CONSTANTS(:,231)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,231)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,132), 2.00000);
ALGEBRAIC(:,151) = (( ALGEBRAIC(:,143).*power(10.0000, CONSTANTS(:,236)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,236)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,143), 2.00000);
ALGEBRAIC(:,167) = (( ALGEBRAIC(:,160).*power(10.0000, CONSTANTS(:,243)))./CONSTANTS(:,135) - ( (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, 2.00000.*CONSTANTS(:,243)))./CONSTANTS(:,135))./power(ALGEBRAIC(:,160), 2.00000);
ALGEBRAIC(:,297) = - 1.00000 - ( ALGEBRAIC(:,21).*STATES(:,11)+ ALGEBRAIC(:,33).*STATES(:,10)+ ALGEBRAIC(:,43).*STATES(:,6)+ ALGEBRAIC(:,10).*STATES(:,4)+ ALGEBRAIC(:,54).*STATES(:,24)+ ALGEBRAIC(:,68).*STATES(:,3)+ ALGEBRAIC(:,91).*STATES(:,9)+ ALGEBRAIC(:,107).*STATES(:,14)+ ALGEBRAIC(:,117).*STATES(:,12)+ ALGEBRAIC(:,151).*STATES(:,20)+ ALGEBRAIC(:,140).*STATES(:,19)+ ALGEBRAIC(:,167).*STATES(:,22));
ALGEBRAIC(:,298) = ( ALGEBRAIC(:,22).*STATES(:,11)+ ALGEBRAIC(:,34).*STATES(:,10)+ ALGEBRAIC(:,44).*STATES(:,6)+ ALGEBRAIC(:,11).*STATES(:,4)+ ALGEBRAIC(:,55).*STATES(:,24)+ ALGEBRAIC(:,69).*STATES(:,3)+ ALGEBRAIC(:,92).*STATES(:,9)+ ALGEBRAIC(:,108).*STATES(:,14)+ ALGEBRAIC(:,118).*STATES(:,12)+ ALGEBRAIC(:,152).*STATES(:,20)+ ALGEBRAIC(:,141).*STATES(:,19)+ ALGEBRAIC(:,168).*STATES(:,22))./ALGEBRAIC(:,297);
ALGEBRAIC(:,6) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,151)).*ALGEBRAIC(:,3);
ALGEBRAIC(:,50) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,188)).*ALGEBRAIC(:,47);
ALGEBRAIC(:,64) = (( ALGEBRAIC(:,62).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,197));
ALGEBRAIC(:,87) = (( ALGEBRAIC(:,84).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,208));
ALGEBRAIC(:,103) = (( ALGEBRAIC(:,101).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,214));
ALGEBRAIC(:,113) = (( ALGEBRAIC(:,111).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,221));
ALGEBRAIC(:,136) = (( ALGEBRAIC(:,133).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,231));
ALGEBRAIC(:,147) = (( ALGEBRAIC(:,144).*STATES(:,2))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,236));
ALGEBRAIC(:,163) = (STATES(:,2)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,243)).*ALGEBRAIC(:,161);
ALGEBRAIC(:,312) = ( ALGEBRAIC(:,16).*ALGEBRAIC(:,301)+ ALGEBRAIC(:,27).*ALGEBRAIC(:,302)+ ALGEBRAIC(:,39).*ALGEBRAIC(:,303)+ ALGEBRAIC(:,6).*ALGEBRAIC(:,304)+ ALGEBRAIC(:,50).*ALGEBRAIC(:,300)+ ALGEBRAIC(:,64).*ALGEBRAIC(:,305)+ ALGEBRAIC(:,87).*ALGEBRAIC(:,306)+ ALGEBRAIC(:,103).*ALGEBRAIC(:,308)+ ALGEBRAIC(:,113).*ALGEBRAIC(:,307)+ ALGEBRAIC(:,147).*ALGEBRAIC(:,310)+ ALGEBRAIC(:,136).*ALGEBRAIC(:,309)+ ALGEBRAIC(:,163).*ALGEBRAIC(:,311))./ALGEBRAIC(:,297);
ALGEBRAIC(:,5) = (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,145)).*ALGEBRAIC(:,3);
ALGEBRAIC(:,31) = ((ALGEBRAIC(:,29)+ALGEBRAIC(:,7)) - ALGEBRAIC(:,18))+(((CONSTANTS(:,172)+CONSTANTS(:,152)) - CONSTANTS(:,163)) - CONSTANTS(:,254));
ALGEBRAIC(:,48) = (CONSTANTS(:,134)./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,182)).*ALGEBRAIC(:,47);
ALGEBRAIC(:,135) = (( ALGEBRAIC(:,133).*CONSTANTS(:,134))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,232));
ALGEBRAIC(:,146) = (( ALGEBRAIC(:,144).*CONSTANTS(:,134))./CONSTANTS(:,135)).*power(10.0000, CONSTANTS(:,237));
ALGEBRAIC(:,176) = ((ALGEBRAIC(:,35)+ALGEBRAIC(:,12)+ALGEBRAIC(:,174)) - ALGEBRAIC(:,173)) - ALGEBRAIC(:,23);
ALGEBRAIC(:,178) = ( exp( - ALGEBRAIC(:,176)./( CONSTANTS(:,1).*CONSTANTS(:,131))).*ALGEBRAIC(:,24).*ALGEBRAIC(:,2))./( ALGEBRAIC(:,13).*power(10.0000, - ALGEBRAIC(:,1)));
ALGEBRAIC(:,281) = ALGEBRAIC(:,185).*ALGEBRAIC(:,239)+ ALGEBRAIC(:,182).*(ALGEBRAIC(:,227)+ALGEBRAIC(:,235))+ ALGEBRAIC(:,188).*ALGEBRAIC(:,242)+ ALGEBRAIC(:,191).*ALGEBRAIC(:,249)+ ALGEBRAIC(:,194).*ALGEBRAIC(:,253)+ ALGEBRAIC(:,197).*ALGEBRAIC(:,255)+ ALGEBRAIC(:,200).*ALGEBRAIC(:,262)+ ALGEBRAIC(:,206).*ALGEBRAIC(:,266)+ ALGEBRAIC(:,209).*ALGEBRAIC(:,270)+ ALGEBRAIC(:,212).*ALGEBRAIC(:,273)+ ALGEBRAIC(:,215).*ALGEBRAIC(:,276)+ ALGEBRAIC(:,218).*ALGEBRAIC(:,279)+ ALGEBRAIC(:,203).*ALGEBRAIC(:,258);
ALGEBRAIC(:,289) = ALGEBRAIC(:,175).* - ALGEBRAIC(:,287);
ALGEBRAIC(:,299) = 1.00000 - ALGEBRAIC(:,298).*ALGEBRAIC(:,296);
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
% Compute a logarithm to any base" +
function x = arbitrary_log(a, base)
x = log(a) ./ log(base);
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