C
C There are a total of 69 entries in the algebraic variable array.
C There are a total of 22 entries in each of the rate and state variable arrays.
C There are a total of 63 entries in the constant variable array.
C
C
C VOI is time in component Time (time_units).
C CONSTS(1) is T in component Environment (Temperature_units).
C CONSTS(2) is F in component Environment (F_units).
C CONSTS(3) is R in component Environment (R_units).
C CONSTS(4) is Ca_o in component Environment (millimolar).
C CONSTS(5) is Na_o in component Environment (millimolar).
C CONSTS(6) is K_o in component Environment (millimolar).
C CONSTS(58) is FoRT in component Environment (inverse_voltage_units).
C CONSTS(59) is RToF in component Environment (voltage_units).
C CONSTS(7) is Cm in component membrane (capacitance_units).
C STATES(1) is Vm in component membrane (voltage_units).
C STATES(2) is Ca_i_peripheral in component membrane (millimolar).
C STATES(3) is Ca_i_bulk in component membrane (millimolar).
C STATES(4) is Ca_sr in component membrane (millimolar).
C STATES(5) is Na_i in component membrane (millimolar).
C STATES(6) is K_i in component membrane (millimolar).
C ALGBRC(44) is i_k1 in component i_k1 (current_units).
C ALGBRC(34) is i_to_fast in component i_to_fast (current_units).
C ALGBRC(37) is i_to_sustained in component i_to_sustained (current_units).
C ALGBRC(47) is i_kr in component i_kr (current_units).
C ALGBRC(49) is i_ks in component i_ks (current_units).
C ALGBRC(55) is i_kb in component i_kb (current_units).
C ALGBRC(53) is i_nak in component i_nak (current_units).
C ALGBRC(32) is i_cal in component i_cal (current_units).
C ALGBRC(26) is i_cat in component i_cat (current_units).
C ALGBRC(39) is i_na_fast in component i_na_fast (current_units).
C ALGBRC(41) is i_na_late in component i_na_late (current_units).
C ALGBRC(60) is i_pmca in component i_pmca (current_units).
C ALGBRC(59) is i_cab in component i_cab (current_units).
C ALGBRC(67) is i_f_na in component i_f_na (current_units).
C ALGBRC(63) is i_f_k in component i_f_k (current_units).
C ALGBRC(50) is i_naca in component i_naca (current_units).
C ALGBRC(57) is i_nab in component i_nab (current_units).
C ALGBRC(1) is i_stim in component membrane (current_units).
C ALGBRC(62) is i_rel_per in component i_rel_per (current_units).
C ALGBRC(64) is i_serca_per in component i_serca_per (current_units).
C ALGBRC(66) is i_serca_bulk in component i_serca_bulk (current_units).
C ALGBRC(69) is i_diff in component i_diff (current_units).
C ALGBRC(68) is i_leak in component i_leak (current_units).
C CONSTS(61) is peripheral_volume in component membrane (volume_units).
C CONSTS(60) is bulk_volume in component membrane (volume_units).
C CONSTS(63) is diffusable_volume in component membrane (volume_units).
C CONSTS(62) is sr_volume in component membrane (volume_units).
C CONSTS(8) is cell_volume in component membrane (volume_units).
C CONSTS(9) is bulk_fraction in component membrane (dimensionless).
C CONSTS(10) is periphery_fraction in component membrane (dimensionless).
C CONSTS(11) is sr_fraction in component membrane (dimensionless).
C CONSTS(12) is stim_start in component membrane (time_units).
C CONSTS(13) is stim_end in component membrane (time_units).
C CONSTS(14) is stim_period in component membrane (time_units).
C CONSTS(15) is stim_duration in component membrane (time_units).
C CONSTS(16) is stim_amplitude in component membrane (current_units).
C ALGBRC(2) is x_inf_Ttype in component x_Ttype (dimensionless).
C CONSTS(17) is tau_x_Ttype in component x_Ttype (time_units).
C STATES(7) is x_Ttype in component x_Ttype (dimensionless).
C ALGBRC(3) is y_inf_Ttype in component y_Ttype (dimensionless).
C ALGBRC(18) is tau_y_Ttype in component y_Ttype (time_units).
C STATES(8) is y_Ttype in component y_Ttype (dimensionless).
C ALGBRC(17) is E_Ca in component i_cat (voltage_units).
C CONSTS(18) is G_max_Ttype in component i_cat (conductance_units).
C ALGBRC(4) is x_inf_Ltype in component x_Ltype (dimensionless).
C CONSTS(19) is tau_x_Ltype in component x_Ltype (time_units).
C STATES(9) is x_Ltype in component x_Ltype (dimensionless).
C ALGBRC(5) is y_inf_Ltype in component y_Ltype (dimensionless).
C ALGBRC(19) is tau_y_Ltype in component y_Ltype (time_units).
C STATES(10) is y_Ltype in component y_Ltype (dimensionless).
C ALGBRC(6) is y_ca_inf_Ltype in component y_ca_Ltype (dimensionless).
C ALGBRC(20) is tau_y_ca_Ltype in component y_ca_Ltype (time_units).
C STATES(11) is y_ca_Ltype in component y_ca_Ltype (dimensionless).
C ALGBRC(29) is E_Ca in component i_cal (voltage_units).
C CONSTS(20) is G_max_Ltype in component i_cal (conductance_units).
C ALGBRC(7) is x_inf_to_fast in component x_to_fast (dimensionless).
C CONSTS(21) is tau_x_to_fast in component x_to_fast (time_units).
C STATES(12) is x_to_fast in component x_to_fast (dimensionless).
C ALGBRC(8) is y_inf_to_fast in component y_to_fast (dimensionless).
C CONSTS(22) is tau_y_to_fast in component y_to_fast (time_units).
C STATES(13) is y_to_fast in component y_to_fast (dimensionless).
C ALGBRC(33) is E_k in component i_to_fast (voltage_units).
C CONSTS(23) is G_max_to_fast in component i_to_fast (conductance_units).
C ALGBRC(35) is x_to_sustained in component x_to_sustained (dimensionless).
C ALGBRC(36) is E_k in component i_to_sustained (voltage_units).
C CONSTS(24) is G_max_to_sustained in component i_to_sustained (conductance_units).
C ALGBRC(9) is x_inf_na_fast in component x_na_fast (dimensionless).
C CONSTS(25) is tau_x_na_fast in component x_na_fast (time_units).
C STATES(14) is x_na_fast in component x_na_fast (dimensionless).
C ALGBRC(10) is y_inf_na_fast in component y_na_fast (dimensionless).
C CONSTS(26) is tau_y_na_fast in component y_na_fast (time_units).
C STATES(15) is y_na_fast in component y_na_fast (dimensionless).
C ALGBRC(38) is E_na in component i_na_fast (voltage_units).
C CONSTS(27) is G_max_na_fast in component i_na_fast (conductance_units).
C ALGBRC(11) is x_inf_na_late in component x_na_late (dimensionless).
C CONSTS(28) is tau_x_na_late in component x_na_late (time_units).
C STATES(16) is x_na_late in component x_na_late (dimensionless).
C ALGBRC(12) is y_inf_na_late in component y_na_late (dimensionless).
C ALGBRC(21) is tau_y_na_late in component y_na_late (time_units).
C STATES(17) is y_na_late in component y_na_late (dimensionless).
C ALGBRC(40) is E_na in component i_na_late (voltage_units).
C CONSTS(29) is G_max_na_late in component i_na_late (conductance_units).
C ALGBRC(42) is x_k1 in component x_k1 (dimensionless).
C ALGBRC(43) is E_k in component i_k1 (voltage_units).
C CONSTS(30) is G_max_k1 in component i_k1 (conductance_units).
C ALGBRC(45) is x_kr in component x_kr (dimensionless).
C ALGBRC(13) is y_inf_kr in component y_kr (dimensionless).
C ALGBRC(30) is tau_y_kr in component y_kr (time_units).
C STATES(18) is y_kr in component y_kr (dimensionless).
C ALGBRC(22) is ykrv1 in component y_kr (rate_constants_units).
C ALGBRC(27) is ykrv2 in component y_kr (rate_constants_units).
C ALGBRC(46) is E_k in component i_kr (voltage_units).
C CONSTS(31) is G_max_kr in component i_kr (conductance_units).
C ALGBRC(14) is x_inf_ks in component x_ks (dimensionless).
C ALGBRC(23) is tau_x_ks in component x_ks (time_units).
C STATES(19) is x_ks in component x_ks (dimensionless).
C ALGBRC(28) is y_inf_ks in component y_ks (dimensionless).
C ALGBRC(31) is tau_y_ks in component y_ks (time_units).
C STATES(20) is y_ks in component y_ks (dimensionless).
C ALGBRC(48) is E_k in component i_ks (voltage_units).
C CONSTS(32) is G_max_ks in component i_ks (conductance_units).
C CONSTS(33) is n_NaCa in component i_naca (dimensionless).
C CONSTS(34) is g_NaCa in component i_naca (current_units).
C CONSTS(35) is d_NaCa in component i_naca (dimensionless).
C CONSTS(36) is gamma in component i_naca (dimensionless).
C ALGBRC(51) is x_nak in component x_nak (dimensionless).
C ALGBRC(52) is y_nak in component y_nak (dimensionless).
C CONSTS(37) is g_nak in component i_nak (current_units).
C ALGBRC(54) is E_k in component i_kb (voltage_units).
C CONSTS(38) is G_max_kb in component i_kb (conductance_units).
C ALGBRC(56) is E_na in component i_nab (voltage_units).
C CONSTS(39) is G_max_nab in component i_nab (conductance_units).
C ALGBRC(58) is E_ca in component i_cab (voltage_units).
C CONSTS(40) is G_max_cab in component i_cab (conductance_units).
C CONSTS(41) is PMCA_max in component i_pmca (current_units).
C CONSTS(42) is Kpmca in component i_pmca (millimolar).
C CONSTS(43) is Hpmca in component i_pmca (dimensionless).
C CONSTS(44) is G_f_k in component i_f_k (conductance_units).
C ALGBRC(61) is E_k in component i_f_k (voltage_units).
C STATES(21) is y_gate_f_k in component y_gate_f_k (dimensionless).
C CONSTS(45) is G_f_na in component i_f_na (conductance_units).
C ALGBRC(65) is E_na in component i_f_na (voltage_units).
C STATES(22) is y_gate_f_na in component y_gate_f_na (dimensionless).
C ALGBRC(15) is y_inf_f_gate in component y_gate_f_k (dimensionless).
C ALGBRC(24) is tau_y_f_gate in component y_gate_f_k (time_units).
C ALGBRC(16) is y_inf_f_gate in component y_gate_f_na (dimensionless).
C ALGBRC(25) is tau_y_f_gate in component y_gate_f_na (time_units).
C CONSTS(46) is REL_max in component i_rel_per (current_per_millimolar_units).
C CONSTS(47) is Krel in component i_rel_per (millimolar).
C CONSTS(48) is SERCA_max in component i_serca_per (current_units).
C CONSTS(49) is Kmf in component i_serca_per (millimolar).
C CONSTS(50) is Kmr in component i_serca_per (millimolar).
C CONSTS(51) is H in component i_serca_per (dimensionless).
C CONSTS(52) is SERCA_max in component i_serca_bulk (current_units).
C CONSTS(53) is Kmf in component i_serca_bulk (millimolar).
C CONSTS(54) is Kmr in component i_serca_bulk (millimolar).
C CONSTS(55) is H in component i_serca_bulk (dimensionless).
C CONSTS(56) is LEAK_max in component i_leak (current_per_millimolar_units).
C CONSTS(57) is DIFF_max in component i_diff (current_per_millimolar_units).
C RATES(1) is d/dt Vm in component membrane (voltage_units).
C RATES(2) is d/dt Ca_i_peripheral in component membrane (millimolar).
C RATES(3) is d/dt Ca_i_bulk in component membrane (millimolar).
C RATES(4) is d/dt Ca_sr in component membrane (millimolar).
C RATES(6) is d/dt K_i in component membrane (millimolar).
C RATES(5) is d/dt Na_i in component membrane (millimolar).
C RATES(7) is d/dt x_Ttype in component x_Ttype (dimensionless).
C RATES(8) is d/dt y_Ttype in component y_Ttype (dimensionless).
C RATES(9) is d/dt x_Ltype in component x_Ltype (dimensionless).
C RATES(10) is d/dt y_Ltype in component y_Ltype (dimensionless).
C RATES(11) is d/dt y_ca_Ltype in component y_ca_Ltype (dimensionless).
C RATES(12) is d/dt x_to_fast in component x_to_fast (dimensionless).
C RATES(13) is d/dt y_to_fast in component y_to_fast (dimensionless).
C RATES(14) is d/dt x_na_fast in component x_na_fast (dimensionless).
C RATES(15) is d/dt y_na_fast in component y_na_fast (dimensionless).
C RATES(16) is d/dt x_na_late in component x_na_late (dimensionless).
C RATES(17) is d/dt y_na_late in component y_na_late (dimensionless).
C RATES(18) is d/dt y_kr in component y_kr (dimensionless).
C RATES(19) is d/dt x_ks in component x_ks (dimensionless).
C RATES(20) is d/dt y_ks in component y_ks (dimensionless).
C RATES(21) is d/dt y_gate_f_k in component y_gate_f_k (dimensionless).
C RATES(22) is d/dt y_gate_f_na in component y_gate_f_na (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 310.0
      CONSTS(2) = 96485.3415
      CONSTS(3) = 8314.472
      CONSTS(4) = 2.0
      CONSTS(5) = 140
      CONSTS(6) = 5.4
      CONSTS(7) = 69.0
      STATES(1) = -88.34
      STATES(2) = 0.00001
      STATES(3) = 0.000032
      STATES(4) = 0.17
      STATES(5) = 6.7
      STATES(6) = 140
      CONSTS(8) = 13266.5
      CONSTS(9) = 0.6
      CONSTS(10) = 0.2
      CONSTS(11) = 0.06
      CONSTS(12) = 700
      CONSTS(13) = 9000000
      CONSTS(14) = 500
      CONSTS(15) = 0.5
      CONSTS(16) = -4320.0
      CONSTS(17) = 1.0
      STATES(7) = 0.001337
      STATES(8) = 0.01
      CONSTS(18) = 0.9
      CONSTS(19) = 0.7
      STATES(9) = 0.000003
      STATES(10) = 0.1
      STATES(11) = 0.7
      CONSTS(20) = 5.4
      CONSTS(21) = 5.0
      STATES(12) = 0.0
      CONSTS(22) = 350.0
      STATES(13) = 0.7
      CONSTS(23) = 10.0
      CONSTS(24) = 3.0
      CONSTS(25) = 0.005
      STATES(14) = 0.000007
      CONSTS(26) = 2.0
      STATES(15) = 0.978861
      CONSTS(27) = 1140.0
      CONSTS(28) = 15.0
      STATES(16) = 0.000012
      STATES(17) = 0.864489
      CONSTS(29) = 2.0
      CONSTS(30) = 20.0
      STATES(18) = 0.25
      CONSTS(31) = 1.5
      STATES(19) = 1.0
      STATES(20) = 0.0
      CONSTS(32) = 3.0
      CONSTS(33) = 3
      CONSTS(34) = 0.001
      CONSTS(35) = 0.001
      CONSTS(36) = 0.5
      CONSTS(37) = 442.2
      CONSTS(38) = 0.01
      CONSTS(39) = 0.01
      CONSTS(40) = 0.0001
      CONSTS(41) = 5.0
      CONSTS(42) = 0.0001
      CONSTS(43) = 1.5
      CONSTS(44) = 0.188709677
      STATES(21) = 0.011099
      CONSTS(45) = 0.045290323
      STATES(22) = 0.011099
      CONSTS(46) = 2500.0
      CONSTS(47) = 0.001
      CONSTS(48) = 120.0
      CONSTS(49) = 0.000246
      CONSTS(50) = 1.7
      CONSTS(51) = 1.6
      CONSTS(52) = 120.0
      CONSTS(53) = 0.000246
      CONSTS(54) = 1.7
      CONSTS(55) = 1.6
      CONSTS(56) = 10.0
      CONSTS(57) = 5000.0
      CONSTS(58) = CONSTS(2)/( CONSTS(3)*CONSTS(1))
      CONSTS(59) = ( CONSTS(3)*CONSTS(1))/CONSTS(2)
      CONSTS(60) =  CONSTS(9)*CONSTS(8)
      CONSTS(61) =  CONSTS(10)*CONSTS(8)
      CONSTS(62) =  CONSTS(11)*CONSTS(8)
      CONSTS(63) =  (CONSTS(9)+CONSTS(10))*CONSTS(8)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+47.8000)/-5.50000))
      RATES(7) = (ALGBRC(2) - STATES(7))/CONSTS(17)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+14.6000)/-5.50000))
      RATES(9) = (ALGBRC(4) - STATES(9))/CONSTS(19)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+-7.00000)/-9.00000))
      RATES(12) = (ALGBRC(7) - STATES(12))/CONSTS(21)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+27.5000)/8.00000))
      RATES(13) = (ALGBRC(8) - STATES(13))/CONSTS(22)
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+25.0000)/-5.00000))
      RATES(14) = (ALGBRC(9) - STATES(14))/CONSTS(25)
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+69.0000)/3.96000))
      RATES(15) = (ALGBRC(10) - STATES(15))/CONSTS(26)
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/-5.00000))
      RATES(16) = (ALGBRC(11) - STATES(16))/CONSTS(28)
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+67.9000)/3.87000))
      ALGBRC(18) =  1.42271*EXP( -0.0511900*STATES(1))
      RATES(8) = (ALGBRC(3) - STATES(8))/ALGBRC(18)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+31.0000)/5.54000))
      ALGBRC(19) = 25.1000/(0.0400000+ 0.700000*EXP( -1.00000* 0.0280000*(STATES(1)+14.5000) ** 2.00000))
      RATES(10) = (ALGBRC(5) - STATES(10))/ALGBRC(19)
      ALGBRC(6) = 0.400000+0.600000/(1.00000+STATES(2)/0.000100000 ** 2.00000)
      ALGBRC(20) = 2.00000+80.0000/(1.00000+STATES(2)/0.000100000 ** 2.00000)
      RATES(11) = (ALGBRC(6) - STATES(11))/ALGBRC(20)
      ALGBRC(12) = 0.100000+0.900000/(1.00000+EXP((STATES(1)+75.6000)/6.30000))
      ALGBRC(21) = 120.000+ 1.00000*EXP((STATES(1)+100.000)/25.0000)
      RATES(17) = (ALGBRC(12) - STATES(17))/ALGBRC(21)
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1) - 1.50000)/-16.7000))
      ALGBRC(23) = TERNRY(ABS(STATES(1)+30.0000).LT.0.0145000, 417.946, 1.00000/(( 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(19) = (ALGBRC(14) - STATES(19))/ALGBRC(23)
      ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(1)+109.000)/10.0000))
      ALGBRC(24) = 6000.00/(EXP( -1.00000*(2.90000+ 0.0400000*STATES(1)))+EXP( 1.00000*(3.60000+ 0.110000*STATES(1))))
      RATES(21) = (ALGBRC(15) - STATES(21))/ALGBRC(24)
      ALGBRC(16) = 1.00000/(1.00000+EXP((STATES(1)+109.000)/10.0000))
      ALGBRC(25) = 6000.00/(EXP( -1.00000*(2.90000+ 0.0400000*STATES(1)))+EXP( 1.00000*(3.60000+ 0.110000*STATES(1))))
      RATES(22) = (ALGBRC(16) - STATES(22))/ALGBRC(25)
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+50.0000)/-7.50000))
      ALGBRC(22) = TERNRY(ABS(STATES(1)+7.00000).GT.0.00100000, ( 0.00138000*1.00000*(STATES(1)+7.00000))/(1.00000 - EXP( -0.123000*(STATES(1)+7.00000))), 0.00138000/0.123000)
      ALGBRC(27) = TERNRY(ABS(STATES(1)+10.0000).GT.0.00100000, ( 6.10000e-05*1.00000*(STATES(1)+10.0000))/(EXP( 0.145000*(STATES(1)+10.0000)) - 1.00000), 0.000610000/0.145000)
      ALGBRC(30) = 1.00000/(ALGBRC(22)+ALGBRC(27))
      RATES(18) = (ALGBRC(13) - STATES(18))/ALGBRC(30)
      ALGBRC(28) = ALGBRC(14)
      ALGBRC(31) =  4.00000*ALGBRC(23)
      RATES(20) = (ALGBRC(28) - STATES(20))/ALGBRC(31)
      ALGBRC(42) = 1.00000/(1.00000+EXP((92.0000+STATES(1))/10.0000))
      ALGBRC(43) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(44) =  CONSTS(30)*CONSTS(6)/5.40000 ** 0.800000*ALGBRC(42)*(STATES(1) - ALGBRC(43))
      ALGBRC(33) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(34) =  CONSTS(23)*STATES(12)*STATES(13)*(STATES(1) - ALGBRC(33))
      ALGBRC(35) = 1.00000/(1.00000+EXP((5.00000 - STATES(1))/17.0000))
      ALGBRC(36) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(37) =  CONSTS(24)*ALGBRC(35)*(STATES(1) - ALGBRC(36))
      ALGBRC(45) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
      ALGBRC(46) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(47) =  CONSTS(31)*CONSTS(6)/5.40000 ** 1.00000*ALGBRC(45)*STATES(18)*(STATES(1) - ALGBRC(46))
      ALGBRC(48) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(49) =  CONSTS(32)*STATES(19)*STATES(20)*(STATES(1) - ALGBRC(48))
      ALGBRC(54) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(55) =  CONSTS(38)*(STATES(1) - ALGBRC(54))
      ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(1)+80.0000)/-45.0000))
      ALGBRC(52) = 1.00000/(1.00000+EXP((STATES(1)+0.00000)/125.000))
      ALGBRC(53) =  CONSTS(37)*ALGBRC(51)*ALGBRC(52)*(1.00000/(1.00000+1.90000/CONSTS(6) ** 1.45000))*(1.00000/(1.00000+31.9800/STATES(5) ** 1.00000))
      ALGBRC(61) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(63) =  CONSTS(44)*STATES(21)*(STATES(1) - ALGBRC(61))
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(12).AND.VOI.LE.CONSTS(13).AND.(VOI - CONSTS(12)) -  INT((VOI - CONSTS(12))/CONSTS(14))*CONSTS(14).LE.CONSTS(15), CONSTS(16), 0.00000)
      RATES(6) = ( ( -1.00000*ALGBRC(34)+ -1.00000*ALGBRC(37)+ -1.00000*ALGBRC(47)+ -1.00000*ALGBRC(49)+ -1.00000*ALGBRC(44)+ -1.00000*ALGBRC(55)+ -1.00000*ALGBRC(63)+ -1.00000*ALGBRC(1)+ 2.00000*ALGBRC(53))*1000.00)/( CONSTS(2)*CONSTS(63))
      ALGBRC(29) =  0.500000*CONSTS(59)*log(CONSTS(4)/STATES(2))
      ALGBRC(32) =  CONSTS(20)*STATES(9)*STATES(10)*STATES(11)*(STATES(1) - ALGBRC(29))
      ALGBRC(17) =  0.500000*CONSTS(59)*log(CONSTS(4)/STATES(2))
      ALGBRC(26) =  CONSTS(18)*STATES(7)*STATES(8)*(STATES(1) - ALGBRC(17))
      ALGBRC(38) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(39) =  CONSTS(27)*STATES(14)*STATES(15)*(STATES(1) - ALGBRC(38))
      ALGBRC(40) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(41) =  CONSTS(29)*STATES(16)*STATES(17)*(STATES(1) - ALGBRC(40))
      ALGBRC(60) =  CONSTS(41)*(1.00000/(1.00000+CONSTS(42)/STATES(2) ** CONSTS(43)))
      ALGBRC(58) =  0.500000*CONSTS(59)*log(CONSTS(4)/STATES(2))
      ALGBRC(59) =  CONSTS(40)*(STATES(1) - ALGBRC(58))
      ALGBRC(65) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(67) =  CONSTS(45)*STATES(22)*(STATES(1) - ALGBRC(65))
      ALGBRC(50) = ( 512.009*CONSTS(34)*( EXP(( CONSTS(36)*(CONSTS(33) - 2.00000)*STATES(1))/CONSTS(59))*STATES(5) ** CONSTS(33)*CONSTS(4) -  STATES(2)*EXP(( (CONSTS(36) - 1.00000)*(CONSTS(33) - 2.00000)*STATES(1))/CONSTS(59))*CONSTS(5) ** CONSTS(33)))/( (1.00000+ CONSTS(35)*( STATES(2)*CONSTS(5) ** CONSTS(33)+ CONSTS(4)*STATES(5) ** CONSTS(33)))*(1.00000+STATES(2)/0.00690000))
      ALGBRC(56) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(57) =  CONSTS(39)*(STATES(1) - ALGBRC(56))
      RATES(1) =  (( - 1.00000*1.00000)/CONSTS(7))*(ALGBRC(44)+ALGBRC(34)+ALGBRC(37)+ALGBRC(47)+ALGBRC(49)+ALGBRC(55)+ALGBRC(53)+ALGBRC(32)+ALGBRC(26)+ALGBRC(39)+ALGBRC(41)+ALGBRC(60)+ALGBRC(59)+ALGBRC(67)+ALGBRC(63)+ALGBRC(50)+ALGBRC(57)+ALGBRC(1))
      ALGBRC(62) =  CONSTS(46)*(STATES(4)/(1.00000+CONSTS(47)/STATES(2) ** 2.00000))
      ALGBRC(64) =  CONSTS(48)*((STATES(2)/CONSTS(49) ** CONSTS(51) - STATES(4)/CONSTS(50) ** CONSTS(51))/(1.00000+STATES(2)/CONSTS(49) ** CONSTS(51)+STATES(4)/CONSTS(50) ** CONSTS(51)))
      ALGBRC(66) =  CONSTS(52)*((STATES(3)/CONSTS(53) ** CONSTS(55) - STATES(4)/CONSTS(54) ** CONSTS(55))/(1.00000+STATES(3)/CONSTS(53) ** CONSTS(55)+STATES(4)/CONSTS(54) ** CONSTS(55)))
      ALGBRC(68) =  CONSTS(56)*(STATES(4) - STATES(3))
      RATES(4) = ( ( -1.00000*ALGBRC(62)+ -1.00000*ALGBRC(68)+ALGBRC(66)+ALGBRC(64))*1000.00)/( 2.00000*CONSTS(2)*CONSTS(62))
      RATES(5) = ( ( -1.00000*ALGBRC(39)+ -1.00000*ALGBRC(41)+ -3.00000*ALGBRC(53)+ -3.00000*ALGBRC(50)+ -1.00000*ALGBRC(67)+ -1.00000*ALGBRC(57))*1000.00)/( CONSTS(2)*CONSTS(63))
      ALGBRC(69) =  CONSTS(57)*(STATES(2) - STATES(3))
      RATES(2) = ( ( -1.00000*ALGBRC(32)+ -1.00000*ALGBRC(26)+ -1.00000*ALGBRC(59)+ 1.00000*ALGBRC(62)+ -1.00000*ALGBRC(64)+ -1.00000*ALGBRC(69)+ 2.00000*ALGBRC(50)+ -1.00000*ALGBRC(60))*1000.00)/( 2.00000*CONSTS(2)*CONSTS(61))
      RATES(3) = ( ( -1.00000*ALGBRC(66)+ALGBRC(69)+ALGBRC(68))*1000.00)/( 2.00000*CONSTS(2)*CONSTS(60))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+47.8000)/-5.50000))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+14.6000)/-5.50000))
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+-7.00000)/-9.00000))
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+27.5000)/8.00000))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+25.0000)/-5.00000))
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+69.0000)/3.96000))
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/-5.00000))
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+67.9000)/3.87000))
      ALGBRC(18) =  1.42271*EXP( -0.0511900*STATES(1))
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+31.0000)/5.54000))
      ALGBRC(19) = 25.1000/(0.0400000+ 0.700000*EXP( -1.00000* 0.0280000*(STATES(1)+14.5000) ** 2.00000))
      ALGBRC(6) = 0.400000+0.600000/(1.00000+STATES(2)/0.000100000 ** 2.00000)
      ALGBRC(20) = 2.00000+80.0000/(1.00000+STATES(2)/0.000100000 ** 2.00000)
      ALGBRC(12) = 0.100000+0.900000/(1.00000+EXP((STATES(1)+75.6000)/6.30000))
      ALGBRC(21) = 120.000+ 1.00000*EXP((STATES(1)+100.000)/25.0000)
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1) - 1.50000)/-16.7000))
      ALGBRC(23) = TERNRY(ABS(STATES(1)+30.0000).LT.0.0145000, 417.946, 1.00000/(( 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(15) = 1.00000/(1.00000+EXP((STATES(1)+109.000)/10.0000))
      ALGBRC(24) = 6000.00/(EXP( -1.00000*(2.90000+ 0.0400000*STATES(1)))+EXP( 1.00000*(3.60000+ 0.110000*STATES(1))))
      ALGBRC(16) = 1.00000/(1.00000+EXP((STATES(1)+109.000)/10.0000))
      ALGBRC(25) = 6000.00/(EXP( -1.00000*(2.90000+ 0.0400000*STATES(1)))+EXP( 1.00000*(3.60000+ 0.110000*STATES(1))))
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+50.0000)/-7.50000))
      ALGBRC(22) = TERNRY(ABS(STATES(1)+7.00000).GT.0.00100000, ( 0.00138000*1.00000*(STATES(1)+7.00000))/(1.00000 - EXP( -0.123000*(STATES(1)+7.00000))), 0.00138000/0.123000)
      ALGBRC(27) = TERNRY(ABS(STATES(1)+10.0000).GT.0.00100000, ( 6.10000e-05*1.00000*(STATES(1)+10.0000))/(EXP( 0.145000*(STATES(1)+10.0000)) - 1.00000), 0.000610000/0.145000)
      ALGBRC(30) = 1.00000/(ALGBRC(22)+ALGBRC(27))
      ALGBRC(28) = ALGBRC(14)
      ALGBRC(31) =  4.00000*ALGBRC(23)
      ALGBRC(42) = 1.00000/(1.00000+EXP((92.0000+STATES(1))/10.0000))
      ALGBRC(43) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(44) =  CONSTS(30)*CONSTS(6)/5.40000 ** 0.800000*ALGBRC(42)*(STATES(1) - ALGBRC(43))
      ALGBRC(33) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(34) =  CONSTS(23)*STATES(12)*STATES(13)*(STATES(1) - ALGBRC(33))
      ALGBRC(35) = 1.00000/(1.00000+EXP((5.00000 - STATES(1))/17.0000))
      ALGBRC(36) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(37) =  CONSTS(24)*ALGBRC(35)*(STATES(1) - ALGBRC(36))
      ALGBRC(45) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
      ALGBRC(46) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(47) =  CONSTS(31)*CONSTS(6)/5.40000 ** 1.00000*ALGBRC(45)*STATES(18)*(STATES(1) - ALGBRC(46))
      ALGBRC(48) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(49) =  CONSTS(32)*STATES(19)*STATES(20)*(STATES(1) - ALGBRC(48))
      ALGBRC(54) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(55) =  CONSTS(38)*(STATES(1) - ALGBRC(54))
      ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(1)+80.0000)/-45.0000))
      ALGBRC(52) = 1.00000/(1.00000+EXP((STATES(1)+0.00000)/125.000))
      ALGBRC(53) =  CONSTS(37)*ALGBRC(51)*ALGBRC(52)*(1.00000/(1.00000+1.90000/CONSTS(6) ** 1.45000))*(1.00000/(1.00000+31.9800/STATES(5) ** 1.00000))
      ALGBRC(61) =  CONSTS(59)*log(CONSTS(6)/STATES(6))
      ALGBRC(63) =  CONSTS(44)*STATES(21)*(STATES(1) - ALGBRC(61))
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(12).AND.VOI.LE.CONSTS(13).AND.(VOI - CONSTS(12)) -  INT((VOI - CONSTS(12))/CONSTS(14))*CONSTS(14).LE.CONSTS(15), CONSTS(16), 0.00000)
      ALGBRC(29) =  0.500000*CONSTS(59)*log(CONSTS(4)/STATES(2))
      ALGBRC(32) =  CONSTS(20)*STATES(9)*STATES(10)*STATES(11)*(STATES(1) - ALGBRC(29))
      ALGBRC(17) =  0.500000*CONSTS(59)*log(CONSTS(4)/STATES(2))
      ALGBRC(26) =  CONSTS(18)*STATES(7)*STATES(8)*(STATES(1) - ALGBRC(17))
      ALGBRC(38) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(39) =  CONSTS(27)*STATES(14)*STATES(15)*(STATES(1) - ALGBRC(38))
      ALGBRC(40) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(41) =  CONSTS(29)*STATES(16)*STATES(17)*(STATES(1) - ALGBRC(40))
      ALGBRC(60) =  CONSTS(41)*(1.00000/(1.00000+CONSTS(42)/STATES(2) ** CONSTS(43)))
      ALGBRC(58) =  0.500000*CONSTS(59)*log(CONSTS(4)/STATES(2))
      ALGBRC(59) =  CONSTS(40)*(STATES(1) - ALGBRC(58))
      ALGBRC(65) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(67) =  CONSTS(45)*STATES(22)*(STATES(1) - ALGBRC(65))
      ALGBRC(50) = ( 512.009*CONSTS(34)*( EXP(( CONSTS(36)*(CONSTS(33) - 2.00000)*STATES(1))/CONSTS(59))*STATES(5) ** CONSTS(33)*CONSTS(4) -  STATES(2)*EXP(( (CONSTS(36) - 1.00000)*(CONSTS(33) - 2.00000)*STATES(1))/CONSTS(59))*CONSTS(5) ** CONSTS(33)))/( (1.00000+ CONSTS(35)*( STATES(2)*CONSTS(5) ** CONSTS(33)+ CONSTS(4)*STATES(5) ** CONSTS(33)))*(1.00000+STATES(2)/0.00690000))
      ALGBRC(56) =  CONSTS(59)*log(CONSTS(5)/STATES(5))
      ALGBRC(57) =  CONSTS(39)*(STATES(1) - ALGBRC(56))
      ALGBRC(62) =  CONSTS(46)*(STATES(4)/(1.00000+CONSTS(47)/STATES(2) ** 2.00000))
      ALGBRC(64) =  CONSTS(48)*((STATES(2)/CONSTS(49) ** CONSTS(51) - STATES(4)/CONSTS(50) ** CONSTS(51))/(1.00000+STATES(2)/CONSTS(49) ** CONSTS(51)+STATES(4)/CONSTS(50) ** CONSTS(51)))
      ALGBRC(66) =  CONSTS(52)*((STATES(3)/CONSTS(53) ** CONSTS(55) - STATES(4)/CONSTS(54) ** CONSTS(55))/(1.00000+STATES(3)/CONSTS(53) ** CONSTS(55)+STATES(4)/CONSTS(54) ** CONSTS(55)))
      ALGBRC(68) =  CONSTS(56)*(STATES(4) - STATES(3))
      ALGBRC(69) =  CONSTS(57)*(STATES(2) - STATES(3))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END