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