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 =65; end % There are a total of 34 entries in each of the rate and state variable arrays. % There are a total of 155 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_CONSTANTS(:,1) = strpad('t_ss in component Clamp_parameters (ms)'); LEGEND_CONSTANTS(:,2) = strpad('V_actHolding in component Clamp_parameters (mV)'); LEGEND_CONSTANTS(:,3) = strpad('V_actTest in component Clamp_parameters (mV)'); LEGEND_CONSTANTS(:,4) = strpad('Cai in component Clamp_parameters (mM)'); LEGEND_CONSTANTS(:,5) = strpad('Ki in component Clamp_parameters (mM)'); LEGEND_CONSTANTS(:,6) = strpad('Nai in component Clamp_parameters (mM)'); LEGEND_VOI = strpad('time in component time (ms)'); LEGEND_ALGEBRAIC(:,1) = strpad('V in component Patch_clamp_protocol (mV)'); LEGEND_ALGEBRAIC(:,51) = strpad('I_ion in component Clamped_current_Xi (pA)'); LEGEND_ALGEBRAIC(:,23) = strpad('I in component Ionic_currents (pA)'); LEGEND_ALGEBRAIC(:,27) = strpad('I in component Ionic_currents (pA)'); LEGEND_ALGEBRAIC(:,31) = strpad('I in component Ionic_currents (pA)'); LEGEND_ALGEBRAIC(:,35) = strpad('I in component Ionic_currents (pA)'); LEGEND_ALGEBRAIC(:,39) = strpad('I in component Ionic_currents (pA)'); LEGEND_ALGEBRAIC(:,15) = strpad('I_NCX in component INCX (pA)'); LEGEND_ALGEBRAIC(:,19) = strpad('I_NaK in component INaK (pA)'); LEGEND_ALGEBRAIC(:,43) = strpad('I in component Ionic_currents (pA)'); LEGEND_ALGEBRAIC(:,47) = strpad('I in component Ionic_currents (pA)'); LEGEND_STATES(:,1) = strpad('O_CaL in component ICaL_channel_states (dimensionless)'); LEGEND_ALGEBRAIC(:,2) = strpad('P_CaT in component ICaT (dimensionless)'); LEGEND_STATES(:,2) = strpad('O4 in component IBK_channel_states (dimensionless)'); LEGEND_ALGEBRAIC(:,10) = strpad('P_Kv in component IKv (dimensionless)'); LEGEND_STATES(:,3) = strpad('O_Na in component INa_channel_states (dimensionless)'); LEGEND_STATES(:,4) = strpad('X in component gating_kinetics (dimensionless)'); LEGEND_STATES(:,5) = strpad('X in component gating_kinetics (dimensionless)'); LEGEND_STATES(:,6) = strpad('X in component gating_kinetics (dimensionless)'); LEGEND_STATES(:,7) = strpad('X in component gating_kinetics (dimensionless)'); LEGEND_CONSTANTS(:,7) = strpad('g_CaL in component model_parameters (nS)'); LEGEND_CONSTANTS(:,8) = strpad('R in component model_parameters (J_per_K_mol)'); LEGEND_CONSTANTS(:,9) = strpad('T in component model_parameters (kelvin)'); LEGEND_CONSTANTS(:,10) = strpad('F in component model_parameters (C_per_mmol)'); LEGEND_CONSTANTS(:,11) = strpad('T0CaL in component model_parameters (kelvin)'); LEGEND_CONSTANTS(:,12) = strpad('T0CaT in component model_parameters (kelvin)'); LEGEND_CONSTANTS(:,13) = strpad('Q10Ca in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,14) = strpad('Cao in component model_parameters (mM)'); LEGEND_CONSTANTS(:,15) = strpad('C0_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,16) = strpad('C1_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,17) = strpad('C2_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,18) = strpad('C3_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,19) = strpad('C0Ca_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,20) = strpad('C1Ca_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,21) = strpad('C2Ca_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,22) = strpad('C3Ca_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,23) = strpad('ICa_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,24) = strpad('IVs_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,25) = strpad('IVf_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,26) = strpad('IVfCa_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,27) = strpad('IVsCa_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,28) = strpad('O_CaL_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,29) = strpad('g_CaT in component model_parameters (nS)'); LEGEND_CONSTANTS(:,30) = strpad('d_CaT_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,31) = strpad('f_CaT_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,32) = strpad('tau_dCaT in component model_parameters (ms)'); LEGEND_CONSTANTS(:,33) = strpad('g_Kv in component model_parameters (nS)'); LEGEND_CONSTANTS(:,34) = strpad('T0K in component model_parameters (kelvin)'); LEGEND_CONSTANTS(:,35) = strpad('Q10K in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,36) = strpad('Ko in component model_parameters (mM)'); LEGEND_CONSTANTS(:,37) = strpad('x_Kv_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,38) = strpad('y_Kv_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,39) = strpad('tau_xKv in component model_parameters (ms)'); LEGEND_CONSTANTS(:,40) = strpad('tau_yKv in component model_parameters (ms)'); LEGEND_CONSTANTS(:,41) = strpad('g_BK in component model_parameters (nS)'); LEGEND_CONSTANTS(:,42) = strpad('BK_C0_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,43) = strpad('BK_C1_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,44) = strpad('BK_C2_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,45) = strpad('BK_C3_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,46) = strpad('BK_C4_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,47) = strpad('O0_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,48) = strpad('O1_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,49) = strpad('O2_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,50) = strpad('O3_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,51) = strpad('O4_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,52) = strpad('k_on in component model_parameters (per_ms)'); LEGEND_CONSTANTS(:,53) = strpad('k_off_C in component model_parameters (per_ms)'); LEGEND_CONSTANTS(:,54) = strpad('k_off_O in component model_parameters (per_ms)'); LEGEND_CONSTANTS(:,55) = strpad('g_Na in component model_parameters (nS)'); LEGEND_CONSTANTS(:,56) = strpad('T0Na in component model_parameters (kelvin)'); LEGEND_CONSTANTS(:,57) = strpad('Q10Na in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,58) = strpad('Nao in component model_parameters (mM)'); LEGEND_CONSTANTS(:,59) = strpad('Na_C1_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,60) = strpad('Na_C2_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,61) = strpad('Na_C3_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,62) = strpad('I1_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,63) = strpad('I2_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,64) = strpad('O_Na_init in component initial_conditions (dimensionless)'); LEGEND_CONSTANTS(:,65) = strpad('K_mCa in component model_parameters (mM)'); LEGEND_CONSTANTS(:,66) = strpad('K_mNai in component model_parameters (mM)'); LEGEND_CONSTANTS(:,67) = strpad('k_sat in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,68) = strpad('gamma in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,69) = strpad('P_NCX in component model_parameters (pA)'); LEGEND_CONSTANTS(:,70) = strpad('K_mK in component model_parameters (mM)'); LEGEND_CONSTANTS(:,71) = strpad('K_mNa in component model_parameters (mM)'); LEGEND_CONSTANTS(:,72) = strpad('P_NaK in component model_parameters (pA)'); LEGEND_CONSTANTS(:,73) = strpad('g_NsNa in component model_parameters (nS)'); LEGEND_CONSTANTS(:,74) = strpad('g_NsK in component model_parameters (nS)'); LEGEND_CONSTANTS(:,75) = strpad('Cm in component model_parameters (pF)'); LEGEND_CONSTANTS(:,76) = strpad('V_cell in component model_parameters (litre)'); LEGEND_CONSTANTS(:,77) = strpad('n_CRT in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,78) = strpad('n_CaM in component model_parameters (dimensionless)'); LEGEND_CONSTANTS(:,79) = strpad('CRT_total in component model_parameters (mM)'); LEGEND_CONSTANTS(:,80) = strpad('CaM_total in component model_parameters (mM)'); LEGEND_CONSTANTS(:,81) = strpad('K_D_CRT in component model_parameters (mM)'); LEGEND_CONSTANTS(:,82) = strpad('K_D_CaM in component model_parameters (mM4)'); LEGEND_CONSTANTS(:,83) = strpad('V_init in component initial_conditions (mV)'); LEGEND_CONSTANTS(:,84) = strpad('Cai_total_init in component initial_conditions (mM)'); LEGEND_CONSTANTS(:,85) = strpad('Cai_init in component initial_conditions (mM)'); LEGEND_CONSTANTS(:,86) = strpad('Ki_init in component initial_conditions (mM)'); LEGEND_CONSTANTS(:,87) = strpad('Nai_init in component initial_conditions (mM)'); LEGEND_ALGEBRAIC(:,55) = strpad('norm in component IBK_channel_states (dimensionless)'); LEGEND_ALGEBRAIC(:,56) = strpad('norm in component INa_channel_states (dimensionless)'); LEGEND_CONSTANTS(:,88) = strpad('PO in component INS_Na (dimensionless)'); LEGEND_CONSTANTS(:,89) = strpad('PO in component INS_K (dimensionless)'); LEGEND_CONSTANTS(:,126) = strpad('E in component Nernst_potential (mV)'); LEGEND_STATES(:,8) = strpad('C0 in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,9) = strpad('C1 in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,10) = strpad('C2 in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,11) = strpad('C3 in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,12) = strpad('C0Ca in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,13) = strpad('C1Ca in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,14) = strpad('C2Ca in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,15) = strpad('C3Ca in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,16) = strpad('ICa in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,17) = strpad('IVs in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,18) = strpad('IVf in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,19) = strpad('IVfCa in component ICaL_channel_states (dimensionless)'); LEGEND_STATES(:,20) = strpad('IVsCa in component ICaL_channel_states (dimensionless)'); LEGEND_ALGEBRAIC(:,3) = strpad('alpha in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,11) = strpad('beta in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,16) = strpad('alpha_0 in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,20) = strpad('alpha_1 in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,24) = strpad('alpha_2 in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,28) = strpad('alpha_3 in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,33) = strpad('beta_0 in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,38) = strpad('beta_1 in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,42) = strpad('beta_2 in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,46) = strpad('beta_3 in component ICaL_channel_states (per_ms)'); LEGEND_CONSTANTS(:,127) = strpad('delta in component ICaL_channel_states (per_ms)'); LEGEND_CONSTANTS(:,128) = strpad('theta in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,54) = strpad('phi_s in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,50) = strpad('phi_f in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,62) = strpad('omega_s in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,61) = strpad('omega_f in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,64) = strpad('omega_fs in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,63) = strpad('omega_sf in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,60) = strpad('psi_s in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,59) = strpad('psi_f in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,58) = strpad('xi_s in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,57) = strpad('xi_f in component ICaL_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,65) = strpad('norm in component ICaL_channel_states (dimensionless)'); LEGEND_CONSTANTS(:,125) = strpad('phi in component temperature_factor (dimensionless)'); LEGEND_CONSTANTS(:,151) = strpad('E in component Nernst_potential (mV)'); LEGEND_ALGEBRAIC(:,4) = strpad('X_inf in component d_inf (dimensionless)'); LEGEND_CONSTANTS(:,129) = strpad('phi in component temperature_factor (dimensionless)'); LEGEND_ALGEBRAIC(:,5) = strpad('X_inf in component f_inf (dimensionless)'); LEGEND_ALGEBRAIC(:,12) = strpad('tau in component tau_f (ms)'); LEGEND_CONSTANTS(:,130) = strpad('phi in component temperature_factor (dimensionless)'); LEGEND_CONSTANTS(:,152) = strpad('E in component Nernst_potential (mV)'); LEGEND_ALGEBRAIC(:,6) = strpad('X_inf in component x_inf (dimensionless)'); LEGEND_CONSTANTS(:,131) = strpad('phi in component temperature_factor (dimensionless)'); LEGEND_ALGEBRAIC(:,7) = strpad('X_inf in component y_inf (dimensionless)'); LEGEND_CONSTANTS(:,132) = strpad('phi in component temperature_factor (dimensionless)'); LEGEND_CONSTANTS(:,133) = strpad('E in component Nernst_potential (mV)'); LEGEND_STATES(:,21) = strpad('C0 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,22) = strpad('C1 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,23) = strpad('C2 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,24) = strpad('C3 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,25) = strpad('C4 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,26) = strpad('O0 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,27) = strpad('O1 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,28) = strpad('O2 in component IBK_channel_states (dimensionless)'); LEGEND_STATES(:,29) = strpad('O3 in component IBK_channel_states (dimensionless)'); LEGEND_ALGEBRAIC(:,8) = strpad('alpha in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,13) = strpad('beta in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,17) = strpad('k_C0O0 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,21) = strpad('k_C1O1 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,25) = strpad('k_C2O2 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,29) = strpad('k_C3O3 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,34) = strpad('k_C4O4 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,36) = strpad('k_O0C0 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,40) = strpad('k_O1C1 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,44) = strpad('k_O2C2 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,48) = strpad('k_O3C3 in component IBK_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,52) = strpad('k_O4C4 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,134) = strpad('k_C0C1 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,135) = strpad('k_C1C2 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,136) = strpad('k_C2C3 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,137) = strpad('k_C3C4 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,138) = strpad('k_C4C3 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,139) = strpad('k_C3C2 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,140) = strpad('k_C2C1 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,141) = strpad('k_C1C0 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,142) = strpad('k_O0O1 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,143) = strpad('k_O1O2 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,144) = strpad('k_O2O3 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,145) = strpad('k_O3O4 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,146) = strpad('k_O4O3 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,147) = strpad('k_O3O2 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,148) = strpad('k_O2O1 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,149) = strpad('k_O1O0 in component IBK_channel_states (per_ms)'); LEGEND_CONSTANTS(:,150) = strpad('E in component Nernst_potential (mV)'); LEGEND_STATES(:,30) = strpad('C1 in component INa_channel_states (dimensionless)'); LEGEND_STATES(:,31) = strpad('C2 in component INa_channel_states (dimensionless)'); LEGEND_STATES(:,32) = strpad('C3 in component INa_channel_states (dimensionless)'); LEGEND_STATES(:,33) = strpad('I1 in component INa_channel_states (dimensionless)'); LEGEND_STATES(:,34) = strpad('I2 in component INa_channel_states (dimensionless)'); LEGEND_ALGEBRAIC(:,9) = strpad('k_I2I1 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,14) = strpad('k_I1O in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,18) = strpad('k_OC1 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,22) = strpad('k_C1C2 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,26) = strpad('k_C2C3 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,30) = strpad('k_C3C2 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,32) = strpad('k_C2C1 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,37) = strpad('k_C1O in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,41) = strpad('k_OI1 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,45) = strpad('k_I1I2 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,49) = strpad('k_I1C1 in component INa_channel_states (per_ms)'); LEGEND_ALGEBRAIC(:,53) = strpad('k_C1I1 in component INa_channel_states (per_ms)'); LEGEND_CONSTANTS(:,153) = strpad('phi in component temperature_factor (dimensionless)'); LEGEND_CONSTANTS(:,154) = strpad('E in component Nernst_potential (mV)'); LEGEND_CONSTANTS(:,155) = strpad('E in component Nernst_potential (mV)'); LEGEND_CONSTANTS(:,90) = strpad('R in component constants (J_per_K_mol)'); LEGEND_CONSTANTS(:,91) = strpad('F in component constants (C_per_mmol)'); LEGEND_CONSTANTS(:,92) = strpad('z_Ca in component constants (dimensionless)'); LEGEND_CONSTANTS(:,93) = strpad('z_Na in component constants (dimensionless)'); LEGEND_CONSTANTS(:,94) = strpad('z_K in component constants (dimensionless)'); LEGEND_CONSTANTS(:,95) = strpad('R in component constants (J_per_K_mol)'); LEGEND_CONSTANTS(:,96) = strpad('F in component constants (C_per_mmol)'); LEGEND_CONSTANTS(:,97) = strpad('z_Ca in component constants (dimensionless)'); LEGEND_CONSTANTS(:,98) = strpad('z_Na in component constants (dimensionless)'); LEGEND_CONSTANTS(:,99) = strpad('z_K in component constants (dimensionless)'); LEGEND_CONSTANTS(:,100) = strpad('R in component constants (J_per_K_mol)'); LEGEND_CONSTANTS(:,101) = strpad('F in component constants (C_per_mmol)'); LEGEND_CONSTANTS(:,102) = strpad('z_K in component constants (dimensionless)'); LEGEND_CONSTANTS(:,103) = strpad('z_Na in component constants (dimensionless)'); LEGEND_CONSTANTS(:,104) = strpad('z_Ca in component constants (dimensionless)'); LEGEND_CONSTANTS(:,105) = strpad('R in component constants (J_per_K_mol)'); LEGEND_CONSTANTS(:,106) = strpad('F in component constants (C_per_mmol)'); LEGEND_CONSTANTS(:,107) = strpad('z_K in component constants (dimensionless)'); LEGEND_CONSTANTS(:,108) = strpad('z_Na in component constants (dimensionless)'); LEGEND_CONSTANTS(:,109) = strpad('z_Ca in component constants (dimensionless)'); LEGEND_CONSTANTS(:,110) = strpad('R in component constants (J_per_K_mol)'); LEGEND_CONSTANTS(:,111) = strpad('F in component constants (C_per_mmol)'); LEGEND_CONSTANTS(:,112) = strpad('z_Na in component constants (dimensionless)'); LEGEND_CONSTANTS(:,113) = strpad('z_Ca in component constants (dimensionless)'); LEGEND_CONSTANTS(:,114) = strpad('z_K in component constants (dimensionless)'); LEGEND_CONSTANTS(:,115) = strpad('R in component constants (J_per_K_mol)'); LEGEND_CONSTANTS(:,116) = strpad('F in component constants (C_per_mmol)'); LEGEND_CONSTANTS(:,117) = strpad('z_Na in component constants (dimensionless)'); LEGEND_CONSTANTS(:,118) = strpad('z_Ca in component constants (dimensionless)'); LEGEND_CONSTANTS(:,119) = strpad('z_K in component constants (dimensionless)'); LEGEND_CONSTANTS(:,120) = strpad('R in component constants (J_per_K_mol)'); LEGEND_CONSTANTS(:,121) = strpad('F in component constants (C_per_mmol)'); LEGEND_CONSTANTS(:,122) = strpad('z_K in component constants (dimensionless)'); LEGEND_CONSTANTS(:,123) = strpad('z_Na in component constants (dimensionless)'); LEGEND_CONSTANTS(:,124) = strpad('z_Ca in component constants (dimensionless)'); LEGEND_RATES(:,11) = strpad('d/dt C3 in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,10) = strpad('d/dt C2 in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,9) = strpad('d/dt C1 in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,8) = strpad('d/dt C0 in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,15) = strpad('d/dt C3Ca in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,14) = strpad('d/dt C2Ca in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,13) = strpad('d/dt C1Ca in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,12) = strpad('d/dt C0Ca in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,1) = strpad('d/dt O_CaL in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,16) = strpad('d/dt ICa in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,17) = strpad('d/dt IVs in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,18) = strpad('d/dt IVf in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,20) = strpad('d/dt IVsCa in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,19) = strpad('d/dt IVfCa in component ICaL_channel_states (dimensionless)'); LEGEND_RATES(:,4) = strpad('d/dt X in component gating_kinetics (dimensionless)'); LEGEND_RATES(:,5) = strpad('d/dt X in component gating_kinetics (dimensionless)'); LEGEND_RATES(:,6) = strpad('d/dt X in component gating_kinetics (dimensionless)'); LEGEND_RATES(:,7) = strpad('d/dt X in component gating_kinetics (dimensionless)'); LEGEND_RATES(:,25) = strpad('d/dt C4 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,24) = strpad('d/dt C3 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,23) = strpad('d/dt C2 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,22) = strpad('d/dt C1 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,21) = strpad('d/dt C0 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,2) = strpad('d/dt O4 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,29) = strpad('d/dt O3 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,28) = strpad('d/dt O2 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,27) = strpad('d/dt O1 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,26) = strpad('d/dt O0 in component IBK_channel_states (dimensionless)'); LEGEND_RATES(:,32) = strpad('d/dt C3 in component INa_channel_states (dimensionless)'); LEGEND_RATES(:,31) = strpad('d/dt C2 in component INa_channel_states (dimensionless)'); LEGEND_RATES(:,30) = strpad('d/dt C1 in component INa_channel_states (dimensionless)'); LEGEND_RATES(:,3) = strpad('d/dt O_Na in component INa_channel_states (dimensionless)'); LEGEND_RATES(:,34) = strpad('d/dt I2 in component INa_channel_states (dimensionless)'); LEGEND_RATES(:,33) = strpad('d/dt I1 in component INa_channel_states (dimensionless)'); 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) = 1000; CONSTANTS(:,2) = -100; CONSTANTS(:,3) = -90; CONSTANTS(:,4) = 5.38843941249284e-5; CONSTANTS(:,5) = 153.604280337996; CONSTANTS(:,6) = 10.5731241425458; CONSTANTS(:,7) = 1.44; CONSTANTS(:,8) = 8.314; CONSTANTS(:,9) = 310; CONSTANTS(:,10) = 96.48534; CONSTANTS(:,11) = 310; CONSTANTS(:,12) = 297; CONSTANTS(:,13) = 2.1; CONSTANTS(:,14) = 2; CONSTANTS(:,15) = 0.815464741971086; CONSTANTS(:,16) = 0.152399266235657; CONSTANTS(:,17) = 0.0106805060777161; CONSTANTS(:,18) = 0.000332673548872087; CONSTANTS(:,19) = 0.0175888495282545; CONSTANTS(:,20) = 0.00328711668724504; CONSTANTS(:,21) = 0.000230369020877669; CONSTANTS(:,22) = 7.1754726923539e-6; CONSTANTS(:,23) = 8.38123983500905e-8; CONSTANTS(:,24) = 1.1193313274705e-6; CONSTANTS(:,25) = 4.0998751301597e-6; CONSTANTS(:,26) = 8.84306615061238e-8; CONSTANTS(:,27) = 2.41429816075123e-8; CONSTANTS(:,28) = 3.88576045134351e-6; CONSTANTS(:,29) = 0.0425; CONSTANTS(:,30) = 0.0791635737410974; CONSTANTS(:,31) = 0.377831534375835; CONSTANTS(:,32) = 1.9508; CONSTANTS(:,33) = 1.0217; CONSTANTS(:,34) = 297; CONSTANTS(:,35) = 3.1; CONSTANTS(:,36) = 5.4; CONSTANTS(:,37) = 0.14714161078933; CONSTANTS(:,38) = 0.99994773314105; CONSTANTS(:,39) = 4.7803; CONSTANTS(:,40) = 763.7564; CONSTANTS(:,41) = 80; CONSTANTS(:,42) = 0.48379087935899; CONSTANTS(:,43) = 0.385183559520031; CONSTANTS(:,44) = 0.115002824567753; CONSTANTS(:,45) = 0.0152602714149774; CONSTANTS(:,46) = 0.000759264410974374; CONSTANTS(:,47) = 6.94960798375172e-7; CONSTANTS(:,48) = 5.55636826398253e-8; CONSTANTS(:,49) = 2.85143702125325e-8; CONSTANTS(:,50) = 1.59832806123435e-6; CONSTANTS(:,51) = 1.82113764497095e-6; CONSTANTS(:,52) = 40633; CONSTANTS(:,53) = 11; CONSTANTS(:,54) = 1.1; CONSTANTS(:,55) = 25.1; CONSTANTS(:,56) = 297; CONSTANTS(:,57) = 2.45; CONSTANTS(:,58) = 140; CONSTANTS(:,59) = 0.0119443135223679; CONSTANTS(:,60) = 0.0109545368437155; CONSTANTS(:,61) = 0.973782548650071; CONSTANTS(:,62) = 0.000126882921013389; CONSTANTS(:,63) = 0.00318975045717667; CONSTANTS(:,64) = 1.96760342050475e-6; CONSTANTS(:,65) = 1.38; CONSTANTS(:,66) = 87.5; CONSTANTS(:,67) = 0.1; CONSTANTS(:,68) = 0.35; CONSTANTS(:,69) = 1992.1865; CONSTANTS(:,70) = 1; CONSTANTS(:,71) = 40; CONSTANTS(:,72) = 16.197; CONSTANTS(:,73) = 0.022488; CONSTANTS(:,74) = 0.017512; CONSTANTS(:,75) = 50; CONSTANTS(:,76) = 3.5e-12; CONSTANTS(:,77) = 1; CONSTANTS(:,78) = 4; CONSTANTS(:,79) = 0.034; CONSTANTS(:,80) = 0.012; CONSTANTS(:,81) = 0.0009; CONSTANTS(:,82) = 0.0001; CONSTANTS(:,83) = -73.5049651455872; CONSTANTS(:,84) = 0.004914; CONSTANTS(:,85) = 5.38843941249284e-5; CONSTANTS(:,86) = 153.604280337996; CONSTANTS(:,87) = 10.5731241425458; CONSTANTS(:,88) = 1; CONSTANTS(:,89) = 1; CONSTANTS(:,90) = 8.314; CONSTANTS(:,91) = 96.48534; CONSTANTS(:,92) = 2; CONSTANTS(:,93) = 1; CONSTANTS(:,94) = 1; CONSTANTS(:,95) = 8.314; CONSTANTS(:,96) = 96.48534; CONSTANTS(:,97) = 2; CONSTANTS(:,98) = 1; CONSTANTS(:,99) = 1; CONSTANTS(:,100) = 8.314; CONSTANTS(:,101) = 96.48534; CONSTANTS(:,102) = 1; CONSTANTS(:,103) = 1; CONSTANTS(:,104) = 2; CONSTANTS(:,105) = 8.314; CONSTANTS(:,106) = 96.48534; CONSTANTS(:,107) = 1; CONSTANTS(:,108) = 1; CONSTANTS(:,109) = 2; CONSTANTS(:,110) = 8.314; CONSTANTS(:,111) = 96.48534; CONSTANTS(:,112) = 1; CONSTANTS(:,113) = 2; CONSTANTS(:,114) = 1; CONSTANTS(:,115) = 8.314; CONSTANTS(:,116) = 96.48534; CONSTANTS(:,117) = 1; CONSTANTS(:,118) = 2; CONSTANTS(:,119) = 1; CONSTANTS(:,120) = 8.314; CONSTANTS(:,121) = 96.48534; CONSTANTS(:,122) = 1; CONSTANTS(:,123) = 1; CONSTANTS(:,124) = 2; CONSTANTS(:,125) = power(CONSTANTS(:,13), (CONSTANTS(:,9) - CONSTANTS(:,11))./10.0000); CONSTANTS(:,126) = (( CONSTANTS(:,90).*CONSTANTS(:,9))./( CONSTANTS(:,92).*CONSTANTS(:,91))).*log(CONSTANTS(:,14)./CONSTANTS(:,4)); CONSTANTS(:,127) = 0.00000; CONSTANTS(:,128) = 0.00000; CONSTANTS(:,129) = power(CONSTANTS(:,13), (CONSTANTS(:,9) - CONSTANTS(:,12))./10.0000); CONSTANTS(:,130) = power(CONSTANTS(:,13), (CONSTANTS(:,9) - CONSTANTS(:,12))./10.0000); CONSTANTS(:,131) = power(CONSTANTS(:,35), (CONSTANTS(:,9) - CONSTANTS(:,34))./10.0000); CONSTANTS(:,132) = power(CONSTANTS(:,35), (CONSTANTS(:,9) - CONSTANTS(:,34))./10.0000); CONSTANTS(:,133) = (( CONSTANTS(:,105).*CONSTANTS(:,9))./( CONSTANTS(:,107).*CONSTANTS(:,106))).*log(CONSTANTS(:,36)./CONSTANTS(:,5)); CONSTANTS(:,134) = 4.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,135) = 3.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,136) = 2.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,137) = 1.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,138) = 4.00000.*CONSTANTS(:,53); CONSTANTS(:,139) = 3.00000.*CONSTANTS(:,53); CONSTANTS(:,140) = 2.00000.*CONSTANTS(:,53); CONSTANTS(:,141) = 1.00000.*CONSTANTS(:,53); CONSTANTS(:,142) = 4.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,143) = 3.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,144) = 2.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,145) = 1.00000.*CONSTANTS(:,52).*CONSTANTS(:,4); CONSTANTS(:,146) = 4.00000.*CONSTANTS(:,54); CONSTANTS(:,147) = 3.00000.*CONSTANTS(:,54); CONSTANTS(:,148) = 2.00000.*CONSTANTS(:,54); CONSTANTS(:,149) = 1.00000.*CONSTANTS(:,54); CONSTANTS(:,150) = (( CONSTANTS(:,110).*CONSTANTS(:,9))./( CONSTANTS(:,112).*CONSTANTS(:,111))).*log(CONSTANTS(:,58)./CONSTANTS(:,6)); CONSTANTS(:,151) = (( CONSTANTS(:,95).*CONSTANTS(:,9))./( CONSTANTS(:,97).*CONSTANTS(:,96))).*log(CONSTANTS(:,14)./CONSTANTS(:,4)); CONSTANTS(:,152) = (( CONSTANTS(:,100).*CONSTANTS(:,9))./( CONSTANTS(:,102).*CONSTANTS(:,101))).*log(CONSTANTS(:,36)./CONSTANTS(:,5)); CONSTANTS(:,153) = power(CONSTANTS(:,57), (CONSTANTS(:,9) - CONSTANTS(:,56))./10.0000); CONSTANTS(:,154) = (( CONSTANTS(:,115).*CONSTANTS(:,9))./( CONSTANTS(:,117).*CONSTANTS(:,116))).*log(CONSTANTS(:,58)./CONSTANTS(:,6)); CONSTANTS(:,155) = (( CONSTANTS(:,120).*CONSTANTS(:,9))./( CONSTANTS(:,122).*CONSTANTS(:,121))).*log(CONSTANTS(:,36)./CONSTANTS(:,5)); STATES(:,1) = CONSTANTS(:,28); STATES(:,2) = CONSTANTS(:,51); STATES(:,3) = CONSTANTS(:,64); STATES(:,4) = CONSTANTS(:,30); STATES(:,5) = CONSTANTS(:,31); STATES(:,6) = CONSTANTS(:,37); STATES(:,7) = CONSTANTS(:,38); STATES(:,8) = CONSTANTS(:,15); STATES(:,9) = CONSTANTS(:,16); STATES(:,10) = CONSTANTS(:,17); STATES(:,11) = CONSTANTS(:,18); STATES(:,12) = CONSTANTS(:,19); STATES(:,13) = CONSTANTS(:,20); STATES(:,14) = CONSTANTS(:,21); STATES(:,15) = CONSTANTS(:,22); STATES(:,16) = CONSTANTS(:,23); STATES(:,17) = CONSTANTS(:,24); STATES(:,18) = CONSTANTS(:,25); STATES(:,19) = CONSTANTS(:,26); STATES(:,20) = CONSTANTS(:,26); STATES(:,21) = CONSTANTS(:,42); STATES(:,22) = CONSTANTS(:,43); STATES(:,23) = CONSTANTS(:,44); STATES(:,24) = CONSTANTS(:,45); STATES(:,25) = CONSTANTS(:,46); STATES(:,26) = CONSTANTS(:,47); STATES(:,27) = CONSTANTS(:,48); STATES(:,28) = CONSTANTS(:,49); STATES(:,29) = CONSTANTS(:,50); STATES(:,30) = CONSTANTS(:,59); STATES(:,31) = CONSTANTS(:,60); STATES(:,32) = CONSTANTS(:,61); STATES(:,33) = CONSTANTS(:,62); STATES(:,34) = CONSTANTS(:,63); 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 ALGEBRAIC(:,1) = piecewise({VOI<CONSTANTS(:,1), CONSTANTS(:,2) }, CONSTANTS(:,3)); ALGEBRAIC(:,4) = 1.00000./(1.00000+exp( - (ALGEBRAIC(:,1)+60.5000)./5.30000)); RATES(:,4) = ((ALGEBRAIC(:,4) - STATES(:,4))./CONSTANTS(:,32)).*CONSTANTS(:,129); ALGEBRAIC(:,6) = 1.00000./(1.00000+exp( - (ALGEBRAIC(:,1)+43.0000)./17.3600)); RATES(:,6) = ((ALGEBRAIC(:,6) - STATES(:,6))./CONSTANTS(:,39)).*CONSTANTS(:,131); ALGEBRAIC(:,7) = 1.00000./(1.00000+exp((ALGEBRAIC(:,1) - 44.9000)./12.0096)); RATES(:,7) = ((ALGEBRAIC(:,7) - STATES(:,7))./CONSTANTS(:,40)).*CONSTANTS(:,132); ALGEBRAIC(:,5) = 1.00000./(1.00000+exp((ALGEBRAIC(:,1)+75.5000)./4.00000)); ALGEBRAIC(:,12) = 0.381170.*(8.60000+ 14.7000.*exp( - power(ALGEBRAIC(:,1)+50.0000, 2.00000)./( 900.000.*1.00000))); RATES(:,5) = ((ALGEBRAIC(:,5) - STATES(:,5))./ALGEBRAIC(:,12)).*CONSTANTS(:,130); ALGEBRAIC(:,55) = STATES(:,21)+STATES(:,22)+STATES(:,23)+STATES(:,24)+STATES(:,25)+STATES(:,26)+STATES(:,27)+STATES(:,28)+STATES(:,29)+STATES(:,2); ALGEBRAIC(:,8) = 1.00000.*exp(( 8.47188.*ALGEBRAIC(:,1))./( 1.00000.*CONSTANTS(:,9))); ALGEBRAIC(:,34) = 0.0443240.*ALGEBRAIC(:,8); ALGEBRAIC(:,13) = 1.00000.*exp(( - 7.77556.*ALGEBRAIC(:,1))./( 1.00000.*CONSTANTS(:,9))); ALGEBRAIC(:,52) = 0.000106000.*ALGEBRAIC(:,13); RATES(:,25) = ( - (CONSTANTS(:,138)+ALGEBRAIC(:,34)).*STATES(:,25))./ALGEBRAIC(:,55)+( CONSTANTS(:,137).*STATES(:,24))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,52).*STATES(:,2))./ALGEBRAIC(:,55); ALGEBRAIC(:,29) = 0.000781000.*ALGEBRAIC(:,8); ALGEBRAIC(:,48) = 0.0953120.*ALGEBRAIC(:,13); RATES(:,24) = ( - (CONSTANTS(:,139)+ALGEBRAIC(:,29)+CONSTANTS(:,137)).*STATES(:,24))./ALGEBRAIC(:,55)+( CONSTANTS(:,136).*STATES(:,23))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,48).*STATES(:,29))./ALGEBRAIC(:,55)+( CONSTANTS(:,138).*STATES(:,25))./ALGEBRAIC(:,55); ALGEBRAIC(:,25) = 2.81000e-05.*ALGEBRAIC(:,8); ALGEBRAIC(:,44) = 32.6594.*ALGEBRAIC(:,13); RATES(:,23) = ( - (CONSTANTS(:,140)+ALGEBRAIC(:,25)+CONSTANTS(:,136)).*STATES(:,23))./ALGEBRAIC(:,55)+( CONSTANTS(:,135).*STATES(:,22))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,44).*STATES(:,28))./ALGEBRAIC(:,55)+( CONSTANTS(:,139).*STATES(:,24))./ALGEBRAIC(:,55); ALGEBRAIC(:,21) = 0.000869000.*ALGEBRAIC(:,8); ALGEBRAIC(:,40) = 144.174.*ALGEBRAIC(:,13); RATES(:,22) = ( - (CONSTANTS(:,141)+ALGEBRAIC(:,21)+CONSTANTS(:,135)).*STATES(:,22))./ALGEBRAIC(:,55)+( CONSTANTS(:,134).*STATES(:,21))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,40).*STATES(:,27))./ALGEBRAIC(:,55)+( CONSTANTS(:,140).*STATES(:,23))./ALGEBRAIC(:,55); ALGEBRAIC(:,17) = 0.0216200.*ALGEBRAIC(:,8); ALGEBRAIC(:,36) = 318.108.*ALGEBRAIC(:,13); RATES(:,21) = ( - (CONSTANTS(:,134)+ALGEBRAIC(:,17)).*STATES(:,21))./ALGEBRAIC(:,55)+( CONSTANTS(:,141).*STATES(:,22))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,36).*STATES(:,26))./ALGEBRAIC(:,55); RATES(:,2) = ( - (CONSTANTS(:,146)+ALGEBRAIC(:,52)).*STATES(:,2))./ALGEBRAIC(:,55)+( CONSTANTS(:,145).*STATES(:,29))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,34).*STATES(:,25))./ALGEBRAIC(:,55); RATES(:,29) = ( - (CONSTANTS(:,147)+ALGEBRAIC(:,48)+CONSTANTS(:,145)).*STATES(:,29))./ALGEBRAIC(:,55)+( CONSTANTS(:,144).*STATES(:,28))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,29).*STATES(:,24))./ALGEBRAIC(:,55)+( CONSTANTS(:,146).*STATES(:,2))./ALGEBRAIC(:,55); RATES(:,28) = ( - (CONSTANTS(:,148)+ALGEBRAIC(:,44)+CONSTANTS(:,144)).*STATES(:,28))./ALGEBRAIC(:,55)+( CONSTANTS(:,143).*STATES(:,27))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,25).*STATES(:,23))./ALGEBRAIC(:,55)+( CONSTANTS(:,147).*STATES(:,29))./ALGEBRAIC(:,55); RATES(:,27) = ( - (CONSTANTS(:,149)+ALGEBRAIC(:,40)+CONSTANTS(:,143)).*STATES(:,27))./ALGEBRAIC(:,55)+( CONSTANTS(:,142).*STATES(:,26))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,21).*STATES(:,22))./ALGEBRAIC(:,55)+( CONSTANTS(:,148).*STATES(:,28))./ALGEBRAIC(:,55); RATES(:,26) = ( - (CONSTANTS(:,142)+ALGEBRAIC(:,36)).*STATES(:,26))./ALGEBRAIC(:,55)+( CONSTANTS(:,149).*STATES(:,27))./ALGEBRAIC(:,55)+( ALGEBRAIC(:,17).*STATES(:,21))./ALGEBRAIC(:,55); ALGEBRAIC(:,56) = STATES(:,3)+STATES(:,30)+STATES(:,31)+STATES(:,32)+STATES(:,33)+STATES(:,34); ALGEBRAIC(:,26) = CONSTANTS(:,153).*0.554320.*exp( - 0.0990740+ 0.0364410.*ALGEBRAIC(:,1)); ALGEBRAIC(:,30) = CONSTANTS(:,153).*0.000525480.*exp( - 0.0691020+ 0.00319450.*ALGEBRAIC(:,1)); RATES(:,32) = ( - ALGEBRAIC(:,30).*STATES(:,32))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,26).*STATES(:,31))./ALGEBRAIC(:,56); ALGEBRAIC(:,22) = CONSTANTS(:,153).*3.15660.*exp(0.363520+ 0.0771930.*ALGEBRAIC(:,1)); ALGEBRAIC(:,32) = CONSTANTS(:,153).*1.44960.*exp( - 0.156600+ 0.0583530.*ALGEBRAIC(:,1)); RATES(:,31) = ( - (ALGEBRAIC(:,32)+ALGEBRAIC(:,26)).*STATES(:,31))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,22).*STATES(:,30))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,30).*STATES(:,32))./ALGEBRAIC(:,56); ALGEBRAIC(:,18) = CONSTANTS(:,153).*2.39100.*exp( - 13.3350 - 0.252890.*ALGEBRAIC(:,1)); ALGEBRAIC(:,37) = CONSTANTS(:,153).*1.53290.*exp(0.00931930+ 0.0410750.*ALGEBRAIC(:,1)); ALGEBRAIC(:,49) = CONSTANTS(:,153).*1.90460.*exp( - 2.48400+ 0.0204060.*ALGEBRAIC(:,1)); ALGEBRAIC(:,53) = CONSTANTS(:,153).*0.000216880.*exp( - 0.0634380+ 0.00466830.*ALGEBRAIC(:,1)); RATES(:,30) = ( - (ALGEBRAIC(:,22)+ALGEBRAIC(:,37)+ALGEBRAIC(:,53)).*STATES(:,30))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,18).*STATES(:,3))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,32).*STATES(:,31))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,49).*STATES(:,33))./ALGEBRAIC(:,56); ALGEBRAIC(:,14) = CONSTANTS(:,153).*0.120520.*exp( - 9.60280+ 0.0830250.*ALGEBRAIC(:,1)); ALGEBRAIC(:,41) = CONSTANTS(:,153).*1.61640.*exp(0.307630+ 0.00605350.*ALGEBRAIC(:,1)); RATES(:,3) = ( - (ALGEBRAIC(:,18)+ALGEBRAIC(:,41)).*STATES(:,3))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,37).*STATES(:,30))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,14).*STATES(:,33))./ALGEBRAIC(:,56); ALGEBRAIC(:,9) = CONSTANTS(:,153).*0.00392390.*exp(2.67930+ 0.00614680.*ALGEBRAIC(:,1)); ALGEBRAIC(:,45) = CONSTANTS(:,153).*0.0277350.*exp(0.0514900 - 0.0468650.*ALGEBRAIC(:,1)); RATES(:,34) = ( - ALGEBRAIC(:,9).*STATES(:,34))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,45).*STATES(:,33))./ALGEBRAIC(:,56); RATES(:,33) = ( - (ALGEBRAIC(:,14)+ALGEBRAIC(:,45)+ALGEBRAIC(:,49)).*STATES(:,33))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,9).*STATES(:,34))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,53).*STATES(:,30))./ALGEBRAIC(:,56)+( ALGEBRAIC(:,41).*STATES(:,3))./ALGEBRAIC(:,56); ALGEBRAIC(:,3) = CONSTANTS(:,125).*0.731000.*exp(ALGEBRAIC(:,1)./30.0000); ALGEBRAIC(:,24) = 2.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,28) = 1.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,11) = CONSTANTS(:,125).*0.214900.*exp( - ALGEBRAIC(:,1)./40.0000); ALGEBRAIC(:,42) = 3.00000.*ALGEBRAIC(:,11); ALGEBRAIC(:,46) = 4.00000.*ALGEBRAIC(:,11); ALGEBRAIC(:,54) = CONSTANTS(:,125).*0.00595600.*exp( - ALGEBRAIC(:,1)./40.0000); ALGEBRAIC(:,50) = CONSTANTS(:,125).*0.474200.*exp(ALGEBRAIC(:,1)./10.0000); ALGEBRAIC(:,60) = CONSTANTS(:,125).*0.00232000.*exp( - ALGEBRAIC(:,1)./280.000); ALGEBRAIC(:,58) = CONSTANTS(:,125).*0.0121300.*exp(ALGEBRAIC(:,1)./500.000); ALGEBRAIC(:,62) = ( ALGEBRAIC(:,46).*ALGEBRAIC(:,58).*ALGEBRAIC(:,54))./( ALGEBRAIC(:,28).*ALGEBRAIC(:,60)); ALGEBRAIC(:,59) = CONSTANTS(:,125).*0.0219700.*exp(ALGEBRAIC(:,1)./500.000); ALGEBRAIC(:,57) = CONSTANTS(:,125).*0.0140700.*exp( - ALGEBRAIC(:,1)./300.000); ALGEBRAIC(:,61) = ( ALGEBRAIC(:,46).*ALGEBRAIC(:,57).*ALGEBRAIC(:,50))./( ALGEBRAIC(:,28).*ALGEBRAIC(:,59)); ALGEBRAIC(:,65) = STATES(:,11)+STATES(:,10)+STATES(:,9)+STATES(:,8)+STATES(:,15)+STATES(:,14)+STATES(:,13)+STATES(:,12)+STATES(:,1)+STATES(:,16)+STATES(:,17)+STATES(:,18)+STATES(:,20)+STATES(:,19); RATES(:,11) = (( ALGEBRAIC(:,24).*STATES(:,10))./ALGEBRAIC(:,65) - ( (ALGEBRAIC(:,28)+ALGEBRAIC(:,42)+ALGEBRAIC(:,54)+ALGEBRAIC(:,50)+CONSTANTS(:,128)).*STATES(:,11))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,46).*STATES(:,1))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,61).*STATES(:,18))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,62).*STATES(:,17))./ALGEBRAIC(:,65)+( CONSTANTS(:,127).*STATES(:,15))./ALGEBRAIC(:,65); ALGEBRAIC(:,20) = 3.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,38) = 2.00000.*ALGEBRAIC(:,11); RATES(:,10) = (( ALGEBRAIC(:,20).*STATES(:,9))./ALGEBRAIC(:,65) - ( (ALGEBRAIC(:,24)+ALGEBRAIC(:,38)+CONSTANTS(:,128)).*STATES(:,10))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,42).*STATES(:,11))./ALGEBRAIC(:,65)+( CONSTANTS(:,127).*STATES(:,14))./ALGEBRAIC(:,65); ALGEBRAIC(:,16) = 4.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,33) = 1.00000.*ALGEBRAIC(:,11); RATES(:,9) = (( ALGEBRAIC(:,16).*STATES(:,8))./ALGEBRAIC(:,65) - ( (ALGEBRAIC(:,20)+ALGEBRAIC(:,33)+CONSTANTS(:,128)).*STATES(:,9))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,38).*STATES(:,10))./ALGEBRAIC(:,65)+( CONSTANTS(:,127).*STATES(:,13))./ALGEBRAIC(:,65); RATES(:,8) = ( - (ALGEBRAIC(:,16)+CONSTANTS(:,128)).*STATES(:,8))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,33).*STATES(:,9))./ALGEBRAIC(:,65)+( CONSTANTS(:,127).*STATES(:,12))./ALGEBRAIC(:,65); RATES(:,15) = ((( CONSTANTS(:,128).*STATES(:,11))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,24).*STATES(:,14))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,28)+ALGEBRAIC(:,42)+ALGEBRAIC(:,50)+ALGEBRAIC(:,54)+CONSTANTS(:,127)).*STATES(:,15))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,46).*STATES(:,16))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,61).*STATES(:,19))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,62).*STATES(:,20))./ALGEBRAIC(:,65); RATES(:,14) = ((( CONSTANTS(:,128).*STATES(:,10))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,20).*STATES(:,13))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,24)+ALGEBRAIC(:,38)+CONSTANTS(:,127)).*STATES(:,14))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,42).*STATES(:,15))./ALGEBRAIC(:,65); RATES(:,13) = ((( CONSTANTS(:,128).*STATES(:,9))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,16).*STATES(:,12))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,20)+ALGEBRAIC(:,33)+CONSTANTS(:,127)).*STATES(:,13))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,38).*STATES(:,14))./ALGEBRAIC(:,65); RATES(:,12) = (( CONSTANTS(:,128).*STATES(:,8))./ALGEBRAIC(:,65) - ( (ALGEBRAIC(:,16)+CONSTANTS(:,127)).*STATES(:,12))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,33).*STATES(:,13))./ALGEBRAIC(:,65); RATES(:,1) = (( ALGEBRAIC(:,28).*STATES(:,11))./ALGEBRAIC(:,65) - ( (ALGEBRAIC(:,46)+ALGEBRAIC(:,59)+ALGEBRAIC(:,60)+CONSTANTS(:,128)).*STATES(:,1))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,57).*STATES(:,18))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,58).*STATES(:,17))./ALGEBRAIC(:,65)+( CONSTANTS(:,127).*STATES(:,16))./ALGEBRAIC(:,65); RATES(:,16) = ((( CONSTANTS(:,128).*STATES(:,1))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,28).*STATES(:,15))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,46)+ALGEBRAIC(:,59)+ALGEBRAIC(:,60)+CONSTANTS(:,127)).*STATES(:,16))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,57).*STATES(:,19))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,58).*STATES(:,20))./ALGEBRAIC(:,65); ALGEBRAIC(:,64) = ALGEBRAIC(:,60); ALGEBRAIC(:,63) = ( ALGEBRAIC(:,58).*ALGEBRAIC(:,59))./ALGEBRAIC(:,57); RATES(:,17) = ((( ALGEBRAIC(:,54).*STATES(:,11))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,60).*STATES(:,1))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,64).*STATES(:,18))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,62)+ALGEBRAIC(:,58)+ALGEBRAIC(:,63)+CONSTANTS(:,128)).*STATES(:,17))./ALGEBRAIC(:,65))+( CONSTANTS(:,127).*STATES(:,20))./ALGEBRAIC(:,65); RATES(:,18) = ((( ALGEBRAIC(:,50).*STATES(:,11))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,59).*STATES(:,1))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,61)+ALGEBRAIC(:,57)+ALGEBRAIC(:,64)+CONSTANTS(:,128)).*STATES(:,18))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,63).*STATES(:,17))./ALGEBRAIC(:,65)+( CONSTANTS(:,127).*STATES(:,19))./ALGEBRAIC(:,65); RATES(:,20) = (( CONSTANTS(:,128).*STATES(:,17))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,54).*STATES(:,15))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,60).*STATES(:,16))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,64).*STATES(:,19))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,62)+ALGEBRAIC(:,58)+ALGEBRAIC(:,63)+CONSTANTS(:,127)).*STATES(:,20))./ALGEBRAIC(:,65); RATES(:,19) = ((( CONSTANTS(:,128).*STATES(:,18))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,50).*STATES(:,15))./ALGEBRAIC(:,65)+( ALGEBRAIC(:,59).*STATES(:,16))./ALGEBRAIC(:,65)) - ( (ALGEBRAIC(:,61)+ALGEBRAIC(:,57)+ALGEBRAIC(:,64)+CONSTANTS(:,127)).*STATES(:,19))./ALGEBRAIC(:,65))+( ALGEBRAIC(:,63).*STATES(:,20))./ALGEBRAIC(:,65); RATES = RATES'; end % Calculate algebraic variables function ALGEBRAIC = computeAlgebraic(ALGEBRAIC, CONSTANTS, STATES, VOI) statesSize = size(STATES); statesColumnCount = statesSize(2); if ( statesColumnCount == 1) STATES = STATES'; utilOnes = 1; else statesRowCount = statesSize(1); utilOnes = ones(statesRowCount, 1); end ALGEBRAIC(:,1) = piecewise({VOI<CONSTANTS(:,1), CONSTANTS(:,2) }, CONSTANTS(:,3)); ALGEBRAIC(:,4) = 1.00000./(1.00000+exp( - (ALGEBRAIC(:,1)+60.5000)./5.30000)); ALGEBRAIC(:,6) = 1.00000./(1.00000+exp( - (ALGEBRAIC(:,1)+43.0000)./17.3600)); ALGEBRAIC(:,7) = 1.00000./(1.00000+exp((ALGEBRAIC(:,1) - 44.9000)./12.0096)); ALGEBRAIC(:,5) = 1.00000./(1.00000+exp((ALGEBRAIC(:,1)+75.5000)./4.00000)); ALGEBRAIC(:,12) = 0.381170.*(8.60000+ 14.7000.*exp( - power(ALGEBRAIC(:,1)+50.0000, 2.00000)./( 900.000.*1.00000))); ALGEBRAIC(:,55) = STATES(:,21)+STATES(:,22)+STATES(:,23)+STATES(:,24)+STATES(:,25)+STATES(:,26)+STATES(:,27)+STATES(:,28)+STATES(:,29)+STATES(:,2); ALGEBRAIC(:,8) = 1.00000.*exp(( 8.47188.*ALGEBRAIC(:,1))./( 1.00000.*CONSTANTS(:,9))); ALGEBRAIC(:,34) = 0.0443240.*ALGEBRAIC(:,8); ALGEBRAIC(:,13) = 1.00000.*exp(( - 7.77556.*ALGEBRAIC(:,1))./( 1.00000.*CONSTANTS(:,9))); ALGEBRAIC(:,52) = 0.000106000.*ALGEBRAIC(:,13); ALGEBRAIC(:,29) = 0.000781000.*ALGEBRAIC(:,8); ALGEBRAIC(:,48) = 0.0953120.*ALGEBRAIC(:,13); ALGEBRAIC(:,25) = 2.81000e-05.*ALGEBRAIC(:,8); ALGEBRAIC(:,44) = 32.6594.*ALGEBRAIC(:,13); ALGEBRAIC(:,21) = 0.000869000.*ALGEBRAIC(:,8); ALGEBRAIC(:,40) = 144.174.*ALGEBRAIC(:,13); ALGEBRAIC(:,17) = 0.0216200.*ALGEBRAIC(:,8); ALGEBRAIC(:,36) = 318.108.*ALGEBRAIC(:,13); ALGEBRAIC(:,56) = STATES(:,3)+STATES(:,30)+STATES(:,31)+STATES(:,32)+STATES(:,33)+STATES(:,34); ALGEBRAIC(:,26) = CONSTANTS(:,153).*0.554320.*exp( - 0.0990740+ 0.0364410.*ALGEBRAIC(:,1)); ALGEBRAIC(:,30) = CONSTANTS(:,153).*0.000525480.*exp( - 0.0691020+ 0.00319450.*ALGEBRAIC(:,1)); ALGEBRAIC(:,22) = CONSTANTS(:,153).*3.15660.*exp(0.363520+ 0.0771930.*ALGEBRAIC(:,1)); ALGEBRAIC(:,32) = CONSTANTS(:,153).*1.44960.*exp( - 0.156600+ 0.0583530.*ALGEBRAIC(:,1)); ALGEBRAIC(:,18) = CONSTANTS(:,153).*2.39100.*exp( - 13.3350 - 0.252890.*ALGEBRAIC(:,1)); ALGEBRAIC(:,37) = CONSTANTS(:,153).*1.53290.*exp(0.00931930+ 0.0410750.*ALGEBRAIC(:,1)); ALGEBRAIC(:,49) = CONSTANTS(:,153).*1.90460.*exp( - 2.48400+ 0.0204060.*ALGEBRAIC(:,1)); ALGEBRAIC(:,53) = CONSTANTS(:,153).*0.000216880.*exp( - 0.0634380+ 0.00466830.*ALGEBRAIC(:,1)); ALGEBRAIC(:,14) = CONSTANTS(:,153).*0.120520.*exp( - 9.60280+ 0.0830250.*ALGEBRAIC(:,1)); ALGEBRAIC(:,41) = CONSTANTS(:,153).*1.61640.*exp(0.307630+ 0.00605350.*ALGEBRAIC(:,1)); ALGEBRAIC(:,9) = CONSTANTS(:,153).*0.00392390.*exp(2.67930+ 0.00614680.*ALGEBRAIC(:,1)); ALGEBRAIC(:,45) = CONSTANTS(:,153).*0.0277350.*exp(0.0514900 - 0.0468650.*ALGEBRAIC(:,1)); ALGEBRAIC(:,3) = CONSTANTS(:,125).*0.731000.*exp(ALGEBRAIC(:,1)./30.0000); ALGEBRAIC(:,24) = 2.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,28) = 1.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,11) = CONSTANTS(:,125).*0.214900.*exp( - ALGEBRAIC(:,1)./40.0000); ALGEBRAIC(:,42) = 3.00000.*ALGEBRAIC(:,11); ALGEBRAIC(:,46) = 4.00000.*ALGEBRAIC(:,11); ALGEBRAIC(:,54) = CONSTANTS(:,125).*0.00595600.*exp( - ALGEBRAIC(:,1)./40.0000); ALGEBRAIC(:,50) = CONSTANTS(:,125).*0.474200.*exp(ALGEBRAIC(:,1)./10.0000); ALGEBRAIC(:,60) = CONSTANTS(:,125).*0.00232000.*exp( - ALGEBRAIC(:,1)./280.000); ALGEBRAIC(:,58) = CONSTANTS(:,125).*0.0121300.*exp(ALGEBRAIC(:,1)./500.000); ALGEBRAIC(:,62) = ( ALGEBRAIC(:,46).*ALGEBRAIC(:,58).*ALGEBRAIC(:,54))./( ALGEBRAIC(:,28).*ALGEBRAIC(:,60)); ALGEBRAIC(:,59) = CONSTANTS(:,125).*0.0219700.*exp(ALGEBRAIC(:,1)./500.000); ALGEBRAIC(:,57) = CONSTANTS(:,125).*0.0140700.*exp( - ALGEBRAIC(:,1)./300.000); ALGEBRAIC(:,61) = ( ALGEBRAIC(:,46).*ALGEBRAIC(:,57).*ALGEBRAIC(:,50))./( ALGEBRAIC(:,28).*ALGEBRAIC(:,59)); ALGEBRAIC(:,65) = STATES(:,11)+STATES(:,10)+STATES(:,9)+STATES(:,8)+STATES(:,15)+STATES(:,14)+STATES(:,13)+STATES(:,12)+STATES(:,1)+STATES(:,16)+STATES(:,17)+STATES(:,18)+STATES(:,20)+STATES(:,19); ALGEBRAIC(:,20) = 3.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,38) = 2.00000.*ALGEBRAIC(:,11); ALGEBRAIC(:,16) = 4.00000.*ALGEBRAIC(:,3); ALGEBRAIC(:,33) = 1.00000.*ALGEBRAIC(:,11); ALGEBRAIC(:,64) = ALGEBRAIC(:,60); ALGEBRAIC(:,63) = ( ALGEBRAIC(:,58).*ALGEBRAIC(:,59))./ALGEBRAIC(:,57); ALGEBRAIC(:,2) = STATES(:,4).*STATES(:,5); ALGEBRAIC(:,10) = STATES(:,6).*STATES(:,7); ALGEBRAIC(:,15) = ( CONSTANTS(:,69).*( exp(( CONSTANTS(:,68).*ALGEBRAIC(:,1).*CONSTANTS(:,10))./( CONSTANTS(:,8).*CONSTANTS(:,9))).*power(CONSTANTS(:,6), 3.00000).*CONSTANTS(:,14) - 2.50000.*exp(( (CONSTANTS(:,68) - 1.00000).*ALGEBRAIC(:,1).*CONSTANTS(:,10))./( CONSTANTS(:,8).*CONSTANTS(:,9))).*power(CONSTANTS(:,58), 3.00000).*CONSTANTS(:,4)))./( (power(CONSTANTS(:,66), 3.00000)+power(CONSTANTS(:,58), 3.00000)).*(CONSTANTS(:,65)+CONSTANTS(:,14)).*(1.00000+ CONSTANTS(:,67).*exp(( (CONSTANTS(:,68) - 1.00000).*ALGEBRAIC(:,1).*CONSTANTS(:,10))./( CONSTANTS(:,8).*CONSTANTS(:,9))))); ALGEBRAIC(:,19) = ( CONSTANTS(:,72).*CONSTANTS(:,36).*CONSTANTS(:,6))./( (CONSTANTS(:,36)+CONSTANTS(:,70)).*(CONSTANTS(:,71)+CONSTANTS(:,6)).*(1.00000+ 0.124500.*exp(( - 0.100000.*ALGEBRAIC(:,1).*CONSTANTS(:,10))./( CONSTANTS(:,8).*CONSTANTS(:,9)))+ 0.0353000.*exp(( - ALGEBRAIC(:,1).*CONSTANTS(:,10))./( CONSTANTS(:,8).*CONSTANTS(:,9))))); ALGEBRAIC(:,23) = CONSTANTS(:,7).*STATES(:,1).*(ALGEBRAIC(:,1) - CONSTANTS(:,126)); ALGEBRAIC(:,27) = CONSTANTS(:,29).*ALGEBRAIC(:,2).*(ALGEBRAIC(:,1) - CONSTANTS(:,151)); ALGEBRAIC(:,31) = CONSTANTS(:,33).*ALGEBRAIC(:,10).*(ALGEBRAIC(:,1) - CONSTANTS(:,152)); ALGEBRAIC(:,35) = CONSTANTS(:,41).*STATES(:,2).*(ALGEBRAIC(:,1) - CONSTANTS(:,133)); ALGEBRAIC(:,39) = CONSTANTS(:,55).*STATES(:,3).*(ALGEBRAIC(:,1) - CONSTANTS(:,150)); ALGEBRAIC(:,43) = CONSTANTS(:,73).*CONSTANTS(:,88).*(ALGEBRAIC(:,1) - CONSTANTS(:,154)); ALGEBRAIC(:,47) = CONSTANTS(:,74).*CONSTANTS(:,89).*(ALGEBRAIC(:,1) - CONSTANTS(:,155)); ALGEBRAIC(:,51) = ALGEBRAIC(:,23)+ALGEBRAIC(:,27)+ALGEBRAIC(:,31)+ALGEBRAIC(:,35)+ALGEBRAIC(:,39)+ALGEBRAIC(:,15)+ALGEBRAIC(:,19)+ALGEBRAIC(:,43)+ALGEBRAIC(:,47); end % Compute result of a piecewise function function x = piecewise(cases, default) set = [0]; for i = 1:2:length(cases) if (length(cases{i+1}) == 1) x(cases{i} & ~set,:) = cases{i+1}; else x(cases{i} & ~set,:) = cases{i+1}(cases{i} & ~set); end set = set | cases{i}; if(set), break, end end if (length(default) == 1) x(~set,:) = default; else x(~set,:) = default(~set); 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