C
C There are a total of 119 entries in the algebraic variable array.
C There are a total of 31 entries in each of the rate and state variable arrays.
C There are a total of 130 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is L_totmax in component b1_AR_Gs_parameters (uM).
C CONSTS(2) is sum_b1_AR in component b1_AR_Gs_parameters (uM).
C CONSTS(3) is Gs_tot in component b1_AR_Gs_parameters (uM).
C CONSTS(4) is Kl in component b1_AR_Gs_parameters (uM).
C CONSTS(5) is Kr in component b1_AR_Gs_parameters (uM).
C CONSTS(6) is Kc in component b1_AR_Gs_parameters (uM).
C CONSTS(7) is k_bar_kp in component b1_AR_Gs_parameters (per_sec).
C CONSTS(8) is k_bar_km in component b1_AR_Gs_parameters (per_sec).
C CONSTS(9) is k_p_kap in component b1_AR_Gs_parameters (per_uM_per_sec).
C CONSTS(10) is k_p_kam in component b1_AR_Gs_parameters (per_sec).
C CONSTS(11) is k_g_act in component b1_AR_Gs_parameters (per_sec).
C CONSTS(12) is k_hyd in component b1_AR_Gs_parameters (per_sec).
C CONSTS(13) is k_reassoc in component b1_AR_Gs_parameters (per_uM_per_sec).
C CONSTS(14) is AC_tot in component cAMP_parameters (uM).
C CONSTS(15) is ATP in component cAMP_parameters (uM).
C CONSTS(16) is PDE_tot in component cAMP_parameters (uM).
C CONSTS(17) is IBMX_tot in component cAMP_parameters (uM).
C CONSTS(18) is Fsk_tot in component cAMP_parameters (uM).
C CONSTS(19) is k_ac_basal in component cAMP_parameters (per_sec).
C CONSTS(20) is k_ac_gsa in component cAMP_parameters (per_sec).
C CONSTS(21) is k_ac_fsk in component cAMP_parameters (per_sec).
C CONSTS(22) is k_pde in component cAMP_parameters (per_sec).
C CONSTS(23) is Km_basal in component cAMP_parameters (uM).
C CONSTS(24) is Km_gsa in component cAMP_parameters (uM).
C CONSTS(25) is Km_fsk in component cAMP_parameters (uM).
C CONSTS(26) is Km_pde in component cAMP_parameters (uM).
C CONSTS(27) is K_gsa in component cAMP_parameters (uM).
C CONSTS(28) is K_fsk in component cAMP_parameters (uM).
C CONSTS(29) is Ki_ibmx in component cAMP_parameters (uM).
C CONSTS(30) is PKAI_tot in component PKA_parameters (uM).
C CONSTS(31) is PKAII_tot in component PKA_parameters (uM).
C CONSTS(32) is PKI_tot in component PKA_parameters (uM).
C CONSTS(33) is K_a in component PKA_parameters (uM).
C CONSTS(34) is K_b in component PKA_parameters (uM).
C CONSTS(35) is K_d in component PKA_parameters (uM).
C CONSTS(36) is Ki_pki in component PKA_parameters (uM).
C CONSTS(37) is PLB_tot in component PLB_parameters (uM).
C CONSTS(38) is PP1_tot in component PLB_parameters (uM).
C CONSTS(39) is Inhib1_tot in component PLB_parameters (uM).
C CONSTS(40) is k_pka_plb in component PLB_parameters (per_sec).
C CONSTS(41) is Km_pka_plb in component PLB_parameters (uM).
C CONSTS(42) is k_pp1_plb in component PLB_parameters (per_sec).
C CONSTS(43) is Km_pp1_plb in component PLB_parameters (uM).
C CONSTS(44) is k_pka_i1 in component PLB_parameters (per_sec).
C CONSTS(45) is Km_pka_i1 in component PLB_parameters (uM).
C CONSTS(46) is Vmax_pp2a_i1 in component PLB_parameters (uM_per_sec).
C CONSTS(47) is Km_pp2a_i1 in component PLB_parameters (uM).
C CONSTS(48) is Ki_inhib1 in component PLB_parameters (uM).
C CONSTS(49) is epsilon in component LCC_parameters (dimensionless).
C CONSTS(50) is LCC_tot in component LCC_parameters (uM).
C CONSTS(51) is PP1_lcc_tot in component LCC_parameters (uM).
C CONSTS(52) is PP2A_lcc_tot in component LCC_parameters (uM).
C CONSTS(53) is k_pka_lcc in component LCC_parameters (per_sec).
C CONSTS(54) is Km_pka_lcc in component LCC_parameters (uM).
C CONSTS(55) is k_pp1_lcc in component LCC_parameters (per_sec).
C CONSTS(56) is Km_pp1_lcc in component LCC_parameters (uM).
C CONSTS(57) is k_pp2a_lcc in component LCC_parameters (per_sec).
C CONSTS(58) is Km_pp2a_lcc in component LCC_parameters (uM).
C CONSTS(59) is V_myo in component EC_Coupling_Parameters (uL).
C CONSTS(60) is V_nsr in component EC_Coupling_Parameters (uL).
C CONSTS(61) is V_jsr in component EC_Coupling_Parameters (uL).
C CONSTS(62) is A_Cap in component EC_Coupling_Parameters (cm2).
C CONSTS(63) is Temp in component EC_Coupling_Parameters (kelvin).
C CONSTS(64) is Na_ext in component EC_Coupling_Parameters (mM).
C CONSTS(65) is K_ext in component EC_Coupling_Parameters (mM).
C CONSTS(66) is Ca_ext in component EC_Coupling_Parameters (mM).
C CONSTS(67) is G_Na in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(68) is G_to in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(69) is G_ss in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(70) is G_Ki_bar in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(71) is G_Kp in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(72) is f in component EC_Coupling_Parameters (per_sec).
C CONSTS(73) is g in component EC_Coupling_Parameters (per_sec).
C CONSTS(74) is gamma_o in component EC_Coupling_Parameters (per_mM_per_sec).
C CONSTS(75) is omega in component EC_Coupling_Parameters (per_sec).
C CONSTS(76) is p_Ca in component EC_Coupling_Parameters (cm_per_sec).
C CONSTS(77) is p_K in component EC_Coupling_Parameters (cm_per_sec).
C CONSTS(78) is N_lcc in component EC_Coupling_Parameters (dimensionless).
C CONSTS(79) is I_Ca05 in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(80) is k_NaCa in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(81) is Km_Na in component EC_Coupling_Parameters (mM).
C CONSTS(82) is Km_Ca in component EC_Coupling_Parameters (mM).
C CONSTS(83) is k_sat in component EC_Coupling_Parameters (dimensionless).
C CONSTS(84) is eta in component EC_Coupling_Parameters (dimensionless).
C CONSTS(85) is I_bar_NaK in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(86) is Km_Nai in component EC_Coupling_Parameters (mM).
C CONSTS(87) is Km_Ko in component EC_Coupling_Parameters (mM).
C CONSTS(88) is I_bar_PCa in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(89) is Km_PCa in component EC_Coupling_Parameters (mM).
C CONSTS(90) is G_CaB in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(91) is G_NaB in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(92) is Pns in component EC_Coupling_Parameters (dimensionless).
C CONSTS(93) is Km_NS in component EC_Coupling_Parameters (mM).
C CONSTS(94) is I_up_bar in component EC_Coupling_Parameters (mM_per_sec).
C CONSTS(95) is Km_up0 in component EC_Coupling_Parameters (mM).
C CONSTS(96) is NSR_bar in component EC_Coupling_Parameters (mM).
C CONSTS(97) is tau_on in component EC_Coupling_Parameters (second).
C CONSTS(98) is tau_off in component EC_Coupling_Parameters (second).
C CONSTS(99) is G_max_rel in component EC_Coupling_Parameters (mM_per_sec).
C CONSTS(100) is d_Cai_th in component EC_Coupling_Parameters (mM).
C CONSTS(101) is Km_rel in component EC_Coupling_Parameters (mM).
C CONSTS(102) is CSQN_th in component EC_Coupling_Parameters (mM).
C CONSTS(103) is CSQN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(104) is Km_CSQN in component EC_Coupling_Parameters (mM).
C CONSTS(105) is tau_tr in component EC_Coupling_Parameters (second).
C CONSTS(106) is TRPN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(107) is CMDN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(108) is INDO_bar in component EC_Coupling_Parameters (mM).
C CONSTS(109) is Km_TRPN in component EC_Coupling_Parameters (mM).
C CONSTS(110) is Km_CMDN in component EC_Coupling_Parameters (mM).
C CONSTS(111) is Km_INDO in component EC_Coupling_Parameters (mM).
C ALGBRC(63) is LR in component b1_AR_module (uM).
C ALGBRC(64) is LRG in component b1_AR_module (uM).
C ALGBRC(65) is RG in component b1_AR_module (uM).
C ALGBRC(81) is BARK_DESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(12) is BARK_RESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(89) is PKA_DESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(24) is PKA_RESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(82) is G_ACT in component b1_AR_module (uM_per_sec).
C ALGBRC(26) is HYD in component b1_AR_module (uM_per_sec).
C ALGBRC(28) is REASSOC in component b1_AR_module (per_sec).
C ALGBRC(66) is L in component b1_AR_module (uM).
C ALGBRC(67) is R in component b1_AR_module (uM).
C ALGBRC(68) is Gs in component b1_AR_module (uM).
C STATES(1) is b1_AR_tot in component b1_AR_module (uM).
C STATES(2) is b1_AR_d in component b1_AR_module (uM).
C STATES(3) is b1_AR_p in component b1_AR_module (uM).
C STATES(4) is Gs_agtp_tot in component b1_AR_module (uM).
C STATES(5) is Gs_agdp in component b1_AR_module (uM).
C STATES(6) is Gs_bg in component b1_AR_module (uM).
C ALGBRC(69) is PKAC_I in component PKA_module (uM).
C ALGBRC(70) is cAMP in component PKA_module (uM).
C ALGBRC(43) is Gsa_GTP in component cAMP_module (uM).
C ALGBRC(51) is Fsk in component cAMP_module (uM).
C ALGBRC(44) is AC in component cAMP_module (uM).
C CONSTS(128) is PDE in component cAMP_module (uM).
C CONSTS(129) is IBMX in component cAMP_module (uM).
C STATES(7) is cAMP_tot in component cAMP_module (uM).
C ALGBRC(45) is Gsa_GTP_AC in component cAMP_module (uM).
C ALGBRC(52) is Fsk_AC in component cAMP_module (uM).
C ALGBRC(49) is AC_ACT_GSA in component cAMP_module (uM).
C ALGBRC(47) is AC_ACT_BASAL in component cAMP_module (uM).
C ALGBRC(54) is AC_ACT_FSK in component cAMP_module (uM).
C ALGBRC(83) is PDE_ACT in component cAMP_module (uM).
C CONSTS(130) is PDE_IBMX in component cAMP_module (uM).
C ALGBRC(71) is PKI in component PKA_module (uM).
C ALGBRC(72) is A2RC_I in component PKA_module (uM).
C ALGBRC(73) is A2R_I in component PKA_module (uM).
C ALGBRC(74) is A2RC_II in component PKA_module (uM).
C ALGBRC(75) is A2R_II in component PKA_module (uM).
C ALGBRC(76) is ARC_I in component PKA_module (uM).
C ALGBRC(77) is ARC_II in component PKA_module (uM).
C ALGBRC(78) is PKA_temp in component PKA_module (uM).
C ALGBRC(79) is PKAC_II in component PKA_module (uM).
C ALGBRC(30) is PLB in component PLB_module (uM).
C ALGBRC(84) is PLB_PHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(61) is PLB_DEPHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(31) is Inhib1 in component PLB_module (uM).
C ALGBRC(56) is Inhib1p_PP1 in component PLB_module (uM).
C ALGBRC(85) is Inhib1_PHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(33) is Inhib1_DEPHOSPH in component PLB_module (uM_per_sec).
C STATES(8) is PLB_p in component PLB_module (uM).
C STATES(9) is Inhib1p_tot in component PLB_module (uM).
C ALGBRC(57) is Inhib1p in component PLB_module (uM).
C ALGBRC(58) is PP1 in component PLB_module (uM).
C ALGBRC(1) is frac_PLB_p in component PLB_module (dimensionless).
C ALGBRC(32) is frac_PLB in component PLB_module (dimensionless).
C CONSTS(120) is frac_PLB_o in component PLB_module (dimensionless).
C ALGBRC(35) is LCCa in component LCC_module (uM).
C ALGBRC(86) is LCCa_PHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(37) is LCCa_DEPHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(39) is LCCb in component LCC_module (uM).
C ALGBRC(87) is LCCb_PHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(41) is LCCb_DEPHOSPH in component LCC_module (uM_per_sec).
C STATES(10) is LCCa_p in component LCC_module (uM).
C STATES(11) is LCCb_p in component LCC_module (uM).
C ALGBRC(2) is frac_LCCa_p in component LCC_module (dimensionless).
C CONSTS(121) is frac_LCCa_po in component LCC_module (dimensionless).
C ALGBRC(34) is frac_LCCb_p in component LCC_module (dimensionless).
C CONSTS(122) is frac_LCCb_po in component LCC_module (dimensionless).
C ALGBRC(36) is E_Na in component Nernst_Potentials (mV).
C ALGBRC(38) is E_K in component Nernst_Potentials (mV).
C ALGBRC(40) is E_Ca in component Nernst_Potentials (mV).
C CONSTS(123) is E_Cl in component Nernst_Potentials (mV).
C CONSTS(112) is R in component Nernst_Potentials (joules_per_mole_kelvin).
C CONSTS(113) is Frdy in component Nernst_Potentials (coulombs_per_mole).
C CONSTS(124) is FoRT in component Nernst_Potentials (per_mV).
C CONSTS(114) is z_Na in component Nernst_Potentials (dimensionless).
C CONSTS(115) is z_K in component Nernst_Potentials (dimensionless).
C CONSTS(116) is z_Ca in component Nernst_Potentials (dimensionless).
C STATES(12) is Na_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C STATES(13) is K_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C STATES(14) is Ca_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C ALGBRC(3) is am in component Fast_Na_Current (per_sec).
C ALGBRC(14) is bm in component Fast_Na_Current (per_sec).
C ALGBRC(4) is ah in component Fast_Na_Current (per_sec).
C ALGBRC(5) is aj in component Fast_Na_Current (per_sec).
C ALGBRC(15) is bh in component Fast_Na_Current (per_sec).
C ALGBRC(16) is bj in component Fast_Na_Current (per_sec).
C STATES(15) is m in component Fast_Na_Current (dimensionless).
C STATES(16) is h in component Fast_Na_Current (dimensionless).
C STATES(17) is j in component Fast_Na_Current (dimensionless).
C ALGBRC(42) is I_Na in component Fast_Na_Current (uA_per_uF).
C STATES(18) is V_m in component Ion_Concentrations_and_Membrane_Potential (mV).
C ALGBRC(6) is a_lcc in component L_Type_Calcium_Current (per_sec).
C ALGBRC(17) is b_lcc in component L_Type_Calcium_Current (per_sec).
C ALGBRC(19) is f_lcc in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(7) is y_lcc_inf in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(18) is tau_y_lcc in component L_Type_Calcium_Current (second).
C ALGBRC(25) is gamma in component L_Type_Calcium_Current (per_mM_per_sec).
C ALGBRC(27) is v_gamma in component L_Type_Calcium_Current (per_mM_per_sec).
C ALGBRC(29) is v_omega in component L_Type_Calcium_Current (dimensionless).
C STATES(19) is v in component L_Type_Calcium_Current (dimensionless).
C STATES(20) is w in component L_Type_Calcium_Current (dimensionless).
C STATES(21) is x in component L_Type_Calcium_Current (dimensionless).
C STATES(22) is y in component L_Type_Calcium_Current (dimensionless).
C STATES(23) is z in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(46) is i_bar_Ca in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(48) is i_bar_K in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(50) is f_avail in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(53) is I_Ca in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(55) is I_CaK in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(59) is I_Ca_tot in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(8) is r_toss in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(9) is s_toss in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(20) is tau_r_to in component Transient_Outward_K_Current (second).
C ALGBRC(21) is tau_s_to in component Transient_Outward_K_Current (second).
C ALGBRC(22) is tau_ss_to in component Transient_Outward_K_Current (second).
C STATES(24) is r_to in component Transient_Outward_K_Current (dimensionless).
C STATES(25) is s_to in component Transient_Outward_K_Current (dimensionless).
C STATES(26) is ss_to in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(60) is I_to in component Transient_Outward_K_Current (uA_per_uF).
C ALGBRC(10) is r_ss_inf in component Steady_State_K_Current (dimensionless).
C ALGBRC(23) is tau_r_ss in component Steady_State_K_Current (second).
C ALGBRC(11) is s_ss_inf in component Steady_State_K_Current (dimensionless).
C CONSTS(125) is tau_s_ss in component Steady_State_K_Current (second).
C STATES(27) is r_ss in component Steady_State_K_Current (dimensionless).
C STATES(28) is s_ss in component Steady_State_K_Current (dimensionless).
C ALGBRC(62) is I_ss in component Steady_State_K_Current (uA_per_uF).
C ALGBRC(80) is a_Ki in component Time_Independent_K_Current (dimensionless).
C ALGBRC(88) is b_Ki in component Time_Independent_K_Current (dimensionless).
C ALGBRC(90) is Ki_ss in component Time_Independent_K_Current (dimensionless).
C ALGBRC(91) is I_Ki in component Time_Independent_K_Current (uA_per_uF).
C ALGBRC(92) is Kp in component Plateau_K_Current (dimensionless).
C ALGBRC(93) is I_Kp in component Plateau_K_Current (uA_per_uF).
C ALGBRC(94) is s4 in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(95) is s5 in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(96) is I_NCX in component Na_Ca_Exchanger_Current (uA_per_uF).
C CONSTS(126) is sigma in component Na_K_Pump_Current (dimensionless).
C ALGBRC(97) is f_NaK in component Na_K_Pump_Current (dimensionless).
C ALGBRC(98) is I_NaK in component Na_K_Pump_Current (uA_per_uF).
C ALGBRC(99) is I_PCa in component Sarcolemmal_Ca_Pump_Current (uA_per_uF).
C ALGBRC(100) is I_CaB in component Ca_Background_Current (uA_per_uF).
C ALGBRC(101) is I_NaB in component Na_Background_Current (uA_per_uF).
C ALGBRC(102) is I_Na_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(103) is I_K_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(104) is I_Ca_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(105) is t_rel in component Calcium_Induced_Calcium_Release (second).
C ALGBRC(106) is ryr_on in component Calcium_Induced_Calcium_Release (dimensionless).
C ALGBRC(107) is ryr_off in component Calcium_Induced_Calcium_Release (dimensionless).
C ALGBRC(108) is g_rel in component Calcium_Induced_Calcium_Release (per_sec).
C ALGBRC(109) is I_rel in component Calcium_Induced_Calcium_Release (mM_per_sec).
C STATES(29) is Ca_jsr in component Other_SR_Fluxes_and_Concentrations (mM).
C STATES(30) is trel in component Ion_Concentrations_and_Membrane_Potential (second).
C ALGBRC(110) is Km_up in component Other_SR_Fluxes_and_Concentrations (mM).
C ALGBRC(111) is I_up in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(112) is I_leak in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(113) is I_tr in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(115) is B_jsr in component Other_SR_Fluxes_and_Concentrations (dimensionless).
C STATES(31) is Ca_nsr in component Other_SR_Fluxes_and_Concentrations (mM).
C ALGBRC(117) is SR_content in component Other_SR_Fluxes_and_Concentrations (uM).
C ALGBRC(114) is b_trpn in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(116) is b_cmdn in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(118) is b_indo in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(119) is B_myo in component Cytoplasmic_Calcium_Buffering (dimensionless).
C CONSTS(127) is I_app in component Ion_Concentrations_and_Membrane_Potential (uA_per_uF).
C CONSTS(117) is V_hold in component Ion_Concentrations_and_Membrane_Potential (mV).
C CONSTS(118) is V_test in component Ion_Concentrations_and_Membrane_Potential (mV).
C ALGBRC(13) is V_clamp in component Ion_Concentrations_and_Membrane_Potential (mV).
C CONSTS(119) is R_clamp in component Ion_Concentrations_and_Membrane_Potential (dimensionless).
C RATES(1) is d/dt b1_AR_tot in component b1_AR_module (uM).
C RATES(2) is d/dt b1_AR_d in component b1_AR_module (uM).
C RATES(3) is d/dt b1_AR_p in component b1_AR_module (uM).
C RATES(4) is d/dt Gs_agtp_tot in component b1_AR_module (uM).
C RATES(5) is d/dt Gs_agdp in component b1_AR_module (uM).
C RATES(6) is d/dt Gs_bg in component b1_AR_module (uM).
C RATES(7) is d/dt cAMP_tot in component cAMP_module (uM).
C RATES(8) is d/dt PLB_p in component PLB_module (uM).
C RATES(9) is d/dt Inhib1p_tot in component PLB_module (uM).
C RATES(10) is d/dt LCCa_p in component LCC_module (uM).
C RATES(11) is d/dt LCCb_p in component LCC_module (uM).
C RATES(15) is d/dt m in component Fast_Na_Current (dimensionless).
C RATES(16) is d/dt h in component Fast_Na_Current (dimensionless).
C RATES(17) is d/dt j in component Fast_Na_Current (dimensionless).
C RATES(19) is d/dt v in component L_Type_Calcium_Current (dimensionless).
C RATES(20) is d/dt w in component L_Type_Calcium_Current (dimensionless).
C RATES(21) is d/dt x in component L_Type_Calcium_Current (dimensionless).
C RATES(22) is d/dt y in component L_Type_Calcium_Current (dimensionless).
C RATES(23) is d/dt z in component L_Type_Calcium_Current (dimensionless).
C RATES(24) is d/dt r_to in component Transient_Outward_K_Current (dimensionless).
C RATES(25) is d/dt s_to in component Transient_Outward_K_Current (dimensionless).
C RATES(26) is d/dt ss_to in component Transient_Outward_K_Current (dimensionless).
C RATES(27) is d/dt r_ss in component Steady_State_K_Current (dimensionless).
C RATES(28) is d/dt s_ss in component Steady_State_K_Current (dimensionless).
C RATES(31) is d/dt Ca_nsr in component Other_SR_Fluxes_and_Concentrations (mM).
C RATES(29) is d/dt Ca_jsr in component Other_SR_Fluxes_and_Concentrations (mM).
C RATES(12) is d/dt Na_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C RATES(13) is d/dt K_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C RATES(14) is d/dt Ca_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C RATES(18) is d/dt V_m in component Ion_Concentrations_and_Membrane_Potential (mV).
C RATES(30) is d/dt trel in component Ion_Concentrations_and_Membrane_Potential (second).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      CONSTS(2) = 0.0132
      CONSTS(3) = 3.83
      CONSTS(4) = 0.285
      CONSTS(5) = 0.062
      CONSTS(6) = 33
      CONSTS(7) = 1.1e-3
      CONSTS(8) = 2.2e-3
      CONSTS(9) = 3.6e-3
      CONSTS(10) = 2.2e-3
      CONSTS(11) = 16
      CONSTS(12) = 0.8
      CONSTS(13) = 1.21e3
      CONSTS(14) = 49.7e-3
      CONSTS(15) = 5e3
      CONSTS(16) = 38.9e-3
      CONSTS(17) = 0
      CONSTS(18) = 0
      CONSTS(19) = 0.2
      CONSTS(20) = 8.5
      CONSTS(21) = 7.3
      CONSTS(22) = 5
      CONSTS(23) = 1.03e3
      CONSTS(24) = 315
      CONSTS(25) = 860
      CONSTS(26) = 1.3
      CONSTS(27) = 0.4
      CONSTS(28) = 44
      CONSTS(29) = 30
      CONSTS(30) = 0.59
      CONSTS(31) = 0.025
      CONSTS(32) = 0.18
      CONSTS(33) = 9.14
      CONSTS(34) = 1.64
      CONSTS(35) = 4.375
      CONSTS(36) = 0.2e-3
      CONSTS(37) = 106
      CONSTS(38) = 0.89
      CONSTS(39) = 0.3
      CONSTS(40) = 54
      CONSTS(41) = 21
      CONSTS(42) = 8.5
      CONSTS(43) = 7
      CONSTS(44) = 60
      CONSTS(45) = 1
      CONSTS(46) = 14
      CONSTS(47) = 1
      CONSTS(48) = 1e-3
      CONSTS(49) = 10
      CONSTS(50) = 25e-3
      CONSTS(51) = 25e-3
      CONSTS(52) = 25e-3
      CONSTS(53) = 54
      CONSTS(54) = 21
      CONSTS(55) = 8.52
      CONSTS(56) = 3
      CONSTS(57) = 10.1
      CONSTS(58) = 3
      CONSTS(59) = 20.8e-6
      CONSTS(60) = 9.88e-7
      CONSTS(61) = 9.3e-8
      CONSTS(62) = 1.534e-4
      CONSTS(63) = 310
      CONSTS(64) = 140
      CONSTS(65) = 5.4
      CONSTS(66) = 1.8
      CONSTS(67) = 8
      CONSTS(68) = 0.35
      CONSTS(69) = 0.07
      CONSTS(70) = 0.24
      CONSTS(71) = 0.008
      CONSTS(72) = 300
      CONSTS(73) = 2e3
      CONSTS(74) = 5187.5
      CONSTS(75) = 10
      CONSTS(76) = 1.7469e-8
      CONSTS(77) = 3.234e-11
      CONSTS(78) = 3e5
      CONSTS(79) = -0.458
      CONSTS(80) = 1483
      CONSTS(81) = 87.5
      CONSTS(82) = 1.38
      CONSTS(83) = 0.1
      CONSTS(84) = 0.35
      CONSTS(85) = 1.1
      CONSTS(86) = 10
      CONSTS(87) = 1.5
      CONSTS(88) = 1.15
      CONSTS(89) = 0.5e-3
      CONSTS(90) = 2.8e-3
      CONSTS(91) = 1.18e-3
      CONSTS(92) = 0
      CONSTS(93) = 1.2e-3
      CONSTS(94) = 4.7
      CONSTS(95) = 3e-4
      CONSTS(96) = 15
      CONSTS(97) = 2e-3
      CONSTS(98) = 2e-3
      CONSTS(99) = 60e3
      CONSTS(100) = 0.18e-3
      CONSTS(101) = 0.8e-3
      CONSTS(102) = 8.75
      CONSTS(103) = 15
      CONSTS(104) = 0.8
      CONSTS(105) = 5.7e-4
      CONSTS(106) = 0.07
      CONSTS(107) = 0.05
      CONSTS(108) = 0.07
      CONSTS(109) = 0.5128e-3
      CONSTS(110) = 2.38e-3
      CONSTS(111) = 8.44e-4
      STATES(1) = 0.01205
      STATES(2) = 0
      STATES(3) = 1.154e-3
      STATES(4) = 0.02505
      STATES(5) = 6.446e-4
      STATES(6) = 0.02569
      STATES(7) = 0.8453
      STATES(8) = 4.105
      STATES(9) = 0.0526
      STATES(10) = 5.103e-3
      STATES(11) = 5.841e-3
      CONSTS(112) = 8314
      CONSTS(113) = 96485
      CONSTS(114) = 1
      CONSTS(115) = 1
      CONSTS(116) = 2
      STATES(12) = 16
      STATES(13) = 145
      STATES(14) = 1.58e-4
      STATES(15) = 1.4e-3
      STATES(16) = 0.99
      STATES(17) = 0.99
      STATES(18) = -85.66
      STATES(19) = 0
      STATES(20) = 0
      STATES(21) = 0.13
      STATES(22) = 0.96
      STATES(23) = 0.92
      STATES(24) = 1.4e-3
      STATES(25) = 1
      STATES(26) = 0.613
      STATES(27) = 198e-3
      STATES(28) = 0.43
      STATES(29) = 1.92
      STATES(30) = 0.9
      STATES(31) = 1.92
      CONSTS(117) = -40
      CONSTS(118) = -10
      CONSTS(119) = 0.02
      CONSTS(120) = 0.961300
      CONSTS(121) = 0.204100
      CONSTS(122) = 0.233600
      CONSTS(123) = - 40.0000
      CONSTS(124) = (CONSTS(113)/CONSTS(112))/CONSTS(63)
      CONSTS(125) = 2.10000
      CONSTS(126) = (EXP(CONSTS(64)/67.3000) - 1.00000)/7.00000
      CONSTS(127) = 0.00000
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(18)+87.5000)/10.3000))
      RATES(28) = (ALGBRC(11) - STATES(28))/CONSTS(125)
      ALGBRC(3) = ( 0.320000*(STATES(18)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(18)+47.1300)))
      ALGBRC(14) =  0.0800000*EXP(- STATES(18)/11.0000)
      RATES(15) =  1000.00*( ALGBRC(3)*(1.00000 - STATES(15)) -  ALGBRC(14)*STATES(15))
      ALGBRC(4) = TERNRY(STATES(18).GE.- 40.0000, 0.00000,  0.135000*EXP((80.0000+STATES(18))/- 6.80000))
      ALGBRC(15) = TERNRY(STATES(18).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP(- (STATES(18)+10.6600)/11.1000))),  3.56000*EXP( 0.0790000*STATES(18))+ 310000.*EXP( 0.350000*STATES(18)))
      RATES(16) =  1000.00*( ALGBRC(4)*(1.00000 - STATES(16)) -  ALGBRC(15)*STATES(16))
      ALGBRC(5) = TERNRY(STATES(18).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(18)) -  3.47400e-05*EXP( - 0.0439100*STATES(18)))*(STATES(18)+37.7800))/(1.00000+EXP( 0.311000*(STATES(18)+79.2300))))
      ALGBRC(16) = TERNRY(STATES(18).GE.- 40.0000, ( 0.300000*EXP( - 2.57500e-07*STATES(18)))/(1.00000+EXP( - 0.100000*(STATES(18)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(18)))/(1.00000+EXP( - 0.137800*(STATES(18)+40.1400))))
      RATES(17) =  1000.00*( ALGBRC(5)*(1.00000 - STATES(17)) -  ALGBRC(16)*STATES(17))
      ALGBRC(6) =  400.000*EXP((STATES(18)+2.00000)/10.0000)
      ALGBRC(17) =  50.0000*EXP(( - 1.00000*(STATES(18)+2.00000))/13.0000)
      RATES(19) =  ALGBRC(6)*(1.00000 - STATES(19)) -  ALGBRC(17)*STATES(19)
      RATES(20) =  2.00000*ALGBRC(6)*(1.00000 - STATES(20)) -  (ALGBRC(17)/2.00000)*STATES(20)
      ALGBRC(2) = STATES(10)/CONSTS(50)
      ALGBRC(19) =  CONSTS(72)*(( 0.375000*ALGBRC(2))/CONSTS(121)+0.625000)
      RATES(21) =  ALGBRC(19)*(1.00000 - STATES(21)) -  CONSTS(73)*STATES(21)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(18)+55.0000)/7.50000))+0.100000/(1.00000+EXP((- STATES(18)+21.0000)/6.00000))
      ALGBRC(18) = 0.0200000+0.300000/(1.00000+EXP((STATES(18)+30.0000)/9.50000))
      RATES(22) = (ALGBRC(7) - STATES(22))/ALGBRC(18)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(18)+10.6000)/- 11.4200))
      ALGBRC(20) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
      RATES(24) = (ALGBRC(8) - STATES(24))/ALGBRC(20)
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(18)+43.5000)/6.88410))
      ALGBRC(21) =  0.350000*EXP( - 1.00000*(STATES(18)+70.0000)/15.0000 ** 2.00000)+0.0350000
      RATES(25) = (ALGBRC(9) - STATES(25))/ALGBRC(21)
      ALGBRC(22) =  3.70000*EXP( - 1.00000*(STATES(18)+70.0000)/30.0000 ** 2.00000)+0.0350000
      RATES(26) = (ALGBRC(9) - STATES(26))/ALGBRC(22)
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(18)+11.5000)/- 11.8200))
      ALGBRC(23) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
      RATES(27) = (ALGBRC(10) - STATES(27))/ALGBRC(23)
      ALGBRC(26) =  CONSTS(12)*STATES(4)
      ALGBRC(28) =  CONSTS(13)*STATES(5)*STATES(6)
      RATES(5) = ALGBRC(26) - ALGBRC(28)
      ALGBRC(25) =  CONSTS(74)*STATES(14)
      ALGBRC(27) =  ALGBRC(25)*(1.00000 - STATES(19) ** 4.00000+ 2.00000*STATES(19)*1.00000 - STATES(19) ** 3.00000+ 4.00000*STATES(19) ** 2.00000*1.00000 - STATES(19) ** 2.00000+ 8.00000*STATES(19) ** 3.00000*(1.00000 - STATES(19))+ 16.0000*STATES(19) ** 4.00000*(1.00000 - ALGBRC(19)/CONSTS(73)))
      ALGBRC(29) =  CONSTS(75)*(1.00000 - STATES(20) ** 4.00000+ 0.500000*STATES(20)*1.00000 - STATES(20) ** 3.00000+ 0.250000*STATES(20) ** 2.00000*1.00000 - STATES(20) ** 2.00000+ 0.125000*STATES(20) ** 3.00000*(1.00000 - STATES(20))+ 0.0625000*STATES(20) ** 4.00000)
      RATES(23) =  ALGBRC(29)*(1.00000 - STATES(23)) -  ALGBRC(27)*STATES(23)
