C
C There are a total of 51 entries in the algebraic variable array.
C There are a total of 13 entries in each of the rate and state variable arrays.
C There are a total of 58 entries in the constant variable array.
C
C
C VOI is t in component environment (ms).
C CONSTS(1) is R in component environment (mJ_per_mole_kelvin).
C CONSTS(2) is T in component environment (kelvin).
C CONSTS(3) is F in component environment (coulomb_per_mole).
C CONSTS(49) is V_tau in component environment (mV).
C CONSTS(4) is Ca_o in component environment (mM).
C CONSTS(5) is Na_o in component environment (mM).
C CONSTS(6) is K_o in component environment (mM).
C CONSTS(50) is vol_cyt in component environment (pl).
C CONSTS(7) is vol_pmu in component environment (pl).
C CONSTS(8) is fr_cyt in component environment (dimensionless).
C STATES(1) is V in component Membrane (mV).
C ALGBRC(1) is VD in component Membrane (dimensionless).
C CONSTS(9) is C_sp in component Membrane (pF_per_sqcm).
C STATES(2) is Ca_i in component Cytosol (mM).
C STATES(3) is Na_i in component Cytosol (mM).
C STATES(4) is K_i in component Cytosol (mM).
C ALGBRC(3) is V_Ca in component Membrane (dimensionless).
C ALGBRC(5) is V_Na in component Membrane (dimensionless).
C ALGBRC(7) is V_K in component Membrane (dimensionless).
C CONSTS(56) is A_pmu in component Membrane (sqcm).
C CONSTS(10) is SVR_pmu in component Membrane (per_cm).
C ALGBRC(41) is J_Ca in component calcium_dynamics (mM_per_ms).
C ALGBRC(49) is J_Na in component sodium_dynamics (mM_per_ms).
C ALGBRC(51) is J_K in component potassium_dynamics (mM_per_ms).
C CONSTS(11) is atp in component Cytosol (mM).
C ALGBRC(40) is J_ca in component calcium_dynamics (mM_per_ms).
C ALGBRC(20) is I_CaL in component L_type_Ca_channel (pA).
C ALGBRC(35) is I_pmca in component PMCA (pA).
C ALGBRC(39) is I_xm in component NaCa (pA).
C ALGBRC(16) is J_calb in component calcium_buffer_dynamics (mM_per_ms).
C ALGBRC(18) is J_cam in component calcium_buffer_dynamics (mM_per_ms).
C ALGBRC(22) is I_Na in component transient_Na_channel (pA).
C ALGBRC(23) is I_Nalk in component Leak_Na_channel (pA).
C ALGBRC(24) is I_NaHCN in component HCN_channel (pA).
C ALGBRC(47) is I_nk in component sodium_pump (pA).
C ALGBRC(31) is I_K in component potassium_dynamics (pA).
C ALGBRC(26) is I_Ksk in component SK_K_channel (pA).
C ALGBRC(28) is I_Kdr in component DR_K_channel (pA).
C ALGBRC(30) is I_Kir in component IR_K_channel (pA).
C STATES(5) is Calb in component calcium_buffer_dynamics (mM).
C STATES(6) is Cam in component calcium_buffer_dynamics (mM).
C CONSTS(12) is Calbtot in component calcium_buffer_dynamics (mM).
C CONSTS(13) is Camtot in component calcium_buffer_dynamics (mM).
C ALGBRC(15) is CaCalb in component calcium_buffer_dynamics (mM).
C ALGBRC(17) is CaCam in component calcium_buffer_dynamics (mM).
C CONSTS(14) is kcal_1 in component calcium_buffer_dynamics (per_mM_ms).
C CONSTS(15) is kcal_2 in component calcium_buffer_dynamics (per_ms).
C ALGBRC(9) is kcam_cb in component calcium_buffer_dynamics (per_ms).
C CONSTS(16) is kcam_cd in component calcium_buffer_dynamics (per_ms).
C ALGBRC(11) is kcam_nb in component calcium_buffer_dynamics (per_ms).
C CONSTS(17) is kcam_nd in component calcium_buffer_dynamics (per_ms).
C ALGBRC(13) is alpha_cam in component calcium_buffer_dynamics (per_ms).
C ALGBRC(14) is beta_cam in component calcium_buffer_dynamics (per_ms).
C STATES(7) is m_cal in component L_type_Ca_channel (dimensionless).
C ALGBRC(19) is h_cal in component L_type_Ca_channel (dimensionless).
C CONSTS(18) is g_cal in component L_type_Ca_channel (pA_per_mM).
C CONSTS(19) is g_na in component transient_Na_channel (pA_per_mM).
C ALGBRC(21) is O_na in component transient_Na_channel (dimensionless).
C STATES(8) is m_na in component transient_Na_channel (dimensionless).
C STATES(9) is h_na in component transient_Na_channel (dimensionless).
C CONSTS(20) is A_mna in component transient_Na_channel (per_ms).
C CONSTS(21) is B_mna in component transient_Na_channel (per_ms).
C CONSTS(22) is A_hna in component transient_Na_channel (per_ms).
C CONSTS(23) is B_hna in component transient_Na_channel (per_ms).
C CONSTS(24) is za_mna in component transient_Na_channel (dimensionless).
C CONSTS(25) is zb_mna in component transient_Na_channel (dimensionless).
C CONSTS(26) is za_hna in component transient_Na_channel (dimensionless).
C CONSTS(27) is zb_hna in component transient_Na_channel (dimensionless).
C CONSTS(28) is g_nalk in component Leak_Na_channel (pA_per_mM).
C CONSTS(29) is g_nahcn in component HCN_channel (pA_per_mM).
C STATES(10) is O_hcn in component HCN_channel (dimensionless).
C ALGBRC(8) is kf_hcn in component HCN_channel (per_ms).
C ALGBRC(12) is kr_hcn in component HCN_channel (per_ms).
C ALGBRC(2) is kf_free in component HCN_channel (per_ms).
C ALGBRC(4) is kr_free in component HCN_channel (per_ms).
C ALGBRC(6) is kf_bnd in component HCN_channel (per_ms).
C ALGBRC(10) is kr_bnd in component HCN_channel (per_ms).
C CONSTS(51) is P_c in component HCN_channel (dimensionless).
C CONSTS(52) is P_o in component HCN_channel (dimensionless).
C CONSTS(30) is cAMP in component HCN_channel (mM).
C ALGBRC(25) is O_sk in component SK_K_channel (dimensionless).
C CONSTS(31) is g_ksk in component SK_K_channel (pA_per_mM).
C ALGBRC(27) is O_kdr in component DR_K_channel (dimensionless).
C STATES(11) is m_kdr in component DR_K_channel (dimensionless).
C CONSTS(32) is g_kdr in component DR_K_channel (nS).
C ALGBRC(29) is O_kir in component IR_K_channel (dimensionless).
C CONSTS(33) is g_kir in component IR_K_channel (nS).
C STATES(12) is y_pc in component PMCA (dimensionless).
C ALGBRC(32) is K_pmca in component PMCA (pA).
C CONSTS(53) is k_1pc in component PMCA (per_ms).
C CONSTS(34) is k_2pc in component PMCA (per_ms).
C CONSTS(35) is k_3pc in component PMCA (per_ms).
C CONSTS(36) is k_4pc in component PMCA (per_ms).
C ALGBRC(34) is P_E1Spc in component PMCA (dimensionless).
C CONSTS(54) is P_E2Spc in component PMCA (dimensionless).
C ALGBRC(36) is P_E1pc in component PMCA (dimensionless).
C CONSTS(57) is P_E2pc in component PMCA (dimensionless).
C ALGBRC(38) is alpha_pc in component PMCA (per_ms).
C CONSTS(58) is beta_pc in component PMCA (per_ms).
C ALGBRC(33) is K_pci in component PMCA (mM).
C CONSTS(37) is K_pco in component PMCA (mM).
C CONSTS(38) is k_pmca in component PMCA (dimensionless).
C CONSTS(39) is del in component NaCa (dimensionless).
C CONSTS(40) is k_xm in component NaCa (pA).
C ALGBRC(37) is Dr in component NaCa (mM4).
C ALGBRC(42) is Na_eff in component sodium_pump (mM).
C STATES(13) is y_nk in component sodium_pump (dimensionless).
C ALGBRC(45) is alpha_nk in component sodium_pump (per_ms).
C ALGBRC(50) is beta_nk in component sodium_pump (per_ms).
C ALGBRC(43) is P_E1Snk in component sodium_pump (dimensionless).
C ALGBRC(46) is P_E2Snk in component sodium_pump (dimensionless).
C ALGBRC(44) is P_E1Dnk in component sodium_pump (dimensionless).
C ALGBRC(48) is P_E2Dnk in component sodium_pump (dimensionless).
C CONSTS(55) is k_1nk in component sodium_pump (per_ms).
C CONSTS(41) is k_2nk in component sodium_pump (per_ms).
C CONSTS(42) is k_3nk in component sodium_pump (per_ms).
C CONSTS(43) is k_4nk in component sodium_pump (per_ms).
C CONSTS(44) is K_nknai in component sodium_pump (mM).
C CONSTS(45) is K_nknao in component sodium_pump (mM).
C CONSTS(46) is K_nkki in component sodium_pump (mM).
C CONSTS(47) is K_nkko in component sodium_pump (mM).
C CONSTS(48) is k_nk in component sodium_pump (pA).
C RATES(1) is d/dt V in component Membrane (mV).
C RATES(2) is d/dt Ca_i in component Cytosol (mM).
C RATES(3) is d/dt Na_i in component Cytosol (mM).
C RATES(4) is d/dt K_i in component Cytosol (mM).
C RATES(5) is d/dt Calb in component calcium_buffer_dynamics (mM).
C RATES(6) is d/dt Cam in component calcium_buffer_dynamics (mM).
C RATES(7) is d/dt m_cal in component L_type_Ca_channel (dimensionless).
C RATES(8) is d/dt m_na in component transient_Na_channel (dimensionless).
C RATES(9) is d/dt h_na in component transient_Na_channel (dimensionless).
C RATES(10) is d/dt O_hcn in component HCN_channel (dimensionless).
C RATES(11) is d/dt m_kdr in component DR_K_channel (dimensionless).
C RATES(12) is d/dt y_pc in component PMCA (dimensionless).
C RATES(13) is d/dt y_nk in component sodium_pump (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 8314.472
      CONSTS(2) = 310.15
      CONSTS(3) = 96485.30929
      CONSTS(4) = 1.8
      CONSTS(5) = 137
      CONSTS(6) = 5.4
      CONSTS(7) = 5
      CONSTS(8) = 0.5
      STATES(1) = -49.42
      CONSTS(9) = 0.9e6
      STATES(2) = 0.000188
      STATES(3) = 4.6876
      STATES(4) = 126.05893
      CONSTS(10) = 1.6667e4
      CONSTS(11) = 2
      STATES(5) = 0.0026
      STATES(6) = 0.0222
      CONSTS(12) = 0.005
      CONSTS(13) = 0.0235
      CONSTS(14) = 10
      CONSTS(15) = 2e-3
      CONSTS(16) = 0.003
      CONSTS(17) = 3
      STATES(7) = 0.006271
      CONSTS(18) = 2101.2
      CONSTS(19) = 907.68
      STATES(8) = 0.0952
      STATES(9) = 0.1848
      CONSTS(20) = 1.9651
      CONSTS(21) = 0.0424
      CONSTS(22) = 9.566e-5
      CONSTS(23) = 0.5296
      CONSTS(24) = 1.7127
      CONSTS(25) = 1.5581
      CONSTS(26) = -2.4317
      CONSTS(27) = -1.1868
      CONSTS(28) = 0.0053
      CONSTS(29) = 51.1
      STATES(10) = 0.003
      CONSTS(30) = 1e-5
      CONSTS(31) = 2.2515
      STATES(11) = 0.0932
      CONSTS(32) = 31.237
      CONSTS(33) = 13.816
      STATES(12) = 0.483
      CONSTS(34) = 0.001
      CONSTS(35) = 0.001
      CONSTS(36) = 1
      CONSTS(37) = 2
      CONSTS(38) = 2.233
      CONSTS(39) = 0.35
      CONSTS(40) = 0.0166
      STATES(13) = 0.6213
      CONSTS(41) = 0.04
      CONSTS(42) = 0.01
      CONSTS(43) = 0.165
      CONSTS(44) = 4.05
      CONSTS(45) = 69.8
      CONSTS(46) = 32.88
      CONSTS(47) = 0.258
      CONSTS(48) = 1085.7
      CONSTS(49) = ( CONSTS(1)*CONSTS(2))/CONSTS(3)
      CONSTS(50) =  CONSTS(8)*CONSTS(7)
      CONSTS(51) = 1.00000/(1.00000+CONSTS(30)/0.00116300)
      CONSTS(52) = 1.00000/(1.00000+CONSTS(30)/1.45000e-05)
      CONSTS(53) = 1.00000/(1.00000+0.100000/CONSTS(11))
      CONSTS(54) = 1.00000/(1.00000+CONSTS(37)/CONSTS(4))
      CONSTS(55) = 0.370000/(1.00000+0.0940000/CONSTS(11))
      CONSTS(56) = ( CONSTS(10)*CONSTS(7)*0.00100000*0.00100000*0.00100000)/1.00000
      CONSTS(57) = 1.00000 - CONSTS(54)
      CONSTS(58) =  CONSTS(34)*CONSTS(54)+ CONSTS(36)*CONSTS(57)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(7) = (1.00000/(1.00000+EXP(- (STATES(1)+15.0000)/7.00000)) - STATES(7))/( 7.68000*EXP(- (STATES(1)+65.0000)/17.3300 ** 2.00000)+0.723100)
      RATES(11) = (1.00000/(1.00000+EXP(- (STATES(1)+25.0000)/12.0000)) - STATES(11))/(18.0000/(1.00000+EXP((STATES(1)+39.0000)/8.00000))+1.00000)
      ALGBRC(1) = STATES(1)/CONSTS(49)
      RATES(8) =  CONSTS(20)*EXP( CONSTS(24)*ALGBRC(1))*(1.00000 - STATES(8)) -  CONSTS(21)*EXP( - CONSTS(25)*ALGBRC(1))*STATES(8)
      RATES(9) =  CONSTS(22)*EXP( CONSTS(26)*ALGBRC(1))*(1.00000 - STATES(9)) -  CONSTS(23)*EXP( - CONSTS(27)*ALGBRC(1))*STATES(9)
      ALGBRC(2) = 0.00600000/(1.00000+EXP((STATES(1)+87.7000)/6.45000))
      ALGBRC(6) = 0.0268000/(1.00000+EXP((STATES(1)+94.2000)/13.3000))
      ALGBRC(8) =  ALGBRC(2)*CONSTS(51)+ ALGBRC(6)*(1.00000 - CONSTS(51))
      ALGBRC(4) = 0.0800000/(1.00000+EXP(- (STATES(1)+51.7000)/7.00000))
      ALGBRC(10) = 0.0800000/(1.00000+EXP(- (STATES(1)+35.5000)/7.00000))
      ALGBRC(12) =  ALGBRC(4)*CONSTS(52)+ ALGBRC(10)*(1.00000 - CONSTS(52))
      RATES(10) =  ALGBRC(8)*(1.00000 - STATES(10)) -  ALGBRC(12)*STATES(10)
      ALGBRC(15) = CONSTS(12) - STATES(5)
      ALGBRC(16) =  CONSTS(14)*STATES(5)*STATES(2) -  CONSTS(15)*ALGBRC(15)
      RATES(5) = - ALGBRC(16)
      ALGBRC(17) = CONSTS(13) - STATES(6)
      ALGBRC(9) =  12000.0*STATES(2) ** 2.00000
      ALGBRC(11) =  3.70000e+06*STATES(2) ** 2.00000
      ALGBRC(13) =  ALGBRC(9)*ALGBRC(11)*(1.00000/(ALGBRC(9)+CONSTS(17))+1.00000/(CONSTS(16)+CONSTS(17)))
      ALGBRC(14) =  CONSTS(16)*CONSTS(17)*(1.00000/(ALGBRC(9)+CONSTS(17))+1.00000/(CONSTS(16)+CONSTS(17)))
      ALGBRC(18) =  ALGBRC(13)*STATES(6) -  ALGBRC(14)*ALGBRC(17)
      RATES(6) = - ALGBRC(18)
      ALGBRC(33) =  (173.600/(1.00000+ALGBRC(17)/5.00000e-05)+6.40000)*1.00000e-05
      ALGBRC(34) = 1.00000/(1.00000+ALGBRC(33)/STATES(2))
      ALGBRC(36) = 1.00000 - ALGBRC(34)
      ALGBRC(38) =  CONSTS(53)*ALGBRC(34)+ CONSTS(35)*ALGBRC(36)
      RATES(12) =  CONSTS(58)*(1.00000 - STATES(12)) -  ALGBRC(38)*STATES(12)
      ALGBRC(3) =  0.500000*log(CONSTS(4)/STATES(2))
      ALGBRC(19) = 0.000450000/(0.000450000+STATES(2))
      ALGBRC(20) = ( CONSTS(18)*STATES(7)*ALGBRC(19)* ( STATES(2)*CONSTS(4)) ** (1.0 / 2)* sinh(ALGBRC(1) - ALGBRC(3)))/( sinh(ALGBRC(1))/ALGBRC(1))
      ALGBRC(32) =  CONSTS(38)*(( 10.5600*ALGBRC(17))/(ALGBRC(17)+5.00000e-05)+1.20000)
      ALGBRC(35) =  ALGBRC(32)*( CONSTS(53)*ALGBRC(34)*STATES(12) -  CONSTS(34)*CONSTS(54)*(1.00000 - STATES(12)))*1.00000
      ALGBRC(37) =  (1.00000+ 0.00100000*( STATES(3) ** 3.00000*CONSTS(4)+ CONSTS(5) ** 3.00000*STATES(2)))*(1.00000+STATES(2)/0.00690000)
      ALGBRC(39) = ( CONSTS(40)*( STATES(3) ** 3.00000*CONSTS(4)*EXP( CONSTS(39)*ALGBRC(1)) -  CONSTS(5) ** 3.00000*STATES(2)*EXP( (CONSTS(39) - 1.00000)*ALGBRC(1))))/ALGBRC(37)
      ALGBRC(40) =  (- 1.00000/( 2.00000*CONSTS(3)*CONSTS(50)))*((ALGBRC(20)+ 2.00000*ALGBRC(35)) -  2.00000*ALGBRC(39))
      ALGBRC(41) = ALGBRC(40) - (ALGBRC(16)+ 4.00000*ALGBRC(18))
      RATES(2) = ALGBRC(41)
      ALGBRC(5) = log(CONSTS(5)/STATES(3))
      ALGBRC(21) =  STATES(8) ** 3.00000*STATES(9)
      ALGBRC(22) = ( CONSTS(19)*ALGBRC(21)* ( STATES(3)*CONSTS(5)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(5))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(23) = ( CONSTS(28)* ( STATES(3)*CONSTS(5)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(5))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(24) = ( CONSTS(29)*STATES(10)* ( STATES(3)*CONSTS(5)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(5))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(43) = 1.00000/(1.00000+ (CONSTS(44)/STATES(3))*(1.00000+STATES(4)/CONSTS(46)))
      ALGBRC(42) =  CONSTS(5)*EXP( - 0.820000*ALGBRC(1))
      ALGBRC(46) = 1.00000/(1.00000+ (CONSTS(45)/ALGBRC(42))*(1.00000+CONSTS(6)/CONSTS(47)))
      ALGBRC(47) =  CONSTS(48)*( CONSTS(55)*ALGBRC(43)*STATES(13) -  CONSTS(41)*ALGBRC(46)*(1.00000 - STATES(13)))*1.00000
      ALGBRC(49) =  (- 1.00000/( CONSTS(3)*CONSTS(50)))*( 3.00000*ALGBRC(47)+ 3.00000*ALGBRC(39)+ALGBRC(22)+ALGBRC(23)+ALGBRC(24))
      RATES(3) = ALGBRC(49)
      ALGBRC(44) = 1.00000/(1.00000+ (CONSTS(46)/STATES(4))*(1.00000+STATES(3)/CONSTS(44)))
      ALGBRC(45) =  CONSTS(55)*ALGBRC(43)+ CONSTS(42)*ALGBRC(44)
      ALGBRC(48) = 1.00000/(1.00000+ (CONSTS(47)/CONSTS(6))*(1.00000+ALGBRC(42)/CONSTS(45)))
      ALGBRC(50) =  CONSTS(41)*ALGBRC(46)+ CONSTS(43)*ALGBRC(48)
      RATES(13) =  ALGBRC(50)*(1.00000 - STATES(13)) -  ALGBRC(45)*STATES(13)
      ALGBRC(7) = log(CONSTS(6)/STATES(4))
      ALGBRC(25) = STATES(2) ** 4.20000/(0.000350000 ** 4.20000+STATES(2) ** 4.20000)
      ALGBRC(26) = ( CONSTS(31)*ALGBRC(25)* ( STATES(4)*CONSTS(6)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(7))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(27) = STATES(11) ** 3.00000
      ALGBRC(28) =  CONSTS(32)*ALGBRC(27)*(STATES(1) -  ALGBRC(7)*CONSTS(49))
      ALGBRC(29) = 1.00000/(1.00000+EXP((STATES(1)+85.0000)/12.1000))
      ALGBRC(30) =  CONSTS(33)*ALGBRC(29)*(STATES(1) -  ALGBRC(7)*CONSTS(49))
      ALGBRC(31) = ALGBRC(26)+ALGBRC(28)+ALGBRC(30)
      ALGBRC(51) =  (- 1.00000/( CONSTS(3)*CONSTS(50)))*(ALGBRC(31) -  2.00000*ALGBRC(47))
      RATES(1) =  (( CONSTS(3)*CONSTS(50))/( CONSTS(9)*CONSTS(56)))*(ALGBRC(49)+ALGBRC(51)+ 2.00000*ALGBRC(41))
      RATES(4) = ALGBRC(51)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = STATES(1)/CONSTS(49)
      ALGBRC(2) = 0.00600000/(1.00000+EXP((STATES(1)+87.7000)/6.45000))
      ALGBRC(6) = 0.0268000/(1.00000+EXP((STATES(1)+94.2000)/13.3000))
      ALGBRC(8) =  ALGBRC(2)*CONSTS(51)+ ALGBRC(6)*(1.00000 - CONSTS(51))
      ALGBRC(4) = 0.0800000/(1.00000+EXP(- (STATES(1)+51.7000)/7.00000))
      ALGBRC(10) = 0.0800000/(1.00000+EXP(- (STATES(1)+35.5000)/7.00000))
      ALGBRC(12) =  ALGBRC(4)*CONSTS(52)+ ALGBRC(10)*(1.00000 - CONSTS(52))
      ALGBRC(15) = CONSTS(12) - STATES(5)
      ALGBRC(16) =  CONSTS(14)*STATES(5)*STATES(2) -  CONSTS(15)*ALGBRC(15)
      ALGBRC(17) = CONSTS(13) - STATES(6)
      ALGBRC(9) =  12000.0*STATES(2) ** 2.00000
      ALGBRC(11) =  3.70000e+06*STATES(2) ** 2.00000
      ALGBRC(13) =  ALGBRC(9)*ALGBRC(11)*(1.00000/(ALGBRC(9)+CONSTS(17))+1.00000/(CONSTS(16)+CONSTS(17)))
      ALGBRC(14) =  CONSTS(16)*CONSTS(17)*(1.00000/(ALGBRC(9)+CONSTS(17))+1.00000/(CONSTS(16)+CONSTS(17)))
      ALGBRC(18) =  ALGBRC(13)*STATES(6) -  ALGBRC(14)*ALGBRC(17)
      ALGBRC(33) =  (173.600/(1.00000+ALGBRC(17)/5.00000e-05)+6.40000)*1.00000e-05
      ALGBRC(34) = 1.00000/(1.00000+ALGBRC(33)/STATES(2))
      ALGBRC(36) = 1.00000 - ALGBRC(34)
      ALGBRC(38) =  CONSTS(53)*ALGBRC(34)+ CONSTS(35)*ALGBRC(36)
      ALGBRC(3) =  0.500000*log(CONSTS(4)/STATES(2))
      ALGBRC(19) = 0.000450000/(0.000450000+STATES(2))
      ALGBRC(20) = ( CONSTS(18)*STATES(7)*ALGBRC(19)* ( STATES(2)*CONSTS(4)) ** (1.0 / 2)* sinh(ALGBRC(1) - ALGBRC(3)))/( sinh(ALGBRC(1))/ALGBRC(1))
      ALGBRC(32) =  CONSTS(38)*(( 10.5600*ALGBRC(17))/(ALGBRC(17)+5.00000e-05)+1.20000)
      ALGBRC(35) =  ALGBRC(32)*( CONSTS(53)*ALGBRC(34)*STATES(12) -  CONSTS(34)*CONSTS(54)*(1.00000 - STATES(12)))*1.00000
      ALGBRC(37) =  (1.00000+ 0.00100000*( STATES(3) ** 3.00000*CONSTS(4)+ CONSTS(5) ** 3.00000*STATES(2)))*(1.00000+STATES(2)/0.00690000)
      ALGBRC(39) = ( CONSTS(40)*( STATES(3) ** 3.00000*CONSTS(4)*EXP( CONSTS(39)*ALGBRC(1)) -  CONSTS(5) ** 3.00000*STATES(2)*EXP( (CONSTS(39) - 1.00000)*ALGBRC(1))))/ALGBRC(37)
      ALGBRC(40) =  (- 1.00000/( 2.00000*CONSTS(3)*CONSTS(50)))*((ALGBRC(20)+ 2.00000*ALGBRC(35)) -  2.00000*ALGBRC(39))
      ALGBRC(41) = ALGBRC(40) - (ALGBRC(16)+ 4.00000*ALGBRC(18))
      ALGBRC(5) = log(CONSTS(5)/STATES(3))
      ALGBRC(21) =  STATES(8) ** 3.00000*STATES(9)
      ALGBRC(22) = ( CONSTS(19)*ALGBRC(21)* ( STATES(3)*CONSTS(5)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(5))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(23) = ( CONSTS(28)* ( STATES(3)*CONSTS(5)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(5))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(24) = ( CONSTS(29)*STATES(10)* ( STATES(3)*CONSTS(5)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(5))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(43) = 1.00000/(1.00000+ (CONSTS(44)/STATES(3))*(1.00000+STATES(4)/CONSTS(46)))
      ALGBRC(42) =  CONSTS(5)*EXP( - 0.820000*ALGBRC(1))
      ALGBRC(46) = 1.00000/(1.00000+ (CONSTS(45)/ALGBRC(42))*(1.00000+CONSTS(6)/CONSTS(47)))
      ALGBRC(47) =  CONSTS(48)*( CONSTS(55)*ALGBRC(43)*STATES(13) -  CONSTS(41)*ALGBRC(46)*(1.00000 - STATES(13)))*1.00000
      ALGBRC(49) =  (- 1.00000/( CONSTS(3)*CONSTS(50)))*( 3.00000*ALGBRC(47)+ 3.00000*ALGBRC(39)+ALGBRC(22)+ALGBRC(23)+ALGBRC(24))
      ALGBRC(44) = 1.00000/(1.00000+ (CONSTS(46)/STATES(4))*(1.00000+STATES(3)/CONSTS(44)))
      ALGBRC(45) =  CONSTS(55)*ALGBRC(43)+ CONSTS(42)*ALGBRC(44)
      ALGBRC(48) = 1.00000/(1.00000+ (CONSTS(47)/CONSTS(6))*(1.00000+ALGBRC(42)/CONSTS(45)))
      ALGBRC(50) =  CONSTS(41)*ALGBRC(46)+ CONSTS(43)*ALGBRC(48)
      ALGBRC(7) = log(CONSTS(6)/STATES(4))
      ALGBRC(25) = STATES(2) ** 4.20000/(0.000350000 ** 4.20000+STATES(2) ** 4.20000)
      ALGBRC(26) = ( CONSTS(31)*ALGBRC(25)* ( STATES(4)*CONSTS(6)) ** (1.0 / 2)* sinh( 0.500000*(ALGBRC(1) - ALGBRC(7))))/( sinh( 0.500000*ALGBRC(1))/( 0.500000*ALGBRC(1)))
      ALGBRC(27) = STATES(11) ** 3.00000
      ALGBRC(28) =  CONSTS(32)*ALGBRC(27)*(STATES(1) -  ALGBRC(7)*CONSTS(49))
      ALGBRC(29) = 1.00000/(1.00000+EXP((STATES(1)+85.0000)/12.1000))
      ALGBRC(30) =  CONSTS(33)*ALGBRC(29)*(STATES(1) -  ALGBRC(7)*CONSTS(49))
      ALGBRC(31) = ALGBRC(26)+ALGBRC(28)+ALGBRC(30)
      ALGBRC(51) =  (- 1.00000/( CONSTS(3)*CONSTS(50)))*(ALGBRC(31) -  2.00000*ALGBRC(47))
      RETURN
      END