Generated Code

The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)

The raw code is available.

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