CALL minimize(minfunc_1, CONSTS, VARIABLES, <VAR>)
      ALGBRC(81) =  CONSTS(7)*(ALGBRC(63)+ALGBRC(64))
      ALGBRC(12) =  CONSTS(8)*STATES(2)
      RATES(2) = ALGBRC(81) - ALGBRC(12)
      ALGBRC(82) =  CONSTS(11)*(ALGBRC(65)+ALGBRC(64))
      RATES(4) = ALGBRC(82) - ALGBRC(26)
      RATES(6) = ALGBRC(82) - ALGBRC(28)
CALL minimize(minfunc_2, CONSTS, VARIABLES, <VAR>)
      ALGBRC(49) = ( CONSTS(20)*ALGBRC(45)*CONSTS(15))/(CONSTS(24)+CONSTS(15))
      ALGBRC(47) = ( CONSTS(19)*ALGBRC(44)*CONSTS(15))/(CONSTS(23)+CONSTS(15))
CALL minimize(minfunc_3, CONSTS, VARIABLES, <VAR>)
      ALGBRC(54) = ( CONSTS(21)*ALGBRC(52)*CONSTS(15))/(CONSTS(25)+CONSTS(15))
      ALGBRC(83) = ( CONSTS(22)*CONSTS(128)*ALGBRC(70))/(CONSTS(26)+ALGBRC(70))
      RATES(7) = (ALGBRC(47)+ALGBRC(49)+ALGBRC(54)) - ALGBRC(83)
      ALGBRC(30) = CONSTS(37) - STATES(8)
      ALGBRC(84) = ( CONSTS(40)*ALGBRC(69)*ALGBRC(30))/(CONSTS(41)+ALGBRC(30))
