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