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 =369; end % There are a total of 145 entries in each of the rate and state variable arrays. % There are a total of 410 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(:,325) = strpad('AC47_cyt in component ac (uM)'); LEGEND_CONSTANTS(:,330) = strpad('AC47_eca in component ac (uM)'); LEGEND_CONSTANTS(:,321) = strpad('AC56_cav in component ac (uM)'); LEGEND_CONSTANTS(:,329) = strpad('AC56_cyt in component ac (uM)'); LEGEND_CONSTANTS(:,297) = strpad('AC_tot in component ac (uM)'); LEGEND_CONSTANTS(:,1) = strpad('ATP in component ac (uM)'); LEGEND_CONSTANTS(:,2) = strpad('KmATP in component ac (uM)'); LEGEND_CONSTANTS(:,3) = strpad('KmGiAC56 in component ac (uM)'); LEGEND_CONSTANTS(:,4) = strpad('KmGsAC47 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,5) = strpad('KmGsAC56 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,6) = strpad('KmGsGiAC56 in component ac (dimensionless)'); LEGEND_ALGEBRAIC(:,1) = strpad('ac_kAC47_cyt_gsa in component ac (dimensionless)'); LEGEND_ALGEBRAIC(:,2) = strpad('ac_kAC47_eca_gsa in component ac (dimensionless)'); LEGEND_ALGEBRAIC(:,3) = strpad('ac_kAC56_cav_gsa in component ac (dimensionless)'); LEGEND_ALGEBRAIC(:,31) = strpad('ac_kAC56_cyt_gsa in component ac (dimensionless)'); LEGEND_CONSTANTS(:,7) = strpad('afAC47 in component ac (hertz)'); LEGEND_CONSTANTS(:,8) = strpad('afAC56 in component ac (hertz)'); LEGEND_CONSTANTS(:,9) = strpad('basalAC47 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,10) = strpad('basalAC56 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,295) = strpad('R_b1_tot in component beta (uM)'); LEGEND_STATES(:,1) = strpad('Gi_bg in component beta_cav (uM)'); LEGEND_STATES(:,2) = strpad('Gs_aGTP in component beta_cav (uM)'); LEGEND_STATES(:,3) = strpad('Gs_aGTP in component beta_cyt (uM)'); LEGEND_STATES(:,4) = strpad('Gs_aGTP in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,69) = strpad('dcAMP_AC47_cyt in component ac (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,68) = strpad('dcAMP_AC47_eca in component ac (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,83) = strpad('dcAMP_AC56_cav in component ac (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,107) = strpad('dcAMP_AC56_cyt in component ac (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_CONSTANTS(:,228) = strpad('fATP in component ac (dimensionless)'); LEGEND_CONSTANTS(:,11) = strpad('f_AC47_eca in component ac (dimensionless)'); LEGEND_CONSTANTS(:,261) = strpad('f_AC56_AC47 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,12) = strpad('f_AC56_cav in component ac (dimensionless)'); LEGEND_ALGEBRAIC(:,33) = strpad('gsi in component ac (dimensionless)'); LEGEND_CONSTANTS(:,13) = strpad('hGsAC47 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,14) = strpad('hGsAC56 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,15) = strpad('hGsGiAC56 in component ac (dimensionless)'); LEGEND_ALGEBRAIC(:,56) = strpad('kAC47_cyt in component ac (hertz)'); LEGEND_ALGEBRAIC(:,32) = strpad('kAC47_eca in component ac (hertz)'); LEGEND_ALGEBRAIC(:,57) = strpad('kAC56_cav in component ac (hertz)'); LEGEND_ALGEBRAIC(:,85) = strpad('kAC56_cyt in component ac (hertz)'); LEGEND_CONSTANTS(:,16) = strpad('vGsGiAC56 in component ac (dimensionless)'); LEGEND_CONSTANTS(:,320) = strpad('vr_cav in component cell (dimensionless)'); LEGEND_CONSTANTS(:,322) = strpad('vr_cyt in component cell (dimensionless)'); LEGEND_CONSTANTS(:,327) = strpad('vr_eca in component cell (dimensionless)'); LEGEND_CONSTANTS(:,17) = strpad('ICaL_akap in component akap_sig (uM)'); LEGEND_CONSTANTS(:,392) = strpad('ICaL_akapf in component akap_sig (uM)'); LEGEND_CONSTANTS(:,393) = strpad('ICaL_arn in component akap_sig (uM)'); LEGEND_CONSTANTS(:,394) = strpad('ICaL_arp in component akap_sig (uM)'); LEGEND_CONSTANTS(:,18) = strpad('ICaL_tot in component akap_sig (uM)'); LEGEND_CONSTANTS(:,282) = strpad('ICaLf in component akap_sig (uM)'); LEGEND_STATES(:,5) = strpad('ICaLp in component akap_sig (uM)'); LEGEND_CONSTANTS(:,19) = strpad('Ka_ical in component akap_sig (uM)'); LEGEND_CONSTANTS(:,20) = strpad('Ka_ryr in component akap_sig (uM)'); LEGEND_CONSTANTS(:,21) = strpad('Ki in component akap_sig (uM)'); LEGEND_CONSTANTS(:,22) = strpad('Kp_ical in component akap_sig (uM)'); LEGEND_CONSTANTS(:,23) = strpad('Kp_ryr in component akap_sig (uM)'); LEGEND_CONSTANTS(:,24) = strpad('Kr in component akap_sig (uM)'); LEGEND_CONSTANTS(:,25) = strpad('Li in component akap_sig (uM)'); LEGEND_CONSTANTS(:,26) = strpad('Lr in component akap_sig (uM)'); LEGEND_CONSTANTS(:,27) = strpad('Mi in component akap_sig (uM)'); LEGEND_CONSTANTS(:,28) = strpad('Mr in component akap_sig (uM)'); LEGEND_CONSTANTS(:,381) = strpad('PKA_cav in component pka (uM)'); LEGEND_CONSTANTS(:,391) = strpad('PKAf in component akap_sig (uM)'); LEGEND_CONSTANTS(:,29) = strpad('PP1_cav in component pp1 (uM)'); LEGEND_CONSTANTS(:,292) = strpad('PP1f_cav in component akap_sig (uM)'); LEGEND_CONSTANTS(:,30) = strpad('RyR_akap in component akap_sig (uM)'); LEGEND_CONSTANTS(:,395) = strpad('RyR_akapf in component akap_sig (uM)'); LEGEND_CONSTANTS(:,396) = strpad('RyR_arn in component akap_sig (uM)'); LEGEND_CONSTANTS(:,397) = strpad('RyR_arp in component akap_sig (uM)'); LEGEND_CONSTANTS(:,31) = strpad('RyR_tot in component akap_sig (uM)'); LEGEND_CONSTANTS(:,294) = strpad('RyRf in component akap_sig (uM)'); LEGEND_STATES(:,6) = strpad('RyRp in component akap_sig (uM)'); LEGEND_CONSTANTS(:,274) = strpad('akap_sig_ICaLf_sum in component akap_sig (dimensionless)'); LEGEND_ALGEBRAIC(:,4) = strpad('akap_sig_ICaLp_dif in component akap_sig (uM)'); LEGEND_CONSTANTS(:,388) = strpad('akap_sig_PKAf_arg in component akap_sig (dimensionless)'); LEGEND_CONSTANTS(:,382) = strpad('akap_sig_PKAf_b in component akap_sig (uM)'); LEGEND_CONSTANTS(:,383) = strpad('akap_sig_PKAf_c in component akap_sig (mM2_times_1e_minus_6)'); LEGEND_CONSTANTS(:,384) = strpad('akap_sig_PKAf_d in component akap_sig (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,389) = strpad('akap_sig_PKAf_mag in component akap_sig (uM)'); LEGEND_CONSTANTS(:,385) = strpad('akap_sig_PKAf_rr in component akap_sig (mol6_per_m18_times_1e_minus_18)'); LEGEND_CONSTANTS(:,390) = strpad('akap_sig_PKAf_x in component akap_sig (dimensionless)'); LEGEND_CONSTANTS(:,386) = strpad('akap_sig_PKAf_yi in component akap_sig (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,387) = strpad('akap_sig_PKAf_yr in component akap_sig (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,289) = strpad('akap_sig_PP1f_cav_arg in component akap_sig (dimensionless)'); LEGEND_CONSTANTS(:,283) = strpad('akap_sig_PP1f_cav_b in component akap_sig (uM)'); LEGEND_CONSTANTS(:,284) = strpad('akap_sig_PP1f_cav_c in component akap_sig (mM2_times_1e_minus_6)'); LEGEND_CONSTANTS(:,285) = strpad('akap_sig_PP1f_cav_d in component akap_sig (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,290) = strpad('akap_sig_PP1f_cav_mag in component akap_sig (uM)'); LEGEND_CONSTANTS(:,286) = strpad('akap_sig_PP1f_cav_rr in component akap_sig (mol6_per_m18_times_1e_minus_18)'); LEGEND_CONSTANTS(:,291) = strpad('akap_sig_PP1f_cav_x in component akap_sig (dimensionless)'); LEGEND_CONSTANTS(:,287) = strpad('akap_sig_PP1f_cav_yi in component akap_sig (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,288) = strpad('akap_sig_PP1f_cav_yr in component akap_sig (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,293) = strpad('akap_sig_RyRf_sum in component akap_sig (dimensionless)'); LEGEND_ALGEBRAIC(:,5) = strpad('akap_sig_RyRp_dif in component akap_sig (uM)'); LEGEND_ALGEBRAIC(:,11) = strpad('fp_ICaL in component akap_sig (dimensionless)'); LEGEND_ALGEBRAIC(:,39) = strpad('fp_RyR in component akap_sig (dimensionless)'); LEGEND_CONSTANTS(:,32) = strpad('ka_ical in component akap_sig (hertz)'); LEGEND_CONSTANTS(:,33) = strpad('ka_ryr in component akap_sig (hertz)'); LEGEND_CONSTANTS(:,34) = strpad('kp_ical in component akap_sig (hertz)'); LEGEND_CONSTANTS(:,35) = strpad('kp_ryr in component akap_sig (hertz)'); LEGEND_STATES(:,7) = strpad('C in component pka_cav (uM)'); LEGEND_VOI = strpad('time in component engine (ms)'); LEGEND_CONSTANTS(:,36) = strpad('Gi_tot in component beta (uM)'); LEGEND_CONSTANTS(:,296) = strpad('Gs_tot in component beta (uM)'); LEGEND_CONSTANTS(:,298) = strpad('R_b2_tot in component beta (uM)'); LEGEND_CONSTANTS(:,37) = strpad('f_Gi_cav in component beta (dimensionless)'); LEGEND_CONSTANTS(:,216) = strpad('f_Gi_eca in component beta (dimensionless)'); LEGEND_CONSTANTS(:,38) = strpad('f_Gs_cav in component beta (dimensionless)'); LEGEND_CONSTANTS(:,217) = strpad('f_Gs_cyt in component beta (dimensionless)'); LEGEND_CONSTANTS(:,39) = strpad('f_Gs_eca in component beta (dimensionless)'); LEGEND_CONSTANTS(:,40) = strpad('f_Rb1_cav in component beta (dimensionless)'); LEGEND_CONSTANTS(:,299) = strpad('f_Rb1_cyt in component beta (dimensionless)'); LEGEND_CONSTANTS(:,41) = strpad('f_Rb1_eca in component beta (dimensionless)'); LEGEND_CONSTANTS(:,42) = strpad('f_Rb2_cav in component beta (dimensionless)'); LEGEND_CONSTANTS(:,300) = strpad('f_Rb2_eca in component beta (dimensionless)'); LEGEND_CONSTANTS(:,43) = strpad('k_act1_Gi in component beta (hertz)'); LEGEND_CONSTANTS(:,44) = strpad('k_act1_Gs in component beta (hertz)'); LEGEND_CONSTANTS(:,45) = strpad('k_act2_Gi in component beta (hertz)'); LEGEND_CONSTANTS(:,46) = strpad('k_act2_Gs in component beta (hertz)'); LEGEND_CONSTANTS(:,47) = strpad('k_b1_c in component beta (uM)'); LEGEND_CONSTANTS(:,48) = strpad('k_b1_h in component beta (uM)'); LEGEND_CONSTANTS(:,49) = strpad('k_b1_l in component beta (uM)'); LEGEND_CONSTANTS(:,50) = strpad('k_b2_a in component beta (uM)'); LEGEND_CONSTANTS(:,51) = strpad('k_b2_c in component beta (uM)'); LEGEND_CONSTANTS(:,52) = strpad('k_b2_f in component beta (uM)'); LEGEND_CONSTANTS(:,53) = strpad('k_b2_h in component beta (uM)'); LEGEND_CONSTANTS(:,54) = strpad('k_b2_l in component beta (uM)'); LEGEND_CONSTANTS(:,55) = strpad('k_b2_n in component beta (uM)'); LEGEND_CONSTANTS(:,218) = strpad('k_grk_dp in component beta (hertz)'); LEGEND_CONSTANTS(:,219) = strpad('k_grk_p in component beta (hertz)'); LEGEND_CONSTANTS(:,220) = strpad('k_hydr_Gi in component beta (hertz)'); LEGEND_CONSTANTS(:,56) = strpad('k_hydr_Gs in component beta (hertz)'); LEGEND_CONSTANTS(:,258) = strpad('k_pka_dp in component beta (hertz)'); LEGEND_CONSTANTS(:,221) = strpad('k_pka_p in component beta (per_mM_per_ms)'); LEGEND_CONSTANTS(:,222) = strpad('k_reas_Gi in component beta (per_mM_per_ms)'); LEGEND_CONSTANTS(:,57) = strpad('k_reas_Gs in component beta (per_mM_per_ms)'); LEGEND_CONSTANTS(:,58) = strpad('rate_bds in component beta (per_mM_per_ms)'); LEGEND_CONSTANTS(:,59) = strpad('GRK in component beta_cav (dimensionless)'); LEGEND_STATES(:,8) = strpad('Gi_aGDP in component beta_cav (uM)'); LEGEND_STATES(:,9) = strpad('Gi_aGTP in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,6) = strpad('Gi_abg in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,87) = strpad('Gi_f in component beta_cav (uM)'); LEGEND_STATES(:,10) = strpad('Gs_aGDP in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,7) = strpad('Gs_abg in component beta_cav (uM)'); LEGEND_STATES(:,11) = strpad('Gs_bg in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,157) = strpad('Gs_f in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,189) = strpad('LRGs_tot in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,165) = strpad('LRb1 in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,169) = strpad('LRb1Gs in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,181) = strpad('LRb2 in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,109) = strpad('LRb2Gi in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,185) = strpad('LRb2Gs in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,197) = strpad('RGs_tot in component beta_cav (uM)'); LEGEND_CONSTANTS(:,323) = strpad('R_b1_tot in component beta_cav (uM)'); LEGEND_CONSTANTS(:,324) = strpad('R_b2_tot in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,173) = strpad('Rb1Gs in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,161) = strpad('Rb1_f in component beta_cav (uM)'); LEGEND_STATES(:,12) = strpad('Rb1_grk_tot in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,34) = strpad('Rb1_np_tot in component beta_cav (uM)'); LEGEND_STATES(:,13) = strpad('Rb1_pka_tot in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,99) = strpad('Rb2Gi in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,193) = strpad('Rb2Gs in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,177) = strpad('Rb2_f in component beta_cav (uM)'); LEGEND_STATES(:,14) = strpad('Rb2_grk_tot in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,58) = strpad('Rb2_np_tot in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,71) = strpad('Rb2_pka_f in component beta_cav (uM)'); LEGEND_STATES(:,15) = strpad('Rb2_pka_tot in component beta_cav (uM)'); LEGEND_CONSTANTS(:,224) = strpad('beta_cav_Gs_f_a in component beta_cav (mol4_per_m12_times_1e_minus_12)'); LEGEND_ALGEBRAIC(:,132) = strpad('beta_cav_Gs_f_arg in component beta_cav (dimensionless)'); LEGEND_ALGEBRAIC(:,70) = strpad('beta_cav_Gs_f_b in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,86) = strpad('beta_cav_Gs_f_c in component beta_cav (mM2_times_1e_minus_6)'); LEGEND_CONSTANTS(:,259) = strpad('beta_cav_Gs_f_c11 in component beta_cav (mol5_per_m15_times_1e_minus_15)'); LEGEND_CONSTANTS(:,272) = strpad('beta_cav_Gs_f_c22 in component beta_cav (mol5_per_m15_times_1e_minus_15)'); LEGEND_CONSTANTS(:,280) = strpad('beta_cav_Gs_f_c33 in component beta_cav (mol6_per_m18_times_1e_minus_18)'); LEGEND_ALGEBRAIC(:,98) = strpad('beta_cav_Gs_f_d in component beta_cav (mM3_times_1e_minus_9)'); LEGEND_ALGEBRAIC(:,148) = strpad('beta_cav_Gs_f_i in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,138) = strpad('beta_cav_Gs_f_mag in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,152) = strpad('beta_cav_Gs_f_r in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,110) = strpad('beta_cav_Gs_f_rr in component beta_cav (mol6_per_m18_times_1e_minus_18)'); LEGEND_ALGEBRAIC(:,143) = strpad('beta_cav_Gs_f_x in component beta_cav (dimensionless)'); LEGEND_ALGEBRAIC(:,118) = strpad('beta_cav_Gs_f_yi in component beta_cav (mM3_times_1e_minus_9)'); LEGEND_ALGEBRAIC(:,126) = strpad('beta_cav_Gs_f_yr in component beta_cav (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,223) = strpad('beta_cav_Rb2_pka_f_a in component beta_cav (uM)'); LEGEND_ALGEBRAIC(:,35) = strpad('beta_cav_Rb2_pka_f_b in component beta_cav (mM2_times_1e_minus_6)'); LEGEND_ALGEBRAIC(:,59) = strpad('beta_cav_Rb2_pka_f_c in component beta_cav (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,60) = strpad('L in component iso (uM)'); LEGEND_CONSTANTS(:,61) = strpad('k_GsAct_b2 in component beta_cav (dimensionless)'); LEGEND_CONSTANTS(:,62) = strpad('GRK in component beta_cyt (dimensionless)'); LEGEND_STATES(:,16) = strpad('Gs_aGDP in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,8) = strpad('Gs_abg in component beta_cyt (uM)'); LEGEND_STATES(:,17) = strpad('Gs_bg in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,100) = strpad('Gs_f in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,111) = strpad('LRb1Gs_np in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,119) = strpad('LRb1_np in component beta_cyt (uM)'); LEGEND_CONSTANTS(:,326) = strpad('R_b1_tot in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,120) = strpad('Rb1Gs_np in component beta_cyt (uM)'); LEGEND_STATES(:,18) = strpad('Rb1_grk_tot in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,88) = strpad('Rb1_np_f in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,36) = strpad('Rb1_np_tot in component beta_cyt (uM)'); LEGEND_STATES(:,19) = strpad('Rb1_pka_tot in component beta_cyt (uM)'); LEGEND_CONSTANTS(:,225) = strpad('beta_cyt_Rb1_np_f_a in component beta_cyt (uM)'); LEGEND_ALGEBRAIC(:,60) = strpad('beta_cyt_Rb1_np_f_b in component beta_cyt (mM2_times_1e_minus_6)'); LEGEND_ALGEBRAIC(:,72) = strpad('beta_cyt_Rb1_np_f_c in component beta_cyt (mM3_times_1e_minus_9)'); LEGEND_STATES(:,20) = strpad('C in component pka_cyt (uM)'); LEGEND_CONSTANTS(:,63) = strpad('GRK in component beta_eca (dimensionless)'); LEGEND_STATES(:,21) = strpad('Gi_aGDP in component beta_eca (uM)'); LEGEND_STATES(:,22) = strpad('Gi_aGTP in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,9) = strpad('Gi_abg in component beta_eca (uM)'); LEGEND_STATES(:,23) = strpad('Gi_bg in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,90) = strpad('Gi_f in component beta_eca (uM)'); LEGEND_STATES(:,24) = strpad('Gs_aGDP in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,10) = strpad('Gs_abg in component beta_eca (uM)'); LEGEND_STATES(:,25) = strpad('Gs_bg in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,158) = strpad('Gs_f in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,190) = strpad('LRGs_tot in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,166) = strpad('LRb1 in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,170) = strpad('LRb1Gs in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,182) = strpad('LRb2 in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,112) = strpad('LRb2Gi in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,186) = strpad('LRb2Gs in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,198) = strpad('RGs_tot in component beta_eca (uM)'); LEGEND_CONSTANTS(:,331) = strpad('R_b1_tot in component beta_eca (uM)'); LEGEND_CONSTANTS(:,332) = strpad('R_b2_tot in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,174) = strpad('Rb1Gs in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,162) = strpad('Rb1_f in component beta_eca (uM)'); LEGEND_STATES(:,26) = strpad('Rb1_grk_tot in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,37) = strpad('Rb1_np_tot in component beta_eca (uM)'); LEGEND_STATES(:,27) = strpad('Rb1_pka_tot in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,102) = strpad('Rb2Gi in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,194) = strpad('Rb2Gs in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,178) = strpad('Rb2_f in component beta_eca (uM)'); LEGEND_STATES(:,28) = strpad('Rb2_grk_tot in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,61) = strpad('Rb2_np_tot in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,74) = strpad('Rb2_pka_f in component beta_eca (uM)'); LEGEND_STATES(:,29) = strpad('Rb2_pka_tot in component beta_eca (uM)'); LEGEND_CONSTANTS(:,227) = strpad('beta_eca_Gs_f_a in component beta_eca (mol4_per_m12_times_1e_minus_12)'); LEGEND_ALGEBRAIC(:,133) = strpad('beta_eca_Gs_f_arg in component beta_eca (dimensionless)'); LEGEND_ALGEBRAIC(:,73) = strpad('beta_eca_Gs_f_b in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,89) = strpad('beta_eca_Gs_f_c in component beta_eca (mM2_times_1e_minus_6)'); LEGEND_CONSTANTS(:,260) = strpad('beta_eca_Gs_f_c11 in component beta_eca (mol5_per_m15_times_1e_minus_15)'); LEGEND_CONSTANTS(:,273) = strpad('beta_eca_Gs_f_c22 in component beta_eca (mol5_per_m15_times_1e_minus_15)'); LEGEND_CONSTANTS(:,281) = strpad('beta_eca_Gs_f_c33 in component beta_eca (mol6_per_m18_times_1e_minus_18)'); LEGEND_ALGEBRAIC(:,101) = strpad('beta_eca_Gs_f_d in component beta_eca (mM3_times_1e_minus_9)'); LEGEND_ALGEBRAIC(:,149) = strpad('beta_eca_Gs_f_i in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,139) = strpad('beta_eca_Gs_f_mag in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,153) = strpad('beta_eca_Gs_f_r in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,113) = strpad('beta_eca_Gs_f_rr in component beta_eca (mol6_per_m18_times_1e_minus_18)'); LEGEND_ALGEBRAIC(:,144) = strpad('beta_eca_Gs_f_x in component beta_eca (dimensionless)'); LEGEND_ALGEBRAIC(:,121) = strpad('beta_eca_Gs_f_yi in component beta_eca (mM3_times_1e_minus_9)'); LEGEND_ALGEBRAIC(:,127) = strpad('beta_eca_Gs_f_yr in component beta_eca (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,226) = strpad('beta_eca_Rb2_pka_f_a in component beta_eca (uM)'); LEGEND_ALGEBRAIC(:,38) = strpad('beta_eca_Rb2_pka_f_b in component beta_eca (mM2_times_1e_minus_6)'); LEGEND_ALGEBRAIC(:,62) = strpad('beta_eca_Rb2_pka_f_c in component beta_eca (mM3_times_1e_minus_9)'); LEGEND_CONSTANTS(:,64) = strpad('k_GsAct_b2 in component beta_eca (dimensionless)'); LEGEND_STATES(:,30) = strpad('C in component pka_eca (uM)'); LEGEND_CONSTANTS(:,303) = strpad('AF in component cell (m2_mol_per_s_per_A_times_1e_minus_4)'); LEGEND_ALGEBRAIC(:,187) = strpad('Ca in component calcium (mM)'); LEGEND_ALGEBRAIC(:,103) = strpad('Ca_CaL in component calcium (mM)'); LEGEND_ALGEBRAIC(:,128) = strpad('Ca_jsr in component calcium (mM)'); LEGEND_STATES(:,31) = strpad('Ca_nsr in component calcium (mM)'); LEGEND_ALGEBRAIC(:,150) = strpad('Ca_sr in component calcium (mM)'); LEGEND_ALGEBRAIC(:,306) = strpad('ICa_tot in component calcium (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,218) = strpad('ICab in component icab (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,298) = strpad('INaCa in component inaca (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,292) = strpad('INaCaSR in component inaca (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,202) = strpad('Idiff_Ca in component diff (mM_per_ms)'); LEGEND_ALGEBRAIC(:,210) = strpad('Idiff_sr in component diff (mM_per_ms)'); LEGEND_ALGEBRAIC(:,305) = strpad('IpCa in component ipca (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,320) = strpad('Irel in component irel (mM_per_ms)'); LEGEND_ALGEBRAIC(:,212) = strpad('Itr in component diff (mM_per_ms)'); LEGEND_ALGEBRAIC(:,329) = strpad('Iup in component iup (mM_per_ms)'); LEGEND_CONSTANTS(:,65) = strpad('Ka_tni in component calcium (uM)'); LEGEND_CONSTANTS(:,66) = strpad('Kp_tni in component calcium (uM)'); LEGEND_CONSTANTS(:,67) = strpad('PP2A in component pp1 (uM)'); LEGEND_CONSTANTS(:,229) = strpad('bar_sum in component calcium (mM)'); LEGEND_CONSTANTS(:,68) = strpad('bsl_bar in component calcium (mM)'); LEGEND_CONSTANTS(:,69) = strpad('bsl_km in component calcium (mM)'); LEGEND_CONSTANTS(:,70) = strpad('bsr_bar in component calcium (mM)'); LEGEND_CONSTANTS(:,71) = strpad('bsr_km in component calcium (mM)'); LEGEND_ALGEBRAIC(:,63) = strpad('calcium_Ca_CaL_b in component calcium (mM)'); LEGEND_ALGEBRAIC(:,75) = strpad('calcium_Ca_CaL_c in component calcium (mM2)'); LEGEND_ALGEBRAIC(:,91) = strpad('calcium_Ca_CaL_d in component calcium (mM3)'); LEGEND_ALGEBRAIC(:,179) = strpad('calcium_Ca_b in component calcium (mM)'); LEGEND_ALGEBRAIC(:,183) = strpad('calcium_Ca_c in component calcium (mM2)'); LEGEND_ALGEBRAIC(:,171) = strpad('calcium_Ca_d in component calcium (mM3)'); LEGEND_ALGEBRAIC(:,114) = strpad('calcium_Ca_jsr_b in component calcium (mM)'); LEGEND_ALGEBRAIC(:,122) = strpad('calcium_Ca_jsr_c in component calcium (mM2)'); LEGEND_CONSTANTS(:,313) = strpad('calcium_Ca_nsr_r1 in component calcium (dimensionless)'); LEGEND_ALGEBRAIC(:,134) = strpad('calcium_Ca_sr_b in component calcium (mM)'); LEGEND_ALGEBRAIC(:,140) = strpad('calcium_Ca_sr_c in component calcium (mM2)'); LEGEND_ALGEBRAIC(:,145) = strpad('calcium_Ca_sr_d in component calcium (mM3)'); LEGEND_ALGEBRAIC(:,154) = strpad('calcium_fhat_val in component calcium (dimensionless)'); LEGEND_CONSTANTS(:,306) = strpad('calcium_uCa_CaL_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_CONSTANTS(:,316) = strpad('calcium_uCa_CaL_r2 in component calcium (dimensionless)'); LEGEND_CONSTANTS(:,309) = strpad('calcium_uCa_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_CONSTANTS(:,314) = strpad('calcium_uCa_r2 in component calcium (dimensionless)'); LEGEND_CONSTANTS(:,317) = strpad('calcium_uCa_sr_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_CONSTANTS(:,318) = strpad('calcium_uCa_sr_r2 in component calcium (dimensionless)'); LEGEND_CONSTANTS(:,72) = strpad('cbar in component calcium (mM)'); LEGEND_CONSTANTS(:,73) = strpad('csqn_bar in component calcium (mM)'); LEGEND_CONSTANTS(:,74) = strpad('csqn_km in component calcium (mM)'); LEGEND_STATES(:,32) = strpad('f_tni in component calcium (dimensionless)'); LEGEND_ALGEBRAIC(:,159) = strpad('fhat in component calcium (dimensionless)'); LEGEND_ALGEBRAIC(:,249) = strpad('ICaL in component ical (uA_per_cm2)'); LEGEND_CONSTANTS(:,75) = strpad('ka_tni in component calcium (hertz)'); LEGEND_CONSTANTS(:,76) = strpad('kc in component calcium (mM)'); LEGEND_CONSTANTS(:,231) = strpad('km_pro in component calcium (mM2)'); LEGEND_CONSTANTS(:,262) = strpad('km_sum in component calcium (mM)'); LEGEND_CONSTANTS(:,77) = strpad('kp_tni in component calcium (hertz)'); LEGEND_ALGEBRAIC(:,167) = strpad('kpro in component calcium (mM2)'); LEGEND_ALGEBRAIC(:,175) = strpad('ksum in component calcium (mM)'); LEGEND_ALGEBRAIC(:,163) = strpad('kt in component calcium (mM)'); LEGEND_CONSTANTS(:,78) = strpad('ktn in component calcium (mM)'); LEGEND_CONSTANTS(:,230) = strpad('ktp in component calcium (mM)'); LEGEND_CONSTANTS(:,319) = strpad('r3 in component calcium (dimensionless)'); LEGEND_CONSTANTS(:,263) = strpad('ss_pro in component calcium (mM2)'); LEGEND_CONSTANTS(:,275) = strpad('ss_sum in component calcium (mM)'); LEGEND_CONSTANTS(:,79) = strpad('tbar in component calcium (mM)'); LEGEND_STATES(:,33) = strpad('uCa in component calcium (mM)'); LEGEND_STATES(:,34) = strpad('uCa_CaL in component calcium (mM)'); LEGEND_STATES(:,35) = strpad('uCa_jsr in component calcium (mM)'); LEGEND_STATES(:,36) = strpad('uCa_sr in component calcium (mM)'); LEGEND_CONSTANTS(:,305) = strpad('v_CaL in component cell (uL)'); LEGEND_CONSTANTS(:,308) = strpad('v_cyt in component cell (uL)'); LEGEND_CONSTANTS(:,311) = strpad('v_jsr in component cell (uL)'); LEGEND_CONSTANTS(:,312) = strpad('v_nsr in component cell (uL)'); LEGEND_CONSTANTS(:,315) = strpad('v_sr in component cell (uL)'); LEGEND_CONSTANTS(:,80) = strpad('CaMK0 in component camk (dimensionless)'); LEGEND_CONSTANTS(:,81) = strpad('K in component camk (dimensionless)'); LEGEND_CONSTANTS(:,82) = strpad('Km in component camk (mM)'); LEGEND_CONSTANTS(:,83) = strpad('PP1_eca in component pp1 (uM)'); LEGEND_ALGEBRAIC(:,207) = strpad('PP1_tot in component camk (uM)'); LEGEND_ALGEBRAIC(:,204) = strpad('PP1f_cyt in component pp1 (uM)'); LEGEND_ALGEBRAIC(:,195) = strpad('active in component camk (dimensionless)'); LEGEND_ALGEBRAIC(:,191) = strpad('bound in component camk (dimensionless)'); LEGEND_ALGEBRAIC(:,200) = strpad('c in component camk (dimensionless)'); LEGEND_ALGEBRAIC(:,201) = strpad('camk_f_ryr_d in component camk (dimensionless)'); LEGEND_CONSTANTS(:,84) = strpad('camk_trap_alpha in component camk (mS_per_uF)'); LEGEND_CONSTANTS(:,85) = strpad('camk_trap_beta in component camk (mS_per_uF)'); LEGEND_STATES(:,37) = strpad('f_ical in component camk (dimensionless)'); LEGEND_STATES(:,38) = strpad('f_ik1 in component camk (dimensionless)'); LEGEND_STATES(:,39) = strpad('f_ina in component camk (dimensionless)'); LEGEND_STATES(:,40) = strpad('f_ito in component camk (dimensionless)'); LEGEND_STATES(:,41) = strpad('f_plb in component camk (dimensionless)'); LEGEND_STATES(:,42) = strpad('f_ryr in component camk (dimensionless)'); LEGEND_CONSTANTS(:,232) = strpad('tau_cal in component camk (ms)'); LEGEND_CONSTANTS(:,233) = strpad('tau_ik1 in component camk (ms)'); LEGEND_CONSTANTS(:,234) = strpad('tau_ina in component camk (ms)'); LEGEND_CONSTANTS(:,235) = strpad('tau_ito in component camk (ms)'); LEGEND_CONSTANTS(:,86) = strpad('tau_plb in component camk (ms)'); LEGEND_CONSTANTS(:,87) = strpad('tau_ryr in component camk (ms)'); LEGEND_STATES(:,43) = strpad('trap in component camk (dimensionless)'); LEGEND_STATES(:,44) = strpad('cAMP_cav in component camp (uM)'); LEGEND_STATES(:,45) = strpad('cAMP_cyt in component camp (uM)'); LEGEND_STATES(:,46) = strpad('cAMP_eca in component camp (uM)'); LEGEND_ALGEBRAIC(:,96) = strpad('camp_cAMP_cav_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,106) = strpad('camp_cAMP_cav_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,136) = strpad('camp_cAMP_cav_pde in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,116) = strpad('camp_cAMP_cyt_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,124) = strpad('camp_cAMP_cyt_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,147) = strpad('camp_cAMP_cyt_pde in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,84) = strpad('camp_cAMP_eca_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,97) = strpad('camp_cAMP_eca_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,125) = strpad('camp_cAMP_eca_pde in component camp (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,115) = strpad('dcAMP_PDE2_cav in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,130) = strpad('dcAMP_PDE2_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,108) = strpad('dcAMP_PDE2_eca in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,123) = strpad('dcAMP_PDE3_cav in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,137) = strpad('dcAMP_PDE3_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,129) = strpad('dcAMP_PDE4_cav in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,142) = strpad('dcAMP_PDE4_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,117) = strpad('dcAMP_PDE4_eca in component pde (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_CONSTANTS(:,236) = strpad('j_cav_cyt in component camp (m3_per_s_times_1e_minus_9)'); LEGEND_CONSTANTS(:,237) = strpad('j_cav_eca in component camp (m3_per_s_times_1e_minus_9)'); LEGEND_CONSTANTS(:,238) = strpad('j_eca_cyt in component camp (m3_per_s_times_1e_minus_9)'); LEGEND_ALGEBRAIC(:,146) = strpad('dcAMP in component pka_cav (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,155) = strpad('dcAMP in component pka_cyt (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_ALGEBRAIC(:,135) = strpad('dcAMP in component pka_eca (mol_per_m3_per_s_times_1e_minus_3)'); LEGEND_CONSTANTS(:,307) = strpad('v_cav in component cell (uL)'); LEGEND_CONSTANTS(:,310) = strpad('v_eca in component cell (uL)'); LEGEND_CONSTANTS(:,88) = strpad('F in component phys (C_per_mol)'); LEGEND_CONSTANTS(:,302) = strpad('capArea in component cell (cm2)'); LEGEND_CONSTANTS(:,301) = strpad('geoArea in component cell (cm2)'); LEGEND_CONSTANTS(:,89) = strpad('length in component cell (cm)'); LEGEND_CONSTANTS(:,90) = strpad('pi in component cell (dimensionless)'); LEGEND_CONSTANTS(:,91) = strpad('radius in component cell (cm)'); LEGEND_CONSTANTS(:,304) = strpad('volume in component cell (uL)'); LEGEND_ALGEBRAIC(:,334) = strpad('CTKCl in component ctkcl (mM_per_ms)'); LEGEND_ALGEBRAIC(:,355) = strpad('CTNaCl in component ctnacl (mM_per_ms)'); LEGEND_STATES(:,47) = strpad('Cl in component chloride (mM)'); LEGEND_STATES(:,48) = strpad('Cl_sr in component chloride (mM)'); LEGEND_ALGEBRAIC(:,324) = strpad('IClCa in component iclca (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,331) = strpad('ICl_tot in component chloride (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,330) = strpad('IClb in component iclb (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,205) = strpad('Idiff_Cl in component diff (mM_per_ms)'); LEGEND_CONSTANTS(:,333) = strpad('chloride_Cl_r1 in component chloride (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_CONSTANTS(:,338) = strpad('chloride_Cl_r2 in component chloride (dimensionless)'); LEGEND_CONSTANTS(:,342) = strpad('chloride_Cl_sr_r1 in component chloride (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_ALGEBRAIC(:,328) = strpad('ECl in component nernst (mV)'); LEGEND_ALGEBRAIC(:,332) = strpad('EK in component nernst (mV)'); LEGEND_CONSTANTS(:,92) = strpad('KClBar in component ctkcl (mM_per_ms)'); LEGEND_ALGEBRAIC(:,333) = strpad('ctkcl_CTKCl_z1 in component ctkcl (mV)'); LEGEND_CONSTANTS(:,93) = strpad('ctkcl_CTKCl_z2 in component ctkcl (mV)'); LEGEND_ALGEBRAIC(:,353) = strpad('ENa in component nernst (mV)'); LEGEND_CONSTANTS(:,94) = strpad('NaClBar in component ctnacl (mM_per_ms)'); LEGEND_ALGEBRAIC(:,354) = strpad('ctnacl_CTNaCl_z1 in component ctnacl (g4_m8_per_s12_per_A4)'); LEGEND_CONSTANTS(:,347) = strpad('ctnacl_CTNaCl_z2 in component ctnacl (g4_m8_per_s12_per_A4)'); LEGEND_ALGEBRAIC(:,208) = strpad('Idiff_Na in component diff (mM_per_ms)'); LEGEND_STATES(:,49) = strpad('Na in component sodium (mM)'); LEGEND_STATES(:,50) = strpad('Na_sr in component sodium (mM)'); LEGEND_CONSTANTS(:,95) = strpad('tau in component diff (ms)'); LEGEND_CONSTANTS(:,96) = strpad('tau_sr in component diff (ms)'); LEGEND_CONSTANTS(:,97) = strpad('tau_tr in component diff (ms)'); LEGEND_CONSTANTS(:,98) = strpad('Cao in component extra (mM)'); LEGEND_CONSTANTS(:,99) = strpad('Clo in component extra (mM)'); LEGEND_CONSTANTS(:,100) = strpad('Ko in component extra (mM)'); LEGEND_CONSTANTS(:,101) = strpad('Nao in component extra (mM)'); LEGEND_CONSTANTS(:,264) = strpad('FRT in component phys (per_mV)'); LEGEND_STATES(:,51) = strpad('V in component membrane (mV)'); LEGEND_ALGEBRAIC(:,216) = strpad('efrt in component icab (dimensionless)'); LEGEND_CONSTANTS(:,102) = strpad('pCab in component icab (cm_per_s)'); LEGEND_ALGEBRAIC(:,214) = strpad('vfrt in component icab (dimensionless)'); LEGEND_ALGEBRAIC(:,224) = strpad('f_hat in component ical (dimensionless)'); LEGEND_ALGEBRAIC(:,236) = strpad('ICaL in component ical_camk (uA_per_cm2)'); LEGEND_CONSTANTS(:,398) = strpad('ical_f_hat_ratio in component ical (dimensionless)'); LEGEND_ALGEBRAIC(:,221) = strpad('ical_f_hat_val in component ical (dimensionless)'); LEGEND_ALGEBRAIC(:,247) = strpad('ICaL in component ical_np (uA_per_cm2)'); LEGEND_STATES(:,52) = strpad('C in component ical_camk (dimensionless)'); LEGEND_STATES(:,53) = strpad('CI in component ical_camk (dimensionless)'); LEGEND_STATES(:,54) = strpad('CIs in component ical_camk (dimensionless)'); LEGEND_STATES(:,55) = strpad('Cs in component ical_camk (dimensionless)'); LEGEND_CONSTANTS(:,276) = strpad('FFRT in component phys (s4_A2_per_g_per_m2_per_mol)'); LEGEND_ALGEBRAIC(:,233) = strpad('IBar in component ical_camk (uA_per_cm2)'); LEGEND_STATES(:,56) = strpad('O in component ical_camk (dimensionless)'); LEGEND_STATES(:,57) = strpad('OI in component ical_camk (dimensionless)'); LEGEND_STATES(:,58) = strpad('OIs in component ical_camk (dimensionless)'); LEGEND_STATES(:,59) = strpad('Os in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,227) = strpad('PCa in component ical_camk (cm_per_s)'); LEGEND_ALGEBRAIC(:,156) = strpad('ac_inf in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,196) = strpad('ac_tau in component ical_np (ms)'); LEGEND_ALGEBRAIC(:,203) = strpad('alpha in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,206) = strpad('beta in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,184) = strpad('delta in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,243) = strpad('delta1 in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,160) = strpad('delta_tau in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,230) = strpad('ical_camk_IBar_vv in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,231) = strpad('ical_camk_delta1_xs_cor in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,237) = strpad('ical_camk_delta1_y_cor in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,215) = strpad('in_a in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,217) = strpad('in_b in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,168) = strpad('in_hi_inf in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,219) = strpad('in_hi_tau in component ical_camk (ms)'); LEGEND_ALGEBRAIC(:,164) = strpad('in_inf in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,172) = strpad('in_lo_inf in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,222) = strpad('in_lo_tau in component ical_camk (ms)'); LEGEND_ALGEBRAIC(:,188) = strpad('inca in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,176) = strpad('ss_cal_10 in component ical_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,180) = strpad('ss_cal_4 in component ical_camk (dimensionless)'); LEGEND_CONSTANTS(:,103) = strpad('theta in component ical_camk (mS_per_uF)'); LEGEND_CONSTANTS(:,104) = strpad('theta1 in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,225) = strpad('x in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,228) = strpad('xs in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,234) = strpad('y in component ical_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,240) = strpad('ys in component ical_camk (mS_per_uF)'); LEGEND_STATES(:,60) = strpad('C in component ical_np (dimensionless)'); LEGEND_STATES(:,61) = strpad('CI in component ical_np (dimensionless)'); LEGEND_STATES(:,62) = strpad('CIs in component ical_np (dimensionless)'); LEGEND_STATES(:,63) = strpad('Cs in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,245) = strpad('IBar in component ical_np (uA_per_cm2)'); LEGEND_STATES(:,64) = strpad('O in component ical_np (dimensionless)'); LEGEND_STATES(:,65) = strpad('OI in component ical_np (dimensionless)'); LEGEND_STATES(:,66) = strpad('OIs in component ical_np (dimensionless)'); LEGEND_STATES(:,67) = strpad('Os in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,239) = strpad('PCa in component ical_np (cm_per_s)'); LEGEND_ALGEBRAIC(:,192) = strpad('ac_inf in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,209) = strpad('alpha in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,211) = strpad('beta in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,238) = strpad('delta in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,261) = strpad('delta1 in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,213) = strpad('delta_tau in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,242) = strpad('ical_np_IBar_vv in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,250) = strpad('ical_np_delta1_xs_cor in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,255) = strpad('ical_np_delta1_y_cor in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,223) = strpad('in_hi_inf in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,244) = strpad('in_hi_tau in component ical_np (ms)'); LEGEND_ALGEBRAIC(:,220) = strpad('in_inf in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,226) = strpad('in_lo_inf in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,229) = strpad('in_lo_tau in component ical_np (ms)'); LEGEND_ALGEBRAIC(:,241) = strpad('inca in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,232) = strpad('ss_cal_10 in component ical_np (dimensionless)'); LEGEND_ALGEBRAIC(:,235) = strpad('ss_cal_4 in component ical_np (dimensionless)'); LEGEND_CONSTANTS(:,105) = strpad('theta in component ical_np (mS_per_uF)'); LEGEND_CONSTANTS(:,106) = strpad('theta1 in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,246) = strpad('x in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,248) = strpad('xs in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,252) = strpad('y in component ical_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,258) = strpad('ys in component ical_np (mS_per_uF)'); LEGEND_CONSTANTS(:,107) = strpad('Gbar in component iclb (mS_per_cm2)'); LEGEND_ALGEBRAIC(:,254) = strpad('IClCa_bar in component iclca (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,318) = strpad('Irel_pure in component irel (mM_per_ms)'); LEGEND_ALGEBRAIC(:,321) = strpad('KClCa in component iclca (dimensionless)'); LEGEND_CONSTANTS(:,108) = strpad('PCl in component iclca (cm_per_s)'); LEGEND_STATES(:,68) = strpad('i2 in component iclca (dimensionless)'); LEGEND_ALGEBRAIC(:,12) = strpad('iclca_i2_alpha in component iclca (dimensionless)'); LEGEND_ALGEBRAIC(:,40) = strpad('iclca_i2_beta in component iclca (dimensionless)'); LEGEND_CONSTANTS(:,109) = strpad('kCaCl in component iclca (mM_per_ms)'); LEGEND_CONSTANTS(:,110) = strpad('tau in component iclca (ms)'); LEGEND_ALGEBRAIC(:,251) = strpad('vexp in component iclca (dimensionless)'); LEGEND_CONSTANTS(:,240) = strpad('Gbar in component ik1 (mS_per_cm2)'); LEGEND_ALGEBRAIC(:,340) = strpad('IK1 in component ik1 (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,339) = strpad('IK1_camk in component ik1 (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,338) = strpad('IK1_np in component ik1 (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,336) = strpad('ik1_IK1_np_alpha in component ik1 (dimensionless)'); LEGEND_ALGEBRAIC(:,337) = strpad('ik1_IK1_np_beta in component ik1 (dimensionless)'); LEGEND_ALGEBRAIC(:,335) = strpad('ik1_IK1_np_vv in component ik1 (mV)'); LEGEND_CONSTANTS(:,241) = strpad('GKr in component ikr (mS_per_cm2)'); LEGEND_ALGEBRAIC(:,341) = strpad('IKr in component ikr (uA_per_cm2)'); LEGEND_STATES(:,69) = strpad('ac in component ikr (dimensionless)'); LEGEND_ALGEBRAIC(:,13) = strpad('ikr_ac_tau in component ikr (ms)'); LEGEND_ALGEBRAIC(:,41) = strpad('inf in component ikr (dimensionless)'); LEGEND_ALGEBRAIC(:,257) = strpad('inx in component ikr (dimensionless)'); LEGEND_ALGEBRAIC(:,349) = strpad('EKs in component nernst (mV)'); LEGEND_ALGEBRAIC(:,260) = strpad('G in component iks (mS_per_cm2)'); LEGEND_ALGEBRAIC(:,352) = strpad('IKs in component iks (uA_per_cm2)'); LEGEND_CONSTANTS(:,404) = strpad('IKs_arn in component iks_sig (uM)'); LEGEND_ALGEBRAIC(:,350) = strpad('IKs_np in component iks (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,351) = strpad('IKs_pka in component iks (uA_per_cm2)'); LEGEND_CONSTANTS(:,111) = strpad('IKs_tot in component iks_sig (uM)'); LEGEND_ALGEBRAIC(:,267) = strpad('f_hat in component iks (dimensionless)'); LEGEND_ALGEBRAIC(:,263) = strpad('fp_iks in component iks_sig (dimensionless)'); LEGEND_CONSTANTS(:,405) = strpad('iks_f_hat_ratio in component iks (dimensionless)'); LEGEND_ALGEBRAIC(:,265) = strpad('iks_f_hat_val in component iks (dimensionless)'); LEGEND_STATES(:,70) = strpad('O1 in component iks_np (dimensionless)'); LEGEND_STATES(:,71) = strpad('O2 in component iks_np (dimensionless)'); LEGEND_STATES(:,72) = strpad('O1 in component iks_pka (dimensionless)'); LEGEND_STATES(:,73) = strpad('O2 in component iks_pka (dimensionless)'); LEGEND_STATES(:,74) = strpad('C1 in component iks_np (dimensionless)'); LEGEND_STATES(:,75) = strpad('C10 in component iks_np (dimensionless)'); LEGEND_STATES(:,76) = strpad('C11 in component iks_np (dimensionless)'); LEGEND_STATES(:,77) = strpad('C12 in component iks_np (dimensionless)'); LEGEND_STATES(:,78) = strpad('C13 in component iks_np (dimensionless)'); LEGEND_STATES(:,79) = strpad('C14 in component iks_np (dimensionless)'); LEGEND_STATES(:,80) = strpad('C15 in component iks_np (dimensionless)'); LEGEND_STATES(:,81) = strpad('C2 in component iks_np (dimensionless)'); LEGEND_STATES(:,82) = strpad('C3 in component iks_np (dimensionless)'); LEGEND_STATES(:,83) = strpad('C4 in component iks_np (dimensionless)'); LEGEND_STATES(:,84) = strpad('C5 in component iks_np (dimensionless)'); LEGEND_STATES(:,85) = strpad('C6 in component iks_np (dimensionless)'); LEGEND_STATES(:,86) = strpad('C7 in component iks_np (dimensionless)'); LEGEND_STATES(:,87) = strpad('C8 in component iks_np (dimensionless)'); LEGEND_STATES(:,88) = strpad('C9 in component iks_np (dimensionless)'); LEGEND_ALGEBRAIC(:,14) = strpad('a in component iks_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,42) = strpad('b in component iks_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,64) = strpad('d in component iks_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,76) = strpad('e in component iks_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,92) = strpad('g in component iks_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,93) = strpad('o in component iks_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,104) = strpad('p in component iks_np (mS_per_uF)'); LEGEND_CONSTANTS(:,112) = strpad('t in component iks_np (mS_per_uF)'); LEGEND_STATES(:,89) = strpad('C1 in component iks_pka (dimensionless)'); LEGEND_STATES(:,90) = strpad('C10 in component iks_pka (dimensionless)'); LEGEND_STATES(:,91) = strpad('C11 in component iks_pka (dimensionless)'); LEGEND_STATES(:,92) = strpad('C12 in component iks_pka (dimensionless)'); LEGEND_STATES(:,93) = strpad('C13 in component iks_pka (dimensionless)'); LEGEND_STATES(:,94) = strpad('C14 in component iks_pka (dimensionless)'); LEGEND_STATES(:,95) = strpad('C15 in component iks_pka (dimensionless)'); LEGEND_STATES(:,96) = strpad('C2 in component iks_pka (dimensionless)'); LEGEND_STATES(:,97) = strpad('C3 in component iks_pka (dimensionless)'); LEGEND_STATES(:,98) = strpad('C4 in component iks_pka (dimensionless)'); LEGEND_STATES(:,99) = strpad('C5 in component iks_pka (dimensionless)'); LEGEND_STATES(:,100) = strpad('C6 in component iks_pka (dimensionless)'); LEGEND_STATES(:,101) = strpad('C7 in component iks_pka (dimensionless)'); LEGEND_STATES(:,102) = strpad('C8 in component iks_pka (dimensionless)'); LEGEND_STATES(:,103) = strpad('C9 in component iks_pka (dimensionless)'); LEGEND_ALGEBRAIC(:,15) = strpad('a in component iks_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,43) = strpad('b in component iks_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,65) = strpad('d in component iks_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,77) = strpad('e in component iks_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,94) = strpad('g in component iks_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,95) = strpad('o in component iks_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,105) = strpad('p in component iks_pka (mS_per_uF)'); LEGEND_CONSTANTS(:,113) = strpad('t in component iks_pka (mS_per_uF)'); LEGEND_CONSTANTS(:,406) = strpad('IKs_arp in component iks_sig (uM)'); LEGEND_CONSTANTS(:,357) = strpad('IKsf in component iks_sig (uM)'); LEGEND_STATES(:,104) = strpad('IKsp in component iks_sig (uM)'); LEGEND_CONSTANTS(:,114) = strpad('K in component iks_sig (uM)'); LEGEND_CONSTANTS(:,115) = strpad('Ka_iks in component iks_sig (uM)'); LEGEND_CONSTANTS(:,116) = strpad('Kp_iks in component iks_sig (uM)'); LEGEND_CONSTANTS(:,117) = strpad('L in component iks_sig (uM)'); LEGEND_CONSTANTS(:,118) = strpad('M in component iks_sig (uM)'); LEGEND_CONSTANTS(:,400) = strpad('PKA_eca in component pka (uM)'); LEGEND_CONSTANTS(:,402) = strpad('PKAf in component iks_sig (uM)'); LEGEND_CONSTANTS(:,361) = strpad('PP1f_eca in component iks_sig (uM)'); LEGEND_CONSTANTS(:,119) = strpad('Yotiao in component iks_sig (uM)'); LEGEND_CONSTANTS(:,403) = strpad('Yotiaof in component iks_sig (uM)'); LEGEND_CONSTANTS(:,354) = strpad('iks_sig_IKsf_sum in component iks_sig (dimensionless)'); LEGEND_ALGEBRAIC(:,16) = strpad('iks_sig_IKsp_dif in component iks_sig (uM)'); LEGEND_CONSTANTS(:,401) = strpad('iks_sig_PKAf_sum in component iks_sig (dimensionless)'); LEGEND_CONSTANTS(:,359) = strpad('iks_sig_PP1f_eca_sum in component iks_sig (dimensionless)'); LEGEND_CONSTANTS(:,120) = strpad('ka_iks in component iks_sig (hertz)'); LEGEND_CONSTANTS(:,121) = strpad('kp_iks in component iks_sig (hertz)'); LEGEND_ALGEBRAIC(:,344) = strpad('IKur in component ikur (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,342) = strpad('IKur_np in component ikur (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,343) = strpad('IKur_p in component ikur (uA_per_cm2)'); LEGEND_CONSTANTS(:,122) = strpad('Ka_ikur in component ikur (uM)'); LEGEND_CONSTANTS(:,123) = strpad('Kp_ikur in component ikur (uM)'); LEGEND_STATES(:,105) = strpad('f_ikur in component ikur (dimensionless)'); LEGEND_ALGEBRAIC(:,271) = strpad('fhat in component ikur (dimensionless)'); LEGEND_CONSTANTS(:,124) = strpad('gbar_np in component ikur (mS_per_cm2)'); LEGEND_ALGEBRAIC(:,269) = strpad('ikur_fhat_val in component ikur (dimensionless)'); LEGEND_CONSTANTS(:,125) = strpad('ka_ikur in component ikur (hertz)'); LEGEND_CONSTANTS(:,126) = strpad('kp_ikur in component ikur (hertz)'); LEGEND_ALGEBRAIC(:,360) = strpad('INa in component ina (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,356) = strpad('INa_both in component ina (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,357) = strpad('INa_camk in component ina (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,358) = strpad('INa_np in component ina (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,359) = strpad('INa_pka in component ina (uA_per_cm2)'); LEGEND_CONSTANTS(:,127) = strpad('Ka_ina in component ina (uM)'); LEGEND_CONSTANTS(:,128) = strpad('Kp_ina in component ina (uM)'); LEGEND_ALGEBRAIC(:,276) = strpad('f_both in component ina (dimensionless)'); LEGEND_ALGEBRAIC(:,277) = strpad('f_camk_only in component ina (dimensionless)'); LEGEND_STATES(:,106) = strpad('f_ina in component ina (dimensionless)'); LEGEND_ALGEBRAIC(:,279) = strpad('f_np in component ina (dimensionless)'); LEGEND_ALGEBRAIC(:,275) = strpad('f_pka in component ina (dimensionless)'); LEGEND_ALGEBRAIC(:,278) = strpad('f_pka_only in component ina (dimensionless)'); LEGEND_CONSTANTS(:,242) = strpad('gNaBar in component ina (mS_per_cm2)'); LEGEND_STATES(:,107) = strpad('h in component ina_camk (dimensionless)'); LEGEND_STATES(:,108) = strpad('j in component ina_camk (dimensionless)'); LEGEND_STATES(:,109) = strpad('m in component ina_camk (dimensionless)'); LEGEND_ALGEBRAIC(:,274) = strpad('ina_f_pka_val in component ina (dimensionless)'); LEGEND_STATES(:,110) = strpad('h in component ina_np (dimensionless)'); LEGEND_STATES(:,111) = strpad('j in component ina_np (dimensionless)'); LEGEND_STATES(:,112) = strpad('m in component ina_np (dimensionless)'); LEGEND_STATES(:,113) = strpad('h in component ina_pka (dimensionless)'); LEGEND_STATES(:,114) = strpad('j in component ina_pka (dimensionless)'); LEGEND_STATES(:,115) = strpad('m in component ina_pka (dimensionless)'); LEGEND_CONSTANTS(:,129) = strpad('ka_ina in component ina (hertz)'); LEGEND_CONSTANTS(:,130) = strpad('kp_ina in component ina (hertz)'); LEGEND_CONSTANTS(:,131) = strpad('dVIn in component ina_camk (mV)'); LEGEND_ALGEBRAIC(:,17) = strpad('ina_camk_h_alpha in component ina_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,44) = strpad('ina_camk_h_beta in component ina_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,18) = strpad('ina_camk_j_alpha in component ina_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,45) = strpad('ina_camk_j_beta in component ina_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,19) = strpad('ina_camk_m_alpha in component ina_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,46) = strpad('ina_camk_m_beta in component ina_camk (mS_per_uF)'); LEGEND_ALGEBRAIC(:,20) = strpad('ina_np_h_alpha in component ina_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,47) = strpad('ina_np_h_beta in component ina_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,21) = strpad('ina_np_j_alpha in component ina_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,48) = strpad('ina_np_j_beta in component ina_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,22) = strpad('ina_np_m_alpha in component ina_np (mS_per_uF)'); LEGEND_ALGEBRAIC(:,49) = strpad('ina_np_m_beta in component ina_np (mS_per_uF)'); LEGEND_CONSTANTS(:,132) = strpad('dVAc in component ina_pka (mV)'); LEGEND_CONSTANTS(:,133) = strpad('dVIn in component ina_pka (mV)'); LEGEND_ALGEBRAIC(:,23) = strpad('ina_pka_h_alpha in component ina_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,50) = strpad('ina_pka_h_beta in component ina_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,24) = strpad('ina_pka_j_alpha in component ina_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,51) = strpad('ina_pka_j_beta in component ina_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,25) = strpad('ina_pka_m_alpha in component ina_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,52) = strpad('ina_pka_m_beta in component ina_pka (mS_per_uF)'); LEGEND_ALGEBRAIC(:,282) = strpad('INab in component inab (uA_per_cm2)'); LEGEND_CONSTANTS(:,134) = strpad('P in component inab (cm_per_s)'); LEGEND_ALGEBRAIC(:,281) = strpad('ePhi in component inab (dimensionless)'); LEGEND_ALGEBRAIC(:,280) = strpad('inab_INab_phi in component inab (dimensionless)'); LEGEND_CONSTANTS(:,363) = strpad('KmNai3 in component inaca (mM3)'); LEGEND_CONSTANTS(:,365) = strpad('KmNao3 in component inaca (mM3)'); LEGEND_CONSTANTS(:,135) = strpad('Km_Ca in component inaca (mM)'); LEGEND_CONSTANTS(:,136) = strpad('Km_Cai in component inaca (mM)'); LEGEND_CONSTANTS(:,137) = strpad('Km_Cao in component inaca (mM)'); LEGEND_CONSTANTS(:,138) = strpad('Km_Nai in component inaca (mM)'); LEGEND_CONSTANTS(:,139) = strpad('Km_Nao in component inaca (mM)'); LEGEND_ALGEBRAIC(:,283) = strpad('Na_i3 in component inaca (mM3)'); LEGEND_CONSTANTS(:,367) = strpad('Na_o3 in component inaca (mM3)'); LEGEND_ALGEBRAIC(:,284) = strpad('Na_ss3 in component inaca (mM3)'); LEGEND_CONSTANTS(:,140) = strpad('eta in component inaca (dimensionless)'); LEGEND_ALGEBRAIC(:,285) = strpad('exp1 in component inaca (dimensionless)'); LEGEND_ALGEBRAIC(:,286) = strpad('exp2 in component inaca (dimensionless)'); LEGEND_ALGEBRAIC(:,287) = strpad('inaca_INaCaSR_denom1 in component inaca (dimensionless)'); LEGEND_ALGEBRAIC(:,288) = strpad('inaca_INaCaSR_denom2 in component inaca (dimensionless)'); LEGEND_ALGEBRAIC(:,289) = strpad('inaca_INaCaSR_denom3 in component inaca (mol4_per_m12)'); LEGEND_ALGEBRAIC(:,290) = strpad('inaca_INaCaSR_denom4 in component inaca (mol4_per_m12)'); LEGEND_ALGEBRAIC(:,291) = strpad('inaca_INaCaSR_num in component inaca (A_mol4_per_m14_times_1e_minus_2)'); LEGEND_ALGEBRAIC(:,293) = strpad('inaca_INaCa_denom1 in component inaca (dimensionless)'); LEGEND_ALGEBRAIC(:,294) = strpad('inaca_INaCa_denom2 in component inaca (dimensionless)'); LEGEND_ALGEBRAIC(:,295) = strpad('inaca_INaCa_denom3 in component inaca (mol4_per_m12)'); LEGEND_ALGEBRAIC(:,296) = strpad('inaca_INaCa_denom4 in component inaca (mol4_per_m12)'); LEGEND_ALGEBRAIC(:,297) = strpad('inaca_INaCa_num in component inaca (A_mol4_per_m14_times_1e_minus_2)'); LEGEND_CONSTANTS(:,141) = strpad('kSat in component inaca (dimensionless)'); LEGEND_CONSTANTS(:,142) = strpad('vMax in component inaca (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,304) = strpad('INaK in component inak (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,302) = strpad('INaK_np in component inak (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,303) = strpad('INaK_p in component inak (uA_per_cm2)'); LEGEND_CONSTANTS(:,143) = strpad('Ka_inak in component inak (uM)'); LEGEND_CONSTANTS(:,144) = strpad('Kp_inak in component inak (uM)'); LEGEND_STATES(:,116) = strpad('f_inak in component inak (dimensionless)'); LEGEND_ALGEBRAIC(:,300) = strpad('fhat in component inak (dimensionless)'); LEGEND_CONSTANTS(:,145) = strpad('ibar in component inak (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,299) = strpad('inak_fhat_val in component inak (dimensionless)'); LEGEND_CONSTANTS(:,146) = strpad('ka_inak in component inak (hertz)'); LEGEND_CONSTANTS(:,147) = strpad('km_ko in component inak (mM)'); LEGEND_CONSTANTS(:,148) = strpad('km_np in component inak (mM)'); LEGEND_CONSTANTS(:,149) = strpad('km_p in component inak (mM)'); LEGEND_CONSTANTS(:,150) = strpad('kp_inak in component inak (hertz)'); LEGEND_ALGEBRAIC(:,301) = strpad('phi in component inak (uA_per_cm2)'); LEGEND_CONSTANTS(:,243) = strpad('pk in component inak (dimensionless)'); LEGEND_ALGEBRAIC(:,364) = strpad('INaL in component inal (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,362) = strpad('INaL_camk in component inal (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,363) = strpad('INaL_np in component inal (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,361) = strpad('conductance in component inal (mV)'); LEGEND_STATES(:,117) = strpad('h in component inal (dimensionless)'); LEGEND_ALGEBRAIC(:,26) = strpad('h_inf in component inal (dimensionless)'); LEGEND_ALGEBRAIC(:,27) = strpad('inal_m_alpha in component inal (mS_per_uF)'); LEGEND_ALGEBRAIC(:,53) = strpad('inal_m_beta in component inal (mS_per_uF)'); LEGEND_STATES(:,118) = strpad('m in component inal (dimensionless)'); LEGEND_CONSTANTS(:,151) = strpad('tau_h in component inal (ms)'); LEGEND_CONSTANTS(:,152) = strpad('IpCa_bar in component ipca (uA_per_cm2)'); LEGEND_CONSTANTS(:,153) = strpad('Km_pCa in component ipca (mM)'); LEGEND_ALGEBRAIC(:,316) = strpad('Ileak_ryr in component irel (mM_per_ms)'); LEGEND_ALGEBRAIC(:,308) = strpad('Ileak_ryr_np in component irel (mM_per_ms)'); LEGEND_ALGEBRAIC(:,310) = strpad('Ileak_ryr_p in component irel (mM_per_ms)'); LEGEND_STATES(:,119) = strpad('Irel_np in component irel (mM_per_ms)'); LEGEND_STATES(:,120) = strpad('Irel_p in component irel (mM_per_ms)'); LEGEND_CONSTANTS(:,154) = strpad('Km_ryr_leak_np in component irel (mM)'); LEGEND_CONSTANTS(:,155) = strpad('Km_ryr_leak_p in component irel (mM)'); LEGEND_ALGEBRAIC(:,256) = strpad('alpha_np in component irel (mM_per_ms)'); LEGEND_ALGEBRAIC(:,262) = strpad('alpha_p in component irel (mM_per_ms)'); LEGEND_CONSTANTS(:,244) = strpad('beta_0 in component irel (ms)'); LEGEND_ALGEBRAIC(:,253) = strpad('beta_np in component irel (ms)'); LEGEND_ALGEBRAIC(:,259) = strpad('beta_p in component irel (ms)'); LEGEND_ALGEBRAIC(:,314) = strpad('fhat in component irel (dimensionless)'); LEGEND_CONSTANTS(:,399) = strpad('irel_fhat_ratio in component irel (dimensionless)'); LEGEND_ALGEBRAIC(:,312) = strpad('irel_fhat_val in component irel (dimensionless)'); LEGEND_ALGEBRAIC(:,272) = strpad('irel_inf_np in component irel (mM_per_ms)'); LEGEND_ALGEBRAIC(:,273) = strpad('irel_inf_p in component irel (mM_per_ms)'); LEGEND_ALGEBRAIC(:,266) = strpad('irel_tau_np in component irel (ms)'); LEGEND_ALGEBRAIC(:,268) = strpad('irel_tau_p in component irel (ms)'); LEGEND_CONSTANTS(:,156) = strpad('k_ryr_leak_np in component irel (mS_per_uF)'); LEGEND_CONSTANTS(:,157) = strpad('k_ryr_leak_p in component irel (mS_per_uF)'); LEGEND_ALGEBRAIC(:,264) = strpad('x in component irel (dimensionless)'); LEGEND_ALGEBRAIC(:,270) = strpad('y in component irel (dimensionless)'); LEGEND_CONSTANTS(:,158) = strpad('Gbar in component ito (mS_per_cm2)'); LEGEND_ALGEBRAIC(:,348) = strpad('ITo in component ito (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,346) = strpad('ITo_camk in component ito (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,347) = strpad('ITo_np in component ito (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,326) = strpad('R in component ito (dimensionless)'); LEGEND_ALGEBRAIC(:,28) = strpad('a_inf in component ito (dimensionless)'); LEGEND_STATES(:,121) = strpad('a_np in component ito (dimensionless)'); LEGEND_ALGEBRAIC(:,78) = strpad('a_tau in component ito (ms)'); LEGEND_ALGEBRAIC(:,29) = strpad('alph_if in component ito (dimensionless)'); LEGEND_ALGEBRAIC(:,55) = strpad('alph_is in component ito (dimensionless)'); LEGEND_ALGEBRAIC(:,67) = strpad('beta_i in component ito (mS_per_uF)'); LEGEND_STATES(:,122) = strpad('if_camk in component ito (dimensionless)'); LEGEND_STATES(:,123) = strpad('if_np in component ito (dimensionless)'); LEGEND_STATES(:,124) = strpad('is_camk in component ito (dimensionless)'); LEGEND_STATES(:,125) = strpad('is_np in component ito (dimensionless)'); LEGEND_ALGEBRAIC(:,54) = strpad('ito_a_np_alpha in component ito (dimensionless)'); LEGEND_ALGEBRAIC(:,66) = strpad('ito_a_np_beta in component ito (dimensionless)'); LEGEND_ALGEBRAIC(:,79) = strpad('ito_if_camk_alpha in component ito (mS_per_uF)'); LEGEND_ALGEBRAIC(:,80) = strpad('ito_if_np_alpha in component ito (mS_per_uF)'); LEGEND_ALGEBRAIC(:,81) = strpad('ito_is_camk_alpha in component ito (mS_per_uF)'); LEGEND_ALGEBRAIC(:,82) = strpad('ito_is_np_alpha in component ito (mS_per_uF)'); LEGEND_ALGEBRAIC(:,345) = strpad('x in component ito (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,309) = strpad('Imax in component iup (mM_per_ms)'); LEGEND_CONSTANTS(:,159) = strpad('Ka_plb in component iup (uM)'); LEGEND_CONSTANTS(:,377) = strpad('Km_both in component iup (mM)'); LEGEND_CONSTANTS(:,370) = strpad('Km_camk in component iup (mM)'); LEGEND_CONSTANTS(:,160) = strpad('Km_np in component iup (mM)'); LEGEND_CONSTANTS(:,375) = strpad('Km_pka in component iup (mM)'); LEGEND_ALGEBRAIC(:,323) = strpad('Km_up in component iup (mM)'); LEGEND_CONSTANTS(:,161) = strpad('Kp_plb in component iup (uM)'); LEGEND_ALGEBRAIC(:,307) = strpad('f_SERCA2a in component iup (dimensionless)'); LEGEND_ALGEBRAIC(:,315) = strpad('f_both in component iup (dimensionless)'); LEGEND_ALGEBRAIC(:,317) = strpad('f_camk_only in component iup (dimensionless)'); LEGEND_ALGEBRAIC(:,322) = strpad('f_np in component iup (dimensionless)'); LEGEND_ALGEBRAIC(:,313) = strpad('f_pka in component iup (dimensionless)'); LEGEND_ALGEBRAIC(:,319) = strpad('f_pka_only in component iup (dimensionless)'); LEGEND_STATES(:,126) = strpad('f_plb in component iup (dimensionless)'); LEGEND_ALGEBRAIC(:,311) = strpad('iup_f_pka_val in component iup (dimensionless)'); LEGEND_CONSTANTS(:,162) = strpad('iupmax in component iup (mM_per_ms)'); LEGEND_CONSTANTS(:,380) = strpad('iupmaxCAMK in component iup (mM_per_ms)'); LEGEND_CONSTANTS(:,163) = strpad('ka_plb in component iup (hertz)'); LEGEND_CONSTANTS(:,164) = strpad('kp_plb in component iup (hertz)'); LEGEND_ALGEBRAIC(:,325) = strpad('leak in component iup (mM_per_ms)'); LEGEND_CONSTANTS(:,165) = strpad('nsrmax in component iup (mM)'); LEGEND_ALGEBRAIC(:,327) = strpad('uptake in component iup (mM_per_ms)'); LEGEND_ALGEBRAIC(:,365) = strpad('IK_tot in component potassium (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,367) = strpad('INa_tot in component sodium (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,368) = strpad('i_ion in component membrane (uA_per_cm2)'); LEGEND_ALGEBRAIC(:,369) = strpad('i_stim in component stimulus (uA_per_cm2)'); LEGEND_CONSTANTS(:,166) = strpad('PNaK in component nernst (dimensionless)'); LEGEND_CONSTANTS(:,239) = strpad('RTF in component phys (mV)'); LEGEND_STATES(:,127) = strpad('K in component potassium (mM)'); LEGEND_CONSTANTS(:,167) = strpad('KPDEp in component pde (uM)'); LEGEND_CONSTANTS(:,168) = strpad('KmIbmxPde2 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,169) = strpad('KmIbmxPde3 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,170) = strpad('KmIbmxPde4 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,171) = strpad('KmPDE2 in component pde (uM)'); LEGEND_CONSTANTS(:,172) = strpad('KmPDE3 in component pde (uM)'); LEGEND_CONSTANTS(:,173) = strpad('KmPDE4 in component pde (uM)'); LEGEND_CONSTANTS(:,348) = strpad('PDE2_cav in component pde (uM)'); LEGEND_CONSTANTS(:,349) = strpad('PDE2_cyt in component pde (uM)'); LEGEND_CONSTANTS(:,350) = strpad('PDE2_eca in component pde (uM)'); LEGEND_CONSTANTS(:,174) = strpad('PDE2_tot in component pde (uM)'); LEGEND_STATES(:,128) = strpad('PDE3_P_cav in component pde (uM)'); LEGEND_STATES(:,129) = strpad('PDE3_P_cyt in component pde (uM)'); LEGEND_CONSTANTS(:,378) = strpad('PDE3_cav in component pde (uM)'); LEGEND_CONSTANTS(:,379) = strpad('PDE3_cyt in component pde (uM)'); LEGEND_CONSTANTS(:,364) = strpad('PDE3_tot in component pde (uM)'); LEGEND_STATES(:,130) = strpad('PDE4_P_cav in component pde (uM)'); LEGEND_STATES(:,131) = strpad('PDE4_P_cyt in component pde (uM)'); LEGEND_STATES(:,132) = strpad('PDE4_P_eca in component pde (uM)'); LEGEND_CONSTANTS(:,371) = strpad('PDE4_cav in component pde (uM)'); LEGEND_CONSTANTS(:,372) = strpad('PDE4_cyt in component pde (uM)'); LEGEND_CONSTANTS(:,373) = strpad('PDE4_eca in component pde (uM)'); LEGEND_CONSTANTS(:,366) = strpad('PDE4_tot in component pde (uM)'); LEGEND_CONSTANTS(:,175) = strpad('delta_k_pde34 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,176) = strpad('f_pde2_cav in component pde (dimensionless)'); LEGEND_CONSTANTS(:,334) = strpad('f_pde2_cyt in component pde (dimensionless)'); LEGEND_CONSTANTS(:,177) = strpad('f_pde2_eca in component pde (dimensionless)'); LEGEND_CONSTANTS(:,335) = strpad('f_pde2_part in component pde (dimensionless)'); LEGEND_CONSTANTS(:,368) = strpad('f_pde3_cav in component pde (dimensionless)'); LEGEND_CONSTANTS(:,374) = strpad('f_pde3_cyt in component pde (dimensionless)'); LEGEND_CONSTANTS(:,178) = strpad('f_pde4_cav in component pde (dimensionless)'); LEGEND_CONSTANTS(:,340) = strpad('f_pde4_cyt in component pde (dimensionless)'); LEGEND_CONSTANTS(:,344) = strpad('f_pde4_eca in component pde (dimensionless)'); LEGEND_CONSTANTS(:,179) = strpad('f_pde4_part in component pde (dimensionless)'); LEGEND_CONSTANTS(:,180) = strpad('f_pde_part in component pde (dimensionless)'); LEGEND_CONSTANTS(:,181) = strpad('ff_pde3_cyt in component pde (dimensionless)'); LEGEND_CONSTANTS(:,182) = strpad('h_ibmx_pde2 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,183) = strpad('h_ibmx_pde3 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,184) = strpad('h_ibmx_pde4 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,185) = strpad('ibmx in component pde (uM)'); LEGEND_CONSTANTS(:,343) = strpad('ibmx2 in component pde (uM)'); LEGEND_CONSTANTS(:,376) = strpad('ibmx3 in component pde (uM)'); LEGEND_CONSTANTS(:,369) = strpad('ibmx4 in component pde (uM)'); LEGEND_CONSTANTS(:,339) = strpad('ibmx_h2 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,351) = strpad('ibmx_h3 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,355) = strpad('ibmx_h4 in component pde (dimensionless)'); LEGEND_CONSTANTS(:,186) = strpad('kPDE2 in component pde (hertz)'); LEGEND_CONSTANTS(:,187) = strpad('kPDE3 in component pde (hertz)'); LEGEND_CONSTANTS(:,188) = strpad('kPDE4 in component pde (hertz)'); LEGEND_CONSTANTS(:,245) = strpad('kbPDEp in component pde (hertz)'); LEGEND_CONSTANTS(:,189) = strpad('kfPDEp in component pde (per_mM_per_ms)'); LEGEND_CONSTANTS(:,360) = strpad('pde_PDE3_tot_alpha in component pde (dimensionless)'); LEGEND_CONSTANTS(:,362) = strpad('pde_PDE3_tot_beta in component pde (dimensionless)'); LEGEND_CONSTANTS(:,190) = strpad('r_pde34_frac in component pde (dimensionless)'); LEGEND_CONSTANTS(:,358) = strpad('r_pde3_cyt in component pde (dimensionless)'); LEGEND_CONSTANTS(:,191) = strpad('R in component phys (mJ_per_mol_per_K)'); LEGEND_CONSTANTS(:,192) = strpad('T in component phys (kelvin)'); LEGEND_CONSTANTS(:,246) = strpad('K_pki in component pka (uM)'); LEGEND_CONSTANTS(:,336) = strpad('PKA_cyt in component pka (uM)'); LEGEND_CONSTANTS(:,193) = strpad('PKA_tot in component pka (uM)'); LEGEND_CONSTANTS(:,345) = strpad('PKI_cav in component pka (uM)'); LEGEND_CONSTANTS(:,356) = strpad('PKI_cyt in component pka (uM)'); LEGEND_CONSTANTS(:,352) = strpad('PKI_eca in component pka (uM)'); LEGEND_CONSTANTS(:,337) = strpad('PKI_tot in component pka (uM)'); LEGEND_CONSTANTS(:,265) = strpad('b_pki in component pka (hertz)'); LEGEND_CONSTANTS(:,194) = strpad('f_cav in component pka (dimensionless)'); LEGEND_CONSTANTS(:,328) = strpad('f_cyt in component pka (dimensionless)'); LEGEND_CONSTANTS(:,195) = strpad('f_eca in component pka (dimensionless)'); LEGEND_CONSTANTS(:,196) = strpad('f_pki in component pka (per_mM_per_ms)'); LEGEND_CONSTANTS(:,341) = strpad('f_pki_cav in component pka (dimensionless)'); LEGEND_CONSTANTS(:,353) = strpad('f_pki_cyt in component pka (dimensionless)'); LEGEND_CONSTANTS(:,346) = strpad('f_pki_eca in component pka (dimensionless)'); LEGEND_STATES(:,133) = strpad('A2R in component pka_cav (uM)'); LEGEND_STATES(:,134) = strpad('A2RC in component pka_cav (uM)'); LEGEND_STATES(:,135) = strpad('ARC in component pka_cav (uM)'); LEGEND_CONSTANTS(:,197) = strpad('K1 in component pka_cav (uM)'); LEGEND_CONSTANTS(:,198) = strpad('K2 in component pka_cav (uM)'); LEGEND_CONSTANTS(:,199) = strpad('K3 in component pka_cav (m3_per_mol_times_1e3)'); LEGEND_STATES(:,136) = strpad('PKIC in component pka_cav (uM)'); LEGEND_ALGEBRAIC(:,141) = strpad('RCf in component pka_cav (uM)'); LEGEND_CONSTANTS(:,247) = strpad('b1 in component pka_cav (hertz)'); LEGEND_CONSTANTS(:,248) = strpad('b2 in component pka_cav (hertz)'); LEGEND_CONSTANTS(:,249) = strpad('b3 in component pka_cav (per_mM_per_ms)'); LEGEND_CONSTANTS(:,200) = strpad('f1 in component pka_cav (per_mM_per_ms)'); LEGEND_CONSTANTS(:,201) = strpad('f2 in component pka_cav (per_mM_per_ms)'); LEGEND_CONSTANTS(:,202) = strpad('f3 in component pka_cav (hertz)'); LEGEND_STATES(:,137) = strpad('A2R in component pka_cyt (uM)'); LEGEND_STATES(:,138) = strpad('A2RC in component pka_cyt (uM)'); LEGEND_STATES(:,139) = strpad('ARC in component pka_cyt (uM)'); LEGEND_CONSTANTS(:,203) = strpad('K1 in component pka_cyt (uM)'); LEGEND_CONSTANTS(:,204) = strpad('K2 in component pka_cyt (uM)'); LEGEND_CONSTANTS(:,205) = strpad('K3 in component pka_cyt (m3_per_mol_times_1e3)'); LEGEND_STATES(:,140) = strpad('PKIC in component pka_cyt (uM)'); LEGEND_ALGEBRAIC(:,151) = strpad('RCf in component pka_cyt (uM)'); LEGEND_CONSTANTS(:,266) = strpad('b1 in component pka_cyt (hertz)'); LEGEND_CONSTANTS(:,267) = strpad('b2 in component pka_cyt (hertz)'); LEGEND_CONSTANTS(:,268) = strpad('b3 in component pka_cyt (per_mM_per_ms)'); LEGEND_CONSTANTS(:,250) = strpad('f1 in component pka_cyt (per_mM_per_ms)'); LEGEND_CONSTANTS(:,251) = strpad('f2 in component pka_cyt (per_mM_per_ms)'); LEGEND_CONSTANTS(:,252) = strpad('f3 in component pka_cyt (hertz)'); LEGEND_STATES(:,141) = strpad('A2R in component pka_eca (uM)'); LEGEND_STATES(:,142) = strpad('A2RC in component pka_eca (uM)'); LEGEND_STATES(:,143) = strpad('ARC in component pka_eca (uM)'); LEGEND_CONSTANTS(:,253) = strpad('K1 in component pka_eca (uM)'); LEGEND_CONSTANTS(:,254) = strpad('K2 in component pka_eca (uM)'); LEGEND_CONSTANTS(:,255) = strpad('K3 in component pka_eca (m3_per_mol_times_1e3)'); LEGEND_STATES(:,144) = strpad('PKIC in component pka_eca (uM)'); LEGEND_ALGEBRAIC(:,131) = strpad('RCf in component pka_eca (uM)'); LEGEND_CONSTANTS(:,277) = strpad('b1 in component pka_eca (hertz)'); LEGEND_CONSTANTS(:,278) = strpad('b2 in component pka_eca (hertz)'); LEGEND_CONSTANTS(:,279) = strpad('b3 in component pka_eca (per_mM_per_ms)'); LEGEND_CONSTANTS(:,269) = strpad('f1 in component pka_eca (per_mM_per_ms)'); LEGEND_CONSTANTS(:,270) = strpad('f2 in component pka_eca (per_mM_per_ms)'); LEGEND_CONSTANTS(:,271) = strpad('f3 in component pka_eca (hertz)'); LEGEND_CONSTANTS(:,407) = strpad('potassium_K_r1 in component potassium (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_CONSTANTS(:,206) = strpad('K in component pp1 (uM)'); LEGEND_CONSTANTS(:,207) = strpad('Kdp in component pp1 (uM)'); LEGEND_CONSTANTS(:,208) = strpad('Kp in component pp1 (uM)'); LEGEND_CONSTANTS(:,209) = strpad('PP1_cyt in component pp1 (uM)'); LEGEND_ALGEBRAIC(:,30) = strpad('di in component pp1 (uM)'); LEGEND_CONSTANTS(:,210) = strpad('f in component pp1 (dimensionless)'); LEGEND_STATES(:,145) = strpad('inhib1_p in component pp1 (uM)'); LEGEND_CONSTANTS(:,256) = strpad('inhib1_tot in component pp1 (uM)'); LEGEND_CONSTANTS(:,211) = strpad('kdp in component pp1 (hertz)'); LEGEND_CONSTANTS(:,212) = strpad('kp in component pp1 (hertz)'); LEGEND_ALGEBRAIC(:,199) = strpad('pp1_PP1f_cyt_sum in component pp1 (uM)'); LEGEND_ALGEBRAIC(:,366) = strpad('INa_cyt in component sodium (uA_per_cm2)'); LEGEND_CONSTANTS(:,408) = strpad('sodium_Na_r1 in component sodium (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_CONSTANTS(:,409) = strpad('sodium_Na_r2 in component sodium (dimensionless)'); LEGEND_CONSTANTS(:,410) = strpad('sodium_Na_sr_r1 in component sodium (mol_per_m_per_s_per_A_times_1e5)'); LEGEND_CONSTANTS(:,257) = strpad('amplitude in component stimulus (uA_per_cm2)'); LEGEND_CONSTANTS(:,213) = strpad('duration in component stimulus (ms)'); LEGEND_CONSTANTS(:,214) = strpad('offset in component stimulus (ms)'); LEGEND_CONSTANTS(:,215) = strpad('period in component stimulus (ms)'); LEGEND_RATES(:,5) = strpad('d/dt ICaLp in component akap_sig (uM)'); LEGEND_RATES(:,6) = strpad('d/dt RyRp in component akap_sig (uM)'); LEGEND_RATES(:,8) = strpad('d/dt Gi_aGDP in component beta_cav (uM)'); LEGEND_RATES(:,9) = strpad('d/dt Gi_aGTP in component beta_cav (uM)'); LEGEND_RATES(:,1) = strpad('d/dt Gi_bg in component beta_cav (uM)'); LEGEND_RATES(:,10) = strpad('d/dt Gs_aGDP in component beta_cav (uM)'); LEGEND_RATES(:,2) = strpad('d/dt Gs_aGTP in component beta_cav (uM)'); LEGEND_RATES(:,11) = strpad('d/dt Gs_bg in component beta_cav (uM)'); LEGEND_RATES(:,12) = strpad('d/dt Rb1_grk_tot in component beta_cav (uM)'); LEGEND_RATES(:,13) = strpad('d/dt Rb1_pka_tot in component beta_cav (uM)'); LEGEND_RATES(:,14) = strpad('d/dt Rb2_grk_tot in component beta_cav (uM)'); LEGEND_RATES(:,15) = strpad('d/dt Rb2_pka_tot in component beta_cav (uM)'); LEGEND_RATES(:,16) = strpad('d/dt Gs_aGDP in component beta_cyt (uM)'); LEGEND_RATES(:,3) = strpad('d/dt Gs_aGTP in component beta_cyt (uM)'); LEGEND_RATES(:,17) = strpad('d/dt Gs_bg in component beta_cyt (uM)'); LEGEND_RATES(:,18) = strpad('d/dt Rb1_grk_tot in component beta_cyt (uM)'); LEGEND_RATES(:,19) = strpad('d/dt Rb1_pka_tot in component beta_cyt (uM)'); LEGEND_RATES(:,21) = strpad('d/dt Gi_aGDP in component beta_eca (uM)'); LEGEND_RATES(:,22) = strpad('d/dt Gi_aGTP in component beta_eca (uM)'); LEGEND_RATES(:,23) = strpad('d/dt Gi_bg in component beta_eca (uM)'); LEGEND_RATES(:,24) = strpad('d/dt Gs_aGDP in component beta_eca (uM)'); LEGEND_RATES(:,4) = strpad('d/dt Gs_aGTP in component beta_eca (uM)'); LEGEND_RATES(:,25) = strpad('d/dt Gs_bg in component beta_eca (uM)'); LEGEND_RATES(:,26) = strpad('d/dt Rb1_grk_tot in component beta_eca (uM)'); LEGEND_RATES(:,27) = strpad('d/dt Rb1_pka_tot in component beta_eca (uM)'); LEGEND_RATES(:,28) = strpad('d/dt Rb2_grk_tot in component beta_eca (uM)'); LEGEND_RATES(:,29) = strpad('d/dt Rb2_pka_tot in component beta_eca (uM)'); LEGEND_RATES(:,31) = strpad('d/dt Ca_nsr in component calcium (mM)'); LEGEND_RATES(:,32) = strpad('d/dt f_tni in component calcium (dimensionless)'); LEGEND_RATES(:,33) = strpad('d/dt uCa in component calcium (mM)'); LEGEND_RATES(:,34) = strpad('d/dt uCa_CaL in component calcium (mM)'); LEGEND_RATES(:,35) = strpad('d/dt uCa_jsr in component calcium (mM)'); LEGEND_RATES(:,36) = strpad('d/dt uCa_sr in component calcium (mM)'); LEGEND_RATES(:,37) = strpad('d/dt f_ical in component camk (dimensionless)'); LEGEND_RATES(:,38) = strpad('d/dt f_ik1 in component camk (dimensionless)'); LEGEND_RATES(:,39) = strpad('d/dt f_ina in component camk (dimensionless)'); LEGEND_RATES(:,40) = strpad('d/dt f_ito in component camk (dimensionless)'); LEGEND_RATES(:,41) = strpad('d/dt f_plb in component camk (dimensionless)'); LEGEND_RATES(:,42) = strpad('d/dt f_ryr in component camk (dimensionless)'); LEGEND_RATES(:,43) = strpad('d/dt trap in component camk (dimensionless)'); LEGEND_RATES(:,44) = strpad('d/dt cAMP_cav in component camp (uM)'); LEGEND_RATES(:,45) = strpad('d/dt cAMP_cyt in component camp (uM)'); LEGEND_RATES(:,46) = strpad('d/dt cAMP_eca in component camp (uM)'); LEGEND_RATES(:,47) = strpad('d/dt Cl in component chloride (mM)'); LEGEND_RATES(:,48) = strpad('d/dt Cl_sr in component chloride (mM)'); LEGEND_RATES(:,52) = strpad('d/dt C in component ical_camk (dimensionless)'); LEGEND_RATES(:,53) = strpad('d/dt CI in component ical_camk (dimensionless)'); LEGEND_RATES(:,54) = strpad('d/dt CIs in component ical_camk (dimensionless)'); LEGEND_RATES(:,55) = strpad('d/dt Cs in component ical_camk (dimensionless)'); LEGEND_RATES(:,56) = strpad('d/dt O in component ical_camk (dimensionless)'); LEGEND_RATES(:,57) = strpad('d/dt OI in component ical_camk (dimensionless)'); LEGEND_RATES(:,58) = strpad('d/dt OIs in component ical_camk (dimensionless)'); LEGEND_RATES(:,59) = strpad('d/dt Os in component ical_camk (dimensionless)'); LEGEND_RATES(:,60) = strpad('d/dt C in component ical_np (dimensionless)'); LEGEND_RATES(:,61) = strpad('d/dt CI in component ical_np (dimensionless)'); LEGEND_RATES(:,62) = strpad('d/dt CIs in component ical_np (dimensionless)'); LEGEND_RATES(:,63) = strpad('d/dt Cs in component ical_np (dimensionless)'); LEGEND_RATES(:,64) = strpad('d/dt O in component ical_np (dimensionless)'); LEGEND_RATES(:,65) = strpad('d/dt OI in component ical_np (dimensionless)'); LEGEND_RATES(:,66) = strpad('d/dt OIs in component ical_np (dimensionless)'); LEGEND_RATES(:,67) = strpad('d/dt Os in component ical_np (dimensionless)'); LEGEND_RATES(:,68) = strpad('d/dt i2 in component iclca (dimensionless)'); LEGEND_RATES(:,69) = strpad('d/dt ac in component ikr (dimensionless)'); LEGEND_RATES(:,74) = strpad('d/dt C1 in component iks_np (dimensionless)'); LEGEND_RATES(:,75) = strpad('d/dt C10 in component iks_np (dimensionless)'); LEGEND_RATES(:,76) = strpad('d/dt C11 in component iks_np (dimensionless)'); LEGEND_RATES(:,77) = strpad('d/dt C12 in component iks_np (dimensionless)'); LEGEND_RATES(:,78) = strpad('d/dt C13 in component iks_np (dimensionless)'); LEGEND_RATES(:,79) = strpad('d/dt C14 in component iks_np (dimensionless)'); LEGEND_RATES(:,80) = strpad('d/dt C15 in component iks_np (dimensionless)'); LEGEND_RATES(:,81) = strpad('d/dt C2 in component iks_np (dimensionless)'); LEGEND_RATES(:,82) = strpad('d/dt C3 in component iks_np (dimensionless)'); LEGEND_RATES(:,83) = strpad('d/dt C4 in component iks_np (dimensionless)'); LEGEND_RATES(:,84) = strpad('d/dt C5 in component iks_np (dimensionless)'); LEGEND_RATES(:,85) = strpad('d/dt C6 in component iks_np (dimensionless)'); LEGEND_RATES(:,86) = strpad('d/dt C7 in component iks_np (dimensionless)'); LEGEND_RATES(:,87) = strpad('d/dt C8 in component iks_np (dimensionless)'); LEGEND_RATES(:,88) = strpad('d/dt C9 in component iks_np (dimensionless)'); LEGEND_RATES(:,70) = strpad('d/dt O1 in component iks_np (dimensionless)'); LEGEND_RATES(:,71) = strpad('d/dt O2 in component iks_np (dimensionless)'); LEGEND_RATES(:,89) = strpad('d/dt C1 in component iks_pka (dimensionless)'); LEGEND_RATES(:,90) = strpad('d/dt C10 in component iks_pka (dimensionless)'); LEGEND_RATES(:,91) = strpad('d/dt C11 in component iks_pka (dimensionless)'); LEGEND_RATES(:,92) = strpad('d/dt C12 in component iks_pka (dimensionless)'); LEGEND_RATES(:,93) = strpad('d/dt C13 in component iks_pka (dimensionless)'); LEGEND_RATES(:,94) = strpad('d/dt C14 in component iks_pka (dimensionless)'); LEGEND_RATES(:,95) = strpad('d/dt C15 in component iks_pka (dimensionless)'); LEGEND_RATES(:,96) = strpad('d/dt C2 in component iks_pka (dimensionless)'); LEGEND_RATES(:,97) = strpad('d/dt C3 in component iks_pka (dimensionless)'); LEGEND_RATES(:,98) = strpad('d/dt C4 in component iks_pka (dimensionless)'); LEGEND_RATES(:,99) = strpad('d/dt C5 in component iks_pka (dimensionless)'); LEGEND_RATES(:,100) = strpad('d/dt C6 in component iks_pka (dimensionless)'); LEGEND_RATES(:,101) = strpad('d/dt C7 in component iks_pka (dimensionless)'); LEGEND_RATES(:,102) = strpad('d/dt C8 in component iks_pka (dimensionless)'); LEGEND_RATES(:,103) = strpad('d/dt C9 in component iks_pka (dimensionless)'); LEGEND_RATES(:,72) = strpad('d/dt O1 in component iks_pka (dimensionless)'); LEGEND_RATES(:,73) = strpad('d/dt O2 in component iks_pka (dimensionless)'); LEGEND_RATES(:,104) = strpad('d/dt IKsp in component iks_sig (uM)'); LEGEND_RATES(:,105) = strpad('d/dt f_ikur in component ikur (dimensionless)'); LEGEND_RATES(:,106) = strpad('d/dt f_ina in component ina (dimensionless)'); LEGEND_RATES(:,107) = strpad('d/dt h in component ina_camk (dimensionless)'); LEGEND_RATES(:,108) = strpad('d/dt j in component ina_camk (dimensionless)'); LEGEND_RATES(:,109) = strpad('d/dt m in component ina_camk (dimensionless)'); LEGEND_RATES(:,110) = strpad('d/dt h in component ina_np (dimensionless)'); LEGEND_RATES(:,111) = strpad('d/dt j in component ina_np (dimensionless)'); LEGEND_RATES(:,112) = strpad('d/dt m in component ina_np (dimensionless)'); LEGEND_RATES(:,113) = strpad('d/dt h in component ina_pka (dimensionless)'); LEGEND_RATES(:,114) = strpad('d/dt j in component ina_pka (dimensionless)'); LEGEND_RATES(:,115) = strpad('d/dt m in component ina_pka (dimensionless)'); LEGEND_RATES(:,116) = strpad('d/dt f_inak in component inak (dimensionless)'); LEGEND_RATES(:,117) = strpad('d/dt h in component inal (dimensionless)'); LEGEND_RATES(:,118) = strpad('d/dt m in component inal (dimensionless)'); LEGEND_RATES(:,119) = strpad('d/dt Irel_np in component irel (mM_per_ms)'); LEGEND_RATES(:,120) = strpad('d/dt Irel_p in component irel (mM_per_ms)'); LEGEND_RATES(:,121) = strpad('d/dt a_np in component ito (dimensionless)'); LEGEND_RATES(:,122) = strpad('d/dt if_camk in component ito (dimensionless)'); LEGEND_RATES(:,123) = strpad('d/dt if_np in component ito (dimensionless)'); LEGEND_RATES(:,124) = strpad('d/dt is_camk in component ito (dimensionless)'); LEGEND_RATES(:,125) = strpad('d/dt is_np in component ito (dimensionless)'); LEGEND_RATES(:,126) = strpad('d/dt f_plb in component iup (dimensionless)'); LEGEND_RATES(:,51) = strpad('d/dt V in component membrane (mV)'); LEGEND_RATES(:,128) = strpad('d/dt PDE3_P_cav in component pde (uM)'); LEGEND_RATES(:,129) = strpad('d/dt PDE3_P_cyt in component pde (uM)'); LEGEND_RATES(:,130) = strpad('d/dt PDE4_P_cav in component pde (uM)'); LEGEND_RATES(:,131) = strpad('d/dt PDE4_P_cyt in component pde (uM)'); LEGEND_RATES(:,132) = strpad('d/dt PDE4_P_eca in component pde (uM)'); LEGEND_RATES(:,133) = strpad('d/dt A2R in component pka_cav (uM)'); LEGEND_RATES(:,134) = strpad('d/dt A2RC in component pka_cav (uM)'); LEGEND_RATES(:,135) = strpad('d/dt ARC in component pka_cav (uM)'); LEGEND_RATES(:,7) = strpad('d/dt C in component pka_cav (uM)'); LEGEND_RATES(:,136) = strpad('d/dt PKIC in component pka_cav (uM)'); LEGEND_RATES(:,137) = strpad('d/dt A2R in component pka_cyt (uM)'); LEGEND_RATES(:,138) = strpad('d/dt A2RC in component pka_cyt (uM)'); LEGEND_RATES(:,139) = strpad('d/dt ARC in component pka_cyt (uM)'); LEGEND_RATES(:,20) = strpad('d/dt C in component pka_cyt (uM)'); LEGEND_RATES(:,140) = strpad('d/dt PKIC in component pka_cyt (uM)'); LEGEND_RATES(:,141) = strpad('d/dt A2R in component pka_eca (uM)'); LEGEND_RATES(:,142) = strpad('d/dt A2RC in component pka_eca (uM)'); LEGEND_RATES(:,143) = strpad('d/dt ARC in component pka_eca (uM)'); LEGEND_RATES(:,30) = strpad('d/dt C in component pka_eca (uM)'); LEGEND_RATES(:,144) = strpad('d/dt PKIC in component pka_eca (uM)'); LEGEND_RATES(:,127) = strpad('d/dt K in component potassium (mM)'); LEGEND_RATES(:,145) = strpad('d/dt inhib1_p in component pp1 (uM)'); LEGEND_RATES(:,49) = strpad('d/dt Na in component sodium (mM)'); LEGEND_RATES(:,50) = strpad('d/dt Na_sr in component sodium (mM)'); 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) = 5000.0; CONSTANTS(:,2) = 315.0; CONSTANTS(:,3) = 0.0465; CONSTANTS(:,4) = 0.031544; CONSTANTS(:,5) = 0.0852; CONSTANTS(:,6) = 0.4824; CONSTANTS(:,7) = 3.3757; CONSTANTS(:,8) = 41.32; CONSTANTS(:,9) = 0.03135; CONSTANTS(:,10) = 0.037696; STATES(:,1) = 2.09911481235842013e-3; STATES(:,2) = 6.85041638458664965e-3; STATES(:,3) = 7.31420577213055985e-4; STATES(:,4) = 1.84627603007976003e-2; CONSTANTS(:,11) = 0.16479; CONSTANTS(:,12) = 0.087459; CONSTANTS(:,13) = 1.0043; CONSTANTS(:,14) = 1.3574; CONSTANTS(:,15) = 0.6623; CONSTANTS(:,16) = 0.8569; CONSTANTS(:,17) = 0.025; CONSTANTS(:,18) = 0.025; STATES(:,5) = 6.73713947839316954e-4; CONSTANTS(:,19) = 1.27019999999999993e-6; CONSTANTS(:,20) = 6.62979999999999944e-5; CONSTANTS(:,21) = 0.01; CONSTANTS(:,22) = 0.0063064; CONSTANTS(:,23) = 0.043003; CONSTANTS(:,24) = 0.01; CONSTANTS(:,25) = 0.0001; CONSTANTS(:,26) = 0.0001; CONSTANTS(:,27) = 0.01; CONSTANTS(:,28) = 0.01; CONSTANTS(:,29) = 0.25; CONSTANTS(:,30) = 0.125; CONSTANTS(:,31) = 0.125; STATES(:,6) = 4.10693810508170991e-3; CONSTANTS(:,32) = 5.10090000000000044e-4; CONSTANTS(:,33) = 0.0025548; CONSTANTS(:,34) = 0.0006903; CONSTANTS(:,35) = 0.0038257; STATES(:,7) = 3.26565916584702978e-02; CONSTANTS(:,36) = 0.5; CONSTANTS(:,37) = 0.85; CONSTANTS(:,38) = 0.0011071; CONSTANTS(:,39) = 0.5664; CONSTANTS(:,40) = 0.081161; CONSTANTS(:,41) = 0.48744; CONSTANTS(:,42) = 0.85; CONSTANTS(:,43) = 4.0; CONSTANTS(:,44) = 4.9054; CONSTANTS(:,45) = 0.05; CONSTANTS(:,46) = 0.25945; CONSTANTS(:,47) = 2.449; CONSTANTS(:,48) = 0.062; CONSTANTS(:,49) = 0.567; CONSTANTS(:,50) = 1.6655; CONSTANTS(:,51) = 1.8463; CONSTANTS(:,52) = 0.1; CONSTANTS(:,53) = 0.012; CONSTANTS(:,54) = 1.053; CONSTANTS(:,55) = 1.053; CONSTANTS(:,56) = 0.8; CONSTANTS(:,57) = 1210.0; CONSTANTS(:,58) = 0.35; CONSTANTS(:,59) = 1.0; STATES(:,8) = 5.02792845976641014e-4; STATES(:,9) = 1.59632196638178009e-3; STATES(:,10) = 6.07316088556675987e-4; STATES(:,11) = 7.45773247314215036e-3; STATES(:,12) = 2.49592854373432000e-10; STATES(:,13) = 1.49041813757830998e-2; STATES(:,14) = 8.91799633266019011e-10; STATES(:,15) = 2.75455839709412009e-2; CONSTANTS(:,60) = 0.0; CONSTANTS(:,61) = 1.0; CONSTANTS(:,62) = 1.0; STATES(:,16) = 4.19991861054322011e-4; STATES(:,17) = 1.15141243826746994e-3; STATES(:,18) = 7.07824478944670999e-11; STATES(:,19) = 9.44463350378085993e-3; STATES(:,20) = 3.62113356111495976e-01; CONSTANTS(:,63) = 1.0; STATES(:,21) = 3.41341142614041016e-4; STATES(:,22) = 3.64315164237569004e-4; STATES(:,23) = 7.05656306851923029e-4; STATES(:,24) = 6.39038440072506948e-4; STATES(:,25) = 1.91017987408719017e-2; STATES(:,26) = 1.18055788874765002e-9; STATES(:,27) = 2.03016833596287999e-1; STATES(:,28) = 1.13428924662652000e-10; STATES(:,29) = 1.10248953370551007e-2; CONSTANTS(:,64) = 1.0; STATES(:,30) = 5.67249910261072965e-01; STATES(:,31) = 1.191; CONSTANTS(:,65) = 2.71430000000000008e-5; CONSTANTS(:,66) = 0.26714; CONSTANTS(:,67) = 1.0; CONSTANTS(:,68) = 1.124; CONSTANTS(:,69) = 0.0087; CONSTANTS(:,70) = 0.047; CONSTANTS(:,71) = 0.00087; CONSTANTS(:,72) = 0.05; CONSTANTS(:,73) = 10.0; CONSTANTS(:,74) = 0.8; STATES(:,32) = 6.73518785672381992e-1; CONSTANTS(:,75) = 0.10408; CONSTANTS(:,76) = 0.00238; CONSTANTS(:,77) = 0.052633; CONSTANTS(:,78) = 0.0005; CONSTANTS(:,79) = 0.07; STATES(:,33) = 0.013394; STATES(:,34) = 0.023413; STATES(:,35) = 6.8659; STATES(:,36) = 0.023413; CONSTANTS(:,80) = 0.05; CONSTANTS(:,81) = 0.25; CONSTANTS(:,82) = 0.0015; CONSTANTS(:,83) = 0.1; CONSTANTS(:,84) = 0.05; CONSTANTS(:,85) = 0.00068; STATES(:,37) = 0.0; STATES(:,38) = 0.0; STATES(:,39) = 0.0; STATES(:,40) = 0.0; STATES(:,41) = 0.0; STATES(:,42) = 0.0; CONSTANTS(:,86) = 100000.0; CONSTANTS(:,87) = 10000.0; STATES(:,43) = 0.0017546; STATES(:,44) = 3.47102959606005013e-1; STATES(:,45) = 4.74081735738210996e-1; STATES(:,46) = 9.62359241535767040e+00; CONSTANTS(:,88) = 96487.0; CONSTANTS(:,89) = 0.01; CONSTANTS(:,90) = 3.14159265358979312e+00; CONSTANTS(:,91) = 0.0011; STATES(:,47) = 20.273; STATES(:,48) = 20.273; CONSTANTS(:,92) = 1.77e-5; CONSTANTS(:,93) = 87.8251; CONSTANTS(:,94) = 2.46108000000000016e-5; STATES(:,49) = 6.8909; STATES(:,50) = 6.8909; CONSTANTS(:,95) = 0.2; CONSTANTS(:,96) = 0.02; CONSTANTS(:,97) = 75.0; CONSTANTS(:,98) = 1.8; CONSTANTS(:,99) = 100.0; CONSTANTS(:,100) = 5.4; CONSTANTS(:,101) = 140.0; STATES(:,51) = -87.491; CONSTANTS(:,102) = 1.995e-07; STATES(:,52) = 1.0; STATES(:,53) = 0.0; STATES(:,54) = 0.0; STATES(:,55) = 0.0; STATES(:,56) = 0.0; STATES(:,57) = 0.0; STATES(:,58) = 0.0; STATES(:,59) = 0.0; CONSTANTS(:,103) = 1.0; CONSTANTS(:,104) = 1e-06; STATES(:,60) = 1.0; STATES(:,61) = 0.0; STATES(:,62) = 0.0; STATES(:,63) = 0.0; STATES(:,64) = 0.0; STATES(:,65) = 0.0; STATES(:,66) = 0.0; STATES(:,67) = 0.0; CONSTANTS(:,105) = 1.0; CONSTANTS(:,106) = 1e-06; CONSTANTS(:,107) = 0.000225; CONSTANTS(:,108) = 9e-07; STATES(:,68) = 0.99604; CONSTANTS(:,109) = 0.4; CONSTANTS(:,110) = 8.0; STATES(:,69) = 1.23059999999999995e-08; CONSTANTS(:,111) = 0.025; STATES(:,70) = 9.37220000000000071e-16; STATES(:,71) = 1.65950000000000014e-17; STATES(:,72) = 1.12010000000000000e-15; STATES(:,73) = 1.66129999999999997e-18; STATES(:,74) = 0.91141; STATES(:,75) = 5.36959999999999959e-07; STATES(:,76) = 2.48610000000000001e-08; STATES(:,77) = 2.87759999999999997e-10; STATES(:,78) = 1.12169999999999999e-10; STATES(:,79) = 2.59670000000000003e-12; STATES(:,80) = 8.78740000000000076e-15; STATES(:,81) = 0.084395; STATES(:,82) = 0.0029306; STATES(:,83) = 4.52285000000000024e-05; STATES(:,84) = 2.61750000000000024e-07; STATES(:,85) = 0.0011424; STATES(:,86) = 7.93370000000000029e-05; STATES(:,87) = 1.83659999999999992e-06; STATES(:,88) = 1.41719999999999995e-08; CONSTANTS(:,112) = 0.0027304; STATES(:,89) = 0.95624; STATES(:,90) = 3.19559999999999984e-07; STATES(:,91) = 7.039e-09; STATES(:,92) = 3.87629999999999971e-11; STATES(:,93) = 5.02769999999999999e-11; STATES(:,94) = 5.53740000000000021e-13; STATES(:,95) = 2.96639999999999984e-15; STATES(:,96) = 0.042127; STATES(:,97) = 6.95969999999999965e-04; STATES(:,98) = 5.11010000000000042e-06; STATES(:,99) = 1.407e-08; STATES(:,100) = 9.02690000000000046e-04; STATES(:,101) = 2.98259999999999987e-05; STATES(:,102) = 3.285e-07; STATES(:,103) = 1.206e-09; CONSTANTS(:,113) = 0.0046171; STATES(:,104) = 7.65988420110534033e-04; CONSTANTS(:,114) = 0.01; CONSTANTS(:,115) = 9.97940000000000003e-05; CONSTANTS(:,116) = 1.11470000000000002e-04; CONSTANTS(:,117) = 0.0001; CONSTANTS(:,118) = 0.01; CONSTANTS(:,119) = 0.025; CONSTANTS(:,120) = 0.16305; CONSTANTS(:,121) = 1.0542; CONSTANTS(:,122) = 0.27623; CONSTANTS(:,123) = 0.002331; STATES(:,105) = 5.89379755147717982e-02; CONSTANTS(:,124) = 0.00384; CONSTANTS(:,125) = 0.069537; CONSTANTS(:,126) = 0.317; CONSTANTS(:,127) = 0.10988; CONSTANTS(:,128) = 7.8605; STATES(:,106) = 2.39479458960527997e-01; STATES(:,107) = 0.83805; STATES(:,108) = 0.99281; STATES(:,109) = 6.81269999999999988e-04; STATES(:,110) = 0.0068172; STATES(:,111) = 0.99709; STATES(:,112) = 0.90163; STATES(:,113) = 0.001236; STATES(:,114) = 0.99123; STATES(:,115) = 0.79472; CONSTANTS(:,129) = 0.01368; CONSTANTS(:,130) = 0.052811; CONSTANTS(:,131) = 3.25; CONSTANTS(:,132) = 3.7; CONSTANTS(:,133) = 4.9; CONSTANTS(:,134) = 3.2e-09; CONSTANTS(:,135) = 0.000125; CONSTANTS(:,136) = 0.0036; CONSTANTS(:,137) = 1.3; CONSTANTS(:,138) = 12.3; CONSTANTS(:,139) = 87.5; CONSTANTS(:,140) = 0.27; CONSTANTS(:,141) = 0.32; CONSTANTS(:,142) = 4.5; CONSTANTS(:,143) = 0.0011001; CONSTANTS(:,144) = 5.7392; STATES(:,116) = 1.26345311579565994e-01; CONSTANTS(:,145) = 1.4; CONSTANTS(:,146) = 0.015265; CONSTANTS(:,147) = 1.5; CONSTANTS(:,148) = 2.6; CONSTANTS(:,149) = 1.846; CONSTANTS(:,150) = 0.092455; STATES(:,117) = 0.36003; STATES(:,118) = 0.0007053; CONSTANTS(:,151) = 600.0; CONSTANTS(:,152) = 0.0575; CONSTANTS(:,153) = 0.0005; STATES(:,119) = 3.66750000000000000e-09; STATES(:,120) = 7.30739999999999981e-09; CONSTANTS(:,154) = 20.0; CONSTANTS(:,155) = 1.1; CONSTANTS(:,156) = 0.000175; CONSTANTS(:,157) = 0.0005; CONSTANTS(:,158) = 0.4975; STATES(:,121) = 1.76869999999999985e-05; STATES(:,122) = 1.0; STATES(:,123) = 0.99798; STATES(:,124) = 1.0; STATES(:,125) = 0.98747; CONSTANTS(:,159) = 9.88539999999999992e-04; CONSTANTS(:,160) = 0.00092; CONSTANTS(:,161) = 0.80737; STATES(:,126) = 5.92167467082830967e-01; CONSTANTS(:,162) = 0.004375; CONSTANTS(:,163) = 0.11348; CONSTANTS(:,164) = 0.48302; CONSTANTS(:,165) = 15.0; CONSTANTS(:,166) = 0.01833; STATES(:,127) = 145.62; CONSTANTS(:,167) = 0.52218; CONSTANTS(:,168) = 21.58; CONSTANTS(:,169) = 2.642; CONSTANTS(:,170) = 11.89; CONSTANTS(:,171) = 50.0; CONSTANTS(:,172) = 0.8; CONSTANTS(:,173) = 1.4; CONSTANTS(:,174) = 0.029268; STATES(:,128) = 2.36821659448036986e-02; STATES(:,129) = 1.28402905095187994e-02; STATES(:,130) = 6.37363047239019025e-03; STATES(:,131) = 9.17039986149184062e-03; STATES(:,132) = 4.29171113639321980e-05; CONSTANTS(:,175) = 3.0; CONSTANTS(:,176) = 0.16957; CONSTANTS(:,177) = 2.12570000000000006e-04; CONSTANTS(:,178) = 0.12481; CONSTANTS(:,179) = 0.125; CONSTANTS(:,180) = 0.2; CONSTANTS(:,181) = 0.35; CONSTANTS(:,182) = 1.167; CONSTANTS(:,183) = 0.7629; CONSTANTS(:,184) = 0.9024; CONSTANTS(:,185) = 0.0; CONSTANTS(:,186) = 20.0; CONSTANTS(:,187) = 2.5; CONSTANTS(:,188) = 4.0; CONSTANTS(:,189) = 0.0196; CONSTANTS(:,190) = 3.71; CONSTANTS(:,191) = 8314.0; CONSTANTS(:,192) = 310.0; CONSTANTS(:,193) = 0.5; CONSTANTS(:,194) = 0.0388; CONSTANTS(:,195) = 0.1; CONSTANTS(:,196) = 50.0; STATES(:,133) = 2.25475702283052998e-01; STATES(:,134) = 2.76490711096605019e-03; STATES(:,135) = 9.04820284659604013e-02; CONSTANTS(:,197) = 2.4984; CONSTANTS(:,198) = 11.359; CONSTANTS(:,199) = 0.3755; STATES(:,136) = 1.92819110624504991e-01; CONSTANTS(:,200) = 100.0; CONSTANTS(:,201) = 100.0; CONSTANTS(:,202) = 100.0; STATES(:,137) = 4.89063888619455989e-01; STATES(:,138) = 6.64997605558790977e-02; STATES(:,139) = 6.46928309115710060e-02; CONSTANTS(:,203) = 0.1088; CONSTANTS(:,204) = 0.4612; CONSTANTS(:,205) = 0.3755; STATES(:,140) = 1.26950532507959013e-01; STATES(:,141) = 8.17161796756963987e-01; STATES(:,142) = 1.74057375932567010e-01; STATES(:,143) = 2.05444874210056000e-01; STATES(:,144) = 2.49911886495889995e-01; CONSTANTS(:,206) = 0.001; CONSTANTS(:,207) = 1.95259999999999991e-05; CONSTANTS(:,208) = 0.001469; CONSTANTS(:,209) = 0.2; CONSTANTS(:,210) = 0.3; STATES(:,145) = 2.82662056977524001e-02; CONSTANTS(:,211) = 0.0035731; CONSTANTS(:,212) = 0.010145; CONSTANTS(:,213) = 0.5; CONSTANTS(:,214) = 100.0; CONSTANTS(:,215) = 1000.0; CONSTANTS(:,216) = 1.00000 - CONSTANTS(:,37); CONSTANTS(:,217) = (1.00000 - CONSTANTS(:,38)) - CONSTANTS(:,39); CONSTANTS(:,218) = CONSTANTS(:,58).*0.000983300; CONSTANTS(:,219) = CONSTANTS(:,58).*0.00133000; CONSTANTS(:,220) = CONSTANTS(:,56); CONSTANTS(:,221) = CONSTANTS(:,58).*0.00650000; CONSTANTS(:,222) = CONSTANTS(:,57); CONSTANTS(:,223) = ( (CONSTANTS(:,52)+CONSTANTS(:,60)).*(CONSTANTS(:,55)+CONSTANTS(:,60)))./CONSTANTS(:,55); CONSTANTS(:,224) = ( ( CONSTANTS(:,49).*CONSTANTS(:,54)).*(CONSTANTS(:,48)+CONSTANTS(:,60))).*(CONSTANTS(:,53)+CONSTANTS(:,60)); CONSTANTS(:,225) = ( (CONSTANTS(:,48)+CONSTANTS(:,60)).*(CONSTANTS(:,49)+CONSTANTS(:,60)))./CONSTANTS(:,49); CONSTANTS(:,226) = ( (CONSTANTS(:,52)+CONSTANTS(:,60)).*(CONSTANTS(:,55)+CONSTANTS(:,60)))./CONSTANTS(:,55); CONSTANTS(:,227) = ( ( CONSTANTS(:,49).*CONSTANTS(:,54)).*(CONSTANTS(:,48)+CONSTANTS(:,60))).*(CONSTANTS(:,53)+CONSTANTS(:,60)); CONSTANTS(:,228) = CONSTANTS(:,1)./(CONSTANTS(:,2)+CONSTANTS(:,1)); CONSTANTS(:,229) = CONSTANTS(:,70)+CONSTANTS(:,68); CONSTANTS(:,230) = 1.50000.*CONSTANTS(:,78); CONSTANTS(:,231) = CONSTANTS(:,71).*CONSTANTS(:,69); CONSTANTS(:,232) = CONSTANTS(:,87); CONSTANTS(:,233) = CONSTANTS(:,86); CONSTANTS(:,234) = CONSTANTS(:,86); CONSTANTS(:,235) = CONSTANTS(:,86); CONSTANTS(:,236) = 7.50000e-14.*1.00000e+06; CONSTANTS(:,237) = 5.00000e-15.*1.00000e+06; CONSTANTS(:,238) = 9.00000e-15.*1.00000e+06; CONSTANTS(:,239) = ( CONSTANTS(:,191).*CONSTANTS(:,192))./CONSTANTS(:,88); CONSTANTS(:,240) = 0.500000.*power((CONSTANTS(:,100)./5.40000), 1.0 ./ 2); CONSTANTS(:,241) = 0.0138542.*power((CONSTANTS(:,100)./5.40000), 1.0 ./ 2); CONSTANTS(:,242) = ( 2.15000.*8.25000).*1.10000; CONSTANTS(:,243) = CONSTANTS(:,100)./(CONSTANTS(:,100)+CONSTANTS(:,147)); CONSTANTS(:,244) = 0.666700.*4.75000; CONSTANTS(:,245) = CONSTANTS(:,167).*CONSTANTS(:,189); CONSTANTS(:,246) = 0.0100000./50.0000; CONSTANTS(:,247) = CONSTANTS(:,200).*CONSTANTS(:,197); CONSTANTS(:,248) = CONSTANTS(:,201).*CONSTANTS(:,198); CONSTANTS(:,249) = CONSTANTS(:,202).*CONSTANTS(:,199); CONSTANTS(:,250) = CONSTANTS(:,200); CONSTANTS(:,251) = CONSTANTS(:,201); CONSTANTS(:,252) = CONSTANTS(:,202); CONSTANTS(:,253) = CONSTANTS(:,197); CONSTANTS(:,254) = CONSTANTS(:,198); CONSTANTS(:,255) = CONSTANTS(:,199); CONSTANTS(:,256) = (CONSTANTS(:,210)./(1.00000 - CONSTANTS(:,210))).*CONSTANTS(:,206)+ CONSTANTS(:,210).*CONSTANTS(:,209); CONSTANTS(:,257) = - 80.0000; CONSTANTS(:,258) = 0.156290.*CONSTANTS(:,221); CONSTANTS(:,259) = ( ( ( CONSTANTS(:,47).*CONSTANTS(:,48)).*CONSTANTS(:,54)).*(CONSTANTS(:,53)+CONSTANTS(:,60))).*(CONSTANTS(:,49)+CONSTANTS(:,60)); CONSTANTS(:,260) = ( ( ( CONSTANTS(:,47).*CONSTANTS(:,48)).*CONSTANTS(:,54)).*(CONSTANTS(:,53)+CONSTANTS(:,60))).*(CONSTANTS(:,49)+CONSTANTS(:,60)); CONSTANTS(:,261) = 1.00000./(1.00000+0.350000); CONSTANTS(:,262) = CONSTANTS(:,71)+CONSTANTS(:,69); CONSTANTS(:,263) = ( CONSTANTS(:,70).*CONSTANTS(:,69)+ CONSTANTS(:,68).*CONSTANTS(:,71))+CONSTANTS(:,231); CONSTANTS(:,264) = 1.00000./CONSTANTS(:,239); CONSTANTS(:,265) = CONSTANTS(:,196).*CONSTANTS(:,246); CONSTANTS(:,266) = CONSTANTS(:,250).*CONSTANTS(:,203); CONSTANTS(:,267) = CONSTANTS(:,251).*CONSTANTS(:,204); CONSTANTS(:,268) = CONSTANTS(:,252).*CONSTANTS(:,205); CONSTANTS(:,269) = CONSTANTS(:,200); CONSTANTS(:,270) = CONSTANTS(:,201); CONSTANTS(:,271) = CONSTANTS(:,202); CONSTANTS(:,272) = ( ( ( CONSTANTS(:,51).*CONSTANTS(:,53)).*CONSTANTS(:,49)).*(CONSTANTS(:,48)+CONSTANTS(:,60))).*(CONSTANTS(:,54)+CONSTANTS(:,60)); CONSTANTS(:,273) = ( ( ( CONSTANTS(:,51).*CONSTANTS(:,53)).*CONSTANTS(:,49)).*(CONSTANTS(:,48)+CONSTANTS(:,60))).*(CONSTANTS(:,54)+CONSTANTS(:,60)); CONSTANTS(:,274) = 1.00000+(CONSTANTS(:,17) - CONSTANTS(:,18))./CONSTANTS(:,25); CONSTANTS(:,275) = CONSTANTS(:,229)+CONSTANTS(:,262); CONSTANTS(:,276) = CONSTANTS(:,88).*CONSTANTS(:,264); CONSTANTS(:,277) = CONSTANTS(:,269).*CONSTANTS(:,253); CONSTANTS(:,278) = CONSTANTS(:,270).*CONSTANTS(:,254); CONSTANTS(:,279) = CONSTANTS(:,271).*CONSTANTS(:,255); CONSTANTS(:,280) = ( ( ( ( CONSTANTS(:,47).*CONSTANTS(:,51)).*CONSTANTS(:,48)).*CONSTANTS(:,53)).*(CONSTANTS(:,49)+CONSTANTS(:,60))).*(CONSTANTS(:,54)+CONSTANTS(:,60)); CONSTANTS(:,281) = ( ( ( ( CONSTANTS(:,47).*CONSTANTS(:,51)).*CONSTANTS(:,48)).*CONSTANTS(:,53)).*(CONSTANTS(:,49)+CONSTANTS(:,60))).*(CONSTANTS(:,54)+CONSTANTS(:,60)); CONSTANTS(:,282) = (CONSTANTS(:,25)./2.00000).*(power((power(CONSTANTS(:,274), 2.00000)+( 4.00000.*CONSTANTS(:,18))./CONSTANTS(:,25)), 1.0 ./ 2) - CONSTANTS(:,274)); CONSTANTS(:,283) = (((CONSTANTS(:,17)+CONSTANTS(:,30))+CONSTANTS(:,21))+CONSTANTS(:,24)) - CONSTANTS(:,29); CONSTANTS(:,284) = (( CONSTANTS(:,17).*CONSTANTS(:,24)+ CONSTANTS(:,30).*CONSTANTS(:,21))+ CONSTANTS(:,21).*CONSTANTS(:,24)) - CONSTANTS(:,29).*(CONSTANTS(:,21)+CONSTANTS(:,24)); CONSTANTS(:,285) = ( CONSTANTS(:,29).*CONSTANTS(:,21)).*CONSTANTS(:,24); CONSTANTS(:,286) = ((( ( - CONSTANTS(:,285)./27.0000).*power(CONSTANTS(:,283), 3.00000) - ( ( ( CONSTANTS(:,283).*CONSTANTS(:,283)).*CONSTANTS(:,284)).*CONSTANTS(:,284))./108.000)+( ( CONSTANTS(:,283).*CONSTANTS(:,284)).*CONSTANTS(:,285))./6.00000)+power(CONSTANTS(:,284), 3.00000)./27.0000)+( CONSTANTS(:,285).*CONSTANTS(:,285))./4.00000; CONSTANTS(:,287) = piecewise({CONSTANTS(:,286)<0.00000, power( - CONSTANTS(:,286), 1.0 ./ 2) }, 0.00000); CONSTANTS(:,288) = ((piecewise({CONSTANTS(:,286)>0.00000, power(CONSTANTS(:,286), 1.0 ./ 2) }, 0.00000)+CONSTANTS(:,285)./2.00000)+( CONSTANTS(:,283).*CONSTANTS(:,284))./6.00000) - power(CONSTANTS(:,283), 3.00000)./27.0000; CONSTANTS(:,289) = atan(CONSTANTS(:,287)./CONSTANTS(:,288))./3.00000; CONSTANTS(:,290) = power( CONSTANTS(:,288).*CONSTANTS(:,288)+ CONSTANTS(:,287).*CONSTANTS(:,287), 1.00000./6.00000); CONSTANTS(:,291) = (CONSTANTS(:,284)./3.00000 - ( CONSTANTS(:,283).*CONSTANTS(:,283))./9.00000)./( CONSTANTS(:,290).*CONSTANTS(:,290)); CONSTANTS(:,292) = ( CONSTANTS(:,290).*cos(CONSTANTS(:,289))).*(1.00000 - CONSTANTS(:,291)) - CONSTANTS(:,283)./3.00000; CONSTANTS(:,293) = 1.00000+(CONSTANTS(:,30) - CONSTANTS(:,31))./CONSTANTS(:,26); CONSTANTS(:,294) = (CONSTANTS(:,26)./2.00000).*(power((power(CONSTANTS(:,293), 2.00000)+( 4.00000.*CONSTANTS(:,31))./CONSTANTS(:,26)), 1.0 ./ 2) - CONSTANTS(:,293)); CONSTANTS(:,295) = 0.850000.*0.0250000; CONSTANTS(:,296) = 224.000.*CONSTANTS(:,295); CONSTANTS(:,297) = 3.00000.*CONSTANTS(:,295); CONSTANTS(:,298) = 0.150000.*0.0250000; CONSTANTS(:,299) = (1.00000 - CONSTANTS(:,40)) - CONSTANTS(:,41); CONSTANTS(:,300) = 1.00000 - CONSTANTS(:,42); CONSTANTS(:,301) = ( ( 2.00000.*CONSTANTS(:,90)).*CONSTANTS(:,91)).*(CONSTANTS(:,91)+CONSTANTS(:,89)); CONSTANTS(:,302) = 2.00000.*CONSTANTS(:,301); CONSTANTS(:,303) = CONSTANTS(:,302)./CONSTANTS(:,88); CONSTANTS(:,304) = ( ( ( 1000.00.*CONSTANTS(:,90)).*CONSTANTS(:,91)).*CONSTANTS(:,91)).*CONSTANTS(:,89); CONSTANTS(:,305) = CONSTANTS(:,304).*0.00200000; CONSTANTS(:,306) = CONSTANTS(:,303)./( 2.00000.*CONSTANTS(:,305)); CONSTANTS(:,307) = 0.0200000.*CONSTANTS(:,304); CONSTANTS(:,308) = CONSTANTS(:,304).*0.678000; CONSTANTS(:,309) = CONSTANTS(:,303)./( 2.00000.*CONSTANTS(:,308)); CONSTANTS(:,310) = 0.0400000.*CONSTANTS(:,304); CONSTANTS(:,311) = CONSTANTS(:,304).*0.00480000; CONSTANTS(:,312) = CONSTANTS(:,304).*0.0552000; CONSTANTS(:,313) = CONSTANTS(:,311)./CONSTANTS(:,312); CONSTANTS(:,314) = CONSTANTS(:,312)./CONSTANTS(:,308); CONSTANTS(:,315) = CONSTANTS(:,304).*0.0200000; CONSTANTS(:,316) = CONSTANTS(:,315)./CONSTANTS(:,305); CONSTANTS(:,317) = - CONSTANTS(:,303)./CONSTANTS(:,315); CONSTANTS(:,318) = - CONSTANTS(:,311)./CONSTANTS(:,315); CONSTANTS(:,319) = CONSTANTS(:,315)./CONSTANTS(:,308); CONSTANTS(:,320) = CONSTANTS(:,304)./CONSTANTS(:,307); CONSTANTS(:,321) = ( ( CONSTANTS(:,12).*CONSTANTS(:,261)).*CONSTANTS(:,297)).*CONSTANTS(:,320); CONSTANTS(:,322) = CONSTANTS(:,304)./CONSTANTS(:,308); CONSTANTS(:,323) = ( CONSTANTS(:,40).*CONSTANTS(:,295)).*CONSTANTS(:,320); CONSTANTS(:,324) = ( CONSTANTS(:,42).*CONSTANTS(:,298)).*CONSTANTS(:,320); CONSTANTS(:,325) = ( ( (1.00000 - CONSTANTS(:,11)).*(1.00000 - CONSTANTS(:,261))).*CONSTANTS(:,297)).*CONSTANTS(:,322); CONSTANTS(:,326) = ( CONSTANTS(:,299).*CONSTANTS(:,295)).*CONSTANTS(:,322); CONSTANTS(:,327) = CONSTANTS(:,304)./CONSTANTS(:,310); CONSTANTS(:,328) = (1.00000 - CONSTANTS(:,194)) - CONSTANTS(:,195); CONSTANTS(:,329) = ( ( (1.00000 - CONSTANTS(:,12)).*CONSTANTS(:,261)).*CONSTANTS(:,297)).*CONSTANTS(:,322); CONSTANTS(:,330) = ( ( CONSTANTS(:,11).*(1.00000 - CONSTANTS(:,261))).*CONSTANTS(:,297)).*CONSTANTS(:,327); CONSTANTS(:,331) = ( CONSTANTS(:,41).*CONSTANTS(:,295)).*CONSTANTS(:,327); CONSTANTS(:,332) = ( CONSTANTS(:,300).*CONSTANTS(:,298)).*CONSTANTS(:,327); CONSTANTS(:,333) = CONSTANTS(:,303)./CONSTANTS(:,308); CONSTANTS(:,334) = (1.00000 - CONSTANTS(:,176)) - CONSTANTS(:,177); CONSTANTS(:,335) = CONSTANTS(:,176)+CONSTANTS(:,177); CONSTANTS(:,336) = ( CONSTANTS(:,328).*CONSTANTS(:,193)).*CONSTANTS(:,322); CONSTANTS(:,337) = 0.200000.*CONSTANTS(:,193); CONSTANTS(:,338) = CONSTANTS(:,315)./CONSTANTS(:,308); CONSTANTS(:,339) = power( CONSTANTS(:,185).*1.00000, CONSTANTS(:,182)); CONSTANTS(:,340) = 1.00000 - CONSTANTS(:,179); CONSTANTS(:,341) = CONSTANTS(:,194); CONSTANTS(:,342) = CONSTANTS(:,303)./CONSTANTS(:,315); CONSTANTS(:,343) = (1.00000 - CONSTANTS(:,339)./(CONSTANTS(:,168)+CONSTANTS(:,339))).*CONSTANTS(:,174); CONSTANTS(:,344) = CONSTANTS(:,179) - CONSTANTS(:,178); CONSTANTS(:,345) = ( CONSTANTS(:,341).*CONSTANTS(:,337)).*CONSTANTS(:,320); CONSTANTS(:,346) = CONSTANTS(:,195); CONSTANTS(:,347) = power(87.8251, 4.00000); CONSTANTS(:,348) = ( CONSTANTS(:,343).*CONSTANTS(:,176)).*CONSTANTS(:,320); CONSTANTS(:,349) = ( CONSTANTS(:,343).*CONSTANTS(:,334)).*CONSTANTS(:,322); CONSTANTS(:,350) = ( CONSTANTS(:,343).*CONSTANTS(:,177)).*CONSTANTS(:,327); CONSTANTS(:,351) = power( CONSTANTS(:,185).*1.00000, CONSTANTS(:,183)); CONSTANTS(:,352) = ( CONSTANTS(:,346).*CONSTANTS(:,337)).*CONSTANTS(:,327); CONSTANTS(:,353) = (1.00000 - CONSTANTS(:,341)) - CONSTANTS(:,346); CONSTANTS(:,354) = 1.00000+(CONSTANTS(:,119) - CONSTANTS(:,111))./CONSTANTS(:,117); CONSTANTS(:,355) = power( CONSTANTS(:,185).*1.00000, CONSTANTS(:,184)); CONSTANTS(:,356) = ( CONSTANTS(:,353).*CONSTANTS(:,337)).*CONSTANTS(:,322); CONSTANTS(:,357) = (CONSTANTS(:,117)./2.00000).*(power((power(CONSTANTS(:,354), 2.00000)+( 4.00000.*CONSTANTS(:,111))./CONSTANTS(:,117)), 1.0 ./ 2) - CONSTANTS(:,354)); CONSTANTS(:,358) = CONSTANTS(:,181)./(1.00000 - CONSTANTS(:,181)); CONSTANTS(:,359) = 1.00000+(CONSTANTS(:,119) - CONSTANTS(:,83))./CONSTANTS(:,114); CONSTANTS(:,360) = CONSTANTS(:,358).*( CONSTANTS(:,179).*(((1.00000+CONSTANTS(:,190)) - CONSTANTS(:,190).*CONSTANTS(:,335)) - CONSTANTS(:,180))+ CONSTANTS(:,335).*(CONSTANTS(:,180) - 1.00000))+ ( CONSTANTS(:,190).*CONSTANTS(:,179)).*(CONSTANTS(:,180) - CONSTANTS(:,335)); CONSTANTS(:,361) = (CONSTANTS(:,114)./2.00000).*(power((power(CONSTANTS(:,359), 2.00000)+( 4.00000.*CONSTANTS(:,83))./CONSTANTS(:,114)), 1.0 ./ 2) - CONSTANTS(:,359)); CONSTANTS(:,362) = CONSTANTS(:,179).*((1.00000+CONSTANTS(:,190))+ CONSTANTS(:,180).*(CONSTANTS(:,358) - CONSTANTS(:,190))) - CONSTANTS(:,180).*(1.00000+CONSTANTS(:,358)); CONSTANTS(:,363) = power(CONSTANTS(:,138), 3.00000); CONSTANTS(:,364) = (CONSTANTS(:,360)./CONSTANTS(:,362)).*CONSTANTS(:,174); CONSTANTS(:,365) = power(CONSTANTS(:,139), 3.00000); CONSTANTS(:,366) = ( (CONSTANTS(:,180) - CONSTANTS(:,335)).*CONSTANTS(:,174)+ CONSTANTS(:,180).*CONSTANTS(:,364))./( (1.00000+CONSTANTS(:,190)).*CONSTANTS(:,179) - CONSTANTS(:,180)); CONSTANTS(:,367) = power(CONSTANTS(:,101), 3.00000); CONSTANTS(:,368) = ( ( CONSTANTS(:,190).*CONSTANTS(:,179)).*CONSTANTS(:,366))./CONSTANTS(:,364); CONSTANTS(:,369) = (1.00000 - CONSTANTS(:,355)./(CONSTANTS(:,170)+CONSTANTS(:,355))).*CONSTANTS(:,366); CONSTANTS(:,370) = CONSTANTS(:,160) - 0.000170000; CONSTANTS(:,371) = ( CONSTANTS(:,369).*CONSTANTS(:,178)).*CONSTANTS(:,320); CONSTANTS(:,372) = ( CONSTANTS(:,369).*CONSTANTS(:,340)).*CONSTANTS(:,322); CONSTANTS(:,373) = ( CONSTANTS(:,369).*CONSTANTS(:,344)).*CONSTANTS(:,327); CONSTANTS(:,374) = 1.00000 - CONSTANTS(:,368); CONSTANTS(:,375) = CONSTANTS(:,160).*(1.00000 - 0.460000); CONSTANTS(:,376) = (1.00000 - CONSTANTS(:,351)./(CONSTANTS(:,169)+CONSTANTS(:,351))).*CONSTANTS(:,364); CONSTANTS(:,377) = CONSTANTS(:,375); CONSTANTS(:,378) = ( CONSTANTS(:,376).*CONSTANTS(:,368)).*CONSTANTS(:,320); CONSTANTS(:,379) = ( CONSTANTS(:,376).*CONSTANTS(:,374)).*CONSTANTS(:,322); CONSTANTS(:,380) = 3.25000.*CONSTANTS(:,162); CONSTANTS(:,381) = ( CONSTANTS(:,194).*CONSTANTS(:,193)).*CONSTANTS(:,320); CONSTANTS(:,382) = (((CONSTANTS(:,17)+CONSTANTS(:,30))+CONSTANTS(:,27))+CONSTANTS(:,28)) - CONSTANTS(:,381); CONSTANTS(:,383) = (( CONSTANTS(:,17).*CONSTANTS(:,28)+ CONSTANTS(:,30).*CONSTANTS(:,27))+ CONSTANTS(:,27).*CONSTANTS(:,28)) - CONSTANTS(:,381).*(CONSTANTS(:,27)+CONSTANTS(:,28)); CONSTANTS(:,384) = ( CONSTANTS(:,381).*CONSTANTS(:,27)).*CONSTANTS(:,28); CONSTANTS(:,385) = ((( ( - CONSTANTS(:,384)./27.0000).*power(CONSTANTS(:,382), 3.00000) - ( ( ( CONSTANTS(:,382).*CONSTANTS(:,382)).*CONSTANTS(:,383)).*CONSTANTS(:,383))./108.000)+( ( CONSTANTS(:,382).*CONSTANTS(:,383)).*CONSTANTS(:,384))./6.00000)+power(CONSTANTS(:,383), 3.00000)./27.0000)+( CONSTANTS(:,384).*CONSTANTS(:,384))./4.00000; CONSTANTS(:,386) = piecewise({CONSTANTS(:,385)<0.00000, power( - CONSTANTS(:,385), 1.0 ./ 2) }, 0.00000); CONSTANTS(:,387) = ((piecewise({CONSTANTS(:,385)>0.00000, power(CONSTANTS(:,385), 1.0 ./ 2) }, 0.00000)+CONSTANTS(:,384)./2.00000)+( CONSTANTS(:,382).*CONSTANTS(:,383))./6.00000) - power(CONSTANTS(:,382), 3.00000)./27.0000; CONSTANTS(:,388) = atan(CONSTANTS(:,386)./CONSTANTS(:,387))./3.00000; CONSTANTS(:,389) = power( CONSTANTS(:,387).*CONSTANTS(:,387)+ CONSTANTS(:,386).*CONSTANTS(:,386), 1.00000./6.00000); CONSTANTS(:,390) = (CONSTANTS(:,383)./3.00000 - ( CONSTANTS(:,382).*CONSTANTS(:,382))./9.00000)./( CONSTANTS(:,389).*CONSTANTS(:,389)); CONSTANTS(:,391) = ( CONSTANTS(:,389).*cos(CONSTANTS(:,388))).*(1.00000 - CONSTANTS(:,390)) - CONSTANTS(:,382)./3.00000; CONSTANTS(:,392) = ((CONSTANTS(:,17) - CONSTANTS(:,18))+CONSTANTS(:,282))./( (CONSTANTS(:,292)./CONSTANTS(:,21)+1.00000).*(CONSTANTS(:,391)./CONSTANTS(:,27)+1.00000)); CONSTANTS(:,393) = ( ( CONSTANTS(:,282).*CONSTANTS(:,392)).*CONSTANTS(:,391))./( CONSTANTS(:,25).*CONSTANTS(:,27)); CONSTANTS(:,394) = ( CONSTANTS(:,393).*CONSTANTS(:,292))./CONSTANTS(:,21); CONSTANTS(:,395) = ((CONSTANTS(:,30) - CONSTANTS(:,31))+CONSTANTS(:,294))./( (CONSTANTS(:,292)./CONSTANTS(:,24)+1.00000).*(CONSTANTS(:,391)./CONSTANTS(:,28)+1.00000)); CONSTANTS(:,396) = ( ( CONSTANTS(:,294).*CONSTANTS(:,395)).*CONSTANTS(:,391))./( CONSTANTS(:,26).*CONSTANTS(:,28)); CONSTANTS(:,397) = ( CONSTANTS(:,396).*CONSTANTS(:,292))./CONSTANTS(:,24); CONSTANTS(:,398) = 0.0269000+CONSTANTS(:,393)./CONSTANTS(:,18); CONSTANTS(:,399) = 0.0329000+CONSTANTS(:,396)./CONSTANTS(:,31); CONSTANTS(:,400) = ( CONSTANTS(:,195).*CONSTANTS(:,193)).*CONSTANTS(:,327); CONSTANTS(:,401) = 1.00000+(CONSTANTS(:,119) - CONSTANTS(:,400))./CONSTANTS(:,118); CONSTANTS(:,402) = (CONSTANTS(:,118)./2.00000).*(power((power(CONSTANTS(:,401), 2.00000)+( 4.00000.*CONSTANTS(:,400))./CONSTANTS(:,118)), 1.0 ./ 2) - CONSTANTS(:,401)); CONSTANTS(:,403) = ((CONSTANTS(:,119) - CONSTANTS(:,111))+CONSTANTS(:,357))./( (1.00000+CONSTANTS(:,361)./CONSTANTS(:,114)).*(1.00000+CONSTANTS(:,402)./CONSTANTS(:,118))); CONSTANTS(:,404) = ( ( CONSTANTS(:,357).*CONSTANTS(:,403)).*CONSTANTS(:,402))./( CONSTANTS(:,117).*CONSTANTS(:,118)); CONSTANTS(:,405) = 0.0306000+CONSTANTS(:,404)./CONSTANTS(:,111); CONSTANTS(:,406) = ( CONSTANTS(:,404).*CONSTANTS(:,361))./CONSTANTS(:,114); CONSTANTS(:,407) = - CONSTANTS(:,303)./CONSTANTS(:,308); CONSTANTS(:,408) = - CONSTANTS(:,303)./CONSTANTS(:,308); CONSTANTS(:,409) = CONSTANTS(:,315)./CONSTANTS(:,308); CONSTANTS(:,410) = ( 3.00000.*CONSTANTS(:,303))./CONSTANTS(:,315); 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(:,8) = 0.00100000.*( CONSTANTS(:,220).*STATES(:,9) - ( CONSTANTS(:,222).*STATES(:,1)).*STATES(:,8)); RATES(:,10) = 0.00100000.*( CONSTANTS(:,56).*STATES(:,2) - ( CONSTANTS(:,57).*STATES(:,11)).*STATES(:,10)); RATES(:,16) = 0.00100000.*( CONSTANTS(:,56).*STATES(:,3) - ( CONSTANTS(:,57).*STATES(:,17)).*STATES(:,16)); RATES(:,21) = 0.00100000.*( CONSTANTS(:,220).*STATES(:,22) - ( CONSTANTS(:,222).*STATES(:,23)).*STATES(:,21)); RATES(:,24) = 0.00100000.*( CONSTANTS(:,56).*STATES(:,4) - ( CONSTANTS(:,57).*STATES(:,25)).*STATES(:,24)); RATES(:,32) = 0.00100000.*(( ( CONSTANTS(:,75).*STATES(:,20)).*(1.00000 - STATES(:,32)))./(CONSTANTS(:,65)+ (1.00000 - STATES(:,32)).*1.00000) - ( ( CONSTANTS(:,77).*CONSTANTS(:,67)).*STATES(:,32))./(CONSTANTS(:,66)+ STATES(:,32).*1.00000)); RATES(:,105) = 0.00100000.*(( ( CONSTANTS(:,125).*STATES(:,30)).*(1.00000 - STATES(:,105)))./(CONSTANTS(:,122)+ (1.00000 - STATES(:,105)).*1.00000) - ( ( CONSTANTS(:,126).*CONSTANTS(:,83)).*STATES(:,105))./(CONSTANTS(:,123)+ STATES(:,105).*1.00000)); RATES(:,106) = 0.00100000.*(( ( CONSTANTS(:,129).*STATES(:,7)).*(1.00000 - STATES(:,106)))./(CONSTANTS(:,127)+ (1.00000 - STATES(:,106)).*1.00000) - ( ( CONSTANTS(:,130).*CONSTANTS(:,29)).*STATES(:,106))./(CONSTANTS(:,128)+ STATES(:,106).*1.00000)); RATES(:,116) = 0.00100000.*(( ( CONSTANTS(:,146).*STATES(:,7)).*(1.00000 - STATES(:,116)))./(CONSTANTS(:,143)+ (1.00000 - STATES(:,116)).*1.00000) - ( ( CONSTANTS(:,150).*CONSTANTS(:,29)).*STATES(:,116))./(CONSTANTS(:,144)+ STATES(:,116).*1.00000)); RATES(:,128) = 0.00100000.*( ( CONSTANTS(:,189).*STATES(:,7)).*(CONSTANTS(:,378) - STATES(:,128)) - CONSTANTS(:,245).*STATES(:,128)); RATES(:,129) = 0.00100000.*( ( CONSTANTS(:,189).*STATES(:,20)).*(CONSTANTS(:,379) - STATES(:,129)) - CONSTANTS(:,245).*STATES(:,129)); RATES(:,130) = 0.00100000.*( ( CONSTANTS(:,189).*STATES(:,7)).*(CONSTANTS(:,371) - STATES(:,130)) - CONSTANTS(:,245).*STATES(:,130)); RATES(:,131) = 0.00100000.*( ( CONSTANTS(:,189).*STATES(:,20)).*(CONSTANTS(:,372) - STATES(:,131)) - CONSTANTS(:,245).*STATES(:,131)); RATES(:,132) = 0.00100000.*( ( CONSTANTS(:,189).*STATES(:,30)).*(CONSTANTS(:,373) - STATES(:,132)) - CONSTANTS(:,245).*STATES(:,132)); RATES(:,133) = 0.00100000.*( CONSTANTS(:,202).*STATES(:,134) - ( CONSTANTS(:,249).*STATES(:,133)).*STATES(:,7)); RATES(:,134) = 0.00100000.*(( ( CONSTANTS(:,201).*STATES(:,135)).*STATES(:,44) - (CONSTANTS(:,248)+CONSTANTS(:,202)).*STATES(:,134))+ ( CONSTANTS(:,249).*STATES(:,133)).*STATES(:,7)); RATES(:,7) = 0.00100000.*((( CONSTANTS(:,202).*STATES(:,134) - ( CONSTANTS(:,249).*STATES(:,133)).*STATES(:,7))+ CONSTANTS(:,265).*STATES(:,136)) - ( CONSTANTS(:,196).*(CONSTANTS(:,345) - STATES(:,136))).*STATES(:,7)); RATES(:,136) = 0.00100000.*( ( CONSTANTS(:,196).*(CONSTANTS(:,345) - STATES(:,136))).*STATES(:,7) - CONSTANTS(:,265).*STATES(:,136)); RATES(:,137) = 0.00100000.*( CONSTANTS(:,252).*STATES(:,138) - ( CONSTANTS(:,268).*STATES(:,137)).*STATES(:,20)); RATES(:,138) = 0.00100000.*(( ( CONSTANTS(:,251).*STATES(:,139)).*STATES(:,45) - (CONSTANTS(:,267)+CONSTANTS(:,252)).*STATES(:,138))+ ( CONSTANTS(:,268).*STATES(:,137)).*STATES(:,20)); RATES(:,20) = 0.00100000.*((( CONSTANTS(:,252).*STATES(:,138) - ( CONSTANTS(:,268).*STATES(:,137)).*STATES(:,20))+ CONSTANTS(:,265).*STATES(:,140)) - ( CONSTANTS(:,196).*(CONSTANTS(:,356) - STATES(:,140))).*STATES(:,20)); RATES(:,140) = 0.00100000.*( ( CONSTANTS(:,196).*(CONSTANTS(:,356) - STATES(:,140))).*STATES(:,20) - CONSTANTS(:,265).*STATES(:,140)); RATES(:,141) = 0.00100000.*( CONSTANTS(:,271).*STATES(:,142) - ( CONSTANTS(:,279).*STATES(:,141)).*STATES(:,30)); RATES(:,142) = 0.00100000.*(( ( CONSTANTS(:,270).*STATES(:,143)).*STATES(:,46) - (CONSTANTS(:,278)+CONSTANTS(:,271)).*STATES(:,142))+ ( CONSTANTS(:,279).*STATES(:,141)).*STATES(:,30)); RATES(:,30) = 0.00100000.*((( CONSTANTS(:,271).*STATES(:,142) - ( CONSTANTS(:,279).*STATES(:,141)).*STATES(:,30))+ CONSTANTS(:,265).*STATES(:,144)) - ( CONSTANTS(:,196).*(CONSTANTS(:,352) - STATES(:,144))).*STATES(:,30)); RATES(:,144) = 0.00100000.*( ( CONSTANTS(:,196).*(CONSTANTS(:,352) - STATES(:,144))).*STATES(:,30) - CONSTANTS(:,265).*STATES(:,144)); ALGEBRAIC(:,4) = CONSTANTS(:,394) - STATES(:,5); RATES(:,5) = 0.00100000.*(( ( CONSTANTS(:,32).*STATES(:,7)).*ALGEBRAIC(:,4))./(CONSTANTS(:,19)+ALGEBRAIC(:,4)) - ( ( CONSTANTS(:,34).*CONSTANTS(:,29)).*STATES(:,5))./(CONSTANTS(:,22)+STATES(:,5))); ALGEBRAIC(:,5) = CONSTANTS(:,397) - STATES(:,6); RATES(:,6) = 0.00100000.*(( ( CONSTANTS(:,33).*STATES(:,7)).*ALGEBRAIC(:,5))./(CONSTANTS(:,20)+ALGEBRAIC(:,5)) - ( ( CONSTANTS(:,35).*CONSTANTS(:,29)).*STATES(:,6))./(CONSTANTS(:,23)+STATES(:,6))); ALGEBRAIC(:,16) = CONSTANTS(:,406) - STATES(:,104); RATES(:,104) = 0.00100000.*(( ( CONSTANTS(:,120).*STATES(:,30)).*ALGEBRAIC(:,16))./(CONSTANTS(:,115)+ALGEBRAIC(:,16)) - ( ( CONSTANTS(:,121).*CONSTANTS(:,83)).*STATES(:,104))./(CONSTANTS(:,116)+STATES(:,104))); ALGEBRAIC(:,26) = 1.00000./(1.00000+exp((STATES(:,51)+91.0000)./6.10000)); RATES(:,117) = (ALGEBRAIC(:,26) - STATES(:,117))./CONSTANTS(:,151); ALGEBRAIC(:,30) = CONSTANTS(:,256) - STATES(:,145); RATES(:,145) = 0.00100000.*(( ( CONSTANTS(:,212).*STATES(:,20)).*ALGEBRAIC(:,30))./(CONSTANTS(:,208)+ALGEBRAIC(:,30)) - ( ( CONSTANTS(:,211).*CONSTANTS(:,67)).*STATES(:,145))./(CONSTANTS(:,207)+STATES(:,145))); ALGEBRAIC(:,34) = (CONSTANTS(:,323) - STATES(:,13)) - STATES(:,12); RATES(:,13) = 0.00100000.*( ( CONSTANTS(:,221).*STATES(:,7)).*ALGEBRAIC(:,34) - CONSTANTS(:,258).*STATES(:,13)); ALGEBRAIC(:,36) = (CONSTANTS(:,326) - STATES(:,19)) - STATES(:,18); RATES(:,19) = 0.00100000.*( ( CONSTANTS(:,221).*STATES(:,20)).*ALGEBRAIC(:,36) - CONSTANTS(:,258).*STATES(:,19)); ALGEBRAIC(:,37) = (CONSTANTS(:,331) - STATES(:,27)) - STATES(:,26); RATES(:,27) = 0.00100000.*( ( CONSTANTS(:,221).*STATES(:,30)).*ALGEBRAIC(:,37) - CONSTANTS(:,258).*STATES(:,27)); ALGEBRAIC(:,12) = 0.0250000./(1.00000+exp((STATES(:,51)+58.0000)./5.00000)); ALGEBRAIC(:,40) = 0.200000./(1.00000+exp((STATES(:,51)+19.0000)./ - 9.00000)); RATES(:,68) = (ALGEBRAIC(:,12)./(ALGEBRAIC(:,12)+ALGEBRAIC(:,40)) - STATES(:,68))./CONSTANTS(:,110); ALGEBRAIC(:,13) = 1.00000./(( 0.000600000.*(STATES(:,51) - 1.73840))./(1.00000 - exp( - 0.136000.*(STATES(:,51) - 1.73840))) - ( 0.000300000.*(STATES(:,51)+38.3608))./(1.00000 - exp( 0.152200.*(STATES(:,51)+38.3608)))); ALGEBRAIC(:,41) = 1.00000./(1.00000+exp((STATES(:,51)+10.0850)./ - 4.25000)); RATES(:,69) = (ALGEBRAIC(:,41) - STATES(:,69))./ALGEBRAIC(:,13); ALGEBRAIC(:,14) = 0.00739900./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0311960))./ - 0.800190)); ALGEBRAIC(:,42) = 0.00569920./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0415200))./1.34890)); RATES(:,74) = ALGEBRAIC(:,42).*STATES(:,81) - STATES(:,74).*( 4.00000.*ALGEBRAIC(:,14)); ALGEBRAIC(:,15) = 0.00994150./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0448090))./ - 0.581720)); ALGEBRAIC(:,43) = 0.00332010./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0942170))./0.953640)); RATES(:,89) = ALGEBRAIC(:,43).*STATES(:,96) - STATES(:,89).*( 4.00000.*ALGEBRAIC(:,15)); ALGEBRAIC(:,17) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, 0.00000 }, 0.135000.*exp(((87.0000+STATES(:,51))+CONSTANTS(:,131))./ - 6.80000)); ALGEBRAIC(:,44) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, 1.00000./( 0.130000.*(1.00000+exp(((STATES(:,51)+CONSTANTS(:,131))+27.4034)./ - 11.1000))) }, 3.56000.*exp( 0.0790000.*((STATES(:,51)+CONSTANTS(:,131))+7.00000))+ 310000..*exp( 0.350000.*((STATES(:,51)+CONSTANTS(:,131))+7.00000))); RATES(:,107) = ALGEBRAIC(:,17).*(1.00000 - STATES(:,107)) - ALGEBRAIC(:,44).*STATES(:,107); ALGEBRAIC(:,18) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, 0.00000 }, ( ( - 127140..*exp( 0.244400.*(STATES(:,51)+CONSTANTS(:,131)))+ - 6.94800e-05.*exp( - 0.0439100.*(STATES(:,51)+CONSTANTS(:,131)))).*((STATES(:,51)+CONSTANTS(:,131))+37.7800))./(1.00000+exp( 0.311000.*((STATES(:,51)+CONSTANTS(:,131))+79.2300)))); ALGEBRAIC(:,45) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, ( 0.300000.*exp( - 2.53500e-07.*(STATES(:,51)+CONSTANTS(:,131))))./(1.00000+exp( - 0.100000.*((STATES(:,51)+CONSTANTS(:,131))+32.0000))) }, ( 0.121200.*exp( - 0.0105200.*(STATES(:,51)+CONSTANTS(:,131))))./(1.00000+exp( - 0.137800.*((STATES(:,51)+CONSTANTS(:,131))+40.1400)))); RATES(:,108) = ALGEBRAIC(:,18).*(1.00000 - STATES(:,108)) - ALGEBRAIC(:,45).*STATES(:,108); ALGEBRAIC(:,19) = ( 0.320000.*(STATES(:,51)+58.4729))./(1.00000 - exp( - 0.100000.*(STATES(:,51)+58.4729))); ALGEBRAIC(:,46) = 0.0800000.*exp((13.7299 - STATES(:,51))./11.0000); RATES(:,109) = ALGEBRAIC(:,19).*(1.00000 - STATES(:,109)) - ALGEBRAIC(:,46).*STATES(:,109); ALGEBRAIC(:,20) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, 0.135000.*exp((87.0000+STATES(:,51))./ - 6.80000)); ALGEBRAIC(:,47) = piecewise({STATES(:,51)>= - 40.0000, 1.00000./( 0.130000.*(1.00000+exp((STATES(:,51)+27.4034)./ - 11.1000))) }, 3.56000.*exp( 0.0790000.*(STATES(:,51)+7.00000))+ 310000..*exp( 0.350000.*(STATES(:,51)+7.00000))); RATES(:,110) = ALGEBRAIC(:,20).*(1.00000 - STATES(:,110)) - ALGEBRAIC(:,47).*STATES(:,110); ALGEBRAIC(:,21) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, ( ( - 127140..*exp( 0.244400.*STATES(:,51)) - 6.94800e-05.*exp( - 0.0439100.*STATES(:,51))).*(STATES(:,51)+37.7800))./(1.00000+exp( 0.311000.*(STATES(:,51)+79.2300)))); ALGEBRAIC(:,48) = piecewise({STATES(:,51)>= - 40.0000, ( 0.300000.*exp( - 2.53500e-07.*STATES(:,51)))./(1.00000+exp( - 0.100000.*(STATES(:,51)+32.0000))) }, ( 0.121200.*exp( - 0.0105200.*STATES(:,51)))./(1.00000+exp( - 0.137800.*(STATES(:,51)+40.1400)))); RATES(:,111) = ALGEBRAIC(:,21).*(1.00000 - STATES(:,111)) - ALGEBRAIC(:,48).*STATES(:,111); ALGEBRAIC(:,22) = ( 0.320000.*(STATES(:,51)+58.4729))./(1.00000 - exp( - 0.100000.*(STATES(:,51)+58.4729))); ALGEBRAIC(:,49) = 0.0800000.*exp((13.7299 - STATES(:,51))./11.0000); RATES(:,112) = ALGEBRAIC(:,22).*(1.00000 - STATES(:,112)) - ALGEBRAIC(:,49).*STATES(:,112); ALGEBRAIC(:,23) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, 0.135000.*exp(((87.0000+STATES(:,51))+CONSTANTS(:,133))./ - 6.80000)); ALGEBRAIC(:,50) = piecewise({STATES(:,51)>= - 40.0000, (1.00000./0.130000)./(1.00000+exp(((STATES(:,51)+CONSTANTS(:,133))+27.4034)./ - 11.1000)) }, 3.56000.*exp( 0.0790000.*((STATES(:,51)+CONSTANTS(:,133))+7.00000))+ 310000..*exp( 0.350000.*((STATES(:,51)+CONSTANTS(:,133))+7.00000))); RATES(:,113) = ALGEBRAIC(:,23).*(1.00000 - STATES(:,113)) - ALGEBRAIC(:,50).*STATES(:,113); ALGEBRAIC(:,24) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, ( ( - 127140..*exp( 0.244400.*(STATES(:,51)+CONSTANTS(:,133))) - 6.94800e-05.*exp( - 0.0439100.*(STATES(:,51)+CONSTANTS(:,133)))).*((STATES(:,51)+CONSTANTS(:,133))+37.7800))./(1.00000+exp( 0.311000.*((STATES(:,51)+CONSTANTS(:,133))+79.2300)))); ALGEBRAIC(:,51) = piecewise({STATES(:,51)>= - 40.0000, ( 0.300000.*exp( - 2.53500e-07.*(STATES(:,51)+CONSTANTS(:,133))))./(1.00000+exp( - 0.100000.*((STATES(:,51)+CONSTANTS(:,133))+32.0000))) }, ( 0.121200.*exp( - 0.0105200.*(STATES(:,51)+CONSTANTS(:,133))))./(1.00000+exp( - 0.137800.*((STATES(:,51)+CONSTANTS(:,133))+40.1400)))); RATES(:,114) = ALGEBRAIC(:,24).*(1.00000 - STATES(:,114)) - ALGEBRAIC(:,51).*STATES(:,114); ALGEBRAIC(:,25) = ( 0.320000.*((STATES(:,51)+CONSTANTS(:,132))+58.4729))./(1.00000 - exp( - 0.100000.*((STATES(:,51)+CONSTANTS(:,132))+58.4729))); ALGEBRAIC(:,52) = 0.0800000.*exp(((STATES(:,51)+CONSTANTS(:,132)) - 13.7299)./ - 11.0000); RATES(:,115) = ALGEBRAIC(:,25).*(1.00000 - STATES(:,115)) - ALGEBRAIC(:,52).*STATES(:,115); ALGEBRAIC(:,27) = ( 0.320000.*(STATES(:,51)+47.1300))./(1.00000 - exp( - 0.100000.*(STATES(:,51)+47.1300))); ALGEBRAIC(:,53) = 0.0800000.*exp(STATES(:,51)./ - 11.0000); RATES(:,118) = ALGEBRAIC(:,27).*(1.00000 - STATES(:,118)) - ALGEBRAIC(:,53).*STATES(:,118); ALGEBRAIC(:,58) = (CONSTANTS(:,324) - STATES(:,15)) - STATES(:,14); RATES(:,15) = 0.00100000.*( ( CONSTANTS(:,221).*STATES(:,7)).*ALGEBRAIC(:,58) - CONSTANTS(:,258).*STATES(:,15)); ALGEBRAIC(:,61) = (CONSTANTS(:,332) - STATES(:,29)) - STATES(:,28); RATES(:,29) = 0.00100000.*( ( CONSTANTS(:,221).*STATES(:,30)).*ALGEBRAIC(:,61) - CONSTANTS(:,258).*STATES(:,29)); ALGEBRAIC(:,28) = 1.00000./(1.00000+exp((STATES(:,51)+9.43700)./ - 7.13300)); ALGEBRAIC(:,54) = 1.00000./(1.00000+exp((STATES(:,51) - 18.4099)./ - 29.3814)); ALGEBRAIC(:,66) = 1.00000./(1.00000+exp((STATES(:,51)+100.000)./29.3814)); ALGEBRAIC(:,78) = 1.00000./(ALGEBRAIC(:,54)./1.20890+ 3.50000.*ALGEBRAIC(:,66)); RATES(:,121) = (ALGEBRAIC(:,28) - STATES(:,121))./ALGEBRAIC(:,78); ALGEBRAIC(:,67) = ((1.00000./(1.00000+exp((STATES(:,51)+19.0000)./ - 9.00000)))./0.500000)./9.79530; ALGEBRAIC(:,29) = 1.00000./(1.00000+exp((STATES(:,51)+58.0000)./5.00000)); ALGEBRAIC(:,79) = 0.0479600.*ALGEBRAIC(:,29); RATES(:,122) = ALGEBRAIC(:,79).*(1.00000 - STATES(:,122)) - ALGEBRAIC(:,67).*STATES(:,122); ALGEBRAIC(:,80) = 0.0214400.*ALGEBRAIC(:,29); RATES(:,123) = ALGEBRAIC(:,80).*(1.00000 - STATES(:,123)) - ALGEBRAIC(:,67).*STATES(:,123); ALGEBRAIC(:,55) = (1.00000./(1.00000+exp((STATES(:,51)+60.0000)./5.00000)))./250.000; ALGEBRAIC(:,81) = 2.46000.*ALGEBRAIC(:,55); RATES(:,124) = ALGEBRAIC(:,81).*(1.00000 - STATES(:,124)) - ALGEBRAIC(:,67).*STATES(:,124); ALGEBRAIC(:,82) = 0.560340.*ALGEBRAIC(:,55); RATES(:,125) = ALGEBRAIC(:,82).*(1.00000 - STATES(:,125)) - ALGEBRAIC(:,67).*STATES(:,125); ALGEBRAIC(:,64) = 0.0906540.*exp( ( - 0.111570.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,92) = 0.388390./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.150190))./ - 0.606930)); RATES(:,75) = ( ALGEBRAIC(:,42).*STATES(:,76)+ ALGEBRAIC(:,92).*STATES(:,86)) - STATES(:,75).*( 2.00000.*ALGEBRAIC(:,14)+ 2.00000.*ALGEBRAIC(:,64)); RATES(:,76) = ((( ( 2.00000.*ALGEBRAIC(:,14)).*STATES(:,75)+ ( 2.00000.*ALGEBRAIC(:,42)).*STATES(:,77))+ ( 2.00000.*ALGEBRAIC(:,92)).*STATES(:,87))+ ( 3.00000.*ALGEBRAIC(:,64)).*STATES(:,78)) - STATES(:,76).*(((ALGEBRAIC(:,14)+ALGEBRAIC(:,42))+ALGEBRAIC(:,92))+ 2.00000.*ALGEBRAIC(:,64)); RATES(:,77) = (( ALGEBRAIC(:,14).*STATES(:,76)+ ( 3.00000.*ALGEBRAIC(:,92)).*STATES(:,88))+ ( 3.00000.*ALGEBRAIC(:,64)).*STATES(:,79)) - STATES(:,77).*(( 2.00000.*ALGEBRAIC(:,42)+ 2.00000.*ALGEBRAIC(:,92))+ 2.00000.*ALGEBRAIC(:,64)); RATES(:,78) = ( ALGEBRAIC(:,42).*STATES(:,79)+ ALGEBRAIC(:,92).*STATES(:,76)) - STATES(:,78).*(ALGEBRAIC(:,14)+ 3.00000.*ALGEBRAIC(:,64)); RATES(:,79) = (( ALGEBRAIC(:,14).*STATES(:,78)+ ( 2.00000.*ALGEBRAIC(:,92)).*STATES(:,77))+ ( 4.00000.*ALGEBRAIC(:,64)).*STATES(:,80)) - STATES(:,79).*((ALGEBRAIC(:,42)+ALGEBRAIC(:,92))+ 3.00000.*ALGEBRAIC(:,64)); ALGEBRAIC(:,76) = 0.00311240+(0.0283300 - 0.00311240)./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.0516600))./1.55220)); RATES(:,80) = ( ALGEBRAIC(:,92).*STATES(:,79) - STATES(:,80).*( 4.00000.*ALGEBRAIC(:,64)+CONSTANTS(:,112)))+ ALGEBRAIC(:,76).*STATES(:,70); RATES(:,81) = (( ( 4.00000.*ALGEBRAIC(:,14)).*STATES(:,74)+ ( 2.00000.*ALGEBRAIC(:,42)).*STATES(:,82))+ ALGEBRAIC(:,64).*STATES(:,85)) - STATES(:,81).*(( 3.00000.*ALGEBRAIC(:,14)+ALGEBRAIC(:,42))+ALGEBRAIC(:,92)); RATES(:,82) = (( ( 3.00000.*ALGEBRAIC(:,14)).*STATES(:,81)+ ( 3.00000.*ALGEBRAIC(:,42)).*STATES(:,83))+ ALGEBRAIC(:,64).*STATES(:,86)) - STATES(:,82).*(( 2.00000.*ALGEBRAIC(:,14)+ 2.00000.*ALGEBRAIC(:,42))+ 2.00000.*ALGEBRAIC(:,92)); RATES(:,83) = (( ( 2.00000.*ALGEBRAIC(:,14)).*STATES(:,82)+ ( 4.00000.*ALGEBRAIC(:,42)).*STATES(:,84))+ ALGEBRAIC(:,64).*STATES(:,87)) - STATES(:,83).*((ALGEBRAIC(:,14)+ 3.00000.*ALGEBRAIC(:,42))+ 3.00000.*ALGEBRAIC(:,92)); RATES(:,84) = ( ALGEBRAIC(:,14).*STATES(:,83)+ ALGEBRAIC(:,64).*STATES(:,88)) - STATES(:,84).*( 4.00000.*ALGEBRAIC(:,42)+ 4.00000.*ALGEBRAIC(:,92)); RATES(:,85) = ( ALGEBRAIC(:,42).*STATES(:,86)+ ALGEBRAIC(:,92).*STATES(:,81)) - STATES(:,85).*( 3.00000.*ALGEBRAIC(:,14)+ALGEBRAIC(:,64)); RATES(:,86) = ((( ( 3.00000.*ALGEBRAIC(:,14)).*STATES(:,85)+ ( 2.00000.*ALGEBRAIC(:,42)).*STATES(:,87))+ ( 2.00000.*ALGEBRAIC(:,92)).*STATES(:,82))+ ( 2.00000.*ALGEBRAIC(:,64)).*STATES(:,75)) - STATES(:,86).*((( 2.00000.*ALGEBRAIC(:,14)+ALGEBRAIC(:,42))+ALGEBRAIC(:,92))+ALGEBRAIC(:,64)); RATES(:,87) = ((( ( 2.00000.*ALGEBRAIC(:,14)).*STATES(:,86)+ ( 3.00000.*ALGEBRAIC(:,42)).*STATES(:,88))+ ( 3.00000.*ALGEBRAIC(:,92)).*STATES(:,83))+ ( 2.00000.*ALGEBRAIC(:,64)).*STATES(:,76)) - STATES(:,87).*(((ALGEBRAIC(:,14)+ 2.00000.*ALGEBRAIC(:,42))+ 2.00000.*ALGEBRAIC(:,92))+ALGEBRAIC(:,64)); RATES(:,88) = (( ALGEBRAIC(:,14).*STATES(:,87)+ ( 4.00000.*ALGEBRAIC(:,92)).*STATES(:,84))+ ( 2.00000.*ALGEBRAIC(:,64)).*STATES(:,77)) - STATES(:,88).*(( 3.00000.*ALGEBRAIC(:,42)+ 3.00000.*ALGEBRAIC(:,92))+ALGEBRAIC(:,64)); ALGEBRAIC(:,65) = 0.0657000.*exp( ( - 0.118990.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,94) = 0.563560./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.179860))./ - 0.583810)); RATES(:,90) = ( ALGEBRAIC(:,43).*STATES(:,91)+ ALGEBRAIC(:,94).*STATES(:,101)) - STATES(:,90).*( 2.00000.*ALGEBRAIC(:,15)+ 2.00000.*ALGEBRAIC(:,65)); RATES(:,91) = ((( ( 2.00000.*ALGEBRAIC(:,15)).*STATES(:,90)+ ( 2.00000.*ALGEBRAIC(:,43)).*STATES(:,92))+ ( 2.00000.*ALGEBRAIC(:,94)).*STATES(:,102))+ ( 3.00000.*ALGEBRAIC(:,65)).*STATES(:,93)) - STATES(:,91).*(((ALGEBRAIC(:,15)+ALGEBRAIC(:,43))+ALGEBRAIC(:,94))+ 2.00000.*ALGEBRAIC(:,65)); RATES(:,92) = (( ALGEBRAIC(:,15).*STATES(:,91)+ ( 3.00000.*ALGEBRAIC(:,94)).*STATES(:,103))+ ( 3.00000.*ALGEBRAIC(:,65)).*STATES(:,94)) - STATES(:,92).*(( 2.00000.*ALGEBRAIC(:,43)+ 2.00000.*ALGEBRAIC(:,94))+ 2.00000.*ALGEBRAIC(:,65)); RATES(:,93) = ( ALGEBRAIC(:,43).*STATES(:,94)+ ALGEBRAIC(:,94).*STATES(:,91)) - STATES(:,93).*(ALGEBRAIC(:,15)+ 3.00000.*ALGEBRAIC(:,65)); RATES(:,94) = (( ALGEBRAIC(:,15).*STATES(:,93)+ ( 2.00000.*ALGEBRAIC(:,94)).*STATES(:,92))+ ( 4.00000.*ALGEBRAIC(:,65)).*STATES(:,95)) - STATES(:,94).*((ALGEBRAIC(:,43)+ALGEBRAIC(:,94))+ 3.00000.*ALGEBRAIC(:,65)); ALGEBRAIC(:,77) = 0.000385250+(0.0124060 - 0.000385250)./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.0641180))./0.779920)); RATES(:,95) = ( ALGEBRAIC(:,94).*STATES(:,94) - STATES(:,95).*( 4.00000.*ALGEBRAIC(:,65)+CONSTANTS(:,113)))+ ALGEBRAIC(:,77).*STATES(:,72); RATES(:,96) = (( ( 4.00000.*ALGEBRAIC(:,15)).*STATES(:,89)+ ( 2.00000.*ALGEBRAIC(:,43)).*STATES(:,97))+ ALGEBRAIC(:,65).*STATES(:,100)) - STATES(:,96).*(( 3.00000.*ALGEBRAIC(:,15)+ALGEBRAIC(:,43))+ALGEBRAIC(:,94)); RATES(:,97) = (( ( 3.00000.*ALGEBRAIC(:,15)).*STATES(:,96)+ ( 3.00000.*ALGEBRAIC(:,43)).*STATES(:,98))+ ALGEBRAIC(:,65).*STATES(:,101)) - STATES(:,97).*(( 2.00000.*ALGEBRAIC(:,15)+ 2.00000.*ALGEBRAIC(:,43))+ 2.00000.*ALGEBRAIC(:,94)); RATES(:,98) = (( ( 2.00000.*ALGEBRAIC(:,15)).*STATES(:,97)+ ( 4.00000.*ALGEBRAIC(:,43)).*STATES(:,99))+ ALGEBRAIC(:,65).*STATES(:,102)) - STATES(:,98).*((ALGEBRAIC(:,15)+ 3.00000.*ALGEBRAIC(:,43))+ 3.00000.*ALGEBRAIC(:,94)); RATES(:,99) = ( ALGEBRAIC(:,15).*STATES(:,98)+ ALGEBRAIC(:,65).*STATES(:,103)) - STATES(:,99).*( 4.00000.*ALGEBRAIC(:,43)+ 4.00000.*ALGEBRAIC(:,94)); RATES(:,100) = ( ALGEBRAIC(:,43).*STATES(:,101)+ ALGEBRAIC(:,94).*STATES(:,96)) - STATES(:,100).*( 3.00000.*ALGEBRAIC(:,15)+ALGEBRAIC(:,65)); RATES(:,101) = ((( ( 3.00000.*ALGEBRAIC(:,15)).*STATES(:,100)+ ( 2.00000.*ALGEBRAIC(:,43)).*STATES(:,102))+ ( 2.00000.*ALGEBRAIC(:,94)).*STATES(:,97))+ ( 2.00000.*ALGEBRAIC(:,65)).*STATES(:,90)) - STATES(:,101).*((( 2.00000.*ALGEBRAIC(:,15)+ALGEBRAIC(:,43))+ALGEBRAIC(:,94))+ALGEBRAIC(:,65)); RATES(:,102) = ((( ( 2.00000.*ALGEBRAIC(:,15)).*STATES(:,101)+ ( 3.00000.*ALGEBRAIC(:,43)).*STATES(:,103))+ ( 3.00000.*ALGEBRAIC(:,94)).*STATES(:,98))+ ( 2.00000.*ALGEBRAIC(:,65)).*STATES(:,91)) - STATES(:,102).*(((ALGEBRAIC(:,15)+ 2.00000.*ALGEBRAIC(:,43))+ 2.00000.*ALGEBRAIC(:,94))+ALGEBRAIC(:,65)); RATES(:,103) = (( ALGEBRAIC(:,15).*STATES(:,102)+ ( 4.00000.*ALGEBRAIC(:,94)).*STATES(:,99))+ ( 2.00000.*ALGEBRAIC(:,65)).*STATES(:,92)) - STATES(:,103).*(( 3.00000.*ALGEBRAIC(:,43)+ 3.00000.*ALGEBRAIC(:,94))+ALGEBRAIC(:,65)); ALGEBRAIC(:,93) = 0.000441980.*exp( ( - 1.20220.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,104) = 0.000401730.*exp( ( 0.000208730.*STATES(:,51)).*CONSTANTS(:,264)); RATES(:,70) = ( - (ALGEBRAIC(:,76)+ALGEBRAIC(:,104)).*STATES(:,70)+ ALGEBRAIC(:,93).*STATES(:,71))+ CONSTANTS(:,112).*STATES(:,80); RATES(:,71) = ALGEBRAIC(:,104).*STATES(:,70) - ALGEBRAIC(:,93).*STATES(:,71); ALGEBRAIC(:,95) = 0.000237300.*exp( ( - 1.97420.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,105) = 0.000226520.*exp( ( 0.000246900.*STATES(:,51)).*CONSTANTS(:,264)); RATES(:,72) = ( - (ALGEBRAIC(:,77)+ALGEBRAIC(:,105)).*STATES(:,72)+ ALGEBRAIC(:,95).*STATES(:,73))+ CONSTANTS(:,113).*STATES(:,95); RATES(:,73) = ALGEBRAIC(:,105).*STATES(:,72) - ALGEBRAIC(:,95).*STATES(:,73); ALGEBRAIC(:,6) = ( ( CONSTANTS(:,37).*CONSTANTS(:,36)).*CONSTANTS(:,320) - STATES(:,9)) - STATES(:,8); ALGEBRAIC(:,35) = ( ALGEBRAIC(:,6).*(CONSTANTS(:,60)+CONSTANTS(:,52)) - STATES(:,15).*(CONSTANTS(:,52)+CONSTANTS(:,60)))+ ( CONSTANTS(:,50).*CONSTANTS(:,52)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,55)); ALGEBRAIC(:,59) = ( - STATES(:,15).*CONSTANTS(:,50)).*CONSTANTS(:,52); ALGEBRAIC(:,71) = ( - ALGEBRAIC(:,35)+power(( ALGEBRAIC(:,35).*ALGEBRAIC(:,35) - ( 4.00000.*CONSTANTS(:,223)).*ALGEBRAIC(:,59)), 1.0 ./ 2))./( 2.00000.*CONSTANTS(:,223)); ALGEBRAIC(:,87) = ALGEBRAIC(:,6)./(1.00000+ (ALGEBRAIC(:,71)./CONSTANTS(:,50)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,52))); ALGEBRAIC(:,99) = ( ALGEBRAIC(:,71).*ALGEBRAIC(:,87))./CONSTANTS(:,50); ALGEBRAIC(:,109) = ( ALGEBRAIC(:,99).*CONSTANTS(:,60))./CONSTANTS(:,52); RATES(:,9) = 0.00100000.*(( CONSTANTS(:,45).*ALGEBRAIC(:,99)+ CONSTANTS(:,43).*ALGEBRAIC(:,109)) - CONSTANTS(:,220).*STATES(:,9)); RATES(:,1) = 0.00100000.*(( CONSTANTS(:,45).*ALGEBRAIC(:,99)+ CONSTANTS(:,43).*ALGEBRAIC(:,109)) - ( CONSTANTS(:,222).*STATES(:,1)).*STATES(:,8)); ALGEBRAIC(:,9) = ( ( CONSTANTS(:,216).*CONSTANTS(:,36)).*CONSTANTS(:,327) - STATES(:,22)) - STATES(:,21); ALGEBRAIC(:,38) = ( ALGEBRAIC(:,9).*(CONSTANTS(:,60)+CONSTANTS(:,52)) - STATES(:,29).*(CONSTANTS(:,52)+CONSTANTS(:,60)))+ ( CONSTANTS(:,50).*CONSTANTS(:,52)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,55)); ALGEBRAIC(:,62) = ( - STATES(:,29).*CONSTANTS(:,50)).*CONSTANTS(:,52); ALGEBRAIC(:,74) = ( - ALGEBRAIC(:,38)+power(( ALGEBRAIC(:,38).*ALGEBRAIC(:,38) - ( 4.00000.*CONSTANTS(:,226)).*ALGEBRAIC(:,62)), 1.0 ./ 2))./( 2.00000.*CONSTANTS(:,226)); ALGEBRAIC(:,90) = ALGEBRAIC(:,9)./(1.00000+ (ALGEBRAIC(:,74)./CONSTANTS(:,50)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,52))); ALGEBRAIC(:,102) = ( ALGEBRAIC(:,74).*ALGEBRAIC(:,90))./CONSTANTS(:,50); ALGEBRAIC(:,112) = (CONSTANTS(:,60)./CONSTANTS(:,52)).*ALGEBRAIC(:,102); RATES(:,22) = 0.00100000.*(( CONSTANTS(:,45).*ALGEBRAIC(:,102)+ CONSTANTS(:,43).*ALGEBRAIC(:,112)) - CONSTANTS(:,220).*STATES(:,22)); RATES(:,23) = 0.00100000.*(( CONSTANTS(:,45).*ALGEBRAIC(:,102)+ CONSTANTS(:,43).*ALGEBRAIC(:,112)) - ( CONSTANTS(:,222).*STATES(:,23)).*STATES(:,21)); ALGEBRAIC(:,8) = ( ( CONSTANTS(:,217).*CONSTANTS(:,296)).*CONSTANTS(:,322) - STATES(:,3)) - STATES(:,16); ALGEBRAIC(:,60) = ( ALGEBRAIC(:,8).*(CONSTANTS(:,48)+CONSTANTS(:,60)) - ALGEBRAIC(:,36).*(CONSTANTS(:,48)+CONSTANTS(:,60)))+ ( CONSTANTS(:,47).*CONSTANTS(:,48)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,49)); ALGEBRAIC(:,72) = ( - ALGEBRAIC(:,36).*CONSTANTS(:,48)).*CONSTANTS(:,47); ALGEBRAIC(:,88) = ( - ALGEBRAIC(:,60)+power(( ALGEBRAIC(:,60).*ALGEBRAIC(:,60) - ( 4.00000.*CONSTANTS(:,225)).*ALGEBRAIC(:,72)), 1.0 ./ 2))./( 2.00000.*CONSTANTS(:,225)); ALGEBRAIC(:,100) = ALGEBRAIC(:,8)./(1.00000+ (ALGEBRAIC(:,88)./CONSTANTS(:,47)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,48))); ALGEBRAIC(:,111) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,88)).*ALGEBRAIC(:,100))./( CONSTANTS(:,47).*CONSTANTS(:,48)); ALGEBRAIC(:,120) = ( ALGEBRAIC(:,100).*ALGEBRAIC(:,88))./CONSTANTS(:,47); RATES(:,3) = 0.00100000.*(( CONSTANTS(:,46).*ALGEBRAIC(:,120)+ CONSTANTS(:,44).*ALGEBRAIC(:,111)) - CONSTANTS(:,56).*STATES(:,3)); RATES(:,17) = 0.00100000.*(( CONSTANTS(:,46).*ALGEBRAIC(:,120)+ CONSTANTS(:,44).*ALGEBRAIC(:,111)) - ( CONSTANTS(:,57).*STATES(:,17)).*STATES(:,16)); ALGEBRAIC(:,119) = ( CONSTANTS(:,60).*ALGEBRAIC(:,88))./CONSTANTS(:,49); RATES(:,18) = 0.00100000.*( ( CONSTANTS(:,219).*CONSTANTS(:,62)).*(ALGEBRAIC(:,119)+ALGEBRAIC(:,111)) - CONSTANTS(:,218).*STATES(:,18)); ALGEBRAIC(:,131) = ((CONSTANTS(:,400) - STATES(:,143)) - STATES(:,142)) - STATES(:,141); RATES(:,143) = 0.00100000.*((( ( CONSTANTS(:,269).*ALGEBRAIC(:,131)).*STATES(:,46) - CONSTANTS(:,277).*STATES(:,143)) - ( CONSTANTS(:,270).*STATES(:,143)).*STATES(:,46))+ CONSTANTS(:,278).*STATES(:,142)); ALGEBRAIC(:,2) = power( STATES(:,4).*1.00000, CONSTANTS(:,13)); ALGEBRAIC(:,32) = CONSTANTS(:,7).*(CONSTANTS(:,9)+ALGEBRAIC(:,2)./(CONSTANTS(:,4)+ALGEBRAIC(:,2))); ALGEBRAIC(:,68) = ( ALGEBRAIC(:,32).*CONSTANTS(:,330)).*CONSTANTS(:,228); ALGEBRAIC(:,84) = ( CONSTANTS(:,237).*(STATES(:,44) - STATES(:,46)))./CONSTANTS(:,310); ALGEBRAIC(:,97) = ( CONSTANTS(:,238).*(STATES(:,46) - STATES(:,45)))./CONSTANTS(:,310); ALGEBRAIC(:,108) = ( CONSTANTS(:,350).*CONSTANTS(:,186))./(1.00000+CONSTANTS(:,171)./STATES(:,46)); ALGEBRAIC(:,117) = ( (CONSTANTS(:,373)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,132)).*CONSTANTS(:,188))./(1.00000+CONSTANTS(:,173)./STATES(:,46)); ALGEBRAIC(:,125) = ALGEBRAIC(:,108)+ALGEBRAIC(:,117); ALGEBRAIC(:,135) = (( ( - CONSTANTS(:,269).*ALGEBRAIC(:,131)).*STATES(:,46)+ CONSTANTS(:,277).*STATES(:,143)) - ( CONSTANTS(:,270).*STATES(:,143)).*STATES(:,46))+ CONSTANTS(:,278).*STATES(:,142); RATES(:,46) = 0.00100000.*((((ALGEBRAIC(:,135)+ALGEBRAIC(:,68)) - ALGEBRAIC(:,125))+ALGEBRAIC(:,84)) - ALGEBRAIC(:,97)); ALGEBRAIC(:,141) = ((CONSTANTS(:,381) - STATES(:,135)) - STATES(:,134)) - STATES(:,133); RATES(:,135) = 0.00100000.*((( ( CONSTANTS(:,200).*ALGEBRAIC(:,141)).*STATES(:,44) - CONSTANTS(:,247).*STATES(:,135)) - ( CONSTANTS(:,201).*STATES(:,135)).*STATES(:,44))+ CONSTANTS(:,248).*STATES(:,134)); ALGEBRAIC(:,3) = power( STATES(:,2).*1.00000, CONSTANTS(:,14)); ALGEBRAIC(:,33) = power( STATES(:,2).*1.00000, CONSTANTS(:,15)); ALGEBRAIC(:,57) = ( CONSTANTS(:,8).*(CONSTANTS(:,10)+ALGEBRAIC(:,3)./(CONSTANTS(:,5)+ALGEBRAIC(:,3)))).*(1.00000 - ( (1.00000 - ( CONSTANTS(:,16).*ALGEBRAIC(:,33))./(CONSTANTS(:,6)+ALGEBRAIC(:,33))).*STATES(:,1))./(CONSTANTS(:,3)+STATES(:,1))); ALGEBRAIC(:,83) = ( ALGEBRAIC(:,57).*CONSTANTS(:,321)).*CONSTANTS(:,228); ALGEBRAIC(:,96) = ( CONSTANTS(:,237).*(STATES(:,44) - STATES(:,46)))./CONSTANTS(:,307); ALGEBRAIC(:,106) = ( CONSTANTS(:,236).*(STATES(:,44) - STATES(:,45)))./CONSTANTS(:,307); ALGEBRAIC(:,115) = ( CONSTANTS(:,348).*CONSTANTS(:,186))./(1.00000+CONSTANTS(:,171)./STATES(:,44)); ALGEBRAIC(:,123) = ( (CONSTANTS(:,378)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,128)).*CONSTANTS(:,187))./(1.00000+CONSTANTS(:,172)./STATES(:,44)); ALGEBRAIC(:,129) = ( (CONSTANTS(:,371)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,130)).*CONSTANTS(:,188))./(1.00000+CONSTANTS(:,173)./STATES(:,44)); ALGEBRAIC(:,136) = (ALGEBRAIC(:,115)+ALGEBRAIC(:,123))+ALGEBRAIC(:,129); ALGEBRAIC(:,146) = (( ( - CONSTANTS(:,200).*ALGEBRAIC(:,141)).*STATES(:,44)+ CONSTANTS(:,247).*STATES(:,135)) - ( CONSTANTS(:,201).*STATES(:,135)).*STATES(:,44))+ CONSTANTS(:,248).*STATES(:,134); RATES(:,44) = 0.00100000.*((((ALGEBRAIC(:,146)+ALGEBRAIC(:,83)) - ALGEBRAIC(:,136)) - ALGEBRAIC(:,96)) - ALGEBRAIC(:,106)); ALGEBRAIC(:,151) = ((CONSTANTS(:,336) - STATES(:,139)) - STATES(:,138)) - STATES(:,137); RATES(:,139) = 0.00100000.*((( ( CONSTANTS(:,250).*ALGEBRAIC(:,151)).*STATES(:,45) - CONSTANTS(:,266).*STATES(:,139)) - ( CONSTANTS(:,251).*STATES(:,139)).*STATES(:,45))+ CONSTANTS(:,267).*STATES(:,138)); ALGEBRAIC(:,1) = power( STATES(:,3).*1.00000, CONSTANTS(:,13)); ALGEBRAIC(:,56) = CONSTANTS(:,7).*(CONSTANTS(:,9)+ALGEBRAIC(:,1)./(CONSTANTS(:,4)+ALGEBRAIC(:,1))); ALGEBRAIC(:,69) = ( ALGEBRAIC(:,56).*CONSTANTS(:,325)).*CONSTANTS(:,228); ALGEBRAIC(:,31) = power( STATES(:,3).*1.00000, CONSTANTS(:,14)); ALGEBRAIC(:,85) = CONSTANTS(:,8).*(CONSTANTS(:,10)+ALGEBRAIC(:,31)./(CONSTANTS(:,5)+ALGEBRAIC(:,31))); ALGEBRAIC(:,107) = ( ALGEBRAIC(:,85).*CONSTANTS(:,329)).*CONSTANTS(:,228); ALGEBRAIC(:,116) = ( CONSTANTS(:,236).*(STATES(:,44) - STATES(:,45)))./CONSTANTS(:,308); ALGEBRAIC(:,124) = ( CONSTANTS(:,238).*(STATES(:,46) - STATES(:,45)))./CONSTANTS(:,308); ALGEBRAIC(:,130) = ( CONSTANTS(:,349).*CONSTANTS(:,186))./(1.00000+CONSTANTS(:,171)./STATES(:,45)); ALGEBRAIC(:,137) = ( (CONSTANTS(:,379)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,129)).*CONSTANTS(:,187))./(1.00000+CONSTANTS(:,172)./STATES(:,45)); ALGEBRAIC(:,142) = ( (CONSTANTS(:,372)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,131)).*CONSTANTS(:,188))./(1.00000+CONSTANTS(:,173)./STATES(:,45)); ALGEBRAIC(:,147) = (ALGEBRAIC(:,130)+ALGEBRAIC(:,137))+ALGEBRAIC(:,142); ALGEBRAIC(:,155) = (( ( - CONSTANTS(:,250).*ALGEBRAIC(:,151)).*STATES(:,45)+ CONSTANTS(:,266).*STATES(:,139)) - ( CONSTANTS(:,251).*STATES(:,139)).*STATES(:,45))+ CONSTANTS(:,267).*STATES(:,138); RATES(:,45) = 0.00100000.*(((((ALGEBRAIC(:,155)+ALGEBRAIC(:,69))+ALGEBRAIC(:,107)) - ALGEBRAIC(:,147))+ALGEBRAIC(:,116))+ALGEBRAIC(:,124)); ALGEBRAIC(:,7) = ( ( CONSTANTS(:,38).*CONSTANTS(:,296)).*CONSTANTS(:,320) - STATES(:,2)) - STATES(:,10); ALGEBRAIC(:,70) = (((CONSTANTS(:,259)+CONSTANTS(:,272))./CONSTANTS(:,224)+ALGEBRAIC(:,34))+ALGEBRAIC(:,58)) - ALGEBRAIC(:,7); ALGEBRAIC(:,86) = (( CONSTANTS(:,272).*(ALGEBRAIC(:,34) - ALGEBRAIC(:,7))+ CONSTANTS(:,259).*(ALGEBRAIC(:,58) - ALGEBRAIC(:,7)))+CONSTANTS(:,280))./CONSTANTS(:,224); ALGEBRAIC(:,98) = ( ALGEBRAIC(:,7).*CONSTANTS(:,280))./CONSTANTS(:,224); ALGEBRAIC(:,110) = ((( ( - ALGEBRAIC(:,98)./27.0000).*power(ALGEBRAIC(:,70), 3.00000) - ( ( ( ALGEBRAIC(:,70).*ALGEBRAIC(:,70)).*ALGEBRAIC(:,86)).*ALGEBRAIC(:,86))./108.000)+( ( ALGEBRAIC(:,70).*ALGEBRAIC(:,86)).*ALGEBRAIC(:,98))./6.00000)+power(ALGEBRAIC(:,86), 3.00000)./27.0000)+( ALGEBRAIC(:,98).*ALGEBRAIC(:,98))./4.00000; ALGEBRAIC(:,118) = piecewise({ALGEBRAIC(:,110)<0.00000, power( - ALGEBRAIC(:,110), 1.0 ./ 2) }, 0.00000); ALGEBRAIC(:,126) = ((piecewise({ALGEBRAIC(:,110)>0.00000, power(ALGEBRAIC(:,110), 1.0 ./ 2) }, 0.00000)+ALGEBRAIC(:,98)./2.00000)+( ALGEBRAIC(:,70).*ALGEBRAIC(:,86))./6.00000) - power(ALGEBRAIC(:,70), 3.00000)./27.0000; ALGEBRAIC(:,132) = atan(ALGEBRAIC(:,118)./ALGEBRAIC(:,126))./3.00000; ALGEBRAIC(:,138) = power( ALGEBRAIC(:,126).*ALGEBRAIC(:,126)+ ALGEBRAIC(:,118).*ALGEBRAIC(:,118), 1.00000./6.00000); ALGEBRAIC(:,143) = (ALGEBRAIC(:,86)./3.00000 - ( ALGEBRAIC(:,70).*ALGEBRAIC(:,70))./9.00000)./( ALGEBRAIC(:,138).*ALGEBRAIC(:,138)); ALGEBRAIC(:,148) = ( ALGEBRAIC(:,138).* sin(ALGEBRAIC(:,132))).*(1.00000+ALGEBRAIC(:,143)); ALGEBRAIC(:,152) = ( ALGEBRAIC(:,138).*cos(ALGEBRAIC(:,132))).*(1.00000 - ALGEBRAIC(:,143)) - ALGEBRAIC(:,70)./3.00000; ALGEBRAIC(:,157) = power(( ALGEBRAIC(:,152).*ALGEBRAIC(:,152)+ ALGEBRAIC(:,148).*ALGEBRAIC(:,148)), 1.0 ./ 2); ALGEBRAIC(:,161) = ALGEBRAIC(:,34)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,49))+( ALGEBRAIC(:,157).*(CONSTANTS(:,48)+CONSTANTS(:,60)))./( CONSTANTS(:,47).*CONSTANTS(:,48))); ALGEBRAIC(:,165) = ( CONSTANTS(:,60).*ALGEBRAIC(:,161))./CONSTANTS(:,49); ALGEBRAIC(:,169) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,161)).*ALGEBRAIC(:,157))./( CONSTANTS(:,47).*CONSTANTS(:,48)); RATES(:,12) = 0.00100000.*( ( CONSTANTS(:,219).*CONSTANTS(:,59)).*(ALGEBRAIC(:,165)+ALGEBRAIC(:,169)) - CONSTANTS(:,218).*STATES(:,12)); ALGEBRAIC(:,10) = ( ( CONSTANTS(:,39).*CONSTANTS(:,296)).*CONSTANTS(:,327) - STATES(:,4)) - STATES(:,24); ALGEBRAIC(:,73) = (((CONSTANTS(:,260)+CONSTANTS(:,273))./CONSTANTS(:,227)+ALGEBRAIC(:,37))+ALGEBRAIC(:,61)) - ALGEBRAIC(:,10); ALGEBRAIC(:,89) = (( CONSTANTS(:,273).*(ALGEBRAIC(:,37) - ALGEBRAIC(:,10))+ CONSTANTS(:,260).*(ALGEBRAIC(:,61) - ALGEBRAIC(:,10)))+CONSTANTS(:,281))./CONSTANTS(:,227); ALGEBRAIC(:,101) = ( ALGEBRAIC(:,10).*CONSTANTS(:,281))./CONSTANTS(:,227); ALGEBRAIC(:,113) = ((( ( - ALGEBRAIC(:,101)./27.0000).*power(ALGEBRAIC(:,73), 3.00000) - ( ( ( ALGEBRAIC(:,73).*ALGEBRAIC(:,73)).*ALGEBRAIC(:,89)).*ALGEBRAIC(:,89))./108.000)+( ( ALGEBRAIC(:,73).*ALGEBRAIC(:,89)).*ALGEBRAIC(:,101))./6.00000)+power(ALGEBRAIC(:,89), 3.00000)./27.0000)+( ALGEBRAIC(:,101).*ALGEBRAIC(:,101))./4.00000; ALGEBRAIC(:,121) = piecewise({ALGEBRAIC(:,113)<0.00000, power( - ALGEBRAIC(:,113), 1.0 ./ 2) }, 0.00000); ALGEBRAIC(:,127) = ((piecewise({ALGEBRAIC(:,113)>0.00000, power(ALGEBRAIC(:,113), 1.0 ./ 2) }, 0.00000)+ALGEBRAIC(:,101)./2.00000)+( ALGEBRAIC(:,73).*ALGEBRAIC(:,89))./6.00000) - power(ALGEBRAIC(:,73), 3.00000)./27.0000; ALGEBRAIC(:,133) = atan(ALGEBRAIC(:,121)./ALGEBRAIC(:,127))./3.00000; ALGEBRAIC(:,139) = power( ALGEBRAIC(:,127).*ALGEBRAIC(:,127)+ ALGEBRAIC(:,121).*ALGEBRAIC(:,121), 1.00000./6.00000); ALGEBRAIC(:,144) = (ALGEBRAIC(:,89)./3.00000 - ( ALGEBRAIC(:,73).*ALGEBRAIC(:,73))./9.00000)./( ALGEBRAIC(:,139).*ALGEBRAIC(:,139)); ALGEBRAIC(:,149) = ( ALGEBRAIC(:,139).* sin(ALGEBRAIC(:,133))).*(1.00000+ALGEBRAIC(:,144)); ALGEBRAIC(:,153) = ( ALGEBRAIC(:,139).*cos(ALGEBRAIC(:,133))).*(1.00000 - ALGEBRAIC(:,144)) - ALGEBRAIC(:,73)./3.00000; ALGEBRAIC(:,158) = power(( ALGEBRAIC(:,153).*ALGEBRAIC(:,153)+ ALGEBRAIC(:,149).*ALGEBRAIC(:,149)), 1.0 ./ 2); ALGEBRAIC(:,162) = ALGEBRAIC(:,37)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,49))+( ALGEBRAIC(:,158).*(CONSTANTS(:,48)+CONSTANTS(:,60)))./( CONSTANTS(:,47).*CONSTANTS(:,48))); ALGEBRAIC(:,166) = ( CONSTANTS(:,60).*ALGEBRAIC(:,162))./CONSTANTS(:,49); ALGEBRAIC(:,170) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,162)).*ALGEBRAIC(:,158))./( CONSTANTS(:,47).*CONSTANTS(:,48)); RATES(:,26) = 0.00100000.*( ( CONSTANTS(:,219).*CONSTANTS(:,63)).*(ALGEBRAIC(:,166)+ALGEBRAIC(:,170)) - CONSTANTS(:,218).*STATES(:,26)); ALGEBRAIC(:,177) = ALGEBRAIC(:,58)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,54))+( ALGEBRAIC(:,157).*(CONSTANTS(:,53)+CONSTANTS(:,60)))./( CONSTANTS(:,51).*CONSTANTS(:,53))); ALGEBRAIC(:,181) = ( CONSTANTS(:,60).*ALGEBRAIC(:,177))./CONSTANTS(:,54); ALGEBRAIC(:,185) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,177)).*ALGEBRAIC(:,157))./( CONSTANTS(:,51).*CONSTANTS(:,53)); RATES(:,14) = 0.00100000.*( ( CONSTANTS(:,219).*CONSTANTS(:,59)).*(ALGEBRAIC(:,181)+ALGEBRAIC(:,185)) - CONSTANTS(:,218).*STATES(:,14)); ALGEBRAIC(:,178) = ALGEBRAIC(:,61)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,54))+( ALGEBRAIC(:,158).*(CONSTANTS(:,53)+CONSTANTS(:,60)))./( CONSTANTS(:,51).*CONSTANTS(:,53))); ALGEBRAIC(:,182) = ( CONSTANTS(:,60).*ALGEBRAIC(:,178))./CONSTANTS(:,54); ALGEBRAIC(:,186) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,178)).*ALGEBRAIC(:,158))./( CONSTANTS(:,51).*CONSTANTS(:,53)); RATES(:,28) = 0.00100000.*( ( CONSTANTS(:,219).*CONSTANTS(:,63)).*(ALGEBRAIC(:,182)+ALGEBRAIC(:,186)) - CONSTANTS(:,218).*STATES(:,28)); ALGEBRAIC(:,189) = ALGEBRAIC(:,169)+ CONSTANTS(:,61).*ALGEBRAIC(:,185); ALGEBRAIC(:,173) = ( ALGEBRAIC(:,161).*ALGEBRAIC(:,157))./CONSTANTS(:,47); ALGEBRAIC(:,193) = ( ALGEBRAIC(:,177).*ALGEBRAIC(:,157))./CONSTANTS(:,51); ALGEBRAIC(:,197) = ALGEBRAIC(:,173)+ CONSTANTS(:,61).*ALGEBRAIC(:,193); RATES(:,2) = 0.00100000.*(( CONSTANTS(:,46).*ALGEBRAIC(:,197)+ CONSTANTS(:,44).*ALGEBRAIC(:,189)) - CONSTANTS(:,56).*STATES(:,2)); RATES(:,11) = 0.00100000.*(( CONSTANTS(:,46).*ALGEBRAIC(:,197)+ CONSTANTS(:,44).*ALGEBRAIC(:,189)) - ( CONSTANTS(:,57).*STATES(:,11)).*STATES(:,10)); ALGEBRAIC(:,190) = ALGEBRAIC(:,170)+ CONSTANTS(:,64).*ALGEBRAIC(:,186); ALGEBRAIC(:,174) = ( ALGEBRAIC(:,162).*ALGEBRAIC(:,158))./CONSTANTS(:,47); ALGEBRAIC(:,194) = ( ALGEBRAIC(:,178).*ALGEBRAIC(:,158))./CONSTANTS(:,51); ALGEBRAIC(:,198) = ALGEBRAIC(:,174)+ CONSTANTS(:,64).*ALGEBRAIC(:,194); RATES(:,4) = 0.00100000.*(( CONSTANTS(:,46).*ALGEBRAIC(:,198)+ CONSTANTS(:,44).*ALGEBRAIC(:,190)) - CONSTANTS(:,56).*STATES(:,4)); RATES(:,25) = 0.00100000.*(( CONSTANTS(:,46).*ALGEBRAIC(:,198)+ CONSTANTS(:,44).*ALGEBRAIC(:,190)) - ( CONSTANTS(:,57).*STATES(:,25)).*STATES(:,24)); ALGEBRAIC(:,134) = CONSTANTS(:,275) - STATES(:,36); ALGEBRAIC(:,140) = CONSTANTS(:,263) - STATES(:,36).*CONSTANTS(:,262); ALGEBRAIC(:,145) = - CONSTANTS(:,231).*STATES(:,36); ALGEBRAIC(:,150) = - ALGEBRAIC(:,134)./3.00000+ ( (2.00000./3.00000).*power(( ALGEBRAIC(:,134).*ALGEBRAIC(:,134) - 3.00000.*ALGEBRAIC(:,140)), 1.0 ./ 2)).*cos(acos((( ( 9.00000.*ALGEBRAIC(:,134)).*ALGEBRAIC(:,140) - ( ( 2.00000.*ALGEBRAIC(:,134)).*ALGEBRAIC(:,134)).*ALGEBRAIC(:,134)) - 27.0000.*ALGEBRAIC(:,145))./( 2.00000.*power( ALGEBRAIC(:,134).*ALGEBRAIC(:,134) - 3.00000.*ALGEBRAIC(:,140), 1.50000)))./3.00000); ALGEBRAIC(:,191) = ( CONSTANTS(:,80).*(1.00000 - STATES(:,43)))./(1.00000+CONSTANTS(:,82)./ALGEBRAIC(:,150)); ALGEBRAIC(:,195) = ALGEBRAIC(:,191)+STATES(:,43); ALGEBRAIC(:,200) = ALGEBRAIC(:,195)./(ALGEBRAIC(:,195)+CONSTANTS(:,81)); RATES(:,37) = (ALGEBRAIC(:,200) - STATES(:,37))./CONSTANTS(:,232); RATES(:,38) = (ALGEBRAIC(:,200) - STATES(:,38))./CONSTANTS(:,233); RATES(:,39) = (ALGEBRAIC(:,200) - STATES(:,39))./CONSTANTS(:,234); RATES(:,40) = (ALGEBRAIC(:,200) - STATES(:,40))./CONSTANTS(:,235); RATES(:,41) = (ALGEBRAIC(:,200) - STATES(:,41))./CONSTANTS(:,86); ALGEBRAIC(:,201) = 1.00000./(1.00000+power(CONSTANTS(:,81)./ALGEBRAIC(:,195), 2.00000)); RATES(:,42) = (ALGEBRAIC(:,201) - STATES(:,42))./CONSTANTS(:,87); ALGEBRAIC(:,199) = (CONSTANTS(:,206) - CONSTANTS(:,209))+STATES(:,145); ALGEBRAIC(:,204) = 0.500000.*(power((power(ALGEBRAIC(:,199), 2.00000)+ ( 4.00000.*CONSTANTS(:,206)).*CONSTANTS(:,209)), 1.0 ./ 2) - ALGEBRAIC(:,199)); RATES(:,126) = 0.00100000.*(( ( CONSTANTS(:,163).*STATES(:,20)).*(1.00000 - STATES(:,126)))./(CONSTANTS(:,159)+ (1.00000 - STATES(:,126)).*1.00000) - ( ( CONSTANTS(:,164).*ALGEBRAIC(:,204)).*STATES(:,126))./(CONSTANTS(:,161)+ STATES(:,126).*1.00000)); ALGEBRAIC(:,207) = (CONSTANTS(:,29)./CONSTANTS(:,320)+CONSTANTS(:,83)./CONSTANTS(:,327))+ALGEBRAIC(:,204)./CONSTANTS(:,322); RATES(:,43) = ( CONSTANTS(:,84).*ALGEBRAIC(:,191)).*ALGEBRAIC(:,195) - ( CONSTANTS(:,85).*STATES(:,43)).*(0.100000+( 0.900000.*ALGEBRAIC(:,207))./0.137100); ALGEBRAIC(:,156) = 1.00000./( (1.00000+exp((4.79800+STATES(:,51))./ - 7.56990)).*(1.00000+exp((25.0000+STATES(:,51))./ - 5.00000))); ALGEBRAIC(:,196) = 0.590000+( 0.800000.*exp( 0.0520000.*(STATES(:,51)+13.0000)))./(1.00000+exp( 0.132000.*(STATES(:,51)+13.0000))); ALGEBRAIC(:,203) = ALGEBRAIC(:,156)./ALGEBRAIC(:,196); ALGEBRAIC(:,206) = (1.00000 - ALGEBRAIC(:,156))./ALGEBRAIC(:,196); ALGEBRAIC(:,180) = 1.00000+power(0.00200000./ALGEBRAIC(:,150), 4.00000); ALGEBRAIC(:,184) = 6.00000./ALGEBRAIC(:,180); ALGEBRAIC(:,164) = 1.00000./(1.00000+exp((29.9790+STATES(:,51))./3.17750)); ALGEBRAIC(:,172) = (0.100000+ALGEBRAIC(:,164))./1.10000; ALGEBRAIC(:,215) = 1.00000./( ( 70.0000.*(1.00000 - 0.500000.*STATES(:,37))).*(1.00000+exp((STATES(:,51)+49.1000)./10.3490))); ALGEBRAIC(:,217) = 1.00000./(1.00000+exp((STATES(:,51)+0.213000)./ - 10.8070)); ALGEBRAIC(:,222) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./38.4940); ALGEBRAIC(:,225) = ALGEBRAIC(:,172)./ALGEBRAIC(:,222); ALGEBRAIC(:,234) = (1.00000 - ALGEBRAIC(:,172))./ALGEBRAIC(:,222); RATES(:,52) = (( - ((ALGEBRAIC(:,203)+ALGEBRAIC(:,184))+ALGEBRAIC(:,234)).*STATES(:,52)+ ALGEBRAIC(:,206).*STATES(:,56))+ CONSTANTS(:,103).*STATES(:,55))+ ALGEBRAIC(:,225).*STATES(:,53); RATES(:,56) = (( - ((ALGEBRAIC(:,206)+ALGEBRAIC(:,184))+ALGEBRAIC(:,234)).*STATES(:,56)+ ALGEBRAIC(:,203).*STATES(:,52))+ CONSTANTS(:,103).*STATES(:,59))+ ALGEBRAIC(:,225).*STATES(:,57); ALGEBRAIC(:,168) = (0.000100000+ALGEBRAIC(:,164))./1.00010; ALGEBRAIC(:,160) = 0.100000.*STATES(:,37); ALGEBRAIC(:,176) = 1.00000+power(0.0100000./ALGEBRAIC(:,150), 10.0000); ALGEBRAIC(:,188) = (32.5000 - (18.0000 - ALGEBRAIC(:,160))./ALGEBRAIC(:,180)) - 10.0000./ALGEBRAIC(:,176); ALGEBRAIC(:,219) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./ALGEBRAIC(:,188)); ALGEBRAIC(:,228) = ALGEBRAIC(:,168)./ALGEBRAIC(:,219); ALGEBRAIC(:,240) = (1.00000 - ALGEBRAIC(:,168))./ALGEBRAIC(:,219); RATES(:,55) = (( - ((ALGEBRAIC(:,203)+CONSTANTS(:,103))+ALGEBRAIC(:,240)).*STATES(:,55)+ ALGEBRAIC(:,184).*STATES(:,52))+ ALGEBRAIC(:,206).*STATES(:,59))+ ALGEBRAIC(:,228).*STATES(:,54); RATES(:,59) = (( - ((ALGEBRAIC(:,206)+CONSTANTS(:,103))+ALGEBRAIC(:,240)).*STATES(:,59)+ ALGEBRAIC(:,184).*STATES(:,56))+ ALGEBRAIC(:,203).*STATES(:,55))+ ALGEBRAIC(:,228).*STATES(:,58); ALGEBRAIC(:,231) = piecewise({abs(ALGEBRAIC(:,228))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,228)); ALGEBRAIC(:,237) = piecewise({abs(ALGEBRAIC(:,234))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,234)); ALGEBRAIC(:,243) = ( CONSTANTS(:,104).*( ( ALGEBRAIC(:,225).*ALGEBRAIC(:,240)).*ALGEBRAIC(:,184)))./( ( ALGEBRAIC(:,237).*ALGEBRAIC(:,231)).*CONSTANTS(:,103)); RATES(:,53) = (( - ((ALGEBRAIC(:,203)+ALGEBRAIC(:,243))+ALGEBRAIC(:,225)).*STATES(:,53)+ ALGEBRAIC(:,234).*STATES(:,52))+ CONSTANTS(:,104).*STATES(:,54))+ ALGEBRAIC(:,206).*STATES(:,57); RATES(:,54) = (( - ((ALGEBRAIC(:,203)+CONSTANTS(:,104))+ALGEBRAIC(:,228)).*STATES(:,54)+ ALGEBRAIC(:,240).*STATES(:,55))+ ALGEBRAIC(:,243).*STATES(:,53))+ ALGEBRAIC(:,206).*STATES(:,58); RATES(:,57) = (( - ((ALGEBRAIC(:,206)+ALGEBRAIC(:,243))+ALGEBRAIC(:,225)).*STATES(:,57)+ ALGEBRAIC(:,234).*STATES(:,56))+ CONSTANTS(:,104).*STATES(:,58))+ ALGEBRAIC(:,203).*STATES(:,53); RATES(:,58) = (( - ((ALGEBRAIC(:,206)+CONSTANTS(:,104))+ALGEBRAIC(:,228)).*STATES(:,58)+ ALGEBRAIC(:,240).*STATES(:,59))+ ALGEBRAIC(:,243).*STATES(:,57))+ ALGEBRAIC(:,203).*STATES(:,54); ALGEBRAIC(:,63) = CONSTANTS(:,275) - STATES(:,34); ALGEBRAIC(:,75) = CONSTANTS(:,263) - STATES(:,34).*CONSTANTS(:,262); ALGEBRAIC(:,91) = - CONSTANTS(:,231).*STATES(:,34); ALGEBRAIC(:,103) = - ALGEBRAIC(:,63)./3.00000+ ( (2.00000./3.00000).*power(( ALGEBRAIC(:,63).*ALGEBRAIC(:,63) - 3.00000.*ALGEBRAIC(:,75)), 1.0 ./ 2)).*cos(acos((( ( 9.00000.*ALGEBRAIC(:,63)).*ALGEBRAIC(:,75) - ( ( 2.00000.*ALGEBRAIC(:,63)).*ALGEBRAIC(:,63)).*ALGEBRAIC(:,63)) - 27.0000.*ALGEBRAIC(:,91))./( 2.00000.*power( ALGEBRAIC(:,63).*ALGEBRAIC(:,63) - 3.00000.*ALGEBRAIC(:,75), 1.50000)))./3.00000); ALGEBRAIC(:,210) = (ALGEBRAIC(:,150) - ALGEBRAIC(:,103))./CONSTANTS(:,96); ALGEBRAIC(:,11) = (STATES(:,5)+CONSTANTS(:,393))./CONSTANTS(:,18); ALGEBRAIC(:,221) = (ALGEBRAIC(:,11) - CONSTANTS(:,398))./(0.927300 - CONSTANTS(:,398)); ALGEBRAIC(:,224) = piecewise({ALGEBRAIC(:,221)<0.00000, 0.00000 }, ALGEBRAIC(:,221)); ALGEBRAIC(:,227) = 0.000257900.*(1.00000+ 0.100000.*STATES(:,37)); ALGEBRAIC(:,230) = exp( ( 2.00000.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,233) = ( ( ( ( ALGEBRAIC(:,227).*4.00000).*STATES(:,51)).*CONSTANTS(:,276)).*( ALGEBRAIC(:,150).*ALGEBRAIC(:,230) - 0.341000.*CONSTANTS(:,98)))./(ALGEBRAIC(:,230) - 1.00000); ALGEBRAIC(:,236) = ALGEBRAIC(:,233).*(STATES(:,56)+STATES(:,59)); ALGEBRAIC(:,239) = 0.000155200.*(1.00000+ 0.400000.*STATES(:,37)); ALGEBRAIC(:,242) = exp( ( 2.00000.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,245) = ( ( ( ( ALGEBRAIC(:,239).*4.00000).*STATES(:,51)).*CONSTANTS(:,276)).*( ALGEBRAIC(:,150).*ALGEBRAIC(:,242) - 0.341000.*CONSTANTS(:,98)))./(ALGEBRAIC(:,242) - 1.00000); ALGEBRAIC(:,247) = ALGEBRAIC(:,245).*(STATES(:,64)+STATES(:,67)); ALGEBRAIC(:,249) = (1.00000 - ALGEBRAIC(:,224)).*ALGEBRAIC(:,247)+ ALGEBRAIC(:,224).*ALGEBRAIC(:,236); RATES(:,34) = - CONSTANTS(:,306).*ALGEBRAIC(:,249)+ CONSTANTS(:,316).*ALGEBRAIC(:,210); ALGEBRAIC(:,192) = 1.00000./( (1.00000+exp((13.5600 - STATES(:,51))./9.45000)).*(1.00000+exp((25.0000+STATES(:,51))./ - 5.00000))); ALGEBRAIC(:,209) = ALGEBRAIC(:,192)./ALGEBRAIC(:,196); ALGEBRAIC(:,211) = (1.00000 - ALGEBRAIC(:,192))./ALGEBRAIC(:,196); ALGEBRAIC(:,235) = 1.00000+power(0.00110000./ALGEBRAIC(:,150), 4.00000); ALGEBRAIC(:,238) = 14.9186./ALGEBRAIC(:,235); ALGEBRAIC(:,220) = 1.00000./(1.00000+exp((17.5000+STATES(:,51))./3.00000)); ALGEBRAIC(:,226) = (0.247400+ALGEBRAIC(:,220))./1.24740; ALGEBRAIC(:,229) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./26.5530); ALGEBRAIC(:,246) = ALGEBRAIC(:,226)./ALGEBRAIC(:,229); ALGEBRAIC(:,252) = (1.00000 - ALGEBRAIC(:,226))./ALGEBRAIC(:,229); RATES(:,60) = (( - ((ALGEBRAIC(:,209)+ALGEBRAIC(:,238))+ALGEBRAIC(:,252)).*STATES(:,60)+ ALGEBRAIC(:,211).*STATES(:,64))+ CONSTANTS(:,105).*STATES(:,63))+ ALGEBRAIC(:,246).*STATES(:,61); RATES(:,64) = (( - ((ALGEBRAIC(:,211)+ALGEBRAIC(:,238))+ALGEBRAIC(:,252)).*STATES(:,64)+ ALGEBRAIC(:,209).*STATES(:,60))+ CONSTANTS(:,105).*STATES(:,67))+ ALGEBRAIC(:,246).*STATES(:,65); ALGEBRAIC(:,223) = (0.00100000+ALGEBRAIC(:,220))./1.00100; ALGEBRAIC(:,213) = 5.00000.*STATES(:,37); ALGEBRAIC(:,232) = 1.00000+power(0.0120000./ALGEBRAIC(:,150), 10.0000); ALGEBRAIC(:,241) = (13.8250 - (6.38360 - ALGEBRAIC(:,213))./ALGEBRAIC(:,235)) - 3.36960./ALGEBRAIC(:,232); ALGEBRAIC(:,244) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./ALGEBRAIC(:,241)); ALGEBRAIC(:,248) = ALGEBRAIC(:,223)./ALGEBRAIC(:,244); ALGEBRAIC(:,258) = (1.00000 - ALGEBRAIC(:,223))./ALGEBRAIC(:,244); RATES(:,63) = (( - ((ALGEBRAIC(:,209)+CONSTANTS(:,105))+ALGEBRAIC(:,258)).*STATES(:,63)+ ALGEBRAIC(:,238).*STATES(:,60))+ ALGEBRAIC(:,211).*STATES(:,67))+ ALGEBRAIC(:,248).*STATES(:,62); RATES(:,67) = (( - ((ALGEBRAIC(:,211)+CONSTANTS(:,105))+ALGEBRAIC(:,258)).*STATES(:,67)+ ALGEBRAIC(:,238).*STATES(:,64))+ ALGEBRAIC(:,209).*STATES(:,63))+ ALGEBRAIC(:,248).*STATES(:,66); ALGEBRAIC(:,250) = piecewise({abs(ALGEBRAIC(:,248))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,248)); ALGEBRAIC(:,255) = piecewise({abs(ALGEBRAIC(:,252))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,252)); ALGEBRAIC(:,261) = ( CONSTANTS(:,106).*( ( ALGEBRAIC(:,246).*ALGEBRAIC(:,258)).*ALGEBRAIC(:,238)))./( ( ALGEBRAIC(:,255).*ALGEBRAIC(:,250)).*CONSTANTS(:,105)); RATES(:,61) = (( - ((ALGEBRAIC(:,209)+ALGEBRAIC(:,261))+ALGEBRAIC(:,246)).*STATES(:,61)+ ALGEBRAIC(:,252).*STATES(:,60))+ CONSTANTS(:,106).*STATES(:,62))+ ALGEBRAIC(:,211).*STATES(:,65); RATES(:,62) = (( - ((ALGEBRAIC(:,209)+CONSTANTS(:,106))+ALGEBRAIC(:,248)).*STATES(:,62)+ ALGEBRAIC(:,258).*STATES(:,63))+ ALGEBRAIC(:,261).*STATES(:,61))+ ALGEBRAIC(:,211).*STATES(:,66); RATES(:,65) = (( - ((ALGEBRAIC(:,211)+ALGEBRAIC(:,261))+ALGEBRAIC(:,246)).*STATES(:,65)+ ALGEBRAIC(:,252).*STATES(:,64))+ CONSTANTS(:,106).*STATES(:,66))+ ALGEBRAIC(:,209).*STATES(:,61); RATES(:,66) = (( - ((ALGEBRAIC(:,211)+CONSTANTS(:,106))+ALGEBRAIC(:,248)).*STATES(:,66)+ ALGEBRAIC(:,258).*STATES(:,67))+ ALGEBRAIC(:,261).*STATES(:,65))+ ALGEBRAIC(:,209).*STATES(:,62); ALGEBRAIC(:,253) = CONSTANTS(:,244).*(1.00000+ 2.00000.*STATES(:,42)); ALGEBRAIC(:,256) = 0.112500.*ALGEBRAIC(:,253); ALGEBRAIC(:,114) = (CONSTANTS(:,73)+CONSTANTS(:,74)) - STATES(:,35); ALGEBRAIC(:,122) = STATES(:,35).*CONSTANTS(:,74); ALGEBRAIC(:,128) = (power(( ALGEBRAIC(:,114).*ALGEBRAIC(:,114)+ 4.00000.*ALGEBRAIC(:,122)), 1.0 ./ 2) - ALGEBRAIC(:,114))./2.00000; ALGEBRAIC(:,270) = ( ALGEBRAIC(:,249).*1.00000)./(1.00000+power(1.00000./ALGEBRAIC(:,128), 8.00000)); ALGEBRAIC(:,272) = ALGEBRAIC(:,256).*ALGEBRAIC(:,270); ALGEBRAIC(:,264) = 1.00000+0.0123000./ALGEBRAIC(:,128); ALGEBRAIC(:,266) = ALGEBRAIC(:,253)./ALGEBRAIC(:,264); RATES(:,119) = - (ALGEBRAIC(:,272)+STATES(:,119))./ALGEBRAIC(:,266); ALGEBRAIC(:,259) = CONSTANTS(:,244).*(1.00000+ 0.00000.*STATES(:,42)); ALGEBRAIC(:,262) = 0.112500.*ALGEBRAIC(:,259); ALGEBRAIC(:,273) = ( 1.99250.*ALGEBRAIC(:,262)).*ALGEBRAIC(:,270); ALGEBRAIC(:,268) = ( 0.535700.*ALGEBRAIC(:,259))./ALGEBRAIC(:,264); RATES(:,120) = - (ALGEBRAIC(:,273)+STATES(:,120))./ALGEBRAIC(:,268); ALGEBRAIC(:,287) = 1.00000+power(CONSTANTS(:,135)./ALGEBRAIC(:,150), 2.00000); ALGEBRAIC(:,286) = exp( ( (CONSTANTS(:,140) - 1.00000).*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,288) = 1.00000+ CONSTANTS(:,141).*ALGEBRAIC(:,286); ALGEBRAIC(:,284) = power(STATES(:,50), 3.00000); ALGEBRAIC(:,289) = ( CONSTANTS(:,137).*ALGEBRAIC(:,284)+ CONSTANTS(:,365).*ALGEBRAIC(:,150))+ ( CONSTANTS(:,363).*CONSTANTS(:,98)).*(1.00000+ALGEBRAIC(:,150)./CONSTANTS(:,136)); ALGEBRAIC(:,290) = ( ( CONSTANTS(:,136).*CONSTANTS(:,367)).*(1.00000+ALGEBRAIC(:,284)./CONSTANTS(:,363))+ ALGEBRAIC(:,284).*CONSTANTS(:,98))+ CONSTANTS(:,367).*ALGEBRAIC(:,150); ALGEBRAIC(:,285) = exp( ( CONSTANTS(:,140).*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,291) = ( 0.200000.*CONSTANTS(:,142)).*( ( ALGEBRAIC(:,284).*CONSTANTS(:,98)).*ALGEBRAIC(:,285) - ( CONSTANTS(:,367).*ALGEBRAIC(:,150)).*ALGEBRAIC(:,286)); ALGEBRAIC(:,292) = ALGEBRAIC(:,291)./( ( ALGEBRAIC(:,287).*ALGEBRAIC(:,288)).*(ALGEBRAIC(:,289)+ALGEBRAIC(:,290))); ALGEBRAIC(:,208) = (STATES(:,50) - STATES(:,49))./CONSTANTS(:,95); RATES(:,50) = - ( CONSTANTS(:,410).*ALGEBRAIC(:,292)+ALGEBRAIC(:,208)); ALGEBRAIC(:,39) = (STATES(:,6)+CONSTANTS(:,396))./CONSTANTS(:,31); ALGEBRAIC(:,312) = (ALGEBRAIC(:,39) - CONSTANTS(:,399))./(0.958600 - CONSTANTS(:,399)); ALGEBRAIC(:,314) = piecewise({ALGEBRAIC(:,312)<0.00000, 0.00000 }, ALGEBRAIC(:,312)); ALGEBRAIC(:,318) = (1.00000 - ALGEBRAIC(:,314)).*STATES(:,119)+ ALGEBRAIC(:,314).*STATES(:,120); ALGEBRAIC(:,308) = ( CONSTANTS(:,156).*exp(ALGEBRAIC(:,128)./CONSTANTS(:,154))).*(ALGEBRAIC(:,128) - ALGEBRAIC(:,150)); ALGEBRAIC(:,310) = ( CONSTANTS(:,157).*exp(ALGEBRAIC(:,128)./CONSTANTS(:,155))).*(ALGEBRAIC(:,128) - ALGEBRAIC(:,150)); ALGEBRAIC(:,316) = (1.00000 - ALGEBRAIC(:,314)).*ALGEBRAIC(:,308)+ ALGEBRAIC(:,314).*ALGEBRAIC(:,310); ALGEBRAIC(:,320) = ALGEBRAIC(:,316)+ALGEBRAIC(:,318); ALGEBRAIC(:,212) = (STATES(:,31) - ALGEBRAIC(:,128))./CONSTANTS(:,97); RATES(:,35) = ALGEBRAIC(:,212) - ALGEBRAIC(:,320); ALGEBRAIC(:,154) = (STATES(:,32) - 0.673519)./(0.999180 - 0.673519); ALGEBRAIC(:,159) = piecewise({ALGEBRAIC(:,154)<0.00000, 0.00000 }, ALGEBRAIC(:,154)); ALGEBRAIC(:,163) = (1.00000 - ALGEBRAIC(:,159)).*CONSTANTS(:,78)+ ALGEBRAIC(:,159).*CONSTANTS(:,230); ALGEBRAIC(:,175) = ALGEBRAIC(:,163)+CONSTANTS(:,76); ALGEBRAIC(:,179) = ((ALGEBRAIC(:,175) - STATES(:,33))+CONSTANTS(:,72))+CONSTANTS(:,79); ALGEBRAIC(:,167) = ALGEBRAIC(:,163).*CONSTANTS(:,76); ALGEBRAIC(:,183) = ((ALGEBRAIC(:,167) - STATES(:,33).*ALGEBRAIC(:,175))+ CONSTANTS(:,79).*CONSTANTS(:,76))+ CONSTANTS(:,72).*ALGEBRAIC(:,163); ALGEBRAIC(:,171) = - ALGEBRAIC(:,167).*STATES(:,33); ALGEBRAIC(:,187) = - ALGEBRAIC(:,179)./3.00000+ ( (2.00000./3.00000).*power(( ALGEBRAIC(:,179).*ALGEBRAIC(:,179) - 3.00000.*ALGEBRAIC(:,183)), 1.0 ./ 2)).*cos(acos((( ( 9.00000.*ALGEBRAIC(:,179)).*ALGEBRAIC(:,183) - ( ( 2.00000.*ALGEBRAIC(:,179)).*ALGEBRAIC(:,179)).*ALGEBRAIC(:,179)) - 27.0000.*ALGEBRAIC(:,171))./( 2.00000.*power( ALGEBRAIC(:,179).*ALGEBRAIC(:,179) - 3.00000.*ALGEBRAIC(:,183), 1.50000)))./3.00000); ALGEBRAIC(:,202) = (ALGEBRAIC(:,150) - ALGEBRAIC(:,187))./CONSTANTS(:,95); RATES(:,36) = - ((( CONSTANTS(:,317).*ALGEBRAIC(:,292)+ CONSTANTS(:,318).*ALGEBRAIC(:,320))+ALGEBRAIC(:,202))+ALGEBRAIC(:,210)); ALGEBRAIC(:,251) = exp( STATES(:,51).*CONSTANTS(:,264)); ALGEBRAIC(:,254) = ( ( ( CONSTANTS(:,108).*STATES(:,51)).*CONSTANTS(:,276)).*(STATES(:,47) - CONSTANTS(:,99).*ALGEBRAIC(:,251)))./(1.00000 - ALGEBRAIC(:,251)); ALGEBRAIC(:,321) = 1.00000 - 1.00000./(1.00000+power(ALGEBRAIC(:,318)./CONSTANTS(:,109), 2.00000)); ALGEBRAIC(:,324) = ( ALGEBRAIC(:,254).*ALGEBRAIC(:,321)).*STATES(:,68); ALGEBRAIC(:,205) = (STATES(:,48) - STATES(:,47))./CONSTANTS(:,95); RATES(:,48) = CONSTANTS(:,342).*ALGEBRAIC(:,324) - ALGEBRAIC(:,205); ALGEBRAIC(:,307) = 1.00000./(1.00000+power(0.0300000./ALGEBRAIC(:,195), 2.00000)); ALGEBRAIC(:,309) = (1.00000 - ALGEBRAIC(:,307)).*CONSTANTS(:,162)+ ALGEBRAIC(:,307).*CONSTANTS(:,380); ALGEBRAIC(:,325) = ( ALGEBRAIC(:,309).*STATES(:,31))./CONSTANTS(:,165); ALGEBRAIC(:,311) = (STATES(:,126) - 0.659100)./(0.994500 - 0.659100); ALGEBRAIC(:,313) = piecewise({ALGEBRAIC(:,311)<0.00000, 0.00000 }, ALGEBRAIC(:,311)); ALGEBRAIC(:,315) = ALGEBRAIC(:,313).*STATES(:,41); ALGEBRAIC(:,317) = STATES(:,41) - ALGEBRAIC(:,315); ALGEBRAIC(:,319) = ALGEBRAIC(:,313) - ALGEBRAIC(:,315); ALGEBRAIC(:,322) = ((1.00000 - ALGEBRAIC(:,319)) - ALGEBRAIC(:,317)) - ALGEBRAIC(:,315); ALGEBRAIC(:,323) = (( ALGEBRAIC(:,322).*CONSTANTS(:,160)+ ALGEBRAIC(:,319).*CONSTANTS(:,375))+ ALGEBRAIC(:,317).*CONSTANTS(:,370))+ ALGEBRAIC(:,315).*CONSTANTS(:,377); ALGEBRAIC(:,327) = ( ALGEBRAIC(:,309).*ALGEBRAIC(:,187))./(ALGEBRAIC(:,187)+ALGEBRAIC(:,323)); ALGEBRAIC(:,329) = ALGEBRAIC(:,327) - ALGEBRAIC(:,325); RATES(:,31) = ALGEBRAIC(:,329) - CONSTANTS(:,313).*ALGEBRAIC(:,212); ALGEBRAIC(:,214) = ( 2.00000.*STATES(:,51)).*CONSTANTS(:,264); ALGEBRAIC(:,216) = exp(ALGEBRAIC(:,214)); ALGEBRAIC(:,218) = ( ( ( ( CONSTANTS(:,102).*2.00000).*CONSTANTS(:,88)).*ALGEBRAIC(:,214)).*( ALGEBRAIC(:,187).*ALGEBRAIC(:,216) - 0.341000.*CONSTANTS(:,98)))./(ALGEBRAIC(:,216) - 1.00000); ALGEBRAIC(:,293) = 1.00000+power(CONSTANTS(:,135)./ALGEBRAIC(:,187), 2.00000); ALGEBRAIC(:,294) = 1.00000+ CONSTANTS(:,141).*ALGEBRAIC(:,286); ALGEBRAIC(:,283) = power(STATES(:,49), 3.00000); ALGEBRAIC(:,295) = ( CONSTANTS(:,137).*ALGEBRAIC(:,283)+ CONSTANTS(:,365).*ALGEBRAIC(:,187))+ ( CONSTANTS(:,363).*CONSTANTS(:,98)).*(1.00000+ALGEBRAIC(:,187)./CONSTANTS(:,136)); ALGEBRAIC(:,296) = ( ( CONSTANTS(:,136).*CONSTANTS(:,367)).*(1.00000+ALGEBRAIC(:,283)./CONSTANTS(:,363))+ ALGEBRAIC(:,283).*CONSTANTS(:,98))+ CONSTANTS(:,367).*ALGEBRAIC(:,187); ALGEBRAIC(:,297) = ( 0.800000.*CONSTANTS(:,142)).*( ( ALGEBRAIC(:,283).*CONSTANTS(:,98)).*ALGEBRAIC(:,285) - ( CONSTANTS(:,367).*ALGEBRAIC(:,187)).*ALGEBRAIC(:,286)); ALGEBRAIC(:,298) = ALGEBRAIC(:,297)./( ( ALGEBRAIC(:,293).*ALGEBRAIC(:,294)).*(ALGEBRAIC(:,295)+ALGEBRAIC(:,296))); ALGEBRAIC(:,305) = ( CONSTANTS(:,152).*ALGEBRAIC(:,187))./(CONSTANTS(:,153)+ALGEBRAIC(:,187)); RATES(:,33) = ( - CONSTANTS(:,309).*((ALGEBRAIC(:,218)+ALGEBRAIC(:,305)) - 2.00000.*ALGEBRAIC(:,298)) - CONSTANTS(:,314).*ALGEBRAIC(:,329))+ CONSTANTS(:,319).*ALGEBRAIC(:,202); ALGEBRAIC(:,328) = - CONSTANTS(:,239).*log(CONSTANTS(:,99)./STATES(:,47)); ALGEBRAIC(:,332) = CONSTANTS(:,239).*log(CONSTANTS(:,100)./STATES(:,127)); ALGEBRAIC(:,333) = ALGEBRAIC(:,332) - ALGEBRAIC(:,328); ALGEBRAIC(:,334) = ( CONSTANTS(:,92).*ALGEBRAIC(:,333))./(ALGEBRAIC(:,333)+CONSTANTS(:,93)); ALGEBRAIC(:,353) = CONSTANTS(:,239).*log(CONSTANTS(:,101)./STATES(:,49)); ALGEBRAIC(:,354) = power(ALGEBRAIC(:,353) - ALGEBRAIC(:,328), 4.00000); ALGEBRAIC(:,355) = ( CONSTANTS(:,94).*ALGEBRAIC(:,354))./(ALGEBRAIC(:,354)+CONSTANTS(:,347)); ALGEBRAIC(:,330) = CONSTANTS(:,107).*(STATES(:,51) - ALGEBRAIC(:,328)); RATES(:,47) = (( CONSTANTS(:,333).*ALGEBRAIC(:,330)+ALGEBRAIC(:,355))+ALGEBRAIC(:,334))+ CONSTANTS(:,338).*ALGEBRAIC(:,205); ALGEBRAIC(:,356) = ( ( ( CONSTANTS(:,242).*power(STATES(:,115), 3.00000)).*STATES(:,107)).*STATES(:,108)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,357) = ( ( ( CONSTANTS(:,242).*power(STATES(:,109), 3.00000)).*STATES(:,107)).*STATES(:,108)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,358) = ( ( ( CONSTANTS(:,242).*power(STATES(:,112), 3.00000)).*STATES(:,110)).*STATES(:,111)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,359) = ( ( ( ( CONSTANTS(:,242).*power(STATES(:,115), 3.00000)).*STATES(:,113)).*STATES(:,114)).*(STATES(:,51) - ALGEBRAIC(:,353))).*1.25000; ALGEBRAIC(:,274) = (STATES(:,106) - 0.239480)./(0.950143 - 0.239480); ALGEBRAIC(:,275) = piecewise({ALGEBRAIC(:,274)<0.00000, 0.00000 }, ALGEBRAIC(:,274)); ALGEBRAIC(:,276) = ALGEBRAIC(:,275).*STATES(:,39); ALGEBRAIC(:,277) = STATES(:,39) - ALGEBRAIC(:,276); ALGEBRAIC(:,278) = ALGEBRAIC(:,275) - ALGEBRAIC(:,276); ALGEBRAIC(:,279) = ((1.00000 - ALGEBRAIC(:,278)) - ALGEBRAIC(:,277)) - ALGEBRAIC(:,276); ALGEBRAIC(:,360) = (( ALGEBRAIC(:,279).*ALGEBRAIC(:,358)+ ALGEBRAIC(:,278).*ALGEBRAIC(:,359))+ ALGEBRAIC(:,277).*ALGEBRAIC(:,357))+ ALGEBRAIC(:,276).*ALGEBRAIC(:,356); ALGEBRAIC(:,280) = STATES(:,51).*CONSTANTS(:,264); ALGEBRAIC(:,281) = exp(ALGEBRAIC(:,280)); ALGEBRAIC(:,282) = ( ( ( CONSTANTS(:,134).*CONSTANTS(:,88)).*ALGEBRAIC(:,280)).*( STATES(:,49).*ALGEBRAIC(:,281) - CONSTANTS(:,101)))./(ALGEBRAIC(:,281) - 1.00000); ALGEBRAIC(:,301) = ( CONSTANTS(:,145).*CONSTANTS(:,243))./(1.00000+exp( - (STATES(:,51)+92.0000).*CONSTANTS(:,264))); ALGEBRAIC(:,302) = ALGEBRAIC(:,301).*power(STATES(:,49)./(STATES(:,49)+CONSTANTS(:,148)), 3.00000); ALGEBRAIC(:,303) = ALGEBRAIC(:,301).*power(STATES(:,49)./(STATES(:,49)+CONSTANTS(:,149)), 3.00000); ALGEBRAIC(:,299) = (STATES(:,116) - 0.126345)./(0.998014 - 0.126345); ALGEBRAIC(:,300) = piecewise({ALGEBRAIC(:,299)<0.00000, 0.00000 }, ALGEBRAIC(:,299)); ALGEBRAIC(:,304) = (1.00000 - ALGEBRAIC(:,300)).*ALGEBRAIC(:,302)+ ALGEBRAIC(:,300).*ALGEBRAIC(:,303); ALGEBRAIC(:,361) = ( power(STATES(:,118), 3.00000).*STATES(:,117)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,362) = 0.0160000.*ALGEBRAIC(:,361); ALGEBRAIC(:,363) = 0.00650000.*ALGEBRAIC(:,361); ALGEBRAIC(:,364) = (1.00000 - STATES(:,39)).*ALGEBRAIC(:,363)+ STATES(:,39).*ALGEBRAIC(:,362); ALGEBRAIC(:,366) = (((ALGEBRAIC(:,360)+ALGEBRAIC(:,282))+ALGEBRAIC(:,364))+ ALGEBRAIC(:,304).*3.00000)+ ALGEBRAIC(:,298).*3.00000; RATES(:,49) = ( CONSTANTS(:,408).*ALGEBRAIC(:,366)+ CONSTANTS(:,409).*ALGEBRAIC(:,208))+ALGEBRAIC(:,355); ALGEBRAIC(:,306) = (((ALGEBRAIC(:,249)+ALGEBRAIC(:,218))+ALGEBRAIC(:,305)) - ALGEBRAIC(:,298).*2.00000) - ALGEBRAIC(:,292).*2.00000; ALGEBRAIC(:,331) = ALGEBRAIC(:,330)+ALGEBRAIC(:,324); ALGEBRAIC(:,335) = STATES(:,51) - ALGEBRAIC(:,332); ALGEBRAIC(:,336) = 1.02000./(1.00000+exp( 0.238500.*(ALGEBRAIC(:,335) - 59.2150))); ALGEBRAIC(:,337) = ( 0.491240.*exp( 0.0803200.*(ALGEBRAIC(:,335)+5.47600))+exp( 0.0617500.*(ALGEBRAIC(:,335) - 594.310)))./(1.00000+exp( - 0.514300.*(ALGEBRAIC(:,335)+4.75300))); ALGEBRAIC(:,338) = ( CONSTANTS(:,240).*(ALGEBRAIC(:,336)./(ALGEBRAIC(:,336)+ALGEBRAIC(:,337)))).*ALGEBRAIC(:,335); ALGEBRAIC(:,339) = ALGEBRAIC(:,338).*1.20000; ALGEBRAIC(:,340) = (1.00000 - STATES(:,38)).*ALGEBRAIC(:,338)+ STATES(:,38).*ALGEBRAIC(:,339); ALGEBRAIC(:,257) = 1.00000./(1.00000+exp((STATES(:,51)+10.0000)./15.4000)); ALGEBRAIC(:,341) = ( ( CONSTANTS(:,241).*STATES(:,69)).*ALGEBRAIC(:,257)).*(STATES(:,51) - ALGEBRAIC(:,332)); ALGEBRAIC(:,349) = CONSTANTS(:,239).*log((CONSTANTS(:,100)+ CONSTANTS(:,166).*CONSTANTS(:,101))./(STATES(:,127)+ CONSTANTS(:,166).*STATES(:,49))); ALGEBRAIC(:,260) = 0.195610.*(1.00000+0.600000./(1.00000+power(3.80000e-05./ALGEBRAIC(:,187), 1.40000))); ALGEBRAIC(:,350) = ( ALGEBRAIC(:,260).*(STATES(:,70)+STATES(:,71))).*(STATES(:,51) - ALGEBRAIC(:,349)); ALGEBRAIC(:,351) = ( ALGEBRAIC(:,260).*(STATES(:,72)+STATES(:,73))).*(STATES(:,51) - ALGEBRAIC(:,349)); ALGEBRAIC(:,263) = (STATES(:,104)+CONSTANTS(:,404))./CONSTANTS(:,111); ALGEBRAIC(:,265) = (ALGEBRAIC(:,263) - CONSTANTS(:,405))./(0.785000 - CONSTANTS(:,405)); ALGEBRAIC(:,267) = piecewise({ALGEBRAIC(:,265)<0.00000, 0.00000 }, ALGEBRAIC(:,265)); ALGEBRAIC(:,352) = ALGEBRAIC(:,267).*ALGEBRAIC(:,351)+ (1.00000 - ALGEBRAIC(:,267)).*ALGEBRAIC(:,350); ALGEBRAIC(:,342) = ( CONSTANTS(:,124).*(STATES(:,51) - ALGEBRAIC(:,332)))./(1.00000+exp((15.0000 - STATES(:,51))./17.0000)); ALGEBRAIC(:,343) = (( CONSTANTS(:,124).*(STATES(:,51) - ALGEBRAIC(:,332)))./(1.00000+exp((36.0000 - STATES(:,51))./17.0000))).*3.62000; ALGEBRAIC(:,269) = (STATES(:,105) - 0.0589380)./(0.393747 - 0.0589380); ALGEBRAIC(:,271) = piecewise({ALGEBRAIC(:,269)<0.00000, 0.00000 }, ALGEBRAIC(:,269)); ALGEBRAIC(:,344) = (1.00000 - ALGEBRAIC(:,271)).*ALGEBRAIC(:,342)+ ALGEBRAIC(:,271).*ALGEBRAIC(:,343); ALGEBRAIC(:,326) = exp(STATES(:,51)./550.000); ALGEBRAIC(:,345) = ( ( CONSTANTS(:,158).*power(STATES(:,121), 3.00000)).*ALGEBRAIC(:,326)).*(STATES(:,51) - ALGEBRAIC(:,332)); ALGEBRAIC(:,346) = ALGEBRAIC(:,345).*( 0.735600.*STATES(:,122)+ 0.264400.*STATES(:,124)); ALGEBRAIC(:,347) = ALGEBRAIC(:,345).*( 0.735600.*STATES(:,123)+ 0.264400.*STATES(:,125)); ALGEBRAIC(:,348) = (1.00000 - STATES(:,40)).*ALGEBRAIC(:,347)+ STATES(:,40).*ALGEBRAIC(:,346); ALGEBRAIC(:,365) = ((((ALGEBRAIC(:,340)+ALGEBRAIC(:,341))+ALGEBRAIC(:,352))+ALGEBRAIC(:,344))+ALGEBRAIC(:,348)) - 2.00000.*ALGEBRAIC(:,304); ALGEBRAIC(:,367) = ALGEBRAIC(:,366)+ ALGEBRAIC(:,292).*3.00000; ALGEBRAIC(:,368) = ((ALGEBRAIC(:,367)+ALGEBRAIC(:,365))+ALGEBRAIC(:,306))+ALGEBRAIC(:,331); ALGEBRAIC(:,369) = piecewise({(VOI - CONSTANTS(:,214)) - CONSTANTS(:,215).*floor((VOI - CONSTANTS(:,214))./CONSTANTS(:,215))<CONSTANTS(:,213), 1.00000 }, 0.00000).*CONSTANTS(:,257); RATES(:,51) = - (ALGEBRAIC(:,368)+ALGEBRAIC(:,369))./1.00000; RATES(:,127) = CONSTANTS(:,407).*(ALGEBRAIC(:,365)+ALGEBRAIC(:,369))+ALGEBRAIC(:,334); 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(:,4) = CONSTANTS(:,394) - STATES(:,5); ALGEBRAIC(:,5) = CONSTANTS(:,397) - STATES(:,6); ALGEBRAIC(:,16) = CONSTANTS(:,406) - STATES(:,104); ALGEBRAIC(:,26) = 1.00000./(1.00000+exp((STATES(:,51)+91.0000)./6.10000)); ALGEBRAIC(:,30) = CONSTANTS(:,256) - STATES(:,145); ALGEBRAIC(:,34) = (CONSTANTS(:,323) - STATES(:,13)) - STATES(:,12); ALGEBRAIC(:,36) = (CONSTANTS(:,326) - STATES(:,19)) - STATES(:,18); ALGEBRAIC(:,37) = (CONSTANTS(:,331) - STATES(:,27)) - STATES(:,26); ALGEBRAIC(:,12) = 0.0250000./(1.00000+exp((STATES(:,51)+58.0000)./5.00000)); ALGEBRAIC(:,40) = 0.200000./(1.00000+exp((STATES(:,51)+19.0000)./ - 9.00000)); ALGEBRAIC(:,13) = 1.00000./(( 0.000600000.*(STATES(:,51) - 1.73840))./(1.00000 - exp( - 0.136000.*(STATES(:,51) - 1.73840))) - ( 0.000300000.*(STATES(:,51)+38.3608))./(1.00000 - exp( 0.152200.*(STATES(:,51)+38.3608)))); ALGEBRAIC(:,41) = 1.00000./(1.00000+exp((STATES(:,51)+10.0850)./ - 4.25000)); ALGEBRAIC(:,14) = 0.00739900./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0311960))./ - 0.800190)); ALGEBRAIC(:,42) = 0.00569920./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0415200))./1.34890)); ALGEBRAIC(:,15) = 0.00994150./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0448090))./ - 0.581720)); ALGEBRAIC(:,43) = 0.00332010./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51) - 0.0942170))./0.953640)); ALGEBRAIC(:,17) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, 0.00000 }, 0.135000.*exp(((87.0000+STATES(:,51))+CONSTANTS(:,131))./ - 6.80000)); ALGEBRAIC(:,44) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, 1.00000./( 0.130000.*(1.00000+exp(((STATES(:,51)+CONSTANTS(:,131))+27.4034)./ - 11.1000))) }, 3.56000.*exp( 0.0790000.*((STATES(:,51)+CONSTANTS(:,131))+7.00000))+ 310000..*exp( 0.350000.*((STATES(:,51)+CONSTANTS(:,131))+7.00000))); ALGEBRAIC(:,18) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, 0.00000 }, ( ( - 127140..*exp( 0.244400.*(STATES(:,51)+CONSTANTS(:,131)))+ - 6.94800e-05.*exp( - 0.0439100.*(STATES(:,51)+CONSTANTS(:,131)))).*((STATES(:,51)+CONSTANTS(:,131))+37.7800))./(1.00000+exp( 0.311000.*((STATES(:,51)+CONSTANTS(:,131))+79.2300)))); ALGEBRAIC(:,45) = piecewise({STATES(:,51)+CONSTANTS(:,131)>= - 40.0000, ( 0.300000.*exp( - 2.53500e-07.*(STATES(:,51)+CONSTANTS(:,131))))./(1.00000+exp( - 0.100000.*((STATES(:,51)+CONSTANTS(:,131))+32.0000))) }, ( 0.121200.*exp( - 0.0105200.*(STATES(:,51)+CONSTANTS(:,131))))./(1.00000+exp( - 0.137800.*((STATES(:,51)+CONSTANTS(:,131))+40.1400)))); ALGEBRAIC(:,19) = ( 0.320000.*(STATES(:,51)+58.4729))./(1.00000 - exp( - 0.100000.*(STATES(:,51)+58.4729))); ALGEBRAIC(:,46) = 0.0800000.*exp((13.7299 - STATES(:,51))./11.0000); ALGEBRAIC(:,20) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, 0.135000.*exp((87.0000+STATES(:,51))./ - 6.80000)); ALGEBRAIC(:,47) = piecewise({STATES(:,51)>= - 40.0000, 1.00000./( 0.130000.*(1.00000+exp((STATES(:,51)+27.4034)./ - 11.1000))) }, 3.56000.*exp( 0.0790000.*(STATES(:,51)+7.00000))+ 310000..*exp( 0.350000.*(STATES(:,51)+7.00000))); ALGEBRAIC(:,21) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, ( ( - 127140..*exp( 0.244400.*STATES(:,51)) - 6.94800e-05.*exp( - 0.0439100.*STATES(:,51))).*(STATES(:,51)+37.7800))./(1.00000+exp( 0.311000.*(STATES(:,51)+79.2300)))); ALGEBRAIC(:,48) = piecewise({STATES(:,51)>= - 40.0000, ( 0.300000.*exp( - 2.53500e-07.*STATES(:,51)))./(1.00000+exp( - 0.100000.*(STATES(:,51)+32.0000))) }, ( 0.121200.*exp( - 0.0105200.*STATES(:,51)))./(1.00000+exp( - 0.137800.*(STATES(:,51)+40.1400)))); ALGEBRAIC(:,22) = ( 0.320000.*(STATES(:,51)+58.4729))./(1.00000 - exp( - 0.100000.*(STATES(:,51)+58.4729))); ALGEBRAIC(:,49) = 0.0800000.*exp((13.7299 - STATES(:,51))./11.0000); ALGEBRAIC(:,23) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, 0.135000.*exp(((87.0000+STATES(:,51))+CONSTANTS(:,133))./ - 6.80000)); ALGEBRAIC(:,50) = piecewise({STATES(:,51)>= - 40.0000, (1.00000./0.130000)./(1.00000+exp(((STATES(:,51)+CONSTANTS(:,133))+27.4034)./ - 11.1000)) }, 3.56000.*exp( 0.0790000.*((STATES(:,51)+CONSTANTS(:,133))+7.00000))+ 310000..*exp( 0.350000.*((STATES(:,51)+CONSTANTS(:,133))+7.00000))); ALGEBRAIC(:,24) = piecewise({STATES(:,51)>= - 40.0000, 0.00000 }, ( ( - 127140..*exp( 0.244400.*(STATES(:,51)+CONSTANTS(:,133))) - 6.94800e-05.*exp( - 0.0439100.*(STATES(:,51)+CONSTANTS(:,133)))).*((STATES(:,51)+CONSTANTS(:,133))+37.7800))./(1.00000+exp( 0.311000.*((STATES(:,51)+CONSTANTS(:,133))+79.2300)))); ALGEBRAIC(:,51) = piecewise({STATES(:,51)>= - 40.0000, ( 0.300000.*exp( - 2.53500e-07.*(STATES(:,51)+CONSTANTS(:,133))))./(1.00000+exp( - 0.100000.*((STATES(:,51)+CONSTANTS(:,133))+32.0000))) }, ( 0.121200.*exp( - 0.0105200.*(STATES(:,51)+CONSTANTS(:,133))))./(1.00000+exp( - 0.137800.*((STATES(:,51)+CONSTANTS(:,133))+40.1400)))); ALGEBRAIC(:,25) = ( 0.320000.*((STATES(:,51)+CONSTANTS(:,132))+58.4729))./(1.00000 - exp( - 0.100000.*((STATES(:,51)+CONSTANTS(:,132))+58.4729))); ALGEBRAIC(:,52) = 0.0800000.*exp(((STATES(:,51)+CONSTANTS(:,132)) - 13.7299)./ - 11.0000); ALGEBRAIC(:,27) = ( 0.320000.*(STATES(:,51)+47.1300))./(1.00000 - exp( - 0.100000.*(STATES(:,51)+47.1300))); ALGEBRAIC(:,53) = 0.0800000.*exp(STATES(:,51)./ - 11.0000); ALGEBRAIC(:,58) = (CONSTANTS(:,324) - STATES(:,15)) - STATES(:,14); ALGEBRAIC(:,61) = (CONSTANTS(:,332) - STATES(:,29)) - STATES(:,28); ALGEBRAIC(:,28) = 1.00000./(1.00000+exp((STATES(:,51)+9.43700)./ - 7.13300)); ALGEBRAIC(:,54) = 1.00000./(1.00000+exp((STATES(:,51) - 18.4099)./ - 29.3814)); ALGEBRAIC(:,66) = 1.00000./(1.00000+exp((STATES(:,51)+100.000)./29.3814)); ALGEBRAIC(:,78) = 1.00000./(ALGEBRAIC(:,54)./1.20890+ 3.50000.*ALGEBRAIC(:,66)); ALGEBRAIC(:,67) = ((1.00000./(1.00000+exp((STATES(:,51)+19.0000)./ - 9.00000)))./0.500000)./9.79530; ALGEBRAIC(:,29) = 1.00000./(1.00000+exp((STATES(:,51)+58.0000)./5.00000)); ALGEBRAIC(:,79) = 0.0479600.*ALGEBRAIC(:,29); ALGEBRAIC(:,80) = 0.0214400.*ALGEBRAIC(:,29); ALGEBRAIC(:,55) = (1.00000./(1.00000+exp((STATES(:,51)+60.0000)./5.00000)))./250.000; ALGEBRAIC(:,81) = 2.46000.*ALGEBRAIC(:,55); ALGEBRAIC(:,82) = 0.560340.*ALGEBRAIC(:,55); ALGEBRAIC(:,64) = 0.0906540.*exp( ( - 0.111570.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,92) = 0.388390./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.150190))./ - 0.606930)); ALGEBRAIC(:,76) = 0.00311240+(0.0283300 - 0.00311240)./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.0516600))./1.55220)); ALGEBRAIC(:,65) = 0.0657000.*exp( ( - 0.118990.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,94) = 0.563560./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.179860))./ - 0.583810)); ALGEBRAIC(:,77) = 0.000385250+(0.0124060 - 0.000385250)./(1.00000+exp(( CONSTANTS(:,264).*(STATES(:,51)+0.0641180))./0.779920)); ALGEBRAIC(:,93) = 0.000441980.*exp( ( - 1.20220.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,104) = 0.000401730.*exp( ( 0.000208730.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,95) = 0.000237300.*exp( ( - 1.97420.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,105) = 0.000226520.*exp( ( 0.000246900.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,6) = ( ( CONSTANTS(:,37).*CONSTANTS(:,36)).*CONSTANTS(:,320) - STATES(:,9)) - STATES(:,8); ALGEBRAIC(:,35) = ( ALGEBRAIC(:,6).*(CONSTANTS(:,60)+CONSTANTS(:,52)) - STATES(:,15).*(CONSTANTS(:,52)+CONSTANTS(:,60)))+ ( CONSTANTS(:,50).*CONSTANTS(:,52)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,55)); ALGEBRAIC(:,59) = ( - STATES(:,15).*CONSTANTS(:,50)).*CONSTANTS(:,52); ALGEBRAIC(:,71) = ( - ALGEBRAIC(:,35)+power(( ALGEBRAIC(:,35).*ALGEBRAIC(:,35) - ( 4.00000.*CONSTANTS(:,223)).*ALGEBRAIC(:,59)), 1.0 ./ 2))./( 2.00000.*CONSTANTS(:,223)); ALGEBRAIC(:,87) = ALGEBRAIC(:,6)./(1.00000+ (ALGEBRAIC(:,71)./CONSTANTS(:,50)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,52))); ALGEBRAIC(:,99) = ( ALGEBRAIC(:,71).*ALGEBRAIC(:,87))./CONSTANTS(:,50); ALGEBRAIC(:,109) = ( ALGEBRAIC(:,99).*CONSTANTS(:,60))./CONSTANTS(:,52); ALGEBRAIC(:,9) = ( ( CONSTANTS(:,216).*CONSTANTS(:,36)).*CONSTANTS(:,327) - STATES(:,22)) - STATES(:,21); ALGEBRAIC(:,38) = ( ALGEBRAIC(:,9).*(CONSTANTS(:,60)+CONSTANTS(:,52)) - STATES(:,29).*(CONSTANTS(:,52)+CONSTANTS(:,60)))+ ( CONSTANTS(:,50).*CONSTANTS(:,52)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,55)); ALGEBRAIC(:,62) = ( - STATES(:,29).*CONSTANTS(:,50)).*CONSTANTS(:,52); ALGEBRAIC(:,74) = ( - ALGEBRAIC(:,38)+power(( ALGEBRAIC(:,38).*ALGEBRAIC(:,38) - ( 4.00000.*CONSTANTS(:,226)).*ALGEBRAIC(:,62)), 1.0 ./ 2))./( 2.00000.*CONSTANTS(:,226)); ALGEBRAIC(:,90) = ALGEBRAIC(:,9)./(1.00000+ (ALGEBRAIC(:,74)./CONSTANTS(:,50)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,52))); ALGEBRAIC(:,102) = ( ALGEBRAIC(:,74).*ALGEBRAIC(:,90))./CONSTANTS(:,50); ALGEBRAIC(:,112) = (CONSTANTS(:,60)./CONSTANTS(:,52)).*ALGEBRAIC(:,102); ALGEBRAIC(:,8) = ( ( CONSTANTS(:,217).*CONSTANTS(:,296)).*CONSTANTS(:,322) - STATES(:,3)) - STATES(:,16); ALGEBRAIC(:,60) = ( ALGEBRAIC(:,8).*(CONSTANTS(:,48)+CONSTANTS(:,60)) - ALGEBRAIC(:,36).*(CONSTANTS(:,48)+CONSTANTS(:,60)))+ ( CONSTANTS(:,47).*CONSTANTS(:,48)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,49)); ALGEBRAIC(:,72) = ( - ALGEBRAIC(:,36).*CONSTANTS(:,48)).*CONSTANTS(:,47); ALGEBRAIC(:,88) = ( - ALGEBRAIC(:,60)+power(( ALGEBRAIC(:,60).*ALGEBRAIC(:,60) - ( 4.00000.*CONSTANTS(:,225)).*ALGEBRAIC(:,72)), 1.0 ./ 2))./( 2.00000.*CONSTANTS(:,225)); ALGEBRAIC(:,100) = ALGEBRAIC(:,8)./(1.00000+ (ALGEBRAIC(:,88)./CONSTANTS(:,47)).*(1.00000+CONSTANTS(:,60)./CONSTANTS(:,48))); ALGEBRAIC(:,111) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,88)).*ALGEBRAIC(:,100))./( CONSTANTS(:,47).*CONSTANTS(:,48)); ALGEBRAIC(:,120) = ( ALGEBRAIC(:,100).*ALGEBRAIC(:,88))./CONSTANTS(:,47); ALGEBRAIC(:,119) = ( CONSTANTS(:,60).*ALGEBRAIC(:,88))./CONSTANTS(:,49); ALGEBRAIC(:,131) = ((CONSTANTS(:,400) - STATES(:,143)) - STATES(:,142)) - STATES(:,141); ALGEBRAIC(:,2) = power( STATES(:,4).*1.00000, CONSTANTS(:,13)); ALGEBRAIC(:,32) = CONSTANTS(:,7).*(CONSTANTS(:,9)+ALGEBRAIC(:,2)./(CONSTANTS(:,4)+ALGEBRAIC(:,2))); ALGEBRAIC(:,68) = ( ALGEBRAIC(:,32).*CONSTANTS(:,330)).*CONSTANTS(:,228); ALGEBRAIC(:,84) = ( CONSTANTS(:,237).*(STATES(:,44) - STATES(:,46)))./CONSTANTS(:,310); ALGEBRAIC(:,97) = ( CONSTANTS(:,238).*(STATES(:,46) - STATES(:,45)))./CONSTANTS(:,310); ALGEBRAIC(:,108) = ( CONSTANTS(:,350).*CONSTANTS(:,186))./(1.00000+CONSTANTS(:,171)./STATES(:,46)); ALGEBRAIC(:,117) = ( (CONSTANTS(:,373)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,132)).*CONSTANTS(:,188))./(1.00000+CONSTANTS(:,173)./STATES(:,46)); ALGEBRAIC(:,125) = ALGEBRAIC(:,108)+ALGEBRAIC(:,117); ALGEBRAIC(:,135) = (( ( - CONSTANTS(:,269).*ALGEBRAIC(:,131)).*STATES(:,46)+ CONSTANTS(:,277).*STATES(:,143)) - ( CONSTANTS(:,270).*STATES(:,143)).*STATES(:,46))+ CONSTANTS(:,278).*STATES(:,142); ALGEBRAIC(:,141) = ((CONSTANTS(:,381) - STATES(:,135)) - STATES(:,134)) - STATES(:,133); ALGEBRAIC(:,3) = power( STATES(:,2).*1.00000, CONSTANTS(:,14)); ALGEBRAIC(:,33) = power( STATES(:,2).*1.00000, CONSTANTS(:,15)); ALGEBRAIC(:,57) = ( CONSTANTS(:,8).*(CONSTANTS(:,10)+ALGEBRAIC(:,3)./(CONSTANTS(:,5)+ALGEBRAIC(:,3)))).*(1.00000 - ( (1.00000 - ( CONSTANTS(:,16).*ALGEBRAIC(:,33))./(CONSTANTS(:,6)+ALGEBRAIC(:,33))).*STATES(:,1))./(CONSTANTS(:,3)+STATES(:,1))); ALGEBRAIC(:,83) = ( ALGEBRAIC(:,57).*CONSTANTS(:,321)).*CONSTANTS(:,228); ALGEBRAIC(:,96) = ( CONSTANTS(:,237).*(STATES(:,44) - STATES(:,46)))./CONSTANTS(:,307); ALGEBRAIC(:,106) = ( CONSTANTS(:,236).*(STATES(:,44) - STATES(:,45)))./CONSTANTS(:,307); ALGEBRAIC(:,115) = ( CONSTANTS(:,348).*CONSTANTS(:,186))./(1.00000+CONSTANTS(:,171)./STATES(:,44)); ALGEBRAIC(:,123) = ( (CONSTANTS(:,378)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,128)).*CONSTANTS(:,187))./(1.00000+CONSTANTS(:,172)./STATES(:,44)); ALGEBRAIC(:,129) = ( (CONSTANTS(:,371)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,130)).*CONSTANTS(:,188))./(1.00000+CONSTANTS(:,173)./STATES(:,44)); ALGEBRAIC(:,136) = (ALGEBRAIC(:,115)+ALGEBRAIC(:,123))+ALGEBRAIC(:,129); ALGEBRAIC(:,146) = (( ( - CONSTANTS(:,200).*ALGEBRAIC(:,141)).*STATES(:,44)+ CONSTANTS(:,247).*STATES(:,135)) - ( CONSTANTS(:,201).*STATES(:,135)).*STATES(:,44))+ CONSTANTS(:,248).*STATES(:,134); ALGEBRAIC(:,151) = ((CONSTANTS(:,336) - STATES(:,139)) - STATES(:,138)) - STATES(:,137); ALGEBRAIC(:,1) = power( STATES(:,3).*1.00000, CONSTANTS(:,13)); ALGEBRAIC(:,56) = CONSTANTS(:,7).*(CONSTANTS(:,9)+ALGEBRAIC(:,1)./(CONSTANTS(:,4)+ALGEBRAIC(:,1))); ALGEBRAIC(:,69) = ( ALGEBRAIC(:,56).*CONSTANTS(:,325)).*CONSTANTS(:,228); ALGEBRAIC(:,31) = power( STATES(:,3).*1.00000, CONSTANTS(:,14)); ALGEBRAIC(:,85) = CONSTANTS(:,8).*(CONSTANTS(:,10)+ALGEBRAIC(:,31)./(CONSTANTS(:,5)+ALGEBRAIC(:,31))); ALGEBRAIC(:,107) = ( ALGEBRAIC(:,85).*CONSTANTS(:,329)).*CONSTANTS(:,228); ALGEBRAIC(:,116) = ( CONSTANTS(:,236).*(STATES(:,44) - STATES(:,45)))./CONSTANTS(:,308); ALGEBRAIC(:,124) = ( CONSTANTS(:,238).*(STATES(:,46) - STATES(:,45)))./CONSTANTS(:,308); ALGEBRAIC(:,130) = ( CONSTANTS(:,349).*CONSTANTS(:,186))./(1.00000+CONSTANTS(:,171)./STATES(:,45)); ALGEBRAIC(:,137) = ( (CONSTANTS(:,379)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,129)).*CONSTANTS(:,187))./(1.00000+CONSTANTS(:,172)./STATES(:,45)); ALGEBRAIC(:,142) = ( (CONSTANTS(:,372)+ (CONSTANTS(:,175) - 1.00000).*STATES(:,131)).*CONSTANTS(:,188))./(1.00000+CONSTANTS(:,173)./STATES(:,45)); ALGEBRAIC(:,147) = (ALGEBRAIC(:,130)+ALGEBRAIC(:,137))+ALGEBRAIC(:,142); ALGEBRAIC(:,155) = (( ( - CONSTANTS(:,250).*ALGEBRAIC(:,151)).*STATES(:,45)+ CONSTANTS(:,266).*STATES(:,139)) - ( CONSTANTS(:,251).*STATES(:,139)).*STATES(:,45))+ CONSTANTS(:,267).*STATES(:,138); ALGEBRAIC(:,7) = ( ( CONSTANTS(:,38).*CONSTANTS(:,296)).*CONSTANTS(:,320) - STATES(:,2)) - STATES(:,10); ALGEBRAIC(:,70) = (((CONSTANTS(:,259)+CONSTANTS(:,272))./CONSTANTS(:,224)+ALGEBRAIC(:,34))+ALGEBRAIC(:,58)) - ALGEBRAIC(:,7); ALGEBRAIC(:,86) = (( CONSTANTS(:,272).*(ALGEBRAIC(:,34) - ALGEBRAIC(:,7))+ CONSTANTS(:,259).*(ALGEBRAIC(:,58) - ALGEBRAIC(:,7)))+CONSTANTS(:,280))./CONSTANTS(:,224); ALGEBRAIC(:,98) = ( ALGEBRAIC(:,7).*CONSTANTS(:,280))./CONSTANTS(:,224); ALGEBRAIC(:,110) = ((( ( - ALGEBRAIC(:,98)./27.0000).*power(ALGEBRAIC(:,70), 3.00000) - ( ( ( ALGEBRAIC(:,70).*ALGEBRAIC(:,70)).*ALGEBRAIC(:,86)).*ALGEBRAIC(:,86))./108.000)+( ( ALGEBRAIC(:,70).*ALGEBRAIC(:,86)).*ALGEBRAIC(:,98))./6.00000)+power(ALGEBRAIC(:,86), 3.00000)./27.0000)+( ALGEBRAIC(:,98).*ALGEBRAIC(:,98))./4.00000; ALGEBRAIC(:,118) = piecewise({ALGEBRAIC(:,110)<0.00000, power( - ALGEBRAIC(:,110), 1.0 ./ 2) }, 0.00000); ALGEBRAIC(:,126) = ((piecewise({ALGEBRAIC(:,110)>0.00000, power(ALGEBRAIC(:,110), 1.0 ./ 2) }, 0.00000)+ALGEBRAIC(:,98)./2.00000)+( ALGEBRAIC(:,70).*ALGEBRAIC(:,86))./6.00000) - power(ALGEBRAIC(:,70), 3.00000)./27.0000; ALGEBRAIC(:,132) = atan(ALGEBRAIC(:,118)./ALGEBRAIC(:,126))./3.00000; ALGEBRAIC(:,138) = power( ALGEBRAIC(:,126).*ALGEBRAIC(:,126)+ ALGEBRAIC(:,118).*ALGEBRAIC(:,118), 1.00000./6.00000); ALGEBRAIC(:,143) = (ALGEBRAIC(:,86)./3.00000 - ( ALGEBRAIC(:,70).*ALGEBRAIC(:,70))./9.00000)./( ALGEBRAIC(:,138).*ALGEBRAIC(:,138)); ALGEBRAIC(:,148) = ( ALGEBRAIC(:,138).* sin(ALGEBRAIC(:,132))).*(1.00000+ALGEBRAIC(:,143)); ALGEBRAIC(:,152) = ( ALGEBRAIC(:,138).*cos(ALGEBRAIC(:,132))).*(1.00000 - ALGEBRAIC(:,143)) - ALGEBRAIC(:,70)./3.00000; ALGEBRAIC(:,157) = power(( ALGEBRAIC(:,152).*ALGEBRAIC(:,152)+ ALGEBRAIC(:,148).*ALGEBRAIC(:,148)), 1.0 ./ 2); ALGEBRAIC(:,161) = ALGEBRAIC(:,34)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,49))+( ALGEBRAIC(:,157).*(CONSTANTS(:,48)+CONSTANTS(:,60)))./( CONSTANTS(:,47).*CONSTANTS(:,48))); ALGEBRAIC(:,165) = ( CONSTANTS(:,60).*ALGEBRAIC(:,161))./CONSTANTS(:,49); ALGEBRAIC(:,169) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,161)).*ALGEBRAIC(:,157))./( CONSTANTS(:,47).*CONSTANTS(:,48)); ALGEBRAIC(:,10) = ( ( CONSTANTS(:,39).*CONSTANTS(:,296)).*CONSTANTS(:,327) - STATES(:,4)) - STATES(:,24); ALGEBRAIC(:,73) = (((CONSTANTS(:,260)+CONSTANTS(:,273))./CONSTANTS(:,227)+ALGEBRAIC(:,37))+ALGEBRAIC(:,61)) - ALGEBRAIC(:,10); ALGEBRAIC(:,89) = (( CONSTANTS(:,273).*(ALGEBRAIC(:,37) - ALGEBRAIC(:,10))+ CONSTANTS(:,260).*(ALGEBRAIC(:,61) - ALGEBRAIC(:,10)))+CONSTANTS(:,281))./CONSTANTS(:,227); ALGEBRAIC(:,101) = ( ALGEBRAIC(:,10).*CONSTANTS(:,281))./CONSTANTS(:,227); ALGEBRAIC(:,113) = ((( ( - ALGEBRAIC(:,101)./27.0000).*power(ALGEBRAIC(:,73), 3.00000) - ( ( ( ALGEBRAIC(:,73).*ALGEBRAIC(:,73)).*ALGEBRAIC(:,89)).*ALGEBRAIC(:,89))./108.000)+( ( ALGEBRAIC(:,73).*ALGEBRAIC(:,89)).*ALGEBRAIC(:,101))./6.00000)+power(ALGEBRAIC(:,89), 3.00000)./27.0000)+( ALGEBRAIC(:,101).*ALGEBRAIC(:,101))./4.00000; ALGEBRAIC(:,121) = piecewise({ALGEBRAIC(:,113)<0.00000, power( - ALGEBRAIC(:,113), 1.0 ./ 2) }, 0.00000); ALGEBRAIC(:,127) = ((piecewise({ALGEBRAIC(:,113)>0.00000, power(ALGEBRAIC(:,113), 1.0 ./ 2) }, 0.00000)+ALGEBRAIC(:,101)./2.00000)+( ALGEBRAIC(:,73).*ALGEBRAIC(:,89))./6.00000) - power(ALGEBRAIC(:,73), 3.00000)./27.0000; ALGEBRAIC(:,133) = atan(ALGEBRAIC(:,121)./ALGEBRAIC(:,127))./3.00000; ALGEBRAIC(:,139) = power( ALGEBRAIC(:,127).*ALGEBRAIC(:,127)+ ALGEBRAIC(:,121).*ALGEBRAIC(:,121), 1.00000./6.00000); ALGEBRAIC(:,144) = (ALGEBRAIC(:,89)./3.00000 - ( ALGEBRAIC(:,73).*ALGEBRAIC(:,73))./9.00000)./( ALGEBRAIC(:,139).*ALGEBRAIC(:,139)); ALGEBRAIC(:,149) = ( ALGEBRAIC(:,139).* sin(ALGEBRAIC(:,133))).*(1.00000+ALGEBRAIC(:,144)); ALGEBRAIC(:,153) = ( ALGEBRAIC(:,139).*cos(ALGEBRAIC(:,133))).*(1.00000 - ALGEBRAIC(:,144)) - ALGEBRAIC(:,73)./3.00000; ALGEBRAIC(:,158) = power(( ALGEBRAIC(:,153).*ALGEBRAIC(:,153)+ ALGEBRAIC(:,149).*ALGEBRAIC(:,149)), 1.0 ./ 2); ALGEBRAIC(:,162) = ALGEBRAIC(:,37)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,49))+( ALGEBRAIC(:,158).*(CONSTANTS(:,48)+CONSTANTS(:,60)))./( CONSTANTS(:,47).*CONSTANTS(:,48))); ALGEBRAIC(:,166) = ( CONSTANTS(:,60).*ALGEBRAIC(:,162))./CONSTANTS(:,49); ALGEBRAIC(:,170) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,162)).*ALGEBRAIC(:,158))./( CONSTANTS(:,47).*CONSTANTS(:,48)); ALGEBRAIC(:,177) = ALGEBRAIC(:,58)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,54))+( ALGEBRAIC(:,157).*(CONSTANTS(:,53)+CONSTANTS(:,60)))./( CONSTANTS(:,51).*CONSTANTS(:,53))); ALGEBRAIC(:,181) = ( CONSTANTS(:,60).*ALGEBRAIC(:,177))./CONSTANTS(:,54); ALGEBRAIC(:,185) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,177)).*ALGEBRAIC(:,157))./( CONSTANTS(:,51).*CONSTANTS(:,53)); ALGEBRAIC(:,178) = ALGEBRAIC(:,61)./((1.00000+CONSTANTS(:,60)./CONSTANTS(:,54))+( ALGEBRAIC(:,158).*(CONSTANTS(:,53)+CONSTANTS(:,60)))./( CONSTANTS(:,51).*CONSTANTS(:,53))); ALGEBRAIC(:,182) = ( CONSTANTS(:,60).*ALGEBRAIC(:,178))./CONSTANTS(:,54); ALGEBRAIC(:,186) = ( ( CONSTANTS(:,60).*ALGEBRAIC(:,178)).*ALGEBRAIC(:,158))./( CONSTANTS(:,51).*CONSTANTS(:,53)); ALGEBRAIC(:,189) = ALGEBRAIC(:,169)+ CONSTANTS(:,61).*ALGEBRAIC(:,185); ALGEBRAIC(:,173) = ( ALGEBRAIC(:,161).*ALGEBRAIC(:,157))./CONSTANTS(:,47); ALGEBRAIC(:,193) = ( ALGEBRAIC(:,177).*ALGEBRAIC(:,157))./CONSTANTS(:,51); ALGEBRAIC(:,197) = ALGEBRAIC(:,173)+ CONSTANTS(:,61).*ALGEBRAIC(:,193); ALGEBRAIC(:,190) = ALGEBRAIC(:,170)+ CONSTANTS(:,64).*ALGEBRAIC(:,186); ALGEBRAIC(:,174) = ( ALGEBRAIC(:,162).*ALGEBRAIC(:,158))./CONSTANTS(:,47); ALGEBRAIC(:,194) = ( ALGEBRAIC(:,178).*ALGEBRAIC(:,158))./CONSTANTS(:,51); ALGEBRAIC(:,198) = ALGEBRAIC(:,174)+ CONSTANTS(:,64).*ALGEBRAIC(:,194); ALGEBRAIC(:,134) = CONSTANTS(:,275) - STATES(:,36); ALGEBRAIC(:,140) = CONSTANTS(:,263) - STATES(:,36).*CONSTANTS(:,262); ALGEBRAIC(:,145) = - CONSTANTS(:,231).*STATES(:,36); ALGEBRAIC(:,150) = - ALGEBRAIC(:,134)./3.00000+ ( (2.00000./3.00000).*power(( ALGEBRAIC(:,134).*ALGEBRAIC(:,134) - 3.00000.*ALGEBRAIC(:,140)), 1.0 ./ 2)).*cos(acos((( ( 9.00000.*ALGEBRAIC(:,134)).*ALGEBRAIC(:,140) - ( ( 2.00000.*ALGEBRAIC(:,134)).*ALGEBRAIC(:,134)).*ALGEBRAIC(:,134)) - 27.0000.*ALGEBRAIC(:,145))./( 2.00000.*power( ALGEBRAIC(:,134).*ALGEBRAIC(:,134) - 3.00000.*ALGEBRAIC(:,140), 1.50000)))./3.00000); ALGEBRAIC(:,191) = ( CONSTANTS(:,80).*(1.00000 - STATES(:,43)))./(1.00000+CONSTANTS(:,82)./ALGEBRAIC(:,150)); ALGEBRAIC(:,195) = ALGEBRAIC(:,191)+STATES(:,43); ALGEBRAIC(:,200) = ALGEBRAIC(:,195)./(ALGEBRAIC(:,195)+CONSTANTS(:,81)); ALGEBRAIC(:,201) = 1.00000./(1.00000+power(CONSTANTS(:,81)./ALGEBRAIC(:,195), 2.00000)); ALGEBRAIC(:,199) = (CONSTANTS(:,206) - CONSTANTS(:,209))+STATES(:,145); ALGEBRAIC(:,204) = 0.500000.*(power((power(ALGEBRAIC(:,199), 2.00000)+ ( 4.00000.*CONSTANTS(:,206)).*CONSTANTS(:,209)), 1.0 ./ 2) - ALGEBRAIC(:,199)); ALGEBRAIC(:,207) = (CONSTANTS(:,29)./CONSTANTS(:,320)+CONSTANTS(:,83)./CONSTANTS(:,327))+ALGEBRAIC(:,204)./CONSTANTS(:,322); ALGEBRAIC(:,156) = 1.00000./( (1.00000+exp((4.79800+STATES(:,51))./ - 7.56990)).*(1.00000+exp((25.0000+STATES(:,51))./ - 5.00000))); ALGEBRAIC(:,196) = 0.590000+( 0.800000.*exp( 0.0520000.*(STATES(:,51)+13.0000)))./(1.00000+exp( 0.132000.*(STATES(:,51)+13.0000))); ALGEBRAIC(:,203) = ALGEBRAIC(:,156)./ALGEBRAIC(:,196); ALGEBRAIC(:,206) = (1.00000 - ALGEBRAIC(:,156))./ALGEBRAIC(:,196); ALGEBRAIC(:,180) = 1.00000+power(0.00200000./ALGEBRAIC(:,150), 4.00000); ALGEBRAIC(:,184) = 6.00000./ALGEBRAIC(:,180); ALGEBRAIC(:,164) = 1.00000./(1.00000+exp((29.9790+STATES(:,51))./3.17750)); ALGEBRAIC(:,172) = (0.100000+ALGEBRAIC(:,164))./1.10000; ALGEBRAIC(:,215) = 1.00000./( ( 70.0000.*(1.00000 - 0.500000.*STATES(:,37))).*(1.00000+exp((STATES(:,51)+49.1000)./10.3490))); ALGEBRAIC(:,217) = 1.00000./(1.00000+exp((STATES(:,51)+0.213000)./ - 10.8070)); ALGEBRAIC(:,222) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./38.4940); ALGEBRAIC(:,225) = ALGEBRAIC(:,172)./ALGEBRAIC(:,222); ALGEBRAIC(:,234) = (1.00000 - ALGEBRAIC(:,172))./ALGEBRAIC(:,222); ALGEBRAIC(:,168) = (0.000100000+ALGEBRAIC(:,164))./1.00010; ALGEBRAIC(:,160) = 0.100000.*STATES(:,37); ALGEBRAIC(:,176) = 1.00000+power(0.0100000./ALGEBRAIC(:,150), 10.0000); ALGEBRAIC(:,188) = (32.5000 - (18.0000 - ALGEBRAIC(:,160))./ALGEBRAIC(:,180)) - 10.0000./ALGEBRAIC(:,176); ALGEBRAIC(:,219) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./ALGEBRAIC(:,188)); ALGEBRAIC(:,228) = ALGEBRAIC(:,168)./ALGEBRAIC(:,219); ALGEBRAIC(:,240) = (1.00000 - ALGEBRAIC(:,168))./ALGEBRAIC(:,219); ALGEBRAIC(:,231) = piecewise({abs(ALGEBRAIC(:,228))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,228)); ALGEBRAIC(:,237) = piecewise({abs(ALGEBRAIC(:,234))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,234)); ALGEBRAIC(:,243) = ( CONSTANTS(:,104).*( ( ALGEBRAIC(:,225).*ALGEBRAIC(:,240)).*ALGEBRAIC(:,184)))./( ( ALGEBRAIC(:,237).*ALGEBRAIC(:,231)).*CONSTANTS(:,103)); ALGEBRAIC(:,63) = CONSTANTS(:,275) - STATES(:,34); ALGEBRAIC(:,75) = CONSTANTS(:,263) - STATES(:,34).*CONSTANTS(:,262); ALGEBRAIC(:,91) = - CONSTANTS(:,231).*STATES(:,34); ALGEBRAIC(:,103) = - ALGEBRAIC(:,63)./3.00000+ ( (2.00000./3.00000).*power(( ALGEBRAIC(:,63).*ALGEBRAIC(:,63) - 3.00000.*ALGEBRAIC(:,75)), 1.0 ./ 2)).*cos(acos((( ( 9.00000.*ALGEBRAIC(:,63)).*ALGEBRAIC(:,75) - ( ( 2.00000.*ALGEBRAIC(:,63)).*ALGEBRAIC(:,63)).*ALGEBRAIC(:,63)) - 27.0000.*ALGEBRAIC(:,91))./( 2.00000.*power( ALGEBRAIC(:,63).*ALGEBRAIC(:,63) - 3.00000.*ALGEBRAIC(:,75), 1.50000)))./3.00000); ALGEBRAIC(:,210) = (ALGEBRAIC(:,150) - ALGEBRAIC(:,103))./CONSTANTS(:,96); ALGEBRAIC(:,11) = (STATES(:,5)+CONSTANTS(:,393))./CONSTANTS(:,18); ALGEBRAIC(:,221) = (ALGEBRAIC(:,11) - CONSTANTS(:,398))./(0.927300 - CONSTANTS(:,398)); ALGEBRAIC(:,224) = piecewise({ALGEBRAIC(:,221)<0.00000, 0.00000 }, ALGEBRAIC(:,221)); ALGEBRAIC(:,227) = 0.000257900.*(1.00000+ 0.100000.*STATES(:,37)); ALGEBRAIC(:,230) = exp( ( 2.00000.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,233) = ( ( ( ( ALGEBRAIC(:,227).*4.00000).*STATES(:,51)).*CONSTANTS(:,276)).*( ALGEBRAIC(:,150).*ALGEBRAIC(:,230) - 0.341000.*CONSTANTS(:,98)))./(ALGEBRAIC(:,230) - 1.00000); ALGEBRAIC(:,236) = ALGEBRAIC(:,233).*(STATES(:,56)+STATES(:,59)); ALGEBRAIC(:,239) = 0.000155200.*(1.00000+ 0.400000.*STATES(:,37)); ALGEBRAIC(:,242) = exp( ( 2.00000.*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,245) = ( ( ( ( ALGEBRAIC(:,239).*4.00000).*STATES(:,51)).*CONSTANTS(:,276)).*( ALGEBRAIC(:,150).*ALGEBRAIC(:,242) - 0.341000.*CONSTANTS(:,98)))./(ALGEBRAIC(:,242) - 1.00000); ALGEBRAIC(:,247) = ALGEBRAIC(:,245).*(STATES(:,64)+STATES(:,67)); ALGEBRAIC(:,249) = (1.00000 - ALGEBRAIC(:,224)).*ALGEBRAIC(:,247)+ ALGEBRAIC(:,224).*ALGEBRAIC(:,236); ALGEBRAIC(:,192) = 1.00000./( (1.00000+exp((13.5600 - STATES(:,51))./9.45000)).*(1.00000+exp((25.0000+STATES(:,51))./ - 5.00000))); ALGEBRAIC(:,209) = ALGEBRAIC(:,192)./ALGEBRAIC(:,196); ALGEBRAIC(:,211) = (1.00000 - ALGEBRAIC(:,192))./ALGEBRAIC(:,196); ALGEBRAIC(:,235) = 1.00000+power(0.00110000./ALGEBRAIC(:,150), 4.00000); ALGEBRAIC(:,238) = 14.9186./ALGEBRAIC(:,235); ALGEBRAIC(:,220) = 1.00000./(1.00000+exp((17.5000+STATES(:,51))./3.00000)); ALGEBRAIC(:,226) = (0.247400+ALGEBRAIC(:,220))./1.24740; ALGEBRAIC(:,229) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./26.5530); ALGEBRAIC(:,246) = ALGEBRAIC(:,226)./ALGEBRAIC(:,229); ALGEBRAIC(:,252) = (1.00000 - ALGEBRAIC(:,226))./ALGEBRAIC(:,229); ALGEBRAIC(:,223) = (0.00100000+ALGEBRAIC(:,220))./1.00100; ALGEBRAIC(:,213) = 5.00000.*STATES(:,37); ALGEBRAIC(:,232) = 1.00000+power(0.0120000./ALGEBRAIC(:,150), 10.0000); ALGEBRAIC(:,241) = (13.8250 - (6.38360 - ALGEBRAIC(:,213))./ALGEBRAIC(:,235)) - 3.36960./ALGEBRAIC(:,232); ALGEBRAIC(:,244) = 1.00000./(ALGEBRAIC(:,215)+ALGEBRAIC(:,217)./ALGEBRAIC(:,241)); ALGEBRAIC(:,248) = ALGEBRAIC(:,223)./ALGEBRAIC(:,244); ALGEBRAIC(:,258) = (1.00000 - ALGEBRAIC(:,223))./ALGEBRAIC(:,244); ALGEBRAIC(:,250) = piecewise({abs(ALGEBRAIC(:,248))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,248)); ALGEBRAIC(:,255) = piecewise({abs(ALGEBRAIC(:,252))<1.00000e-12, 1.00000e-12 }, ALGEBRAIC(:,252)); ALGEBRAIC(:,261) = ( CONSTANTS(:,106).*( ( ALGEBRAIC(:,246).*ALGEBRAIC(:,258)).*ALGEBRAIC(:,238)))./( ( ALGEBRAIC(:,255).*ALGEBRAIC(:,250)).*CONSTANTS(:,105)); ALGEBRAIC(:,253) = CONSTANTS(:,244).*(1.00000+ 2.00000.*STATES(:,42)); ALGEBRAIC(:,256) = 0.112500.*ALGEBRAIC(:,253); ALGEBRAIC(:,114) = (CONSTANTS(:,73)+CONSTANTS(:,74)) - STATES(:,35); ALGEBRAIC(:,122) = STATES(:,35).*CONSTANTS(:,74); ALGEBRAIC(:,128) = (power(( ALGEBRAIC(:,114).*ALGEBRAIC(:,114)+ 4.00000.*ALGEBRAIC(:,122)), 1.0 ./ 2) - ALGEBRAIC(:,114))./2.00000; ALGEBRAIC(:,270) = ( ALGEBRAIC(:,249).*1.00000)./(1.00000+power(1.00000./ALGEBRAIC(:,128), 8.00000)); ALGEBRAIC(:,272) = ALGEBRAIC(:,256).*ALGEBRAIC(:,270); ALGEBRAIC(:,264) = 1.00000+0.0123000./ALGEBRAIC(:,128); ALGEBRAIC(:,266) = ALGEBRAIC(:,253)./ALGEBRAIC(:,264); ALGEBRAIC(:,259) = CONSTANTS(:,244).*(1.00000+ 0.00000.*STATES(:,42)); ALGEBRAIC(:,262) = 0.112500.*ALGEBRAIC(:,259); ALGEBRAIC(:,273) = ( 1.99250.*ALGEBRAIC(:,262)).*ALGEBRAIC(:,270); ALGEBRAIC(:,268) = ( 0.535700.*ALGEBRAIC(:,259))./ALGEBRAIC(:,264); ALGEBRAIC(:,287) = 1.00000+power(CONSTANTS(:,135)./ALGEBRAIC(:,150), 2.00000); ALGEBRAIC(:,286) = exp( ( (CONSTANTS(:,140) - 1.00000).*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,288) = 1.00000+ CONSTANTS(:,141).*ALGEBRAIC(:,286); ALGEBRAIC(:,284) = power(STATES(:,50), 3.00000); ALGEBRAIC(:,289) = ( CONSTANTS(:,137).*ALGEBRAIC(:,284)+ CONSTANTS(:,365).*ALGEBRAIC(:,150))+ ( CONSTANTS(:,363).*CONSTANTS(:,98)).*(1.00000+ALGEBRAIC(:,150)./CONSTANTS(:,136)); ALGEBRAIC(:,290) = ( ( CONSTANTS(:,136).*CONSTANTS(:,367)).*(1.00000+ALGEBRAIC(:,284)./CONSTANTS(:,363))+ ALGEBRAIC(:,284).*CONSTANTS(:,98))+ CONSTANTS(:,367).*ALGEBRAIC(:,150); ALGEBRAIC(:,285) = exp( ( CONSTANTS(:,140).*STATES(:,51)).*CONSTANTS(:,264)); ALGEBRAIC(:,291) = ( 0.200000.*CONSTANTS(:,142)).*( ( ALGEBRAIC(:,284).*CONSTANTS(:,98)).*ALGEBRAIC(:,285) - ( CONSTANTS(:,367).*ALGEBRAIC(:,150)).*ALGEBRAIC(:,286)); ALGEBRAIC(:,292) = ALGEBRAIC(:,291)./( ( ALGEBRAIC(:,287).*ALGEBRAIC(:,288)).*(ALGEBRAIC(:,289)+ALGEBRAIC(:,290))); ALGEBRAIC(:,208) = (STATES(:,50) - STATES(:,49))./CONSTANTS(:,95); ALGEBRAIC(:,39) = (STATES(:,6)+CONSTANTS(:,396))./CONSTANTS(:,31); ALGEBRAIC(:,312) = (ALGEBRAIC(:,39) - CONSTANTS(:,399))./(0.958600 - CONSTANTS(:,399)); ALGEBRAIC(:,314) = piecewise({ALGEBRAIC(:,312)<0.00000, 0.00000 }, ALGEBRAIC(:,312)); ALGEBRAIC(:,318) = (1.00000 - ALGEBRAIC(:,314)).*STATES(:,119)+ ALGEBRAIC(:,314).*STATES(:,120); ALGEBRAIC(:,308) = ( CONSTANTS(:,156).*exp(ALGEBRAIC(:,128)./CONSTANTS(:,154))).*(ALGEBRAIC(:,128) - ALGEBRAIC(:,150)); ALGEBRAIC(:,310) = ( CONSTANTS(:,157).*exp(ALGEBRAIC(:,128)./CONSTANTS(:,155))).*(ALGEBRAIC(:,128) - ALGEBRAIC(:,150)); ALGEBRAIC(:,316) = (1.00000 - ALGEBRAIC(:,314)).*ALGEBRAIC(:,308)+ ALGEBRAIC(:,314).*ALGEBRAIC(:,310); ALGEBRAIC(:,320) = ALGEBRAIC(:,316)+ALGEBRAIC(:,318); ALGEBRAIC(:,212) = (STATES(:,31) - ALGEBRAIC(:,128))./CONSTANTS(:,97); ALGEBRAIC(:,154) = (STATES(:,32) - 0.673519)./(0.999180 - 0.673519); ALGEBRAIC(:,159) = piecewise({ALGEBRAIC(:,154)<0.00000, 0.00000 }, ALGEBRAIC(:,154)); ALGEBRAIC(:,163) = (1.00000 - ALGEBRAIC(:,159)).*CONSTANTS(:,78)+ ALGEBRAIC(:,159).*CONSTANTS(:,230); ALGEBRAIC(:,175) = ALGEBRAIC(:,163)+CONSTANTS(:,76); ALGEBRAIC(:,179) = ((ALGEBRAIC(:,175) - STATES(:,33))+CONSTANTS(:,72))+CONSTANTS(:,79); ALGEBRAIC(:,167) = ALGEBRAIC(:,163).*CONSTANTS(:,76); ALGEBRAIC(:,183) = ((ALGEBRAIC(:,167) - STATES(:,33).*ALGEBRAIC(:,175))+ CONSTANTS(:,79).*CONSTANTS(:,76))+ CONSTANTS(:,72).*ALGEBRAIC(:,163); ALGEBRAIC(:,171) = - ALGEBRAIC(:,167).*STATES(:,33); ALGEBRAIC(:,187) = - ALGEBRAIC(:,179)./3.00000+ ( (2.00000./3.00000).*power(( ALGEBRAIC(:,179).*ALGEBRAIC(:,179) - 3.00000.*ALGEBRAIC(:,183)), 1.0 ./ 2)).*cos(acos((( ( 9.00000.*ALGEBRAIC(:,179)).*ALGEBRAIC(:,183) - ( ( 2.00000.*ALGEBRAIC(:,179)).*ALGEBRAIC(:,179)).*ALGEBRAIC(:,179)) - 27.0000.*ALGEBRAIC(:,171))./( 2.00000.*power( ALGEBRAIC(:,179).*ALGEBRAIC(:,179) - 3.00000.*ALGEBRAIC(:,183), 1.50000)))./3.00000); ALGEBRAIC(:,202) = (ALGEBRAIC(:,150) - ALGEBRAIC(:,187))./CONSTANTS(:,95); ALGEBRAIC(:,251) = exp( STATES(:,51).*CONSTANTS(:,264)); ALGEBRAIC(:,254) = ( ( ( CONSTANTS(:,108).*STATES(:,51)).*CONSTANTS(:,276)).*(STATES(:,47) - CONSTANTS(:,99).*ALGEBRAIC(:,251)))./(1.00000 - ALGEBRAIC(:,251)); ALGEBRAIC(:,321) = 1.00000 - 1.00000./(1.00000+power(ALGEBRAIC(:,318)./CONSTANTS(:,109), 2.00000)); ALGEBRAIC(:,324) = ( ALGEBRAIC(:,254).*ALGEBRAIC(:,321)).*STATES(:,68); ALGEBRAIC(:,205) = (STATES(:,48) - STATES(:,47))./CONSTANTS(:,95); ALGEBRAIC(:,307) = 1.00000./(1.00000+power(0.0300000./ALGEBRAIC(:,195), 2.00000)); ALGEBRAIC(:,309) = (1.00000 - ALGEBRAIC(:,307)).*CONSTANTS(:,162)+ ALGEBRAIC(:,307).*CONSTANTS(:,380); ALGEBRAIC(:,325) = ( ALGEBRAIC(:,309).*STATES(:,31))./CONSTANTS(:,165); ALGEBRAIC(:,311) = (STATES(:,126) - 0.659100)./(0.994500 - 0.659100); ALGEBRAIC(:,313) = piecewise({ALGEBRAIC(:,311)<0.00000, 0.00000 }, ALGEBRAIC(:,311)); ALGEBRAIC(:,315) = ALGEBRAIC(:,313).*STATES(:,41); ALGEBRAIC(:,317) = STATES(:,41) - ALGEBRAIC(:,315); ALGEBRAIC(:,319) = ALGEBRAIC(:,313) - ALGEBRAIC(:,315); ALGEBRAIC(:,322) = ((1.00000 - ALGEBRAIC(:,319)) - ALGEBRAIC(:,317)) - ALGEBRAIC(:,315); ALGEBRAIC(:,323) = (( ALGEBRAIC(:,322).*CONSTANTS(:,160)+ ALGEBRAIC(:,319).*CONSTANTS(:,375))+ ALGEBRAIC(:,317).*CONSTANTS(:,370))+ ALGEBRAIC(:,315).*CONSTANTS(:,377); ALGEBRAIC(:,327) = ( ALGEBRAIC(:,309).*ALGEBRAIC(:,187))./(ALGEBRAIC(:,187)+ALGEBRAIC(:,323)); ALGEBRAIC(:,329) = ALGEBRAIC(:,327) - ALGEBRAIC(:,325); ALGEBRAIC(:,214) = ( 2.00000.*STATES(:,51)).*CONSTANTS(:,264); ALGEBRAIC(:,216) = exp(ALGEBRAIC(:,214)); ALGEBRAIC(:,218) = ( ( ( ( CONSTANTS(:,102).*2.00000).*CONSTANTS(:,88)).*ALGEBRAIC(:,214)).*( ALGEBRAIC(:,187).*ALGEBRAIC(:,216) - 0.341000.*CONSTANTS(:,98)))./(ALGEBRAIC(:,216) - 1.00000); ALGEBRAIC(:,293) = 1.00000+power(CONSTANTS(:,135)./ALGEBRAIC(:,187), 2.00000); ALGEBRAIC(:,294) = 1.00000+ CONSTANTS(:,141).*ALGEBRAIC(:,286); ALGEBRAIC(:,283) = power(STATES(:,49), 3.00000); ALGEBRAIC(:,295) = ( CONSTANTS(:,137).*ALGEBRAIC(:,283)+ CONSTANTS(:,365).*ALGEBRAIC(:,187))+ ( CONSTANTS(:,363).*CONSTANTS(:,98)).*(1.00000+ALGEBRAIC(:,187)./CONSTANTS(:,136)); ALGEBRAIC(:,296) = ( ( CONSTANTS(:,136).*CONSTANTS(:,367)).*(1.00000+ALGEBRAIC(:,283)./CONSTANTS(:,363))+ ALGEBRAIC(:,283).*CONSTANTS(:,98))+ CONSTANTS(:,367).*ALGEBRAIC(:,187); ALGEBRAIC(:,297) = ( 0.800000.*CONSTANTS(:,142)).*( ( ALGEBRAIC(:,283).*CONSTANTS(:,98)).*ALGEBRAIC(:,285) - ( CONSTANTS(:,367).*ALGEBRAIC(:,187)).*ALGEBRAIC(:,286)); ALGEBRAIC(:,298) = ALGEBRAIC(:,297)./( ( ALGEBRAIC(:,293).*ALGEBRAIC(:,294)).*(ALGEBRAIC(:,295)+ALGEBRAIC(:,296))); ALGEBRAIC(:,305) = ( CONSTANTS(:,152).*ALGEBRAIC(:,187))./(CONSTANTS(:,153)+ALGEBRAIC(:,187)); ALGEBRAIC(:,328) = - CONSTANTS(:,239).*log(CONSTANTS(:,99)./STATES(:,47)); ALGEBRAIC(:,332) = CONSTANTS(:,239).*log(CONSTANTS(:,100)./STATES(:,127)); ALGEBRAIC(:,333) = ALGEBRAIC(:,332) - ALGEBRAIC(:,328); ALGEBRAIC(:,334) = ( CONSTANTS(:,92).*ALGEBRAIC(:,333))./(ALGEBRAIC(:,333)+CONSTANTS(:,93)); ALGEBRAIC(:,353) = CONSTANTS(:,239).*log(CONSTANTS(:,101)./STATES(:,49)); ALGEBRAIC(:,354) = power(ALGEBRAIC(:,353) - ALGEBRAIC(:,328), 4.00000); ALGEBRAIC(:,355) = ( CONSTANTS(:,94).*ALGEBRAIC(:,354))./(ALGEBRAIC(:,354)+CONSTANTS(:,347)); ALGEBRAIC(:,330) = CONSTANTS(:,107).*(STATES(:,51) - ALGEBRAIC(:,328)); ALGEBRAIC(:,356) = ( ( ( CONSTANTS(:,242).*power(STATES(:,115), 3.00000)).*STATES(:,107)).*STATES(:,108)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,357) = ( ( ( CONSTANTS(:,242).*power(STATES(:,109), 3.00000)).*STATES(:,107)).*STATES(:,108)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,358) = ( ( ( CONSTANTS(:,242).*power(STATES(:,112), 3.00000)).*STATES(:,110)).*STATES(:,111)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,359) = ( ( ( ( CONSTANTS(:,242).*power(STATES(:,115), 3.00000)).*STATES(:,113)).*STATES(:,114)).*(STATES(:,51) - ALGEBRAIC(:,353))).*1.25000; ALGEBRAIC(:,274) = (STATES(:,106) - 0.239480)./(0.950143 - 0.239480); ALGEBRAIC(:,275) = piecewise({ALGEBRAIC(:,274)<0.00000, 0.00000 }, ALGEBRAIC(:,274)); ALGEBRAIC(:,276) = ALGEBRAIC(:,275).*STATES(:,39); ALGEBRAIC(:,277) = STATES(:,39) - ALGEBRAIC(:,276); ALGEBRAIC(:,278) = ALGEBRAIC(:,275) - ALGEBRAIC(:,276); ALGEBRAIC(:,279) = ((1.00000 - ALGEBRAIC(:,278)) - ALGEBRAIC(:,277)) - ALGEBRAIC(:,276); ALGEBRAIC(:,360) = (( ALGEBRAIC(:,279).*ALGEBRAIC(:,358)+ ALGEBRAIC(:,278).*ALGEBRAIC(:,359))+ ALGEBRAIC(:,277).*ALGEBRAIC(:,357))+ ALGEBRAIC(:,276).*ALGEBRAIC(:,356); ALGEBRAIC(:,280) = STATES(:,51).*CONSTANTS(:,264); ALGEBRAIC(:,281) = exp(ALGEBRAIC(:,280)); ALGEBRAIC(:,282) = ( ( ( CONSTANTS(:,134).*CONSTANTS(:,88)).*ALGEBRAIC(:,280)).*( STATES(:,49).*ALGEBRAIC(:,281) - CONSTANTS(:,101)))./(ALGEBRAIC(:,281) - 1.00000); ALGEBRAIC(:,301) = ( CONSTANTS(:,145).*CONSTANTS(:,243))./(1.00000+exp( - (STATES(:,51)+92.0000).*CONSTANTS(:,264))); ALGEBRAIC(:,302) = ALGEBRAIC(:,301).*power(STATES(:,49)./(STATES(:,49)+CONSTANTS(:,148)), 3.00000); ALGEBRAIC(:,303) = ALGEBRAIC(:,301).*power(STATES(:,49)./(STATES(:,49)+CONSTANTS(:,149)), 3.00000); ALGEBRAIC(:,299) = (STATES(:,116) - 0.126345)./(0.998014 - 0.126345); ALGEBRAIC(:,300) = piecewise({ALGEBRAIC(:,299)<0.00000, 0.00000 }, ALGEBRAIC(:,299)); ALGEBRAIC(:,304) = (1.00000 - ALGEBRAIC(:,300)).*ALGEBRAIC(:,302)+ ALGEBRAIC(:,300).*ALGEBRAIC(:,303); ALGEBRAIC(:,361) = ( power(STATES(:,118), 3.00000).*STATES(:,117)).*(STATES(:,51) - ALGEBRAIC(:,353)); ALGEBRAIC(:,362) = 0.0160000.*ALGEBRAIC(:,361); ALGEBRAIC(:,363) = 0.00650000.*ALGEBRAIC(:,361); ALGEBRAIC(:,364) = (1.00000 - STATES(:,39)).*ALGEBRAIC(:,363)+ STATES(:,39).*ALGEBRAIC(:,362); ALGEBRAIC(:,366) = (((ALGEBRAIC(:,360)+ALGEBRAIC(:,282))+ALGEBRAIC(:,364))+ ALGEBRAIC(:,304).*3.00000)+ ALGEBRAIC(:,298).*3.00000; ALGEBRAIC(:,306) = (((ALGEBRAIC(:,249)+ALGEBRAIC(:,218))+ALGEBRAIC(:,305)) - ALGEBRAIC(:,298).*2.00000) - ALGEBRAIC(:,292).*2.00000; ALGEBRAIC(:,331) = ALGEBRAIC(:,330)+ALGEBRAIC(:,324); ALGEBRAIC(:,335) = STATES(:,51) - ALGEBRAIC(:,332); ALGEBRAIC(:,336) = 1.02000./(1.00000+exp( 0.238500.*(ALGEBRAIC(:,335) - 59.2150))); ALGEBRAIC(:,337) = ( 0.491240.*exp( 0.0803200.*(ALGEBRAIC(:,335)+5.47600))+exp( 0.0617500.*(ALGEBRAIC(:,335) - 594.310)))./(1.00000+exp( - 0.514300.*(ALGEBRAIC(:,335)+4.75300))); ALGEBRAIC(:,338) = ( CONSTANTS(:,240).*(ALGEBRAIC(:,336)./(ALGEBRAIC(:,336)+ALGEBRAIC(:,337)))).*ALGEBRAIC(:,335); ALGEBRAIC(:,339) = ALGEBRAIC(:,338).*1.20000; ALGEBRAIC(:,340) = (1.00000 - STATES(:,38)).*ALGEBRAIC(:,338)+ STATES(:,38).*ALGEBRAIC(:,339); ALGEBRAIC(:,257) = 1.00000./(1.00000+exp((STATES(:,51)+10.0000)./15.4000)); ALGEBRAIC(:,341) = ( ( CONSTANTS(:,241).*STATES(:,69)).*ALGEBRAIC(:,257)).*(STATES(:,51) - ALGEBRAIC(:,332)); ALGEBRAIC(:,349) = CONSTANTS(:,239).*log((CONSTANTS(:,100)+ CONSTANTS(:,166).*CONSTANTS(:,101))./(STATES(:,127)+ CONSTANTS(:,166).*STATES(:,49))); ALGEBRAIC(:,260) = 0.195610.*(1.00000+0.600000./(1.00000+power(3.80000e-05./ALGEBRAIC(:,187), 1.40000))); ALGEBRAIC(:,350) = ( ALGEBRAIC(:,260).*(STATES(:,70)+STATES(:,71))).*(STATES(:,51) - ALGEBRAIC(:,349)); ALGEBRAIC(:,351) = ( ALGEBRAIC(:,260).*(STATES(:,72)+STATES(:,73))).*(STATES(:,51) - ALGEBRAIC(:,349)); ALGEBRAIC(:,263) = (STATES(:,104)+CONSTANTS(:,404))./CONSTANTS(:,111); ALGEBRAIC(:,265) = (ALGEBRAIC(:,263) - CONSTANTS(:,405))./(0.785000 - CONSTANTS(:,405)); ALGEBRAIC(:,267) = piecewise({ALGEBRAIC(:,265)<0.00000, 0.00000 }, ALGEBRAIC(:,265)); ALGEBRAIC(:,352) = ALGEBRAIC(:,267).*ALGEBRAIC(:,351)+ (1.00000 - ALGEBRAIC(:,267)).*ALGEBRAIC(:,350); ALGEBRAIC(:,342) = ( CONSTANTS(:,124).*(STATES(:,51) - ALGEBRAIC(:,332)))./(1.00000+exp((15.0000 - STATES(:,51))./17.0000)); ALGEBRAIC(:,343) = (( CONSTANTS(:,124).*(STATES(:,51) - ALGEBRAIC(:,332)))./(1.00000+exp((36.0000 - STATES(:,51))./17.0000))).*3.62000; ALGEBRAIC(:,269) = (STATES(:,105) - 0.0589380)./(0.393747 - 0.0589380); ALGEBRAIC(:,271) = piecewise({ALGEBRAIC(:,269)<0.00000, 0.00000 }, ALGEBRAIC(:,269)); ALGEBRAIC(:,344) = (1.00000 - ALGEBRAIC(:,271)).*ALGEBRAIC(:,342)+ ALGEBRAIC(:,271).*ALGEBRAIC(:,343); ALGEBRAIC(:,326) = exp(STATES(:,51)./550.000); ALGEBRAIC(:,345) = ( ( CONSTANTS(:,158).*power(STATES(:,121), 3.00000)).*ALGEBRAIC(:,326)).*(STATES(:,51) - ALGEBRAIC(:,332)); ALGEBRAIC(:,346) = ALGEBRAIC(:,345).*( 0.735600.*STATES(:,122)+ 0.264400.*STATES(:,124)); ALGEBRAIC(:,347) = ALGEBRAIC(:,345).*( 0.735600.*STATES(:,123)+ 0.264400.*STATES(:,125)); ALGEBRAIC(:,348) = (1.00000 - STATES(:,40)).*ALGEBRAIC(:,347)+ STATES(:,40).*ALGEBRAIC(:,346); ALGEBRAIC(:,365) = ((((ALGEBRAIC(:,340)+ALGEBRAIC(:,341))+ALGEBRAIC(:,352))+ALGEBRAIC(:,344))+ALGEBRAIC(:,348)) - 2.00000.*ALGEBRAIC(:,304); ALGEBRAIC(:,367) = ALGEBRAIC(:,366)+ ALGEBRAIC(:,292).*3.00000; ALGEBRAIC(:,368) = ((ALGEBRAIC(:,367)+ALGEBRAIC(:,365))+ALGEBRAIC(:,306))+ALGEBRAIC(:,331); ALGEBRAIC(:,369) = piecewise({(VOI - CONSTANTS(:,214)) - CONSTANTS(:,215).*floor((VOI - CONSTANTS(:,214))./CONSTANTS(:,215))<CONSTANTS(:,213), 1.00000 }, 0.00000).*CONSTANTS(:,257); 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