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 =296;
end
% There are a total of 80 entries in each of the rate and state variable arrays.
% There are a total of 770 entries in the constant variable array.
%
function [VOI, STATES, ALGEBRAIC, CONSTANTS] = solveModel()
% Create ALGEBRAIC of correct size
global algebraicVariableCount; algebraicVariableCount = getAlgebraicVariableCount();
% Initialise constants and state variables
[INIT_STATES, CONSTANTS] = initConsts;
% Set timespan to solve over
tspan = [0, 10];
% Set numerical accuracy options for ODE solver
options = odeset('RelTol', 1e-06, 'AbsTol', 1e-06, 'MaxStep', 1);
% Solve model with ODE solver
[VOI, STATES] = ode15s(@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS), tspan, INIT_STATES, options);
% Compute algebraic variables
[RATES, ALGEBRAIC] = computeRates(VOI, STATES, CONSTANTS);
ALGEBRAIC = computeAlgebraic(ALGEBRAIC, CONSTANTS, STATES, VOI);
% Plot state variables against variable of integration
[LEGEND_STATES, LEGEND_ALGEBRAIC, LEGEND_VOI, LEGEND_CONSTANTS] = createLegends();
figure();
plot(VOI, STATES);
xlabel(LEGEND_VOI);
l = legend(LEGEND_STATES);
set(l,'Interpreter','none');
end
function [LEGEND_STATES, LEGEND_ALGEBRAIC, LEGEND_VOI, LEGEND_CONSTANTS] = createLegends()
LEGEND_STATES = ''; LEGEND_ALGEBRAIC = ''; LEGEND_VOI = ''; LEGEND_CONSTANTS = '';
LEGEND_VOI = strpad('time in component environment (second)');
LEGEND_ALGEBRAIC(:,8) = strpad('I_CaT in component I_CaT (picoampere)');
LEGEND_CONSTANTS(:,1) = strpad('g_CaT in component I_CaT (picosiemens)');
LEGEND_CONSTANTS(:,2) = strpad('ECaT in component I_CaT (millivolt)');
LEGEND_CONSTANTS(:,3) = strpad('alpha_OT in component I_CaT (first_order_rate_constant)');
LEGEND_CONSTANTS(:,4) = strpad('beta_OT in component I_CaT (first_order_rate_constant)');
LEGEND_STATES(:,1) = strpad('OCaT in component I_CaT (dimensionless)');
LEGEND_STATES(:,2) = strpad('Vm in component membrane (millivolt)');
LEGEND_STATES(:,3) = strpad('dT in component I_CaT_d_gate (dimensionless)');
LEGEND_STATES(:,4) = strpad('fT in component I_CaT_f_gate (dimensionless)');
LEGEND_ALGEBRAIC(:,1) = strpad('infinity_dT in component I_CaT_d_gate (dimensionless)');
LEGEND_CONSTANTS(:,688) = strpad('tau_dT1 in component I_CaT_d_gate (second)');
LEGEND_CONSTANTS(:,5) = strpad('A_dT1 in component I_CaT_d_gate (second)');
LEGEND_CONSTANTS(:,6) = strpad('k_dT in component I_CaT_d_gate (per_millivolt)');
LEGEND_CONSTANTS(:,7) = strpad('V_dT in component I_CaT_d_gate (millivolt)');
LEGEND_ALGEBRAIC(:,2) = strpad('infinity_fT in component I_CaT_f_gate (dimensionless)');
LEGEND_ALGEBRAIC(:,20) = strpad('tau_fT1 in component I_CaT_f_gate (second)');
LEGEND_CONSTANTS(:,8) = strpad('A_fT1 in component I_CaT_f_gate (second)');
LEGEND_CONSTANTS(:,9) = strpad('A_fT2 in component I_CaT_f_gate (second)');
LEGEND_CONSTANTS(:,10) = strpad('A_fT3 in component I_CaT_f_gate (per_millivolt)');
LEGEND_CONSTANTS(:,11) = strpad('A_fT4 in component I_CaT_f_gate (millivolt)');
LEGEND_CONSTANTS(:,12) = strpad('k_fT in component I_CaT_f_gate (per_millivolt)');
LEGEND_CONSTANTS(:,13) = strpad('V_fT in component I_CaT_f_gate (millivolt)');
LEGEND_ALGEBRAIC(:,26) = strpad('I_CaExt in component I_CaExt (picoampere)');
LEGEND_CONSTANTS(:,14) = strpad('gCaExt in component I_CaExt (femtoampere)');
LEGEND_CONSTANTS(:,15) = strpad('KCaExt in component I_CaExt (micromolar)');
LEGEND_STATES(:,5) = strpad('CCy in component CCy (micromolar)');
LEGEND_ALGEBRAIC(:,30) = strpad('I_Kv1_1 in component I_Kv1_1 (picoampere)');
LEGEND_CONSTANTS(:,16) = strpad('g_Kv1_1 in component I_Kv1_1 (picosiemens)');
LEGEND_CONSTANTS(:,689) = strpad('EK in component Nerst_potentials (millivolt)');
LEGEND_STATES(:,6) = strpad('dv1_1 in component I_Kv1_1_d_gate (dimensionless)');
LEGEND_STATES(:,7) = strpad('fv1_1 in component I_Kv1_1_f_gate (dimensionless)');
LEGEND_ALGEBRAIC(:,21) = strpad('alpha_dv1_1 in component I_Kv1_1_d_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,3) = strpad('A_dv1_11 in component I_Kv1_1_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,17) = strpad('A_dv1_11a in component I_Kv1_1_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,18) = strpad('A_dv1_11b in component I_Kv1_1_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,19) = strpad('A_dv1_11c in component I_Kv1_1_d_gate (per_millivolt)');
LEGEND_CONSTANTS(:,20) = strpad('A_dv1_11d in component I_Kv1_1_d_gate (millivolt)');
LEGEND_CONSTANTS(:,21) = strpad('A_dv1_12 in component I_Kv1_1_d_gate (dimensionless)');
LEGEND_CONSTANTS(:,22) = strpad('A_dv1_13 in component I_Kv1_1_d_gate (per_millivolt)');
LEGEND_CONSTANTS(:,23) = strpad('A_dv1_14 in component I_Kv1_1_d_gate (millivolt)');
LEGEND_ALGEBRAIC(:,27) = strpad('beta_dv1_1 in component I_Kv1_1_d_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,4) = strpad('alpha_fv1_1 in component I_Kv1_1_f_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,24) = strpad('A_fv1_11 in component I_Kv1_1_f_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,25) = strpad('A_fv1_12 in component I_Kv1_1_f_gate (dimensionless)');
LEGEND_CONSTANTS(:,26) = strpad('A_fv1_13 in component I_Kv1_1_f_gate (per_millivolt)');
LEGEND_CONSTANTS(:,27) = strpad('A_fv1_14 in component I_Kv1_1_f_gate (millivolt)');
LEGEND_ALGEBRAIC(:,22) = strpad('beta_fv1_1 in component I_Kv1_1_f_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,31) = strpad('I_KERG in component I_KERG (picoampere)');
LEGEND_CONSTANTS(:,28) = strpad('g_KERG in component I_KERG (picosiemens)');
LEGEND_STATES(:,8) = strpad('dERG in component I_KERG_d_gate (dimensionless)');
LEGEND_ALGEBRAIC(:,5) = strpad('alpha_dERG in component I_KERG_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,29) = strpad('A_dERG1 in component I_KERG_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,30) = strpad('A_dERG2 in component I_KERG_d_gate (dimensionless)');
LEGEND_CONSTANTS(:,31) = strpad('A_dERG3 in component I_KERG_d_gate (per_millivolt)');
LEGEND_CONSTANTS(:,32) = strpad('A_dERG4 in component I_KERG_d_gate (millivolt)');
LEGEND_ALGEBRAIC(:,23) = strpad('beta_dERG in component I_KERG_d_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,32) = strpad('I_KB in component I_KB (picoampere)');
LEGEND_CONSTANTS(:,33) = strpad('g_KB in component I_KB (picosiemens)');
LEGEND_ALGEBRAIC(:,33) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_CONSTANTS(:,34) = strpad('g_Na in component I_Na (picosiemens)');
LEGEND_CONSTANTS(:,690) = strpad('ENa_Cy in component Nerst_potentials (millivolt)');
LEGEND_STATES(:,9) = strpad('dNa in component I_Na_d_gate (dimensionless)');
LEGEND_STATES(:,10) = strpad('fNa in component I_Na_f_gate (dimensionless)');
LEGEND_ALGEBRAIC(:,24) = strpad('alpha_dNa in component I_Na_d_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,6) = strpad('A_dNa1 in component I_Na_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,35) = strpad('A_dNa1a in component I_Na_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,36) = strpad('A_dNa1b in component I_Na_d_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,37) = strpad('A_dNa1c in component I_Na_d_gate (per_millivolt)');
LEGEND_CONSTANTS(:,38) = strpad('A_dNa1d in component I_Na_d_gate (millivolt)');
LEGEND_CONSTANTS(:,39) = strpad('A_dNa2 in component I_Na_d_gate (dimensionless)');
LEGEND_CONSTANTS(:,40) = strpad('A_dNa3 in component I_Na_d_gate (per_millivolt)');
LEGEND_CONSTANTS(:,41) = strpad('A_dNa4 in component I_Na_d_gate (millivolt)');
LEGEND_ALGEBRAIC(:,28) = strpad('beta_dNa in component I_Na_d_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,25) = strpad('alpha_fNa in component I_Na_f_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,7) = strpad('A_fNa1 in component I_Na_f_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,42) = strpad('A_fNa1a in component I_Na_f_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,43) = strpad('A_fNa1b in component I_Na_f_gate (first_order_rate_constant)');
LEGEND_CONSTANTS(:,44) = strpad('A_fNa1c in component I_Na_f_gate (per_millivolt)');
LEGEND_CONSTANTS(:,45) = strpad('A_fNa1d in component I_Na_f_gate (millivolt)');
LEGEND_CONSTANTS(:,46) = strpad('A_fNa2 in component I_Na_f_gate (dimensionless)');
LEGEND_CONSTANTS(:,47) = strpad('A_fNa3 in component I_Na_f_gate (per_millivolt)');
LEGEND_CONSTANTS(:,48) = strpad('A_fNa4 in component I_Na_f_gate (millivolt)');
LEGEND_ALGEBRAIC(:,29) = strpad('beta_fNa in component I_Na_f_gate (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,34) = strpad('I_L in component I_L (picoampere)');
LEGEND_CONSTANTS(:,49) = strpad('g_L in component I_L (picosiemens)');
LEGEND_CONSTANTS(:,50) = strpad('EL in component I_L (millivolt)');
LEGEND_ALGEBRAIC(:,35) = strpad('I_ion_Cy in component membrane (picoampere)');
LEGEND_CONSTANTS(:,51) = strpad('Cm in component membrane (picofarad)');
LEGEND_ALGEBRAIC(:,73) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,97) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,121) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,145) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,169) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,193) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,217) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,241) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,265) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,288) = strpad('I_ion_PU in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,9) = strpad('ECa_Cy in component Nerst_potentials (millivolt)');
LEGEND_CONSTANTS(:,52) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,53) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,54) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,55) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,56) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,57) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,58) = strpad('NO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,59) = strpad('Ni in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,60) = strpad('lambda_S2_Cy in component CCy (dimensionless)');
LEGEND_CONSTANTS(:,61) = strpad('delta_SCy in component CCy (micromolar_per_coulomb)');
LEGEND_ALGEBRAIC(:,46) = strpad('JCy in component JCy (flux)');
LEGEND_ALGEBRAIC(:,36) = strpad('JS2Cy1 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,37) = strpad('JS2Cy2 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,38) = strpad('JS2Cy3 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,39) = strpad('JS2Cy4 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,40) = strpad('JS2Cy5 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,41) = strpad('JS2Cy6 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,42) = strpad('JS2Cy7 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,43) = strpad('JS2Cy8 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,44) = strpad('JS2Cy9 in component JS2Cy (flux)');
LEGEND_ALGEBRAIC(:,45) = strpad('JS2Cy10 in component JS2Cy (flux)');
LEGEND_CONSTANTS(:,62) = strpad('mu_S2Cy1 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,63) = strpad('mu_S2Cy2 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,64) = strpad('mu_S2Cy3 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,65) = strpad('mu_S2Cy4 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,66) = strpad('mu_S2Cy5 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,67) = strpad('mu_S2Cy6 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,68) = strpad('mu_S2Cy7 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,69) = strpad('mu_S2Cy8 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,70) = strpad('mu_S2Cy9 in component JS2Cy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,71) = strpad('mu_S2Cy10 in component JS2Cy (first_order_rate_constant)');
LEGEND_STATES(:,11) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,12) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,13) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,14) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,15) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,16) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,17) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,18) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,19) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_STATES(:,20) = strpad('CS2 in component CS2 (micromolar)');
LEGEND_ALGEBRAIC(:,47) = strpad('JS1S2_1 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,48) = strpad('JS1S2_2 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,49) = strpad('JS1S2_3 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,50) = strpad('JS1S2_4 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,51) = strpad('JS1S2_5 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,52) = strpad('JS1S2_6 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,53) = strpad('JS1S2_7 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,54) = strpad('JS1S2_8 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,55) = strpad('JS1S2_9 in component JS1S2 (flux)');
LEGEND_ALGEBRAIC(:,56) = strpad('JS1S2_10 in component JS1S2 (flux)');
LEGEND_CONSTANTS(:,692) = strpad('mu_S1S2_1 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,693) = strpad('mu_S1S2_2 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,694) = strpad('mu_S1S2_3 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,695) = strpad('mu_S1S2_4 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,696) = strpad('mu_S1S2_5 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,697) = strpad('mu_S1S2_6 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,698) = strpad('mu_S1S2_7 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,699) = strpad('mu_S1S2_8 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,700) = strpad('mu_S1S2_9 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,701) = strpad('mu_S1S2_10 in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,72) = strpad('mu_A in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,73) = strpad('mu_B in component JS1S2 (first_order_rate_constant)');
LEGEND_CONSTANTS(:,74) = strpad('n_PU in component model_parameters (dimensionless)');
LEGEND_STATES(:,21) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,22) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,23) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,24) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,25) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,26) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,27) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,28) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,29) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_STATES(:,30) = strpad('CS1 in component CS1 (micromolar)');
LEGEND_CONSTANTS(:,75) = strpad('mu_Cy in component JCy (first_order_rate_constant)');
LEGEND_CONSTANTS(:,76) = strpad('C_infinity in component JCy (micromolar)');
LEGEND_CONSTANTS(:,77) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,67) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,71) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,716) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,59) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,61) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,69) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,63) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,66) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,65) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,78) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,79) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,80) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,81) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,57) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,58) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,82) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,83) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,84) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,85) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,86) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,60) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,87) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,31) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,88) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,62) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,89) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,90) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,91) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,92) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,64) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,93) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,94) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,95) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,96) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,97) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,98) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,68) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,99) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,100) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,101) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,102) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,103) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,32) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,72) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,104) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,105) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,70) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,106) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,107) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,33) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,80) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,108) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,109) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,110) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,111) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,112) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,113) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,114) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,115) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,116) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,74) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,77) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,81) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,83) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,34) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,117) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,118) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,119) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,120) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,702) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,10) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,35) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,121) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,122) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,123) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,124) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,125) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,126) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,75) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,127) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,128) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,691) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,712) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,129) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,130) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,131) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,132) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,713) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,714) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,715) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,133) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,78) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,134) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,135) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,136) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,137) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,138) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,91) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,95) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,722) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,82) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,85) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,93) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,87) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,90) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,89) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,139) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,140) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,141) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,142) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,76) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,79) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,143) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,144) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,145) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,146) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,147) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,84) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,148) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,36) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,149) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,86) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,150) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,151) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,152) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,153) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,88) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,154) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,155) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,156) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,157) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,158) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,159) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,92) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,160) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,161) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,162) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,163) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,164) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,37) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,96) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,165) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,166) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,94) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,167) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,168) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,38) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,104) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,169) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,170) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,171) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,172) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,173) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,174) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,175) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,176) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,177) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,98) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,101) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,105) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,107) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,39) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,178) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,179) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,180) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,181) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,703) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,11) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,40) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,182) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,183) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,184) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,185) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,186) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,187) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,99) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,188) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,189) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,717) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,718) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,190) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,191) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,192) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,193) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,719) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,720) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,721) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,194) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,102) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,195) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,196) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,197) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,198) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,199) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,115) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,119) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,728) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,106) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,109) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,117) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,111) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,114) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,113) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,200) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,201) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,202) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,203) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,100) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,103) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,204) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,205) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,206) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,207) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,208) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,108) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,209) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,41) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,210) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,110) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,211) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,212) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,213) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,214) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,112) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,215) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,216) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,217) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,218) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,219) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,220) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,116) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,221) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,222) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,223) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,224) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,225) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,42) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,120) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,226) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,227) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,118) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,228) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,229) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,43) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,128) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,230) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,231) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,232) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,233) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,234) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,235) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,236) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,237) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,238) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,122) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,125) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,129) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,131) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,44) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,239) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,240) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,241) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,242) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,704) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,12) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,45) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,243) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,244) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,245) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,246) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,247) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,248) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,123) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,249) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,250) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,723) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,724) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,251) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,252) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,253) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,254) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,725) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,726) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,727) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,255) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,126) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,256) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,257) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,258) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,259) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,260) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,139) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,143) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,734) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,130) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,133) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,141) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,135) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,138) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,137) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,261) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,262) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,263) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,264) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,124) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,127) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,265) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,266) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,267) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,268) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,269) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,132) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,270) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,46) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,271) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,134) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,272) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,273) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,274) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,275) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,136) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,276) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,277) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,278) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,279) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,280) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,281) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,140) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,282) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,283) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,284) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,285) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,286) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,47) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,144) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,287) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,288) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,142) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,289) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,290) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,48) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,152) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,291) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,292) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,293) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,294) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,295) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,296) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,297) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,298) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,299) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,146) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,149) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,153) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,155) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,49) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,300) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,301) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,302) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,303) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,705) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,13) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,50) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,304) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,305) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,306) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,307) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,308) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,309) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,147) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,310) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,311) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,729) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,730) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,312) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,313) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,314) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,315) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,731) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,732) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,733) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,316) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,150) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,317) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,318) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,319) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,320) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,321) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,163) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,167) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,740) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,154) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,157) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,165) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,159) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,162) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,161) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,322) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,323) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,324) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,325) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,148) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,151) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,326) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,327) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,328) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,329) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,330) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,156) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,331) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,51) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,332) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,158) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,333) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,334) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,335) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,336) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,160) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,337) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,338) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,339) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,340) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,341) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,342) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,164) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,343) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,344) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,345) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,346) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,347) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,52) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,168) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,348) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,349) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,166) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,350) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,351) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,53) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,176) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,352) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,353) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,354) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,355) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,356) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,357) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,358) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,359) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,360) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,170) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,173) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,177) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,179) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,54) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,361) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,362) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,363) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,364) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,706) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,14) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,55) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,365) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,366) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,367) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,368) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,369) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,370) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,171) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,371) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,372) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,735) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,736) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,373) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,374) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,375) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,376) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,737) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,738) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,739) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,377) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,174) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,378) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,379) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,380) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,381) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,382) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,187) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,191) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,746) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,178) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,181) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,189) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,183) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,186) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,185) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,383) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,384) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,385) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,386) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,172) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,175) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,387) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,388) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,389) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,390) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,391) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,180) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,392) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,56) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,393) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,182) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,394) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,395) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,396) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,397) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,184) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,398) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,399) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,400) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,401) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,402) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,403) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,188) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,404) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,405) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,406) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,407) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,408) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,57) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,192) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,409) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,410) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,190) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,411) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,412) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,58) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,200) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,413) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,414) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,415) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,416) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,417) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,418) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,419) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,420) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,421) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,194) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,197) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,201) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,203) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,59) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,422) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,423) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,424) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,425) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,707) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,15) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,60) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,426) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,427) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,428) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,429) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,430) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,431) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,195) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,432) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,433) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,741) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,742) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,434) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,435) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,436) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,437) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,743) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,744) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,745) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,438) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,198) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,439) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,440) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,441) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,442) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,443) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,211) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,215) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,752) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,202) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,205) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,213) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,207) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,210) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,209) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,444) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,445) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,446) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,447) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,196) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,199) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,448) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,449) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,450) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,451) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,452) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,204) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,453) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,61) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,454) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,206) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,455) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,456) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,457) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,458) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,208) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,459) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,460) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,461) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,462) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,463) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,464) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,212) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,465) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,466) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,467) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,468) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,469) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,62) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,216) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,470) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,471) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,214) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,472) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,473) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,63) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,224) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,474) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,475) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,476) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,477) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,478) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,479) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,480) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,481) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,482) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,218) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,221) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,225) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,227) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,64) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,483) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,484) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,485) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,486) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,708) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,16) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,65) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,487) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,488) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,489) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,490) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,491) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,492) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,219) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,493) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,494) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,747) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,748) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,495) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,496) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,497) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,498) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,749) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,750) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,751) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,499) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,222) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,500) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,501) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,502) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,503) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,504) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,235) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,239) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,758) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,226) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,229) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,237) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,231) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,234) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,233) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,505) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,506) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,507) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,508) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,220) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,223) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,509) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,510) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,511) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,512) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,513) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,228) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,514) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,66) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,515) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,230) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,516) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,517) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,518) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,519) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,232) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,520) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,521) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,522) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,523) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,524) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,525) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,236) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,526) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,527) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,528) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,529) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,530) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,67) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,240) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,531) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,532) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,238) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,533) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,534) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,68) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,248) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,535) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,536) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,537) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,538) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,539) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,540) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,541) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,542) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,543) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,242) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,245) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,249) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,251) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,69) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,544) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,545) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,546) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,547) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,709) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,17) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,70) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,548) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,549) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,550) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,551) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,552) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,553) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,243) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,554) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,555) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,753) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,754) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,556) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,557) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,558) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,559) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,755) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,756) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,757) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,560) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,246) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,561) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,562) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,563) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,564) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,565) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,259) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,263) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,764) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,250) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,253) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,261) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,255) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,258) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,257) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,566) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,567) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,568) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,569) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,244) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,247) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,570) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,571) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,572) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,573) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,574) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,252) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,575) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,71) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,576) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,254) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,577) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,578) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,579) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,580) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,256) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,581) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,582) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,583) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,584) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,585) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,586) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,260) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,587) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,588) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,589) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,590) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,591) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,72) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,264) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,592) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,593) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,262) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,594) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,595) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,73) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,272) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,596) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,597) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,598) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,599) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,600) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,601) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,602) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,603) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,604) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,266) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,269) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,273) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,275) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,74) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,605) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,606) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,607) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,608) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,710) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,18) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,75) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,609) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,610) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,611) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,612) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,613) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,614) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,267) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,615) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,616) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,759) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,760) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,617) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,618) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,619) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,620) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,761) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,762) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,763) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,621) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,270) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,622) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,623) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,624) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,625) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,626) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,283) = strpad('I_iCa in component PU_membrane (picoampere)');
LEGEND_ALGEBRAIC(:,286) = strpad('I_iNa in component PU_membrane (picoampere)');
LEGEND_CONSTANTS(:,770) = strpad('alpha_scale in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,274) = strpad('I_Ca in component I_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,277) = strpad('I_Na in component I_Na (picoampere)');
LEGEND_ALGEBRAIC(:,284) = strpad('I_NaP in component I_NaP (picoampere)');
LEGEND_ALGEBRAIC(:,279) = strpad('I_NSCC_Ca in component I_NSCC_Ca (picoampere)');
LEGEND_ALGEBRAIC(:,282) = strpad('I_PM in component I_PM (picoampere)');
LEGEND_ALGEBRAIC(:,281) = strpad('I_NSCC_Na in component I_NSCC_Na (picoampere)');
LEGEND_CONSTANTS(:,627) = strpad('gCa_ in component I_Ca (picosiemens)');
LEGEND_CONSTANTS(:,628) = strpad('kCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,629) = strpad('kVCa in component I_Ca (per_millivolt)');
LEGEND_CONSTANTS(:,630) = strpad('VhCa in component I_Ca (millivolt)');
LEGEND_ALGEBRAIC(:,268) = strpad('gCa in component I_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,271) = strpad('ECa_PU in component I_Ca (millivolt)');
LEGEND_CONSTANTS(:,631) = strpad('T in component model_parameters (kelvin)');
LEGEND_CONSTANTS(:,632) = strpad('R in component model_parameters (attojoule_per_zeptomole_kelvin)');
LEGEND_CONSTANTS(:,633) = strpad('F in component model_parameters (femtocoulomb_per_zeptomole)');
LEGEND_CONSTANTS(:,634) = strpad('CO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,635) = strpad('gNa in component I_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,276) = strpad('ENa_PU in component I_Na (millivolt)');
LEGEND_CONSTANTS(:,636) = strpad('NO in component model_parameters (micromolar)');
LEGEND_STATES(:,76) = strpad('NS1 in component NS1 (micromolar)');
LEGEND_CONSTANTS(:,637) = strpad('gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)');
LEGEND_ALGEBRAIC(:,278) = strpad('gNSCC_Ca in component I_NSCC_Ca (picosiemens)');
LEGEND_CONSTANTS(:,638) = strpad('hNSCC in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,639) = strpad('ENSCC in component model_parameters (millivolt)');
LEGEND_CONSTANTS(:,640) = strpad('KNSCC in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,641) = strpad('gNSCC_Na_ in component I_NSCC_Na (picosiemens)');
LEGEND_ALGEBRAIC(:,280) = strpad('gNSCC_Na in component I_NSCC_Na (picosiemens)');
LEGEND_CONSTANTS(:,642) = strpad('gPM in component I_PM (femtoampere)');
LEGEND_CONSTANTS(:,643) = strpad('KPM in component I_PM (micromolar)');
LEGEND_CONSTANTS(:,644) = strpad('gNaP in component I_NaP (picosiemens)');
LEGEND_CONSTANTS(:,645) = strpad('KNaP in component I_NaP (micromolar)');
LEGEND_CONSTANTS(:,646) = strpad('hNaP in component I_NaP (dimensionless)');
LEGEND_CONSTANTS(:,647) = strpad('ENaP in component I_NaP (millivolt)');
LEGEND_ALGEBRAIC(:,285) = strpad('JSERCA in component JSERCA (flux)');
LEGEND_CONSTANTS(:,648) = strpad('VSERCA in component JSERCA (first_order_rate_constant)');
LEGEND_CONSTANTS(:,649) = strpad('A2 in component JSERCA (dimensionless)');
LEGEND_CONSTANTS(:,650) = strpad('A4 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,651) = strpad('A5 in component JSERCA (per_micromolar)');
LEGEND_CONSTANTS(:,652) = strpad('A6 in component JSERCA (per_micromolar2)');
LEGEND_STATES(:,77) = strpad('CER in component CER (micromolar)');
LEGEND_ALGEBRAIC(:,289) = strpad('JMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,653) = strpad('VMCU in component JMCU (flux)');
LEGEND_CONSTANTS(:,654) = strpad('KMCU in component JMCU (micromolar)');
LEGEND_ALGEBRAIC(:,287) = strpad('epsilon_INH in component JMCU (dimensionless)');
LEGEND_CONSTANTS(:,655) = strpad('KINH in component JMCU (micromolar)');
LEGEND_CONSTANTS(:,656) = strpad('hINH in component JMCU (dimensionless)');
LEGEND_STATES(:,78) = strpad('CMT in component CMT (micromolar)');
LEGEND_ALGEBRAIC(:,294) = strpad('JIPR in component JIPR (flux)');
LEGEND_CONSTANTS(:,657) = strpad('kIPR in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,658) = strpad('k_1 in component JIPR (flux)');
LEGEND_CONSTANTS(:,659) = strpad('k1 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,660) = strpad('k2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,661) = strpad('r2 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,662) = strpad('r_2 in component JIPR (flux)');
LEGEND_CONSTANTS(:,663) = strpad('r4 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,664) = strpad('R1 in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,665) = strpad('R3 in component JIPR (micromolar)');
LEGEND_ALGEBRAIC(:,290) = strpad('phi1 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,292) = strpad('phi_1 in component JIPR (flux)');
LEGEND_ALGEBRAIC(:,295) = strpad('phi2 in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,296) = strpad('phi3 in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,79) = strpad('H in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,666) = strpad('g_beta in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,667) = strpad('h_beta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,668) = strpad('g_alpha in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,669) = strpad('K_beta in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,711) = strpad('alpha_zeta in component JIPR (first_order_rate_constant)');
LEGEND_ALGEBRAIC(:,19) = strpad('beta_zeta in component JIPR (first_order_rate_constant)');
LEGEND_STATES(:,80) = strpad('zeta in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,670) = strpad('K_phi3_act in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,671) = strpad('K_phi3_inh in component JIPR (micromolar)');
LEGEND_CONSTANTS(:,672) = strpad('h_phi3_act in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,673) = strpad('h_phi3_inh in component JIPR (dimensionless)');
LEGEND_CONSTANTS(:,674) = strpad('g_phi3 in component JIPR (first_order_rate_constant)');
LEGEND_CONSTANTS(:,675) = strpad('P in component model_parameters (micromolar)');
LEGEND_ALGEBRAIC(:,291) = strpad('JNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,676) = strpad('VNCX in component JNCX (flux)');
LEGEND_CONSTANTS(:,677) = strpad('KNCX in component JNCX (micromolar)');
LEGEND_CONSTANTS(:,765) = strpad('lambda_MT_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,766) = strpad('lambda_ER_S1 in component CS1 (dimensionless)');
LEGEND_CONSTANTS(:,678) = strpad('delta_SPU in component model_parameters (micromolar_per_coulomb)');
LEGEND_CONSTANTS(:,679) = strpad('gamma_S1 in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,680) = strpad('gamma_MT in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,681) = strpad('gamma_ER in component model_parameters (dimensionless)');
LEGEND_CONSTANTS(:,767) = strpad('lambda_MT_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,768) = strpad('lambda_ER_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,769) = strpad('lambda_S1_S2 in component CS2 (dimensionless)');
LEGEND_CONSTANTS(:,682) = strpad('gamma_S2 in component model_parameters (dimensionless)');
LEGEND_ALGEBRAIC(:,293) = strpad('fm in component CMT (dimensionless)');
LEGEND_CONSTANTS(:,683) = strpad('Km in component CMT (micromolar)');
LEGEND_CONSTANTS(:,684) = strpad('Bm in component CMT (micromolar)');
LEGEND_CONSTANTS(:,685) = strpad('KO in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,686) = strpad('Ki in component model_parameters (micromolar)');
LEGEND_CONSTANTS(:,687) = strpad('n_PU_base in component model_parameters (dimensionless)');
LEGEND_RATES(:,1) = strpad('d/dt OCaT in component I_CaT (dimensionless)');
LEGEND_RATES(:,3) = strpad('d/dt dT in component I_CaT_d_gate (dimensionless)');
LEGEND_RATES(:,4) = strpad('d/dt fT in component I_CaT_f_gate (dimensionless)');
LEGEND_RATES(:,6) = strpad('d/dt dv1_1 in component I_Kv1_1_d_gate (dimensionless)');
LEGEND_RATES(:,7) = strpad('d/dt fv1_1 in component I_Kv1_1_f_gate (dimensionless)');
LEGEND_RATES(:,8) = strpad('d/dt dERG in component I_KERG_d_gate (dimensionless)');
LEGEND_RATES(:,9) = strpad('d/dt dNa in component I_Na_d_gate (dimensionless)');
LEGEND_RATES(:,10) = strpad('d/dt fNa in component I_Na_f_gate (dimensionless)');
LEGEND_RATES(:,2) = strpad('d/dt Vm in component membrane (millivolt)');
LEGEND_RATES(:,5) = strpad('d/dt CCy in component CCy (micromolar)');
LEGEND_RATES(:,34) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,35) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,21) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,11) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,32) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,33) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,31) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,39) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,40) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,22) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,12) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,37) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,38) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,36) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,44) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,45) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,23) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,13) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,42) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,43) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,41) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,49) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,50) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,24) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,14) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,47) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,48) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,46) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,54) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,55) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,25) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,15) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,52) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,53) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,51) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,59) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,60) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,26) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,16) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,57) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,58) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,56) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,64) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,65) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,27) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,17) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,62) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,63) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,61) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,69) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,70) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,28) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,18) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,67) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,68) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,66) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,74) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,75) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,29) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,19) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,72) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,73) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,71) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_RATES(:,79) = strpad('d/dt H in component JIPR (dimensionless)');
LEGEND_RATES(:,80) = strpad('d/dt zeta in component JIPR (dimensionless)');
LEGEND_RATES(:,30) = strpad('d/dt CS1 in component CS1 (micromolar)');
LEGEND_RATES(:,20) = strpad('d/dt CS2 in component CS2 (micromolar)');
LEGEND_RATES(:,77) = strpad('d/dt CER in component CER (micromolar)');
LEGEND_RATES(:,78) = strpad('d/dt CMT in component CMT (micromolar)');
LEGEND_RATES(:,76) = strpad('d/dt NS1 in component NS1 (micromolar)');
LEGEND_STATES = LEGEND_STATES';
LEGEND_ALGEBRAIC = LEGEND_ALGEBRAIC';
LEGEND_RATES = LEGEND_RATES';
LEGEND_CONSTANTS = LEGEND_CONSTANTS';
end
function [STATES, CONSTANTS] = initConsts()
VOI = 0; CONSTANTS = []; STATES = []; ALGEBRAIC = [];
CONSTANTS(:,1) = 700.0;
CONSTANTS(:,2) = 17.0;
CONSTANTS(:,3) = 800.0;
CONSTANTS(:,4) = 72;
STATES(:,1) = 0;
STATES(:,2) = -70.0;
STATES(:,3) = 0.010;
STATES(:,4) = 0.001;
CONSTANTS(:,5) = 0.0025;
CONSTANTS(:,6) = -0.60;
CONSTANTS(:,7) = -53;
CONSTANTS(:,8) = 0.019;
CONSTANTS(:,9) = 6.75;
CONSTANTS(:,10) = 2;
CONSTANTS(:,11) = -40;
CONSTANTS(:,12) = 1;
CONSTANTS(:,13) = -65;
CONSTANTS(:,14) = 100.0;
CONSTANTS(:,15) = 1.0;
STATES(:,5) = 0.12;
CONSTANTS(:,16) = 10.0;
STATES(:,6) = 0.000;
STATES(:,7) = 1.0;
CONSTANTS(:,17) = 0.0;
CONSTANTS(:,18) = 250.0;
CONSTANTS(:,19) = -0.077;
CONSTANTS(:,20) = -4.8;
CONSTANTS(:,21) = 0.80;
CONSTANTS(:,22) = -0.13;
CONSTANTS(:,23) = -25.0;
CONSTANTS(:,24) = 333.0;
CONSTANTS(:,25) = 0.10;
CONSTANTS(:,26) = 0.23;
CONSTANTS(:,27) = -44.8;
CONSTANTS(:,28) = 10.0;
STATES(:,8) = 0.000;
CONSTANTS(:,29) = 1000.0;
CONSTANTS(:,30) = 0.70;
CONSTANTS(:,31) = -0.56;
CONSTANTS(:,32) = 30;
CONSTANTS(:,33) = 6.0;
CONSTANTS(:,34) = 2000.0;
STATES(:,9) = 0.000;
STATES(:,10) = 0.100;
CONSTANTS(:,35) = 250.0;
CONSTANTS(:,36) = 1000.0;
CONSTANTS(:,37) = -0.2;
CONSTANTS(:,38) = -15.0;
CONSTANTS(:,39) = 1.0;
CONSTANTS(:,40) = -0.147;
CONSTANTS(:,41) = -35.3;
CONSTANTS(:,42) = 10.0;
CONSTANTS(:,43) = 50.0;
CONSTANTS(:,44) = -0.036;
CONSTANTS(:,45) = -38.0;
CONSTANTS(:,46) = 1.00;
CONSTANTS(:,47) = 0.167;
CONSTANTS(:,48) = -83.0;
CONSTANTS(:,49) = 0.20;
CONSTANTS(:,50) = 0.0;
CONSTANTS(:,51) = 20.0;
CONSTANTS(:,52) = 310.16;
CONSTANTS(:,53) = 8.314E-3;
CONSTANTS(:,54) = 0.09649;
CONSTANTS(:,55) = 1.8E3;
CONSTANTS(:,56) = 5400.0;
CONSTANTS(:,57) = 145E3;
CONSTANTS(:,58) = 140E3;
CONSTANTS(:,59) = 10E3;
CONSTANTS(:,60) = 1.0E-3;
CONSTANTS(:,61) = 0.002;
CONSTANTS(:,62) = 0.000;
CONSTANTS(:,63) = 0.3361;
CONSTANTS(:,64) = 0.3658;
CONSTANTS(:,65) = 0.4357;
CONSTANTS(:,66) = 0.5351;
CONSTANTS(:,67) = 0.6366;
CONSTANTS(:,68) = 0.7386;
CONSTANTS(:,69) = 0.8361;
CONSTANTS(:,70) = 0.9320;
CONSTANTS(:,71) = 1.0185;
STATES(:,11) = 0.023;
STATES(:,12) = 0.023;
STATES(:,13) = 0.023;
STATES(:,14) = 0.023;
STATES(:,15) = 0.023;
STATES(:,16) = 0.023;
STATES(:,17) = 0.023;
STATES(:,18) = 0.023;
STATES(:,19) = 0.023;
STATES(:,20) = 0.023;
CONSTANTS(:,72) = 0.30;
CONSTANTS(:,73) = 0.24;
CONSTANTS(:,74) = 10.0;
STATES(:,21) = 0.120;
STATES(:,22) = 0.120;
STATES(:,23) = 0.120;
STATES(:,24) = 0.120;
STATES(:,25) = 0.120;
STATES(:,26) = 0.120;
STATES(:,27) = 0.120;
STATES(:,28) = 0.120;
STATES(:,29) = 0.120;
STATES(:,30) = 0.120;
CONSTANTS(:,75) = 1.30;
CONSTANTS(:,76) = 0.12;
CONSTANTS(:,77) = 1.0;
CONSTANTS(:,78) = 0.074;
CONSTANTS(:,79) = 0.013;
CONSTANTS(:,80) = -0.20;
CONSTANTS(:,81) = -85.0;
CONSTANTS(:,82) = 310.16;
CONSTANTS(:,83) = 8.314E-3;
CONSTANTS(:,84) = 0.09649;
CONSTANTS(:,85) = 1.8E3;
CONSTANTS(:,86) = 13.5;
CONSTANTS(:,87) = 140.0E3;
STATES(:,31) = 1.01E4;
CONSTANTS(:,88) = 0.10;
CONSTANTS(:,89) = 3.0;
CONSTANTS(:,90) = 0.0;
CONSTANTS(:,91) = 0.12;
CONSTANTS(:,92) = 160.0;
CONSTANTS(:,93) = 675.0;
CONSTANTS(:,94) = 1.0;
CONSTANTS(:,95) = 187.5;
CONSTANTS(:,96) = 1.0E4;
CONSTANTS(:,97) = 4.0;
CONSTANTS(:,98) = 10.0;
CONSTANTS(:,99) = 1.0E5;
CONSTANTS(:,100) = 6E-4;
CONSTANTS(:,101) = 3.57;
CONSTANTS(:,102) = 2.7E-5;
CONSTANTS(:,103) = 2.31E-5;
STATES(:,32) = 200.0;
CONSTANTS(:,104) = 800.0;
CONSTANTS(:,105) = 10.0;
CONSTANTS(:,106) = 10.0;
CONSTANTS(:,107) = 4.0;
STATES(:,33) = 0.200;
CONSTANTS(:,108) = 2000.0;
CONSTANTS(:,109) = 6.4;
CONSTANTS(:,110) = 0.0;
CONSTANTS(:,111) = 4.0;
CONSTANTS(:,112) = 250.0;
CONSTANTS(:,113) = 0.0;
CONSTANTS(:,114) = 750.0;
CONSTANTS(:,115) = 36.0;
CONSTANTS(:,116) = 300.0;
STATES(:,34) = 0.200;
CONSTANTS(:,117) = 1500.0;
CONSTANTS(:,118) = 4.0;
CONSTANTS(:,119) = 0.85;
CONSTANTS(:,120) = 0.35;
STATES(:,35) = 0.300;
CONSTANTS(:,121) = 0.1;
CONSTANTS(:,122) = 0.5;
CONSTANTS(:,123) = 3.0;
CONSTANTS(:,124) = 3.0;
CONSTANTS(:,125) = 4.5;
CONSTANTS(:,126) = 1.0;
CONSTANTS(:,127) = 3.0;
CONSTANTS(:,128) = 0.3;
CONSTANTS(:,129) = 18.5;
CONSTANTS(:,130) = 100.0;
CONSTANTS(:,131) = 200.0;
CONSTANTS(:,132) = 20.0;
CONSTANTS(:,133) = 1.0;
CONSTANTS(:,134) = 0.01;
CONSTANTS(:,135) = 100.0;
CONSTANTS(:,136) = 5400.0;
CONSTANTS(:,137) = 145E3;
CONSTANTS(:,138) = 50.0;
CONSTANTS(:,139) = 0.074;
CONSTANTS(:,140) = 0.013;
CONSTANTS(:,141) = -0.20;
CONSTANTS(:,142) = -85.0;
CONSTANTS(:,143) = 310.16;
CONSTANTS(:,144) = 8.314E-3;
CONSTANTS(:,145) = 0.09649;
CONSTANTS(:,146) = 1.8E3;
CONSTANTS(:,147) = 13.5;
CONSTANTS(:,148) = 140.0E3;
STATES(:,36) = 1.01E4;
CONSTANTS(:,149) = 0.10;
CONSTANTS(:,150) = 3.0;
CONSTANTS(:,151) = 0.0;
CONSTANTS(:,152) = 0.12;
CONSTANTS(:,153) = 160.0;
CONSTANTS(:,154) = 675.0;
CONSTANTS(:,155) = 1.0;
CONSTANTS(:,156) = 187.5;
CONSTANTS(:,157) = 1.0E4;
CONSTANTS(:,158) = 4.0;
CONSTANTS(:,159) = 10.0;
CONSTANTS(:,160) = 1.0E5;
CONSTANTS(:,161) = 6E-4;
CONSTANTS(:,162) = 3.57;
CONSTANTS(:,163) = 2.7E-5;
CONSTANTS(:,164) = 2.31E-5;
STATES(:,37) = 200.0;
CONSTANTS(:,165) = 800.0;
CONSTANTS(:,166) = 10.0;
CONSTANTS(:,167) = 10.0;
CONSTANTS(:,168) = 4.0;
STATES(:,38) = 0.200;
CONSTANTS(:,169) = 2000.0;
CONSTANTS(:,170) = 6.4;
CONSTANTS(:,171) = 0.0;
CONSTANTS(:,172) = 4.0;
CONSTANTS(:,173) = 250.0;
CONSTANTS(:,174) = 0.0;
CONSTANTS(:,175) = 750.0;
CONSTANTS(:,176) = 36.0;
CONSTANTS(:,177) = 300.0;
STATES(:,39) = 0.200;
CONSTANTS(:,178) = 1500.0;
CONSTANTS(:,179) = 4.0;
CONSTANTS(:,180) = 0.85;
CONSTANTS(:,181) = 0.35;
STATES(:,40) = 0.300;
CONSTANTS(:,182) = 0.1;
CONSTANTS(:,183) = 0.5;
CONSTANTS(:,184) = 3.0;
CONSTANTS(:,185) = 3.0;
CONSTANTS(:,186) = 4.5;
CONSTANTS(:,187) = 1.0;
CONSTANTS(:,188) = 3.0;
CONSTANTS(:,189) = 0.3;
CONSTANTS(:,190) = 18.5;
CONSTANTS(:,191) = 100.0;
CONSTANTS(:,192) = 200.0;
CONSTANTS(:,193) = 20.0;
CONSTANTS(:,194) = 1.0;
CONSTANTS(:,195) = 0.01;
CONSTANTS(:,196) = 100.0;
CONSTANTS(:,197) = 5400.0;
CONSTANTS(:,198) = 145E3;
CONSTANTS(:,199) = 50.0;
CONSTANTS(:,200) = 0.074;
CONSTANTS(:,201) = 0.013;
CONSTANTS(:,202) = -0.20;
CONSTANTS(:,203) = -85.0;
CONSTANTS(:,204) = 310.16;
CONSTANTS(:,205) = 8.314E-3;
CONSTANTS(:,206) = 0.09649;
CONSTANTS(:,207) = 1.8E3;
CONSTANTS(:,208) = 13.5;
CONSTANTS(:,209) = 140.0E3;
STATES(:,41) = 1.01E4;
CONSTANTS(:,210) = 0.10;
CONSTANTS(:,211) = 3.0;
CONSTANTS(:,212) = 0.0;
CONSTANTS(:,213) = 0.12;
CONSTANTS(:,214) = 160.0;
CONSTANTS(:,215) = 675.0;
CONSTANTS(:,216) = 1.0;
CONSTANTS(:,217) = 187.5;
CONSTANTS(:,218) = 1.0E4;
CONSTANTS(:,219) = 4.0;
CONSTANTS(:,220) = 10.0;
CONSTANTS(:,221) = 1.0E5;
CONSTANTS(:,222) = 6E-4;
CONSTANTS(:,223) = 3.57;
CONSTANTS(:,224) = 2.7E-5;
CONSTANTS(:,225) = 2.31E-5;
STATES(:,42) = 200.0;
CONSTANTS(:,226) = 800.0;
CONSTANTS(:,227) = 10.0;
CONSTANTS(:,228) = 10.0;
CONSTANTS(:,229) = 4.0;
STATES(:,43) = 0.200;
CONSTANTS(:,230) = 2000.0;
CONSTANTS(:,231) = 6.4;
CONSTANTS(:,232) = 0.0;
CONSTANTS(:,233) = 4.0;
CONSTANTS(:,234) = 250.0;
CONSTANTS(:,235) = 0.0;
CONSTANTS(:,236) = 750.0;
CONSTANTS(:,237) = 36.0;
CONSTANTS(:,238) = 300.0;
STATES(:,44) = 0.200;
CONSTANTS(:,239) = 1500.0;
CONSTANTS(:,240) = 4.0;
CONSTANTS(:,241) = 0.85;
CONSTANTS(:,242) = 0.35;
STATES(:,45) = 0.300;
CONSTANTS(:,243) = 0.1;
CONSTANTS(:,244) = 0.5;
CONSTANTS(:,245) = 3.0;
CONSTANTS(:,246) = 3.0;
CONSTANTS(:,247) = 4.5;
CONSTANTS(:,248) = 1.0;
CONSTANTS(:,249) = 3.0;
CONSTANTS(:,250) = 0.3;
CONSTANTS(:,251) = 18.5;
CONSTANTS(:,252) = 100.0;
CONSTANTS(:,253) = 200.0;
CONSTANTS(:,254) = 20.0;
CONSTANTS(:,255) = 1.0;
CONSTANTS(:,256) = 0.01;
CONSTANTS(:,257) = 100.0;
CONSTANTS(:,258) = 5400.0;
CONSTANTS(:,259) = 145E3;
CONSTANTS(:,260) = 50.0;
CONSTANTS(:,261) = 0.074;
CONSTANTS(:,262) = 0.013;
CONSTANTS(:,263) = -0.20;
CONSTANTS(:,264) = -85.0;
CONSTANTS(:,265) = 310.16;
CONSTANTS(:,266) = 8.314E-3;
CONSTANTS(:,267) = 0.09649;
CONSTANTS(:,268) = 1.8E3;
CONSTANTS(:,269) = 13.5;
CONSTANTS(:,270) = 140.0E3;
STATES(:,46) = 1.01E4;
CONSTANTS(:,271) = 0.10;
CONSTANTS(:,272) = 3.0;
CONSTANTS(:,273) = 0.0;
CONSTANTS(:,274) = 0.12;
CONSTANTS(:,275) = 160.0;
CONSTANTS(:,276) = 675.0;
CONSTANTS(:,277) = 1.0;
CONSTANTS(:,278) = 187.5;
CONSTANTS(:,279) = 1.0E4;
CONSTANTS(:,280) = 4.0;
CONSTANTS(:,281) = 10.0;
CONSTANTS(:,282) = 1.0E5;
CONSTANTS(:,283) = 6E-4;
CONSTANTS(:,284) = 3.57;
CONSTANTS(:,285) = 2.7E-5;
CONSTANTS(:,286) = 2.31E-5;
STATES(:,47) = 200.0;
CONSTANTS(:,287) = 800.0;
CONSTANTS(:,288) = 10.0;
CONSTANTS(:,289) = 10.0;
CONSTANTS(:,290) = 4.0;
STATES(:,48) = 0.200;
CONSTANTS(:,291) = 2000.0;
CONSTANTS(:,292) = 6.4;
CONSTANTS(:,293) = 0.0;
CONSTANTS(:,294) = 4.0;
CONSTANTS(:,295) = 250.0;
CONSTANTS(:,296) = 0.0;
CONSTANTS(:,297) = 750.0;
CONSTANTS(:,298) = 36.0;
CONSTANTS(:,299) = 300.0;
STATES(:,49) = 0.200;
CONSTANTS(:,300) = 1500.0;
CONSTANTS(:,301) = 4.0;
CONSTANTS(:,302) = 0.85;
CONSTANTS(:,303) = 0.35;
STATES(:,50) = 0.300;
CONSTANTS(:,304) = 0.1;
CONSTANTS(:,305) = 0.5;
CONSTANTS(:,306) = 3.0;
CONSTANTS(:,307) = 3.0;
CONSTANTS(:,308) = 4.5;
CONSTANTS(:,309) = 1.0;
CONSTANTS(:,310) = 3.0;
CONSTANTS(:,311) = 0.3;
CONSTANTS(:,312) = 18.5;
CONSTANTS(:,313) = 100.0;
CONSTANTS(:,314) = 200.0;
CONSTANTS(:,315) = 20.0;
CONSTANTS(:,316) = 1.0;
CONSTANTS(:,317) = 0.01;
CONSTANTS(:,318) = 100.0;
CONSTANTS(:,319) = 5400.0;
CONSTANTS(:,320) = 145E3;
CONSTANTS(:,321) = 50.0;
CONSTANTS(:,322) = 0.074;
CONSTANTS(:,323) = 0.013;
CONSTANTS(:,324) = -0.20;
CONSTANTS(:,325) = -85.0;
CONSTANTS(:,326) = 310.16;
CONSTANTS(:,327) = 8.314E-3;
CONSTANTS(:,328) = 0.09649;
CONSTANTS(:,329) = 1.8E3;
CONSTANTS(:,330) = 13.5;
CONSTANTS(:,331) = 140.0E3;
STATES(:,51) = 1.01E4;
CONSTANTS(:,332) = 0.10;
CONSTANTS(:,333) = 3.0;
CONSTANTS(:,334) = 0.0;
CONSTANTS(:,335) = 0.12;
CONSTANTS(:,336) = 160.0;
CONSTANTS(:,337) = 675.0;
CONSTANTS(:,338) = 1.0;
CONSTANTS(:,339) = 187.5;
CONSTANTS(:,340) = 1.0E4;
CONSTANTS(:,341) = 4.0;
CONSTANTS(:,342) = 10.0;
CONSTANTS(:,343) = 1.0E5;
CONSTANTS(:,344) = 6E-4;
CONSTANTS(:,345) = 3.57;
CONSTANTS(:,346) = 2.7E-5;
CONSTANTS(:,347) = 2.31E-5;
STATES(:,52) = 200.0;
CONSTANTS(:,348) = 800.0;
CONSTANTS(:,349) = 10.0;
CONSTANTS(:,350) = 10.0;
CONSTANTS(:,351) = 4.0;
STATES(:,53) = 0.200;
CONSTANTS(:,352) = 2000.0;
CONSTANTS(:,353) = 6.4;
CONSTANTS(:,354) = 0.0;
CONSTANTS(:,355) = 4.0;
CONSTANTS(:,356) = 250.0;
CONSTANTS(:,357) = 0.0;
CONSTANTS(:,358) = 750.0;
CONSTANTS(:,359) = 36.0;
CONSTANTS(:,360) = 300.0;
STATES(:,54) = 0.200;
CONSTANTS(:,361) = 1500.0;
CONSTANTS(:,362) = 4.0;
CONSTANTS(:,363) = 0.85;
CONSTANTS(:,364) = 0.35;
STATES(:,55) = 0.300;
CONSTANTS(:,365) = 0.1;
CONSTANTS(:,366) = 0.5;
CONSTANTS(:,367) = 3.0;
CONSTANTS(:,368) = 3.0;
CONSTANTS(:,369) = 4.5;
CONSTANTS(:,370) = 1.0;
CONSTANTS(:,371) = 3.0;
CONSTANTS(:,372) = 0.3;
CONSTANTS(:,373) = 18.5;
CONSTANTS(:,374) = 100.0;
CONSTANTS(:,375) = 200.0;
CONSTANTS(:,376) = 20.0;
CONSTANTS(:,377) = 1.0;
CONSTANTS(:,378) = 0.01;
CONSTANTS(:,379) = 100.0;
CONSTANTS(:,380) = 5400.0;
CONSTANTS(:,381) = 145E3;
CONSTANTS(:,382) = 50.0;
CONSTANTS(:,383) = 0.074;
CONSTANTS(:,384) = 0.013;
CONSTANTS(:,385) = -0.20;
CONSTANTS(:,386) = -85.0;
CONSTANTS(:,387) = 310.16;
CONSTANTS(:,388) = 8.314E-3;
CONSTANTS(:,389) = 0.09649;
CONSTANTS(:,390) = 1.8E3;
CONSTANTS(:,391) = 13.5;
CONSTANTS(:,392) = 140.0E3;
STATES(:,56) = 1.01E4;
CONSTANTS(:,393) = 0.10;
CONSTANTS(:,394) = 3.0;
CONSTANTS(:,395) = 0.0;
CONSTANTS(:,396) = 0.12;
CONSTANTS(:,397) = 160.0;
CONSTANTS(:,398) = 675.0;
CONSTANTS(:,399) = 1.0;
CONSTANTS(:,400) = 187.5;
CONSTANTS(:,401) = 1.0E4;
CONSTANTS(:,402) = 4.0;
CONSTANTS(:,403) = 10.0;
CONSTANTS(:,404) = 1.0E5;
CONSTANTS(:,405) = 6E-4;
CONSTANTS(:,406) = 3.57;
CONSTANTS(:,407) = 2.7E-5;
CONSTANTS(:,408) = 2.31E-5;
STATES(:,57) = 200.0;
CONSTANTS(:,409) = 800.0;
CONSTANTS(:,410) = 10.0;
CONSTANTS(:,411) = 10.0;
CONSTANTS(:,412) = 4.0;
STATES(:,58) = 0.200;
CONSTANTS(:,413) = 2000.0;
CONSTANTS(:,414) = 6.4;
CONSTANTS(:,415) = 0.0;
CONSTANTS(:,416) = 4.0;
CONSTANTS(:,417) = 250.0;
CONSTANTS(:,418) = 0.0;
CONSTANTS(:,419) = 750.0;
CONSTANTS(:,420) = 36.0;
CONSTANTS(:,421) = 300.0;
STATES(:,59) = 0.200;
CONSTANTS(:,422) = 1500.0;
CONSTANTS(:,423) = 4.0;
CONSTANTS(:,424) = 0.85;
CONSTANTS(:,425) = 0.35;
STATES(:,60) = 0.300;
CONSTANTS(:,426) = 0.1;
CONSTANTS(:,427) = 0.5;
CONSTANTS(:,428) = 3.0;
CONSTANTS(:,429) = 3.0;
CONSTANTS(:,430) = 4.5;
CONSTANTS(:,431) = 1.0;
CONSTANTS(:,432) = 3.0;
CONSTANTS(:,433) = 0.3;
CONSTANTS(:,434) = 18.5;
CONSTANTS(:,435) = 100.0;
CONSTANTS(:,436) = 200.0;
CONSTANTS(:,437) = 20.0;
CONSTANTS(:,438) = 1.0;
CONSTANTS(:,439) = 0.01;
CONSTANTS(:,440) = 100.0;
CONSTANTS(:,441) = 5400.0;
CONSTANTS(:,442) = 145E3;
CONSTANTS(:,443) = 50.0;
CONSTANTS(:,444) = 0.074;
CONSTANTS(:,445) = 0.013;
CONSTANTS(:,446) = -0.20;
CONSTANTS(:,447) = -85.0;
CONSTANTS(:,448) = 310.16;
CONSTANTS(:,449) = 8.314E-3;
CONSTANTS(:,450) = 0.09649;
CONSTANTS(:,451) = 1.8E3;
CONSTANTS(:,452) = 13.5;
CONSTANTS(:,453) = 140.0E3;
STATES(:,61) = 1.01E4;
CONSTANTS(:,454) = 0.10;
CONSTANTS(:,455) = 3.0;
CONSTANTS(:,456) = 0.0;
CONSTANTS(:,457) = 0.12;
CONSTANTS(:,458) = 160.0;
CONSTANTS(:,459) = 675.0;
CONSTANTS(:,460) = 1.0;
CONSTANTS(:,461) = 187.5;
CONSTANTS(:,462) = 1.0E4;
CONSTANTS(:,463) = 4.0;
CONSTANTS(:,464) = 10.0;
CONSTANTS(:,465) = 1.0E5;
CONSTANTS(:,466) = 6E-4;
CONSTANTS(:,467) = 3.57;
CONSTANTS(:,468) = 2.7E-5;
CONSTANTS(:,469) = 2.31E-5;
STATES(:,62) = 200.0;
CONSTANTS(:,470) = 800.0;
CONSTANTS(:,471) = 10.0;
CONSTANTS(:,472) = 10.0;
CONSTANTS(:,473) = 4.0;
STATES(:,63) = 0.200;
CONSTANTS(:,474) = 2000.0;
CONSTANTS(:,475) = 6.4;
CONSTANTS(:,476) = 0.0;
CONSTANTS(:,477) = 4.0;
CONSTANTS(:,478) = 250.0;
CONSTANTS(:,479) = 0.0;
CONSTANTS(:,480) = 750.0;
CONSTANTS(:,481) = 36.0;
CONSTANTS(:,482) = 300.0;
STATES(:,64) = 0.200;
CONSTANTS(:,483) = 1500.0;
CONSTANTS(:,484) = 4.0;
CONSTANTS(:,485) = 0.85;
CONSTANTS(:,486) = 0.35;
STATES(:,65) = 0.300;
CONSTANTS(:,487) = 0.1;
CONSTANTS(:,488) = 0.5;
CONSTANTS(:,489) = 3.0;
CONSTANTS(:,490) = 3.0;
CONSTANTS(:,491) = 4.5;
CONSTANTS(:,492) = 1.0;
CONSTANTS(:,493) = 3.0;
CONSTANTS(:,494) = 0.3;
CONSTANTS(:,495) = 18.5;
CONSTANTS(:,496) = 100.0;
CONSTANTS(:,497) = 200.0;
CONSTANTS(:,498) = 20.0;
CONSTANTS(:,499) = 1.0;
CONSTANTS(:,500) = 0.01;
CONSTANTS(:,501) = 100.0;
CONSTANTS(:,502) = 5400.0;
CONSTANTS(:,503) = 145E3;
CONSTANTS(:,504) = 50.0;
CONSTANTS(:,505) = 0.074;
CONSTANTS(:,506) = 0.013;
CONSTANTS(:,507) = -0.20;
CONSTANTS(:,508) = -85.0;
CONSTANTS(:,509) = 310.16;
CONSTANTS(:,510) = 8.314E-3;
CONSTANTS(:,511) = 0.09649;
CONSTANTS(:,512) = 1.8E3;
CONSTANTS(:,513) = 13.5;
CONSTANTS(:,514) = 140.0E3;
STATES(:,66) = 1.01E4;
CONSTANTS(:,515) = 0.10;
CONSTANTS(:,516) = 3.0;
CONSTANTS(:,517) = 0.0;
CONSTANTS(:,518) = 0.12;
CONSTANTS(:,519) = 160.0;
CONSTANTS(:,520) = 675.0;
CONSTANTS(:,521) = 1.0;
CONSTANTS(:,522) = 187.5;
CONSTANTS(:,523) = 1.0E4;
CONSTANTS(:,524) = 4.0;
CONSTANTS(:,525) = 10.0;
CONSTANTS(:,526) = 1.0E5;
CONSTANTS(:,527) = 6E-4;
CONSTANTS(:,528) = 3.57;
CONSTANTS(:,529) = 2.7E-5;
CONSTANTS(:,530) = 2.31E-5;
STATES(:,67) = 200.0;
CONSTANTS(:,531) = 800.0;
CONSTANTS(:,532) = 10.0;
CONSTANTS(:,533) = 10.0;
CONSTANTS(:,534) = 4.0;
STATES(:,68) = 0.200;
CONSTANTS(:,535) = 2000.0;
CONSTANTS(:,536) = 6.4;
CONSTANTS(:,537) = 0.0;
CONSTANTS(:,538) = 4.0;
CONSTANTS(:,539) = 250.0;
CONSTANTS(:,540) = 0.0;
CONSTANTS(:,541) = 750.0;
CONSTANTS(:,542) = 36.0;
CONSTANTS(:,543) = 300.0;
STATES(:,69) = 0.200;
CONSTANTS(:,544) = 1500.0;
CONSTANTS(:,545) = 4.0;
CONSTANTS(:,546) = 0.85;
CONSTANTS(:,547) = 0.35;
STATES(:,70) = 0.300;
CONSTANTS(:,548) = 0.1;
CONSTANTS(:,549) = 0.5;
CONSTANTS(:,550) = 3.0;
CONSTANTS(:,551) = 3.0;
CONSTANTS(:,552) = 4.5;
CONSTANTS(:,553) = 1.0;
CONSTANTS(:,554) = 3.0;
CONSTANTS(:,555) = 0.3;
CONSTANTS(:,556) = 18.5;
CONSTANTS(:,557) = 100.0;
CONSTANTS(:,558) = 200.0;
CONSTANTS(:,559) = 20.0;
CONSTANTS(:,560) = 1.0;
CONSTANTS(:,561) = 0.01;
CONSTANTS(:,562) = 100.0;
CONSTANTS(:,563) = 5400.0;
CONSTANTS(:,564) = 145E3;
CONSTANTS(:,565) = 50.0;
CONSTANTS(:,566) = 0.074;
CONSTANTS(:,567) = 0.013;
CONSTANTS(:,568) = -0.20;
CONSTANTS(:,569) = -85.0;
CONSTANTS(:,570) = 310.16;
CONSTANTS(:,571) = 8.314E-3;
CONSTANTS(:,572) = 0.09649;
CONSTANTS(:,573) = 1.8E3;
CONSTANTS(:,574) = 13.5;
CONSTANTS(:,575) = 140.0E3;
STATES(:,71) = 1.01E4;
CONSTANTS(:,576) = 0.10;
CONSTANTS(:,577) = 3.0;
CONSTANTS(:,578) = 0.0;
CONSTANTS(:,579) = 0.12;
CONSTANTS(:,580) = 160.0;
CONSTANTS(:,581) = 675.0;
CONSTANTS(:,582) = 1.0;
CONSTANTS(:,583) = 187.5;
CONSTANTS(:,584) = 1.0E4;
CONSTANTS(:,585) = 4.0;
CONSTANTS(:,586) = 10.0;
CONSTANTS(:,587) = 1.0E5;
CONSTANTS(:,588) = 6E-4;
CONSTANTS(:,589) = 3.57;
CONSTANTS(:,590) = 2.7E-5;
CONSTANTS(:,591) = 2.31E-5;
STATES(:,72) = 200.0;
CONSTANTS(:,592) = 800.0;
CONSTANTS(:,593) = 10.0;
CONSTANTS(:,594) = 10.0;
CONSTANTS(:,595) = 4.0;
STATES(:,73) = 0.200;
CONSTANTS(:,596) = 2000.0;
CONSTANTS(:,597) = 6.4;
CONSTANTS(:,598) = 0.0;
CONSTANTS(:,599) = 4.0;
CONSTANTS(:,600) = 250.0;
CONSTANTS(:,601) = 0.0;
CONSTANTS(:,602) = 750.0;
CONSTANTS(:,603) = 36.0;
CONSTANTS(:,604) = 300.0;
STATES(:,74) = 0.200;
CONSTANTS(:,605) = 1500.0;
CONSTANTS(:,606) = 4.0;
CONSTANTS(:,607) = 0.85;
CONSTANTS(:,608) = 0.35;
STATES(:,75) = 0.300;
CONSTANTS(:,609) = 0.1;
CONSTANTS(:,610) = 0.5;
CONSTANTS(:,611) = 3.0;
CONSTANTS(:,612) = 3.0;
CONSTANTS(:,613) = 4.5;
CONSTANTS(:,614) = 1.0;
CONSTANTS(:,615) = 3.0;
CONSTANTS(:,616) = 0.3;
CONSTANTS(:,617) = 18.5;
CONSTANTS(:,618) = 100.0;
CONSTANTS(:,619) = 200.0;
CONSTANTS(:,620) = 20.0;
CONSTANTS(:,621) = 1.0;
CONSTANTS(:,622) = 0.01;
CONSTANTS(:,623) = 100.0;
CONSTANTS(:,624) = 5400.0;
CONSTANTS(:,625) = 145E3;
CONSTANTS(:,626) = 50.0;
CONSTANTS(:,627) = 0.074;
CONSTANTS(:,628) = 0.013;
CONSTANTS(:,629) = -0.20;
CONSTANTS(:,630) = -85.0;
CONSTANTS(:,631) = 310.16;
CONSTANTS(:,632) = 8.314E-3;
CONSTANTS(:,633) = 0.09649;
CONSTANTS(:,634) = 1.8E3;
CONSTANTS(:,635) = 13.5;
CONSTANTS(:,636) = 140.0E3;
STATES(:,76) = 1.01E4;
CONSTANTS(:,637) = 0.10;
CONSTANTS(:,638) = 3.0;
CONSTANTS(:,639) = 0.0;
CONSTANTS(:,640) = 0.12;
CONSTANTS(:,641) = 160.0;
CONSTANTS(:,642) = 675.0;
CONSTANTS(:,643) = 1.0;
CONSTANTS(:,644) = 187.5;
CONSTANTS(:,645) = 1.0E4;
CONSTANTS(:,646) = 4.0;
CONSTANTS(:,647) = 10.0;
CONSTANTS(:,648) = 1.0E5;
CONSTANTS(:,649) = 6E-4;
CONSTANTS(:,650) = 3.57;
CONSTANTS(:,651) = 2.7E-5;
CONSTANTS(:,652) = 2.31E-5;
STATES(:,77) = 200.0;
CONSTANTS(:,653) = 800.0;
CONSTANTS(:,654) = 10.0;
CONSTANTS(:,655) = 10.0;
CONSTANTS(:,656) = 4.0;
STATES(:,78) = 0.200;
CONSTANTS(:,657) = 2000.0;
CONSTANTS(:,658) = 6.4;
CONSTANTS(:,659) = 0.0;
CONSTANTS(:,660) = 4.0;
CONSTANTS(:,661) = 250.0;
CONSTANTS(:,662) = 0.0;
CONSTANTS(:,663) = 750.0;
CONSTANTS(:,664) = 36.0;
CONSTANTS(:,665) = 300.0;
STATES(:,79) = 0.200;
CONSTANTS(:,666) = 1500.0;
CONSTANTS(:,667) = 4.0;
CONSTANTS(:,668) = 0.85;
CONSTANTS(:,669) = 0.35;
STATES(:,80) = 0.300;
CONSTANTS(:,670) = 0.1;
CONSTANTS(:,671) = 0.5;
CONSTANTS(:,672) = 3.0;
CONSTANTS(:,673) = 3.0;
CONSTANTS(:,674) = 4.5;
CONSTANTS(:,675) = 1.0;
CONSTANTS(:,676) = 3.0;
CONSTANTS(:,677) = 0.3;
CONSTANTS(:,678) = 18.5;
CONSTANTS(:,679) = 100.0;
CONSTANTS(:,680) = 200.0;
CONSTANTS(:,681) = 20.0;
CONSTANTS(:,682) = 1.0;
CONSTANTS(:,683) = 0.01;
CONSTANTS(:,684) = 100.0;
CONSTANTS(:,685) = 5400.0;
CONSTANTS(:,686) = 145E3;
CONSTANTS(:,687) = 50.0;
CONSTANTS(:,688) = CONSTANTS(:,5);
CONSTANTS(:,689) = (( CONSTANTS(:,53).*CONSTANTS(:,52))./CONSTANTS(:,54)).*log(CONSTANTS(:,56)./CONSTANTS(:,57));
CONSTANTS(:,690) = (( CONSTANTS(:,53).*CONSTANTS(:,52))./CONSTANTS(:,54)).*log(CONSTANTS(:,58)./CONSTANTS(:,59));
CONSTANTS(:,691) = CONSTANTS(:,131)./CONSTANTS(:,130);
CONSTANTS(:,692) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((1.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,693) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((2.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,694) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((3.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,695) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((4.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,696) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((5.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,697) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((6.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,698) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((7.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,699) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((8.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,700) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((9.00000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,701) = CONSTANTS(:,72)+ (CONSTANTS(:,73) - CONSTANTS(:,72)).*((10.0000 - 1.00000)./(CONSTANTS(:,74) - 1.00000));
CONSTANTS(:,702) = CONSTANTS(:,119);
CONSTANTS(:,703) = CONSTANTS(:,180);
CONSTANTS(:,704) = CONSTANTS(:,241);
CONSTANTS(:,705) = CONSTANTS(:,302);
CONSTANTS(:,706) = CONSTANTS(:,363);
CONSTANTS(:,707) = CONSTANTS(:,424);
CONSTANTS(:,708) = CONSTANTS(:,485);
CONSTANTS(:,709) = CONSTANTS(:,546);
CONSTANTS(:,710) = CONSTANTS(:,607);
CONSTANTS(:,711) = CONSTANTS(:,668);
CONSTANTS(:,712) = CONSTANTS(:,132)./CONSTANTS(:,130);
CONSTANTS(:,713) = CONSTANTS(:,131)./CONSTANTS(:,133);
CONSTANTS(:,714) = CONSTANTS(:,132)./CONSTANTS(:,133);
CONSTANTS(:,715) = CONSTANTS(:,130)./CONSTANTS(:,133);
CONSTANTS(:,716) = CONSTANTS(:,138)./CONSTANTS(:,74);
CONSTANTS(:,717) = CONSTANTS(:,192)./CONSTANTS(:,191);
CONSTANTS(:,718) = CONSTANTS(:,193)./CONSTANTS(:,191);
CONSTANTS(:,719) = CONSTANTS(:,192)./CONSTANTS(:,194);
CONSTANTS(:,720) = CONSTANTS(:,193)./CONSTANTS(:,194);
CONSTANTS(:,721) = CONSTANTS(:,191)./CONSTANTS(:,194);
CONSTANTS(:,722) = CONSTANTS(:,199)./CONSTANTS(:,74);
CONSTANTS(:,723) = CONSTANTS(:,253)./CONSTANTS(:,252);
CONSTANTS(:,724) = CONSTANTS(:,254)./CONSTANTS(:,252);
CONSTANTS(:,725) = CONSTANTS(:,253)./CONSTANTS(:,255);
CONSTANTS(:,726) = CONSTANTS(:,254)./CONSTANTS(:,255);
CONSTANTS(:,727) = CONSTANTS(:,252)./CONSTANTS(:,255);
CONSTANTS(:,728) = CONSTANTS(:,260)./CONSTANTS(:,74);
CONSTANTS(:,729) = CONSTANTS(:,314)./CONSTANTS(:,313);
CONSTANTS(:,730) = CONSTANTS(:,315)./CONSTANTS(:,313);
CONSTANTS(:,731) = CONSTANTS(:,314)./CONSTANTS(:,316);
CONSTANTS(:,732) = CONSTANTS(:,315)./CONSTANTS(:,316);
CONSTANTS(:,733) = CONSTANTS(:,313)./CONSTANTS(:,316);
CONSTANTS(:,734) = CONSTANTS(:,321)./CONSTANTS(:,74);
CONSTANTS(:,735) = CONSTANTS(:,375)./CONSTANTS(:,374);
CONSTANTS(:,736) = CONSTANTS(:,376)./CONSTANTS(:,374);
CONSTANTS(:,737) = CONSTANTS(:,375)./CONSTANTS(:,377);
CONSTANTS(:,738) = CONSTANTS(:,376)./CONSTANTS(:,377);
CONSTANTS(:,739) = CONSTANTS(:,374)./CONSTANTS(:,377);
CONSTANTS(:,740) = CONSTANTS(:,382)./CONSTANTS(:,74);
CONSTANTS(:,741) = CONSTANTS(:,436)./CONSTANTS(:,435);
CONSTANTS(:,742) = CONSTANTS(:,437)./CONSTANTS(:,435);
CONSTANTS(:,743) = CONSTANTS(:,436)./CONSTANTS(:,438);
CONSTANTS(:,744) = CONSTANTS(:,437)./CONSTANTS(:,438);
CONSTANTS(:,745) = CONSTANTS(:,435)./CONSTANTS(:,438);
CONSTANTS(:,746) = CONSTANTS(:,443)./CONSTANTS(:,74);
CONSTANTS(:,747) = CONSTANTS(:,497)./CONSTANTS(:,496);
CONSTANTS(:,748) = CONSTANTS(:,498)./CONSTANTS(:,496);
CONSTANTS(:,749) = CONSTANTS(:,497)./CONSTANTS(:,499);
CONSTANTS(:,750) = CONSTANTS(:,498)./CONSTANTS(:,499);
CONSTANTS(:,751) = CONSTANTS(:,496)./CONSTANTS(:,499);
CONSTANTS(:,752) = CONSTANTS(:,504)./CONSTANTS(:,74);
CONSTANTS(:,753) = CONSTANTS(:,558)./CONSTANTS(:,557);
CONSTANTS(:,754) = CONSTANTS(:,559)./CONSTANTS(:,557);
CONSTANTS(:,755) = CONSTANTS(:,558)./CONSTANTS(:,560);
CONSTANTS(:,756) = CONSTANTS(:,559)./CONSTANTS(:,560);
CONSTANTS(:,757) = CONSTANTS(:,557)./CONSTANTS(:,560);
CONSTANTS(:,758) = CONSTANTS(:,565)./CONSTANTS(:,74);
CONSTANTS(:,759) = CONSTANTS(:,619)./CONSTANTS(:,618);
CONSTANTS(:,760) = CONSTANTS(:,620)./CONSTANTS(:,618);
CONSTANTS(:,761) = CONSTANTS(:,619)./CONSTANTS(:,621);
CONSTANTS(:,762) = CONSTANTS(:,620)./CONSTANTS(:,621);
CONSTANTS(:,763) = CONSTANTS(:,618)./CONSTANTS(:,621);
CONSTANTS(:,764) = CONSTANTS(:,626)./CONSTANTS(:,74);
CONSTANTS(:,765) = CONSTANTS(:,680)./CONSTANTS(:,679);
CONSTANTS(:,766) = CONSTANTS(:,681)./CONSTANTS(:,679);
CONSTANTS(:,767) = CONSTANTS(:,680)./CONSTANTS(:,682);
CONSTANTS(:,768) = CONSTANTS(:,681)./CONSTANTS(:,682);
CONSTANTS(:,769) = CONSTANTS(:,679)./CONSTANTS(:,682);
CONSTANTS(:,770) = CONSTANTS(:,687)./CONSTANTS(:,74);
if (isempty(STATES)), warning('Initial values for states not set');, end
end
function [RATES, ALGEBRAIC] = computeRates(VOI, STATES, CONSTANTS)
global algebraicVariableCount;
statesSize = size(STATES);
statesColumnCount = statesSize(2);
if ( statesColumnCount == 1)
STATES = STATES';
ALGEBRAIC = zeros(1, algebraicVariableCount);
utilOnes = 1;
else
statesRowCount = statesSize(1);
ALGEBRAIC = zeros(statesRowCount, algebraicVariableCount);
RATES = zeros(statesRowCount, statesColumnCount);
utilOnes = ones(statesRowCount, 1);
end
RATES(:,1) = CONSTANTS(:,3).*STATES(:,3).*STATES(:,4).*(1.00000 - STATES(:,1)) - CONSTANTS(:,4).*STATES(:,1);
ALGEBRAIC(:,1) = 1.00000./(1.00000+exp( CONSTANTS(:,6).*(STATES(:,2) - CONSTANTS(:,7))));
RATES(:,3) = (ALGEBRAIC(:,1) - STATES(:,3))./CONSTANTS(:,688);
ALGEBRAIC(:,10) = CONSTANTS(:,117)./(1.00000+power(CONSTANTS(:,120)./STATES(:,11), CONSTANTS(:,118)));
RATES(:,35) = CONSTANTS(:,702).*(1.00000 - STATES(:,35)) - ALGEBRAIC(:,10).*STATES(:,35);
ALGEBRAIC(:,11) = CONSTANTS(:,178)./(1.00000+power(CONSTANTS(:,181)./STATES(:,12), CONSTANTS(:,179)));
RATES(:,40) = CONSTANTS(:,703).*(1.00000 - STATES(:,40)) - ALGEBRAIC(:,11).*STATES(:,40);
ALGEBRAIC(:,12) = CONSTANTS(:,239)./(1.00000+power(CONSTANTS(:,242)./STATES(:,13), CONSTANTS(:,240)));
RATES(:,45) = CONSTANTS(:,704).*(1.00000 - STATES(:,45)) - ALGEBRAIC(:,12).*STATES(:,45);
ALGEBRAIC(:,13) = CONSTANTS(:,300)./(1.00000+power(CONSTANTS(:,303)./STATES(:,14), CONSTANTS(:,301)));
RATES(:,50) = CONSTANTS(:,705).*(1.00000 - STATES(:,50)) - ALGEBRAIC(:,13).*STATES(:,50);
ALGEBRAIC(:,14) = CONSTANTS(:,361)./(1.00000+power(CONSTANTS(:,364)./STATES(:,15), CONSTANTS(:,362)));
RATES(:,55) = CONSTANTS(:,706).*(1.00000 - STATES(:,55)) - ALGEBRAIC(:,14).*STATES(:,55);
ALGEBRAIC(:,15) = CONSTANTS(:,422)./(1.00000+power(CONSTANTS(:,425)./STATES(:,16), CONSTANTS(:,423)));
RATES(:,60) = CONSTANTS(:,707).*(1.00000 - STATES(:,60)) - ALGEBRAIC(:,15).*STATES(:,60);
ALGEBRAIC(:,16) = CONSTANTS(:,483)./(1.00000+power(CONSTANTS(:,486)./STATES(:,17), CONSTANTS(:,484)));
RATES(:,65) = CONSTANTS(:,708).*(1.00000 - STATES(:,65)) - ALGEBRAIC(:,16).*STATES(:,65);
ALGEBRAIC(:,17) = CONSTANTS(:,544)./(1.00000+power(CONSTANTS(:,547)./STATES(:,18), CONSTANTS(:,545)));
RATES(:,70) = CONSTANTS(:,709).*(1.00000 - STATES(:,70)) - ALGEBRAIC(:,17).*STATES(:,70);
ALGEBRAIC(:,18) = CONSTANTS(:,605)./(1.00000+power(CONSTANTS(:,608)./STATES(:,19), CONSTANTS(:,606)));
RATES(:,75) = CONSTANTS(:,710).*(1.00000 - STATES(:,75)) - ALGEBRAIC(:,18).*STATES(:,75);
ALGEBRAIC(:,19) = CONSTANTS(:,666)./(1.00000+power(CONSTANTS(:,669)./STATES(:,20), CONSTANTS(:,667)));
RATES(:,80) = CONSTANTS(:,711).*(1.00000 - STATES(:,80)) - ALGEBRAIC(:,19).*STATES(:,80);
ALGEBRAIC(:,2) = 1.00000./(1.00000+exp( CONSTANTS(:,12).*(STATES(:,2) - CONSTANTS(:,13))));
ALGEBRAIC(:,20) = CONSTANTS(:,8)+(CONSTANTS(:,9) - CONSTANTS(:,8))./(1.00000+exp( CONSTANTS(:,10).*(STATES(:,2) - CONSTANTS(:,11))));
RATES(:,4) = (ALGEBRAIC(:,2) - STATES(:,4))./ALGEBRAIC(:,20);
ALGEBRAIC(:,4) = CONSTANTS(:,24).*(CONSTANTS(:,25)./(1.00000+exp( CONSTANTS(:,26).*(STATES(:,2) - CONSTANTS(:,27))))+(1.00000 - CONSTANTS(:,25)));
ALGEBRAIC(:,22) = CONSTANTS(:,24).*(CONSTANTS(:,25) - CONSTANTS(:,25)./(1.00000+exp( CONSTANTS(:,26).*(STATES(:,2) - CONSTANTS(:,27)))));
RATES(:,7) = ALGEBRAIC(:,4).*(1.00000 - STATES(:,7)) - ALGEBRAIC(:,22).*STATES(:,7);
ALGEBRAIC(:,5) = CONSTANTS(:,29).*(CONSTANTS(:,30)./(1.00000+exp( CONSTANTS(:,31).*(STATES(:,2) - CONSTANTS(:,32))))+(1.00000 - CONSTANTS(:,30)));
ALGEBRAIC(:,23) = CONSTANTS(:,29).*(CONSTANTS(:,30) - CONSTANTS(:,30)./(1.00000+exp( CONSTANTS(:,31).*(STATES(:,2) - CONSTANTS(:,32)))));
RATES(:,8) = ALGEBRAIC(:,5).*(1.00000 - STATES(:,8)) - ALGEBRAIC(:,23).*STATES(:,8);
ALGEBRAIC(:,3) = CONSTANTS(:,17)+(CONSTANTS(:,18) - CONSTANTS(:,17))./(1.00000+exp( CONSTANTS(:,19).*(STATES(:,2) - CONSTANTS(:,20))));
ALGEBRAIC(:,21) = ALGEBRAIC(:,3).*(CONSTANTS(:,21)./(1.00000+exp( CONSTANTS(:,22).*(STATES(:,2) - CONSTANTS(:,23))))+(1.00000 - CONSTANTS(:,21)));
ALGEBRAIC(:,27) = ALGEBRAIC(:,3).*(CONSTANTS(:,21) - CONSTANTS(:,21)./(1.00000+exp( CONSTANTS(:,22).*(STATES(:,2) - CONSTANTS(:,23)))));
RATES(:,6) = ALGEBRAIC(:,21).*(1.00000 - STATES(:,6)) - ALGEBRAIC(:,27).*STATES(:,6);
ALGEBRAIC(:,6) = CONSTANTS(:,35)+(CONSTANTS(:,36) - CONSTANTS(:,35))./(1.00000+exp( CONSTANTS(:,37).*(STATES(:,2) - CONSTANTS(:,38))));
ALGEBRAIC(:,24) = ALGEBRAIC(:,6).*(CONSTANTS(:,39)./(1.00000+exp( CONSTANTS(:,40).*(STATES(:,2) - CONSTANTS(:,41))))+(1.00000 - CONSTANTS(:,39)));
ALGEBRAIC(:,28) = ALGEBRAIC(:,6).*(CONSTANTS(:,39) - CONSTANTS(:,39)./(1.00000+exp( CONSTANTS(:,40).*(STATES(:,2) - CONSTANTS(:,41)))));
RATES(:,9) = ALGEBRAIC(:,24).*(1.00000 - STATES(:,9)) - ALGEBRAIC(:,28).*STATES(:,9);
ALGEBRAIC(:,7) = CONSTANTS(:,42)+(CONSTANTS(:,43) - CONSTANTS(:,42))./(1.00000+exp( CONSTANTS(:,44).*(STATES(:,2) - CONSTANTS(:,45))));
ALGEBRAIC(:,25) = ALGEBRAIC(:,7).*(CONSTANTS(:,46)./(1.00000+exp( CONSTANTS(:,47).*(STATES(:,2) - CONSTANTS(:,48))))+(1.00000 - CONSTANTS(:,46)));
ALGEBRAIC(:,29) = ALGEBRAIC(:,7).*(CONSTANTS(:,46) - CONSTANTS(:,46)./(1.00000+exp( CONSTANTS(:,47).*(STATES(:,2) - CONSTANTS(:,48)))));
RATES(:,10) = ALGEBRAIC(:,25).*(1.00000 - STATES(:,10)) - ALGEBRAIC(:,29).*STATES(:,10);
ALGEBRAIC(:,8) = CONSTANTS(:,1).*STATES(:,1).*(STATES(:,2) - CONSTANTS(:,2));
ALGEBRAIC(:,26) = CONSTANTS(:,14).*(STATES(:,5)./(CONSTANTS(:,15)+STATES(:,5)));
ALGEBRAIC(:,46) = CONSTANTS(:,75).*(CONSTANTS(:,76) - STATES(:,5));
ALGEBRAIC(:,36) = CONSTANTS(:,62).*(STATES(:,5) - STATES(:,11));
ALGEBRAIC(:,37) = CONSTANTS(:,63).*(STATES(:,5) - STATES(:,12));
ALGEBRAIC(:,38) = CONSTANTS(:,64).*(STATES(:,5) - STATES(:,13));
ALGEBRAIC(:,39) = CONSTANTS(:,65).*(STATES(:,5) - STATES(:,14));
ALGEBRAIC(:,40) = CONSTANTS(:,66).*(STATES(:,5) - STATES(:,15));
ALGEBRAIC(:,41) = CONSTANTS(:,67).*(STATES(:,5) - STATES(:,16));
ALGEBRAIC(:,42) = CONSTANTS(:,68).*(STATES(:,5) - STATES(:,17));
ALGEBRAIC(:,43) = CONSTANTS(:,69).*(STATES(:,5) - STATES(:,18));
ALGEBRAIC(:,44) = CONSTANTS(:,70).*(STATES(:,5) - STATES(:,19));
ALGEBRAIC(:,45) = CONSTANTS(:,71).*(STATES(:,5) - STATES(:,20));
RATES(:,5) = ALGEBRAIC(:,46) - ( CONSTANTS(:,60).*ALGEBRAIC(:,36)+ CONSTANTS(:,60).*ALGEBRAIC(:,37)+ CONSTANTS(:,60).*ALGEBRAIC(:,38)+ CONSTANTS(:,60).*ALGEBRAIC(:,39)+ CONSTANTS(:,60).*ALGEBRAIC(:,40)+ CONSTANTS(:,60).*ALGEBRAIC(:,41)+ CONSTANTS(:,60).*ALGEBRAIC(:,42)+ CONSTANTS(:,60).*ALGEBRAIC(:,43)+ CONSTANTS(:,60).*ALGEBRAIC(:,44)+ CONSTANTS(:,60).*ALGEBRAIC(:,45)+ (CONSTANTS(:,61)./2.00000).*(ALGEBRAIC(:,8)+ALGEBRAIC(:,26)));
ALGEBRAIC(:,60) = (( CONSTANTS(:,83).*CONSTANTS(:,82))./CONSTANTS(:,84)).*log(CONSTANTS(:,87)./STATES(:,31));
ALGEBRAIC(:,61) = CONSTANTS(:,86).*(STATES(:,2) - ALGEBRAIC(:,60));
ALGEBRAIC(:,69) = CONSTANTS(:,95).*(power(STATES(:,31), CONSTANTS(:,97))./(power(CONSTANTS(:,96), CONSTANTS(:,97))+power(STATES(:,31), CONSTANTS(:,97)))).*(CONSTANTS(:,98) - STATES(:,2));
ALGEBRAIC(:,64) = CONSTANTS(:,92).*(power(CONSTANTS(:,91), CONSTANTS(:,89))./(power(CONSTANTS(:,91), CONSTANTS(:,89))+power(STATES(:,21), CONSTANTS(:,89))));
ALGEBRAIC(:,65) = ALGEBRAIC(:,64).*(STATES(:,2) - CONSTANTS(:,90));
ALGEBRAIC(:,71) = CONSTANTS(:,716).*(ALGEBRAIC(:,65)+ALGEBRAIC(:,61)+ALGEBRAIC(:,69));
RATES(:,31) = - (CONSTANTS(:,129)./( 1.00000.*CONSTANTS(:,716))).*ALGEBRAIC(:,71);
ALGEBRAIC(:,47) = CONSTANTS(:,692).*(STATES(:,11) - STATES(:,21));
ALGEBRAIC(:,57) = ( CONSTANTS(:,78).*exp( CONSTANTS(:,79).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,80).*(STATES(:,2) - CONSTANTS(:,81))));
ALGEBRAIC(:,58) = (( CONSTANTS(:,83).*CONSTANTS(:,82))./( 2.00000.*CONSTANTS(:,84))).*log(CONSTANTS(:,85)./STATES(:,21));
ALGEBRAIC(:,59) = ALGEBRAIC(:,57).*(STATES(:,2) - ALGEBRAIC(:,58));
ALGEBRAIC(:,62) = CONSTANTS(:,88).*(power(CONSTANTS(:,91), CONSTANTS(:,89))./(power(CONSTANTS(:,91), CONSTANTS(:,89))+power(STATES(:,21), CONSTANTS(:,89))));
ALGEBRAIC(:,63) = ALGEBRAIC(:,62).*(STATES(:,2) - CONSTANTS(:,90));
ALGEBRAIC(:,66) = CONSTANTS(:,93).*(power(STATES(:,21), 2.00000)./(power(CONSTANTS(:,94), 2.00000)+power(STATES(:,21), 2.00000)));
ALGEBRAIC(:,67) = CONSTANTS(:,716).*(ALGEBRAIC(:,59)+ALGEBRAIC(:,63)+ALGEBRAIC(:,66));
ALGEBRAIC(:,68) = ( CONSTANTS(:,99).*(STATES(:,21) - CONSTANTS(:,100).*STATES(:,32)))./(1.00000+ CONSTANTS(:,101).*STATES(:,21)+ CONSTANTS(:,102).*STATES(:,32)+ CONSTANTS(:,103).*STATES(:,21).*STATES(:,32));
ALGEBRAIC(:,75) = CONSTANTS(:,127).*(STATES(:,33)./(STATES(:,33)+CONSTANTS(:,128)));
RATES(:,21) = (ALGEBRAIC(:,47)+ CONSTANTS(:,691).*ALGEBRAIC(:,75)) - ( (CONSTANTS(:,129)./( CONSTANTS(:,716).*2.00000)).*ALGEBRAIC(:,67)+ CONSTANTS(:,712).*ALGEBRAIC(:,68));
ALGEBRAIC(:,70) = power(CONSTANTS(:,106), CONSTANTS(:,107))./(power(CONSTANTS(:,106), CONSTANTS(:,107))+power(STATES(:,33), CONSTANTS(:,107)));
ALGEBRAIC(:,72) = CONSTANTS(:,104).*(power(STATES(:,11), 2.00000)./(power(CONSTANTS(:,105), 2.00000)+power(STATES(:,11), 2.00000))).*ALGEBRAIC(:,70);
ALGEBRAIC(:,78) = 1.00000./(1.00000+( CONSTANTS(:,134).*CONSTANTS(:,135))./power(CONSTANTS(:,134)+STATES(:,33), 2.00000));
RATES(:,33) = ALGEBRAIC(:,78).*(ALGEBRAIC(:,72) - ALGEBRAIC(:,75));
ALGEBRAIC(:,74) = ( CONSTANTS(:,110).*CONSTANTS(:,115)+ CONSTANTS(:,112).*STATES(:,11))./(CONSTANTS(:,115)+STATES(:,11));
ALGEBRAIC(:,77) = ( (CONSTANTS(:,109)+CONSTANTS(:,113)).*CONSTANTS(:,116))./(CONSTANTS(:,116)+STATES(:,11));
ALGEBRAIC(:,80) = CONSTANTS(:,108).*power(( CONSTANTS(:,126).*ALGEBRAIC(:,74).*STATES(:,34))./( CONSTANTS(:,126).*ALGEBRAIC(:,74)+ALGEBRAIC(:,77)), 4.00000).*(STATES(:,32) - STATES(:,11));
RATES(:,11) = (ALGEBRAIC(:,36)+ CONSTANTS(:,714).*ALGEBRAIC(:,80)) - ( CONSTANTS(:,715).*ALGEBRAIC(:,47)+ CONSTANTS(:,713).*ALGEBRAIC(:,72));
RATES(:,32) = ALGEBRAIC(:,68) - ALGEBRAIC(:,80);
ALGEBRAIC(:,81) = ( CONSTANTS(:,111).*CONSTANTS(:,116)+ CONSTANTS(:,114).*STATES(:,11))./(CONSTANTS(:,116)+STATES(:,11));
ALGEBRAIC(:,83) = ( CONSTANTS(:,125).*STATES(:,35))./( (1.00000+power(CONSTANTS(:,121)./STATES(:,11), CONSTANTS(:,123))).*(1.00000+power(STATES(:,11)./CONSTANTS(:,122), CONSTANTS(:,124))));
RATES(:,34) = ALGEBRAIC(:,83).*(1.00000 - STATES(:,34)) - (( CONSTANTS(:,126).*ALGEBRAIC(:,74).*ALGEBRAIC(:,81))./( CONSTANTS(:,126).*ALGEBRAIC(:,74)+ALGEBRAIC(:,77))).*STATES(:,34);
ALGEBRAIC(:,84) = (( CONSTANTS(:,144).*CONSTANTS(:,143))./CONSTANTS(:,145)).*log(CONSTANTS(:,148)./STATES(:,36));
ALGEBRAIC(:,85) = CONSTANTS(:,147).*(STATES(:,2) - ALGEBRAIC(:,84));
ALGEBRAIC(:,93) = CONSTANTS(:,156).*(power(STATES(:,36), CONSTANTS(:,158))./(power(CONSTANTS(:,157), CONSTANTS(:,158))+power(STATES(:,36), CONSTANTS(:,158)))).*(CONSTANTS(:,159) - STATES(:,2));
ALGEBRAIC(:,88) = CONSTANTS(:,153).*(power(CONSTANTS(:,152), CONSTANTS(:,150))./(power(CONSTANTS(:,152), CONSTANTS(:,150))+power(STATES(:,22), CONSTANTS(:,150))));
ALGEBRAIC(:,89) = ALGEBRAIC(:,88).*(STATES(:,2) - CONSTANTS(:,151));
ALGEBRAIC(:,95) = CONSTANTS(:,722).*(ALGEBRAIC(:,89)+ALGEBRAIC(:,85)+ALGEBRAIC(:,93));
RATES(:,36) = - (CONSTANTS(:,190)./( 1.00000.*CONSTANTS(:,722))).*ALGEBRAIC(:,95);
ALGEBRAIC(:,48) = CONSTANTS(:,693).*(STATES(:,12) - STATES(:,22));
ALGEBRAIC(:,76) = ( CONSTANTS(:,139).*exp( CONSTANTS(:,140).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,141).*(STATES(:,2) - CONSTANTS(:,142))));
ALGEBRAIC(:,79) = (( CONSTANTS(:,144).*CONSTANTS(:,143))./( 2.00000.*CONSTANTS(:,145))).*log(CONSTANTS(:,146)./STATES(:,22));
ALGEBRAIC(:,82) = ALGEBRAIC(:,76).*(STATES(:,2) - ALGEBRAIC(:,79));
ALGEBRAIC(:,86) = CONSTANTS(:,149).*(power(CONSTANTS(:,152), CONSTANTS(:,150))./(power(CONSTANTS(:,152), CONSTANTS(:,150))+power(STATES(:,22), CONSTANTS(:,150))));
ALGEBRAIC(:,87) = ALGEBRAIC(:,86).*(STATES(:,2) - CONSTANTS(:,151));
ALGEBRAIC(:,90) = CONSTANTS(:,154).*(power(STATES(:,22), 2.00000)./(power(CONSTANTS(:,155), 2.00000)+power(STATES(:,22), 2.00000)));
ALGEBRAIC(:,91) = CONSTANTS(:,722).*(ALGEBRAIC(:,82)+ALGEBRAIC(:,87)+ALGEBRAIC(:,90));
ALGEBRAIC(:,92) = ( CONSTANTS(:,160).*(STATES(:,22) - CONSTANTS(:,161).*STATES(:,37)))./(1.00000+ CONSTANTS(:,162).*STATES(:,22)+ CONSTANTS(:,163).*STATES(:,37)+ CONSTANTS(:,164).*STATES(:,22).*STATES(:,37));
ALGEBRAIC(:,99) = CONSTANTS(:,188).*(STATES(:,38)./(STATES(:,38)+CONSTANTS(:,189)));
RATES(:,22) = (ALGEBRAIC(:,48)+ CONSTANTS(:,717).*ALGEBRAIC(:,99)) - ( (CONSTANTS(:,190)./( CONSTANTS(:,722).*2.00000)).*ALGEBRAIC(:,91)+ CONSTANTS(:,718).*ALGEBRAIC(:,92));
ALGEBRAIC(:,94) = power(CONSTANTS(:,167), CONSTANTS(:,168))./(power(CONSTANTS(:,167), CONSTANTS(:,168))+power(STATES(:,38), CONSTANTS(:,168)));
ALGEBRAIC(:,96) = CONSTANTS(:,165).*(power(STATES(:,12), 2.00000)./(power(CONSTANTS(:,166), 2.00000)+power(STATES(:,12), 2.00000))).*ALGEBRAIC(:,94);
ALGEBRAIC(:,102) = 1.00000./(1.00000+( CONSTANTS(:,195).*CONSTANTS(:,196))./power(CONSTANTS(:,195)+STATES(:,38), 2.00000));
RATES(:,38) = ALGEBRAIC(:,102).*(ALGEBRAIC(:,96) - ALGEBRAIC(:,99));
ALGEBRAIC(:,98) = ( CONSTANTS(:,171).*CONSTANTS(:,176)+ CONSTANTS(:,173).*STATES(:,12))./(CONSTANTS(:,176)+STATES(:,12));
ALGEBRAIC(:,101) = ( (CONSTANTS(:,170)+CONSTANTS(:,174)).*CONSTANTS(:,177))./(CONSTANTS(:,177)+STATES(:,12));
ALGEBRAIC(:,104) = CONSTANTS(:,169).*power(( CONSTANTS(:,187).*ALGEBRAIC(:,98).*STATES(:,39))./( CONSTANTS(:,187).*ALGEBRAIC(:,98)+ALGEBRAIC(:,101)), 4.00000).*(STATES(:,37) - STATES(:,12));
RATES(:,12) = (ALGEBRAIC(:,37)+ CONSTANTS(:,720).*ALGEBRAIC(:,104)) - ( CONSTANTS(:,721).*ALGEBRAIC(:,48)+ CONSTANTS(:,719).*ALGEBRAIC(:,96));
RATES(:,37) = ALGEBRAIC(:,92) - ALGEBRAIC(:,104);
ALGEBRAIC(:,105) = ( CONSTANTS(:,172).*CONSTANTS(:,177)+ CONSTANTS(:,175).*STATES(:,12))./(CONSTANTS(:,177)+STATES(:,12));
ALGEBRAIC(:,107) = ( CONSTANTS(:,186).*STATES(:,40))./( (1.00000+power(CONSTANTS(:,182)./STATES(:,12), CONSTANTS(:,184))).*(1.00000+power(STATES(:,12)./CONSTANTS(:,183), CONSTANTS(:,185))));
RATES(:,39) = ALGEBRAIC(:,107).*(1.00000 - STATES(:,39)) - (( CONSTANTS(:,187).*ALGEBRAIC(:,98).*ALGEBRAIC(:,105))./( CONSTANTS(:,187).*ALGEBRAIC(:,98)+ALGEBRAIC(:,101))).*STATES(:,39);
ALGEBRAIC(:,108) = (( CONSTANTS(:,205).*CONSTANTS(:,204))./CONSTANTS(:,206)).*log(CONSTANTS(:,209)./STATES(:,41));
ALGEBRAIC(:,109) = CONSTANTS(:,208).*(STATES(:,2) - ALGEBRAIC(:,108));
ALGEBRAIC(:,117) = CONSTANTS(:,217).*(power(STATES(:,41), CONSTANTS(:,219))./(power(CONSTANTS(:,218), CONSTANTS(:,219))+power(STATES(:,41), CONSTANTS(:,219)))).*(CONSTANTS(:,220) - STATES(:,2));
ALGEBRAIC(:,112) = CONSTANTS(:,214).*(power(CONSTANTS(:,213), CONSTANTS(:,211))./(power(CONSTANTS(:,213), CONSTANTS(:,211))+power(STATES(:,23), CONSTANTS(:,211))));
ALGEBRAIC(:,113) = ALGEBRAIC(:,112).*(STATES(:,2) - CONSTANTS(:,212));
ALGEBRAIC(:,119) = CONSTANTS(:,728).*(ALGEBRAIC(:,113)+ALGEBRAIC(:,109)+ALGEBRAIC(:,117));
RATES(:,41) = - (CONSTANTS(:,251)./( 1.00000.*CONSTANTS(:,728))).*ALGEBRAIC(:,119);
ALGEBRAIC(:,49) = CONSTANTS(:,694).*(STATES(:,13) - STATES(:,23));
ALGEBRAIC(:,100) = ( CONSTANTS(:,200).*exp( CONSTANTS(:,201).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,202).*(STATES(:,2) - CONSTANTS(:,203))));
ALGEBRAIC(:,103) = (( CONSTANTS(:,205).*CONSTANTS(:,204))./( 2.00000.*CONSTANTS(:,206))).*log(CONSTANTS(:,207)./STATES(:,23));
ALGEBRAIC(:,106) = ALGEBRAIC(:,100).*(STATES(:,2) - ALGEBRAIC(:,103));
ALGEBRAIC(:,110) = CONSTANTS(:,210).*(power(CONSTANTS(:,213), CONSTANTS(:,211))./(power(CONSTANTS(:,213), CONSTANTS(:,211))+power(STATES(:,23), CONSTANTS(:,211))));
ALGEBRAIC(:,111) = ALGEBRAIC(:,110).*(STATES(:,2) - CONSTANTS(:,212));
ALGEBRAIC(:,114) = CONSTANTS(:,215).*(power(STATES(:,23), 2.00000)./(power(CONSTANTS(:,216), 2.00000)+power(STATES(:,23), 2.00000)));
ALGEBRAIC(:,115) = CONSTANTS(:,728).*(ALGEBRAIC(:,106)+ALGEBRAIC(:,111)+ALGEBRAIC(:,114));
ALGEBRAIC(:,116) = ( CONSTANTS(:,221).*(STATES(:,23) - CONSTANTS(:,222).*STATES(:,42)))./(1.00000+ CONSTANTS(:,223).*STATES(:,23)+ CONSTANTS(:,224).*STATES(:,42)+ CONSTANTS(:,225).*STATES(:,23).*STATES(:,42));
ALGEBRAIC(:,123) = CONSTANTS(:,249).*(STATES(:,43)./(STATES(:,43)+CONSTANTS(:,250)));
RATES(:,23) = (ALGEBRAIC(:,49)+ CONSTANTS(:,723).*ALGEBRAIC(:,123)) - ( (CONSTANTS(:,251)./( CONSTANTS(:,728).*2.00000)).*ALGEBRAIC(:,115)+ CONSTANTS(:,724).*ALGEBRAIC(:,116));
ALGEBRAIC(:,118) = power(CONSTANTS(:,228), CONSTANTS(:,229))./(power(CONSTANTS(:,228), CONSTANTS(:,229))+power(STATES(:,43), CONSTANTS(:,229)));
ALGEBRAIC(:,120) = CONSTANTS(:,226).*(power(STATES(:,13), 2.00000)./(power(CONSTANTS(:,227), 2.00000)+power(STATES(:,13), 2.00000))).*ALGEBRAIC(:,118);
ALGEBRAIC(:,126) = 1.00000./(1.00000+( CONSTANTS(:,256).*CONSTANTS(:,257))./power(CONSTANTS(:,256)+STATES(:,43), 2.00000));
RATES(:,43) = ALGEBRAIC(:,126).*(ALGEBRAIC(:,120) - ALGEBRAIC(:,123));
ALGEBRAIC(:,122) = ( CONSTANTS(:,232).*CONSTANTS(:,237)+ CONSTANTS(:,234).*STATES(:,13))./(CONSTANTS(:,237)+STATES(:,13));
ALGEBRAIC(:,125) = ( (CONSTANTS(:,231)+CONSTANTS(:,235)).*CONSTANTS(:,238))./(CONSTANTS(:,238)+STATES(:,13));
ALGEBRAIC(:,128) = CONSTANTS(:,230).*power(( CONSTANTS(:,248).*ALGEBRAIC(:,122).*STATES(:,44))./( CONSTANTS(:,248).*ALGEBRAIC(:,122)+ALGEBRAIC(:,125)), 4.00000).*(STATES(:,42) - STATES(:,13));
RATES(:,13) = (ALGEBRAIC(:,38)+ CONSTANTS(:,726).*ALGEBRAIC(:,128)) - ( CONSTANTS(:,727).*ALGEBRAIC(:,49)+ CONSTANTS(:,725).*ALGEBRAIC(:,120));
RATES(:,42) = ALGEBRAIC(:,116) - ALGEBRAIC(:,128);
ALGEBRAIC(:,129) = ( CONSTANTS(:,233).*CONSTANTS(:,238)+ CONSTANTS(:,236).*STATES(:,13))./(CONSTANTS(:,238)+STATES(:,13));
ALGEBRAIC(:,131) = ( CONSTANTS(:,247).*STATES(:,45))./( (1.00000+power(CONSTANTS(:,243)./STATES(:,13), CONSTANTS(:,245))).*(1.00000+power(STATES(:,13)./CONSTANTS(:,244), CONSTANTS(:,246))));
RATES(:,44) = ALGEBRAIC(:,131).*(1.00000 - STATES(:,44)) - (( CONSTANTS(:,248).*ALGEBRAIC(:,122).*ALGEBRAIC(:,129))./( CONSTANTS(:,248).*ALGEBRAIC(:,122)+ALGEBRAIC(:,125))).*STATES(:,44);
ALGEBRAIC(:,132) = (( CONSTANTS(:,266).*CONSTANTS(:,265))./CONSTANTS(:,267)).*log(CONSTANTS(:,270)./STATES(:,46));
ALGEBRAIC(:,133) = CONSTANTS(:,269).*(STATES(:,2) - ALGEBRAIC(:,132));
ALGEBRAIC(:,141) = CONSTANTS(:,278).*(power(STATES(:,46), CONSTANTS(:,280))./(power(CONSTANTS(:,279), CONSTANTS(:,280))+power(STATES(:,46), CONSTANTS(:,280)))).*(CONSTANTS(:,281) - STATES(:,2));
ALGEBRAIC(:,136) = CONSTANTS(:,275).*(power(CONSTANTS(:,274), CONSTANTS(:,272))./(power(CONSTANTS(:,274), CONSTANTS(:,272))+power(STATES(:,24), CONSTANTS(:,272))));
ALGEBRAIC(:,137) = ALGEBRAIC(:,136).*(STATES(:,2) - CONSTANTS(:,273));
ALGEBRAIC(:,143) = CONSTANTS(:,734).*(ALGEBRAIC(:,137)+ALGEBRAIC(:,133)+ALGEBRAIC(:,141));
RATES(:,46) = - (CONSTANTS(:,312)./( 1.00000.*CONSTANTS(:,734))).*ALGEBRAIC(:,143);
ALGEBRAIC(:,50) = CONSTANTS(:,695).*(STATES(:,14) - STATES(:,24));
ALGEBRAIC(:,124) = ( CONSTANTS(:,261).*exp( CONSTANTS(:,262).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,263).*(STATES(:,2) - CONSTANTS(:,264))));
ALGEBRAIC(:,127) = (( CONSTANTS(:,266).*CONSTANTS(:,265))./( 2.00000.*CONSTANTS(:,267))).*log(CONSTANTS(:,268)./STATES(:,24));
ALGEBRAIC(:,130) = ALGEBRAIC(:,124).*(STATES(:,2) - ALGEBRAIC(:,127));
ALGEBRAIC(:,134) = CONSTANTS(:,271).*(power(CONSTANTS(:,274), CONSTANTS(:,272))./(power(CONSTANTS(:,274), CONSTANTS(:,272))+power(STATES(:,24), CONSTANTS(:,272))));
ALGEBRAIC(:,135) = ALGEBRAIC(:,134).*(STATES(:,2) - CONSTANTS(:,273));
ALGEBRAIC(:,138) = CONSTANTS(:,276).*(power(STATES(:,24), 2.00000)./(power(CONSTANTS(:,277), 2.00000)+power(STATES(:,24), 2.00000)));
ALGEBRAIC(:,139) = CONSTANTS(:,734).*(ALGEBRAIC(:,130)+ALGEBRAIC(:,135)+ALGEBRAIC(:,138));
ALGEBRAIC(:,140) = ( CONSTANTS(:,282).*(STATES(:,24) - CONSTANTS(:,283).*STATES(:,47)))./(1.00000+ CONSTANTS(:,284).*STATES(:,24)+ CONSTANTS(:,285).*STATES(:,47)+ CONSTANTS(:,286).*STATES(:,24).*STATES(:,47));
ALGEBRAIC(:,147) = CONSTANTS(:,310).*(STATES(:,48)./(STATES(:,48)+CONSTANTS(:,311)));
RATES(:,24) = (ALGEBRAIC(:,50)+ CONSTANTS(:,729).*ALGEBRAIC(:,147)) - ( (CONSTANTS(:,312)./( CONSTANTS(:,734).*2.00000)).*ALGEBRAIC(:,139)+ CONSTANTS(:,730).*ALGEBRAIC(:,140));
ALGEBRAIC(:,142) = power(CONSTANTS(:,289), CONSTANTS(:,290))./(power(CONSTANTS(:,289), CONSTANTS(:,290))+power(STATES(:,48), CONSTANTS(:,290)));
ALGEBRAIC(:,144) = CONSTANTS(:,287).*(power(STATES(:,14), 2.00000)./(power(CONSTANTS(:,288), 2.00000)+power(STATES(:,14), 2.00000))).*ALGEBRAIC(:,142);
ALGEBRAIC(:,150) = 1.00000./(1.00000+( CONSTANTS(:,317).*CONSTANTS(:,318))./power(CONSTANTS(:,317)+STATES(:,48), 2.00000));
RATES(:,48) = ALGEBRAIC(:,150).*(ALGEBRAIC(:,144) - ALGEBRAIC(:,147));
ALGEBRAIC(:,146) = ( CONSTANTS(:,293).*CONSTANTS(:,298)+ CONSTANTS(:,295).*STATES(:,14))./(CONSTANTS(:,298)+STATES(:,14));
ALGEBRAIC(:,149) = ( (CONSTANTS(:,292)+CONSTANTS(:,296)).*CONSTANTS(:,299))./(CONSTANTS(:,299)+STATES(:,14));
ALGEBRAIC(:,152) = CONSTANTS(:,291).*power(( CONSTANTS(:,309).*ALGEBRAIC(:,146).*STATES(:,49))./( CONSTANTS(:,309).*ALGEBRAIC(:,146)+ALGEBRAIC(:,149)), 4.00000).*(STATES(:,47) - STATES(:,14));
RATES(:,14) = (ALGEBRAIC(:,39)+ CONSTANTS(:,732).*ALGEBRAIC(:,152)) - ( CONSTANTS(:,733).*ALGEBRAIC(:,50)+ CONSTANTS(:,731).*ALGEBRAIC(:,144));
RATES(:,47) = ALGEBRAIC(:,140) - ALGEBRAIC(:,152);
ALGEBRAIC(:,153) = ( CONSTANTS(:,294).*CONSTANTS(:,299)+ CONSTANTS(:,297).*STATES(:,14))./(CONSTANTS(:,299)+STATES(:,14));
ALGEBRAIC(:,155) = ( CONSTANTS(:,308).*STATES(:,50))./( (1.00000+power(CONSTANTS(:,304)./STATES(:,14), CONSTANTS(:,306))).*(1.00000+power(STATES(:,14)./CONSTANTS(:,305), CONSTANTS(:,307))));
RATES(:,49) = ALGEBRAIC(:,155).*(1.00000 - STATES(:,49)) - (( CONSTANTS(:,309).*ALGEBRAIC(:,146).*ALGEBRAIC(:,153))./( CONSTANTS(:,309).*ALGEBRAIC(:,146)+ALGEBRAIC(:,149))).*STATES(:,49);
ALGEBRAIC(:,156) = (( CONSTANTS(:,327).*CONSTANTS(:,326))./CONSTANTS(:,328)).*log(CONSTANTS(:,331)./STATES(:,51));
ALGEBRAIC(:,157) = CONSTANTS(:,330).*(STATES(:,2) - ALGEBRAIC(:,156));
ALGEBRAIC(:,165) = CONSTANTS(:,339).*(power(STATES(:,51), CONSTANTS(:,341))./(power(CONSTANTS(:,340), CONSTANTS(:,341))+power(STATES(:,51), CONSTANTS(:,341)))).*(CONSTANTS(:,342) - STATES(:,2));
ALGEBRAIC(:,160) = CONSTANTS(:,336).*(power(CONSTANTS(:,335), CONSTANTS(:,333))./(power(CONSTANTS(:,335), CONSTANTS(:,333))+power(STATES(:,25), CONSTANTS(:,333))));
ALGEBRAIC(:,161) = ALGEBRAIC(:,160).*(STATES(:,2) - CONSTANTS(:,334));
ALGEBRAIC(:,167) = CONSTANTS(:,740).*(ALGEBRAIC(:,161)+ALGEBRAIC(:,157)+ALGEBRAIC(:,165));
RATES(:,51) = - (CONSTANTS(:,373)./( 1.00000.*CONSTANTS(:,740))).*ALGEBRAIC(:,167);
ALGEBRAIC(:,51) = CONSTANTS(:,696).*(STATES(:,15) - STATES(:,25));
ALGEBRAIC(:,148) = ( CONSTANTS(:,322).*exp( CONSTANTS(:,323).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,324).*(STATES(:,2) - CONSTANTS(:,325))));
ALGEBRAIC(:,151) = (( CONSTANTS(:,327).*CONSTANTS(:,326))./( 2.00000.*CONSTANTS(:,328))).*log(CONSTANTS(:,329)./STATES(:,25));
ALGEBRAIC(:,154) = ALGEBRAIC(:,148).*(STATES(:,2) - ALGEBRAIC(:,151));
ALGEBRAIC(:,158) = CONSTANTS(:,332).*(power(CONSTANTS(:,335), CONSTANTS(:,333))./(power(CONSTANTS(:,335), CONSTANTS(:,333))+power(STATES(:,25), CONSTANTS(:,333))));
ALGEBRAIC(:,159) = ALGEBRAIC(:,158).*(STATES(:,2) - CONSTANTS(:,334));
ALGEBRAIC(:,162) = CONSTANTS(:,337).*(power(STATES(:,25), 2.00000)./(power(CONSTANTS(:,338), 2.00000)+power(STATES(:,25), 2.00000)));
ALGEBRAIC(:,163) = CONSTANTS(:,740).*(ALGEBRAIC(:,154)+ALGEBRAIC(:,159)+ALGEBRAIC(:,162));
ALGEBRAIC(:,164) = ( CONSTANTS(:,343).*(STATES(:,25) - CONSTANTS(:,344).*STATES(:,52)))./(1.00000+ CONSTANTS(:,345).*STATES(:,25)+ CONSTANTS(:,346).*STATES(:,52)+ CONSTANTS(:,347).*STATES(:,25).*STATES(:,52));
ALGEBRAIC(:,171) = CONSTANTS(:,371).*(STATES(:,53)./(STATES(:,53)+CONSTANTS(:,372)));
RATES(:,25) = (ALGEBRAIC(:,51)+ CONSTANTS(:,735).*ALGEBRAIC(:,171)) - ( (CONSTANTS(:,373)./( CONSTANTS(:,740).*2.00000)).*ALGEBRAIC(:,163)+ CONSTANTS(:,736).*ALGEBRAIC(:,164));
ALGEBRAIC(:,166) = power(CONSTANTS(:,350), CONSTANTS(:,351))./(power(CONSTANTS(:,350), CONSTANTS(:,351))+power(STATES(:,53), CONSTANTS(:,351)));
ALGEBRAIC(:,168) = CONSTANTS(:,348).*(power(STATES(:,15), 2.00000)./(power(CONSTANTS(:,349), 2.00000)+power(STATES(:,15), 2.00000))).*ALGEBRAIC(:,166);
ALGEBRAIC(:,174) = 1.00000./(1.00000+( CONSTANTS(:,378).*CONSTANTS(:,379))./power(CONSTANTS(:,378)+STATES(:,53), 2.00000));
RATES(:,53) = ALGEBRAIC(:,174).*(ALGEBRAIC(:,168) - ALGEBRAIC(:,171));
ALGEBRAIC(:,170) = ( CONSTANTS(:,354).*CONSTANTS(:,359)+ CONSTANTS(:,356).*STATES(:,15))./(CONSTANTS(:,359)+STATES(:,15));
ALGEBRAIC(:,173) = ( (CONSTANTS(:,353)+CONSTANTS(:,357)).*CONSTANTS(:,360))./(CONSTANTS(:,360)+STATES(:,15));
ALGEBRAIC(:,176) = CONSTANTS(:,352).*power(( CONSTANTS(:,370).*ALGEBRAIC(:,170).*STATES(:,54))./( CONSTANTS(:,370).*ALGEBRAIC(:,170)+ALGEBRAIC(:,173)), 4.00000).*(STATES(:,52) - STATES(:,15));
RATES(:,15) = (ALGEBRAIC(:,40)+ CONSTANTS(:,738).*ALGEBRAIC(:,176)) - ( CONSTANTS(:,739).*ALGEBRAIC(:,51)+ CONSTANTS(:,737).*ALGEBRAIC(:,168));
RATES(:,52) = ALGEBRAIC(:,164) - ALGEBRAIC(:,176);
ALGEBRAIC(:,177) = ( CONSTANTS(:,355).*CONSTANTS(:,360)+ CONSTANTS(:,358).*STATES(:,15))./(CONSTANTS(:,360)+STATES(:,15));
ALGEBRAIC(:,179) = ( CONSTANTS(:,369).*STATES(:,55))./( (1.00000+power(CONSTANTS(:,365)./STATES(:,15), CONSTANTS(:,367))).*(1.00000+power(STATES(:,15)./CONSTANTS(:,366), CONSTANTS(:,368))));
RATES(:,54) = ALGEBRAIC(:,179).*(1.00000 - STATES(:,54)) - (( CONSTANTS(:,370).*ALGEBRAIC(:,170).*ALGEBRAIC(:,177))./( CONSTANTS(:,370).*ALGEBRAIC(:,170)+ALGEBRAIC(:,173))).*STATES(:,54);
ALGEBRAIC(:,180) = (( CONSTANTS(:,388).*CONSTANTS(:,387))./CONSTANTS(:,389)).*log(CONSTANTS(:,392)./STATES(:,56));
ALGEBRAIC(:,181) = CONSTANTS(:,391).*(STATES(:,2) - ALGEBRAIC(:,180));
ALGEBRAIC(:,189) = CONSTANTS(:,400).*(power(STATES(:,56), CONSTANTS(:,402))./(power(CONSTANTS(:,401), CONSTANTS(:,402))+power(STATES(:,56), CONSTANTS(:,402)))).*(CONSTANTS(:,403) - STATES(:,2));
ALGEBRAIC(:,184) = CONSTANTS(:,397).*(power(CONSTANTS(:,396), CONSTANTS(:,394))./(power(CONSTANTS(:,396), CONSTANTS(:,394))+power(STATES(:,26), CONSTANTS(:,394))));
ALGEBRAIC(:,185) = ALGEBRAIC(:,184).*(STATES(:,2) - CONSTANTS(:,395));
ALGEBRAIC(:,191) = CONSTANTS(:,746).*(ALGEBRAIC(:,185)+ALGEBRAIC(:,181)+ALGEBRAIC(:,189));
RATES(:,56) = - (CONSTANTS(:,434)./( 1.00000.*CONSTANTS(:,746))).*ALGEBRAIC(:,191);
ALGEBRAIC(:,52) = CONSTANTS(:,697).*(STATES(:,16) - STATES(:,26));
ALGEBRAIC(:,172) = ( CONSTANTS(:,383).*exp( CONSTANTS(:,384).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,385).*(STATES(:,2) - CONSTANTS(:,386))));
ALGEBRAIC(:,175) = (( CONSTANTS(:,388).*CONSTANTS(:,387))./( 2.00000.*CONSTANTS(:,389))).*log(CONSTANTS(:,390)./STATES(:,26));
ALGEBRAIC(:,178) = ALGEBRAIC(:,172).*(STATES(:,2) - ALGEBRAIC(:,175));
ALGEBRAIC(:,182) = CONSTANTS(:,393).*(power(CONSTANTS(:,396), CONSTANTS(:,394))./(power(CONSTANTS(:,396), CONSTANTS(:,394))+power(STATES(:,26), CONSTANTS(:,394))));
ALGEBRAIC(:,183) = ALGEBRAIC(:,182).*(STATES(:,2) - CONSTANTS(:,395));
ALGEBRAIC(:,186) = CONSTANTS(:,398).*(power(STATES(:,26), 2.00000)./(power(CONSTANTS(:,399), 2.00000)+power(STATES(:,26), 2.00000)));
ALGEBRAIC(:,187) = CONSTANTS(:,746).*(ALGEBRAIC(:,178)+ALGEBRAIC(:,183)+ALGEBRAIC(:,186));
ALGEBRAIC(:,188) = ( CONSTANTS(:,404).*(STATES(:,26) - CONSTANTS(:,405).*STATES(:,57)))./(1.00000+ CONSTANTS(:,406).*STATES(:,26)+ CONSTANTS(:,407).*STATES(:,57)+ CONSTANTS(:,408).*STATES(:,26).*STATES(:,57));
ALGEBRAIC(:,195) = CONSTANTS(:,432).*(STATES(:,58)./(STATES(:,58)+CONSTANTS(:,433)));
RATES(:,26) = (ALGEBRAIC(:,52)+ CONSTANTS(:,741).*ALGEBRAIC(:,195)) - ( (CONSTANTS(:,434)./( CONSTANTS(:,746).*2.00000)).*ALGEBRAIC(:,187)+ CONSTANTS(:,742).*ALGEBRAIC(:,188));
ALGEBRAIC(:,190) = power(CONSTANTS(:,411), CONSTANTS(:,412))./(power(CONSTANTS(:,411), CONSTANTS(:,412))+power(STATES(:,58), CONSTANTS(:,412)));
ALGEBRAIC(:,192) = CONSTANTS(:,409).*(power(STATES(:,16), 2.00000)./(power(CONSTANTS(:,410), 2.00000)+power(STATES(:,16), 2.00000))).*ALGEBRAIC(:,190);
ALGEBRAIC(:,198) = 1.00000./(1.00000+( CONSTANTS(:,439).*CONSTANTS(:,440))./power(CONSTANTS(:,439)+STATES(:,58), 2.00000));
RATES(:,58) = ALGEBRAIC(:,198).*(ALGEBRAIC(:,192) - ALGEBRAIC(:,195));
ALGEBRAIC(:,194) = ( CONSTANTS(:,415).*CONSTANTS(:,420)+ CONSTANTS(:,417).*STATES(:,16))./(CONSTANTS(:,420)+STATES(:,16));
ALGEBRAIC(:,197) = ( (CONSTANTS(:,414)+CONSTANTS(:,418)).*CONSTANTS(:,421))./(CONSTANTS(:,421)+STATES(:,16));
ALGEBRAIC(:,200) = CONSTANTS(:,413).*power(( CONSTANTS(:,431).*ALGEBRAIC(:,194).*STATES(:,59))./( CONSTANTS(:,431).*ALGEBRAIC(:,194)+ALGEBRAIC(:,197)), 4.00000).*(STATES(:,57) - STATES(:,16));
RATES(:,16) = (ALGEBRAIC(:,41)+ CONSTANTS(:,744).*ALGEBRAIC(:,200)) - ( CONSTANTS(:,745).*ALGEBRAIC(:,52)+ CONSTANTS(:,743).*ALGEBRAIC(:,192));
RATES(:,57) = ALGEBRAIC(:,188) - ALGEBRAIC(:,200);
ALGEBRAIC(:,201) = ( CONSTANTS(:,416).*CONSTANTS(:,421)+ CONSTANTS(:,419).*STATES(:,16))./(CONSTANTS(:,421)+STATES(:,16));
ALGEBRAIC(:,203) = ( CONSTANTS(:,430).*STATES(:,60))./( (1.00000+power(CONSTANTS(:,426)./STATES(:,16), CONSTANTS(:,428))).*(1.00000+power(STATES(:,16)./CONSTANTS(:,427), CONSTANTS(:,429))));
RATES(:,59) = ALGEBRAIC(:,203).*(1.00000 - STATES(:,59)) - (( CONSTANTS(:,431).*ALGEBRAIC(:,194).*ALGEBRAIC(:,201))./( CONSTANTS(:,431).*ALGEBRAIC(:,194)+ALGEBRAIC(:,197))).*STATES(:,59);
ALGEBRAIC(:,204) = (( CONSTANTS(:,449).*CONSTANTS(:,448))./CONSTANTS(:,450)).*log(CONSTANTS(:,453)./STATES(:,61));
ALGEBRAIC(:,205) = CONSTANTS(:,452).*(STATES(:,2) - ALGEBRAIC(:,204));
ALGEBRAIC(:,213) = CONSTANTS(:,461).*(power(STATES(:,61), CONSTANTS(:,463))./(power(CONSTANTS(:,462), CONSTANTS(:,463))+power(STATES(:,61), CONSTANTS(:,463)))).*(CONSTANTS(:,464) - STATES(:,2));
ALGEBRAIC(:,208) = CONSTANTS(:,458).*(power(CONSTANTS(:,457), CONSTANTS(:,455))./(power(CONSTANTS(:,457), CONSTANTS(:,455))+power(STATES(:,27), CONSTANTS(:,455))));
ALGEBRAIC(:,209) = ALGEBRAIC(:,208).*(STATES(:,2) - CONSTANTS(:,456));
ALGEBRAIC(:,215) = CONSTANTS(:,752).*(ALGEBRAIC(:,209)+ALGEBRAIC(:,205)+ALGEBRAIC(:,213));
RATES(:,61) = - (CONSTANTS(:,495)./( 1.00000.*CONSTANTS(:,752))).*ALGEBRAIC(:,215);
ALGEBRAIC(:,53) = CONSTANTS(:,698).*(STATES(:,17) - STATES(:,27));
ALGEBRAIC(:,196) = ( CONSTANTS(:,444).*exp( CONSTANTS(:,445).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,446).*(STATES(:,2) - CONSTANTS(:,447))));
ALGEBRAIC(:,199) = (( CONSTANTS(:,449).*CONSTANTS(:,448))./( 2.00000.*CONSTANTS(:,450))).*log(CONSTANTS(:,451)./STATES(:,27));
ALGEBRAIC(:,202) = ALGEBRAIC(:,196).*(STATES(:,2) - ALGEBRAIC(:,199));
ALGEBRAIC(:,206) = CONSTANTS(:,454).*(power(CONSTANTS(:,457), CONSTANTS(:,455))./(power(CONSTANTS(:,457), CONSTANTS(:,455))+power(STATES(:,27), CONSTANTS(:,455))));
ALGEBRAIC(:,207) = ALGEBRAIC(:,206).*(STATES(:,2) - CONSTANTS(:,456));
ALGEBRAIC(:,210) = CONSTANTS(:,459).*(power(STATES(:,27), 2.00000)./(power(CONSTANTS(:,460), 2.00000)+power(STATES(:,27), 2.00000)));
ALGEBRAIC(:,211) = CONSTANTS(:,752).*(ALGEBRAIC(:,202)+ALGEBRAIC(:,207)+ALGEBRAIC(:,210));
ALGEBRAIC(:,212) = ( CONSTANTS(:,465).*(STATES(:,27) - CONSTANTS(:,466).*STATES(:,62)))./(1.00000+ CONSTANTS(:,467).*STATES(:,27)+ CONSTANTS(:,468).*STATES(:,62)+ CONSTANTS(:,469).*STATES(:,27).*STATES(:,62));
ALGEBRAIC(:,219) = CONSTANTS(:,493).*(STATES(:,63)./(STATES(:,63)+CONSTANTS(:,494)));
RATES(:,27) = (ALGEBRAIC(:,53)+ CONSTANTS(:,747).*ALGEBRAIC(:,219)) - ( (CONSTANTS(:,495)./( CONSTANTS(:,752).*2.00000)).*ALGEBRAIC(:,211)+ CONSTANTS(:,748).*ALGEBRAIC(:,212));
ALGEBRAIC(:,214) = power(CONSTANTS(:,472), CONSTANTS(:,473))./(power(CONSTANTS(:,472), CONSTANTS(:,473))+power(STATES(:,63), CONSTANTS(:,473)));
ALGEBRAIC(:,216) = CONSTANTS(:,470).*(power(STATES(:,17), 2.00000)./(power(CONSTANTS(:,471), 2.00000)+power(STATES(:,17), 2.00000))).*ALGEBRAIC(:,214);
ALGEBRAIC(:,222) = 1.00000./(1.00000+( CONSTANTS(:,500).*CONSTANTS(:,501))./power(CONSTANTS(:,500)+STATES(:,63), 2.00000));
RATES(:,63) = ALGEBRAIC(:,222).*(ALGEBRAIC(:,216) - ALGEBRAIC(:,219));
ALGEBRAIC(:,218) = ( CONSTANTS(:,476).*CONSTANTS(:,481)+ CONSTANTS(:,478).*STATES(:,17))./(CONSTANTS(:,481)+STATES(:,17));
ALGEBRAIC(:,221) = ( (CONSTANTS(:,475)+CONSTANTS(:,479)).*CONSTANTS(:,482))./(CONSTANTS(:,482)+STATES(:,17));
ALGEBRAIC(:,224) = CONSTANTS(:,474).*power(( CONSTANTS(:,492).*ALGEBRAIC(:,218).*STATES(:,64))./( CONSTANTS(:,492).*ALGEBRAIC(:,218)+ALGEBRAIC(:,221)), 4.00000).*(STATES(:,62) - STATES(:,17));
RATES(:,17) = (ALGEBRAIC(:,42)+ CONSTANTS(:,750).*ALGEBRAIC(:,224)) - ( CONSTANTS(:,751).*ALGEBRAIC(:,53)+ CONSTANTS(:,749).*ALGEBRAIC(:,216));
RATES(:,62) = ALGEBRAIC(:,212) - ALGEBRAIC(:,224);
ALGEBRAIC(:,225) = ( CONSTANTS(:,477).*CONSTANTS(:,482)+ CONSTANTS(:,480).*STATES(:,17))./(CONSTANTS(:,482)+STATES(:,17));
ALGEBRAIC(:,227) = ( CONSTANTS(:,491).*STATES(:,65))./( (1.00000+power(CONSTANTS(:,487)./STATES(:,17), CONSTANTS(:,489))).*(1.00000+power(STATES(:,17)./CONSTANTS(:,488), CONSTANTS(:,490))));
RATES(:,64) = ALGEBRAIC(:,227).*(1.00000 - STATES(:,64)) - (( CONSTANTS(:,492).*ALGEBRAIC(:,218).*ALGEBRAIC(:,225))./( CONSTANTS(:,492).*ALGEBRAIC(:,218)+ALGEBRAIC(:,221))).*STATES(:,64);
ALGEBRAIC(:,228) = (( CONSTANTS(:,510).*CONSTANTS(:,509))./CONSTANTS(:,511)).*log(CONSTANTS(:,514)./STATES(:,66));
ALGEBRAIC(:,229) = CONSTANTS(:,513).*(STATES(:,2) - ALGEBRAIC(:,228));
ALGEBRAIC(:,237) = CONSTANTS(:,522).*(power(STATES(:,66), CONSTANTS(:,524))./(power(CONSTANTS(:,523), CONSTANTS(:,524))+power(STATES(:,66), CONSTANTS(:,524)))).*(CONSTANTS(:,525) - STATES(:,2));
ALGEBRAIC(:,232) = CONSTANTS(:,519).*(power(CONSTANTS(:,518), CONSTANTS(:,516))./(power(CONSTANTS(:,518), CONSTANTS(:,516))+power(STATES(:,28), CONSTANTS(:,516))));
ALGEBRAIC(:,233) = ALGEBRAIC(:,232).*(STATES(:,2) - CONSTANTS(:,517));
ALGEBRAIC(:,239) = CONSTANTS(:,758).*(ALGEBRAIC(:,233)+ALGEBRAIC(:,229)+ALGEBRAIC(:,237));
RATES(:,66) = - (CONSTANTS(:,556)./( 1.00000.*CONSTANTS(:,758))).*ALGEBRAIC(:,239);
ALGEBRAIC(:,54) = CONSTANTS(:,699).*(STATES(:,18) - STATES(:,28));
ALGEBRAIC(:,220) = ( CONSTANTS(:,505).*exp( CONSTANTS(:,506).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,507).*(STATES(:,2) - CONSTANTS(:,508))));
ALGEBRAIC(:,223) = (( CONSTANTS(:,510).*CONSTANTS(:,509))./( 2.00000.*CONSTANTS(:,511))).*log(CONSTANTS(:,512)./STATES(:,28));
ALGEBRAIC(:,226) = ALGEBRAIC(:,220).*(STATES(:,2) - ALGEBRAIC(:,223));
ALGEBRAIC(:,230) = CONSTANTS(:,515).*(power(CONSTANTS(:,518), CONSTANTS(:,516))./(power(CONSTANTS(:,518), CONSTANTS(:,516))+power(STATES(:,28), CONSTANTS(:,516))));
ALGEBRAIC(:,231) = ALGEBRAIC(:,230).*(STATES(:,2) - CONSTANTS(:,517));
ALGEBRAIC(:,234) = CONSTANTS(:,520).*(power(STATES(:,28), 2.00000)./(power(CONSTANTS(:,521), 2.00000)+power(STATES(:,28), 2.00000)));
ALGEBRAIC(:,235) = CONSTANTS(:,758).*(ALGEBRAIC(:,226)+ALGEBRAIC(:,231)+ALGEBRAIC(:,234));
ALGEBRAIC(:,236) = ( CONSTANTS(:,526).*(STATES(:,28) - CONSTANTS(:,527).*STATES(:,67)))./(1.00000+ CONSTANTS(:,528).*STATES(:,28)+ CONSTANTS(:,529).*STATES(:,67)+ CONSTANTS(:,530).*STATES(:,28).*STATES(:,67));
ALGEBRAIC(:,243) = CONSTANTS(:,554).*(STATES(:,68)./(STATES(:,68)+CONSTANTS(:,555)));
RATES(:,28) = (ALGEBRAIC(:,54)+ CONSTANTS(:,753).*ALGEBRAIC(:,243)) - ( (CONSTANTS(:,556)./( CONSTANTS(:,758).*2.00000)).*ALGEBRAIC(:,235)+ CONSTANTS(:,754).*ALGEBRAIC(:,236));
ALGEBRAIC(:,238) = power(CONSTANTS(:,533), CONSTANTS(:,534))./(power(CONSTANTS(:,533), CONSTANTS(:,534))+power(STATES(:,68), CONSTANTS(:,534)));
ALGEBRAIC(:,240) = CONSTANTS(:,531).*(power(STATES(:,18), 2.00000)./(power(CONSTANTS(:,532), 2.00000)+power(STATES(:,18), 2.00000))).*ALGEBRAIC(:,238);
ALGEBRAIC(:,246) = 1.00000./(1.00000+( CONSTANTS(:,561).*CONSTANTS(:,562))./power(CONSTANTS(:,561)+STATES(:,68), 2.00000));
RATES(:,68) = ALGEBRAIC(:,246).*(ALGEBRAIC(:,240) - ALGEBRAIC(:,243));
ALGEBRAIC(:,242) = ( CONSTANTS(:,537).*CONSTANTS(:,542)+ CONSTANTS(:,539).*STATES(:,18))./(CONSTANTS(:,542)+STATES(:,18));
ALGEBRAIC(:,245) = ( (CONSTANTS(:,536)+CONSTANTS(:,540)).*CONSTANTS(:,543))./(CONSTANTS(:,543)+STATES(:,18));
ALGEBRAIC(:,248) = CONSTANTS(:,535).*power(( CONSTANTS(:,553).*ALGEBRAIC(:,242).*STATES(:,69))./( CONSTANTS(:,553).*ALGEBRAIC(:,242)+ALGEBRAIC(:,245)), 4.00000).*(STATES(:,67) - STATES(:,18));
RATES(:,18) = (ALGEBRAIC(:,43)+ CONSTANTS(:,756).*ALGEBRAIC(:,248)) - ( CONSTANTS(:,757).*ALGEBRAIC(:,54)+ CONSTANTS(:,755).*ALGEBRAIC(:,240));
RATES(:,67) = ALGEBRAIC(:,236) - ALGEBRAIC(:,248);
ALGEBRAIC(:,249) = ( CONSTANTS(:,538).*CONSTANTS(:,543)+ CONSTANTS(:,541).*STATES(:,18))./(CONSTANTS(:,543)+STATES(:,18));
ALGEBRAIC(:,251) = ( CONSTANTS(:,552).*STATES(:,70))./( (1.00000+power(CONSTANTS(:,548)./STATES(:,18), CONSTANTS(:,550))).*(1.00000+power(STATES(:,18)./CONSTANTS(:,549), CONSTANTS(:,551))));
RATES(:,69) = ALGEBRAIC(:,251).*(1.00000 - STATES(:,69)) - (( CONSTANTS(:,553).*ALGEBRAIC(:,242).*ALGEBRAIC(:,249))./( CONSTANTS(:,553).*ALGEBRAIC(:,242)+ALGEBRAIC(:,245))).*STATES(:,69);
ALGEBRAIC(:,252) = (( CONSTANTS(:,571).*CONSTANTS(:,570))./CONSTANTS(:,572)).*log(CONSTANTS(:,575)./STATES(:,71));
ALGEBRAIC(:,253) = CONSTANTS(:,574).*(STATES(:,2) - ALGEBRAIC(:,252));
ALGEBRAIC(:,261) = CONSTANTS(:,583).*(power(STATES(:,71), CONSTANTS(:,585))./(power(CONSTANTS(:,584), CONSTANTS(:,585))+power(STATES(:,71), CONSTANTS(:,585)))).*(CONSTANTS(:,586) - STATES(:,2));
ALGEBRAIC(:,256) = CONSTANTS(:,580).*(power(CONSTANTS(:,579), CONSTANTS(:,577))./(power(CONSTANTS(:,579), CONSTANTS(:,577))+power(STATES(:,29), CONSTANTS(:,577))));
ALGEBRAIC(:,257) = ALGEBRAIC(:,256).*(STATES(:,2) - CONSTANTS(:,578));
ALGEBRAIC(:,263) = CONSTANTS(:,764).*(ALGEBRAIC(:,257)+ALGEBRAIC(:,253)+ALGEBRAIC(:,261));
RATES(:,71) = - (CONSTANTS(:,617)./( 1.00000.*CONSTANTS(:,764))).*ALGEBRAIC(:,263);
ALGEBRAIC(:,55) = CONSTANTS(:,700).*(STATES(:,19) - STATES(:,29));
ALGEBRAIC(:,244) = ( CONSTANTS(:,566).*exp( CONSTANTS(:,567).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,568).*(STATES(:,2) - CONSTANTS(:,569))));
ALGEBRAIC(:,247) = (( CONSTANTS(:,571).*CONSTANTS(:,570))./( 2.00000.*CONSTANTS(:,572))).*log(CONSTANTS(:,573)./STATES(:,29));
ALGEBRAIC(:,250) = ALGEBRAIC(:,244).*(STATES(:,2) - ALGEBRAIC(:,247));
ALGEBRAIC(:,254) = CONSTANTS(:,576).*(power(CONSTANTS(:,579), CONSTANTS(:,577))./(power(CONSTANTS(:,579), CONSTANTS(:,577))+power(STATES(:,29), CONSTANTS(:,577))));
ALGEBRAIC(:,255) = ALGEBRAIC(:,254).*(STATES(:,2) - CONSTANTS(:,578));
ALGEBRAIC(:,258) = CONSTANTS(:,581).*(power(STATES(:,29), 2.00000)./(power(CONSTANTS(:,582), 2.00000)+power(STATES(:,29), 2.00000)));
ALGEBRAIC(:,259) = CONSTANTS(:,764).*(ALGEBRAIC(:,250)+ALGEBRAIC(:,255)+ALGEBRAIC(:,258));
ALGEBRAIC(:,260) = ( CONSTANTS(:,587).*(STATES(:,29) - CONSTANTS(:,588).*STATES(:,72)))./(1.00000+ CONSTANTS(:,589).*STATES(:,29)+ CONSTANTS(:,590).*STATES(:,72)+ CONSTANTS(:,591).*STATES(:,29).*STATES(:,72));
ALGEBRAIC(:,267) = CONSTANTS(:,615).*(STATES(:,73)./(STATES(:,73)+CONSTANTS(:,616)));
RATES(:,29) = (ALGEBRAIC(:,55)+ CONSTANTS(:,759).*ALGEBRAIC(:,267)) - ( (CONSTANTS(:,617)./( CONSTANTS(:,764).*2.00000)).*ALGEBRAIC(:,259)+ CONSTANTS(:,760).*ALGEBRAIC(:,260));
ALGEBRAIC(:,262) = power(CONSTANTS(:,594), CONSTANTS(:,595))./(power(CONSTANTS(:,594), CONSTANTS(:,595))+power(STATES(:,73), CONSTANTS(:,595)));
ALGEBRAIC(:,264) = CONSTANTS(:,592).*(power(STATES(:,19), 2.00000)./(power(CONSTANTS(:,593), 2.00000)+power(STATES(:,19), 2.00000))).*ALGEBRAIC(:,262);
ALGEBRAIC(:,270) = 1.00000./(1.00000+( CONSTANTS(:,622).*CONSTANTS(:,623))./power(CONSTANTS(:,622)+STATES(:,73), 2.00000));
RATES(:,73) = ALGEBRAIC(:,270).*(ALGEBRAIC(:,264) - ALGEBRAIC(:,267));
ALGEBRAIC(:,266) = ( CONSTANTS(:,598).*CONSTANTS(:,603)+ CONSTANTS(:,600).*STATES(:,19))./(CONSTANTS(:,603)+STATES(:,19));
ALGEBRAIC(:,269) = ( (CONSTANTS(:,597)+CONSTANTS(:,601)).*CONSTANTS(:,604))./(CONSTANTS(:,604)+STATES(:,19));
ALGEBRAIC(:,272) = CONSTANTS(:,596).*power(( CONSTANTS(:,614).*ALGEBRAIC(:,266).*STATES(:,74))./( CONSTANTS(:,614).*ALGEBRAIC(:,266)+ALGEBRAIC(:,269)), 4.00000).*(STATES(:,72) - STATES(:,19));
RATES(:,19) = (ALGEBRAIC(:,44)+ CONSTANTS(:,762).*ALGEBRAIC(:,272)) - ( CONSTANTS(:,763).*ALGEBRAIC(:,55)+ CONSTANTS(:,761).*ALGEBRAIC(:,264));
RATES(:,72) = ALGEBRAIC(:,260) - ALGEBRAIC(:,272);
ALGEBRAIC(:,273) = ( CONSTANTS(:,599).*CONSTANTS(:,604)+ CONSTANTS(:,602).*STATES(:,19))./(CONSTANTS(:,604)+STATES(:,19));
ALGEBRAIC(:,275) = ( CONSTANTS(:,613).*STATES(:,75))./( (1.00000+power(CONSTANTS(:,609)./STATES(:,19), CONSTANTS(:,611))).*(1.00000+power(STATES(:,19)./CONSTANTS(:,610), CONSTANTS(:,612))));
RATES(:,74) = ALGEBRAIC(:,275).*(1.00000 - STATES(:,74)) - (( CONSTANTS(:,614).*ALGEBRAIC(:,266).*ALGEBRAIC(:,273))./( CONSTANTS(:,614).*ALGEBRAIC(:,266)+ALGEBRAIC(:,269))).*STATES(:,74);
ALGEBRAIC(:,276) = (( CONSTANTS(:,632).*CONSTANTS(:,631))./CONSTANTS(:,633)).*log(CONSTANTS(:,636)./STATES(:,76));
ALGEBRAIC(:,277) = CONSTANTS(:,635).*(STATES(:,2) - ALGEBRAIC(:,276));
ALGEBRAIC(:,284) = CONSTANTS(:,644).*(power(STATES(:,76), CONSTANTS(:,646))./(power(CONSTANTS(:,645), CONSTANTS(:,646))+power(STATES(:,76), CONSTANTS(:,646)))).*(CONSTANTS(:,647) - STATES(:,2));
ALGEBRAIC(:,280) = CONSTANTS(:,641).*(power(CONSTANTS(:,640), CONSTANTS(:,638))./(power(CONSTANTS(:,640), CONSTANTS(:,638))+power(STATES(:,30), CONSTANTS(:,638))));
ALGEBRAIC(:,281) = ALGEBRAIC(:,280).*(STATES(:,2) - CONSTANTS(:,639));
ALGEBRAIC(:,286) = CONSTANTS(:,770).*(ALGEBRAIC(:,281)+ALGEBRAIC(:,277)+ALGEBRAIC(:,284));
RATES(:,76) = - (CONSTANTS(:,678)./( 1.00000.*CONSTANTS(:,770))).*ALGEBRAIC(:,286);
ALGEBRAIC(:,30) = CONSTANTS(:,16).*STATES(:,6).*STATES(:,7).*(STATES(:,2) - CONSTANTS(:,689));
ALGEBRAIC(:,31) = CONSTANTS(:,28).*STATES(:,8).*(STATES(:,2) - CONSTANTS(:,689));
ALGEBRAIC(:,32) = CONSTANTS(:,33).*(STATES(:,2) - CONSTANTS(:,689));
ALGEBRAIC(:,33) = CONSTANTS(:,34).*STATES(:,9).*STATES(:,10).*(STATES(:,2) - CONSTANTS(:,690));
ALGEBRAIC(:,34) = CONSTANTS(:,49).*(STATES(:,2) - CONSTANTS(:,50));
ALGEBRAIC(:,35) = ALGEBRAIC(:,8)+ALGEBRAIC(:,26)+ALGEBRAIC(:,30)+ALGEBRAIC(:,31)+ALGEBRAIC(:,32)+ALGEBRAIC(:,33)+ALGEBRAIC(:,34);
ALGEBRAIC(:,73) = ALGEBRAIC(:,67)+ALGEBRAIC(:,71);
ALGEBRAIC(:,97) = ALGEBRAIC(:,91)+ALGEBRAIC(:,95);
ALGEBRAIC(:,121) = ALGEBRAIC(:,115)+ALGEBRAIC(:,119);
ALGEBRAIC(:,145) = ALGEBRAIC(:,139)+ALGEBRAIC(:,143);
ALGEBRAIC(:,169) = ALGEBRAIC(:,163)+ALGEBRAIC(:,167);
ALGEBRAIC(:,193) = ALGEBRAIC(:,187)+ALGEBRAIC(:,191);
ALGEBRAIC(:,217) = ALGEBRAIC(:,211)+ALGEBRAIC(:,215);
ALGEBRAIC(:,241) = ALGEBRAIC(:,235)+ALGEBRAIC(:,239);
ALGEBRAIC(:,265) = ALGEBRAIC(:,259)+ALGEBRAIC(:,263);
ALGEBRAIC(:,268) = ( CONSTANTS(:,627).*exp( CONSTANTS(:,628).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,629).*(STATES(:,2) - CONSTANTS(:,630))));
ALGEBRAIC(:,271) = (( CONSTANTS(:,632).*CONSTANTS(:,631))./( 2.00000.*CONSTANTS(:,633))).*log(CONSTANTS(:,634)./STATES(:,30));
ALGEBRAIC(:,274) = ALGEBRAIC(:,268).*(STATES(:,2) - ALGEBRAIC(:,271));
ALGEBRAIC(:,278) = CONSTANTS(:,637).*(power(CONSTANTS(:,640), CONSTANTS(:,638))./(power(CONSTANTS(:,640), CONSTANTS(:,638))+power(STATES(:,30), CONSTANTS(:,638))));
ALGEBRAIC(:,279) = ALGEBRAIC(:,278).*(STATES(:,2) - CONSTANTS(:,639));
ALGEBRAIC(:,282) = CONSTANTS(:,642).*(power(STATES(:,30), 2.00000)./(power(CONSTANTS(:,643), 2.00000)+power(STATES(:,30), 2.00000)));
ALGEBRAIC(:,283) = CONSTANTS(:,770).*(ALGEBRAIC(:,274)+ALGEBRAIC(:,279)+ALGEBRAIC(:,282));
ALGEBRAIC(:,288) = ALGEBRAIC(:,283)+ALGEBRAIC(:,286);
RATES(:,2) = - (1.00000./CONSTANTS(:,51)).*(ALGEBRAIC(:,35)+ALGEBRAIC(:,73)+ALGEBRAIC(:,97)+ALGEBRAIC(:,121)+ALGEBRAIC(:,145)+ALGEBRAIC(:,169)+ALGEBRAIC(:,193)+ALGEBRAIC(:,217)+ALGEBRAIC(:,241)+ALGEBRAIC(:,265)+ALGEBRAIC(:,288));
ALGEBRAIC(:,56) = CONSTANTS(:,701).*(STATES(:,20) - STATES(:,30));
ALGEBRAIC(:,285) = ( CONSTANTS(:,648).*(STATES(:,30) - CONSTANTS(:,649).*STATES(:,77)))./(1.00000+ CONSTANTS(:,650).*STATES(:,30)+ CONSTANTS(:,651).*STATES(:,77)+ CONSTANTS(:,652).*STATES(:,30).*STATES(:,77));
ALGEBRAIC(:,291) = CONSTANTS(:,676).*(STATES(:,78)./(STATES(:,78)+CONSTANTS(:,677)));
RATES(:,30) = (ALGEBRAIC(:,56)+ CONSTANTS(:,765).*ALGEBRAIC(:,291)) - ( (CONSTANTS(:,678)./( CONSTANTS(:,770).*2.00000)).*ALGEBRAIC(:,283)+ CONSTANTS(:,766).*ALGEBRAIC(:,285));
ALGEBRAIC(:,287) = power(CONSTANTS(:,655), CONSTANTS(:,656))./(power(CONSTANTS(:,655), CONSTANTS(:,656))+power(STATES(:,78), CONSTANTS(:,656)));
ALGEBRAIC(:,289) = CONSTANTS(:,653).*(power(STATES(:,20), 2.00000)./(power(CONSTANTS(:,654), 2.00000)+power(STATES(:,20), 2.00000))).*ALGEBRAIC(:,287);
ALGEBRAIC(:,293) = 1.00000./(1.00000+( CONSTANTS(:,683).*CONSTANTS(:,684))./power(CONSTANTS(:,683)+STATES(:,78), 2.00000));
RATES(:,78) = ALGEBRAIC(:,293).*(ALGEBRAIC(:,289) - ALGEBRAIC(:,291));
ALGEBRAIC(:,290) = ( CONSTANTS(:,659).*CONSTANTS(:,664)+ CONSTANTS(:,661).*STATES(:,20))./(CONSTANTS(:,664)+STATES(:,20));
ALGEBRAIC(:,292) = ( (CONSTANTS(:,658)+CONSTANTS(:,662)).*CONSTANTS(:,665))./(CONSTANTS(:,665)+STATES(:,20));
ALGEBRAIC(:,294) = CONSTANTS(:,657).*power(( CONSTANTS(:,675).*ALGEBRAIC(:,290).*STATES(:,79))./( CONSTANTS(:,675).*ALGEBRAIC(:,290)+ALGEBRAIC(:,292)), 4.00000).*(STATES(:,77) - STATES(:,20));
RATES(:,20) = (ALGEBRAIC(:,45)+ CONSTANTS(:,768).*ALGEBRAIC(:,294)) - ( CONSTANTS(:,769).*ALGEBRAIC(:,56)+ CONSTANTS(:,767).*ALGEBRAIC(:,289));
RATES(:,77) = ALGEBRAIC(:,285) - ALGEBRAIC(:,294);
ALGEBRAIC(:,295) = ( CONSTANTS(:,660).*CONSTANTS(:,665)+ CONSTANTS(:,663).*STATES(:,20))./(CONSTANTS(:,665)+STATES(:,20));
ALGEBRAIC(:,296) = ( CONSTANTS(:,674).*STATES(:,80))./( (1.00000+power(CONSTANTS(:,670)./STATES(:,20), CONSTANTS(:,672))).*(1.00000+power(STATES(:,20)./CONSTANTS(:,671), CONSTANTS(:,673))));
RATES(:,79) = ALGEBRAIC(:,296).*(1.00000 - STATES(:,79)) - (( CONSTANTS(:,675).*ALGEBRAIC(:,290).*ALGEBRAIC(:,295))./( CONSTANTS(:,675).*ALGEBRAIC(:,290)+ALGEBRAIC(:,292))).*STATES(:,79);
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(:,1) = 1.00000./(1.00000+exp( CONSTANTS(:,6).*(STATES(:,2) - CONSTANTS(:,7))));
ALGEBRAIC(:,10) = CONSTANTS(:,117)./(1.00000+power(CONSTANTS(:,120)./STATES(:,11), CONSTANTS(:,118)));
ALGEBRAIC(:,11) = CONSTANTS(:,178)./(1.00000+power(CONSTANTS(:,181)./STATES(:,12), CONSTANTS(:,179)));
ALGEBRAIC(:,12) = CONSTANTS(:,239)./(1.00000+power(CONSTANTS(:,242)./STATES(:,13), CONSTANTS(:,240)));
ALGEBRAIC(:,13) = CONSTANTS(:,300)./(1.00000+power(CONSTANTS(:,303)./STATES(:,14), CONSTANTS(:,301)));
ALGEBRAIC(:,14) = CONSTANTS(:,361)./(1.00000+power(CONSTANTS(:,364)./STATES(:,15), CONSTANTS(:,362)));
ALGEBRAIC(:,15) = CONSTANTS(:,422)./(1.00000+power(CONSTANTS(:,425)./STATES(:,16), CONSTANTS(:,423)));
ALGEBRAIC(:,16) = CONSTANTS(:,483)./(1.00000+power(CONSTANTS(:,486)./STATES(:,17), CONSTANTS(:,484)));
ALGEBRAIC(:,17) = CONSTANTS(:,544)./(1.00000+power(CONSTANTS(:,547)./STATES(:,18), CONSTANTS(:,545)));
ALGEBRAIC(:,18) = CONSTANTS(:,605)./(1.00000+power(CONSTANTS(:,608)./STATES(:,19), CONSTANTS(:,606)));
ALGEBRAIC(:,19) = CONSTANTS(:,666)./(1.00000+power(CONSTANTS(:,669)./STATES(:,20), CONSTANTS(:,667)));
ALGEBRAIC(:,2) = 1.00000./(1.00000+exp( CONSTANTS(:,12).*(STATES(:,2) - CONSTANTS(:,13))));
ALGEBRAIC(:,20) = CONSTANTS(:,8)+(CONSTANTS(:,9) - CONSTANTS(:,8))./(1.00000+exp( CONSTANTS(:,10).*(STATES(:,2) - CONSTANTS(:,11))));
ALGEBRAIC(:,4) = CONSTANTS(:,24).*(CONSTANTS(:,25)./(1.00000+exp( CONSTANTS(:,26).*(STATES(:,2) - CONSTANTS(:,27))))+(1.00000 - CONSTANTS(:,25)));
ALGEBRAIC(:,22) = CONSTANTS(:,24).*(CONSTANTS(:,25) - CONSTANTS(:,25)./(1.00000+exp( CONSTANTS(:,26).*(STATES(:,2) - CONSTANTS(:,27)))));
ALGEBRAIC(:,5) = CONSTANTS(:,29).*(CONSTANTS(:,30)./(1.00000+exp( CONSTANTS(:,31).*(STATES(:,2) - CONSTANTS(:,32))))+(1.00000 - CONSTANTS(:,30)));
ALGEBRAIC(:,23) = CONSTANTS(:,29).*(CONSTANTS(:,30) - CONSTANTS(:,30)./(1.00000+exp( CONSTANTS(:,31).*(STATES(:,2) - CONSTANTS(:,32)))));
ALGEBRAIC(:,3) = CONSTANTS(:,17)+(CONSTANTS(:,18) - CONSTANTS(:,17))./(1.00000+exp( CONSTANTS(:,19).*(STATES(:,2) - CONSTANTS(:,20))));
ALGEBRAIC(:,21) = ALGEBRAIC(:,3).*(CONSTANTS(:,21)./(1.00000+exp( CONSTANTS(:,22).*(STATES(:,2) - CONSTANTS(:,23))))+(1.00000 - CONSTANTS(:,21)));
ALGEBRAIC(:,27) = ALGEBRAIC(:,3).*(CONSTANTS(:,21) - CONSTANTS(:,21)./(1.00000+exp( CONSTANTS(:,22).*(STATES(:,2) - CONSTANTS(:,23)))));
ALGEBRAIC(:,6) = CONSTANTS(:,35)+(CONSTANTS(:,36) - CONSTANTS(:,35))./(1.00000+exp( CONSTANTS(:,37).*(STATES(:,2) - CONSTANTS(:,38))));
ALGEBRAIC(:,24) = ALGEBRAIC(:,6).*(CONSTANTS(:,39)./(1.00000+exp( CONSTANTS(:,40).*(STATES(:,2) - CONSTANTS(:,41))))+(1.00000 - CONSTANTS(:,39)));
ALGEBRAIC(:,28) = ALGEBRAIC(:,6).*(CONSTANTS(:,39) - CONSTANTS(:,39)./(1.00000+exp( CONSTANTS(:,40).*(STATES(:,2) - CONSTANTS(:,41)))));
ALGEBRAIC(:,7) = CONSTANTS(:,42)+(CONSTANTS(:,43) - CONSTANTS(:,42))./(1.00000+exp( CONSTANTS(:,44).*(STATES(:,2) - CONSTANTS(:,45))));
ALGEBRAIC(:,25) = ALGEBRAIC(:,7).*(CONSTANTS(:,46)./(1.00000+exp( CONSTANTS(:,47).*(STATES(:,2) - CONSTANTS(:,48))))+(1.00000 - CONSTANTS(:,46)));
ALGEBRAIC(:,29) = ALGEBRAIC(:,7).*(CONSTANTS(:,46) - CONSTANTS(:,46)./(1.00000+exp( CONSTANTS(:,47).*(STATES(:,2) - CONSTANTS(:,48)))));
ALGEBRAIC(:,8) = CONSTANTS(:,1).*STATES(:,1).*(STATES(:,2) - CONSTANTS(:,2));
ALGEBRAIC(:,26) = CONSTANTS(:,14).*(STATES(:,5)./(CONSTANTS(:,15)+STATES(:,5)));
ALGEBRAIC(:,46) = CONSTANTS(:,75).*(CONSTANTS(:,76) - STATES(:,5));
ALGEBRAIC(:,36) = CONSTANTS(:,62).*(STATES(:,5) - STATES(:,11));
ALGEBRAIC(:,37) = CONSTANTS(:,63).*(STATES(:,5) - STATES(:,12));
ALGEBRAIC(:,38) = CONSTANTS(:,64).*(STATES(:,5) - STATES(:,13));
ALGEBRAIC(:,39) = CONSTANTS(:,65).*(STATES(:,5) - STATES(:,14));
ALGEBRAIC(:,40) = CONSTANTS(:,66).*(STATES(:,5) - STATES(:,15));
ALGEBRAIC(:,41) = CONSTANTS(:,67).*(STATES(:,5) - STATES(:,16));
ALGEBRAIC(:,42) = CONSTANTS(:,68).*(STATES(:,5) - STATES(:,17));
ALGEBRAIC(:,43) = CONSTANTS(:,69).*(STATES(:,5) - STATES(:,18));
ALGEBRAIC(:,44) = CONSTANTS(:,70).*(STATES(:,5) - STATES(:,19));
ALGEBRAIC(:,45) = CONSTANTS(:,71).*(STATES(:,5) - STATES(:,20));
ALGEBRAIC(:,60) = (( CONSTANTS(:,83).*CONSTANTS(:,82))./CONSTANTS(:,84)).*log(CONSTANTS(:,87)./STATES(:,31));
ALGEBRAIC(:,61) = CONSTANTS(:,86).*(STATES(:,2) - ALGEBRAIC(:,60));
ALGEBRAIC(:,69) = CONSTANTS(:,95).*(power(STATES(:,31), CONSTANTS(:,97))./(power(CONSTANTS(:,96), CONSTANTS(:,97))+power(STATES(:,31), CONSTANTS(:,97)))).*(CONSTANTS(:,98) - STATES(:,2));
ALGEBRAIC(:,64) = CONSTANTS(:,92).*(power(CONSTANTS(:,91), CONSTANTS(:,89))./(power(CONSTANTS(:,91), CONSTANTS(:,89))+power(STATES(:,21), CONSTANTS(:,89))));
ALGEBRAIC(:,65) = ALGEBRAIC(:,64).*(STATES(:,2) - CONSTANTS(:,90));
ALGEBRAIC(:,71) = CONSTANTS(:,716).*(ALGEBRAIC(:,65)+ALGEBRAIC(:,61)+ALGEBRAIC(:,69));
ALGEBRAIC(:,47) = CONSTANTS(:,692).*(STATES(:,11) - STATES(:,21));
ALGEBRAIC(:,57) = ( CONSTANTS(:,78).*exp( CONSTANTS(:,79).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,80).*(STATES(:,2) - CONSTANTS(:,81))));
ALGEBRAIC(:,58) = (( CONSTANTS(:,83).*CONSTANTS(:,82))./( 2.00000.*CONSTANTS(:,84))).*log(CONSTANTS(:,85)./STATES(:,21));
ALGEBRAIC(:,59) = ALGEBRAIC(:,57).*(STATES(:,2) - ALGEBRAIC(:,58));
ALGEBRAIC(:,62) = CONSTANTS(:,88).*(power(CONSTANTS(:,91), CONSTANTS(:,89))./(power(CONSTANTS(:,91), CONSTANTS(:,89))+power(STATES(:,21), CONSTANTS(:,89))));
ALGEBRAIC(:,63) = ALGEBRAIC(:,62).*(STATES(:,2) - CONSTANTS(:,90));
ALGEBRAIC(:,66) = CONSTANTS(:,93).*(power(STATES(:,21), 2.00000)./(power(CONSTANTS(:,94), 2.00000)+power(STATES(:,21), 2.00000)));
ALGEBRAIC(:,67) = CONSTANTS(:,716).*(ALGEBRAIC(:,59)+ALGEBRAIC(:,63)+ALGEBRAIC(:,66));
ALGEBRAIC(:,68) = ( CONSTANTS(:,99).*(STATES(:,21) - CONSTANTS(:,100).*STATES(:,32)))./(1.00000+ CONSTANTS(:,101).*STATES(:,21)+ CONSTANTS(:,102).*STATES(:,32)+ CONSTANTS(:,103).*STATES(:,21).*STATES(:,32));
ALGEBRAIC(:,75) = CONSTANTS(:,127).*(STATES(:,33)./(STATES(:,33)+CONSTANTS(:,128)));
ALGEBRAIC(:,70) = power(CONSTANTS(:,106), CONSTANTS(:,107))./(power(CONSTANTS(:,106), CONSTANTS(:,107))+power(STATES(:,33), CONSTANTS(:,107)));
ALGEBRAIC(:,72) = CONSTANTS(:,104).*(power(STATES(:,11), 2.00000)./(power(CONSTANTS(:,105), 2.00000)+power(STATES(:,11), 2.00000))).*ALGEBRAIC(:,70);
ALGEBRAIC(:,78) = 1.00000./(1.00000+( CONSTANTS(:,134).*CONSTANTS(:,135))./power(CONSTANTS(:,134)+STATES(:,33), 2.00000));
ALGEBRAIC(:,74) = ( CONSTANTS(:,110).*CONSTANTS(:,115)+ CONSTANTS(:,112).*STATES(:,11))./(CONSTANTS(:,115)+STATES(:,11));
ALGEBRAIC(:,77) = ( (CONSTANTS(:,109)+CONSTANTS(:,113)).*CONSTANTS(:,116))./(CONSTANTS(:,116)+STATES(:,11));
ALGEBRAIC(:,80) = CONSTANTS(:,108).*power(( CONSTANTS(:,126).*ALGEBRAIC(:,74).*STATES(:,34))./( CONSTANTS(:,126).*ALGEBRAIC(:,74)+ALGEBRAIC(:,77)), 4.00000).*(STATES(:,32) - STATES(:,11));
ALGEBRAIC(:,81) = ( CONSTANTS(:,111).*CONSTANTS(:,116)+ CONSTANTS(:,114).*STATES(:,11))./(CONSTANTS(:,116)+STATES(:,11));
ALGEBRAIC(:,83) = ( CONSTANTS(:,125).*STATES(:,35))./( (1.00000+power(CONSTANTS(:,121)./STATES(:,11), CONSTANTS(:,123))).*(1.00000+power(STATES(:,11)./CONSTANTS(:,122), CONSTANTS(:,124))));
ALGEBRAIC(:,84) = (( CONSTANTS(:,144).*CONSTANTS(:,143))./CONSTANTS(:,145)).*log(CONSTANTS(:,148)./STATES(:,36));
ALGEBRAIC(:,85) = CONSTANTS(:,147).*(STATES(:,2) - ALGEBRAIC(:,84));
ALGEBRAIC(:,93) = CONSTANTS(:,156).*(power(STATES(:,36), CONSTANTS(:,158))./(power(CONSTANTS(:,157), CONSTANTS(:,158))+power(STATES(:,36), CONSTANTS(:,158)))).*(CONSTANTS(:,159) - STATES(:,2));
ALGEBRAIC(:,88) = CONSTANTS(:,153).*(power(CONSTANTS(:,152), CONSTANTS(:,150))./(power(CONSTANTS(:,152), CONSTANTS(:,150))+power(STATES(:,22), CONSTANTS(:,150))));
ALGEBRAIC(:,89) = ALGEBRAIC(:,88).*(STATES(:,2) - CONSTANTS(:,151));
ALGEBRAIC(:,95) = CONSTANTS(:,722).*(ALGEBRAIC(:,89)+ALGEBRAIC(:,85)+ALGEBRAIC(:,93));
ALGEBRAIC(:,48) = CONSTANTS(:,693).*(STATES(:,12) - STATES(:,22));
ALGEBRAIC(:,76) = ( CONSTANTS(:,139).*exp( CONSTANTS(:,140).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,141).*(STATES(:,2) - CONSTANTS(:,142))));
ALGEBRAIC(:,79) = (( CONSTANTS(:,144).*CONSTANTS(:,143))./( 2.00000.*CONSTANTS(:,145))).*log(CONSTANTS(:,146)./STATES(:,22));
ALGEBRAIC(:,82) = ALGEBRAIC(:,76).*(STATES(:,2) - ALGEBRAIC(:,79));
ALGEBRAIC(:,86) = CONSTANTS(:,149).*(power(CONSTANTS(:,152), CONSTANTS(:,150))./(power(CONSTANTS(:,152), CONSTANTS(:,150))+power(STATES(:,22), CONSTANTS(:,150))));
ALGEBRAIC(:,87) = ALGEBRAIC(:,86).*(STATES(:,2) - CONSTANTS(:,151));
ALGEBRAIC(:,90) = CONSTANTS(:,154).*(power(STATES(:,22), 2.00000)./(power(CONSTANTS(:,155), 2.00000)+power(STATES(:,22), 2.00000)));
ALGEBRAIC(:,91) = CONSTANTS(:,722).*(ALGEBRAIC(:,82)+ALGEBRAIC(:,87)+ALGEBRAIC(:,90));
ALGEBRAIC(:,92) = ( CONSTANTS(:,160).*(STATES(:,22) - CONSTANTS(:,161).*STATES(:,37)))./(1.00000+ CONSTANTS(:,162).*STATES(:,22)+ CONSTANTS(:,163).*STATES(:,37)+ CONSTANTS(:,164).*STATES(:,22).*STATES(:,37));
ALGEBRAIC(:,99) = CONSTANTS(:,188).*(STATES(:,38)./(STATES(:,38)+CONSTANTS(:,189)));
ALGEBRAIC(:,94) = power(CONSTANTS(:,167), CONSTANTS(:,168))./(power(CONSTANTS(:,167), CONSTANTS(:,168))+power(STATES(:,38), CONSTANTS(:,168)));
ALGEBRAIC(:,96) = CONSTANTS(:,165).*(power(STATES(:,12), 2.00000)./(power(CONSTANTS(:,166), 2.00000)+power(STATES(:,12), 2.00000))).*ALGEBRAIC(:,94);
ALGEBRAIC(:,102) = 1.00000./(1.00000+( CONSTANTS(:,195).*CONSTANTS(:,196))./power(CONSTANTS(:,195)+STATES(:,38), 2.00000));
ALGEBRAIC(:,98) = ( CONSTANTS(:,171).*CONSTANTS(:,176)+ CONSTANTS(:,173).*STATES(:,12))./(CONSTANTS(:,176)+STATES(:,12));
ALGEBRAIC(:,101) = ( (CONSTANTS(:,170)+CONSTANTS(:,174)).*CONSTANTS(:,177))./(CONSTANTS(:,177)+STATES(:,12));
ALGEBRAIC(:,104) = CONSTANTS(:,169).*power(( CONSTANTS(:,187).*ALGEBRAIC(:,98).*STATES(:,39))./( CONSTANTS(:,187).*ALGEBRAIC(:,98)+ALGEBRAIC(:,101)), 4.00000).*(STATES(:,37) - STATES(:,12));
ALGEBRAIC(:,105) = ( CONSTANTS(:,172).*CONSTANTS(:,177)+ CONSTANTS(:,175).*STATES(:,12))./(CONSTANTS(:,177)+STATES(:,12));
ALGEBRAIC(:,107) = ( CONSTANTS(:,186).*STATES(:,40))./( (1.00000+power(CONSTANTS(:,182)./STATES(:,12), CONSTANTS(:,184))).*(1.00000+power(STATES(:,12)./CONSTANTS(:,183), CONSTANTS(:,185))));
ALGEBRAIC(:,108) = (( CONSTANTS(:,205).*CONSTANTS(:,204))./CONSTANTS(:,206)).*log(CONSTANTS(:,209)./STATES(:,41));
ALGEBRAIC(:,109) = CONSTANTS(:,208).*(STATES(:,2) - ALGEBRAIC(:,108));
ALGEBRAIC(:,117) = CONSTANTS(:,217).*(power(STATES(:,41), CONSTANTS(:,219))./(power(CONSTANTS(:,218), CONSTANTS(:,219))+power(STATES(:,41), CONSTANTS(:,219)))).*(CONSTANTS(:,220) - STATES(:,2));
ALGEBRAIC(:,112) = CONSTANTS(:,214).*(power(CONSTANTS(:,213), CONSTANTS(:,211))./(power(CONSTANTS(:,213), CONSTANTS(:,211))+power(STATES(:,23), CONSTANTS(:,211))));
ALGEBRAIC(:,113) = ALGEBRAIC(:,112).*(STATES(:,2) - CONSTANTS(:,212));
ALGEBRAIC(:,119) = CONSTANTS(:,728).*(ALGEBRAIC(:,113)+ALGEBRAIC(:,109)+ALGEBRAIC(:,117));
ALGEBRAIC(:,49) = CONSTANTS(:,694).*(STATES(:,13) - STATES(:,23));
ALGEBRAIC(:,100) = ( CONSTANTS(:,200).*exp( CONSTANTS(:,201).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,202).*(STATES(:,2) - CONSTANTS(:,203))));
ALGEBRAIC(:,103) = (( CONSTANTS(:,205).*CONSTANTS(:,204))./( 2.00000.*CONSTANTS(:,206))).*log(CONSTANTS(:,207)./STATES(:,23));
ALGEBRAIC(:,106) = ALGEBRAIC(:,100).*(STATES(:,2) - ALGEBRAIC(:,103));
ALGEBRAIC(:,110) = CONSTANTS(:,210).*(power(CONSTANTS(:,213), CONSTANTS(:,211))./(power(CONSTANTS(:,213), CONSTANTS(:,211))+power(STATES(:,23), CONSTANTS(:,211))));
ALGEBRAIC(:,111) = ALGEBRAIC(:,110).*(STATES(:,2) - CONSTANTS(:,212));
ALGEBRAIC(:,114) = CONSTANTS(:,215).*(power(STATES(:,23), 2.00000)./(power(CONSTANTS(:,216), 2.00000)+power(STATES(:,23), 2.00000)));
ALGEBRAIC(:,115) = CONSTANTS(:,728).*(ALGEBRAIC(:,106)+ALGEBRAIC(:,111)+ALGEBRAIC(:,114));
ALGEBRAIC(:,116) = ( CONSTANTS(:,221).*(STATES(:,23) - CONSTANTS(:,222).*STATES(:,42)))./(1.00000+ CONSTANTS(:,223).*STATES(:,23)+ CONSTANTS(:,224).*STATES(:,42)+ CONSTANTS(:,225).*STATES(:,23).*STATES(:,42));
ALGEBRAIC(:,123) = CONSTANTS(:,249).*(STATES(:,43)./(STATES(:,43)+CONSTANTS(:,250)));
ALGEBRAIC(:,118) = power(CONSTANTS(:,228), CONSTANTS(:,229))./(power(CONSTANTS(:,228), CONSTANTS(:,229))+power(STATES(:,43), CONSTANTS(:,229)));
ALGEBRAIC(:,120) = CONSTANTS(:,226).*(power(STATES(:,13), 2.00000)./(power(CONSTANTS(:,227), 2.00000)+power(STATES(:,13), 2.00000))).*ALGEBRAIC(:,118);
ALGEBRAIC(:,126) = 1.00000./(1.00000+( CONSTANTS(:,256).*CONSTANTS(:,257))./power(CONSTANTS(:,256)+STATES(:,43), 2.00000));
ALGEBRAIC(:,122) = ( CONSTANTS(:,232).*CONSTANTS(:,237)+ CONSTANTS(:,234).*STATES(:,13))./(CONSTANTS(:,237)+STATES(:,13));
ALGEBRAIC(:,125) = ( (CONSTANTS(:,231)+CONSTANTS(:,235)).*CONSTANTS(:,238))./(CONSTANTS(:,238)+STATES(:,13));
ALGEBRAIC(:,128) = CONSTANTS(:,230).*power(( CONSTANTS(:,248).*ALGEBRAIC(:,122).*STATES(:,44))./( CONSTANTS(:,248).*ALGEBRAIC(:,122)+ALGEBRAIC(:,125)), 4.00000).*(STATES(:,42) - STATES(:,13));
ALGEBRAIC(:,129) = ( CONSTANTS(:,233).*CONSTANTS(:,238)+ CONSTANTS(:,236).*STATES(:,13))./(CONSTANTS(:,238)+STATES(:,13));
ALGEBRAIC(:,131) = ( CONSTANTS(:,247).*STATES(:,45))./( (1.00000+power(CONSTANTS(:,243)./STATES(:,13), CONSTANTS(:,245))).*(1.00000+power(STATES(:,13)./CONSTANTS(:,244), CONSTANTS(:,246))));
ALGEBRAIC(:,132) = (( CONSTANTS(:,266).*CONSTANTS(:,265))./CONSTANTS(:,267)).*log(CONSTANTS(:,270)./STATES(:,46));
ALGEBRAIC(:,133) = CONSTANTS(:,269).*(STATES(:,2) - ALGEBRAIC(:,132));
ALGEBRAIC(:,141) = CONSTANTS(:,278).*(power(STATES(:,46), CONSTANTS(:,280))./(power(CONSTANTS(:,279), CONSTANTS(:,280))+power(STATES(:,46), CONSTANTS(:,280)))).*(CONSTANTS(:,281) - STATES(:,2));
ALGEBRAIC(:,136) = CONSTANTS(:,275).*(power(CONSTANTS(:,274), CONSTANTS(:,272))./(power(CONSTANTS(:,274), CONSTANTS(:,272))+power(STATES(:,24), CONSTANTS(:,272))));
ALGEBRAIC(:,137) = ALGEBRAIC(:,136).*(STATES(:,2) - CONSTANTS(:,273));
ALGEBRAIC(:,143) = CONSTANTS(:,734).*(ALGEBRAIC(:,137)+ALGEBRAIC(:,133)+ALGEBRAIC(:,141));
ALGEBRAIC(:,50) = CONSTANTS(:,695).*(STATES(:,14) - STATES(:,24));
ALGEBRAIC(:,124) = ( CONSTANTS(:,261).*exp( CONSTANTS(:,262).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,263).*(STATES(:,2) - CONSTANTS(:,264))));
ALGEBRAIC(:,127) = (( CONSTANTS(:,266).*CONSTANTS(:,265))./( 2.00000.*CONSTANTS(:,267))).*log(CONSTANTS(:,268)./STATES(:,24));
ALGEBRAIC(:,130) = ALGEBRAIC(:,124).*(STATES(:,2) - ALGEBRAIC(:,127));
ALGEBRAIC(:,134) = CONSTANTS(:,271).*(power(CONSTANTS(:,274), CONSTANTS(:,272))./(power(CONSTANTS(:,274), CONSTANTS(:,272))+power(STATES(:,24), CONSTANTS(:,272))));
ALGEBRAIC(:,135) = ALGEBRAIC(:,134).*(STATES(:,2) - CONSTANTS(:,273));
ALGEBRAIC(:,138) = CONSTANTS(:,276).*(power(STATES(:,24), 2.00000)./(power(CONSTANTS(:,277), 2.00000)+power(STATES(:,24), 2.00000)));
ALGEBRAIC(:,139) = CONSTANTS(:,734).*(ALGEBRAIC(:,130)+ALGEBRAIC(:,135)+ALGEBRAIC(:,138));
ALGEBRAIC(:,140) = ( CONSTANTS(:,282).*(STATES(:,24) - CONSTANTS(:,283).*STATES(:,47)))./(1.00000+ CONSTANTS(:,284).*STATES(:,24)+ CONSTANTS(:,285).*STATES(:,47)+ CONSTANTS(:,286).*STATES(:,24).*STATES(:,47));
ALGEBRAIC(:,147) = CONSTANTS(:,310).*(STATES(:,48)./(STATES(:,48)+CONSTANTS(:,311)));
ALGEBRAIC(:,142) = power(CONSTANTS(:,289), CONSTANTS(:,290))./(power(CONSTANTS(:,289), CONSTANTS(:,290))+power(STATES(:,48), CONSTANTS(:,290)));
ALGEBRAIC(:,144) = CONSTANTS(:,287).*(power(STATES(:,14), 2.00000)./(power(CONSTANTS(:,288), 2.00000)+power(STATES(:,14), 2.00000))).*ALGEBRAIC(:,142);
ALGEBRAIC(:,150) = 1.00000./(1.00000+( CONSTANTS(:,317).*CONSTANTS(:,318))./power(CONSTANTS(:,317)+STATES(:,48), 2.00000));
ALGEBRAIC(:,146) = ( CONSTANTS(:,293).*CONSTANTS(:,298)+ CONSTANTS(:,295).*STATES(:,14))./(CONSTANTS(:,298)+STATES(:,14));
ALGEBRAIC(:,149) = ( (CONSTANTS(:,292)+CONSTANTS(:,296)).*CONSTANTS(:,299))./(CONSTANTS(:,299)+STATES(:,14));
ALGEBRAIC(:,152) = CONSTANTS(:,291).*power(( CONSTANTS(:,309).*ALGEBRAIC(:,146).*STATES(:,49))./( CONSTANTS(:,309).*ALGEBRAIC(:,146)+ALGEBRAIC(:,149)), 4.00000).*(STATES(:,47) - STATES(:,14));
ALGEBRAIC(:,153) = ( CONSTANTS(:,294).*CONSTANTS(:,299)+ CONSTANTS(:,297).*STATES(:,14))./(CONSTANTS(:,299)+STATES(:,14));
ALGEBRAIC(:,155) = ( CONSTANTS(:,308).*STATES(:,50))./( (1.00000+power(CONSTANTS(:,304)./STATES(:,14), CONSTANTS(:,306))).*(1.00000+power(STATES(:,14)./CONSTANTS(:,305), CONSTANTS(:,307))));
ALGEBRAIC(:,156) = (( CONSTANTS(:,327).*CONSTANTS(:,326))./CONSTANTS(:,328)).*log(CONSTANTS(:,331)./STATES(:,51));
ALGEBRAIC(:,157) = CONSTANTS(:,330).*(STATES(:,2) - ALGEBRAIC(:,156));
ALGEBRAIC(:,165) = CONSTANTS(:,339).*(power(STATES(:,51), CONSTANTS(:,341))./(power(CONSTANTS(:,340), CONSTANTS(:,341))+power(STATES(:,51), CONSTANTS(:,341)))).*(CONSTANTS(:,342) - STATES(:,2));
ALGEBRAIC(:,160) = CONSTANTS(:,336).*(power(CONSTANTS(:,335), CONSTANTS(:,333))./(power(CONSTANTS(:,335), CONSTANTS(:,333))+power(STATES(:,25), CONSTANTS(:,333))));
ALGEBRAIC(:,161) = ALGEBRAIC(:,160).*(STATES(:,2) - CONSTANTS(:,334));
ALGEBRAIC(:,167) = CONSTANTS(:,740).*(ALGEBRAIC(:,161)+ALGEBRAIC(:,157)+ALGEBRAIC(:,165));
ALGEBRAIC(:,51) = CONSTANTS(:,696).*(STATES(:,15) - STATES(:,25));
ALGEBRAIC(:,148) = ( CONSTANTS(:,322).*exp( CONSTANTS(:,323).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,324).*(STATES(:,2) - CONSTANTS(:,325))));
ALGEBRAIC(:,151) = (( CONSTANTS(:,327).*CONSTANTS(:,326))./( 2.00000.*CONSTANTS(:,328))).*log(CONSTANTS(:,329)./STATES(:,25));
ALGEBRAIC(:,154) = ALGEBRAIC(:,148).*(STATES(:,2) - ALGEBRAIC(:,151));
ALGEBRAIC(:,158) = CONSTANTS(:,332).*(power(CONSTANTS(:,335), CONSTANTS(:,333))./(power(CONSTANTS(:,335), CONSTANTS(:,333))+power(STATES(:,25), CONSTANTS(:,333))));
ALGEBRAIC(:,159) = ALGEBRAIC(:,158).*(STATES(:,2) - CONSTANTS(:,334));
ALGEBRAIC(:,162) = CONSTANTS(:,337).*(power(STATES(:,25), 2.00000)./(power(CONSTANTS(:,338), 2.00000)+power(STATES(:,25), 2.00000)));
ALGEBRAIC(:,163) = CONSTANTS(:,740).*(ALGEBRAIC(:,154)+ALGEBRAIC(:,159)+ALGEBRAIC(:,162));
ALGEBRAIC(:,164) = ( CONSTANTS(:,343).*(STATES(:,25) - CONSTANTS(:,344).*STATES(:,52)))./(1.00000+ CONSTANTS(:,345).*STATES(:,25)+ CONSTANTS(:,346).*STATES(:,52)+ CONSTANTS(:,347).*STATES(:,25).*STATES(:,52));
ALGEBRAIC(:,171) = CONSTANTS(:,371).*(STATES(:,53)./(STATES(:,53)+CONSTANTS(:,372)));
ALGEBRAIC(:,166) = power(CONSTANTS(:,350), CONSTANTS(:,351))./(power(CONSTANTS(:,350), CONSTANTS(:,351))+power(STATES(:,53), CONSTANTS(:,351)));
ALGEBRAIC(:,168) = CONSTANTS(:,348).*(power(STATES(:,15), 2.00000)./(power(CONSTANTS(:,349), 2.00000)+power(STATES(:,15), 2.00000))).*ALGEBRAIC(:,166);
ALGEBRAIC(:,174) = 1.00000./(1.00000+( CONSTANTS(:,378).*CONSTANTS(:,379))./power(CONSTANTS(:,378)+STATES(:,53), 2.00000));
ALGEBRAIC(:,170) = ( CONSTANTS(:,354).*CONSTANTS(:,359)+ CONSTANTS(:,356).*STATES(:,15))./(CONSTANTS(:,359)+STATES(:,15));
ALGEBRAIC(:,173) = ( (CONSTANTS(:,353)+CONSTANTS(:,357)).*CONSTANTS(:,360))./(CONSTANTS(:,360)+STATES(:,15));
ALGEBRAIC(:,176) = CONSTANTS(:,352).*power(( CONSTANTS(:,370).*ALGEBRAIC(:,170).*STATES(:,54))./( CONSTANTS(:,370).*ALGEBRAIC(:,170)+ALGEBRAIC(:,173)), 4.00000).*(STATES(:,52) - STATES(:,15));
ALGEBRAIC(:,177) = ( CONSTANTS(:,355).*CONSTANTS(:,360)+ CONSTANTS(:,358).*STATES(:,15))./(CONSTANTS(:,360)+STATES(:,15));
ALGEBRAIC(:,179) = ( CONSTANTS(:,369).*STATES(:,55))./( (1.00000+power(CONSTANTS(:,365)./STATES(:,15), CONSTANTS(:,367))).*(1.00000+power(STATES(:,15)./CONSTANTS(:,366), CONSTANTS(:,368))));
ALGEBRAIC(:,180) = (( CONSTANTS(:,388).*CONSTANTS(:,387))./CONSTANTS(:,389)).*log(CONSTANTS(:,392)./STATES(:,56));
ALGEBRAIC(:,181) = CONSTANTS(:,391).*(STATES(:,2) - ALGEBRAIC(:,180));
ALGEBRAIC(:,189) = CONSTANTS(:,400).*(power(STATES(:,56), CONSTANTS(:,402))./(power(CONSTANTS(:,401), CONSTANTS(:,402))+power(STATES(:,56), CONSTANTS(:,402)))).*(CONSTANTS(:,403) - STATES(:,2));
ALGEBRAIC(:,184) = CONSTANTS(:,397).*(power(CONSTANTS(:,396), CONSTANTS(:,394))./(power(CONSTANTS(:,396), CONSTANTS(:,394))+power(STATES(:,26), CONSTANTS(:,394))));
ALGEBRAIC(:,185) = ALGEBRAIC(:,184).*(STATES(:,2) - CONSTANTS(:,395));
ALGEBRAIC(:,191) = CONSTANTS(:,746).*(ALGEBRAIC(:,185)+ALGEBRAIC(:,181)+ALGEBRAIC(:,189));
ALGEBRAIC(:,52) = CONSTANTS(:,697).*(STATES(:,16) - STATES(:,26));
ALGEBRAIC(:,172) = ( CONSTANTS(:,383).*exp( CONSTANTS(:,384).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,385).*(STATES(:,2) - CONSTANTS(:,386))));
ALGEBRAIC(:,175) = (( CONSTANTS(:,388).*CONSTANTS(:,387))./( 2.00000.*CONSTANTS(:,389))).*log(CONSTANTS(:,390)./STATES(:,26));
ALGEBRAIC(:,178) = ALGEBRAIC(:,172).*(STATES(:,2) - ALGEBRAIC(:,175));
ALGEBRAIC(:,182) = CONSTANTS(:,393).*(power(CONSTANTS(:,396), CONSTANTS(:,394))./(power(CONSTANTS(:,396), CONSTANTS(:,394))+power(STATES(:,26), CONSTANTS(:,394))));
ALGEBRAIC(:,183) = ALGEBRAIC(:,182).*(STATES(:,2) - CONSTANTS(:,395));
ALGEBRAIC(:,186) = CONSTANTS(:,398).*(power(STATES(:,26), 2.00000)./(power(CONSTANTS(:,399), 2.00000)+power(STATES(:,26), 2.00000)));
ALGEBRAIC(:,187) = CONSTANTS(:,746).*(ALGEBRAIC(:,178)+ALGEBRAIC(:,183)+ALGEBRAIC(:,186));
ALGEBRAIC(:,188) = ( CONSTANTS(:,404).*(STATES(:,26) - CONSTANTS(:,405).*STATES(:,57)))./(1.00000+ CONSTANTS(:,406).*STATES(:,26)+ CONSTANTS(:,407).*STATES(:,57)+ CONSTANTS(:,408).*STATES(:,26).*STATES(:,57));
ALGEBRAIC(:,195) = CONSTANTS(:,432).*(STATES(:,58)./(STATES(:,58)+CONSTANTS(:,433)));
ALGEBRAIC(:,190) = power(CONSTANTS(:,411), CONSTANTS(:,412))./(power(CONSTANTS(:,411), CONSTANTS(:,412))+power(STATES(:,58), CONSTANTS(:,412)));
ALGEBRAIC(:,192) = CONSTANTS(:,409).*(power(STATES(:,16), 2.00000)./(power(CONSTANTS(:,410), 2.00000)+power(STATES(:,16), 2.00000))).*ALGEBRAIC(:,190);
ALGEBRAIC(:,198) = 1.00000./(1.00000+( CONSTANTS(:,439).*CONSTANTS(:,440))./power(CONSTANTS(:,439)+STATES(:,58), 2.00000));
ALGEBRAIC(:,194) = ( CONSTANTS(:,415).*CONSTANTS(:,420)+ CONSTANTS(:,417).*STATES(:,16))./(CONSTANTS(:,420)+STATES(:,16));
ALGEBRAIC(:,197) = ( (CONSTANTS(:,414)+CONSTANTS(:,418)).*CONSTANTS(:,421))./(CONSTANTS(:,421)+STATES(:,16));
ALGEBRAIC(:,200) = CONSTANTS(:,413).*power(( CONSTANTS(:,431).*ALGEBRAIC(:,194).*STATES(:,59))./( CONSTANTS(:,431).*ALGEBRAIC(:,194)+ALGEBRAIC(:,197)), 4.00000).*(STATES(:,57) - STATES(:,16));
ALGEBRAIC(:,201) = ( CONSTANTS(:,416).*CONSTANTS(:,421)+ CONSTANTS(:,419).*STATES(:,16))./(CONSTANTS(:,421)+STATES(:,16));
ALGEBRAIC(:,203) = ( CONSTANTS(:,430).*STATES(:,60))./( (1.00000+power(CONSTANTS(:,426)./STATES(:,16), CONSTANTS(:,428))).*(1.00000+power(STATES(:,16)./CONSTANTS(:,427), CONSTANTS(:,429))));
ALGEBRAIC(:,204) = (( CONSTANTS(:,449).*CONSTANTS(:,448))./CONSTANTS(:,450)).*log(CONSTANTS(:,453)./STATES(:,61));
ALGEBRAIC(:,205) = CONSTANTS(:,452).*(STATES(:,2) - ALGEBRAIC(:,204));
ALGEBRAIC(:,213) = CONSTANTS(:,461).*(power(STATES(:,61), CONSTANTS(:,463))./(power(CONSTANTS(:,462), CONSTANTS(:,463))+power(STATES(:,61), CONSTANTS(:,463)))).*(CONSTANTS(:,464) - STATES(:,2));
ALGEBRAIC(:,208) = CONSTANTS(:,458).*(power(CONSTANTS(:,457), CONSTANTS(:,455))./(power(CONSTANTS(:,457), CONSTANTS(:,455))+power(STATES(:,27), CONSTANTS(:,455))));
ALGEBRAIC(:,209) = ALGEBRAIC(:,208).*(STATES(:,2) - CONSTANTS(:,456));
ALGEBRAIC(:,215) = CONSTANTS(:,752).*(ALGEBRAIC(:,209)+ALGEBRAIC(:,205)+ALGEBRAIC(:,213));
ALGEBRAIC(:,53) = CONSTANTS(:,698).*(STATES(:,17) - STATES(:,27));
ALGEBRAIC(:,196) = ( CONSTANTS(:,444).*exp( CONSTANTS(:,445).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,446).*(STATES(:,2) - CONSTANTS(:,447))));
ALGEBRAIC(:,199) = (( CONSTANTS(:,449).*CONSTANTS(:,448))./( 2.00000.*CONSTANTS(:,450))).*log(CONSTANTS(:,451)./STATES(:,27));
ALGEBRAIC(:,202) = ALGEBRAIC(:,196).*(STATES(:,2) - ALGEBRAIC(:,199));
ALGEBRAIC(:,206) = CONSTANTS(:,454).*(power(CONSTANTS(:,457), CONSTANTS(:,455))./(power(CONSTANTS(:,457), CONSTANTS(:,455))+power(STATES(:,27), CONSTANTS(:,455))));
ALGEBRAIC(:,207) = ALGEBRAIC(:,206).*(STATES(:,2) - CONSTANTS(:,456));
ALGEBRAIC(:,210) = CONSTANTS(:,459).*(power(STATES(:,27), 2.00000)./(power(CONSTANTS(:,460), 2.00000)+power(STATES(:,27), 2.00000)));
ALGEBRAIC(:,211) = CONSTANTS(:,752).*(ALGEBRAIC(:,202)+ALGEBRAIC(:,207)+ALGEBRAIC(:,210));
ALGEBRAIC(:,212) = ( CONSTANTS(:,465).*(STATES(:,27) - CONSTANTS(:,466).*STATES(:,62)))./(1.00000+ CONSTANTS(:,467).*STATES(:,27)+ CONSTANTS(:,468).*STATES(:,62)+ CONSTANTS(:,469).*STATES(:,27).*STATES(:,62));
ALGEBRAIC(:,219) = CONSTANTS(:,493).*(STATES(:,63)./(STATES(:,63)+CONSTANTS(:,494)));
ALGEBRAIC(:,214) = power(CONSTANTS(:,472), CONSTANTS(:,473))./(power(CONSTANTS(:,472), CONSTANTS(:,473))+power(STATES(:,63), CONSTANTS(:,473)));
ALGEBRAIC(:,216) = CONSTANTS(:,470).*(power(STATES(:,17), 2.00000)./(power(CONSTANTS(:,471), 2.00000)+power(STATES(:,17), 2.00000))).*ALGEBRAIC(:,214);
ALGEBRAIC(:,222) = 1.00000./(1.00000+( CONSTANTS(:,500).*CONSTANTS(:,501))./power(CONSTANTS(:,500)+STATES(:,63), 2.00000));
ALGEBRAIC(:,218) = ( CONSTANTS(:,476).*CONSTANTS(:,481)+ CONSTANTS(:,478).*STATES(:,17))./(CONSTANTS(:,481)+STATES(:,17));
ALGEBRAIC(:,221) = ( (CONSTANTS(:,475)+CONSTANTS(:,479)).*CONSTANTS(:,482))./(CONSTANTS(:,482)+STATES(:,17));
ALGEBRAIC(:,224) = CONSTANTS(:,474).*power(( CONSTANTS(:,492).*ALGEBRAIC(:,218).*STATES(:,64))./( CONSTANTS(:,492).*ALGEBRAIC(:,218)+ALGEBRAIC(:,221)), 4.00000).*(STATES(:,62) - STATES(:,17));
ALGEBRAIC(:,225) = ( CONSTANTS(:,477).*CONSTANTS(:,482)+ CONSTANTS(:,480).*STATES(:,17))./(CONSTANTS(:,482)+STATES(:,17));
ALGEBRAIC(:,227) = ( CONSTANTS(:,491).*STATES(:,65))./( (1.00000+power(CONSTANTS(:,487)./STATES(:,17), CONSTANTS(:,489))).*(1.00000+power(STATES(:,17)./CONSTANTS(:,488), CONSTANTS(:,490))));
ALGEBRAIC(:,228) = (( CONSTANTS(:,510).*CONSTANTS(:,509))./CONSTANTS(:,511)).*log(CONSTANTS(:,514)./STATES(:,66));
ALGEBRAIC(:,229) = CONSTANTS(:,513).*(STATES(:,2) - ALGEBRAIC(:,228));
ALGEBRAIC(:,237) = CONSTANTS(:,522).*(power(STATES(:,66), CONSTANTS(:,524))./(power(CONSTANTS(:,523), CONSTANTS(:,524))+power(STATES(:,66), CONSTANTS(:,524)))).*(CONSTANTS(:,525) - STATES(:,2));
ALGEBRAIC(:,232) = CONSTANTS(:,519).*(power(CONSTANTS(:,518), CONSTANTS(:,516))./(power(CONSTANTS(:,518), CONSTANTS(:,516))+power(STATES(:,28), CONSTANTS(:,516))));
ALGEBRAIC(:,233) = ALGEBRAIC(:,232).*(STATES(:,2) - CONSTANTS(:,517));
ALGEBRAIC(:,239) = CONSTANTS(:,758).*(ALGEBRAIC(:,233)+ALGEBRAIC(:,229)+ALGEBRAIC(:,237));
ALGEBRAIC(:,54) = CONSTANTS(:,699).*(STATES(:,18) - STATES(:,28));
ALGEBRAIC(:,220) = ( CONSTANTS(:,505).*exp( CONSTANTS(:,506).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,507).*(STATES(:,2) - CONSTANTS(:,508))));
ALGEBRAIC(:,223) = (( CONSTANTS(:,510).*CONSTANTS(:,509))./( 2.00000.*CONSTANTS(:,511))).*log(CONSTANTS(:,512)./STATES(:,28));
ALGEBRAIC(:,226) = ALGEBRAIC(:,220).*(STATES(:,2) - ALGEBRAIC(:,223));
ALGEBRAIC(:,230) = CONSTANTS(:,515).*(power(CONSTANTS(:,518), CONSTANTS(:,516))./(power(CONSTANTS(:,518), CONSTANTS(:,516))+power(STATES(:,28), CONSTANTS(:,516))));
ALGEBRAIC(:,231) = ALGEBRAIC(:,230).*(STATES(:,2) - CONSTANTS(:,517));
ALGEBRAIC(:,234) = CONSTANTS(:,520).*(power(STATES(:,28), 2.00000)./(power(CONSTANTS(:,521), 2.00000)+power(STATES(:,28), 2.00000)));
ALGEBRAIC(:,235) = CONSTANTS(:,758).*(ALGEBRAIC(:,226)+ALGEBRAIC(:,231)+ALGEBRAIC(:,234));
ALGEBRAIC(:,236) = ( CONSTANTS(:,526).*(STATES(:,28) - CONSTANTS(:,527).*STATES(:,67)))./(1.00000+ CONSTANTS(:,528).*STATES(:,28)+ CONSTANTS(:,529).*STATES(:,67)+ CONSTANTS(:,530).*STATES(:,28).*STATES(:,67));
ALGEBRAIC(:,243) = CONSTANTS(:,554).*(STATES(:,68)./(STATES(:,68)+CONSTANTS(:,555)));
ALGEBRAIC(:,238) = power(CONSTANTS(:,533), CONSTANTS(:,534))./(power(CONSTANTS(:,533), CONSTANTS(:,534))+power(STATES(:,68), CONSTANTS(:,534)));
ALGEBRAIC(:,240) = CONSTANTS(:,531).*(power(STATES(:,18), 2.00000)./(power(CONSTANTS(:,532), 2.00000)+power(STATES(:,18), 2.00000))).*ALGEBRAIC(:,238);
ALGEBRAIC(:,246) = 1.00000./(1.00000+( CONSTANTS(:,561).*CONSTANTS(:,562))./power(CONSTANTS(:,561)+STATES(:,68), 2.00000));
ALGEBRAIC(:,242) = ( CONSTANTS(:,537).*CONSTANTS(:,542)+ CONSTANTS(:,539).*STATES(:,18))./(CONSTANTS(:,542)+STATES(:,18));
ALGEBRAIC(:,245) = ( (CONSTANTS(:,536)+CONSTANTS(:,540)).*CONSTANTS(:,543))./(CONSTANTS(:,543)+STATES(:,18));
ALGEBRAIC(:,248) = CONSTANTS(:,535).*power(( CONSTANTS(:,553).*ALGEBRAIC(:,242).*STATES(:,69))./( CONSTANTS(:,553).*ALGEBRAIC(:,242)+ALGEBRAIC(:,245)), 4.00000).*(STATES(:,67) - STATES(:,18));
ALGEBRAIC(:,249) = ( CONSTANTS(:,538).*CONSTANTS(:,543)+ CONSTANTS(:,541).*STATES(:,18))./(CONSTANTS(:,543)+STATES(:,18));
ALGEBRAIC(:,251) = ( CONSTANTS(:,552).*STATES(:,70))./( (1.00000+power(CONSTANTS(:,548)./STATES(:,18), CONSTANTS(:,550))).*(1.00000+power(STATES(:,18)./CONSTANTS(:,549), CONSTANTS(:,551))));
ALGEBRAIC(:,252) = (( CONSTANTS(:,571).*CONSTANTS(:,570))./CONSTANTS(:,572)).*log(CONSTANTS(:,575)./STATES(:,71));
ALGEBRAIC(:,253) = CONSTANTS(:,574).*(STATES(:,2) - ALGEBRAIC(:,252));
ALGEBRAIC(:,261) = CONSTANTS(:,583).*(power(STATES(:,71), CONSTANTS(:,585))./(power(CONSTANTS(:,584), CONSTANTS(:,585))+power(STATES(:,71), CONSTANTS(:,585)))).*(CONSTANTS(:,586) - STATES(:,2));
ALGEBRAIC(:,256) = CONSTANTS(:,580).*(power(CONSTANTS(:,579), CONSTANTS(:,577))./(power(CONSTANTS(:,579), CONSTANTS(:,577))+power(STATES(:,29), CONSTANTS(:,577))));
ALGEBRAIC(:,257) = ALGEBRAIC(:,256).*(STATES(:,2) - CONSTANTS(:,578));
ALGEBRAIC(:,263) = CONSTANTS(:,764).*(ALGEBRAIC(:,257)+ALGEBRAIC(:,253)+ALGEBRAIC(:,261));
ALGEBRAIC(:,55) = CONSTANTS(:,700).*(STATES(:,19) - STATES(:,29));
ALGEBRAIC(:,244) = ( CONSTANTS(:,566).*exp( CONSTANTS(:,567).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,568).*(STATES(:,2) - CONSTANTS(:,569))));
ALGEBRAIC(:,247) = (( CONSTANTS(:,571).*CONSTANTS(:,570))./( 2.00000.*CONSTANTS(:,572))).*log(CONSTANTS(:,573)./STATES(:,29));
ALGEBRAIC(:,250) = ALGEBRAIC(:,244).*(STATES(:,2) - ALGEBRAIC(:,247));
ALGEBRAIC(:,254) = CONSTANTS(:,576).*(power(CONSTANTS(:,579), CONSTANTS(:,577))./(power(CONSTANTS(:,579), CONSTANTS(:,577))+power(STATES(:,29), CONSTANTS(:,577))));
ALGEBRAIC(:,255) = ALGEBRAIC(:,254).*(STATES(:,2) - CONSTANTS(:,578));
ALGEBRAIC(:,258) = CONSTANTS(:,581).*(power(STATES(:,29), 2.00000)./(power(CONSTANTS(:,582), 2.00000)+power(STATES(:,29), 2.00000)));
ALGEBRAIC(:,259) = CONSTANTS(:,764).*(ALGEBRAIC(:,250)+ALGEBRAIC(:,255)+ALGEBRAIC(:,258));
ALGEBRAIC(:,260) = ( CONSTANTS(:,587).*(STATES(:,29) - CONSTANTS(:,588).*STATES(:,72)))./(1.00000+ CONSTANTS(:,589).*STATES(:,29)+ CONSTANTS(:,590).*STATES(:,72)+ CONSTANTS(:,591).*STATES(:,29).*STATES(:,72));
ALGEBRAIC(:,267) = CONSTANTS(:,615).*(STATES(:,73)./(STATES(:,73)+CONSTANTS(:,616)));
ALGEBRAIC(:,262) = power(CONSTANTS(:,594), CONSTANTS(:,595))./(power(CONSTANTS(:,594), CONSTANTS(:,595))+power(STATES(:,73), CONSTANTS(:,595)));
ALGEBRAIC(:,264) = CONSTANTS(:,592).*(power(STATES(:,19), 2.00000)./(power(CONSTANTS(:,593), 2.00000)+power(STATES(:,19), 2.00000))).*ALGEBRAIC(:,262);
ALGEBRAIC(:,270) = 1.00000./(1.00000+( CONSTANTS(:,622).*CONSTANTS(:,623))./power(CONSTANTS(:,622)+STATES(:,73), 2.00000));
ALGEBRAIC(:,266) = ( CONSTANTS(:,598).*CONSTANTS(:,603)+ CONSTANTS(:,600).*STATES(:,19))./(CONSTANTS(:,603)+STATES(:,19));
ALGEBRAIC(:,269) = ( (CONSTANTS(:,597)+CONSTANTS(:,601)).*CONSTANTS(:,604))./(CONSTANTS(:,604)+STATES(:,19));
ALGEBRAIC(:,272) = CONSTANTS(:,596).*power(( CONSTANTS(:,614).*ALGEBRAIC(:,266).*STATES(:,74))./( CONSTANTS(:,614).*ALGEBRAIC(:,266)+ALGEBRAIC(:,269)), 4.00000).*(STATES(:,72) - STATES(:,19));
ALGEBRAIC(:,273) = ( CONSTANTS(:,599).*CONSTANTS(:,604)+ CONSTANTS(:,602).*STATES(:,19))./(CONSTANTS(:,604)+STATES(:,19));
ALGEBRAIC(:,275) = ( CONSTANTS(:,613).*STATES(:,75))./( (1.00000+power(CONSTANTS(:,609)./STATES(:,19), CONSTANTS(:,611))).*(1.00000+power(STATES(:,19)./CONSTANTS(:,610), CONSTANTS(:,612))));
ALGEBRAIC(:,276) = (( CONSTANTS(:,632).*CONSTANTS(:,631))./CONSTANTS(:,633)).*log(CONSTANTS(:,636)./STATES(:,76));
ALGEBRAIC(:,277) = CONSTANTS(:,635).*(STATES(:,2) - ALGEBRAIC(:,276));
ALGEBRAIC(:,284) = CONSTANTS(:,644).*(power(STATES(:,76), CONSTANTS(:,646))./(power(CONSTANTS(:,645), CONSTANTS(:,646))+power(STATES(:,76), CONSTANTS(:,646)))).*(CONSTANTS(:,647) - STATES(:,2));
ALGEBRAIC(:,280) = CONSTANTS(:,641).*(power(CONSTANTS(:,640), CONSTANTS(:,638))./(power(CONSTANTS(:,640), CONSTANTS(:,638))+power(STATES(:,30), CONSTANTS(:,638))));
ALGEBRAIC(:,281) = ALGEBRAIC(:,280).*(STATES(:,2) - CONSTANTS(:,639));
ALGEBRAIC(:,286) = CONSTANTS(:,770).*(ALGEBRAIC(:,281)+ALGEBRAIC(:,277)+ALGEBRAIC(:,284));
ALGEBRAIC(:,30) = CONSTANTS(:,16).*STATES(:,6).*STATES(:,7).*(STATES(:,2) - CONSTANTS(:,689));
ALGEBRAIC(:,31) = CONSTANTS(:,28).*STATES(:,8).*(STATES(:,2) - CONSTANTS(:,689));
ALGEBRAIC(:,32) = CONSTANTS(:,33).*(STATES(:,2) - CONSTANTS(:,689));
ALGEBRAIC(:,33) = CONSTANTS(:,34).*STATES(:,9).*STATES(:,10).*(STATES(:,2) - CONSTANTS(:,690));
ALGEBRAIC(:,34) = CONSTANTS(:,49).*(STATES(:,2) - CONSTANTS(:,50));
ALGEBRAIC(:,35) = ALGEBRAIC(:,8)+ALGEBRAIC(:,26)+ALGEBRAIC(:,30)+ALGEBRAIC(:,31)+ALGEBRAIC(:,32)+ALGEBRAIC(:,33)+ALGEBRAIC(:,34);
ALGEBRAIC(:,73) = ALGEBRAIC(:,67)+ALGEBRAIC(:,71);
ALGEBRAIC(:,97) = ALGEBRAIC(:,91)+ALGEBRAIC(:,95);
ALGEBRAIC(:,121) = ALGEBRAIC(:,115)+ALGEBRAIC(:,119);
ALGEBRAIC(:,145) = ALGEBRAIC(:,139)+ALGEBRAIC(:,143);
ALGEBRAIC(:,169) = ALGEBRAIC(:,163)+ALGEBRAIC(:,167);
ALGEBRAIC(:,193) = ALGEBRAIC(:,187)+ALGEBRAIC(:,191);
ALGEBRAIC(:,217) = ALGEBRAIC(:,211)+ALGEBRAIC(:,215);
ALGEBRAIC(:,241) = ALGEBRAIC(:,235)+ALGEBRAIC(:,239);
ALGEBRAIC(:,265) = ALGEBRAIC(:,259)+ALGEBRAIC(:,263);
ALGEBRAIC(:,268) = ( CONSTANTS(:,627).*exp( CONSTANTS(:,628).*STATES(:,2)))./(1.00000+exp( CONSTANTS(:,629).*(STATES(:,2) - CONSTANTS(:,630))));
ALGEBRAIC(:,271) = (( CONSTANTS(:,632).*CONSTANTS(:,631))./( 2.00000.*CONSTANTS(:,633))).*log(CONSTANTS(:,634)./STATES(:,30));
ALGEBRAIC(:,274) = ALGEBRAIC(:,268).*(STATES(:,2) - ALGEBRAIC(:,271));
ALGEBRAIC(:,278) = CONSTANTS(:,637).*(power(CONSTANTS(:,640), CONSTANTS(:,638))./(power(CONSTANTS(:,640), CONSTANTS(:,638))+power(STATES(:,30), CONSTANTS(:,638))));
ALGEBRAIC(:,279) = ALGEBRAIC(:,278).*(STATES(:,2) - CONSTANTS(:,639));
ALGEBRAIC(:,282) = CONSTANTS(:,642).*(power(STATES(:,30), 2.00000)./(power(CONSTANTS(:,643), 2.00000)+power(STATES(:,30), 2.00000)));
ALGEBRAIC(:,283) = CONSTANTS(:,770).*(ALGEBRAIC(:,274)+ALGEBRAIC(:,279)+ALGEBRAIC(:,282));
ALGEBRAIC(:,288) = ALGEBRAIC(:,283)+ALGEBRAIC(:,286);
ALGEBRAIC(:,56) = CONSTANTS(:,701).*(STATES(:,20) - STATES(:,30));
ALGEBRAIC(:,285) = ( CONSTANTS(:,648).*(STATES(:,30) - CONSTANTS(:,649).*STATES(:,77)))./(1.00000+ CONSTANTS(:,650).*STATES(:,30)+ CONSTANTS(:,651).*STATES(:,77)+ CONSTANTS(:,652).*STATES(:,30).*STATES(:,77));
ALGEBRAIC(:,291) = CONSTANTS(:,676).*(STATES(:,78)./(STATES(:,78)+CONSTANTS(:,677)));
ALGEBRAIC(:,287) = power(CONSTANTS(:,655), CONSTANTS(:,656))./(power(CONSTANTS(:,655), CONSTANTS(:,656))+power(STATES(:,78), CONSTANTS(:,656)));
ALGEBRAIC(:,289) = CONSTANTS(:,653).*(power(STATES(:,20), 2.00000)./(power(CONSTANTS(:,654), 2.00000)+power(STATES(:,20), 2.00000))).*ALGEBRAIC(:,287);
ALGEBRAIC(:,293) = 1.00000./(1.00000+( CONSTANTS(:,683).*CONSTANTS(:,684))./power(CONSTANTS(:,683)+STATES(:,78), 2.00000));
ALGEBRAIC(:,290) = ( CONSTANTS(:,659).*CONSTANTS(:,664)+ CONSTANTS(:,661).*STATES(:,20))./(CONSTANTS(:,664)+STATES(:,20));
ALGEBRAIC(:,292) = ( (CONSTANTS(:,658)+CONSTANTS(:,662)).*CONSTANTS(:,665))./(CONSTANTS(:,665)+STATES(:,20));
ALGEBRAIC(:,294) = CONSTANTS(:,657).*power(( CONSTANTS(:,675).*ALGEBRAIC(:,290).*STATES(:,79))./( CONSTANTS(:,675).*ALGEBRAIC(:,290)+ALGEBRAIC(:,292)), 4.00000).*(STATES(:,77) - STATES(:,20));
ALGEBRAIC(:,295) = ( CONSTANTS(:,660).*CONSTANTS(:,665)+ CONSTANTS(:,663).*STATES(:,20))./(CONSTANTS(:,665)+STATES(:,20));
ALGEBRAIC(:,296) = ( CONSTANTS(:,674).*STATES(:,80))./( (1.00000+power(CONSTANTS(:,670)./STATES(:,20), CONSTANTS(:,672))).*(1.00000+power(STATES(:,20)./CONSTANTS(:,671), CONSTANTS(:,673))));
ALGEBRAIC(:,9) = (( CONSTANTS(:,53).*CONSTANTS(:,52))./( 2.00000.*CONSTANTS(:,54))).*log(CONSTANTS(:,55)./STATES(:,5));
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