Generated Code
The following is matlab code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
function [VOI, STATES, ALGEBRAIC, CONSTANTS] = mainFunction()
% This is the "main function". In Matlab, things work best if you rename this function to match the filename.
[VOI, STATES, ALGEBRAIC, CONSTANTS] = solveModel();
end
function [algebraicVariableCount] = getAlgebraicVariableCount()
% Used later when setting a global variable with the number of algebraic variables.
% Note: This is not the "main method".
algebraicVariableCount =0;
end
% There are a total of 77 entries in each of the rate and state variable arrays.
% There are a total of 200 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('t in component environment (s)');
LEGEND_CONSTANTS(:,1) = strpad('Ca in component Calcium (um)');
LEGEND_CONSTANTS(:,2) = strpad('GTP in component undefinedvariables (um)');
LEGEND_CONSTANTS(:,3) = strpad('PKA in component undefinedvariables (um)');
LEGEND_STATES(:,1) = strpad('PKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,2) = strpad('GAP in component B_Model_Modules (um)');
LEGEND_STATES(:,3) = strpad('Raf in component B_Model_Modules (um)');
LEGEND_STATES(:,4) = strpad('GEF_star in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,5) = strpad('CaM_GEF in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,6) = strpad('Gbg_GEF in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,7) = strpad('CaGqPLC in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,8) = strpad('CaPLC in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,9) = strpad('NgCaM in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,10) = strpad('GEF in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,11) = strpad('Ng in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,12) = strpad('Gbg in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,13) = strpad('Ga_GDP in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,14) = strpad('mGluR in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,15) = strpad('PLC in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,16) = strpad('Ca2CaM in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,17) = strpad('Ca3CaM in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,18) = strpad('Ng_star in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,19) = strpad('Gabg in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,20) = strpad('R in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,21) = strpad('Gabg_R in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,22) = strpad('Glu_synapse in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,23) = strpad('Gabg_GluR in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,24) = strpad('Glu in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,25) = strpad('Ga_GTP in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,26) = strpad('GqPLC in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,27) = strpad('GEF_inact in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,28) = strpad('Ca4CAM in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,29) = strpad('CaM in component DGBM_Model_Modules (um)');
LEGEND_CONSTANTS(:,4) = strpad('one in component DGBM_Model_Modules (dimensionless)');
LEGEND_CONSTANTS(:,5) = strpad('kfD1 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,6) = strpad('kbD1 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,7) = strpad('kfD2 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,8) = strpad('kbD2 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,9) = strpad('kfD3 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,10) = strpad('kbD3 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,11) = strpad('kfD4 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,12) = strpad('kbD4 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,13) = strpad('kfD5 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,14) = strpad('kbD5 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,15) = strpad('kfD6 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,16) = strpad('kfD7 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,17) = strpad('kfD8 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,18) = strpad('kfD9 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,19) = strpad('kfG1 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,20) = strpad('kbG1 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,21) = strpad('kfG2 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,22) = strpad('kbG2 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,23) = strpad('kfG3 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,24) = strpad('kbG3 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,25) = strpad('kfG4 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,26) = strpad('kbG4 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,27) = strpad('kfG5 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,28) = strpad('kfB2 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,29) = strpad('kbB2 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,30) = strpad('kmB3 in component DGBM_Model_Modules (um)');
LEGEND_CONSTANTS(:,31) = strpad('kfB5 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,32) = strpad('kbB5 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,33) = strpad('kfB7 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,34) = strpad('kbB7 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,35) = strpad('kfB8 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,36) = strpad('kbB8 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,37) = strpad('kmB1 in component DGBM_Model_Modules (um)');
LEGEND_CONSTANTS(:,38) = strpad('VmaxB1 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,39) = strpad('kmB6 in component DGBM_Model_Modules (um)');
LEGEND_CONSTANTS(:,40) = strpad('VmaxB6 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,41) = strpad('kfM1 in component DGBM_Model_Modules (per_um2_s)');
LEGEND_CONSTANTS(:,42) = strpad('kbM1 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,43) = strpad('kfM2 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,44) = strpad('kbM2 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,45) = strpad('kfM3 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,46) = strpad('kbM3 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,47) = strpad('kfM4 in component DGBM_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,48) = strpad('kbM4 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,49) = strpad('kfM6 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,50) = strpad('kmM5 in component DGBM_Model_Modules (um)');
LEGEND_CONSTANTS(:,51) = strpad('VmaxM5 in component DGBM_Model_Modules (um_per_s)');
LEGEND_CONSTANTS(:,52) = strpad('kmM7 in component DGBM_Model_Modules (um)');
LEGEND_CONSTANTS(:,53) = strpad('VmaxM7 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,54) = strpad('kmM8 in component DGBM_Model_Modules (um)');
LEGEND_CONSTANTS(:,55) = strpad('VmaxM8 in component DGBM_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,56) = strpad('kmH1 in component DGBM_Model_Modules (um)');
LEGEND_STATES(:,30) = strpad('SHCstar_SOS_GRB2 in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,31) = strpad('MAPK_star in component H_Model_Module_MAPK (um)');
LEGEND_STATES(:,32) = strpad('SOS in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,33) = strpad('MAPKK_star_star in component H_Model_Module_MAPKK (um)');
LEGEND_STATES(:,34) = strpad('MAPKK_star in component H_Model_Module_MAPKK (um)');
LEGEND_STATES(:,35) = strpad('Raf_star in component B_Model_Modules (um)');
LEGEND_STATES(:,36) = strpad('GTPRasRaf_star in component B_Model_Modules (um)');
LEGEND_STATES(:,37) = strpad('Raf_star_star in component B_Model_Modules (um)');
LEGEND_STATES(:,38) = strpad('GTP_RAS in component B_Model_Modules (um)');
LEGEND_STATES(:,39) = strpad('GDP_RAS in component B_Model_Modules (um)');
LEGEND_STATES(:,40) = strpad('GAPstar in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,57) = strpad('PP2A in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,58) = strpad('one in component B_Model_Modules (dimensionless)');
LEGEND_CONSTANTS(:,59) = strpad('kfB12 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,60) = strpad('kbB12 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,61) = strpad('kmB9 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,62) = strpad('VmaxB9 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,63) = strpad('kmB10 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,64) = strpad('VmaxB10 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,65) = strpad('kmB11 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,66) = strpad('VmaxB11 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,67) = strpad('kmB13 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,68) = strpad('VmaxB13 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,69) = strpad('kfB4 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,70) = strpad('kbB4 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,71) = strpad('kmB3 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,72) = strpad('VmaxB3 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,73) = strpad('kmA7 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,74) = strpad('VmaxA7 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,75) = strpad('kfH5 in component B_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,76) = strpad('kbH5 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,77) = strpad('kmH1 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,78) = strpad('kmH2 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,79) = strpad('kmH3 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,80) = strpad('kmH4 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,81) = strpad('VmaxH1 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,82) = strpad('VmaxH2 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,83) = strpad('VmaxH3 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,84) = strpad('VmaxH4 in component B_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,85) = strpad('kmM8 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,86) = strpad('kmB6 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,87) = strpad('kmM7 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,88) = strpad('kmA9 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,89) = strpad('kmH8 in component B_Model_Modules (um)');
LEGEND_CONSTANTS(:,90) = strpad('kmH9 in component B_Model_Modules (um)');
LEGEND_STATES(:,41) = strpad('MAPKK in component H_Model_Module_MAPKK (um)');
LEGEND_CONSTANTS(:,91) = strpad('PP2A in component H_Model_Module_MAPKK (um)');
LEGEND_CONSTANTS(:,92) = strpad('one in component H_Model_Module_MAPKK (dimensionless)');
LEGEND_CONSTANTS(:,93) = strpad('kmH6 in component H_Model_Module_MAPKK (um)');
LEGEND_CONSTANTS(:,94) = strpad('kmH7 in component H_Model_Module_MAPKK (um)');
LEGEND_CONSTANTS(:,95) = strpad('kmH8 in component H_Model_Module_MAPKK (um)');
LEGEND_CONSTANTS(:,96) = strpad('kmH9 in component H_Model_Module_MAPKK (um)');
LEGEND_CONSTANTS(:,97) = strpad('VmaxH6 in component H_Model_Module_MAPKK (per_s)');
LEGEND_CONSTANTS(:,98) = strpad('VmaxH7 in component H_Model_Module_MAPKK (per_s)');
LEGEND_CONSTANTS(:,99) = strpad('VmaxH8 in component H_Model_Module_MAPKK (per_s)');
LEGEND_CONSTANTS(:,100) = strpad('VmaxH9 in component H_Model_Module_MAPKK (per_s)');
LEGEND_CONSTANTS(:,101) = strpad('kmH3 in component H_Model_Module_MAPKK (um)');
LEGEND_CONSTANTS(:,102) = strpad('kmH4 in component H_Model_Module_MAPKK (um)');
LEGEND_STATES(:,42) = strpad('MAPK in component H_Model_Module_MAPK (um)');
LEGEND_STATES(:,43) = strpad('MAPK_tyr in component H_Model_Module_MAPK (um)');
LEGEND_CONSTANTS(:,103) = strpad('MKP1 in component H_Model_Module_MAPK (um)');
LEGEND_CONSTANTS(:,104) = strpad('one in component H_Model_Module_MAPK (dimensionless)');
LEGEND_CONSTANTS(:,105) = strpad('kmH10 in component H_Model_Module_MAPK (um)');
LEGEND_CONSTANTS(:,106) = strpad('kmH11 in component H_Model_Module_MAPK (um)');
LEGEND_CONSTANTS(:,107) = strpad('kmH12 in component H_Model_Module_MAPK (um)');
LEGEND_CONSTANTS(:,108) = strpad('kmH13 in component H_Model_Module_MAPK (um)');
LEGEND_CONSTANTS(:,109) = strpad('VmaxH10 in component H_Model_Module_MAPK (per_s)');
LEGEND_CONSTANTS(:,110) = strpad('VmaxH11 in component H_Model_Module_MAPK (per_s)');
LEGEND_CONSTANTS(:,111) = strpad('VmaxH12 in component H_Model_Module_MAPK (per_s)');
LEGEND_CONSTANTS(:,112) = strpad('VmaxH13 in component H_Model_Module_MAPK (per_s)');
LEGEND_STATES(:,44) = strpad('EGF_EGFR in component A_Model_Module_EGFR (um)');
LEGEND_STATES(:,45) = strpad('CaPLCg in component F_Model_Module_PLC (um)');
LEGEND_STATES(:,46) = strpad('SHC in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,47) = strpad('SHCstar in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,48) = strpad('SOSstar_GRB2 in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,49) = strpad('SOS_GRB2 in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,50) = strpad('GRB2 in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,51) = strpad('SOSstar in component A_Model_Module_SHC (um)');
LEGEND_CONSTANTS(:,113) = strpad('one in component A_Model_Module_SHC (dimensionless)');
LEGEND_CONSTANTS(:,114) = strpad('kmA3 in component A_Model_Module_SHC (um)');
LEGEND_CONSTANTS(:,115) = strpad('kfA4 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,116) = strpad('kbA4 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,117) = strpad('kfA5 in component A_Model_Module_SHC (per_um_s)');
LEGEND_CONSTANTS(:,118) = strpad('kbA5 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,119) = strpad('kfA6 in component A_Model_Module_SHC (per_um_s)');
LEGEND_CONSTANTS(:,120) = strpad('kbA6 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,121) = strpad('kfA8 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,122) = strpad('kbA8 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,123) = strpad('kfA10 in component A_Model_Module_SHC (per_um_s)');
LEGEND_CONSTANTS(:,124) = strpad('kbA10 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,125) = strpad('kmA7 in component A_Model_Module_SHC (um)');
LEGEND_CONSTANTS(:,126) = strpad('kmA9 in component A_Model_Module_SHC (um)');
LEGEND_CONSTANTS(:,127) = strpad('VmaxA7 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,128) = strpad('VmaxA3 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,129) = strpad('VmaxA9 in component A_Model_Module_SHC (per_s)');
LEGEND_CONSTANTS(:,130) = strpad('kmH2 in component A_Model_Module_SHC (um)');
LEGEND_CONSTANTS(:,131) = strpad('kmF4 in component A_Model_Module_SHC (um)');
LEGEND_STATES(:,52) = strpad('EGF_EGFR_INTERNAL in component A_Model_Module_EGFR (um)');
LEGEND_STATES(:,53) = strpad('EGF in component A_Model_Module_EGFR (um)');
LEGEND_STATES(:,54) = strpad('EGFR in component A_Model_Module_EGFR (um)');
LEGEND_CONSTANTS(:,132) = strpad('kfA1 in component A_Model_Module_EGFR (per_um_s)');
LEGEND_CONSTANTS(:,133) = strpad('kbA1 in component A_Model_Module_EGFR (per_s)');
LEGEND_CONSTANTS(:,134) = strpad('kfA2 in component A_Model_Module_EGFR (per_s)');
LEGEND_CONSTANTS(:,135) = strpad('kbA2 in component A_Model_Module_EGFR (per_s)');
LEGEND_STATES(:,55) = strpad('CaPLCg_star in component F_Model_Module_PLC (um)');
LEGEND_STATES(:,56) = strpad('PLCg in component F_Model_Module_PLC (um)');
LEGEND_STATES(:,57) = strpad('PLCg_star in component F_Model_Module_PLC (um)');
LEGEND_CONSTANTS(:,136) = strpad('kfF1 in component F_Model_Module_PLC (per_um_s)');
LEGEND_CONSTANTS(:,137) = strpad('kbF1 in component F_Model_Module_PLC (per_s)');
LEGEND_CONSTANTS(:,138) = strpad('kfF3 in component F_Model_Module_PLC (per_s)');
LEGEND_CONSTANTS(:,139) = strpad('kfF5 in component F_Model_Module_PLC (per_um_s)');
LEGEND_CONSTANTS(:,140) = strpad('kbF5 in component F_Model_Module_PLC (per_s)');
LEGEND_CONSTANTS(:,141) = strpad('kmF4 in component F_Model_Module_PLC (um)');
LEGEND_CONSTANTS(:,142) = strpad('VmaxF4 in component F_Model_Module_PLC (per_s)');
LEGEND_CONSTANTS(:,143) = strpad('kmA3 in component F_Model_Module_PLC (um)');
LEGEND_CONSTANTS(:,144) = strpad('one in component F_Model_Module_PLC (dimensionless)');
LEGEND_STATES(:,58) = strpad('PKC_i in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,59) = strpad('DAGPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,60) = strpad('AADAGPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,61) = strpad('DAGCaPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,62) = strpad('CaPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,63) = strpad('AA in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,64) = strpad('APC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,65) = strpad('PIP2_star in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,66) = strpad('PLA2_cyt in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,67) = strpad('PLA2_star in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,68) = strpad('CaPLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,69) = strpad('PIP2CaPLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,70) = strpad('PIP2PLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,71) = strpad('DAGCaPLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,72) = strpad('CaPLA2_star in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,73) = strpad('PIP2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,74) = strpad('DAG in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,75) = strpad('IP3 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,76) = strpad('Inositol in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_STATES(:,77) = strpad('PC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,145) = strpad('kfK1 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,146) = strpad('kbK1 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,147) = strpad('kfK2 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,148) = strpad('kbK2 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,149) = strpad('kfK3 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,150) = strpad('kbK3 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,151) = strpad('kfK4 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,152) = strpad('kbK4 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,153) = strpad('kfK5 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,154) = strpad('kbK5 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,155) = strpad('kfK6 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,156) = strpad('kbK6 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,157) = strpad('kfK7 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,158) = strpad('kbK7 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,159) = strpad('kfK8 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,160) = strpad('kbK8 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,161) = strpad('kfK9 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,162) = strpad('kbK9 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,163) = strpad('kfK10 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,164) = strpad('kbK10 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,165) = strpad('kfE1 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,166) = strpad('kbE1 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,167) = strpad('kfE3 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,168) = strpad('kbE3 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,169) = strpad('kfE5 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,170) = strpad('kbE5 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,171) = strpad('kfE7 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,172) = strpad('kbE7 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,173) = strpad('kfE10 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,174) = strpad('kbE10 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,175) = strpad('kfE11 in component KEL_DAGIP3_Model_Modules (per_um_s)');
LEGEND_CONSTANTS(:,176) = strpad('kbE11 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,177) = strpad('kmE9 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,178) = strpad('VmaxE9 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,179) = strpad('kfE13 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,180) = strpad('kbE13 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,181) = strpad('kmE2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,182) = strpad('kmE4 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,183) = strpad('kmE6 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,184) = strpad('kmE8 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,185) = strpad('kmE12 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,186) = strpad('VmaxE2 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,187) = strpad('VmaxE4 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,188) = strpad('VmaxE6 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,189) = strpad('VmaxE8 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,190) = strpad('VmaxE12 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,191) = strpad('kmF2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,192) = strpad('kmF6 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,193) = strpad('VmaxF2 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,194) = strpad('VmaxF6 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,195) = strpad('kfG8 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,196) = strpad('kfG9 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,197) = strpad('kmG6 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,198) = strpad('VmaxG6 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_CONSTANTS(:,199) = strpad('kmG7 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_CONSTANTS(:,200) = strpad('VmaxG7 in component KEL_DAGIP3_Model_Modules (per_s)');
LEGEND_RATES(:,13) = strpad('d/dt Ga_GDP in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,14) = strpad('d/dt mGluR in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,19) = strpad('d/dt Gabg in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,20) = strpad('d/dt R in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,12) = strpad('d/dt Gbg in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,21) = strpad('d/dt Gabg_R in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,22) = strpad('d/dt Glu_synapse in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,23) = strpad('d/dt Gabg_GluR in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,24) = strpad('d/dt Glu in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,25) = strpad('d/dt Ga_GTP in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,15) = strpad('d/dt PLC in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,26) = strpad('d/dt GqPLC in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,8) = strpad('d/dt CaPLC in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,7) = strpad('d/dt CaGqPLC in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,27) = strpad('d/dt GEF_inact in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,10) = strpad('d/dt GEF in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,5) = strpad('d/dt CaM_GEF in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,28) = strpad('d/dt Ca4CAM in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,6) = strpad('d/dt Gbg_GEF in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,4) = strpad('d/dt GEF_star in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,29) = strpad('d/dt CaM in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,9) = strpad('d/dt NgCaM in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,16) = strpad('d/dt Ca2CaM in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,17) = strpad('d/dt Ca3CaM in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,18) = strpad('d/dt Ng_star in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,11) = strpad('d/dt Ng in component DGBM_Model_Modules (um)');
LEGEND_RATES(:,38) = strpad('d/dt GTP_RAS in component B_Model_Modules (um)');
LEGEND_RATES(:,39) = strpad('d/dt GDP_RAS in component B_Model_Modules (um)');
LEGEND_RATES(:,2) = strpad('d/dt GAP in component B_Model_Modules (um)');
LEGEND_RATES(:,40) = strpad('d/dt GAPstar in component B_Model_Modules (um)');
LEGEND_RATES(:,36) = strpad('d/dt GTPRasRaf_star in component B_Model_Modules (um)');
LEGEND_RATES(:,3) = strpad('d/dt Raf in component B_Model_Modules (um)');
LEGEND_RATES(:,35) = strpad('d/dt Raf_star in component B_Model_Modules (um)');
LEGEND_RATES(:,37) = strpad('d/dt Raf_star_star in component B_Model_Modules (um)');
LEGEND_RATES(:,41) = strpad('d/dt MAPKK in component H_Model_Module_MAPKK (um)');
LEGEND_RATES(:,34) = strpad('d/dt MAPKK_star in component H_Model_Module_MAPKK (um)');
LEGEND_RATES(:,33) = strpad('d/dt MAPKK_star_star in component H_Model_Module_MAPKK (um)');
LEGEND_RATES(:,42) = strpad('d/dt MAPK in component H_Model_Module_MAPK (um)');
LEGEND_RATES(:,43) = strpad('d/dt MAPK_tyr in component H_Model_Module_MAPK (um)');
LEGEND_RATES(:,31) = strpad('d/dt MAPK_star in component H_Model_Module_MAPK (um)');
LEGEND_RATES(:,47) = strpad('d/dt SHCstar in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,46) = strpad('d/dt SHC in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,30) = strpad('d/dt SHCstar_SOS_GRB2 in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,49) = strpad('d/dt SOS_GRB2 in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,32) = strpad('d/dt SOS in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,51) = strpad('d/dt SOSstar in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,48) = strpad('d/dt SOSstar_GRB2 in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,50) = strpad('d/dt GRB2 in component A_Model_Module_SHC (um)');
LEGEND_RATES(:,53) = strpad('d/dt EGF in component A_Model_Module_EGFR (um)');
LEGEND_RATES(:,54) = strpad('d/dt EGFR in component A_Model_Module_EGFR (um)');
LEGEND_RATES(:,44) = strpad('d/dt EGF_EGFR in component A_Model_Module_EGFR (um)');
LEGEND_RATES(:,52) = strpad('d/dt EGF_EGFR_INTERNAL in component A_Model_Module_EGFR (um)');
LEGEND_RATES(:,56) = strpad('d/dt PLCg in component F_Model_Module_PLC (um)');
LEGEND_RATES(:,45) = strpad('d/dt CaPLCg in component F_Model_Module_PLC (um)');
LEGEND_RATES(:,55) = strpad('d/dt CaPLCg_star in component F_Model_Module_PLC (um)');
LEGEND_RATES(:,57) = strpad('d/dt PLCg_star in component F_Model_Module_PLC (um)');
LEGEND_RATES(:,58) = strpad('d/dt PKC_i in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,59) = strpad('d/dt DAGPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,60) = strpad('d/dt AADAGPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,62) = strpad('d/dt CaPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,61) = strpad('d/dt DAGCaPKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,1) = strpad('d/dt PKC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,63) = strpad('d/dt AA in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,66) = strpad('d/dt PLA2_cyt in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,70) = strpad('d/dt PIP2PLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,65) = strpad('d/dt PIP2_star in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,67) = strpad('d/dt PLA2_star in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,72) = strpad('d/dt CaPLA2_star in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,68) = strpad('d/dt CaPLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,69) = strpad('d/dt PIP2CaPLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,71) = strpad('d/dt DAGCaPLA2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,64) = strpad('d/dt APC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,73) = strpad('d/dt PIP2 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,74) = strpad('d/dt DAG in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,75) = strpad('d/dt IP3 in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,77) = strpad('d/dt PC in component KEL_DAGIP3_Model_Modules (um)');
LEGEND_RATES(:,76) = strpad('d/dt Inositol in component KEL_DAGIP3_Model_Modules (um)');
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) = 1;
CONSTANTS(:,2) = 0.5;
CONSTANTS(:,3) = 0.1;
STATES(:,1) = 0;
STATES(:,2) = 0.002;
STATES(:,3) = 0.2;
STATES(:,4) = 0;
STATES(:,5) = 0;
STATES(:,6) = 0;
STATES(:,7) = 0;
STATES(:,8) = 0;
STATES(:,9) = 0;
STATES(:,10) = 0;
STATES(:,11) = 10;
STATES(:,12) = 0.1;
STATES(:,13) = 1;
STATES(:,14) = 0.3;
STATES(:,15) = 0.8;
STATES(:,16) = 0;
STATES(:,17) = 0;
STATES(:,18) = 0;
STATES(:,19) = 0;
STATES(:,20) = 0;
STATES(:,21) = 0;
STATES(:,22) = 0;
STATES(:,23) = 0;
STATES(:,24) = 0;
STATES(:,25) = 0;
STATES(:,26) = 0;
STATES(:,27) = 0.1;
STATES(:,28) = 0;
STATES(:,29) = 20;
CONSTANTS(:,4) = 1;
CONSTANTS(:,5) = 500;
CONSTANTS(:,6) = 1000;
CONSTANTS(:,7) = 0.000028;
CONSTANTS(:,8) = 0.1;
CONSTANTS(:,9) = 0.000028;
CONSTANTS(:,10) = 10;
CONSTANTS(:,11) = 0.000001;
CONSTANTS(:,12) = 1;
CONSTANTS(:,13) = 0.00000001;
CONSTANTS(:,14) = 0.0001;
CONSTANTS(:,15) = 0.01;
CONSTANTS(:,16) = 0.0133;
CONSTANTS(:,17) = 0.00001;
CONSTANTS(:,18) = 0.0001;
CONSTANTS(:,19) = 0.000005;
CONSTANTS(:,20) = 1;
CONSTANTS(:,21) = 0.0000042;
CONSTANTS(:,22) = 1;
CONSTANTS(:,23) = 0.00005;
CONSTANTS(:,24) = 1;
CONSTANTS(:,25) = 0.0000042;
CONSTANTS(:,26) = 1;
CONSTANTS(:,27) = 0.0133;
CONSTANTS(:,28) = 1;
CONSTANTS(:,29) = 0;
CONSTANTS(:,30) = 3.333333333333;
CONSTANTS(:,31) = 0.0001;
CONSTANTS(:,32) = 1;
CONSTANTS(:,33) = 1;
CONSTANTS(:,34) = 0;
CONSTANTS(:,35) = 0.00001;
CONSTANTS(:,36) = 1;
CONSTANTS(:,37) = 7.5;
CONSTANTS(:,38) = 9;
CONSTANTS(:,39) = 3.33333333333;
CONSTANTS(:,40) = 4;
CONSTANTS(:,41) = 0.0000000002;
CONSTANTS(:,42) = 72;
CONSTANTS(:,43) = 0.000006;
CONSTANTS(:,44) = 10;
CONSTANTS(:,45) = 0.000000775;
CONSTANTS(:,46) = 10;
CONSTANTS(:,47) = 0.0000005;
CONSTANTS(:,48) = 1;
CONSTANTS(:,49) = 0.005;
CONSTANTS(:,50) = 10.012;
CONSTANTS(:,51) = 0.67;
CONSTANTS(:,52) = 28.626667;
CONSTANTS(:,53) = 0.58;
CONSTANTS(:,54) = 28.595;
CONSTANTS(:,55) = 0.35;
CONSTANTS(:,56) = 66.666666667;
STATES(:,30) = 0;
STATES(:,31) = 0;
STATES(:,32) = 0.1;
STATES(:,33) = 0;
STATES(:,34) = 0;
STATES(:,35) = 0;
STATES(:,36) = 0;
STATES(:,37) = 0;
STATES(:,38) = 0;
STATES(:,39) = 0.2;
STATES(:,40) = 0;
CONSTANTS(:,57) = 0.224;
CONSTANTS(:,58) = 1;
CONSTANTS(:,59) = 0.0001;
CONSTANTS(:,60) = 0;
CONSTANTS(:,61) = 0.50505;
CONSTANTS(:,62) = 0.02;
CONSTANTS(:,63) = 0.50505;
CONSTANTS(:,64) = 0.02;
CONSTANTS(:,65) = 0.50505;
CONSTANTS(:,66) = 0.02;
CONSTANTS(:,67) = 1.0104;
CONSTANTS(:,68) = 10;
CONSTANTS(:,69) = 0.1;
CONSTANTS(:,70) = 0;
CONSTANTS(:,71) = 3.33333333333;
CONSTANTS(:,72) = 4;
CONSTANTS(:,73) = 0.50505;
CONSTANTS(:,74) = 0.02;
CONSTANTS(:,75) = 0.00004;
CONSTANTS(:,76) = 0.5;
CONSTANTS(:,77) = 66.666666667;
CONSTANTS(:,78) = 25.64166667;
CONSTANTS(:,79) = 15.6565;
CONSTANTS(:,80) = 15.6565;
CONSTANTS(:,81) = 4;
CONSTANTS(:,82) = 10;
CONSTANTS(:,83) = 6;
CONSTANTS(:,84) = 6;
CONSTANTS(:,85) = 28.595;
CONSTANTS(:,86) = 3.33333333333;
CONSTANTS(:,87) = 28.626667;
CONSTANTS(:,88) = 2.564166667;
CONSTANTS(:,89) = 15.6565;
CONSTANTS(:,90) = 15.6565;
STATES(:,41) = 0.18;
CONSTANTS(:,91) = 0.224;
CONSTANTS(:,92) = 1;
CONSTANTS(:,93) = 0.159091667;
CONSTANTS(:,94) = 0.159091667;
CONSTANTS(:,95) = 15.6565;
CONSTANTS(:,96) = 15.6565;
CONSTANTS(:,97) = 0.105;
CONSTANTS(:,98) = 0.105;
CONSTANTS(:,99) = 6;
CONSTANTS(:,100) = 6;
CONSTANTS(:,101) = 15.6565;
CONSTANTS(:,102) = 15.6565;
STATES(:,42) = 0.36;
STATES(:,43) = 0;
CONSTANTS(:,103) = 0.032;
CONSTANTS(:,104) = 1;
CONSTANTS(:,105) = 0.046296667;
CONSTANTS(:,106) = 0.046296667;
CONSTANTS(:,107) = 0.066666667;
CONSTANTS(:,108) = 0.066666667;
CONSTANTS(:,109) = 0.15;
CONSTANTS(:,110) = 0.15;
CONSTANTS(:,111) = 1;
CONSTANTS(:,112) = 1;
STATES(:,44) = 0;
STATES(:,45) = 0;
STATES(:,46) = 0.5;
STATES(:,47) = 0;
STATES(:,48) = 0;
STATES(:,49) = 0;
STATES(:,50) = 1;
STATES(:,51) = 0;
CONSTANTS(:,113) = 1;
CONSTANTS(:,114) = 0.83333333333;
CONSTANTS(:,115) = 0.0016667;
CONSTANTS(:,116) = 0;
CONSTANTS(:,117) = 0.000000041667;
CONSTANTS(:,118) = 0.0168;
CONSTANTS(:,119) = 0.0000008333;
CONSTANTS(:,120) = 0.1;
CONSTANTS(:,121) = 0.001;
CONSTANTS(:,122) = 0;
CONSTANTS(:,123) = 0.000000041667;
CONSTANTS(:,124) = 0.0168;
CONSTANTS(:,125) = 0.40404;
CONSTANTS(:,126) = 2.564166667;
CONSTANTS(:,127) = 0.02;
CONSTANTS(:,128) = 0.2;
CONSTANTS(:,129) = 10;
CONSTANTS(:,130) = 25.64166667;
CONSTANTS(:,131) = 0.33333333333;
STATES(:,52) = 0;
STATES(:,53) = 166.67;
STATES(:,54) = 0.16667;
CONSTANTS(:,132) = 0.000007;
CONSTANTS(:,133) = 0.25;
CONSTANTS(:,134) = 0.002;
CONSTANTS(:,135) = 0.00033;
STATES(:,55) = 0;
STATES(:,56) = 0.82;
STATES(:,57) = 0;
CONSTANTS(:,136) = 0.0003;
CONSTANTS(:,137) = 10;
CONSTANTS(:,138) = 0.05;
CONSTANTS(:,139) = 0.00002;
CONSTANTS(:,140) = 10;
CONSTANTS(:,141) = 0.33333333333;
CONSTANTS(:,142) = 0.2;
CONSTANTS(:,143) = 0.83333333333;
CONSTANTS(:,144) = 1;
STATES(:,58) = 1;
STATES(:,59) = 0;
STATES(:,60) = 0;
STATES(:,61) = 0;
STATES(:,62) = 0;
STATES(:,63) = 0;
STATES(:,64) = 30;
STATES(:,65) = 2.5;
STATES(:,66) = 0.4;
STATES(:,67) = 0;
STATES(:,68) = 0;
STATES(:,69) = 0;
STATES(:,70) = 0;
STATES(:,71) = 0;
STATES(:,72) = 0;
STATES(:,73) = 10;
STATES(:,74) = 0;
STATES(:,75) = 0;
STATES(:,76) = 0;
STATES(:,77) = 0;
CONSTANTS(:,145) = 1;
CONSTANTS(:,146) = 50;
CONSTANTS(:,147) = 0.0000000002;
CONSTANTS(:,148) = 0.1;
CONSTANTS(:,149) = 1.2705;
CONSTANTS(:,150) = 3.5026;
CONSTANTS(:,151) = 0.000000002;
CONSTANTS(:,152) = 0.1;
CONSTANTS(:,153) = 1;
CONSTANTS(:,154) = 0.1;
CONSTANTS(:,155) = 2;
CONSTANTS(:,156) = 0.2;
CONSTANTS(:,157) = 0.000001;
CONSTANTS(:,158) = 0.5;
CONSTANTS(:,159) = 0.000000013333;
CONSTANTS(:,160) = 8.6348;
CONSTANTS(:,161) = 0.000000001;
CONSTANTS(:,162) = 0.1;
CONSTANTS(:,163) = 0.00000003;
CONSTANTS(:,164) = 2;
CONSTANTS(:,165) = 0.000000002;
CONSTANTS(:,166) = 0.5;
CONSTANTS(:,167) = 0.00000016667;
CONSTANTS(:,168) = 0.1;
CONSTANTS(:,169) = 0.00000002;
CONSTANTS(:,170) = 0.1;
CONSTANTS(:,171) = 0.000000005;
CONSTANTS(:,172) = 4;
CONSTANTS(:,173) = 0.17;
CONSTANTS(:,174) = 0;
CONSTANTS(:,175) = 0.00001;
CONSTANTS(:,176) = 0.1;
CONSTANTS(:,177) = 25.64166667;
CONSTANTS(:,178) = 20;
CONSTANTS(:,179) = 0.4;
CONSTANTS(:,180) = 0;
CONSTANTS(:,181) = 20;
CONSTANTS(:,182) = 20;
CONSTANTS(:,183) = 20;
CONSTANTS(:,184) = 20;
CONSTANTS(:,185) = 20;
CONSTANTS(:,186) = 11.04;
CONSTANTS(:,187) = 5.4;
CONSTANTS(:,188) = 36;
CONSTANTS(:,189) = 60;
CONSTANTS(:,190) = 120;
CONSTANTS(:,191) = 97;
CONSTANTS(:,192) = 19.79166667;
CONSTANTS(:,193) = 14;
CONSTANTS(:,194) = 57;
CONSTANTS(:,195) = 0.15;
CONSTANTS(:,196) = 2.5;
CONSTANTS(:,197) = 5;
CONSTANTS(:,198) = 48;
CONSTANTS(:,199) = 19.84166667;
CONSTANTS(:,200) = 10;
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(:,13) = ( STATES(:,25).*CONSTANTS(:,16) - STATES(:,13).*CONSTANTS(:,17).*STATES(:,12))+ STATES(:,7).*CONSTANTS(:,27);
RATES(:,14) = (( CONSTANTS(:,13).*STATES(:,23) - CONSTANTS(:,14).*STATES(:,19).*STATES(:,14)) - ( CONSTANTS(:,9).*STATES(:,14) - CONSTANTS(:,10).*STATES(:,24).*STATES(:,20)))+ CONSTANTS(:,15).*STATES(:,23);
RATES(:,19) = ((( STATES(:,13).*STATES(:,12).*CONSTANTS(:,17) - ( STATES(:,19).*STATES(:,20).*CONSTANTS(:,11) - STATES(:,21).*CONSTANTS(:,12)))+ CONSTANTS(:,13).*STATES(:,23)) - CONSTANTS(:,14).*STATES(:,19).*STATES(:,14)) - CONSTANTS(:,18).*STATES(:,19).*CONSTANTS(:,2);
RATES(:,20) = ( CONSTANTS(:,9).*STATES(:,14) - CONSTANTS(:,10).*STATES(:,24).*STATES(:,20)) - ( STATES(:,19).*STATES(:,20).*CONSTANTS(:,11) - STATES(:,21).*CONSTANTS(:,12));
RATES(:,12) = ( - STATES(:,13).*STATES(:,12).*CONSTANTS(:,17)+ CONSTANTS(:,15).*STATES(:,23)+ CONSTANTS(:,18).*STATES(:,19).*CONSTANTS(:,2)) - ( CONSTANTS(:,35).*STATES(:,10).*STATES(:,12) - CONSTANTS(:,36).*STATES(:,6));
RATES(:,21) = (( STATES(:,19).*STATES(:,20).*CONSTANTS(:,11) - STATES(:,21).*CONSTANTS(:,12))+ CONSTANTS(:,7).*STATES(:,23)) - STATES(:,21).*STATES(:,24).*CONSTANTS(:,8);
RATES(:,22) = - ( CONSTANTS(:,5).*STATES(:,22) - CONSTANTS(:,6).*STATES(:,24));
RATES(:,23) = ( - ( CONSTANTS(:,7).*STATES(:,23) - STATES(:,21).*STATES(:,24).*CONSTANTS(:,8)) - ( CONSTANTS(:,13).*STATES(:,23) - CONSTANTS(:,14).*STATES(:,19).*STATES(:,14))) - CONSTANTS(:,15).*STATES(:,23);
RATES(:,24) = (((( CONSTANTS(:,5).*STATES(:,22) - CONSTANTS(:,6).*STATES(:,24))+ CONSTANTS(:,7).*STATES(:,23)) - STATES(:,21).*STATES(:,24).*CONSTANTS(:,8))+ CONSTANTS(:,9).*STATES(:,14)) - CONSTANTS(:,10).*STATES(:,24).*STATES(:,20);
RATES(:,25) = (( CONSTANTS(:,18).*STATES(:,19).*CONSTANTS(:,2) - STATES(:,25).*CONSTANTS(:,16)) - ( STATES(:,15).*STATES(:,25).*CONSTANTS(:,21) - STATES(:,26).*CONSTANTS(:,22))) - ( STATES(:,7).*CONSTANTS(:,26)+ STATES(:,8).*STATES(:,25).*CONSTANTS(:,25));
RATES(:,15) = - ( STATES(:,15).*CONSTANTS(:,1).*CONSTANTS(:,19) - STATES(:,8).*CONSTANTS(:,20)) - ( STATES(:,15).*STATES(:,25).*CONSTANTS(:,21) - STATES(:,26).*CONSTANTS(:,22));
RATES(:,26) = ( STATES(:,15).*STATES(:,25).*CONSTANTS(:,21) - STATES(:,26).*CONSTANTS(:,22)) - ( CONSTANTS(:,1).*STATES(:,26).*CONSTANTS(:,23) - STATES(:,7).*CONSTANTS(:,24));
RATES(:,8) = (( STATES(:,15).*CONSTANTS(:,1).*CONSTANTS(:,19) - STATES(:,8).*CONSTANTS(:,20))+ STATES(:,7).*CONSTANTS(:,27)) - ( - STATES(:,7).*CONSTANTS(:,26)+ STATES(:,8).*STATES(:,25).*CONSTANTS(:,25));
RATES(:,7) = ((( STATES(:,26).*CONSTANTS(:,1).*CONSTANTS(:,23) - STATES(:,7).*CONSTANTS(:,24))+ STATES(:,8).*STATES(:,25).*CONSTANTS(:,25)) - STATES(:,7).*CONSTANTS(:,26)) - STATES(:,7).*CONSTANTS(:,27);
RATES(:,27) = ( CONSTANTS(:,3).*CONSTANTS(:,38).*STATES(:,10))./(CONSTANTS(:,37)+STATES(:,10)) - CONSTANTS(:,28).*STATES(:,27);
RATES(:,10) = ((( - (((( CONSTANTS(:,3).*CONSTANTS(:,38).*STATES(:,10))./(CONSTANTS(:,37)+STATES(:,10))+( STATES(:,1).*CONSTANTS(:,40).*STATES(:,10))./( CONSTANTS(:,39).*(CONSTANTS(:,4)+STATES(:,2)./CONSTANTS(:,30)+STATES(:,11)./CONSTANTS(:,52)+STATES(:,3)./CONSTANTS(:,56)+STATES(:,9)./CONSTANTS(:,54))+STATES(:,10))) - CONSTANTS(:,28).*STATES(:,27)) - CONSTANTS(:,33).*STATES(:,4)) - CONSTANTS(:,31).*STATES(:,28).*STATES(:,10))+ CONSTANTS(:,32).*STATES(:,5)) - CONSTANTS(:,35).*STATES(:,10).*STATES(:,12))+ CONSTANTS(:,36).*STATES(:,6);
RATES(:,5) = - CONSTANTS(:,32).*STATES(:,5)+ CONSTANTS(:,31).*STATES(:,10).*STATES(:,28);
RATES(:,28) = ( - ( - CONSTANTS(:,32).*STATES(:,5)+ CONSTANTS(:,31).*STATES(:,10).*STATES(:,28))+ CONSTANTS(:,1).*STATES(:,17).*CONSTANTS(:,45)) - CONSTANTS(:,46).*STATES(:,28);
RATES(:,6) = - ( - CONSTANTS(:,35).*STATES(:,10).*STATES(:,12)+ CONSTANTS(:,36).*STATES(:,6));
RATES(:,4) = ( STATES(:,1).*CONSTANTS(:,40).*STATES(:,10))./( CONSTANTS(:,39).*(CONSTANTS(:,4)+STATES(:,2)./CONSTANTS(:,30)+STATES(:,11)./CONSTANTS(:,52)+STATES(:,3)./CONSTANTS(:,56)+STATES(:,9)./CONSTANTS(:,54))+STATES(:,10)) - CONSTANTS(:,33).*STATES(:,4);
RATES(:,29) = ( - ( CONSTANTS(:,41).*CONSTANTS(:,1).*CONSTANTS(:,1).*STATES(:,29) - CONSTANTS(:,42).*STATES(:,16)) - ( CONSTANTS(:,47).*STATES(:,11).*STATES(:,29) - CONSTANTS(:,48).*STATES(:,9)))+( STATES(:,9).*STATES(:,1).*CONSTANTS(:,55))./(STATES(:,9)+CONSTANTS(:,54));
RATES(:,9) = ( CONSTANTS(:,47).*STATES(:,11).*STATES(:,29) - CONSTANTS(:,48).*STATES(:,9)) - ( STATES(:,9).*STATES(:,1).*CONSTANTS(:,55))./(STATES(:,9)+ CONSTANTS(:,54).*(CONSTANTS(:,4)+STATES(:,2)./CONSTANTS(:,30)+STATES(:,10)./CONSTANTS(:,39)+STATES(:,3)./CONSTANTS(:,56)+STATES(:,11)./CONSTANTS(:,52)));
RATES(:,16) = ( CONSTANTS(:,41).*CONSTANTS(:,1).*CONSTANTS(:,1).*STATES(:,29) - CONSTANTS(:,42).*STATES(:,16)) - ( CONSTANTS(:,1).*STATES(:,16).*CONSTANTS(:,43) - CONSTANTS(:,44).*STATES(:,17));
RATES(:,17) = ( CONSTANTS(:,1).*STATES(:,16).*CONSTANTS(:,43) - CONSTANTS(:,44).*STATES(:,17)) - ( CONSTANTS(:,1).*STATES(:,17).*CONSTANTS(:,45) - CONSTANTS(:,46).*STATES(:,28));
RATES(:,18) = ((( STATES(:,9).*STATES(:,1).*CONSTANTS(:,55))./(STATES(:,9)+ CONSTANTS(:,54).*(CONSTANTS(:,4)+STATES(:,2)./CONSTANTS(:,30)+STATES(:,10)./CONSTANTS(:,39)+STATES(:,3)./CONSTANTS(:,56)+STATES(:,11)./CONSTANTS(:,52))) - CONSTANTS(:,49).*STATES(:,18)) - ( CONSTANTS(:,51).*STATES(:,18))./(STATES(:,18)+CONSTANTS(:,50)))+( STATES(:,11).*STATES(:,1).*CONSTANTS(:,53))./(STATES(:,11)+ CONSTANTS(:,52).*(CONSTANTS(:,4)+STATES(:,2)./CONSTANTS(:,30)+STATES(:,10)./CONSTANTS(:,39)+STATES(:,3)./CONSTANTS(:,56)+STATES(:,9)./CONSTANTS(:,54)));
RATES(:,11) = ( CONSTANTS(:,49).*STATES(:,18)+( CONSTANTS(:,51).*STATES(:,18))./(STATES(:,18)+CONSTANTS(:,50))) - ( STATES(:,11).*STATES(:,1).*CONSTANTS(:,53))./(STATES(:,11)+ CONSTANTS(:,52).*(CONSTANTS(:,4)+STATES(:,2)./CONSTANTS(:,30)+STATES(:,10)./CONSTANTS(:,39)+STATES(:,3)./CONSTANTS(:,56)+STATES(:,9)./CONSTANTS(:,54)));
RATES(:,38) = (((( - STATES(:,2).*CONSTANTS(:,68).*STATES(:,38))./(CONSTANTS(:,67)+STATES(:,38)) - CONSTANTS(:,59).*STATES(:,38))+( STATES(:,4).*CONSTANTS(:,64).*STATES(:,39))./(CONSTANTS(:,63)+STATES(:,39))+( STATES(:,6).*CONSTANTS(:,66).*STATES(:,39))./(CONSTANTS(:,65)+STATES(:,39))+( STATES(:,5).*CONSTANTS(:,62).*STATES(:,39))./(CONSTANTS(:,61)+STATES(:,39))+( STATES(:,30).*CONSTANTS(:,74).*STATES(:,39))./(CONSTANTS(:,73)+STATES(:,39))) - STATES(:,38).*STATES(:,35).*CONSTANTS(:,75))+ STATES(:,36).*CONSTANTS(:,76);
RATES(:,39) = - ((( - STATES(:,2).*CONSTANTS(:,68).*STATES(:,38))./(CONSTANTS(:,67)+STATES(:,38)) - CONSTANTS(:,59).*STATES(:,38))+( STATES(:,4).*CONSTANTS(:,64).*STATES(:,39))./(CONSTANTS(:,63)+STATES(:,39))+( STATES(:,6).*CONSTANTS(:,66).*STATES(:,39))./(CONSTANTS(:,65)+STATES(:,39))+( STATES(:,5).*CONSTANTS(:,62).*STATES(:,39))./(CONSTANTS(:,61)+STATES(:,39))) - ( STATES(:,30).*CONSTANTS(:,74).*STATES(:,39))./(CONSTANTS(:,73)+STATES(:,39));
RATES(:,2) = ( - STATES(:,1).*CONSTANTS(:,72).*STATES(:,2))./( CONSTANTS(:,71).*(CONSTANTS(:,58)+STATES(:,3)./CONSTANTS(:,77)+STATES(:,10)./CONSTANTS(:,86)+STATES(:,11)./CONSTANTS(:,87)+STATES(:,9)./CONSTANTS(:,85))+STATES(:,2))+ CONSTANTS(:,69).*STATES(:,40);
RATES(:,40) = ( STATES(:,1).*CONSTANTS(:,72).*STATES(:,2))./( CONSTANTS(:,71).*(CONSTANTS(:,58)+STATES(:,3)./CONSTANTS(:,77)+STATES(:,10)./CONSTANTS(:,86)+STATES(:,11)./CONSTANTS(:,87)+STATES(:,9)./CONSTANTS(:,85))+STATES(:,2)) - CONSTANTS(:,69).*STATES(:,40);
RATES(:,36) = STATES(:,38).*STATES(:,35).*CONSTANTS(:,75) - STATES(:,36).*CONSTANTS(:,76);
RATES(:,3) = ( - STATES(:,1).*STATES(:,3).*CONSTANTS(:,81))./( CONSTANTS(:,77).*(CONSTANTS(:,58)+STATES(:,2)./CONSTANTS(:,71)+STATES(:,10)./CONSTANTS(:,86)+STATES(:,11)./CONSTANTS(:,87)+STATES(:,9)./CONSTANTS(:,85))+STATES(:,3))+( STATES(:,35).*CONSTANTS(:,57).*CONSTANTS(:,83))./( CONSTANTS(:,79).*(CONSTANTS(:,58)+STATES(:,37)./CONSTANTS(:,80)+STATES(:,33)./CONSTANTS(:,90)+STATES(:,34)./CONSTANTS(:,89))+STATES(:,35));
RATES(:,35) = ((( STATES(:,1).*STATES(:,3).*CONSTANTS(:,81))./( CONSTANTS(:,77).*(CONSTANTS(:,58)+STATES(:,2)./CONSTANTS(:,71)+STATES(:,10)./CONSTANTS(:,86)+STATES(:,11)./CONSTANTS(:,87)+STATES(:,9)./CONSTANTS(:,85))+STATES(:,3)) - ( STATES(:,35).*CONSTANTS(:,57).*CONSTANTS(:,83))./( CONSTANTS(:,79).*(CONSTANTS(:,58)+STATES(:,37)./CONSTANTS(:,80)+STATES(:,33)./CONSTANTS(:,90)+STATES(:,34)./CONSTANTS(:,89))+STATES(:,35))) - ( STATES(:,35).*STATES(:,31).*CONSTANTS(:,82))./( CONSTANTS(:,78).*(CONSTANTS(:,58)+STATES(:,32)./CONSTANTS(:,88))+STATES(:,35)))+( STATES(:,37).*CONSTANTS(:,57).*CONSTANTS(:,84))./( CONSTANTS(:,80).*(CONSTANTS(:,58)+STATES(:,35)./CONSTANTS(:,79)+STATES(:,33)./CONSTANTS(:,90)+STATES(:,34)./CONSTANTS(:,89))+STATES(:,37));
RATES(:,37) = ( STATES(:,35).*STATES(:,31).*CONSTANTS(:,82))./( CONSTANTS(:,78).*(CONSTANTS(:,58)+STATES(:,32)./CONSTANTS(:,88))+STATES(:,35)) - ( STATES(:,37).*CONSTANTS(:,57).*CONSTANTS(:,84))./( CONSTANTS(:,80).*(CONSTANTS(:,58)+STATES(:,35)./CONSTANTS(:,79)+STATES(:,33)./CONSTANTS(:,90)+STATES(:,34)./CONSTANTS(:,89))+STATES(:,37));
RATES(:,41) = ( - STATES(:,41).*CONSTANTS(:,97).*STATES(:,36))./( CONSTANTS(:,93).*(CONSTANTS(:,92)+STATES(:,34)./CONSTANTS(:,94))+STATES(:,41))+( STATES(:,34).*CONSTANTS(:,99).*CONSTANTS(:,91))./( CONSTANTS(:,95).*(CONSTANTS(:,92)+STATES(:,35)./CONSTANTS(:,101)+STATES(:,37)./CONSTANTS(:,102)+STATES(:,33)./CONSTANTS(:,96))+STATES(:,34));
RATES(:,34) = ((( STATES(:,41).*CONSTANTS(:,97).*STATES(:,36))./( CONSTANTS(:,93).*(CONSTANTS(:,92)+STATES(:,34)./CONSTANTS(:,94))+STATES(:,41)) - ( STATES(:,34).*CONSTANTS(:,99).*CONSTANTS(:,91))./( CONSTANTS(:,95).*(CONSTANTS(:,92)+STATES(:,35)./CONSTANTS(:,101)+STATES(:,37)./CONSTANTS(:,102)+STATES(:,33)./CONSTANTS(:,96))+STATES(:,34))) - ( STATES(:,34).*CONSTANTS(:,98).*STATES(:,36))./( CONSTANTS(:,94).*(CONSTANTS(:,92)+STATES(:,41)./CONSTANTS(:,93))+STATES(:,34)))+( STATES(:,33).*CONSTANTS(:,100).*CONSTANTS(:,91))./( CONSTANTS(:,96).*(CONSTANTS(:,92)+STATES(:,35)./CONSTANTS(:,101)+STATES(:,37)./CONSTANTS(:,102)+STATES(:,34)./CONSTANTS(:,95))+STATES(:,33));
RATES(:,33) = ( STATES(:,34).*CONSTANTS(:,98).*STATES(:,36))./( CONSTANTS(:,94).*(CONSTANTS(:,92)+STATES(:,41)./CONSTANTS(:,93))+STATES(:,34)) - ( STATES(:,33).*CONSTANTS(:,100).*CONSTANTS(:,91))./( CONSTANTS(:,96).*(CONSTANTS(:,92)+STATES(:,35)./CONSTANTS(:,101)+STATES(:,37)./CONSTANTS(:,102)+STATES(:,34)./CONSTANTS(:,95))+STATES(:,33));
RATES(:,42) = ( - STATES(:,42).*CONSTANTS(:,109).*STATES(:,33))./( CONSTANTS(:,105).*(CONSTANTS(:,104)+STATES(:,43)./CONSTANTS(:,106))+STATES(:,42))+( STATES(:,43).*CONSTANTS(:,111).*CONSTANTS(:,103))./( CONSTANTS(:,107).*(CONSTANTS(:,104)+STATES(:,31)./CONSTANTS(:,108))+STATES(:,43));
RATES(:,43) = ((( STATES(:,42).*CONSTANTS(:,109).*STATES(:,33))./(( CONSTANTS(:,105).*STATES(:,43))./CONSTANTS(:,106)+STATES(:,42)) - ( STATES(:,43).*CONSTANTS(:,111).*CONSTANTS(:,103))./( CONSTANTS(:,107).*(CONSTANTS(:,104)+STATES(:,31)./CONSTANTS(:,108))+STATES(:,43))) - ( STATES(:,43).*CONSTANTS(:,110).*STATES(:,33))./( CONSTANTS(:,106).*(CONSTANTS(:,104)+STATES(:,43)./CONSTANTS(:,107))+STATES(:,43)))+( STATES(:,31).*CONSTANTS(:,112).*CONSTANTS(:,103))./( CONSTANTS(:,108).*(CONSTANTS(:,104)+STATES(:,43)./CONSTANTS(:,107))+STATES(:,31));
RATES(:,31) = ( STATES(:,43).*CONSTANTS(:,110).*STATES(:,33))./( CONSTANTS(:,106).*(CONSTANTS(:,104)+STATES(:,42)./CONSTANTS(:,105))+STATES(:,43)) - ( STATES(:,31).*CONSTANTS(:,112).*CONSTANTS(:,103))./( CONSTANTS(:,108).*(CONSTANTS(:,104)+STATES(:,43)./CONSTANTS(:,107))+STATES(:,31));
RATES(:,47) = ((( STATES(:,44).*CONSTANTS(:,128).*STATES(:,46))./( CONSTANTS(:,114).*(CONSTANTS(:,113)+STATES(:,45)./CONSTANTS(:,131))+STATES(:,46)) - CONSTANTS(:,119).*STATES(:,47).*STATES(:,49))+ CONSTANTS(:,120).*STATES(:,30)) - CONSTANTS(:,115).*STATES(:,47);
RATES(:,46) = ( CONSTANTS(:,115).*STATES(:,47) - CONSTANTS(:,116).*STATES(:,46)) - ( STATES(:,44).*CONSTANTS(:,128).*STATES(:,46))./( CONSTANTS(:,114).*(CONSTANTS(:,113)+STATES(:,45)./CONSTANTS(:,131))+STATES(:,46));
RATES(:,30) = CONSTANTS(:,119).*STATES(:,47).*STATES(:,49) - CONSTANTS(:,120).*STATES(:,30);
RATES(:,49) = (( - CONSTANTS(:,119).*STATES(:,47).*STATES(:,49)+ CONSTANTS(:,120).*STATES(:,30)) - CONSTANTS(:,118).*STATES(:,49))+ CONSTANTS(:,117).*STATES(:,32).*STATES(:,50);
RATES(:,32) = ((( CONSTANTS(:,121).*STATES(:,51) - CONSTANTS(:,122).*STATES(:,32)) - ( STATES(:,31).*CONSTANTS(:,129).*STATES(:,32))./( CONSTANTS(:,126).*(CONSTANTS(:,113)+STATES(:,35)./CONSTANTS(:,130))+STATES(:,32)))+ CONSTANTS(:,118).*STATES(:,49)) - CONSTANTS(:,117).*STATES(:,32).*STATES(:,50);
RATES(:,51) = (( - CONSTANTS(:,121).*STATES(:,51)+ CONSTANTS(:,122).*STATES(:,32)+( STATES(:,31).*CONSTANTS(:,129).*STATES(:,32))./( CONSTANTS(:,126).*(CONSTANTS(:,113)+STATES(:,35)./CONSTANTS(:,130))+STATES(:,32))) - CONSTANTS(:,123).*STATES(:,51).*STATES(:,50))+ CONSTANTS(:,124).*STATES(:,48);
RATES(:,48) = CONSTANTS(:,123).*STATES(:,51).*STATES(:,50) - CONSTANTS(:,124).*STATES(:,48);
RATES(:,50) = (( CONSTANTS(:,118).*STATES(:,49) - CONSTANTS(:,117).*STATES(:,32).*STATES(:,50)) - CONSTANTS(:,123).*STATES(:,51).*STATES(:,50))+ CONSTANTS(:,124).*STATES(:,48);
RATES(:,53) = - CONSTANTS(:,132).*STATES(:,53).*STATES(:,54)+ CONSTANTS(:,133).*STATES(:,44);
RATES(:,54) = - CONSTANTS(:,132).*STATES(:,53).*STATES(:,54)+ CONSTANTS(:,133).*STATES(:,44);
RATES(:,44) = (( CONSTANTS(:,132).*STATES(:,53).*STATES(:,54) - CONSTANTS(:,133).*STATES(:,44)) - CONSTANTS(:,134).*STATES(:,44))+ CONSTANTS(:,135).*STATES(:,52);
RATES(:,52) = CONSTANTS(:,134).*STATES(:,44) - CONSTANTS(:,135).*STATES(:,52);
RATES(:,56) = - CONSTANTS(:,1).*STATES(:,56).*CONSTANTS(:,136)+ STATES(:,45).*CONSTANTS(:,137);
RATES(:,45) = ( - ( - CONSTANTS(:,1).*STATES(:,56).*CONSTANTS(:,136)+ STATES(:,45).*CONSTANTS(:,137))+ CONSTANTS(:,138).*STATES(:,55)) - ( CONSTANTS(:,142).*STATES(:,44).*STATES(:,45))./( CONSTANTS(:,141).*(CONSTANTS(:,144)+STATES(:,46)./CONSTANTS(:,143))+STATES(:,45));
RATES(:,55) = ( - ( CONSTANTS(:,138).*STATES(:,55) - ( CONSTANTS(:,142).*STATES(:,44).*STATES(:,45))./(CONSTANTS(:,141)+STATES(:,45)))+ CONSTANTS(:,1).*STATES(:,57).*CONSTANTS(:,139)) - CONSTANTS(:,140).*STATES(:,55);
RATES(:,57) = - CONSTANTS(:,1).*STATES(:,57).*CONSTANTS(:,139)+ CONSTANTS(:,140).*STATES(:,55);
RATES(:,58) = (((((( - STATES(:,58).*CONSTANTS(:,145)+ STATES(:,1).*CONSTANTS(:,146)) - STATES(:,58).*STATES(:,63).*CONSTANTS(:,147))+ STATES(:,1).*CONSTANTS(:,148)) - STATES(:,58).*STATES(:,74).*CONSTANTS(:,161))+ STATES(:,59).*CONSTANTS(:,162)) - STATES(:,58).*CONSTANTS(:,1).*CONSTANTS(:,157))+ STATES(:,62).*CONSTANTS(:,158);
RATES(:,59) = (( STATES(:,58).*STATES(:,74).*CONSTANTS(:,161) - STATES(:,59).*CONSTANTS(:,162)) - STATES(:,59).*STATES(:,63).*CONSTANTS(:,163))+ STATES(:,60).*CONSTANTS(:,164);
RATES(:,60) = (( STATES(:,59).*STATES(:,63).*CONSTANTS(:,163) - STATES(:,60).*CONSTANTS(:,164)) - STATES(:,60).*CONSTANTS(:,155))+ STATES(:,1).*CONSTANTS(:,156);
RATES(:,62) = (((((( STATES(:,58).*CONSTANTS(:,1).*CONSTANTS(:,157) - STATES(:,62).*CONSTANTS(:,158)) - STATES(:,62).*CONSTANTS(:,149))+ STATES(:,1).*CONSTANTS(:,150)) - STATES(:,62).*STATES(:,63).*CONSTANTS(:,151))+ STATES(:,1).*CONSTANTS(:,152)) - STATES(:,74).*STATES(:,62).*CONSTANTS(:,159))+ STATES(:,61).*CONSTANTS(:,160);
RATES(:,61) = (( STATES(:,74).*STATES(:,62).*CONSTANTS(:,159) - STATES(:,61).*CONSTANTS(:,160)) - STATES(:,61).*CONSTANTS(:,153))+ STATES(:,1).*CONSTANTS(:,154);
RATES(:,1) = (((((((((( STATES(:,58).*CONSTANTS(:,145) - STATES(:,1).*CONSTANTS(:,146))+ STATES(:,58).*STATES(:,63).*CONSTANTS(:,147)) - STATES(:,1).*CONSTANTS(:,148))+ STATES(:,62).*CONSTANTS(:,149)) - STATES(:,1).*CONSTANTS(:,150))+ STATES(:,62).*STATES(:,63).*CONSTANTS(:,151)) - STATES(:,1).*CONSTANTS(:,152))+ STATES(:,61).*CONSTANTS(:,153)) - STATES(:,1).*CONSTANTS(:,154))+ STATES(:,60).*CONSTANTS(:,155)) - STATES(:,1).*CONSTANTS(:,156);
RATES(:,63) = (((((( - STATES(:,58).*STATES(:,63).*CONSTANTS(:,147)+ STATES(:,1).*CONSTANTS(:,148)) - STATES(:,62).*STATES(:,63).*CONSTANTS(:,151))+ STATES(:,1).*CONSTANTS(:,152)) - STATES(:,59).*STATES(:,63).*CONSTANTS(:,163))+ STATES(:,60).*CONSTANTS(:,164)) - STATES(:,63).*CONSTANTS(:,179))+( STATES(:,72).*CONSTANTS(:,190).*STATES(:,64))./(CONSTANTS(:,185)+STATES(:,64))+( STATES(:,68).*CONSTANTS(:,187).*STATES(:,64))./(CONSTANTS(:,182)+STATES(:,64))+( STATES(:,71).*CONSTANTS(:,189).*STATES(:,64))./(CONSTANTS(:,184)+STATES(:,64))+( STATES(:,69).*CONSTANTS(:,188).*STATES(:,64))./(CONSTANTS(:,183)+STATES(:,64))+( STATES(:,70).*CONSTANTS(:,186).*STATES(:,64))./(CONSTANTS(:,181)+STATES(:,64));
RATES(:,66) = (((( STATES(:,67).*CONSTANTS(:,173) - STATES(:,65).*STATES(:,66).*CONSTANTS(:,165))+ STATES(:,70).*CONSTANTS(:,166)) - CONSTANTS(:,1).*STATES(:,66).*CONSTANTS(:,167))+ STATES(:,68).*CONSTANTS(:,168)) - ( STATES(:,42).*CONSTANTS(:,178).*STATES(:,66))./(CONSTANTS(:,177)+STATES(:,66));
RATES(:,70) = STATES(:,65).*STATES(:,66).*CONSTANTS(:,165) - STATES(:,70).*CONSTANTS(:,166);
RATES(:,65) = (( - STATES(:,65).*STATES(:,66).*CONSTANTS(:,165)+ STATES(:,70).*CONSTANTS(:,166)) - STATES(:,65).*STATES(:,68).*CONSTANTS(:,169))+ STATES(:,69).*CONSTANTS(:,170);
RATES(:,67) = (( - STATES(:,67).*CONSTANTS(:,173)+( STATES(:,42).*CONSTANTS(:,178).*STATES(:,66))./(CONSTANTS(:,177)+STATES(:,66))) - CONSTANTS(:,1).*STATES(:,67).*CONSTANTS(:,175))+ STATES(:,72).*CONSTANTS(:,176);
RATES(:,72) = CONSTANTS(:,1).*STATES(:,67).*CONSTANTS(:,175) - STATES(:,72).*CONSTANTS(:,176);
RATES(:,68) = (((( CONSTANTS(:,1).*STATES(:,66).*CONSTANTS(:,167) - STATES(:,68).*CONSTANTS(:,168)) - STATES(:,74).*STATES(:,68).*CONSTANTS(:,171))+ STATES(:,71).*CONSTANTS(:,172)) - STATES(:,65).*STATES(:,68).*CONSTANTS(:,169))+ STATES(:,69).*CONSTANTS(:,170);
RATES(:,69) = STATES(:,65).*STATES(:,68).*CONSTANTS(:,169) - STATES(:,69).*CONSTANTS(:,170);
RATES(:,71) = STATES(:,74).*STATES(:,68).*CONSTANTS(:,171) - STATES(:,71).*CONSTANTS(:,172);
RATES(:,64) = STATES(:,63).*CONSTANTS(:,179) - (( STATES(:,72).*CONSTANTS(:,190).*STATES(:,64))./(CONSTANTS(:,185)+STATES(:,64))+( STATES(:,68).*CONSTANTS(:,187).*STATES(:,64))./(CONSTANTS(:,182)+STATES(:,64))+( STATES(:,71).*CONSTANTS(:,189).*STATES(:,64))./(CONSTANTS(:,184)+STATES(:,64))+( STATES(:,69).*CONSTANTS(:,188).*STATES(:,64))./(CONSTANTS(:,183)+STATES(:,64))+( STATES(:,70).*CONSTANTS(:,186).*STATES(:,64))./(CONSTANTS(:,181)+STATES(:,64)));
RATES(:,73) = ((( - STATES(:,45).*CONSTANTS(:,193).*STATES(:,73))./(CONSTANTS(:,191)+STATES(:,73)) - ( STATES(:,55).*CONSTANTS(:,194).*STATES(:,73))./(CONSTANTS(:,192)+STATES(:,73))) - ( STATES(:,8).*CONSTANTS(:,200).*STATES(:,73))./(CONSTANTS(:,199)+STATES(:,73))) - ( STATES(:,7).*CONSTANTS(:,198).*STATES(:,73))./(CONSTANTS(:,197)+STATES(:,73));
RATES(:,74) = ((((((((( STATES(:,45).*CONSTANTS(:,193).*STATES(:,73))./(CONSTANTS(:,191)+STATES(:,73))+( STATES(:,55).*CONSTANTS(:,194).*STATES(:,73))./(CONSTANTS(:,192)+STATES(:,73))+( STATES(:,8).*CONSTANTS(:,200).*STATES(:,73))./(CONSTANTS(:,199)+STATES(:,73))) - CONSTANTS(:,195).*STATES(:,74))+( STATES(:,7).*CONSTANTS(:,198).*STATES(:,73))./(CONSTANTS(:,197)+STATES(:,73))) - STATES(:,58).*STATES(:,74).*CONSTANTS(:,161))+ STATES(:,59).*CONSTANTS(:,162)) - STATES(:,74).*STATES(:,62).*CONSTANTS(:,159))+ STATES(:,61).*CONSTANTS(:,160)) - STATES(:,74).*STATES(:,68).*CONSTANTS(:,171))+ STATES(:,71).*CONSTANTS(:,172);
RATES(:,75) = ((( STATES(:,45).*CONSTANTS(:,193).*STATES(:,73))./(CONSTANTS(:,191)+STATES(:,73))+( STATES(:,55).*CONSTANTS(:,194).*STATES(:,73))./(CONSTANTS(:,192)+STATES(:,73))+( STATES(:,7).*CONSTANTS(:,198).*STATES(:,73))./(CONSTANTS(:,197)+STATES(:,73))) - CONSTANTS(:,196).*STATES(:,75))+( STATES(:,8).*CONSTANTS(:,200).*STATES(:,73))./(CONSTANTS(:,199)+STATES(:,73));
RATES(:,77) = CONSTANTS(:,195).*STATES(:,74);
RATES(:,76) = CONSTANTS(:,196).*STATES(:,75);
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
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