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 66 entries in each of the rate and state variable arrays. % There are a total of 255 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 (hour)'); LEGEND_STATES(:,1) = strpad('AP1 in component AP1 (micromolar)'); LEGEND_CONSTANTS(:,1) = strpad('v_sap1 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,2) = strpad('GF in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,3) = strpad('K_agf in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,4) = strpad('k_dap1 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,5) = strpad('eps in component model_parameters (dimensionless)'); LEGEND_STATES(:,2) = strpad('pRB in component pRB (micromolar)'); LEGEND_CONSTANTS(:,6) = strpad('v_sprb in component model_parameters (flux)'); LEGEND_CONSTANTS(:,7) = strpad('k_pc1 in component model_parameters (second_order_rate_constant)'); LEGEND_STATES(:,3) = strpad('E2F in component E2F (micromolar)'); LEGEND_CONSTANTS(:,8) = strpad('k_pc2 in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,4) = strpad('pRBc1 in component pRBc1 (micromolar)'); LEGEND_CONSTANTS(:,9) = strpad('V_1 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,10) = strpad('K_1 in component model_parameters (micromolar)'); LEGEND_STATES(:,5) = strpad('Md in component Md (micromolar)'); LEGEND_STATES(:,6) = strpad('Mdp27 in component Mdp27 (micromolar)'); LEGEND_CONSTANTS(:,11) = strpad('V_2 in component model_parameters (flux)'); LEGEND_STATES(:,7) = strpad('pRBp in component pRBp (micromolar)'); LEGEND_CONSTANTS(:,12) = strpad('K_2 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,13) = strpad('k_dprb in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,14) = strpad('V_3 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,15) = strpad('K_3 in component model_parameters (micromolar)'); LEGEND_STATES(:,8) = strpad('Me in component Me (micromolar)'); LEGEND_CONSTANTS(:,16) = strpad('V_4 in component model_parameters (flux)'); LEGEND_STATES(:,9) = strpad('pRBpp in component pRBpp (micromolar)'); LEGEND_CONSTANTS(:,17) = strpad('K_4 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,18) = strpad('k_pc3 in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,19) = strpad('k_pc4 in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,10) = strpad('pRBc2 in component pRBc2 (micromolar)'); LEGEND_CONSTANTS(:,20) = strpad('k_dpRBp in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,21) = strpad('k_dpRBpp in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,22) = strpad('v_se2f in component model_parameters (flux)'); LEGEND_CONSTANTS(:,23) = strpad('V_1e2f in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,11) = strpad('Ma in component Ma (micromolar)'); LEGEND_CONSTANTS(:,24) = strpad('K_1e2f in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,25) = strpad('V_2e2f in component model_parameters (flux)'); LEGEND_STATES(:,12) = strpad('E2Fp in component E2Fp (micromolar)'); LEGEND_CONSTANTS(:,26) = strpad('K_2e2f in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,27) = strpad('k_de2f in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,28) = strpad('k_de2fp in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,13) = strpad('Cd in component Cd (micromolar)'); LEGEND_CONSTANTS(:,29) = strpad('k_cd1 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,30) = strpad('k_cd2 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,31) = strpad('K_i7 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,32) = strpad('K_i8 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,33) = strpad('k_com1 in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,34) = strpad('Cdk4_tot in component model_parameters (micromolar)'); LEGEND_STATES(:,14) = strpad('Mdi in component Mdi (micromolar)'); LEGEND_CONSTANTS(:,35) = strpad('k_decom1 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,36) = strpad('V_dd in component model_parameters (flux)'); LEGEND_CONSTANTS(:,37) = strpad('K_dd in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,38) = strpad('k_ddd in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,39) = strpad('V_m2d in component model_parameters (flux)'); LEGEND_CONSTANTS(:,40) = strpad('K_2d in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,41) = strpad('V_m1d in component model_parameters (flux)'); LEGEND_CONSTANTS(:,42) = strpad('K_1d in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,43) = strpad('k_c1 in component model_parameters (second_order_rate_constant)'); LEGEND_STATES(:,15) = strpad('p27 in component p27 (micromolar)'); LEGEND_CONSTANTS(:,44) = strpad('k_c2 in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,16) = strpad('Ce in component Ce (micromolar)'); LEGEND_CONSTANTS(:,45) = strpad('k_ce in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,46) = strpad('K_i9 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,47) = strpad('K_i10 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,48) = strpad('k_com2 in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,49) = strpad('Cdk2_tot in component model_parameters (micromolar)'); LEGEND_STATES(:,17) = strpad('Mei in component Mei (micromolar)'); LEGEND_STATES(:,18) = strpad('Mep27 in component Mep27 (micromolar)'); LEGEND_STATES(:,19) = strpad('Mai in component Mai (micromolar)'); LEGEND_STATES(:,20) = strpad('Map27 in component Map27 (micromolar)'); LEGEND_CONSTANTS(:,50) = strpad('k_decom2 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,51) = strpad('V_de in component model_parameters (flux)'); LEGEND_STATES(:,21) = strpad('Skp2 in component Skp2 (micromolar)'); LEGEND_CONSTANTS(:,52) = strpad('K_dceskp2 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,53) = strpad('K_de in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,54) = strpad('k_dde in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,55) = strpad('V_m2e in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,22) = strpad('Wee1 in component Wee1 (micromolar)'); LEGEND_CONSTANTS(:,56) = strpad('i_b1 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,57) = strpad('K_2e in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,58) = strpad('V_m1e in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,23) = strpad('Pe in component Pe (micromolar)'); LEGEND_CONSTANTS(:,59) = strpad('K_1e in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,60) = strpad('k_c3 in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,61) = strpad('k_c4 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,62) = strpad('v_sskp2 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,63) = strpad('V_dskp2 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,64) = strpad('K_dskp2 in component model_parameters (micromolar)'); LEGEND_STATES(:,24) = strpad('Cdh1a in component Cdh1a (micromolar)'); LEGEND_CONSTANTS(:,65) = strpad('K_cdh1 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,66) = strpad('k_ddskp2 in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,25) = strpad('Pei in component Pei (micromolar)'); LEGEND_CONSTANTS(:,67) = strpad('v_spei in component model_parameters (flux)'); LEGEND_CONSTANTS(:,68) = strpad('V_6e in component model_parameters (flux)'); LEGEND_CONSTANTS(:,69) = strpad('x_e1 in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,70) = strpad('x_e2 in component model_parameters (per_micromolar)'); LEGEND_STATES(:,26) = strpad('Chk1 in component Chk1 (micromolar)'); LEGEND_CONSTANTS(:,71) = strpad('K_6e in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,72) = strpad('V_m5e in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,73) = strpad('a_e in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,74) = strpad('K_5e in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,75) = strpad('k_dpei in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,76) = strpad('k_dpe in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,27) = strpad('Ca in component Ca (micromolar)'); LEGEND_CONSTANTS(:,77) = strpad('k_ca in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,78) = strpad('K_i11 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,79) = strpad('K_i12 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,80) = strpad('k_com3 in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,81) = strpad('k_decom3 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,82) = strpad('V_da in component model_parameters (flux)'); LEGEND_CONSTANTS(:,83) = strpad('K_da in component model_parameters (micromolar)'); LEGEND_STATES(:,28) = strpad('Cdc20a in component Cdc20a (micromolar)'); LEGEND_CONSTANTS(:,84) = strpad('K_acdc20 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,85) = strpad('k_dda in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,86) = strpad('V_m2a in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,87) = strpad('i_b2 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,88) = strpad('K_2a in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,89) = strpad('V_m1a in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,29) = strpad('Pa in component Pa (micromolar)'); LEGEND_CONSTANTS(:,90) = strpad('K_1a in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,91) = strpad('k_c5 in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,92) = strpad('k_c6 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,93) = strpad('v_s1p27 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,94) = strpad('v_s2p27 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,95) = strpad('K_i13 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,96) = strpad('K_i14 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,97) = strpad('k_c7 in component model_parameters (second_order_rate_constant)'); LEGEND_STATES(:,30) = strpad('Mb in component Mb (micromolar)'); LEGEND_CONSTANTS(:,98) = strpad('k_c8 in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,31) = strpad('Mbp27 in component Mbp27 (micromolar)'); LEGEND_CONSTANTS(:,99) = strpad('V_1p27 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,100) = strpad('K_1p27 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,101) = strpad('V_2p27 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,102) = strpad('K_2p27 in component model_parameters (micromolar)'); LEGEND_STATES(:,32) = strpad('p27p in component p27p (micromolar)'); LEGEND_CONSTANTS(:,103) = strpad('k_ddp27 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,104) = strpad('V_dp27p in component model_parameters (flux)'); LEGEND_CONSTANTS(:,105) = strpad('K_dp27skp2 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,106) = strpad('K_dp27p in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,107) = strpad('k_ddp27p in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,33) = strpad('Cdh1i in component Cdh1i (micromolar)'); LEGEND_CONSTANTS(:,108) = strpad('V_2cdh1 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,109) = strpad('K_2cdh1 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,110) = strpad('V_1cdh1 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,111) = strpad('K_1cdh1 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,112) = strpad('k_dcdh1i in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,113) = strpad('v_scdh1a in component model_parameters (flux)'); LEGEND_CONSTANTS(:,114) = strpad('k_dcdh1a in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,34) = strpad('Pai in component Pai (micromolar)'); LEGEND_CONSTANTS(:,115) = strpad('v_spai in component model_parameters (flux)'); LEGEND_CONSTANTS(:,116) = strpad('V_6a in component model_parameters (flux)'); LEGEND_CONSTANTS(:,117) = strpad('x_a1 in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,118) = strpad('x_a2 in component model_parameters (per_micromolar)'); LEGEND_CONSTANTS(:,119) = strpad('K_6a in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,120) = strpad('V_m5a in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,121) = strpad('a_a in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,122) = strpad('K_5a in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,123) = strpad('k_dpai in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,124) = strpad('k_dpa in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,35) = strpad('Cb in component Cb (micromolar)'); LEGEND_CONSTANTS(:,125) = strpad('v_cb in component model_parameters (flux)'); LEGEND_CONSTANTS(:,126) = strpad('k_com4 in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,127) = strpad('Cdk1_tot in component model_parameters (micromolar)'); LEGEND_STATES(:,36) = strpad('Mbi in component Mbi (micromolar)'); LEGEND_CONSTANTS(:,128) = strpad('k_decom4 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,129) = strpad('V_db in component model_parameters (flux)'); LEGEND_CONSTANTS(:,130) = strpad('K_db in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,131) = strpad('K_dbcdc20 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,132) = strpad('K_dbcdh1 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,133) = strpad('k_ddb in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,134) = strpad('V_m2b in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,135) = strpad('i_b3 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,136) = strpad('K_2b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,137) = strpad('V_m1b in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,37) = strpad('Pb in component Pb (micromolar)'); LEGEND_CONSTANTS(:,138) = strpad('K_1b in component model_parameters (micromolar)'); LEGEND_STATES(:,38) = strpad('Cdc20i in component Cdc20i (micromolar)'); LEGEND_CONSTANTS(:,139) = strpad('v_scdc20i in component model_parameters (flux)'); LEGEND_CONSTANTS(:,140) = strpad('V_m3b in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,141) = strpad('K_3b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,142) = strpad('V_m4b in component model_parameters (flux)'); LEGEND_CONSTANTS(:,143) = strpad('K_4b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,144) = strpad('k_dcdc20i in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,145) = strpad('k_dcdc20a in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,39) = strpad('Pbi in component Pbi (micromolar)'); LEGEND_CONSTANTS(:,146) = strpad('v_spbi in component model_parameters (flux)'); LEGEND_CONSTANTS(:,147) = strpad('V_6b in component model_parameters (flux)'); LEGEND_CONSTANTS(:,148) = strpad('x_b1 in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,149) = strpad('x_b2 in component model_parameters (per_micromolar)'); LEGEND_CONSTANTS(:,150) = strpad('K_6b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,151) = strpad('V_m5b in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,152) = strpad('a_b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,153) = strpad('K_5b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,154) = strpad('k_dpbi in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,155) = strpad('k_dpb in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,156) = strpad('v_swee1 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,157) = strpad('k_sw in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,40) = strpad('Mw in component Mw (micromolar)'); LEGEND_CONSTANTS(:,158) = strpad('V_m7b in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,159) = strpad('i_b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,160) = strpad('K_7b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,161) = strpad('V_m8b in component model_parameters (flux)'); LEGEND_STATES(:,41) = strpad('Wee1p in component Wee1p (micromolar)'); LEGEND_CONSTANTS(:,162) = strpad('K_8b in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,163) = strpad('k_dwee1 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,164) = strpad('k_dwee1p in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,42) = strpad('Cdc45 in component Cdc45 (micromolar)'); LEGEND_CONSTANTS(:,165) = strpad('V_1cdc45 in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,166) = strpad('Cdc45_tot in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,167) = strpad('K_1cdc45 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,168) = strpad('V_2cdc45 in component model_parameters (flux)'); LEGEND_CONSTANTS(:,169) = strpad('K_2cdc45 in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,170) = strpad('k_spol in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,171) = strpad('Pol_tot in component model_parameters (micromolar)'); LEGEND_STATES(:,43) = strpad('Pol in component Pol (micromolar)'); LEGEND_CONSTANTS(:,172) = strpad('k_dpol in component model_parameters (first_order_rate_constant)'); LEGEND_STATES(:,44) = strpad('Primer in component Primer (micromolar)'); LEGEND_CONSTANTS(:,173) = strpad('k_sprim in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,174) = strpad('k_dprim in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,175) = strpad('k_aatr in component model_parameters (second_order_rate_constant)'); LEGEND_CONSTANTS(:,176) = strpad('ATR_tot in component model_parameters (micromolar)'); LEGEND_STATES(:,45) = strpad('ATR in component ATR (micromolar)'); LEGEND_CONSTANTS(:,177) = strpad('k_datr in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,178) = strpad('V_1chk in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,179) = strpad('Chk1_tot in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,180) = strpad('K_1chk in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,181) = strpad('V_2chk in component model_parameters (flux)'); LEGEND_CONSTANTS(:,182) = strpad('K_2chk in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,183) = strpad('v_sw in component model_parameters (flux)'); LEGEND_STATES(:,46) = strpad('BN in component BN (nanomolar)'); LEGEND_CONSTANTS(:,184) = strpad('n_gerard in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,185) = strpad('K_iw in component model_parameters (nanomolar)'); LEGEND_CONSTANTS(:,186) = strpad('v_dw in component model_parameters (flux)'); LEGEND_CONSTANTS(:,187) = strpad('K_dw in component model_parameters (micromolar)'); LEGEND_STATES(:,47) = strpad('X in component X (micromolar)'); LEGEND_CONSTANTS(:,188) = strpad('V_1x in component model_parameters (first_order_rate_constant)'); LEGEND_CONSTANTS(:,189) = strpad('X_tot in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,190) = strpad('K_1x in component model_parameters (micromolar)'); LEGEND_CONSTANTS(:,191) = strpad('V_2x in component model_parameters (flux)'); LEGEND_CONSTANTS(:,192) = strpad('K_2x in component model_parameters (micromolar)'); LEGEND_STATES(:,48) = strpad('CbA in component CbA (micromolar)'); LEGEND_STATES(:,49) = strpad('MP in component MP (nanomolar)'); LEGEND_CONSTANTS(:,193) = strpad('vsP in component MP (nano_flux)'); LEGEND_CONSTANTS(:,194) = strpad('vmP in component MP (nano_flux)'); LEGEND_CONSTANTS(:,195) = strpad('kdmp in component MP (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,196) = strpad('KAP in component MP (nanomolar)'); LEGEND_CONSTANTS(:,197) = strpad('KmP in component MP (nanomolar)'); LEGEND_CONSTANTS(:,198) = strpad('n in component model_parameters (dimensionless)'); LEGEND_STATES(:,50) = strpad('MC in component MC (nanomolar)'); LEGEND_CONSTANTS(:,199) = strpad('vsC in component MC (nano_flux)'); LEGEND_CONSTANTS(:,200) = strpad('vmC in component MC (nano_flux)'); LEGEND_CONSTANTS(:,201) = strpad('kdmc in component MC (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,202) = strpad('KAC in component MC (nanomolar)'); LEGEND_CONSTANTS(:,203) = strpad('KmC in component MC (nanomolar)'); LEGEND_STATES(:,51) = strpad('MB in component MB (nanomolar)'); LEGEND_CONSTANTS(:,204) = strpad('vsB in component MB (nano_flux)'); LEGEND_CONSTANTS(:,205) = strpad('vmB in component MB (nano_flux)'); LEGEND_CONSTANTS(:,206) = strpad('kdmb in component MB (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,207) = strpad('KIB in component MB (nanomolar)'); LEGEND_CONSTANTS(:,208) = strpad('KmB in component MB (nanomolar)'); LEGEND_CONSTANTS(:,209) = strpad('m in component model_parameters (dimensionless)'); LEGEND_STATES(:,52) = strpad('RN in component RN (nanomolar)'); LEGEND_STATES(:,53) = strpad('MR in component MR (nanomolar)'); LEGEND_CONSTANTS(:,210) = strpad('vsR in component MR (nano_flux)'); LEGEND_CONSTANTS(:,211) = strpad('vmR in component MR (nano_flux)'); LEGEND_CONSTANTS(:,212) = strpad('kdmr in component MR (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,213) = strpad('KAR in component MR (nanomolar)'); LEGEND_CONSTANTS(:,214) = strpad('KmR in component MR (nanomolar)'); LEGEND_CONSTANTS(:,215) = strpad('h in component model_parameters (dimensionless)'); LEGEND_STATES(:,54) = strpad('PC in component PC (nanomolar)'); LEGEND_CONSTANTS(:,216) = strpad('ksP in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,217) = strpad('Kp in component model_parameters (nanomolar)'); LEGEND_CONSTANTS(:,218) = strpad('Kdp in component model_parameters (nanomolar)'); LEGEND_CONSTANTS(:,219) = strpad('k3 in component model_parameters (second_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,220) = strpad('k4 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,221) = strpad('kdn in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,222) = strpad('V1P in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,223) = strpad('V2P in component model_parameters (nano_flux)'); LEGEND_STATES(:,55) = strpad('PCP in component PCP (nanomolar)'); LEGEND_STATES(:,56) = strpad('PCC in component PCC (nanomolar)'); LEGEND_STATES(:,57) = strpad('CC in component CC (nanomolar)'); LEGEND_CONSTANTS(:,224) = strpad('ksC in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,225) = strpad('kdnc in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,226) = strpad('V1C in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,227) = strpad('V2C in component model_parameters (nano_flux)'); LEGEND_STATES(:,58) = strpad('CCP in component CCP (nanomolar)'); LEGEND_STATES(:,59) = strpad('RC in component RC (nanomolar)'); LEGEND_CONSTANTS(:,228) = strpad('ksR in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,229) = strpad('Kd in component model_parameters (nanomolar)'); LEGEND_CONSTANTS(:,230) = strpad('k9 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,231) = strpad('k10 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,232) = strpad('vdRC in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,233) = strpad('vdPC in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,234) = strpad('vdCC in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,235) = strpad('k1 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,236) = strpad('k2 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,237) = strpad('V1PC in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,238) = strpad('V2PC in component model_parameters (nano_flux)'); LEGEND_STATES(:,60) = strpad('PCCP in component PCCP (nanomolar)'); LEGEND_STATES(:,61) = strpad('PCN in component PCN (nanomolar)'); LEGEND_CONSTANTS(:,239) = strpad('k7 in component model_parameters (second_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,240) = strpad('k8 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,241) = strpad('V3PC in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,242) = strpad('V4PC in component model_parameters (nano_flux)'); LEGEND_STATES(:,62) = strpad('PCNP in component PCNP (nanomolar)'); LEGEND_STATES(:,63) = strpad('IN in component IN (nanomolar)'); LEGEND_CONSTANTS(:,243) = strpad('vdRN in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,244) = strpad('vdPCC in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,245) = strpad('vdPCN in component model_parameters (nano_flux)'); LEGEND_STATES(:,64) = strpad('BC in component BC (nanomolar)'); LEGEND_CONSTANTS(:,246) = strpad('ksB in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,247) = strpad('k5 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,248) = strpad('k6 in component model_parameters (first_order_rate_constant_nano)'); LEGEND_CONSTANTS(:,249) = strpad('V1B in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,250) = strpad('V2B in component model_parameters (nano_flux)'); LEGEND_STATES(:,65) = strpad('BCP in component BCP (nanomolar)'); LEGEND_CONSTANTS(:,251) = strpad('vdBC in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,252) = strpad('V3B in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,253) = strpad('V4B in component model_parameters (nano_flux)'); LEGEND_STATES(:,66) = strpad('BNP in component BNP (nanomolar)'); LEGEND_CONSTANTS(:,254) = strpad('vdBN in component model_parameters (nano_flux)'); LEGEND_CONSTANTS(:,255) = strpad('vdIN in component model_parameters (nano_flux)'); LEGEND_RATES(:,1) = strpad('d/dt AP1 in component AP1 (micromolar)'); LEGEND_RATES(:,2) = strpad('d/dt pRB in component pRB (micromolar)'); LEGEND_RATES(:,4) = strpad('d/dt pRBc1 in component pRBc1 (micromolar)'); LEGEND_RATES(:,7) = strpad('d/dt pRBp in component pRBp (micromolar)'); LEGEND_RATES(:,10) = strpad('d/dt pRBc2 in component pRBc2 (micromolar)'); LEGEND_RATES(:,9) = strpad('d/dt pRBpp in component pRBpp (micromolar)'); LEGEND_RATES(:,3) = strpad('d/dt E2F in component E2F (micromolar)'); LEGEND_RATES(:,12) = strpad('d/dt E2Fp in component E2Fp (micromolar)'); LEGEND_RATES(:,13) = strpad('d/dt Cd in component Cd (micromolar)'); LEGEND_RATES(:,14) = strpad('d/dt Mdi in component Mdi (micromolar)'); LEGEND_RATES(:,5) = strpad('d/dt Md in component Md (micromolar)'); LEGEND_RATES(:,6) = strpad('d/dt Mdp27 in component Mdp27 (micromolar)'); LEGEND_RATES(:,16) = strpad('d/dt Ce in component Ce (micromolar)'); LEGEND_RATES(:,17) = strpad('d/dt Mei in component Mei (micromolar)'); LEGEND_RATES(:,8) = strpad('d/dt Me in component Me (micromolar)'); LEGEND_RATES(:,21) = strpad('d/dt Skp2 in component Skp2 (micromolar)'); LEGEND_RATES(:,18) = strpad('d/dt Mep27 in component Mep27 (micromolar)'); LEGEND_RATES(:,25) = strpad('d/dt Pei in component Pei (micromolar)'); LEGEND_RATES(:,23) = strpad('d/dt Pe in component Pe (micromolar)'); LEGEND_RATES(:,27) = strpad('d/dt Ca in component Ca (micromolar)'); LEGEND_RATES(:,19) = strpad('d/dt Mai in component Mai (micromolar)'); LEGEND_RATES(:,11) = strpad('d/dt Ma in component Ma (micromolar)'); LEGEND_RATES(:,20) = strpad('d/dt Map27 in component Map27 (micromolar)'); LEGEND_RATES(:,15) = strpad('d/dt p27 in component p27 (micromolar)'); LEGEND_RATES(:,32) = strpad('d/dt p27p in component p27p (micromolar)'); LEGEND_RATES(:,33) = strpad('d/dt Cdh1i in component Cdh1i (micromolar)'); LEGEND_RATES(:,24) = strpad('d/dt Cdh1a in component Cdh1a (micromolar)'); LEGEND_RATES(:,34) = strpad('d/dt Pai in component Pai (micromolar)'); LEGEND_RATES(:,29) = strpad('d/dt Pa in component Pa (micromolar)'); LEGEND_RATES(:,35) = strpad('d/dt Cb in component Cb (micromolar)'); LEGEND_RATES(:,36) = strpad('d/dt Mbi in component Mbi (micromolar)'); LEGEND_RATES(:,30) = strpad('d/dt Mb in component Mb (micromolar)'); LEGEND_RATES(:,31) = strpad('d/dt Mbp27 in component Mbp27 (micromolar)'); LEGEND_RATES(:,38) = strpad('d/dt Cdc20i in component Cdc20i (micromolar)'); LEGEND_RATES(:,28) = strpad('d/dt Cdc20a in component Cdc20a (micromolar)'); LEGEND_RATES(:,39) = strpad('d/dt Pbi in component Pbi (micromolar)'); LEGEND_RATES(:,37) = strpad('d/dt Pb in component Pb (micromolar)'); LEGEND_RATES(:,22) = strpad('d/dt Wee1 in component Wee1 (micromolar)'); LEGEND_RATES(:,41) = strpad('d/dt Wee1p in component Wee1p (micromolar)'); LEGEND_RATES(:,42) = strpad('d/dt Cdc45 in component Cdc45 (micromolar)'); LEGEND_RATES(:,43) = strpad('d/dt Pol in component Pol (micromolar)'); LEGEND_RATES(:,44) = strpad('d/dt Primer in component Primer (micromolar)'); LEGEND_RATES(:,45) = strpad('d/dt ATR in component ATR (micromolar)'); LEGEND_RATES(:,26) = strpad('d/dt Chk1 in component Chk1 (micromolar)'); LEGEND_RATES(:,40) = strpad('d/dt Mw in component Mw (micromolar)'); LEGEND_RATES(:,47) = strpad('d/dt X in component X (micromolar)'); LEGEND_RATES(:,48) = strpad('d/dt CbA in component CbA (micromolar)'); LEGEND_RATES(:,49) = strpad('d/dt MP in component MP (nanomolar)'); LEGEND_RATES(:,50) = strpad('d/dt MC in component MC (nanomolar)'); LEGEND_RATES(:,51) = strpad('d/dt MB in component MB (nanomolar)'); LEGEND_RATES(:,53) = strpad('d/dt MR in component MR (nanomolar)'); LEGEND_RATES(:,54) = strpad('d/dt PC in component PC (nanomolar)'); LEGEND_RATES(:,57) = strpad('d/dt CC in component CC (nanomolar)'); LEGEND_RATES(:,59) = strpad('d/dt RC in component RC (nanomolar)'); LEGEND_RATES(:,55) = strpad('d/dt PCP in component PCP (nanomolar)'); LEGEND_RATES(:,58) = strpad('d/dt CCP in component CCP (nanomolar)'); LEGEND_RATES(:,56) = strpad('d/dt PCC in component PCC (nanomolar)'); LEGEND_RATES(:,61) = strpad('d/dt PCN in component PCN (nanomolar)'); LEGEND_RATES(:,52) = strpad('d/dt RN in component RN (nanomolar)'); LEGEND_RATES(:,60) = strpad('d/dt PCCP in component PCCP (nanomolar)'); LEGEND_RATES(:,62) = strpad('d/dt PCNP in component PCNP (nanomolar)'); LEGEND_RATES(:,64) = strpad('d/dt BC in component BC (nanomolar)'); LEGEND_RATES(:,65) = strpad('d/dt BCP in component BCP (nanomolar)'); LEGEND_RATES(:,46) = strpad('d/dt BN in component BN (nanomolar)'); LEGEND_RATES(:,66) = strpad('d/dt BNP in component BNP (nanomolar)'); LEGEND_RATES(:,63) = strpad('d/dt IN in component IN (nanomolar)'); LEGEND_STATES = LEGEND_STATES'; LEGEND_ALGEBRAIC = LEGEND_ALGEBRAIC'; LEGEND_RATES = LEGEND_RATES'; LEGEND_CONSTANTS = LEGEND_CONSTANTS'; end function [STATES, CONSTANTS] = initConsts() VOI = 0; CONSTANTS = []; STATES = []; ALGEBRAIC = []; STATES(:,1) = 0.01; CONSTANTS(:,1) = 1; CONSTANTS(:,2) = 1; CONSTANTS(:,3) = 0.1; CONSTANTS(:,4) = 0.15; CONSTANTS(:,5) = 17; STATES(:,2) = 1; CONSTANTS(:,6) = 0.8; CONSTANTS(:,7) = 0.05; STATES(:,3) = 0.01; CONSTANTS(:,8) = 0.5; STATES(:,4) = 0.1; CONSTANTS(:,9) = 2.2; CONSTANTS(:,10) = 0.1; STATES(:,5) = 0.01; STATES(:,6) = 0.01; CONSTANTS(:,11) = 2; STATES(:,7) = 0.25; CONSTANTS(:,12) = 0.1; CONSTANTS(:,13) = 0.01; CONSTANTS(:,14) = 1; CONSTANTS(:,15) = 0.1; STATES(:,8) = 0.01; CONSTANTS(:,16) = 2; STATES(:,9) = 0.1; CONSTANTS(:,17) = 0.1; CONSTANTS(:,18) = 0.025; CONSTANTS(:,19) = 0.5; STATES(:,10) = 0.05; CONSTANTS(:,20) = 0.06; CONSTANTS(:,21) = 0.04; CONSTANTS(:,22) = 0.17; CONSTANTS(:,23) = 4; STATES(:,11) = 0.01; CONSTANTS(:,24) = 5; CONSTANTS(:,25) = 0.75; STATES(:,12) = 0.01; CONSTANTS(:,26) = 5; CONSTANTS(:,27) = 0.002; CONSTANTS(:,28) = 1.1; STATES(:,13) = 0.01; CONSTANTS(:,29) = 0.4; CONSTANTS(:,30) = 0.005; CONSTANTS(:,31) = 0.1; CONSTANTS(:,32) = 2; CONSTANTS(:,33) = 0.175; CONSTANTS(:,34) = 1.5; STATES(:,14) = 0.01; CONSTANTS(:,35) = 0.1; CONSTANTS(:,36) = 5; CONSTANTS(:,37) = 0.1; CONSTANTS(:,38) = 0.005; CONSTANTS(:,39) = 0.2; CONSTANTS(:,40) = 0.1; CONSTANTS(:,41) = 1; CONSTANTS(:,42) = 0.1; CONSTANTS(:,43) = 0.15; STATES(:,15) = 0.01; CONSTANTS(:,44) = 0.05; STATES(:,16) = 0.01; CONSTANTS(:,45) = 0.29; CONSTANTS(:,46) = 0.1; CONSTANTS(:,47) = 2; CONSTANTS(:,48) = 0.2; CONSTANTS(:,49) = 2; STATES(:,17) = 0.01; STATES(:,18) = 0.01; STATES(:,19) = 0.01; STATES(:,20) = 0.01; CONSTANTS(:,50) = 0.1; CONSTANTS(:,51) = 3; STATES(:,21) = 0.01; CONSTANTS(:,52) = 2; CONSTANTS(:,53) = 0.1; CONSTANTS(:,54) = 0.005; CONSTANTS(:,55) = 1.4; STATES(:,22) = 0.1; CONSTANTS(:,56) = 0.5; CONSTANTS(:,57) = 0.1; CONSTANTS(:,58) = 2; STATES(:,23) = 0.01; CONSTANTS(:,59) = 0.1; CONSTANTS(:,60) = 0.2; CONSTANTS(:,61) = 0.1; CONSTANTS(:,62) = 0.15; CONSTANTS(:,63) = 1.1; CONSTANTS(:,64) = 0.5; STATES(:,24) = 0.01; CONSTANTS(:,65) = 0.4; CONSTANTS(:,66) = 0.005; STATES(:,25) = 0.01; CONSTANTS(:,67) = 0.13; CONSTANTS(:,68) = 0.8; CONSTANTS(:,69) = 1; CONSTANTS(:,70) = 1; STATES(:,26) = 0.01; CONSTANTS(:,71) = 0.1; CONSTANTS(:,72) = 5; CONSTANTS(:,73) = 0.25; CONSTANTS(:,74) = 0.1; CONSTANTS(:,75) = 0.15; CONSTANTS(:,76) = 0.075; STATES(:,27) = 0.01; CONSTANTS(:,77) = 0.0375; CONSTANTS(:,78) = 0.1; CONSTANTS(:,79) = 2; CONSTANTS(:,80) = 0.2; CONSTANTS(:,81) = 0.1; CONSTANTS(:,82) = 2.5; CONSTANTS(:,83) = 1.1; STATES(:,28) = 0.01; CONSTANTS(:,84) = 2; CONSTANTS(:,85) = 0.005; CONSTANTS(:,86) = 1.85; CONSTANTS(:,87) = 0.5; CONSTANTS(:,88) = 0.1; CONSTANTS(:,89) = 2; STATES(:,29) = 0.01; CONSTANTS(:,90) = 0.1; CONSTANTS(:,91) = 0.15; CONSTANTS(:,92) = 0.125; CONSTANTS(:,93) = 0.8; CONSTANTS(:,94) = 0.1; CONSTANTS(:,95) = 0.1; CONSTANTS(:,96) = 2; CONSTANTS(:,97) = 0.12; STATES(:,30) = 0.01; CONSTANTS(:,98) = 0.2; STATES(:,31) = 0.01; CONSTANTS(:,99) = 100; CONSTANTS(:,100) = 0.5; CONSTANTS(:,101) = 0.1; CONSTANTS(:,102) = 0.5; STATES(:,32) = 0.01; CONSTANTS(:,103) = 0.06; CONSTANTS(:,104) = 5; CONSTANTS(:,105) = 0.1; CONSTANTS(:,106) = 0.1; CONSTANTS(:,107) = 0.01; STATES(:,33) = 0.01; CONSTANTS(:,108) = 8; CONSTANTS(:,109) = 0.01; CONSTANTS(:,110) = 1.25; CONSTANTS(:,111) = 0.01; CONSTANTS(:,112) = 0.2; CONSTANTS(:,113) = 0.11; CONSTANTS(:,114) = 0.1; STATES(:,34) = 0.01; CONSTANTS(:,115) = 0.105; CONSTANTS(:,116) = 1; CONSTANTS(:,117) = 1; CONSTANTS(:,118) = 1; CONSTANTS(:,119) = 0.1; CONSTANTS(:,120) = 4; CONSTANTS(:,121) = 0.2; CONSTANTS(:,122) = 0.1; CONSTANTS(:,123) = 0.15; CONSTANTS(:,124) = 0.075; STATES(:,35) = 0.01; CONSTANTS(:,125) = 0.05; CONSTANTS(:,126) = 0.25; CONSTANTS(:,127) = 0.5; STATES(:,36) = 0.01; CONSTANTS(:,128) = 0.1; CONSTANTS(:,129) = 0.06; CONSTANTS(:,130) = 0.005; CONSTANTS(:,131) = 0.2; CONSTANTS(:,132) = 0.1; CONSTANTS(:,133) = 0.005; CONSTANTS(:,134) = 2.1; CONSTANTS(:,135) = 0.5; CONSTANTS(:,136) = 0.1; CONSTANTS(:,137) = 3.9; STATES(:,37) = 0.01; CONSTANTS(:,138) = 0.1; STATES(:,38) = 0.01; CONSTANTS(:,139) = 0.1; CONSTANTS(:,140) = 8; CONSTANTS(:,141) = 0.1; CONSTANTS(:,142) = 0.7; CONSTANTS(:,143) = 0.1; CONSTANTS(:,144) = 0.14; CONSTANTS(:,145) = 0.05; STATES(:,39) = 0.01; CONSTANTS(:,146) = 0.12; CONSTANTS(:,147) = 1; CONSTANTS(:,148) = 1; CONSTANTS(:,149) = 1; CONSTANTS(:,150) = 0.1; CONSTANTS(:,151) = 5; CONSTANTS(:,152) = 0.11; CONSTANTS(:,153) = 0.1; CONSTANTS(:,154) = 0.2; CONSTANTS(:,155) = 0.1; CONSTANTS(:,156) = 0.06; CONSTANTS(:,157) = 5; STATES(:,40) = 0; CONSTANTS(:,158) = 1.2; CONSTANTS(:,159) = 0.75; CONSTANTS(:,160) = 0.1; CONSTANTS(:,161) = 1; STATES(:,41) = 0.01; CONSTANTS(:,162) = 0.1; CONSTANTS(:,163) = 0.1; CONSTANTS(:,164) = 0.2; STATES(:,42) = 0.01; CONSTANTS(:,165) = 0.8; CONSTANTS(:,166) = 0.5; CONSTANTS(:,167) = 0.02; CONSTANTS(:,168) = 0.12; CONSTANTS(:,169) = 0.02; CONSTANTS(:,170) = 0.8; CONSTANTS(:,171) = 0.5; STATES(:,43) = 0.01; CONSTANTS(:,172) = 0.2; STATES(:,44) = 0.01; CONSTANTS(:,173) = 0.05; CONSTANTS(:,174) = 0.15; CONSTANTS(:,175) = 0.022; CONSTANTS(:,176) = 0.5; STATES(:,45) = 0.01; CONSTANTS(:,177) = 0.15; CONSTANTS(:,178) = 4; CONSTANTS(:,179) = 0.5; CONSTANTS(:,180) = 0.5; CONSTANTS(:,181) = 0.1; CONSTANTS(:,182) = 0.5; CONSTANTS(:,183) = 0; STATES(:,46) = 0.1; CONSTANTS(:,184) = 4; CONSTANTS(:,185) = 0.5; CONSTANTS(:,186) = 0.12; CONSTANTS(:,187) = 0.5; STATES(:,47) = 0.01; CONSTANTS(:,188) = 10; CONSTANTS(:,189) = 1; CONSTANTS(:,190) = 0.1; CONSTANTS(:,191) = 2; CONSTANTS(:,192) = 0.1; STATES(:,48) = 0.01; STATES(:,49) = 0.1; CONSTANTS(:,193) = 2.4; CONSTANTS(:,194) = 2.2; CONSTANTS(:,195) = 0.02; CONSTANTS(:,196) = 0.6; CONSTANTS(:,197) = 0.3; CONSTANTS(:,198) = 2; STATES(:,50) = 1.2; CONSTANTS(:,199) = 2.2; CONSTANTS(:,200) = 2; CONSTANTS(:,201) = 0.02; CONSTANTS(:,202) = 0.6; CONSTANTS(:,203) = 0.4; STATES(:,51) = 9; CONSTANTS(:,204) = 1.8; CONSTANTS(:,205) = 1.3; CONSTANTS(:,206) = 0.02; CONSTANTS(:,207) = 2.2; CONSTANTS(:,208) = 0.4; CONSTANTS(:,209) = 2; STATES(:,52) = 0.1; STATES(:,53) = 1.5; CONSTANTS(:,210) = 1.6; CONSTANTS(:,211) = 1.6; CONSTANTS(:,212) = 0.02; CONSTANTS(:,213) = 0.6; CONSTANTS(:,214) = 0.4; CONSTANTS(:,215) = 2; STATES(:,54) = 0.1; CONSTANTS(:,216) = 1.2; CONSTANTS(:,217) = 1.006; CONSTANTS(:,218) = 0.1; CONSTANTS(:,219) = 0.8; CONSTANTS(:,220) = 0.4; CONSTANTS(:,221) = 0.02; CONSTANTS(:,222) = 9.6; CONSTANTS(:,223) = 0.6; STATES(:,55) = 0.1; STATES(:,56) = 0.1; STATES(:,57) = 0.1; CONSTANTS(:,224) = 3.2; CONSTANTS(:,225) = 0.02; CONSTANTS(:,226) = 1.2; CONSTANTS(:,227) = 0.2; STATES(:,58) = 0.1; STATES(:,59) = 0.1; CONSTANTS(:,228) = 1.7; CONSTANTS(:,229) = 0.3; CONSTANTS(:,230) = 0.8; CONSTANTS(:,231) = 0.4; CONSTANTS(:,232) = 4.4; CONSTANTS(:,233) = 3.4; CONSTANTS(:,234) = 1.4; CONSTANTS(:,235) = 0.8; CONSTANTS(:,236) = 0.4; CONSTANTS(:,237) = 2.4; CONSTANTS(:,238) = 0.2; STATES(:,60) = 0.1; STATES(:,61) = 0.1; CONSTANTS(:,239) = 1; CONSTANTS(:,240) = 0.2; CONSTANTS(:,241) = 2.4; CONSTANTS(:,242) = 0.2; STATES(:,62) = 0.1; STATES(:,63) = 0.1; CONSTANTS(:,243) = 0.8; CONSTANTS(:,244) = 1.4; CONSTANTS(:,245) = 1.4; STATES(:,64) = 0.1; CONSTANTS(:,246) = 0.32; CONSTANTS(:,247) = 0.8; CONSTANTS(:,248) = 0.4; CONSTANTS(:,249) = 1.4; CONSTANTS(:,250) = 0.2; STATES(:,65) = 0.1; CONSTANTS(:,251) = 3; CONSTANTS(:,252) = 1.4; CONSTANTS(:,253) = 0.4; STATES(:,66) = 0.1; CONSTANTS(:,254) = 3; CONSTANTS(:,255) = 1.6; 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(:,1).*CONSTANTS(:,2))./(CONSTANTS(:,3)+CONSTANTS(:,2)) - CONSTANTS(:,4).*STATES(:,1)).*CONSTANTS(:,5); RATES(:,2) = (((((CONSTANTS(:,6) - CONSTANTS(:,7).*STATES(:,2).*STATES(:,3))+ CONSTANTS(:,8).*STATES(:,4)) - (( CONSTANTS(:,9).*STATES(:,2))./(CONSTANTS(:,10)+STATES(:,2))).*(STATES(:,5)+STATES(:,6)))+( CONSTANTS(:,11).*STATES(:,7))./(CONSTANTS(:,12)+STATES(:,7))) - CONSTANTS(:,13).*STATES(:,2)).*CONSTANTS(:,5); RATES(:,4) = ( CONSTANTS(:,7).*STATES(:,2).*STATES(:,3) - CONSTANTS(:,8).*STATES(:,4)).*CONSTANTS(:,5); RATES(:,7) = (((((( (( CONSTANTS(:,9).*STATES(:,2))./(CONSTANTS(:,10)+STATES(:,2))).*(STATES(:,5)+STATES(:,6)) - ( CONSTANTS(:,11).*STATES(:,7))./(CONSTANTS(:,12)+STATES(:,7))) - (( CONSTANTS(:,14).*STATES(:,7))./(CONSTANTS(:,15)+STATES(:,7))).*STATES(:,8))+( CONSTANTS(:,16).*STATES(:,9))./(CONSTANTS(:,17)+STATES(:,9))) - CONSTANTS(:,18).*STATES(:,7).*STATES(:,3))+ CONSTANTS(:,19).*STATES(:,10)) - CONSTANTS(:,20).*STATES(:,7)).*CONSTANTS(:,5); RATES(:,10) = ( CONSTANTS(:,18).*STATES(:,7).*STATES(:,3) - CONSTANTS(:,19).*STATES(:,10)).*CONSTANTS(:,5); RATES(:,9) = (( (( CONSTANTS(:,14).*STATES(:,7))./(CONSTANTS(:,15)+STATES(:,7))).*STATES(:,8) - ( CONSTANTS(:,16).*STATES(:,9))./(CONSTANTS(:,17)+STATES(:,9))) - CONSTANTS(:,21).*STATES(:,9)).*CONSTANTS(:,5); RATES(:,3) = (((((((CONSTANTS(:,22) - CONSTANTS(:,7).*STATES(:,2).*STATES(:,3))+ CONSTANTS(:,8).*STATES(:,4)) - CONSTANTS(:,18).*STATES(:,7).*STATES(:,3))+ CONSTANTS(:,19).*STATES(:,10)) - ( CONSTANTS(:,23).*STATES(:,11).*STATES(:,3))./(CONSTANTS(:,24)+STATES(:,3)))+( CONSTANTS(:,25).*STATES(:,12))./(CONSTANTS(:,26)+STATES(:,12))) - CONSTANTS(:,27).*STATES(:,3)).*CONSTANTS(:,5); RATES(:,12) = ((( CONSTANTS(:,23).*STATES(:,11).*STATES(:,3))./(CONSTANTS(:,24)+STATES(:,3)) - ( CONSTANTS(:,25).*STATES(:,12))./(CONSTANTS(:,26)+STATES(:,12))) - CONSTANTS(:,28).*STATES(:,12)).*CONSTANTS(:,5); RATES(:,13) = ((((( CONSTANTS(:,29).*STATES(:,1)+( (( CONSTANTS(:,30).*STATES(:,3).*CONSTANTS(:,31))./(CONSTANTS(:,31)+STATES(:,2))).*CONSTANTS(:,32))./(CONSTANTS(:,32)+STATES(:,7))) - CONSTANTS(:,33).*STATES(:,13).*(CONSTANTS(:,34) - (STATES(:,14)+STATES(:,5)+STATES(:,6))))+ CONSTANTS(:,35).*STATES(:,14)) - ( CONSTANTS(:,36).*STATES(:,13))./(CONSTANTS(:,37)+STATES(:,13))) - CONSTANTS(:,38).*STATES(:,13)).*CONSTANTS(:,5); RATES(:,14) = ((( CONSTANTS(:,33).*STATES(:,13).*(CONSTANTS(:,34) - (STATES(:,14)+STATES(:,5)+STATES(:,6))) - CONSTANTS(:,35).*STATES(:,14))+( CONSTANTS(:,39).*STATES(:,5))./(CONSTANTS(:,40)+STATES(:,5))) - ( CONSTANTS(:,41).*STATES(:,14))./(CONSTANTS(:,42)+STATES(:,14))).*CONSTANTS(:,5); RATES(:,5) = (((( CONSTANTS(:,41).*STATES(:,14))./(CONSTANTS(:,42)+STATES(:,14)) - ( CONSTANTS(:,39).*STATES(:,5))./(CONSTANTS(:,40)+STATES(:,5))) - CONSTANTS(:,43).*STATES(:,5).*STATES(:,15))+ CONSTANTS(:,44).*STATES(:,6)).*CONSTANTS(:,5); RATES(:,6) = ( CONSTANTS(:,43).*STATES(:,5).*STATES(:,15) - CONSTANTS(:,44).*STATES(:,6)).*CONSTANTS(:,5); RATES(:,16) = ((((( (( CONSTANTS(:,45).*STATES(:,3).*CONSTANTS(:,46))./(CONSTANTS(:,46)+STATES(:,2))).*CONSTANTS(:,47))./(CONSTANTS(:,47)+STATES(:,7)) - CONSTANTS(:,48).*STATES(:,16).*(CONSTANTS(:,49) - (STATES(:,17)+STATES(:,8)+STATES(:,18)+STATES(:,19)+STATES(:,11)+STATES(:,20))))+ CONSTANTS(:,50).*STATES(:,17)) - ( (( CONSTANTS(:,51).*STATES(:,21))./(CONSTANTS(:,52)+STATES(:,21))).*STATES(:,16))./(CONSTANTS(:,53)+STATES(:,16))) - CONSTANTS(:,54).*STATES(:,16)).*CONSTANTS(:,5); RATES(:,17) = ((( CONSTANTS(:,48).*STATES(:,16).*(CONSTANTS(:,49) - (STATES(:,17)+STATES(:,8)+STATES(:,18)+STATES(:,19)+STATES(:,11)+STATES(:,20))) - CONSTANTS(:,50).*STATES(:,17))+( CONSTANTS(:,55).*(STATES(:,22)+CONSTANTS(:,56)).*STATES(:,8))./(CONSTANTS(:,57)+STATES(:,8))) - ( CONSTANTS(:,58).*STATES(:,23).*STATES(:,17))./(CONSTANTS(:,59)+STATES(:,17))).*CONSTANTS(:,5); RATES(:,8) = (((( CONSTANTS(:,58).*STATES(:,23).*STATES(:,17))./(CONSTANTS(:,59)+STATES(:,17)) - ( CONSTANTS(:,55).*(STATES(:,22)+CONSTANTS(:,56)).*STATES(:,8))./(CONSTANTS(:,57)+STATES(:,8))) - CONSTANTS(:,60).*STATES(:,8).*STATES(:,15))+ CONSTANTS(:,61).*STATES(:,18)).*CONSTANTS(:,5); RATES(:,21) = ((CONSTANTS(:,62) - ( (( CONSTANTS(:,63).*STATES(:,21))./(CONSTANTS(:,64)+STATES(:,21))).*STATES(:,24))./(CONSTANTS(:,65)+STATES(:,24))) - CONSTANTS(:,66).*STATES(:,21)).*CONSTANTS(:,5); RATES(:,18) = ( CONSTANTS(:,60).*STATES(:,8).*STATES(:,15) - CONSTANTS(:,61).*STATES(:,18)).*CONSTANTS(:,5); RATES(:,25) = (((CONSTANTS(:,67)+( CONSTANTS(:,68).*(CONSTANTS(:,69)+ CONSTANTS(:,70).*STATES(:,26)).*STATES(:,23))./(CONSTANTS(:,71)+STATES(:,23))) - ( CONSTANTS(:,72).*(STATES(:,8)+CONSTANTS(:,73)).*STATES(:,25))./(CONSTANTS(:,74)+STATES(:,25))) - CONSTANTS(:,75).*STATES(:,25)).*CONSTANTS(:,5); RATES(:,23) = (( CONSTANTS(:,72).*(STATES(:,8)+CONSTANTS(:,73)).*(STATES(:,25)./(CONSTANTS(:,74)+STATES(:,25))) - CONSTANTS(:,68).*(CONSTANTS(:,69)+ CONSTANTS(:,70).*STATES(:,26)).*(STATES(:,23)./(CONSTANTS(:,71)+STATES(:,23)))) - CONSTANTS(:,76).*STATES(:,23)).*CONSTANTS(:,5); RATES(:,27) = ((((( (( CONSTANTS(:,77).*STATES(:,3).*CONSTANTS(:,78))./(CONSTANTS(:,78)+STATES(:,2))).*CONSTANTS(:,79))./(CONSTANTS(:,79)+STATES(:,7)) - CONSTANTS(:,80).*STATES(:,27).*(CONSTANTS(:,49) - (STATES(:,17)+STATES(:,8)+STATES(:,18)+STATES(:,19)+STATES(:,11)+STATES(:,20))))+ CONSTANTS(:,81).*STATES(:,19)) - ( (( CONSTANTS(:,82).*STATES(:,27))./(CONSTANTS(:,83)+STATES(:,27))).*STATES(:,28))./(CONSTANTS(:,84)+STATES(:,28))) - CONSTANTS(:,85).*STATES(:,27)).*CONSTANTS(:,5); RATES(:,19) = ((( CONSTANTS(:,80).*STATES(:,27).*(CONSTANTS(:,49) - (STATES(:,17)+STATES(:,8)+STATES(:,18)+STATES(:,19)+STATES(:,11)+STATES(:,20))) - CONSTANTS(:,81).*STATES(:,19))+( CONSTANTS(:,86).*(STATES(:,22)+CONSTANTS(:,87)).*STATES(:,11))./(CONSTANTS(:,88)+STATES(:,11))) - ( CONSTANTS(:,89).*STATES(:,29).*STATES(:,19))./(CONSTANTS(:,90)+STATES(:,19))).*CONSTANTS(:,5); RATES(:,11) = (((( CONSTANTS(:,89).*STATES(:,29).*STATES(:,19))./(CONSTANTS(:,90)+STATES(:,19)) - ( CONSTANTS(:,86).*(STATES(:,22)+CONSTANTS(:,87)).*STATES(:,11))./(CONSTANTS(:,88)+STATES(:,11))) - CONSTANTS(:,91).*STATES(:,11).*STATES(:,15))+ CONSTANTS(:,92).*STATES(:,20)).*CONSTANTS(:,5); RATES(:,20) = ( CONSTANTS(:,91).*STATES(:,11).*STATES(:,15) - CONSTANTS(:,92).*STATES(:,20)).*CONSTANTS(:,5); RATES(:,15) = ((((((((((((CONSTANTS(:,93)+( (( CONSTANTS(:,94).*STATES(:,3).*CONSTANTS(:,95))./(CONSTANTS(:,95)+STATES(:,2))).*CONSTANTS(:,96))./(CONSTANTS(:,96)+STATES(:,7))) - CONSTANTS(:,43).*STATES(:,5).*STATES(:,15))+ CONSTANTS(:,44).*STATES(:,6)) - CONSTANTS(:,60).*STATES(:,8).*STATES(:,15))+ CONSTANTS(:,61).*STATES(:,18)) - CONSTANTS(:,91).*STATES(:,11).*STATES(:,15))+ CONSTANTS(:,92).*STATES(:,20)) - CONSTANTS(:,97).*STATES(:,30).*STATES(:,15))+ CONSTANTS(:,98).*STATES(:,31)) - ( CONSTANTS(:,99).*STATES(:,8).*STATES(:,15))./(CONSTANTS(:,100)+STATES(:,15)))+( CONSTANTS(:,101).*STATES(:,32))./(CONSTANTS(:,102)+STATES(:,32))) - CONSTANTS(:,103).*STATES(:,15)).*CONSTANTS(:,5); RATES(:,32) = (((( CONSTANTS(:,99).*STATES(:,8).*STATES(:,15))./(CONSTANTS(:,100)+STATES(:,15)) - ( CONSTANTS(:,101).*STATES(:,32))./(CONSTANTS(:,102)+STATES(:,32))) - ( (( CONSTANTS(:,104).*STATES(:,21))./(CONSTANTS(:,105)+STATES(:,21))).*STATES(:,32))./(CONSTANTS(:,106)+STATES(:,32))) - CONSTANTS(:,107).*STATES(:,32)).*CONSTANTS(:,5); RATES(:,33) = (( (( CONSTANTS(:,108).*STATES(:,24))./(CONSTANTS(:,109)+STATES(:,24))).*(STATES(:,11)+STATES(:,30)) - ( CONSTANTS(:,110).*STATES(:,33))./(CONSTANTS(:,111)+STATES(:,33))) - CONSTANTS(:,112).*STATES(:,33)).*CONSTANTS(:,5); RATES(:,24) = (((CONSTANTS(:,113)+( CONSTANTS(:,110).*STATES(:,33))./(CONSTANTS(:,111)+STATES(:,33))) - (( CONSTANTS(:,108).*STATES(:,24))./(CONSTANTS(:,109)+STATES(:,24))).*(STATES(:,11)+STATES(:,30))) - CONSTANTS(:,114).*STATES(:,24)).*CONSTANTS(:,5); RATES(:,34) = (((CONSTANTS(:,115)+( CONSTANTS(:,116).*(CONSTANTS(:,117)+ CONSTANTS(:,118).*STATES(:,26)).*STATES(:,29))./(CONSTANTS(:,119)+STATES(:,29))) - ( CONSTANTS(:,120).*(STATES(:,11)+CONSTANTS(:,121)).*STATES(:,34))./(CONSTANTS(:,122)+STATES(:,34))) - CONSTANTS(:,123).*STATES(:,34)).*CONSTANTS(:,5); RATES(:,29) = ((( CONSTANTS(:,120).*(STATES(:,11)+CONSTANTS(:,121)).*STATES(:,34))./(CONSTANTS(:,122)+STATES(:,34)) - ( CONSTANTS(:,116).*(CONSTANTS(:,117)+ CONSTANTS(:,118).*STATES(:,26)).*STATES(:,29))./(CONSTANTS(:,119)+STATES(:,29))) - CONSTANTS(:,124).*STATES(:,29)).*CONSTANTS(:,5); RATES(:,35) = ((((CONSTANTS(:,125) - CONSTANTS(:,126).*STATES(:,35).*(CONSTANTS(:,127) - (STATES(:,36)+STATES(:,30)+STATES(:,31))))+ CONSTANTS(:,128).*STATES(:,36)) - (( CONSTANTS(:,129).*STATES(:,35))./(CONSTANTS(:,130)+STATES(:,35))).*(STATES(:,28)./(CONSTANTS(:,131)+STATES(:,28))+STATES(:,24)./(CONSTANTS(:,132)+STATES(:,24)))) - CONSTANTS(:,133).*STATES(:,35)).*CONSTANTS(:,5); RATES(:,36) = ((( CONSTANTS(:,126).*STATES(:,35).*(CONSTANTS(:,127) - (STATES(:,36)+STATES(:,30)+STATES(:,31))) - CONSTANTS(:,128).*STATES(:,36))+( CONSTANTS(:,134).*(STATES(:,22)+CONSTANTS(:,135)).*STATES(:,30))./(CONSTANTS(:,136)+STATES(:,30))) - ( CONSTANTS(:,137).*STATES(:,37).*STATES(:,36))./(CONSTANTS(:,138)+STATES(:,36))).*CONSTANTS(:,5); RATES(:,30) = (((( CONSTANTS(:,137).*STATES(:,37).*STATES(:,36))./(CONSTANTS(:,138)+STATES(:,36)) - ( CONSTANTS(:,134).*(STATES(:,22)+CONSTANTS(:,135)).*STATES(:,30))./(CONSTANTS(:,136)+STATES(:,30))) - CONSTANTS(:,97).*STATES(:,30).*STATES(:,15))+ CONSTANTS(:,98).*STATES(:,31)).*CONSTANTS(:,5); RATES(:,31) = ( CONSTANTS(:,97).*STATES(:,30).*STATES(:,15) - CONSTANTS(:,98).*STATES(:,31)).*CONSTANTS(:,5); RATES(:,38) = (((CONSTANTS(:,139) - ( CONSTANTS(:,140).*STATES(:,30).*STATES(:,38))./(CONSTANTS(:,141)+STATES(:,38)))+( CONSTANTS(:,142).*STATES(:,28))./(CONSTANTS(:,143)+STATES(:,28))) - CONSTANTS(:,144).*STATES(:,38)).*CONSTANTS(:,5); RATES(:,28) = ((( CONSTANTS(:,140).*STATES(:,30).*STATES(:,38))./(CONSTANTS(:,141)+STATES(:,38)) - ( CONSTANTS(:,142).*STATES(:,28))./(CONSTANTS(:,143)+STATES(:,28))) - CONSTANTS(:,145).*STATES(:,28)).*CONSTANTS(:,5); RATES(:,39) = (((CONSTANTS(:,146)+( CONSTANTS(:,147).*(CONSTANTS(:,148)+ CONSTANTS(:,149).*STATES(:,26)).*STATES(:,37))./(CONSTANTS(:,150)+STATES(:,37))) - ( CONSTANTS(:,151).*(STATES(:,30)+CONSTANTS(:,152)).*STATES(:,39))./(CONSTANTS(:,153)+STATES(:,39))) - CONSTANTS(:,154).*STATES(:,39)).*CONSTANTS(:,5); RATES(:,37) = ((( CONSTANTS(:,151).*(STATES(:,30)+CONSTANTS(:,152)).*STATES(:,39))./(CONSTANTS(:,153)+STATES(:,39)) - ( CONSTANTS(:,147).*(CONSTANTS(:,148)+ CONSTANTS(:,149).*STATES(:,26)).*STATES(:,37))./(CONSTANTS(:,150)+STATES(:,37))) - CONSTANTS(:,155).*STATES(:,37)).*CONSTANTS(:,5); RATES(:,22) = ((((CONSTANTS(:,156)+ CONSTANTS(:,157).*STATES(:,40)) - ( CONSTANTS(:,158).*(STATES(:,30)+CONSTANTS(:,159)).*STATES(:,22))./(CONSTANTS(:,160)+STATES(:,22)))+( CONSTANTS(:,161).*STATES(:,41))./(CONSTANTS(:,162)+STATES(:,41))) - CONSTANTS(:,163).*STATES(:,22)).*CONSTANTS(:,5); RATES(:,41) = ((( CONSTANTS(:,158).*(STATES(:,30)+CONSTANTS(:,159)).*STATES(:,22))./(CONSTANTS(:,160)+STATES(:,22)) - ( CONSTANTS(:,161).*STATES(:,41))./(CONSTANTS(:,162)+STATES(:,41))) - CONSTANTS(:,164).*STATES(:,41)).*CONSTANTS(:,5); RATES(:,42) = (((( CONSTANTS(:,165).*STATES(:,8).*(CONSTANTS(:,166) - STATES(:,42)))./(CONSTANTS(:,167)+(CONSTANTS(:,166) - STATES(:,42))) - ( CONSTANTS(:,168).*STATES(:,42))./(CONSTANTS(:,169)+STATES(:,42))) - CONSTANTS(:,170).*(CONSTANTS(:,171) - STATES(:,43)).*STATES(:,42))+ CONSTANTS(:,172).*STATES(:,43)).*CONSTANTS(:,5); RATES(:,43) = ( CONSTANTS(:,170).*(CONSTANTS(:,171) - STATES(:,43)).*STATES(:,42) - CONSTANTS(:,172).*STATES(:,43)).*CONSTANTS(:,5); RATES(:,44) = ((( CONSTANTS(:,173).*STATES(:,43) - CONSTANTS(:,174).*STATES(:,44)) - CONSTANTS(:,175).*(CONSTANTS(:,176) - STATES(:,45)).*STATES(:,44))+ CONSTANTS(:,177).*STATES(:,45)).*CONSTANTS(:,5); RATES(:,45) = ( CONSTANTS(:,175).*(CONSTANTS(:,176) - STATES(:,45)).*STATES(:,44) - CONSTANTS(:,177).*STATES(:,45)).*CONSTANTS(:,5); RATES(:,26) = (( CONSTANTS(:,178).*STATES(:,45).*(CONSTANTS(:,179) - STATES(:,26)))./(CONSTANTS(:,180)+(CONSTANTS(:,179) - STATES(:,26))) - ( CONSTANTS(:,181).*STATES(:,26))./(CONSTANTS(:,182)+STATES(:,26))).*CONSTANTS(:,5); RATES(:,40) = ( CONSTANTS(:,183).*power(STATES(:,46), CONSTANTS(:,184)))./(power(CONSTANTS(:,185), CONSTANTS(:,184))+power(STATES(:,46), CONSTANTS(:,184))) - ( CONSTANTS(:,186).*STATES(:,40))./(CONSTANTS(:,187)+STATES(:,40)); RATES(:,47) = ( CONSTANTS(:,188).*STATES(:,11).*((CONSTANTS(:,189) - STATES(:,47))./(CONSTANTS(:,190)+(CONSTANTS(:,189) - STATES(:,47)))) - CONSTANTS(:,191).*(STATES(:,47)./(CONSTANTS(:,192)+STATES(:,47)))).*CONSTANTS(:,5); RATES(:,48) = (((( CONSTANTS(:,125).*STATES(:,47).*1.00000 - CONSTANTS(:,126).*STATES(:,48).*(CONSTANTS(:,127) - (STATES(:,36)+STATES(:,30)+STATES(:,31))))+ CONSTANTS(:,128).*STATES(:,36)) - CONSTANTS(:,129).*(STATES(:,48)./(CONSTANTS(:,130)+STATES(:,48))).*(STATES(:,28)./(CONSTANTS(:,131)+STATES(:,28))+STATES(:,24)./(CONSTANTS(:,132)+STATES(:,24)))) - CONSTANTS(:,133).*STATES(:,48)).*CONSTANTS(:,5); RATES(:,49) = ( CONSTANTS(:,193).*power(STATES(:,46), CONSTANTS(:,198)))./(power(CONSTANTS(:,196), CONSTANTS(:,198))+power(STATES(:,46), CONSTANTS(:,198))) - (( CONSTANTS(:,194).*STATES(:,49))./(CONSTANTS(:,197)+STATES(:,49))+ CONSTANTS(:,195).*STATES(:,49)); RATES(:,50) = ( CONSTANTS(:,199).*power(STATES(:,46), CONSTANTS(:,198)))./(power(CONSTANTS(:,202), CONSTANTS(:,198))+power(STATES(:,46), CONSTANTS(:,198))) - (( CONSTANTS(:,200).*STATES(:,50))./(CONSTANTS(:,203)+STATES(:,50))+ CONSTANTS(:,201).*STATES(:,50)); RATES(:,51) = ( CONSTANTS(:,204).*power(CONSTANTS(:,207), CONSTANTS(:,209)))./(power(CONSTANTS(:,207), CONSTANTS(:,209))+power(STATES(:,52), CONSTANTS(:,209))) - (( CONSTANTS(:,205).*STATES(:,51))./(CONSTANTS(:,208)+STATES(:,51))+ CONSTANTS(:,206).*STATES(:,51)); RATES(:,53) = ( CONSTANTS(:,210).*power(STATES(:,46), CONSTANTS(:,215)))./(power(CONSTANTS(:,213), CONSTANTS(:,215))+power(STATES(:,46), CONSTANTS(:,215))) - (( CONSTANTS(:,211).*STATES(:,53))./(CONSTANTS(:,214)+STATES(:,53))+ CONSTANTS(:,212).*STATES(:,53)); RATES(:,54) = ( CONSTANTS(:,216).*STATES(:,49)+( CONSTANTS(:,223).*STATES(:,55))./(CONSTANTS(:,218)+STATES(:,55))+ CONSTANTS(:,220).*STATES(:,56)) - (( CONSTANTS(:,222).*STATES(:,54))./(CONSTANTS(:,217)+STATES(:,54))+ CONSTANTS(:,219).*STATES(:,54).*STATES(:,57)+ CONSTANTS(:,221).*STATES(:,54)); RATES(:,57) = ( CONSTANTS(:,224).*STATES(:,50)+( CONSTANTS(:,227).*STATES(:,58))./(CONSTANTS(:,218)+STATES(:,58))+ CONSTANTS(:,220).*STATES(:,56)) - (( CONSTANTS(:,226).*STATES(:,57))./(CONSTANTS(:,217)+STATES(:,57))+ CONSTANTS(:,219).*STATES(:,54).*STATES(:,57)+ CONSTANTS(:,225).*STATES(:,57)); RATES(:,59) = ( CONSTANTS(:,228).*STATES(:,53)+ CONSTANTS(:,231).*STATES(:,52)) - ( CONSTANTS(:,230).*STATES(:,59)+( CONSTANTS(:,232).*STATES(:,59))./(CONSTANTS(:,229)+STATES(:,59))+ CONSTANTS(:,221).*STATES(:,59)); RATES(:,55) = ( CONSTANTS(:,222).*STATES(:,54))./(CONSTANTS(:,217)+STATES(:,54)) - (( CONSTANTS(:,223).*STATES(:,55))./(CONSTANTS(:,218)+STATES(:,55))+( CONSTANTS(:,233).*STATES(:,55))./(CONSTANTS(:,229)+STATES(:,55))+ CONSTANTS(:,221).*STATES(:,55)); RATES(:,58) = ( CONSTANTS(:,226).*STATES(:,57))./(CONSTANTS(:,217)+STATES(:,57)) - (( CONSTANTS(:,227).*STATES(:,58))./(CONSTANTS(:,218)+STATES(:,58))+( CONSTANTS(:,234).*STATES(:,58))./(CONSTANTS(:,229)+STATES(:,58))+ CONSTANTS(:,221).*STATES(:,58)); RATES(:,56) = (( CONSTANTS(:,238).*STATES(:,60))./(CONSTANTS(:,218)+STATES(:,60))+ CONSTANTS(:,219).*STATES(:,54).*STATES(:,57)+ CONSTANTS(:,236).*STATES(:,61)) - (( CONSTANTS(:,237).*STATES(:,56))./(CONSTANTS(:,217)+STATES(:,56))+ CONSTANTS(:,220).*STATES(:,56)+ CONSTANTS(:,235).*STATES(:,56)+ CONSTANTS(:,221).*STATES(:,56)); RATES(:,61) = (( CONSTANTS(:,242).*STATES(:,62))./(CONSTANTS(:,218)+STATES(:,62))+ CONSTANTS(:,235).*STATES(:,56)+ CONSTANTS(:,240).*STATES(:,63)) - (( CONSTANTS(:,241).*STATES(:,61))./(CONSTANTS(:,217)+STATES(:,61))+ CONSTANTS(:,236).*STATES(:,61)+ CONSTANTS(:,239).*STATES(:,46).*STATES(:,61)+ CONSTANTS(:,221).*STATES(:,61)); RATES(:,52) = CONSTANTS(:,230).*STATES(:,59) - ( CONSTANTS(:,231).*STATES(:,52)+( CONSTANTS(:,243).*STATES(:,52))./(CONSTANTS(:,229)+STATES(:,52))+ CONSTANTS(:,221).*STATES(:,52)); RATES(:,60) = ( CONSTANTS(:,237).*STATES(:,56))./(CONSTANTS(:,217)+STATES(:,56)) - (( CONSTANTS(:,238).*STATES(:,60))./(CONSTANTS(:,218)+STATES(:,60))+( CONSTANTS(:,244).*STATES(:,60))./(CONSTANTS(:,229)+STATES(:,60))+ CONSTANTS(:,221).*STATES(:,60)); RATES(:,62) = ( CONSTANTS(:,241).*STATES(:,61))./(CONSTANTS(:,217)+STATES(:,61)) - (( CONSTANTS(:,242).*STATES(:,62))./(CONSTANTS(:,218)+STATES(:,62))+( CONSTANTS(:,245).*STATES(:,62))./(CONSTANTS(:,229)+STATES(:,62))+ CONSTANTS(:,221).*STATES(:,62)); RATES(:,64) = (( CONSTANTS(:,250).*STATES(:,65))./(CONSTANTS(:,218)+STATES(:,65))+ CONSTANTS(:,248).*STATES(:,46)+ CONSTANTS(:,246).*STATES(:,51)) - (( CONSTANTS(:,249).*STATES(:,64))./(CONSTANTS(:,217)+STATES(:,64))+ CONSTANTS(:,247).*STATES(:,64)+ CONSTANTS(:,221).*STATES(:,64)); RATES(:,65) = ( CONSTANTS(:,249).*STATES(:,64))./(CONSTANTS(:,217)+STATES(:,64)) - (( CONSTANTS(:,250).*STATES(:,65))./(CONSTANTS(:,218)+STATES(:,65))+( CONSTANTS(:,251).*STATES(:,65))./(CONSTANTS(:,229)+STATES(:,65))+ CONSTANTS(:,221).*STATES(:,65)); RATES(:,46) = (( CONSTANTS(:,253).*STATES(:,66))./(CONSTANTS(:,218)+STATES(:,66))+ CONSTANTS(:,247).*STATES(:,64)+ CONSTANTS(:,240).*STATES(:,63)) - (( CONSTANTS(:,252).*STATES(:,46))./(CONSTANTS(:,217)+STATES(:,46))+ CONSTANTS(:,248).*STATES(:,46)+ CONSTANTS(:,239).*STATES(:,46).*STATES(:,61)+ CONSTANTS(:,221).*STATES(:,46)); RATES(:,66) = ( CONSTANTS(:,252).*STATES(:,46))./(CONSTANTS(:,217)+STATES(:,46)) - (( CONSTANTS(:,253).*STATES(:,66))./(CONSTANTS(:,218)+STATES(:,66))+( CONSTANTS(:,254).*STATES(:,66))./(CONSTANTS(:,229)+STATES(:,66))+ CONSTANTS(:,221).*STATES(:,66)); RATES(:,63) = CONSTANTS(:,239).*STATES(:,46).*STATES(:,61) - ( CONSTANTS(:,240).*STATES(:,63)+( CONSTANTS(:,255).*STATES(:,63))./(CONSTANTS(:,229)+STATES(:,63))+ CONSTANTS(:,221).*STATES(:,63)); 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