C
C There are a total of 108 entries in the algebraic variable array.
C There are a total of 55 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 R in component model_parameters (joule_per_kilomole_kelvin).
C CONSTS(2) is T in component model_parameters (kelvin).
C CONSTS(3) is F in component model_parameters (coulomb_per_mole).
C CONSTS(4) is Na_e in component model_parameters (millimolar).
C CONSTS(5) is Ca_e in component model_parameters (millimolar).
C CONSTS(6) is K_e in component model_parameters (millimolar).
C CONSTS(7) is ATP_i in component model_parameters (millimolar).
C CONSTS(128) is Vt in component model_parameters (cm3).
C CONSTS(85) is Vd in component model_parameters (cm3).
C CONSTS(87) is Vmyo in component model_parameters (cm3).
C CONSTS(81) is Sms in component model_parameters (cm2).
C CONSTS(91) is Smt in component model_parameters (cm2).
C CONSTS(130) is Cms in component model_parameters (microF).
C CONSTS(129) is Cmt in component model_parameters (microF).
C CONSTS(90) is Rst in component model_parameters (ohm).
C CONSTS(88) is VSRup in component model_parameters (cm3).
C CONSTS(86) is VSRrel in component model_parameters (cm3).
C CONSTS(84) is Vc in component model_parameters (cm3).
C CONSTS(89) is pt in component model_parameters (dimensionless).
C CONSTS(8) is Rot in component model_parameters (ohm_cm).
C CONSTS(9) is rt in component model_parameters (cm).
C CONSTS(10) is Lt in component model_parameters (cm).
C CONSTS(11) is rc in component model_parameters (cm).
C CONSTS(12) is Lc in component model_parameters (cm).
C CONSTS(13) is ptcm in component model_parameters (per_cm2).
C ALGBRC(15) is i_circ in component i_circ (microA).
C STATES(1) is Vm_s in component Vm_s (millivolt).
C STATES(2) is Vm_t in component Vm_t (millivolt).
C CONSTS(92) is g_Na_s in component membrane_permeabilities (milliS).
C CONSTS(93) is g_Naps_s in component membrane_permeabilities (milliS).
C CONSTS(94) is q_Kr_s in component membrane_permeabilities (milliS).
C CONSTS(95) is q_Ks_s in component membrane_permeabilities (milliS).
C CONSTS(96) is g_Kp_s in component membrane_permeabilities (milliS).
C CONSTS(97) is g_Kto_s in component membrane_permeabilities (milliS).
C CONSTS(98) is g_K1_s in component membrane_permeabilities (milliS).
C CONSTS(99) is g_KNa_s in component membrane_permeabilities (milliS).
C CONSTS(100) is g_KATP_s in component membrane_permeabilities (milliS).
C CONSTS(101) is g_Nab_s in component membrane_permeabilities (milliS).
C CONSTS(102) is g_Cab_s in component membrane_permeabilities (milliS).
C CONSTS(103) is i_pCa_max_s in component membrane_permeabilities (microA).
C CONSTS(104) is i_NaK_max_s in component membrane_permeabilities (microA).
C CONSTS(105) is i_NaCa_max_s in component membrane_permeabilities (microA).
C CONSTS(106) is P_CaL_s in component membrane_permeabilities (litre_per_second).
C CONSTS(107) is P_KL_s in component membrane_permeabilities (litre_per_second).
C CONSTS(108) is P_nsNa_s in component membrane_permeabilities (litre_per_second).
C CONSTS(109) is P_nsK_s in component membrane_permeabilities (litre_per_second).
C CONSTS(110) is g_Na_t in component membrane_permeabilities (milliS).
C CONSTS(111) is g_Naps_t in component membrane_permeabilities (milliS).
C CONSTS(112) is q_Kr_t in component membrane_permeabilities (milliS).
C CONSTS(113) is q_Ks_t in component membrane_permeabilities (milliS).
C CONSTS(114) is g_Kp_t in component membrane_permeabilities (milliS).
C CONSTS(115) is g_Kto_t in component membrane_permeabilities (milliS).
C CONSTS(116) is g_K1_t in component membrane_permeabilities (milliS).
C CONSTS(117) is g_KNa_t in component membrane_permeabilities (milliS).
C CONSTS(118) is g_KATP_t in component membrane_permeabilities (milliS).
C CONSTS(119) is g_Nab_t in component membrane_permeabilities (milliS).
C CONSTS(120) is g_Cab_t in component membrane_permeabilities (milliS).
C CONSTS(121) is i_pCa_max_t in component membrane_permeabilities (microA).
C CONSTS(122) is i_NaK_max_t in component membrane_permeabilities (microA).
C CONSTS(123) is i_NaCa_max_t in component membrane_permeabilities (microA).
C CONSTS(124) is P_CaL_t in component membrane_permeabilities (litre_per_second).
C CONSTS(125) is P_KL_t in component membrane_permeabilities (litre_per_second).
C CONSTS(126) is P_nsNa_t in component membrane_permeabilities (litre_per_second).
C CONSTS(127) is P_nsK_t in component membrane_permeabilities (litre_per_second).
C CONSTS(14) is g_Na in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(15) is fNat in component membrane_permeabilities (dimensionless).
C CONSTS(16) is g_Naps in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(17) is fNapst in component membrane_permeabilities (dimensionless).
C CONSTS(18) is q_Kr in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(19) is fKrt in component membrane_permeabilities (dimensionless).
C CONSTS(20) is q_Ks in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(21) is fKst in component membrane_permeabilities (dimensionless).
C CONSTS(22) is g_Kp in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(23) is fKpt in component membrane_permeabilities (dimensionless).
C CONSTS(24) is g_Kto in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(25) is fKtot in component membrane_permeabilities (dimensionless).
C CONSTS(26) is g_K1 in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(27) is fK1t in component membrane_permeabilities (dimensionless).
C CONSTS(28) is g_KNa in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(29) is fKNat in component membrane_permeabilities (dimensionless).
C CONSTS(30) is g_KATP in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(31) is fKATPt in component membrane_permeabilities (dimensionless).
C CONSTS(32) is g_Nab in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(33) is fNabt in component membrane_permeabilities (dimensionless).
C CONSTS(34) is g_Cab in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(35) is fCabt in component membrane_permeabilities (dimensionless).
C CONSTS(36) is P_CaL in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(37) is fCaLt in component membrane_permeabilities (dimensionless).
C CONSTS(38) is P_KL in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(39) is P_nsNa in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(40) is fnsNat in component membrane_permeabilities (dimensionless).
C CONSTS(41) is P_nsK in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(42) is fnsKt in component membrane_permeabilities (dimensionless).
C CONSTS(43) is i_NaCa_max in component membrane_permeabilities (microA_per_cm2).
C CONSTS(44) is fNaCat in component membrane_permeabilities (dimensionless).
C CONSTS(45) is i_NaK_max in component membrane_permeabilities (microA_per_cm2).
C CONSTS(46) is fNaKt in component membrane_permeabilities (dimensionless).
C CONSTS(47) is i_pCa_max in component membrane_permeabilities (microA_per_cm2).
C CONSTS(48) is fpCat in component membrane_permeabilities (dimensionless).
C ALGBRC(30) is i_Kext in component i_Kext (microA).
C CONSTS(49) is stim_Period in component i_Kext (second).
C ALGBRC(36) is i_Na_s in component i_Na_s (microA).
C ALGBRC(33) is E_Na_s in component i_Na_s (millivolt).
C STATES(3) is Na_i in component ion_concentrations (millimolar).
C STATES(4) is K_i in component ion_concentrations (millimolar).
C STATES(5) is m in component i_Na_s_m_gate (dimensionless).
C STATES(6) is h in component i_Na_s_h_gate (dimensionless).
C ALGBRC(1) is m_infinity in component i_Na_s_m_gate (dimensionless).
C ALGBRC(16) is tau_m in component i_Na_s_m_gate (second).
C ALGBRC(2) is h_infinity in component i_Na_s_h_gate (dimensionless).
C ALGBRC(17) is tau_h in component i_Na_s_h_gate (second).
C ALGBRC(40) is i_Na_t in component i_Na_t (microA).
C ALGBRC(39) is E_Na_t in component i_Na_t (millivolt).
C STATES(7) is Na_t in component ion_concentrations (millimolar).
C STATES(8) is K_t in component ion_concentrations (millimolar).
C STATES(9) is m in component i_Na_t_m_gate (dimensionless).
C STATES(10) is h in component i_Na_t_h_gate (dimensionless).
C ALGBRC(3) is m_infinity in component i_Na_t_m_gate (dimensionless).
C ALGBRC(18) is tau_m in component i_Na_t_m_gate (second).
C ALGBRC(4) is h_infinity in component i_Na_t_h_gate (dimensionless).
C ALGBRC(19) is tau_h in component i_Na_t_h_gate (second).
C ALGBRC(41) is i_Naps_s in component i_Naps_s (microA).
C ALGBRC(42) is i_Naps_t in component i_Naps_t (microA).
C ALGBRC(43) is i_CaL_s in component i_CaL_s (microA).
C ALGBRC(45) is i_KL_s in component i_CaL_s (microA).
C ALGBRC(5) is alfas in component i_CaL_s (per_second).
C ALGBRC(20) is betas in component i_CaL_s (per_second).
C ALGBRC(31) is gama in component i_CaL_s (per_second).
C CONSTS(50) is omega in component i_CaL_s (per_second).
C CONSTS(51) is a in component i_CaL_s (dimensionless).
C CONSTS(52) is b in component i_CaL_s (dimensionless).
C CONSTS(53) is f in component i_CaL_s (per_second).
C CONSTS(54) is g in component i_CaL_s (per_second).
C CONSTS(55) is f2 in component i_CaL_s (per_second).
C CONSTS(56) is g2 in component i_CaL_s (per_second).
C ALGBRC(34) is alfa2s in component i_CaL_s (per_second).
C ALGBRC(37) is beta2s in component i_CaL_s (per_second).
C STATES(11) is Cst in component i_CaL_s (dimensionless).
C STATES(12) is C1 in component i_CaL_s (dimensionless).
C STATES(13) is C2 in component i_CaL_s (dimensionless).
C STATES(14) is C3 in component i_CaL_s (dimensionless).
C STATES(15) is C4 in component i_CaL_s (dimensionless).
C STATES(16) is Co in component i_CaL_s (dimensionless).
C STATES(17) is Ccast in component i_CaL_s (dimensionless).
C STATES(18) is Cca1 in component i_CaL_s (dimensionless).
C STATES(19) is Cca2 in component i_CaL_s (dimensionless).
C STATES(20) is Cca3 in component i_CaL_s (dimensionless).
C STATES(21) is Cca4 in component i_CaL_s (dimensionless).
C STATES(22) is Ccao in component i_CaL_s (dimensionless).
C STATES(23) is Ca_ss in component ion_concentrations (millimolar).
C ALGBRC(44) is i_CaL_t in component i_CaL_t (microA).
C STATES(24) is y in component i_CaL_s_y_gate (dimensionless).
C ALGBRC(6) is y_infinity in component i_CaL_s_y_gate (dimensionless).
C ALGBRC(21) is tau_y in component i_CaL_s_y_gate (second).
C ALGBRC(46) is i_KL_t in component i_CaL_t (microA).
C ALGBRC(22) is alfat in component i_CaL_t (per_second).
C ALGBRC(32) is betat in component i_CaL_t (per_second).
C ALGBRC(7) is gama in component i_CaL_t (per_second).
C CONSTS(57) is omega in component i_CaL_t (per_second).
C CONSTS(58) is a in component i_CaL_t (dimensionless).
C CONSTS(59) is b in component i_CaL_t (dimensionless).
C CONSTS(60) is f in component i_CaL_t (per_second).
C CONSTS(61) is g in component i_CaL_t (per_second).
C CONSTS(62) is f2 in component i_CaL_t (per_second).
C CONSTS(63) is g2 in component i_CaL_t (per_second).
C ALGBRC(35) is alfa2t in component i_CaL_t (per_second).
C ALGBRC(38) is beta2t in component i_CaL_t (per_second).
C STATES(25) is TCst in component i_CaL_t (dimensionless).
C STATES(26) is TC1 in component i_CaL_t (dimensionless).
C STATES(27) is TC2 in component i_CaL_t (dimensionless).
C STATES(28) is TC3 in component i_CaL_t (dimensionless).
C STATES(29) is TC4 in component i_CaL_t (dimensionless).
C STATES(30) is TCo in component i_CaL_t (dimensionless).
C STATES(31) is TCcast in component i_CaL_t (dimensionless).
C STATES(32) is TCca1 in component i_CaL_t (dimensionless).
C STATES(33) is TCca2 in component i_CaL_t (dimensionless).
C STATES(34) is TCca3 in component i_CaL_t (dimensionless).
C STATES(35) is TCca4 in component i_CaL_t (dimensionless).
C STATES(36) is TCcao in component i_CaL_t (dimensionless).
C STATES(37) is Ca_t in component ion_concentrations (millimolar).
C STATES(38) is y in component i_CaL_t_y_gate (dimensionless).
C ALGBRC(8) is y_infinity in component i_CaL_t_y_gate (dimensionless).
C ALGBRC(23) is tau_y in component i_CaL_t_y_gate (second).
C ALGBRC(49) is i_Kr_s in component i_Kr_s (microA).
C ALGBRC(47) is E_Kr_s in component i_Kr_s (millivolt).
C STATES(39) is xr in component i_Kr_s_xr_gate (dimensionless).
C ALGBRC(48) is xri in component i_Kr_s_xri_gate (dimensionless).
C ALGBRC(9) is xr_infinity in component i_Kr_s_xr_gate (dimensionless).
C ALGBRC(24) is tau_xr in component i_Kr_s_xr_gate (second).
C ALGBRC(52) is i_Kr_t in component i_Kr_t (microA).
C ALGBRC(50) is E_Kr_t in component i_Kr_t (millivolt).
C STATES(40) is xr in component i_Kr_t_xr_gate (dimensionless).
C ALGBRC(51) is xri in component i_Kr_t_xri_gate (dimensionless).
C ALGBRC(10) is xr_infinity in component i_Kr_t_xr_gate (dimensionless).
C ALGBRC(25) is tau_xr in component i_Kr_t_xr_gate (second).
C ALGBRC(54) is i_Ks_s in component i_Ks_s (microA).
C ALGBRC(53) is E_Ks_s in component i_Ks_s (millivolt).
C CONSTS(64) is PRNaK in component i_Ks_s (dimensionless).
C STATES(41) is Ca_i in component ion_concentrations (millimolar).
C STATES(42) is xs in component i_Ks_s_xs_gate (dimensionless).
C ALGBRC(11) is xs_infinity in component i_Ks_s_xs_gate (dimensionless).
C ALGBRC(26) is tau_xs in component i_Ks_s_xs_gate (second).
C ALGBRC(56) is i_Ks_t in component i_Ks_t (microA).
C ALGBRC(55) is E_Ks_t in component i_Ks_t (millivolt).
C CONSTS(65) is PRNaK in component i_Ks_t (dimensionless).
C STATES(43) is xs in component i_Ks_t_xs_gate (dimensionless).
C ALGBRC(12) is xs_infinity in component i_Ks_t_xs_gate (dimensionless).
C ALGBRC(27) is tau_xs in component i_Ks_t_xs_gate (second).
C ALGBRC(60) is i_K1_s in component i_K1_s (microA).
C ALGBRC(57) is E_K1_s in component i_K1_s (millivolt).
C ALGBRC(58) is aK1s in component i_K1_s (dimensionless).
C ALGBRC(59) is bK1s in component i_K1_s (dimensionless).
C ALGBRC(64) is i_K1_t in component i_K1_t (microA).
C ALGBRC(61) is E_K1_t in component i_K1_t (millivolt).
C ALGBRC(62) is aK1t in component i_K1_t (dimensionless).
C ALGBRC(63) is bK1t in component i_K1_t (dimensionless).
C ALGBRC(67) is i_Kp_s in component i_Kp_s (microA).
C ALGBRC(66) is kps in component i_Kp_s (dimensionless).
C ALGBRC(65) is E_Kp_s in component i_Kp_s (millivolt).
C ALGBRC(70) is i_Kp_t in component i_Kp_t (microA).
C ALGBRC(69) is kpt in component i_Kp_t (dimensionless).
C ALGBRC(68) is E_Kp_t in component i_Kp_t (millivolt).
C ALGBRC(72) is i_Kto_s in component i_Kto_s (microA).
C CONSTS(66) is fr2 in component i_Kto_s (dimensionless).
C ALGBRC(71) is E_Kto_s in component i_Kto_s (millivolt).
C STATES(44) is r2 in component i_Kto_s_r2_gate (dimensionless).
C STATES(45) is r3 in component i_Kto_s_r3_gate (dimensionless).
C CONSTS(67) is tau_r2 in component i_Kto_s_r2_gate (second).
C ALGBRC(13) is r2_infinity in component i_Kto_s_r2_gate (dimensionless).
C CONSTS(68) is tau_r3 in component i_Kto_s_r3_gate (second).
C ALGBRC(28) is r3_infinity in component i_Kto_s_r3_gate (dimensionless).
C ALGBRC(74) is i_Kto_t in component i_Kto_t (microA).
C CONSTS(69) is fr2 in component i_Kto_t (dimensionless).
C ALGBRC(73) is E_Kto_t in component i_Kto_t (millivolt).
C STATES(46) is r2 in component i_Kto_t_r2_gate (dimensionless).
C STATES(47) is r3 in component i_Kto_t_r3_gate (dimensionless).
C CONSTS(70) is tau_r2 in component i_Kto_t_r2_gate (second).
C ALGBRC(14) is r2_infinity in component i_Kto_t_r2_gate (dimensionless).
C ALGBRC(29) is r3_infinity in component i_Kto_t_r3_gate (dimensionless).
C CONSTS(71) is tau_r3 in component i_Kto_t_r3_gate (second).
C ALGBRC(76) is i_KNa_s in component i_KNa_s (microA).
C ALGBRC(75) is E_KNa_s in component i_KNa_s (millivolt).
C ALGBRC(78) is i_KNa_t in component i_KNa_t (microA).
C ALGBRC(77) is E_KNa_t in component i_KNa_t (millivolt).
C ALGBRC(79) is i_nsNa_s in component i_nsNa_s (microA).
C ALGBRC(80) is i_nsNa_t in component i_nsNa_t (microA).
C ALGBRC(81) is i_nsK_s in component i_nsK_s (microA).
C ALGBRC(82) is i_nsK_t in component i_nsK_t (microA).
C ALGBRC(83) is i_Nab_s in component i_Nab_s (microA).
C ALGBRC(84) is i_Nab_t in component i_Nab_t (microA).
C ALGBRC(86) is i_Cab_s in component i_Cab_s (microA).
C ALGBRC(85) is E_Ca_s in component i_Cab_s (millivolt).
C ALGBRC(88) is i_Cab_t in component i_Cab_t (microA).
C ALGBRC(87) is E_Ca_t in component i_Cab_t (millivolt).
C ALGBRC(89) is i_NaCa_s in component i_NaCa_s (microA).
C ALGBRC(90) is i_NaCa_t in component i_NaCa_t (microA).
C ALGBRC(91) is i_NaK_s in component i_NaK_s (microA).
C ALGBRC(92) is i_NaK_t in component i_NaK_t (microA).
C ALGBRC(94) is i_pCa_s in component i_pCa_s (microA).
C ALGBRC(95) is i_pCa_t in component i_pCa_t (microA).
C ALGBRC(99) is i_KATP_s in component i_KATP_s (microA).
C ALGBRC(96) is E_KATP_s in component i_KATP_s (millivolt).
C ALGBRC(103) is i_KATP_t in component i_KATP_t (microA).
C ALGBRC(101) is E_KATP_t in component i_KATP_t (millivolt).
C ALGBRC(93) is JteNa in component t_tubular_ion_fluxes (millimolar_per_second).
C ALGBRC(97) is JteCa in component t_tubular_ion_fluxes (millimolar_per_second).
C ALGBRC(105) is JteK in component t_tubular_ion_fluxes (millimolar_per_second).
C CONSTS(72) is tau_Na in component t_tubular_ion_fluxes (second).
C CONSTS(73) is tau_Ca in component t_tubular_ion_fluxes (second).
C CONSTS(74) is tau_K in component t_tubular_ion_fluxes (second).
C ALGBRC(98) is JCaSRup in component JCaSRup (millimolar_per_second).
C ALGBRC(100) is JCaSRleak in component JCaSRleak (millimolar_per_second).
C STATES(48) is CaSRup in component CaSRup (millimolar).
C ALGBRC(102) is Jtr in component Jtr (millimolar_per_second).
C CONSTS(75) is tau_tr in component Jtr (second).
C STATES(49) is CaSRrel in component CaSRrel (millimolar).
C ALGBRC(104) is JCaSRrel in component JCaSRrel (millimolar_per_second).
C CONSTS(82) is kap in component JCaSRrel (per_millimolar4_per_second).
C CONSTS(76) is kam in component JCaSRrel (per_second).
C CONSTS(83) is kbp in component JCaSRrel (per_millimolar3_per_second).
C CONSTS(77) is kbm in component JCaSRrel (per_second).
C CONSTS(78) is kcp in component JCaSRrel (per_second).
C CONSTS(79) is kcm in component JCaSRrel (per_second).
C STATES(50) is F1 in component JCaSRrel (dimensionless).
C STATES(51) is F2 in component JCaSRrel (dimensionless).
C STATES(52) is F3 in component JCaSRrel (dimensionless).
C STATES(53) is F4 in component JCaSRrel (dimensionless).
C ALGBRC(106) is JCad in component JCad (millimolar_per_second).
C CONSTS(80) is tau_d in component JCad (second).
C STATES(54) is BTRH in component ion_concentrations (millimolar).
C STATES(55) is BTRL in component ion_concentrations (millimolar).
C ALGBRC(107) is dBTRH in component ion_concentrations (millimolar_per_second).
C ALGBRC(108) is dBTRL in component ion_concentrations (millimolar_per_second).
C RATES(5) is d/dt m in component i_Na_s_m_gate (dimensionless).
C RATES(6) is d/dt h in component i_Na_s_h_gate (dimensionless).
C RATES(9) is d/dt m in component i_Na_t_m_gate (dimensionless).
C RATES(10) is d/dt h in component i_Na_t_h_gate (dimensionless).
C RATES(11) is d/dt Cst in component i_CaL_s (dimensionless).
C RATES(12) is d/dt C1 in component i_CaL_s (dimensionless).
C RATES(13) is d/dt C2 in component i_CaL_s (dimensionless).
C RATES(14) is d/dt C3 in component i_CaL_s (dimensionless).
C RATES(15) is d/dt C4 in component i_CaL_s (dimensionless).
C RATES(16) is d/dt Co in component i_CaL_s (dimensionless).
C RATES(17) is d/dt Ccast in component i_CaL_s (dimensionless).
C RATES(18) is d/dt Cca1 in component i_CaL_s (dimensionless).
C RATES(19) is d/dt Cca2 in component i_CaL_s (dimensionless).
C RATES(20) is d/dt Cca3 in component i_CaL_s (dimensionless).
C RATES(21) is d/dt Cca4 in component i_CaL_s (dimensionless).
C RATES(22) is d/dt Ccao in component i_CaL_s (dimensionless).
C RATES(24) is d/dt y in component i_CaL_s_y_gate (dimensionless).
C RATES(25) is d/dt TCst in component i_CaL_t (dimensionless).
C RATES(26) is d/dt TC1 in component i_CaL_t (dimensionless).
C RATES(27) is d/dt TC2 in component i_CaL_t (dimensionless).
C RATES(28) is d/dt TC3 in component i_CaL_t (dimensionless).
C RATES(29) is d/dt TC4 in component i_CaL_t (dimensionless).
C RATES(30) is d/dt TCo in component i_CaL_t (dimensionless).
C RATES(31) is d/dt TCcast in component i_CaL_t (dimensionless).
C RATES(32) is d/dt TCca1 in component i_CaL_t (dimensionless).
C RATES(33) is d/dt TCca2 in component i_CaL_t (dimensionless).
C RATES(34) is d/dt TCca3 in component i_CaL_t (dimensionless).
C RATES(35) is d/dt TCca4 in component i_CaL_t (dimensionless).
C RATES(36) is d/dt TCcao in component i_CaL_t (dimensionless).
C RATES(38) is d/dt y in component i_CaL_t_y_gate (dimensionless).
C RATES(39) is d/dt xr in component i_Kr_s_xr_gate (dimensionless).
C RATES(40) is d/dt xr in component i_Kr_t_xr_gate (dimensionless).
C RATES(42) is d/dt xs in component i_Ks_s_xs_gate (dimensionless).
C RATES(43) is d/dt xs in component i_Ks_t_xs_gate (dimensionless).
C RATES(44) is d/dt r2 in component i_Kto_s_r2_gate (dimensionless).
C RATES(45) is d/dt r3 in component i_Kto_s_r3_gate (dimensionless).
C RATES(46) is d/dt r2 in component i_Kto_t_r2_gate (dimensionless).
C RATES(47) is d/dt r3 in component i_Kto_t_r3_gate (dimensionless).
C RATES(50) is d/dt F1 in component JCaSRrel (dimensionless).
C RATES(51) is d/dt F2 in component JCaSRrel (dimensionless).
C RATES(52) is d/dt F3 in component JCaSRrel (dimensionless).
C RATES(53) is d/dt F4 in component JCaSRrel (dimensionless).
C RATES(49) is d/dt CaSRrel in component CaSRrel (millimolar).
C RATES(48) is d/dt CaSRup in component CaSRup (millimolar).
C RATES(23) is d/dt Ca_ss in component ion_concentrations (millimolar).
C RATES(41) is d/dt Ca_i in component ion_concentrations (millimolar).
C RATES(54) is d/dt BTRH in component ion_concentrations (millimolar).
C RATES(55) is d/dt BTRL in component ion_concentrations (millimolar).
C RATES(3) is d/dt Na_i in component ion_concentrations (millimolar).
C RATES(4) is d/dt K_i in component ion_concentrations (millimolar).
C RATES(7) is d/dt Na_t in component ion_concentrations (millimolar).
C RATES(37) is d/dt Ca_t in component ion_concentrations (millimolar).
C RATES(8) is d/dt K_t in component ion_concentrations (millimolar).
C RATES(1) is d/dt Vm_s in component Vm_s (millivolt).
C RATES(2) is d/dt Vm_t in component Vm_t (millivolt).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 8310
      CONSTS(2) = 310
      CONSTS(3) = 96500
      CONSTS(4) = 140
      CONSTS(5) = 1.8
      CONSTS(6) = 5.4
      CONSTS(7) = 6.8
      CONSTS(8) = 83.33
      CONSTS(9) = 148.15e-7
      CONSTS(10) = 5.93e-4
      CONSTS(11) = 5.93e-4
      CONSTS(12) = 0.013
      CONSTS(13) = 21028875
      STATES(1) = -8.5207812e1
      STATES(2) = -8.5208034e1
      CONSTS(14) = 30
      CONSTS(15) = 0.64
      CONSTS(16) = 0.0053
      CONSTS(17) = 0.526
      CONSTS(18) = 0.8
      CONSTS(19) = 0.526
      CONSTS(20) = 0.1975
      CONSTS(21) = 0.526
      CONSTS(22) = 0.006
      CONSTS(23) = 0.526
      CONSTS(24) = 0
      CONSTS(25) = 0.526
      CONSTS(26) = 0.75
      CONSTS(27) = 0.8
      CONSTS(28) = 0.12848
      CONSTS(29) = 0.526
      CONSTS(30) = 1
      CONSTS(31) = 0.526
      CONSTS(32) = 0.00075
      CONSTS(33) = 0.526
      CONSTS(34) = 0.0021
      CONSTS(35) = 0.526
      CONSTS(36) = 48e-4
      CONSTS(37) = 0.64
      CONSTS(38) = 48e-7
      CONSTS(39) = 1.75e-7
      CONSTS(40) = 0.526
      CONSTS(41) = 1.75e-7
      CONSTS(42) = 0.526
      CONSTS(43) = 0.00025
      CONSTS(44) = 0.526
      CONSTS(45) = 1.5
      CONSTS(46) = 0.526
      CONSTS(47) = 1.15
      CONSTS(48) = 0.2
      CONSTS(49) = 1
      STATES(3) = 1.1120279e1
      STATES(4) = 1.3678926e2
      STATES(5) = 3.4130493e-3
      STATES(6) = 8.2699973e-1
      STATES(7) = 1.3993239e2
      STATES(8) = 5.4140321e0
      STATES(9) = 3.4129472e-3
      STATES(10) = 8.2700551e-1
      CONSTS(50) = 10
      CONSTS(51) = 2
      CONSTS(52) = 2
      CONSTS(53) = 300
      CONSTS(54) = 2000
      CONSTS(55) = 5
      CONSTS(56) = 7000
      STATES(11) = 9.9797984e-1
      STATES(12) = 7.5722514e-5
      STATES(13) = 2.1545646e-9
      STATES(14) = 2.7463178e-14
      STATES(15) = -5.7569102e-17
      STATES(16) = -1.3496934e-17
      STATES(17) = 1.9155065e-3
      STATES(18) = 5.8135383e-7
      STATES(19) = 6.615513e-11
      STATES(20) = 3.3426555e-15
      STATES(21) = 4.6640361e-19
      STATES(22) = 3.7986638e-22
      STATES(23) = 9.5977033e-5
      STATES(24) = 7.8134852e-1
      CONSTS(57) = 10
      CONSTS(58) = 2
      CONSTS(59) = 2
      CONSTS(60) = 300
      CONSTS(61) = 2000
      CONSTS(62) = 5
      CONSTS(63) = 7000
      STATES(25) = 9.9798934e-1
      STATES(26) = 7.5720257e-5
      STATES(27) = 2.1544156e-9
      STATES(28) = 2.7459475e-14
      STATES(29) = -5.7367189e-17
      STATES(30) = -1.3449454e-17
      STATES(31) = 1.9155245e-3
      STATES(32) = 5.8133643e-7
      STATES(33) = 6.6150547e-11
      STATES(34) = 3.3423167e-15
      STATES(35) = 4.6165955e-19
      STATES(36) = 3.75926e-22
      STATES(37) = 1.8476402e0
      STATES(38) = 7.813477e-1
      STATES(39) = 2.0469344e-4
      STATES(40) = 2.0469344e-4
      CONSTS(64) = 0.01833
      STATES(41) = 8.8787034e-5
      STATES(42) = 6.1359896e-3
      CONSTS(65) = 0.01833
      STATES(43) = 6.1359896e-3
      CONSTS(66) = 1
      STATES(44) = 9.625701e-1
      STATES(45) = 9.6250049e-1
      CONSTS(67) = 0.02
      CONSTS(68) = 0.08
      CONSTS(69) = 1
      STATES(46) = 9.625701e-1
      STATES(47) = 9.6250049e-1
      CONSTS(70) = 0.02
      CONSTS(71) = 0.08
      CONSTS(72) = 0.2
      CONSTS(73) = 0.24
      CONSTS(74) = 0.2
      STATES(48) = 1.0200296e0
      CONSTS(75) = 0.03448
      STATES(49) = 1.0326252e0
      CONSTS(76) = 576
      CONSTS(77) = 1930
      CONSTS(78) = 18
      CONSTS(79) = 0.8
      STATES(50) = 8.053082e-1
      STATES(51) = 2.6639195e-4
      STATES(52) = 4.9455459e-10
      STATES(53) = 1.9442578e-1
      CONSTS(80) = 0.003125
      STATES(54) = 9.6700747e-1
      STATES(55) = 8.1740868e-2
      CONSTS(81) =  2.00000*3.14159*CONSTS(11) ** 2.00000+ 2.00000*3.14159*CONSTS(11)*CONSTS(12)
      CONSTS(82) =  3.00000*1.21500e+13
      CONSTS(83) = 4.05000e+09
      CONSTS(84) =  3.14159*CONSTS(11) ** 2.00000*CONSTS(12)
      CONSTS(85) =  CONSTS(84)*0.680000*6.00000e-05
      CONSTS(86) =  CONSTS(84)*0.00316000
      CONSTS(87) =  CONSTS(84)*0.680000
      CONSTS(88) =  CONSTS(84)*0.0550000
      CONSTS(89) =  CONSTS(13)*2.00000*3.14159*CONSTS(11)*CONSTS(12)
      CONSTS(90) = (( CONSTS(8)*CONSTS(10))/2.00000)/( 3.14159*CONSTS(9) ** 2.00000*CONSTS(89))
      CONSTS(91) =  2.00000*3.14159*CONSTS(9)*CONSTS(10)*CONSTS(89)
      CONSTS(92) =  CONSTS(14)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(15))
      CONSTS(93) =  CONSTS(16)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(17))
      CONSTS(94) =  CONSTS(18)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(19))
      CONSTS(95) =  CONSTS(20)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(21))
      CONSTS(96) =  CONSTS(22)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(23))
      CONSTS(97) =  CONSTS(24)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(25))
      CONSTS(98) =  CONSTS(26)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(27))
      CONSTS(99) =  CONSTS(28)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(29))
      CONSTS(100) =  CONSTS(30)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(31))
      CONSTS(101) =  CONSTS(32)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(33))
      CONSTS(102) =  CONSTS(34)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(35))
      CONSTS(103) =  CONSTS(47)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(48))
      CONSTS(104) =  CONSTS(45)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(46))
      CONSTS(105) =  CONSTS(43)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(44))
      CONSTS(106) =  CONSTS(36)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(37))
      CONSTS(107) =  CONSTS(38)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(37))
      CONSTS(108) =  CONSTS(39)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(40))
      CONSTS(109) =  CONSTS(41)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(42))
      CONSTS(110) =  CONSTS(14)*(CONSTS(81)+CONSTS(91))*CONSTS(15)
      CONSTS(111) =  CONSTS(16)*(CONSTS(81)+CONSTS(91))*CONSTS(17)
      CONSTS(112) =  CONSTS(18)*(CONSTS(81)+CONSTS(91))*CONSTS(19)
      CONSTS(113) =  CONSTS(20)*(CONSTS(81)+CONSTS(91))*CONSTS(21)
      CONSTS(114) =  CONSTS(22)*(CONSTS(81)+CONSTS(91))*CONSTS(23)
      CONSTS(115) =  CONSTS(24)*(CONSTS(81)+CONSTS(91))*CONSTS(25)
      CONSTS(116) =  CONSTS(26)*(CONSTS(81)+CONSTS(91))*CONSTS(27)
      CONSTS(117) =  CONSTS(28)*(CONSTS(81)+CONSTS(91))*CONSTS(29)
      CONSTS(118) =  CONSTS(30)*(CONSTS(81)+CONSTS(91))*CONSTS(31)
      CONSTS(119) =  CONSTS(32)*(CONSTS(81)+CONSTS(91))*CONSTS(33)
      CONSTS(120) =  CONSTS(34)*(CONSTS(81)+CONSTS(91))*CONSTS(35)
      CONSTS(121) =  CONSTS(47)*(CONSTS(81)+CONSTS(91))*CONSTS(48)
      CONSTS(122) =  CONSTS(45)*(CONSTS(81)+CONSTS(91))*CONSTS(46)
      CONSTS(123) =  CONSTS(43)*(CONSTS(81)+CONSTS(91))*CONSTS(44)
      CONSTS(124) =  CONSTS(36)*(CONSTS(81)+CONSTS(91))*CONSTS(37)
      CONSTS(125) =  CONSTS(38)*(CONSTS(81)+CONSTS(91))*CONSTS(37)
      CONSTS(126) =  CONSTS(39)*(CONSTS(81)+CONSTS(91))*CONSTS(40)
      CONSTS(127) =  CONSTS(41)*(CONSTS(81)+CONSTS(91))*CONSTS(42)
      CONSTS(128) =  3.14159*CONSTS(9) ** 2.00000*CONSTS(10)*CONSTS(89)
      CONSTS(129) =  CONSTS(91)*1.00000
      CONSTS(130) =  CONSTS(81)*1.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(16) =  CONSTS(53)*STATES(15) -  CONSTS(54)*STATES(16)
      RATES(22) =  CONSTS(55)*STATES(21) -  CONSTS(56)*STATES(22)
      RATES(30) =  CONSTS(60)*STATES(29) -  CONSTS(61)*STATES(30)
      RATES(36) =  CONSTS(62)*STATES(35) -  CONSTS(63)*STATES(36)
      RATES(50) =  STATES(51)*CONSTS(76) -  STATES(50)*CONSTS(82)*STATES(23) ** 4.00000
      RATES(51) = ( STATES(50)*CONSTS(82)*STATES(23) ** 4.00000+ STATES(52)*CONSTS(77)+ STATES(53)*CONSTS(79)) -  STATES(51)*(CONSTS(76)+ CONSTS(83)*STATES(23) ** 3.00000+CONSTS(78))
      RATES(52) =  STATES(51)*CONSTS(83)*STATES(23) ** 3.00000 -  STATES(52)*CONSTS(77)
      RATES(53) =  STATES(51)*CONSTS(78) -  STATES(53)*CONSTS(79)
      ALGBRC(13) = 1.00600/(1.00000+EXP((STATES(1)+36.6900+10.0000)/12.4300))
      RATES(44) = (ALGBRC(13) - STATES(44))/CONSTS(67)
      ALGBRC(14) = 1.00600/(1.00000+EXP((STATES(2)+36.6900+10.0000)/12.4300))
      RATES(46) = (ALGBRC(14) - STATES(46))/CONSTS(70)
      ALGBRC(1) = 1.00000/(1.00000+EXP(((STATES(1)+52.2000) - 9.00000)/- 7.40000))
      ALGBRC(16) = (0.00100000/( 101.600*EXP( STATES(1)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(1)))+0.000100000)/4.83000
      RATES(5) = (ALGBRC(1) - STATES(5))/ALGBRC(16)
      ALGBRC(2) = 1.00000/(1.00000+EXP(((STATES(1)+85.6000) - 9.00000)/5.50000))
      ALGBRC(17) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(1))+ 6.53700*EXP( STATES(1)*0.0801600))+0.000500000)/4.83000
      RATES(6) = (ALGBRC(2) - STATES(6))/ALGBRC(17)
      ALGBRC(3) = 1.00000/(1.00000+EXP(((STATES(2)+52.2000) - 9.00000)/- 7.40000))
      ALGBRC(18) = (0.00100000/( 101.600*EXP( STATES(2)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(2)))+0.000100000)/4.83000
      RATES(9) = (ALGBRC(3) - STATES(9))/ALGBRC(18)
      ALGBRC(4) = 1.00000/(1.00000+EXP(((STATES(2)+85.6000) - 9.00000)/5.50000))
      ALGBRC(19) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(2))+ 6.53700*EXP( STATES(2)*0.0801600))+0.000500000)/4.83000
      RATES(10) = (ALGBRC(4) - STATES(10))/ALGBRC(19)
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/6.00000))
      ALGBRC(21) = 0.00100000/(0.0200000+ 0.0197000*EXP(-  (STATES(1)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(1)+40.0000)/9.50000)*4.00000))
      RATES(24) = (ALGBRC(6) - STATES(24))/ALGBRC(21)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(2)+35.0000)/6.00000))
      ALGBRC(23) = 0.00100000/(0.0200000+ 0.0197000*EXP(-  (STATES(2)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(2)+40.0000)/9.50000)*4.00000))
      RATES(38) = (ALGBRC(8) - STATES(38))/ALGBRC(23)
      ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+21.5000)/7.50000))
      ALGBRC(24) = 0.00100000/(( 0.00138000*(STATES(1)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(1)+14.2000)))+( 0.000610000*(STATES(1)+38.9000))/(EXP( 0.145000*(STATES(1)+38.9000)) - 1.00000))
      RATES(39) = (ALGBRC(9) - STATES(39))/ALGBRC(24)
      ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(2)+21.5000)/7.50000))
      ALGBRC(25) = 0.00100000/(( 0.00138000*(STATES(2)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(2)+14.2000)))+( 0.000610000*(STATES(2)+38.9000))/(EXP( 0.145000*(STATES(2)+38.9000)) - 1.00000))
      RATES(40) = (ALGBRC(10) - STATES(40))/ALGBRC(25)
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
      ALGBRC(26) = 0.00100000/(( 7.19000e-05*(STATES(1)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(1)+30.0000)))+( 0.000131000*(STATES(1)+30.0000))/(EXP( 0.0687000*(STATES(1)+30.0000)) - 1.00000))
      RATES(42) = (ALGBRC(11) - STATES(42))/ALGBRC(26)
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(2) - 1.50000)/16.7000))
      ALGBRC(27) = 0.00100000/(( 7.19000e-05*(STATES(2)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(2)+30.0000)))+( 0.000131000*(STATES(2)+30.0000))/(EXP( 0.0687000*(STATES(2)+30.0000)) - 1.00000))
      RATES(43) = (ALGBRC(12) - STATES(43))/ALGBRC(27)
      ALGBRC(28) = ALGBRC(13)
      RATES(45) = (ALGBRC(28) - STATES(45))/CONSTS(68)
      ALGBRC(29) = ALGBRC(14)
      RATES(47) = (ALGBRC(29) - STATES(47))/CONSTS(71)
      ALGBRC(5) =  400.000*EXP((STATES(1)+12.0000)/10.0000)
      ALGBRC(20) =  50.0000*EXP(- (STATES(1)+12.0000)/13.0000)
      ALGBRC(31) =  187.500*STATES(23)
      RATES(11) = ( ALGBRC(20)*STATES(12)+ CONSTS(50)*STATES(17)) -  ( 4.00000*ALGBRC(5)+ALGBRC(31))*STATES(11)
      RATES(12) = ( 4.00000*ALGBRC(5)*STATES(11)+ 2.00000*ALGBRC(20)*STATES(13)+ (CONSTS(50)/CONSTS(52))*STATES(18)) -  (ALGBRC(20)+ 3.00000*ALGBRC(5)+ ALGBRC(31)*CONSTS(51))*STATES(12)
      RATES(13) = ( 3.00000*ALGBRC(5)*STATES(12)+ 3.00000*ALGBRC(20)*STATES(14)+ (CONSTS(50)/CONSTS(52) ** 2.00000)*STATES(19)) -  ( 2.00000*ALGBRC(20)+ 2.00000*ALGBRC(5)+ ALGBRC(31)*CONSTS(51) ** 2.00000)*STATES(13)
      RATES(14) = ( 2.00000*ALGBRC(5)*STATES(13)+ 4.00000*ALGBRC(20)*STATES(15)+ (CONSTS(50)/CONSTS(52) ** 3.00000)*STATES(20)) -  ( 3.00000*ALGBRC(20)+ALGBRC(5)+ ALGBRC(31)*CONSTS(51) ** 3.00000)*STATES(14)
      RATES(15) = ( ALGBRC(5)*STATES(14)+ CONSTS(54)*STATES(16)+ (CONSTS(50)/CONSTS(52) ** 4.00000)*STATES(21)) -  ( 4.00000*ALGBRC(20)+CONSTS(53)+ ALGBRC(31)*CONSTS(51) ** 4.00000)*STATES(15)
      ALGBRC(22) =  400.000*EXP((STATES(2)+12.0000)/10.0000)
      ALGBRC(32) =  50.0000*EXP(- (STATES(2)+12.0000)/13.0000)
      ALGBRC(7) =  187.500*STATES(23)
      RATES(25) = ( ALGBRC(32)*STATES(26)+ CONSTS(57)*STATES(31)) -  ( 4.00000*ALGBRC(22)+ALGBRC(7))*STATES(25)
      RATES(26) = ( 4.00000*ALGBRC(22)*STATES(25)+ 2.00000*ALGBRC(32)*STATES(27)+ (CONSTS(57)/CONSTS(59))*STATES(32)) -  (ALGBRC(32)+ 3.00000*ALGBRC(22)+ ALGBRC(7)*CONSTS(58))*STATES(26)
      RATES(27) = ( 3.00000*ALGBRC(22)*STATES(26)+ 3.00000*ALGBRC(32)*STATES(28)+ (CONSTS(57)/CONSTS(59) ** 2.00000)*STATES(33)) -  ( 2.00000*ALGBRC(32)+ 2.00000*ALGBRC(22)+ ALGBRC(7)*CONSTS(58) ** 2.00000)*STATES(27)
      RATES(28) = ( 2.00000*ALGBRC(22)*STATES(27)+ 4.00000*ALGBRC(32)*STATES(29)+ (CONSTS(57)/CONSTS(59) ** 3.00000)*STATES(34)) -  ( 3.00000*ALGBRC(32)+ALGBRC(22)+ ALGBRC(7)*CONSTS(58) ** 3.00000)*STATES(28)
      RATES(29) = ( ALGBRC(22)*STATES(28)+ CONSTS(61)*STATES(30)+ (CONSTS(57)/CONSTS(59) ** 4.00000)*STATES(35)) -  ( 4.00000*ALGBRC(32)+CONSTS(60)+ ALGBRC(7)*CONSTS(58) ** 4.00000)*STATES(29)
      ALGBRC(34) =  ALGBRC(5)*CONSTS(51)
      ALGBRC(37) = ALGBRC(20)/CONSTS(52)
      RATES(17) = ( ALGBRC(37)*STATES(18)+ ALGBRC(31)*STATES(11)) -  ( 4.00000*ALGBRC(34)+CONSTS(50))*STATES(17)
      RATES(18) = ( 4.00000*ALGBRC(34)*STATES(17)+ 2.00000*ALGBRC(37)*STATES(19)+ ALGBRC(31)*CONSTS(51)*STATES(12)) -  (ALGBRC(37)+ 3.00000*ALGBRC(34)+CONSTS(50)/CONSTS(52))*STATES(18)
      RATES(19) = ( 3.00000*ALGBRC(34)*STATES(18)+ 3.00000*ALGBRC(37)*STATES(20)+ ALGBRC(31)*CONSTS(51) ** 2.00000*STATES(13)) -  ( 2.00000*ALGBRC(37)+ 2.00000*ALGBRC(34)+CONSTS(50)/CONSTS(52) ** 2.00000)*STATES(19)
      RATES(20) = ( 2.00000*ALGBRC(34)*STATES(19)+ 4.00000*ALGBRC(37)*STATES(21)+ ALGBRC(31)*CONSTS(51) ** 3.00000*STATES(14)) -  ( 3.00000*ALGBRC(37)+ALGBRC(34)+CONSTS(50)/CONSTS(52) ** 3.00000)*STATES(20)
      RATES(21) = ( ALGBRC(34)*STATES(20)+ CONSTS(56)*STATES(22)+ ALGBRC(31)*CONSTS(51) ** 4.00000*STATES(15)) -  ( 4.00000*ALGBRC(37)+CONSTS(55)+CONSTS(50)/CONSTS(52) ** 4.00000)*STATES(21)
      ALGBRC(35) =  ALGBRC(22)*CONSTS(58)
      ALGBRC(38) = ALGBRC(32)/CONSTS(59)
      RATES(31) = ( ALGBRC(38)*STATES(32)+ ALGBRC(7)*STATES(25)) -  ( 4.00000*ALGBRC(35)+CONSTS(57))*STATES(31)
      RATES(32) = ( 4.00000*ALGBRC(35)*STATES(31)+ 2.00000*ALGBRC(38)*STATES(33)+ ALGBRC(7)*CONSTS(58)*STATES(26)) -  (ALGBRC(38)+ 3.00000*ALGBRC(35)+CONSTS(57)/CONSTS(59))*STATES(32)
      RATES(33) = ( 3.00000*ALGBRC(35)*STATES(32)+ 3.00000*ALGBRC(38)*STATES(34)+ ALGBRC(7)*CONSTS(58) ** 2.00000*STATES(27)) -  ( 2.00000*ALGBRC(38)+ 2.00000*ALGBRC(35)+CONSTS(57)/CONSTS(59) ** 2.00000)*STATES(33)
      RATES(34) = ( 2.00000*ALGBRC(35)*STATES(33)+ 4.00000*ALGBRC(38)*STATES(35)+ ALGBRC(7)*CONSTS(58) ** 3.00000*STATES(28)) -  ( 3.00000*ALGBRC(38)+ALGBRC(35)+CONSTS(57)/CONSTS(59) ** 3.00000)*STATES(34)
      RATES(35) = ( ALGBRC(35)*STATES(34)+ CONSTS(63)*STATES(36)+ ALGBRC(7)*CONSTS(58) ** 4.00000*STATES(29)) -  ( 4.00000*ALGBRC(38)+CONSTS(62)+CONSTS(57)/CONSTS(59) ** 4.00000)*STATES(35)
      ALGBRC(36) =  CONSTS(92)*STATES(5) ** 3.00000*STATES(6)*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(40) =  CONSTS(110)*STATES(9) ** 3.00000*STATES(10)*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(41) =  (CONSTS(93)/(1.00000+EXP((- 54.0000 - STATES(1))/8.00000)))*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(42) =  (CONSTS(111)/(1.00000+EXP((- 54.0000 - STATES(2))/8.00000)))*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(79) = (( (( 1.00000*CONSTS(108)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*CONSTS(4)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(80) = (( (( 1.00000*CONSTS(126)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*STATES(7)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(33) = ( log(CONSTS(4)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(83) =  CONSTS(101)*(STATES(1) - ALGBRC(33))
      ALGBRC(39) = ( log(STATES(7)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(84) =  CONSTS(119)*(STATES(2) - ALGBRC(39))
      ALGBRC(89) = ( CONSTS(105)*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5) -  CONSTS(4) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5)+ CONSTS(4) ** 3.00000*STATES(41)))
      ALGBRC(90) = ( CONSTS(123)*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37) -  STATES(7) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37)+ STATES(7) ** 3.00000*STATES(41)))
      ALGBRC(91) = ( (( (CONSTS(104)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(CONSTS(4)/67.3000) - 1.00000)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*CONSTS(6))/(CONSTS(6)+1.50000)
      ALGBRC(92) = ( (( (CONSTS(122)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(STATES(7)/67.3000) - 1.00000)*EXP(( - STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*STATES(8))/(STATES(8)+1.50000)
      RATES(3) = - (ALGBRC(36)+ALGBRC(40)+ALGBRC(41)+ALGBRC(42)+ALGBRC(79)+ALGBRC(80)+ALGBRC(83)+ALGBRC(84)+ 3.00000*ALGBRC(89)+ 3.00000*ALGBRC(90)+ 3.00000*ALGBRC(91)+ 3.00000*ALGBRC(92))/( CONSTS(3)*CONSTS(87))
      ALGBRC(93) =  (( 1.00000*CONSTS(128))/CONSTS(72))*(STATES(7) - CONSTS(4))
      RATES(7) = ((ALGBRC(40)+ALGBRC(42)+ALGBRC(80)+ALGBRC(84)+ 3.00000*ALGBRC(90)+ 3.00000*ALGBRC(92))/CONSTS(3) -  1.00000*ALGBRC(93))/CONSTS(128)
      ALGBRC(44) = ( (( 1.00000*CONSTS(124)*4.00000*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.341000*STATES(37)))/(EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(87) = ( log(STATES(37)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
      ALGBRC(88) =  CONSTS(120)*(STATES(2) - ALGBRC(87))
      ALGBRC(95) = ( CONSTS(121)*STATES(41))/(0.000500000+STATES(41))
      ALGBRC(97) =  (( 1.00000*CONSTS(128))/CONSTS(73))*(STATES(37) - CONSTS(5))
      RATES(37) = (( - 2.00000*ALGBRC(90)+ALGBRC(44)+ALGBRC(88)+ALGBRC(95))/( 2.00000*CONSTS(3)) -  1.00000*ALGBRC(97))/CONSTS(128)
      ALGBRC(15) =  1000.00*(STATES(2)/CONSTS(90) - STATES(1)/CONSTS(90))
      ALGBRC(30) = TERNRY(VOI -  INT(VOI/CONSTS(49))*CONSTS(49).GE.0.00000.AND.VOI -  INT(VOI/CONSTS(49))*CONSTS(49).LE.0.00100000,  45.0000*(CONSTS(81)+CONSTS(91)), 0.00000)
      ALGBRC(43) = ( (( 1.00000*CONSTS(106)*4.00000*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.341000*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(45) = ( (( (( 1.00000*CONSTS(107))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(47) = ( log(CONSTS(6)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(48) = 1.00000/(1.00000+EXP((STATES(1)+9.00000)/22.4000))
      ALGBRC(49) =  CONSTS(94)*0.0261400* (CONSTS(6)/5.40000) ** (1.0 / 2)*STATES(39)*ALGBRC(48)*(STATES(1) - ALGBRC(47))
      ALGBRC(53) = ( log((CONSTS(6)+ CONSTS(64)*CONSTS(4))/(STATES(4)+ CONSTS(64)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(54) =  CONSTS(95)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(42) ** 2.00000*(STATES(1) - ALGBRC(53))
      ALGBRC(57) = ALGBRC(47)
      ALGBRC(58) = 1020.00/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(57)) - 59.2150)))
      ALGBRC(59) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(57))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(57)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(57))+4.75300)))
      ALGBRC(60) =  (( CONSTS(98)* (CONSTS(6)/5.40000) ** (1.0 / 2)*ALGBRC(58))/(ALGBRC(58)+ALGBRC(59)))*(STATES(1) - ALGBRC(57))
      ALGBRC(66) = 1.00000/(1.00000+EXP((20.0000 - STATES(1))/5.00000))
      ALGBRC(65) = ALGBRC(47)
      ALGBRC(67) =  CONSTS(96)*ALGBRC(66)*(STATES(1) - ALGBRC(65))
      ALGBRC(71) = ALGBRC(47)
      ALGBRC(72) =  (CONSTS(97)/( (1.00000+EXP((STATES(1)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(1) - 45.8000)/25.8700))))*( STATES(44)*CONSTS(66)+ STATES(45)*(1.00000 - CONSTS(66)))*(STATES(1) - ALGBRC(71))
      ALGBRC(75) = ALGBRC(47)
      ALGBRC(76) =  (( CONSTS(99)*(0.800000 - 0.650000/(1.00000+EXP((STATES(1)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(1) - ALGBRC(75))
      ALGBRC(81) = (( (( 1.00000*CONSTS(109)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(85) = ( log(CONSTS(5)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
      ALGBRC(86) =  CONSTS(102)*(STATES(1) - ALGBRC(85))
      ALGBRC(94) = ( CONSTS(103)*STATES(41))/(0.000500000+STATES(41))
      ALGBRC(96) = ALGBRC(47)
      ALGBRC(99) =  (( CONSTS(100)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*CONSTS(6)/4.00000 ** 0.240000*(STATES(1) - ALGBRC(96))
      RATES(1) = ( 1000.00*(((((((((((((((((((ALGBRC(30)+ALGBRC(15)) - ALGBRC(36)) - ALGBRC(41)) - ALGBRC(43)) - ALGBRC(45)) - ALGBRC(54)) - ALGBRC(49)) - ALGBRC(60)) - ALGBRC(67)) - ALGBRC(76)) - ALGBRC(99)) - ALGBRC(79)) - ALGBRC(81)) - ALGBRC(83)) - ALGBRC(86)) - ALGBRC(89)) - ALGBRC(91)) - ALGBRC(94)) - ALGBRC(72)))/CONSTS(130)
      ALGBRC(98) = ( 1.00000*1000.00*CONSTS(87)*0.00180000*STATES(41) ** 2.00000)/(STATES(41) ** 2.00000+0.000500000 ** 2.00000)
      ALGBRC(100) =  1250.00*CONSTS(87)*5.80000e-05*(STATES(48) - STATES(41))
      ALGBRC(102) =  (( 1.00000*CONSTS(86))/CONSTS(75))*(STATES(48) - STATES(49))
      RATES(48) = ((ALGBRC(98) - ALGBRC(100)) - ALGBRC(102))/( 1.00000*CONSTS(88))
      ALGBRC(104) =  1800.00*CONSTS(86)*(STATES(51)+STATES(52))*(STATES(49) - STATES(23))
      RATES(49) = ( (1.00000/(1.00000+( 15.0000*0.800000)/0.800000+STATES(49) ** 2.00000))*(ALGBRC(102) - ALGBRC(104)))/CONSTS(86)
      ALGBRC(46) = ( (( (( 1.00000*CONSTS(125))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(50) = ( log(STATES(8)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(2)+9.00000)/22.4000))
      ALGBRC(52) =  CONSTS(112)*0.0261400* (STATES(8)/5.40000) ** (1.0 / 2)*STATES(40)*ALGBRC(51)*(STATES(2) - ALGBRC(50))
      ALGBRC(55) = ( log((STATES(8)+ CONSTS(65)*STATES(7))/(STATES(4)+ CONSTS(65)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(56) =  CONSTS(113)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(43) ** 2.00000*(STATES(2) - ALGBRC(55))
      ALGBRC(61) = ALGBRC(50)
      ALGBRC(62) = 1020.00/(1.00000+EXP( 0.238500*((STATES(2) - ALGBRC(61)) - 59.2150)))
      ALGBRC(63) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(2) - ALGBRC(61))+5.47600))+EXP( 0.0617500*((STATES(2) - ALGBRC(61)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(2) - ALGBRC(61))+4.75300)))
      ALGBRC(64) =  (( CONSTS(116)* (STATES(8)/5.40000) ** (1.0 / 2)*ALGBRC(62))/(ALGBRC(62)+ALGBRC(63)))*(STATES(2) - ALGBRC(61))
      ALGBRC(69) = 1.00000/(1.00000+EXP((20.0000 - STATES(2))/5.00000))
      ALGBRC(68) = ALGBRC(50)
      ALGBRC(70) =  CONSTS(114)*ALGBRC(69)*(STATES(2) - ALGBRC(68))
      ALGBRC(73) = ALGBRC(50)
      ALGBRC(74) =  (CONSTS(115)/( (1.00000+EXP((STATES(2)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(2) - 45.8000)/25.8700))))*( STATES(46)*CONSTS(69)+ STATES(47)*(1.00000 - CONSTS(69)))*(STATES(2) - ALGBRC(73))
      ALGBRC(77) = ALGBRC(50)
      ALGBRC(78) =  (( CONSTS(117)*(0.800000 - 0.650000/(1.00000+EXP((STATES(2)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(2) - ALGBRC(77))
      ALGBRC(82) = (( (( 1.00000*CONSTS(127)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(101) = ALGBRC(50)
      ALGBRC(103) =  (( CONSTS(118)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*STATES(8)/4.00000 ** 0.240000*(STATES(2) - ALGBRC(101))
      RATES(4) = - (((( - 0.00000*ALGBRC(30)+ALGBRC(54)+ALGBRC(56)+ALGBRC(49)+ALGBRC(52)+ALGBRC(60)+ALGBRC(64)+ALGBRC(67)+ALGBRC(70)+ALGBRC(45)+ALGBRC(46)+ALGBRC(76)+ALGBRC(78)+ALGBRC(81)+ALGBRC(82)) -  2.00000*ALGBRC(91)) -  2.00000*ALGBRC(92))+ALGBRC(99)+ALGBRC(103)+ALGBRC(72)+ALGBRC(74))/( CONSTS(3)*CONSTS(87))
      RATES(2) = ( 1000.00*((((((((((((((((((- ALGBRC(15) - ALGBRC(40)) - ALGBRC(42)) - ALGBRC(44)) - ALGBRC(46)) - ALGBRC(56)) - ALGBRC(52)) - ALGBRC(64)) - ALGBRC(70)) - ALGBRC(78)) - ALGBRC(103)) - ALGBRC(80)) - ALGBRC(82)) - ALGBRC(84)) - ALGBRC(88)) - ALGBRC(90)) - ALGBRC(92)) - ALGBRC(95)) - ALGBRC(74)))/CONSTS(129)
      ALGBRC(106) = ( 1.00000*(STATES(23) - STATES(41))*CONSTS(87))/CONSTS(80)
      RATES(23) =  (1.00000/(1.00000+( 0.0500000*0.00238000)/0.00238000+STATES(23) ** 2.00000))*((- (ALGBRC(43)+ALGBRC(44))/( 2.00000*CONSTS(3)*CONSTS(85))+( 1.00000*ALGBRC(104))/CONSTS(85)) - ( 1.00000*ALGBRC(106))/CONSTS(85))
      ALGBRC(105) =  (( 1.00000*CONSTS(128))/CONSTS(74))*(STATES(8) - CONSTS(6))
      RATES(8) = ((((ALGBRC(56)+ALGBRC(52)+ALGBRC(64)+ALGBRC(70)+ALGBRC(78)+ALGBRC(46)+ALGBRC(82)) -  2.00000*ALGBRC(92))+ALGBRC(103)+ALGBRC(74))/CONSTS(3) -  1.00000*ALGBRC(105))/CONSTS(128)
      ALGBRC(107) =  20000.0*STATES(41)*(1.00000 - STATES(54)) -  0.0700000*STATES(54)
      RATES(54) = ALGBRC(107)
      ALGBRC(108) =  40000.0*STATES(41)*(1.00000 - STATES(55)) -  STATES(55)*40.0000
      RATES(41) =  (1.00000/(1.00000+( 0.0500000*0.00238000)/0.00238000+STATES(41) ** 2.00000))*(((( - 2.00000*(ALGBRC(89)+ALGBRC(90))+ALGBRC(86)+ALGBRC(88)+ALGBRC(94)+ALGBRC(95))/( - 2.00000*CONSTS(3)*CONSTS(87))+( 1.00000*((ALGBRC(106)+ALGBRC(100)) - ALGBRC(98)))/CONSTS(87)) -  ALGBRC(107)*0.140000) -  ALGBRC(108)*0.0700000)
      RATES(55) = ALGBRC(108)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(13) = 1.00600/(1.00000+EXP((STATES(1)+36.6900+10.0000)/12.4300))
      ALGBRC(14) = 1.00600/(1.00000+EXP((STATES(2)+36.6900+10.0000)/12.4300))
      ALGBRC(1) = 1.00000/(1.00000+EXP(((STATES(1)+52.2000) - 9.00000)/- 7.40000))
      ALGBRC(16) = (0.00100000/( 101.600*EXP( STATES(1)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(1)))+0.000100000)/4.83000
      ALGBRC(2) = 1.00000/(1.00000+EXP(((STATES(1)+85.6000) - 9.00000)/5.50000))
      ALGBRC(17) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(1))+ 6.53700*EXP( STATES(1)*0.0801600))+0.000500000)/4.83000
      ALGBRC(3) = 1.00000/(1.00000+EXP(((STATES(2)+52.2000) - 9.00000)/- 7.40000))
      ALGBRC(18) = (0.00100000/( 101.600*EXP( STATES(2)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(2)))+0.000100000)/4.83000
      ALGBRC(4) = 1.00000/(1.00000+EXP(((STATES(2)+85.6000) - 9.00000)/5.50000))
      ALGBRC(19) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(2))+ 6.53700*EXP( STATES(2)*0.0801600))+0.000500000)/4.83000
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/6.00000))
      ALGBRC(21) = 0.00100000/(0.0200000+ 0.0197000*EXP(-  (STATES(1)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(1)+40.0000)/9.50000)*4.00000))
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(2)+35.0000)/6.00000))
      ALGBRC(23) = 0.00100000/(0.0200000+ 0.0197000*EXP(-  (STATES(2)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(2)+40.0000)/9.50000)*4.00000))
      ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+21.5000)/7.50000))
      ALGBRC(24) = 0.00100000/(( 0.00138000*(STATES(1)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(1)+14.2000)))+( 0.000610000*(STATES(1)+38.9000))/(EXP( 0.145000*(STATES(1)+38.9000)) - 1.00000))
      ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(2)+21.5000)/7.50000))
      ALGBRC(25) = 0.00100000/(( 0.00138000*(STATES(2)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(2)+14.2000)))+( 0.000610000*(STATES(2)+38.9000))/(EXP( 0.145000*(STATES(2)+38.9000)) - 1.00000))
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
      ALGBRC(26) = 0.00100000/(( 7.19000e-05*(STATES(1)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(1)+30.0000)))+( 0.000131000*(STATES(1)+30.0000))/(EXP( 0.0687000*(STATES(1)+30.0000)) - 1.00000))
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(2) - 1.50000)/16.7000))
      ALGBRC(27) = 0.00100000/(( 7.19000e-05*(STATES(2)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(2)+30.0000)))+( 0.000131000*(STATES(2)+30.0000))/(EXP( 0.0687000*(STATES(2)+30.0000)) - 1.00000))
      ALGBRC(28) = ALGBRC(13)
      ALGBRC(29) = ALGBRC(14)
      ALGBRC(5) =  400.000*EXP((STATES(1)+12.0000)/10.0000)
      ALGBRC(20) =  50.0000*EXP(- (STATES(1)+12.0000)/13.0000)
      ALGBRC(31) =  187.500*STATES(23)
      ALGBRC(22) =  400.000*EXP((STATES(2)+12.0000)/10.0000)
      ALGBRC(32) =  50.0000*EXP(- (STATES(2)+12.0000)/13.0000)
      ALGBRC(7) =  187.500*STATES(23)
      ALGBRC(34) =  ALGBRC(5)*CONSTS(51)
      ALGBRC(37) = ALGBRC(20)/CONSTS(52)
      ALGBRC(35) =  ALGBRC(22)*CONSTS(58)
      ALGBRC(38) = ALGBRC(32)/CONSTS(59)
      ALGBRC(36) =  CONSTS(92)*STATES(5) ** 3.00000*STATES(6)*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(40) =  CONSTS(110)*STATES(9) ** 3.00000*STATES(10)*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(41) =  (CONSTS(93)/(1.00000+EXP((- 54.0000 - STATES(1))/8.00000)))*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(42) =  (CONSTS(111)/(1.00000+EXP((- 54.0000 - STATES(2))/8.00000)))*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
      ALGBRC(79) = (( (( 1.00000*CONSTS(108)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*CONSTS(4)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(80) = (( (( 1.00000*CONSTS(126)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*STATES(7)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(33) = ( log(CONSTS(4)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(83) =  CONSTS(101)*(STATES(1) - ALGBRC(33))
      ALGBRC(39) = ( log(STATES(7)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(84) =  CONSTS(119)*(STATES(2) - ALGBRC(39))
      ALGBRC(89) = ( CONSTS(105)*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5) -  CONSTS(4) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5)+ CONSTS(4) ** 3.00000*STATES(41)))
      ALGBRC(90) = ( CONSTS(123)*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37) -  STATES(7) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37)+ STATES(7) ** 3.00000*STATES(41)))
      ALGBRC(91) = ( (( (CONSTS(104)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(CONSTS(4)/67.3000) - 1.00000)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*CONSTS(6))/(CONSTS(6)+1.50000)
      ALGBRC(92) = ( (( (CONSTS(122)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(STATES(7)/67.3000) - 1.00000)*EXP(( - STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*STATES(8))/(STATES(8)+1.50000)
      ALGBRC(93) =  (( 1.00000*CONSTS(128))/CONSTS(72))*(STATES(7) - CONSTS(4))
      ALGBRC(44) = ( (( 1.00000*CONSTS(124)*4.00000*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.341000*STATES(37)))/(EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(87) = ( log(STATES(37)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
      ALGBRC(88) =  CONSTS(120)*(STATES(2) - ALGBRC(87))
      ALGBRC(95) = ( CONSTS(121)*STATES(41))/(0.000500000+STATES(41))
      ALGBRC(97) =  (( 1.00000*CONSTS(128))/CONSTS(73))*(STATES(37) - CONSTS(5))
      ALGBRC(15) =  1000.00*(STATES(2)/CONSTS(90) - STATES(1)/CONSTS(90))
      ALGBRC(30) = TERNRY(VOI -  INT(VOI/CONSTS(49))*CONSTS(49).GE.0.00000.AND.VOI -  INT(VOI/CONSTS(49))*CONSTS(49).LE.0.00100000,  45.0000*(CONSTS(81)+CONSTS(91)), 0.00000)
      ALGBRC(43) = ( (( 1.00000*CONSTS(106)*4.00000*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.341000*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(45) = ( (( (( 1.00000*CONSTS(107))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(47) = ( log(CONSTS(6)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(48) = 1.00000/(1.00000+EXP((STATES(1)+9.00000)/22.4000))
      ALGBRC(49) =  CONSTS(94)*0.0261400* (CONSTS(6)/5.40000) ** (1.0 / 2)*STATES(39)*ALGBRC(48)*(STATES(1) - ALGBRC(47))
      ALGBRC(53) = ( log((CONSTS(6)+ CONSTS(64)*CONSTS(4))/(STATES(4)+ CONSTS(64)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(54) =  CONSTS(95)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(42) ** 2.00000*(STATES(1) - ALGBRC(53))
      ALGBRC(57) = ALGBRC(47)
      ALGBRC(58) = 1020.00/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(57)) - 59.2150)))
      ALGBRC(59) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(57))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(57)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(57))+4.75300)))
      ALGBRC(60) =  (( CONSTS(98)* (CONSTS(6)/5.40000) ** (1.0 / 2)*ALGBRC(58))/(ALGBRC(58)+ALGBRC(59)))*(STATES(1) - ALGBRC(57))
      ALGBRC(66) = 1.00000/(1.00000+EXP((20.0000 - STATES(1))/5.00000))
      ALGBRC(65) = ALGBRC(47)
      ALGBRC(67) =  CONSTS(96)*ALGBRC(66)*(STATES(1) - ALGBRC(65))
      ALGBRC(71) = ALGBRC(47)
      ALGBRC(72) =  (CONSTS(97)/( (1.00000+EXP((STATES(1)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(1) - 45.8000)/25.8700))))*( STATES(44)*CONSTS(66)+ STATES(45)*(1.00000 - CONSTS(66)))*(STATES(1) - ALGBRC(71))
      ALGBRC(75) = ALGBRC(47)
      ALGBRC(76) =  (( CONSTS(99)*(0.800000 - 0.650000/(1.00000+EXP((STATES(1)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(1) - ALGBRC(75))
      ALGBRC(81) = (( (( 1.00000*CONSTS(109)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(85) = ( log(CONSTS(5)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
      ALGBRC(86) =  CONSTS(102)*(STATES(1) - ALGBRC(85))
      ALGBRC(94) = ( CONSTS(103)*STATES(41))/(0.000500000+STATES(41))
      ALGBRC(96) = ALGBRC(47)
      ALGBRC(99) =  (( CONSTS(100)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*CONSTS(6)/4.00000 ** 0.240000*(STATES(1) - ALGBRC(96))
      ALGBRC(98) = ( 1.00000*1000.00*CONSTS(87)*0.00180000*STATES(41) ** 2.00000)/(STATES(41) ** 2.00000+0.000500000 ** 2.00000)
      ALGBRC(100) =  1250.00*CONSTS(87)*5.80000e-05*(STATES(48) - STATES(41))
      ALGBRC(102) =  (( 1.00000*CONSTS(86))/CONSTS(75))*(STATES(48) - STATES(49))
      ALGBRC(104) =  1800.00*CONSTS(86)*(STATES(51)+STATES(52))*(STATES(49) - STATES(23))
      ALGBRC(46) = ( (( (( 1.00000*CONSTS(125))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
      ALGBRC(50) = ( log(STATES(8)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(2)+9.00000)/22.4000))
      ALGBRC(52) =  CONSTS(112)*0.0261400* (STATES(8)/5.40000) ** (1.0 / 2)*STATES(40)*ALGBRC(51)*(STATES(2) - ALGBRC(50))
      ALGBRC(55) = ( log((STATES(8)+ CONSTS(65)*STATES(7))/(STATES(4)+ CONSTS(65)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
      ALGBRC(56) =  CONSTS(113)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(43) ** 2.00000*(STATES(2) - ALGBRC(55))
      ALGBRC(61) = ALGBRC(50)
      ALGBRC(62) = 1020.00/(1.00000+EXP( 0.238500*((STATES(2) - ALGBRC(61)) - 59.2150)))
      ALGBRC(63) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(2) - ALGBRC(61))+5.47600))+EXP( 0.0617500*((STATES(2) - ALGBRC(61)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(2) - ALGBRC(61))+4.75300)))
      ALGBRC(64) =  (( CONSTS(116)* (STATES(8)/5.40000) ** (1.0 / 2)*ALGBRC(62))/(ALGBRC(62)+ALGBRC(63)))*(STATES(2) - ALGBRC(61))
      ALGBRC(69) = 1.00000/(1.00000+EXP((20.0000 - STATES(2))/5.00000))
      ALGBRC(68) = ALGBRC(50)
      ALGBRC(70) =  CONSTS(114)*ALGBRC(69)*(STATES(2) - ALGBRC(68))
      ALGBRC(73) = ALGBRC(50)
      ALGBRC(74) =  (CONSTS(115)/( (1.00000+EXP((STATES(2)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(2) - 45.8000)/25.8700))))*( STATES(46)*CONSTS(69)+ STATES(47)*(1.00000 - CONSTS(69)))*(STATES(2) - ALGBRC(73))
      ALGBRC(77) = ALGBRC(50)
      ALGBRC(78) =  (( CONSTS(117)*(0.800000 - 0.650000/(1.00000+EXP((STATES(2)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(2) - ALGBRC(77))
      ALGBRC(82) = (( (( 1.00000*CONSTS(127)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) -  0.750000*STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
      ALGBRC(101) = ALGBRC(50)
      ALGBRC(103) =  (( CONSTS(118)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*STATES(8)/4.00000 ** 0.240000*(STATES(2) - ALGBRC(101))
      ALGBRC(106) = ( 1.00000*(STATES(23) - STATES(41))*CONSTS(87))/CONSTS(80)
      ALGBRC(105) =  (( 1.00000*CONSTS(128))/CONSTS(74))*(STATES(8) - CONSTS(6))
      ALGBRC(107) =  20000.0*STATES(41)*(1.00000 - STATES(54)) -  0.0700000*STATES(54)
      ALGBRC(108) =  40000.0*STATES(41)*(1.00000 - STATES(55)) -  STATES(55)*40.0000
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END