CALL minimize(minfunc_4, CONSTS, VARIABLES, <VAR>)
      ALGBRC(61) = ( CONSTS(42)*ALGBRC(58)*STATES(8))/(CONSTS(43)+STATES(8))
      RATES(8) = ALGBRC(84) - ALGBRC(61)
      ALGBRC(31) = CONSTS(39) - STATES(9)
      ALGBRC(85) = ( CONSTS(44)*ALGBRC(69)*ALGBRC(31))/(CONSTS(45)+ALGBRC(31))
      ALGBRC(33) = ( CONSTS(46)*STATES(9))/(CONSTS(47)+STATES(9))
      RATES(9) = ALGBRC(85) - ALGBRC(33)
      ALGBRC(35) = CONSTS(50) - STATES(10)
      ALGBRC(86) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(35))/(CONSTS(54)+ CONSTS(49)*ALGBRC(35))
      ALGBRC(37) = ( CONSTS(49)*CONSTS(57)*CONSTS(52)*STATES(10))/(CONSTS(58)+ CONSTS(49)*STATES(10))
      RATES(10) = ALGBRC(86) - ALGBRC(37)
      ALGBRC(39) = CONSTS(50) - STATES(11)
      ALGBRC(87) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(39))/(CONSTS(54)+ CONSTS(49)*ALGBRC(39))
      ALGBRC(41) = ( CONSTS(49)*CONSTS(55)*CONSTS(51)*STATES(11))/(CONSTS(56)+ CONSTS(49)*STATES(11))
      RATES(11) = ALGBRC(87) - ALGBRC(41)
      ALGBRC(89) =  CONSTS(9)*ALGBRC(69)*STATES(1)
      ALGBRC(24) =  CONSTS(10)*STATES(3)
      RATES(1) = (ALGBRC(12) - ALGBRC(81))+(ALGBRC(24) - ALGBRC(89))
      RATES(3) = ALGBRC(89) - ALGBRC(24)
      ALGBRC(36) =  (1.00000/( CONSTS(124)*CONSTS(114)))*log(CONSTS(64)/STATES(12))
      ALGBRC(42) =  CONSTS(67)*STATES(15) ** 3.00000*STATES(16)*STATES(17)*(STATES(18) - ALGBRC(36))
      ALGBRC(94) =  EXP( CONSTS(84)*STATES(18)*CONSTS(124))*STATES(12) ** 3.00000*CONSTS(66)
      ALGBRC(95) =  EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124))*CONSTS(64) ** 3.00000*STATES(14)
      ALGBRC(96) =  (CONSTS(80)/( (CONSTS(81) ** 3.00000+CONSTS(64) ** 3.00000)*(CONSTS(82)+CONSTS(66))*(1.00000+ CONSTS(83)*EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124)))))*(ALGBRC(94) - ALGBRC(95))
      ALGBRC(97) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(18)*CONSTS(124))+ 0.0365000*CONSTS(126)*EXP( - STATES(18)*CONSTS(124)))
      ALGBRC(98) = ( (( CONSTS(85)*ALGBRC(97))/(1.00000+CONSTS(86)/STATES(12) ** 1.50000))*CONSTS(65))/(CONSTS(65)+CONSTS(87))
      ALGBRC(101) =  CONSTS(91)*(STATES(18) - ALGBRC(36))
      ALGBRC(102) = ALGBRC(42)+ALGBRC(101)+ 3.00000*ALGBRC(96)+ 3.00000*ALGBRC(98)
      RATES(12) = ( - 1000.00*ALGBRC(102)*CONSTS(62))/( CONSTS(59)*CONSTS(114)*CONSTS(113))
      ALGBRC(48) = ( CONSTS(77)*STATES(18)*CONSTS(113)*CONSTS(124)*( STATES(13)*EXP( STATES(18)*CONSTS(124)) - CONSTS(65)))/(EXP( STATES(18)*CONSTS(124)) - 1.00000)
      ALGBRC(34) = STATES(11)/CONSTS(50)
      ALGBRC(50) =  0.500000*(( 0.400000*ALGBRC(34))/CONSTS(122)+0.600000)
      ALGBRC(46) = ( CONSTS(76)*4.00000*STATES(18)*CONSTS(113)*CONSTS(124)*( 0.00100000*EXP( 2.00000*STATES(18)*CONSTS(124)) -  0.341000*CONSTS(66)))/(EXP( 2.00000*STATES(18)*CONSTS(124)) - 1.00000)
      ALGBRC(53) =  ALGBRC(46)*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
      ALGBRC(55) =  (ALGBRC(48)/(1.00000+ALGBRC(53)/CONSTS(79)))*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
      ALGBRC(38) =  (1.00000/( CONSTS(124)*CONSTS(115)))*log(CONSTS(65)/STATES(13))
      ALGBRC(60) =  CONSTS(68)*STATES(24)*( 0.886000*STATES(25)+ 0.114000*STATES(26))*(STATES(18) - ALGBRC(38))
      ALGBRC(62) =  CONSTS(69)*STATES(27)*STATES(28)*(STATES(18) - ALGBRC(38))
      ALGBRC(80) = 1.02000/(1.00000+EXP( 0.238500*((STATES(18) - ALGBRC(38)) - 59.2150)))
      ALGBRC(88) = ( 0.491240*EXP( 0.0803200*((STATES(18)+5.47600) - ALGBRC(38)))+EXP( 0.0617500*((STATES(18) - ALGBRC(38)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(18) - ALGBRC(38))+4.75300)))
      ALGBRC(90) = ALGBRC(80)/(ALGBRC(80)+ALGBRC(88))
      ALGBRC(91) =  CONSTS(70)* (CONSTS(65)/5.40000) ** (1.0 / 2)*ALGBRC(90)*(STATES(18) - ALGBRC(38))
      ALGBRC(92) = 1.00000/(1.00000+EXP((7.48800 - STATES(18))/5.98000))
      ALGBRC(93) =  CONSTS(71)*ALGBRC(92)*(STATES(18) - ALGBRC(38))
      ALGBRC(103) = ((ALGBRC(60)+ALGBRC(62)+ALGBRC(91)+ALGBRC(93)) -  2.00000*ALGBRC(98))+ALGBRC(55)
      RATES(13) = ( - 1000.00*ALGBRC(103)*CONSTS(62))/( CONSTS(59)*CONSTS(115)*CONSTS(113))
      ALGBRC(99) = ( CONSTS(88)*STATES(14))/(CONSTS(89)+STATES(14))
      ALGBRC(40) =  (1.00000/( CONSTS(124)*CONSTS(116)))*log(CONSTS(66)/STATES(14))
      ALGBRC(100) =  CONSTS(90)*(STATES(18) - ALGBRC(40))
      ALGBRC(104) = (ALGBRC(53)+ALGBRC(100)+ALGBRC(99)) -  2.00000*ALGBRC(96)
      RATES(18) =  - 1000.00*((ALGBRC(104)+ALGBRC(103)+ALGBRC(102)) - CONSTS(127))
      RATES(30) = TERNRY( - 1000.00*((ALGBRC(104)+ALGBRC(103)+ALGBRC(102)) - CONSTS(127)).GT.30000.0, 1.00000 -  10000.0*STATES(30), 1.00000)
      ALGBRC(32) = ALGBRC(30)/CONSTS(37)
      ALGBRC(110) = ( CONSTS(95)*(1.00000+ 2.00000*ALGBRC(32)))/(1.00000+ 2.00000*CONSTS(120))
      ALGBRC(111) = ( CONSTS(94)*STATES(14) ** 2.00000)/(ALGBRC(110) ** 2.00000+STATES(14) ** 2.00000)
      ALGBRC(112) = ( CONSTS(94)*STATES(31))/CONSTS(96)
      ALGBRC(113) = (STATES(31) - STATES(29))/CONSTS(105)
      RATES(31) = (ALGBRC(111) - ALGBRC(112)) - ( ALGBRC(113)*CONSTS(61))/CONSTS(60)
      ALGBRC(105) = STATES(30)+0.00200000
      ALGBRC(106) = 1.00000 - EXP(- ALGBRC(105)/CONSTS(97))
      ALGBRC(107) = EXP(- ALGBRC(105)/CONSTS(98))
      ALGBRC(108) = CONSTS(99)/(1.00000+EXP((ALGBRC(104)+5.00000)/0.900000))
      ALGBRC(109) =  ALGBRC(108)*ALGBRC(106)*ALGBRC(107)*(STATES(29) - STATES(14))
      ALGBRC(115) = 1.00000/(1.00000+( CONSTS(103)*CONSTS(104))/CONSTS(104)+STATES(29) ** 2.00000)
      RATES(29) =  ALGBRC(115)*(ALGBRC(113) - ALGBRC(109))
      ALGBRC(114) = ( CONSTS(106)*CONSTS(109))/CONSTS(109)+STATES(14) ** 2.00000
      ALGBRC(116) = ( CONSTS(107)*CONSTS(110))/CONSTS(110)+STATES(14) ** 2.00000
      ALGBRC(118) = ( CONSTS(108)*CONSTS(111))/CONSTS(111)+STATES(14) ** 2.00000
      ALGBRC(119) = 1.00000/(1.00000+ALGBRC(116)+ALGBRC(114)+ALGBRC(114)+ALGBRC(118))
      RATES(14) =  - ALGBRC(119)*((( 1000.00*ALGBRC(104)*CONSTS(62))/( CONSTS(59)*CONSTS(116)*CONSTS(113))+( (ALGBRC(111) - ALGBRC(112))*CONSTS(60))/CONSTS(59)) - ( ALGBRC(109)*CONSTS(61))/CONSTS(59))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(18)+87.5000)/10.3000))
      ALGBRC(3) = ( 0.320000*(STATES(18)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(18)+47.1300)))
      ALGBRC(14) =  0.0800000*EXP(- STATES(18)/11.0000)
      ALGBRC(4) = TERNRY(STATES(18).GE.- 40.0000, 0.00000,  0.135000*EXP((80.0000+STATES(18))/- 6.80000))
      ALGBRC(15) = TERNRY(STATES(18).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP(- (STATES(18)+10.6600)/11.1000))),  3.56000*EXP( 0.0790000*STATES(18))+ 310000.*EXP( 0.350000*STATES(18)))
      ALGBRC(5) = TERNRY(STATES(18).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(18)) -  3.47400e-05*EXP( - 0.0439100*STATES(18)))*(STATES(18)+37.7800))/(1.00000+EXP( 0.311000*(STATES(18)+79.2300))))
      ALGBRC(16) = TERNRY(STATES(18).GE.- 40.0000, ( 0.300000*EXP( - 2.57500e-07*STATES(18)))/(1.00000+EXP( - 0.100000*(STATES(18)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(18)))/(1.00000+EXP( - 0.137800*(STATES(18)+40.1400))))
      ALGBRC(6) =  400.000*EXP((STATES(18)+2.00000)/10.0000)
      ALGBRC(17) =  50.0000*EXP(( - 1.00000*(STATES(18)+2.00000))/13.0000)
      ALGBRC(2) = STATES(10)/CONSTS(50)
      ALGBRC(19) =  CONSTS(72)*(( 0.375000*ALGBRC(2))/CONSTS(121)+0.625000)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(18)+55.0000)/7.50000))+0.100000/(1.00000+EXP((- STATES(18)+21.0000)/6.00000))
      ALGBRC(18) = 0.0200000+0.300000/(1.00000+EXP((STATES(18)+30.0000)/9.50000))
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(18)+10.6000)/- 11.4200))
      ALGBRC(20) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(18)+43.5000)/6.88410))
      ALGBRC(21) =  0.350000*EXP( - 1.00000*(STATES(18)+70.0000)/15.0000 ** 2.00000)+0.0350000
      ALGBRC(22) =  3.70000*EXP( - 1.00000*(STATES(18)+70.0000)/30.0000 ** 2.00000)+0.0350000
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(18)+11.5000)/- 11.8200))
      ALGBRC(23) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
      ALGBRC(26) =  CONSTS(12)*STATES(4)
      ALGBRC(28) =  CONSTS(13)*STATES(5)*STATES(6)
      ALGBRC(25) =  CONSTS(74)*STATES(14)
      ALGBRC(27) =  ALGBRC(25)*(1.00000 - STATES(19) ** 4.00000+ 2.00000*STATES(19)*1.00000 - STATES(19) ** 3.00000+ 4.00000*STATES(19) ** 2.00000*1.00000 - STATES(19) ** 2.00000+ 8.00000*STATES(19) ** 3.00000*(1.00000 - STATES(19))+ 16.0000*STATES(19) ** 4.00000*(1.00000 - ALGBRC(19)/CONSTS(73)))
      ALGBRC(29) =  CONSTS(75)*(1.00000 - STATES(20) ** 4.00000+ 0.500000*STATES(20)*1.00000 - STATES(20) ** 3.00000+ 0.250000*STATES(20) ** 2.00000*1.00000 - STATES(20) ** 2.00000+ 0.125000*STATES(20) ** 3.00000*(1.00000 - STATES(20))+ 0.0625000*STATES(20) ** 4.00000)
      ALGBRC(81) =  CONSTS(7)*(ALGBRC(63)+ALGBRC(64))
      ALGBRC(12) =  CONSTS(8)*STATES(2)
      ALGBRC(82) =  CONSTS(11)*(ALGBRC(65)+ALGBRC(64))
      ALGBRC(49) = ( CONSTS(20)*ALGBRC(45)*CONSTS(15))/(CONSTS(24)+CONSTS(15))
      ALGBRC(47) = ( CONSTS(19)*ALGBRC(44)*CONSTS(15))/(CONSTS(23)+CONSTS(15))
      ALGBRC(54) = ( CONSTS(21)*ALGBRC(52)*CONSTS(15))/(CONSTS(25)+CONSTS(15))
      ALGBRC(83) = ( CONSTS(22)*CONSTS(128)*ALGBRC(70))/(CONSTS(26)+ALGBRC(70))
      ALGBRC(30) = CONSTS(37) - STATES(8)
      ALGBRC(84) = ( CONSTS(40)*ALGBRC(69)*ALGBRC(30))/(CONSTS(41)+ALGBRC(30))
      ALGBRC(61) = ( CONSTS(42)*ALGBRC(58)*STATES(8))/(CONSTS(43)+STATES(8))
      ALGBRC(31) = CONSTS(39) - STATES(9)
      ALGBRC(85) = ( CONSTS(44)*ALGBRC(69)*ALGBRC(31))/(CONSTS(45)+ALGBRC(31))
      ALGBRC(33) = ( CONSTS(46)*STATES(9))/(CONSTS(47)+STATES(9))
      ALGBRC(35) = CONSTS(50) - STATES(10)
      ALGBRC(86) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(35))/(CONSTS(54)+ CONSTS(49)*ALGBRC(35))
      ALGBRC(37) = ( CONSTS(49)*CONSTS(57)*CONSTS(52)*STATES(10))/(CONSTS(58)+ CONSTS(49)*STATES(10))
      ALGBRC(39) = CONSTS(50) - STATES(11)
      ALGBRC(87) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(39))/(CONSTS(54)+ CONSTS(49)*ALGBRC(39))
      ALGBRC(41) = ( CONSTS(49)*CONSTS(55)*CONSTS(51)*STATES(11))/(CONSTS(56)+ CONSTS(49)*STATES(11))
      ALGBRC(89) =  CONSTS(9)*ALGBRC(69)*STATES(1)
      ALGBRC(24) =  CONSTS(10)*STATES(3)
      ALGBRC(36) =  (1.00000/( CONSTS(124)*CONSTS(114)))*log(CONSTS(64)/STATES(12))
      ALGBRC(42) =  CONSTS(67)*STATES(15) ** 3.00000*STATES(16)*STATES(17)*(STATES(18) - ALGBRC(36))
      ALGBRC(94) =  EXP( CONSTS(84)*STATES(18)*CONSTS(124))*STATES(12) ** 3.00000*CONSTS(66)
      ALGBRC(95) =  EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124))*CONSTS(64) ** 3.00000*STATES(14)
      ALGBRC(96) =  (CONSTS(80)/( (CONSTS(81) ** 3.00000+CONSTS(64) ** 3.00000)*(CONSTS(82)+CONSTS(66))*(1.00000+ CONSTS(83)*EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124)))))*(ALGBRC(94) - ALGBRC(95))
      ALGBRC(97) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(18)*CONSTS(124))+ 0.0365000*CONSTS(126)*EXP( - STATES(18)*CONSTS(124)))
      ALGBRC(98) = ( (( CONSTS(85)*ALGBRC(97))/(1.00000+CONSTS(86)/STATES(12) ** 1.50000))*CONSTS(65))/(CONSTS(65)+CONSTS(87))
      ALGBRC(101) =  CONSTS(91)*(STATES(18) - ALGBRC(36))
      ALGBRC(102) = ALGBRC(42)+ALGBRC(101)+ 3.00000*ALGBRC(96)+ 3.00000*ALGBRC(98)
      ALGBRC(48) = ( CONSTS(77)*STATES(18)*CONSTS(113)*CONSTS(124)*( STATES(13)*EXP( STATES(18)*CONSTS(124)) - CONSTS(65)))/(EXP( STATES(18)*CONSTS(124)) - 1.00000)
      ALGBRC(34) = STATES(11)/CONSTS(50)
      ALGBRC(50) =  0.500000*(( 0.400000*ALGBRC(34))/CONSTS(122)+0.600000)
      ALGBRC(46) = ( CONSTS(76)*4.00000*STATES(18)*CONSTS(113)*CONSTS(124)*( 0.00100000*EXP( 2.00000*STATES(18)*CONSTS(124)) -  0.341000*CONSTS(66)))/(EXP( 2.00000*STATES(18)*CONSTS(124)) - 1.00000)
      ALGBRC(53) =  ALGBRC(46)*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
      ALGBRC(55) =  (ALGBRC(48)/(1.00000+ALGBRC(53)/CONSTS(79)))*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
      ALGBRC(38) =  (1.00000/( CONSTS(124)*CONSTS(115)))*log(CONSTS(65)/STATES(13))
      ALGBRC(60) =  CONSTS(68)*STATES(24)*( 0.886000*STATES(25)+ 0.114000*STATES(26))*(STATES(18) - ALGBRC(38))
      ALGBRC(62) =  CONSTS(69)*STATES(27)*STATES(28)*(STATES(18) - ALGBRC(38))
      ALGBRC(80) = 1.02000/(1.00000+EXP( 0.238500*((STATES(18) - ALGBRC(38)) - 59.2150)))
      ALGBRC(88) = ( 0.491240*EXP( 0.0803200*((STATES(18)+5.47600) - ALGBRC(38)))+EXP( 0.0617500*((STATES(18) - ALGBRC(38)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(18) - ALGBRC(38))+4.75300)))
      ALGBRC(90) = ALGBRC(80)/(ALGBRC(80)+ALGBRC(88))
      ALGBRC(91) =  CONSTS(70)* (CONSTS(65)/5.40000) ** (1.0 / 2)*ALGBRC(90)*(STATES(18) - ALGBRC(38))
      ALGBRC(92) = 1.00000/(1.00000+EXP((7.48800 - STATES(18))/5.98000))
      ALGBRC(93) =  CONSTS(71)*ALGBRC(92)*(STATES(18) - ALGBRC(38))
      ALGBRC(103) = ((ALGBRC(60)+ALGBRC(62)+ALGBRC(91)+ALGBRC(93)) -  2.00000*ALGBRC(98))+ALGBRC(55)
      ALGBRC(99) = ( CONSTS(88)*STATES(14))/(CONSTS(89)+STATES(14))
      ALGBRC(40) =  (1.00000/( CONSTS(124)*CONSTS(116)))*log(CONSTS(66)/STATES(14))
      ALGBRC(100) =  CONSTS(90)*(STATES(18) - ALGBRC(40))
      ALGBRC(104) = (ALGBRC(53)+ALGBRC(100)+ALGBRC(99)) -  2.00000*ALGBRC(96)
      ALGBRC(32) = ALGBRC(30)/CONSTS(37)
      ALGBRC(110) = ( CONSTS(95)*(1.00000+ 2.00000*ALGBRC(32)))/(1.00000+ 2.00000*CONSTS(120))
      ALGBRC(111) = ( CONSTS(94)*STATES(14) ** 2.00000)/(ALGBRC(110) ** 2.00000+STATES(14) ** 2.00000)
      ALGBRC(112) = ( CONSTS(94)*STATES(31))/CONSTS(96)
      ALGBRC(113) = (STATES(31) - STATES(29))/CONSTS(105)
      ALGBRC(105) = STATES(30)+0.00200000
      ALGBRC(106) = 1.00000 - EXP(- ALGBRC(105)/CONSTS(97))
      ALGBRC(107) = EXP(- ALGBRC(105)/CONSTS(98))
      ALGBRC(108) = CONSTS(99)/(1.00000+EXP((ALGBRC(104)+5.00000)/0.900000))
      ALGBRC(109) =  ALGBRC(108)*ALGBRC(106)*ALGBRC(107)*(STATES(29) - STATES(14))
      ALGBRC(115) = 1.00000/(1.00000+( CONSTS(103)*CONSTS(104))/CONSTS(104)+STATES(29) ** 2.00000)
      ALGBRC(114) = ( CONSTS(106)*CONSTS(109))/CONSTS(109)+STATES(14) ** 2.00000
      ALGBRC(116) = ( CONSTS(107)*CONSTS(110))/CONSTS(110)+STATES(14) ** 2.00000
      ALGBRC(118) = ( CONSTS(108)*CONSTS(111))/CONSTS(111)+STATES(14) ** 2.00000
      ALGBRC(119) = 1.00000/(1.00000+ALGBRC(116)+ALGBRC(114)+ALGBRC(114)+ALGBRC(118))
      ALGBRC(1) = STATES(8)/CONSTS(37)
      ALGBRC(13) = TERNRY(VOI.GT.59.1000.AND.VOI.LT.59.5000, CONSTS(118), CONSTS(117))
      ALGBRC(59) = ALGBRC(53)+ALGBRC(55)
      ALGBRC(117) =  1000.00*(( (STATES(29)+STATES(29)/ALGBRC(115))*CONSTS(61))/CONSTS(59)+( STATES(31)*CONSTS(60))/CONSTS(59))
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(CONSTS(130) - ( CONSTS(128)*CONSTS(129))/CONSTS(29))
      minfunc_0[2] = abs(CONSTS(128) - (CONSTS(16) - CONSTS(130)))
      minfunc_0[3] = abs(CONSTS(129) - (CONSTS(17) - CONSTS(130)))
      RETURN
      END
REAL FUNCTION minfunc_1(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_1[1] = abs(ALGBRC(63) - ( ALGBRC(66)*ALGBRC(67))/CONSTS(4))
      minfunc_1[2] = abs(ALGBRC(64) - ( ALGBRC(63)*ALGBRC(68))/CONSTS(5))
      minfunc_1[3] = abs(ALGBRC(65) - ( ALGBRC(67)*ALGBRC(68))/CONSTS(6))
      minfunc_1[4] = abs(ALGBRC(66) - ((CONSTS(1) - ALGBRC(63)) - ALGBRC(64)))
      minfunc_1[5] = abs(ALGBRC(67) - (((STATES(1) - ALGBRC(63)) - ALGBRC(64)) - ALGBRC(65)))
      minfunc_1[6] = abs(ALGBRC(68) - ((CONSTS(3) - ALGBRC(64)) - ALGBRC(65)))
      minfunc_1[7] = abs(ALGBRC(71) - ( CONSTS(32)*CONSTS(36))/(CONSTS(36)+ALGBRC(69)+ALGBRC(79)))
      minfunc_1[8] = abs(ALGBRC(72) -  (ALGBRC(69)/CONSTS(35))*ALGBRC(69)*(1.00000+ALGBRC(71)/CONSTS(36)))
      minfunc_1[9] = abs(ALGBRC(73) -  ALGBRC(69)*(1.00000+ALGBRC(71)/CONSTS(36)))
      minfunc_1[10] = abs(ALGBRC(74) -  (ALGBRC(79)/CONSTS(35))*ALGBRC(79)*(1.00000+ALGBRC(71)/CONSTS(36)))
      minfunc_1[11] = abs(ALGBRC(75) -  ALGBRC(79)*(1.00000+ALGBRC(71)/CONSTS(36)))
      minfunc_1[12] = abs(ALGBRC(76) -  (CONSTS(33)/ALGBRC(70))*ALGBRC(72))
      minfunc_1[13] = abs(ALGBRC(77) -  (CONSTS(33)/ALGBRC(70))*ALGBRC(74))
      minfunc_1[14] = abs(ALGBRC(78) - (( CONSTS(33)*CONSTS(34))/CONSTS(35)+( CONSTS(33)*ALGBRC(70))/CONSTS(35)+ALGBRC(70) ** 2.00000/CONSTS(35)))
      minfunc_1[15] = abs(ALGBRC(70) - ((STATES(7) - (ALGBRC(76)+ 2.00000*ALGBRC(72)+ 2.00000*ALGBRC(73))) - (ALGBRC(77)+ 2.00000*ALGBRC(74)+ 2.00000*ALGBRC(75))))
      minfunc_1[16] = abs(0.00000 - ( 2.00000*CONSTS(30)*ALGBRC(70) ** 2.00000 -  ALGBRC(69)*(1.00000+ALGBRC(71)/CONSTS(36))*( ALGBRC(78)*ALGBRC(69)+ALGBRC(70) ** 2.00000)))
      minfunc_1[17] = abs(0.00000 - ( 2.00000*CONSTS(31)*ALGBRC(70) ** 2.00000 -  ALGBRC(79)*(1.00000+ALGBRC(71)/CONSTS(36))*( ALGBRC(78)*ALGBRC(79)+ALGBRC(70) ** 2.00000)))
      RETURN
      END
REAL FUNCTION minfunc_2(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_2[1] = abs(ALGBRC(45) - ( ALGBRC(43)*ALGBRC(44))/CONSTS(27))
      minfunc_2[2] = abs(ALGBRC(43) - (STATES(4) - ALGBRC(45)))
      minfunc_2[3] = abs(ALGBRC(44) - (CONSTS(14) - ALGBRC(45)))
      RETURN
      END
REAL FUNCTION minfunc_3(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_3[1] = abs(ALGBRC(52) - ( ALGBRC(51)*ALGBRC(44))/CONSTS(28))
      minfunc_3[2] = abs(ALGBRC(51) - (CONSTS(18) - ALGBRC(52)))
      RETURN
      END
REAL FUNCTION minfunc_4(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_4[1] = abs(ALGBRC(56) - ( ALGBRC(57)*ALGBRC(58))/CONSTS(48))
      minfunc_4[2] = abs(ALGBRC(57) - (STATES(9) - ALGBRC(56)))
      minfunc_4[3] = abs(ALGBRC(58) - (CONSTS(38) - ALGBRC(56)))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END