C
C There are a total of 108 entries in the algebraic variable array.
C There are a total of 26 entries in each of the rate and state variable arrays.
C There are a total of 80 entries in the constant variable array.
C
C
C VOI is time in component Environment (ms).
C CONSTS(1) is R in component Environment (J_per_moleK).
C CONSTS(2) is T in component Environment (kelvin).
C CONSTS(3) is F in component Environment (coulomb_per_mmole).
C CONSTS(4) is K_o in component Environment (mM).
C CONSTS(5) is Ca_o in component Environment (mM).
C CONSTS(6) is Na_o in component Environment (mM).
C CONSTS(76) is FonRT in component Environment (per_mV).
C STATES(1) is V in component cell (mV).
C ALGBRC(50) is xik1 in component IK1 (nA_per_nF).
C ALGBRC(61) is xito in component Ito (nA_per_nF).
C ALGBRC(66) is xiNaK in component INaK (nA_per_nF).
C CONSTS(7) is wca in component cell (mV_per_uM).
C ALGBRC(103) is xiNaCa in component INaCa (nA_per_nF).
C ALGBRC(94) is xica in component ICaL (nA_per_nF).
C ALGBRC(108) is Itotal in component cell (nA_per_nF).
C ALGBRC(107) is xina in component INa (nA_per_nF).
C ALGBRC(52) is xikr in component IKr (nA_per_nF).
C ALGBRC(105) is xiks in component IKs (nA_per_nF).
C ALGBRC(11) is i_Stim in component cell (nA_per_nF).
C CONSTS(8) is stim_offset in component cell (ms).
C CONSTS(9) is stim_period in component cell (ms).
C CONSTS(10) is stim_duration in component cell (ms).
C CONSTS(11) is stim_amplitude in component cell (nA_per_nF).
C ALGBRC(4) is past in component cell (ms).
C ALGBRC(106) is ena in component reversal_potentials (mV).
C STATES(2) is xm in component INa (dimensionless).
C STATES(3) is xh in component INa (dimensionless).
C STATES(4) is xj in component INa (dimensionless).
C CONSTS(12) is gna in component INa (uS_per_nF).
C ALGBRC(1) is am in component INa (per_ms).
C ALGBRC(8) is bm in component INa (per_ms).
C ALGBRC(2) is ah in component INa (per_ms).
C ALGBRC(9) is bh in component INa (per_ms).
C ALGBRC(3) is aj in component INa (per_ms).
C ALGBRC(10) is bj in component INa (per_ms).
C STATES(5) is Ca_dyad in component Ca (uM).
C ALGBRC(89) is csm in component Ca (mM).
C STATES(6) is c1 in component ICaL (dimensionless).
C STATES(7) is c2 in component ICaL (dimensionless).
C STATES(8) is xi1ca in component ICaL (dimensionless).
C STATES(9) is xi1ba in component ICaL (dimensionless).
C STATES(10) is xi2ca in component ICaL (dimensionless).
C STATES(11) is xi2ba in component ICaL (dimensionless).
C CONSTS(13) is gca in component ICaL (mmole_per_coulomb_cm).
C CONSTS(14) is pca in component ICaL (cm_per_s).
C ALGBRC(15) is za in component ICaL (dimensionless).
C ALGBRC(19) is poinf in component ICaL (dimensionless).
C ALGBRC(27) is fca in component ICaL (dimensionless).
C CONSTS(15) is vth in component ICaL (mV).
C CONSTS(16) is s6 in component ICaL (mV).
C CONSTS(17) is vx in component ICaL (mV).
C CONSTS(18) is sx in component ICaL (mV).
C CONSTS(19) is vy in component ICaL (mV).
C CONSTS(20) is sy in component ICaL (mV).
C CONSTS(21) is vyr in component ICaL (mV).
C CONSTS(22) is syr in component ICaL (mV).
C CONSTS(23) is cat in component ICaL (uM).
C CONSTS(24) is cpt in component ICaL (uM).
C ALGBRC(24) is alpha in component ICaL (per_ms).
C ALGBRC(26) is beta in component ICaL (per_ms).
C ALGBRC(29) is k1 in component ICaL (per_ms).
C CONSTS(25) is k2 in component ICaL (per_ms).
C CONSTS(26) is k1t in component ICaL (per_ms).
C CONSTS(27) is k2t in component ICaL (per_ms).
C ALGBRC(32) is k3 in component ICaL (per_ms).
C ALGBRC(33) is k3t in component ICaL (per_ms).
C ALGBRC(40) is k6 in component ICaL (per_ms).
C ALGBRC(41) is k5 in component ICaL (per_ms).
C ALGBRC(42) is k6t in component ICaL (per_ms).
C ALGBRC(43) is k5t in component ICaL (per_ms).
C ALGBRC(44) is k4 in component ICaL (per_ms).
C ALGBRC(45) is k4t in component ICaL (per_ms).
C CONSTS(28) is r1 in component ICaL (per_ms).
C CONSTS(29) is r2 in component ICaL (per_ms).
C ALGBRC(28) is s1 in component ICaL (per_ms).
C CONSTS(30) is s1t in component ICaL (per_ms).
C ALGBRC(30) is s2 in component ICaL (per_ms).
C CONSTS(77) is s2t in component ICaL (per_ms).
C ALGBRC(35) is recov in component ICaL (ms).
C CONSTS(31) is tca in component ICaL (ms).
C ALGBRC(36) is tau_ca in component ICaL (ms).
C ALGBRC(37) is tauca in component ICaL (ms).
C ALGBRC(38) is tauba in component ICaL (ms).
C CONSTS(32) is taupo in component ICaL (ms).
C CONSTS(33) is tau3 in component ICaL (ms).
C ALGBRC(34) is Pr in component ICaL (dimensionless).
C ALGBRC(39) is Ps in component ICaL (dimensionless).
C ALGBRC(31) is poi in component ICaL (dimensionless).
C ALGBRC(46) is po in component ICaL (dimensionless).
C ALGBRC(90) is rxa in component ICaL (mA_per_cm2).
C ALGBRC(91) is jca in component ICaL (uM_per_ms).
C CONSTS(80) is ek in component reversal_potentials (mV).
C CONSTS(34) is gkix in component IK1 (uS_per_nF).
C ALGBRC(47) is aki in component IK1 (per_ms).
C ALGBRC(48) is bki in component IK1 (per_ms).
C ALGBRC(49) is xkin in component IK1 (dimensionless).
C STATES(12) is xr in component IKr (dimensionless).
C CONSTS(35) is gkr in component IKr (uS_per_nF).
C ALGBRC(5) is xkrv1 in component IKr (per_ms).
C ALGBRC(12) is xkrv2 in component IKr (per_ms).
C ALGBRC(16) is taukr in component IKr (ms).
C ALGBRC(20) is xkrinf in component IKr (dimensionless).
C ALGBRC(51) is rg in component IKr (dimensionless).
C STATES(13) is Ca_i in component Ca (uM).
C CONSTS(36) is gks in component IKs (uS_per_nF).
C STATES(14) is xs1 in component IKs (dimensionless).
C STATES(15) is xs2 in component IKs (dimensionless).
C ALGBRC(104) is eks in component reversal_potentials (mV).
C ALGBRC(6) is xs1ss in component IKs (dimensionless).
C ALGBRC(13) is xs2ss in component IKs (dimensionless).
C ALGBRC(17) is tauxs1 in component IKs (ms).
C ALGBRC(21) is tauxs2 in component IKs (ms).
C ALGBRC(53) is gksx in component IKs (dimensionless).
C ALGBRC(57) is xitos in component Ito (nA_per_nF).
C ALGBRC(59) is xitof in component Ito (nA_per_nF).
C STATES(16) is xtos in component Ito (dimensionless).
C STATES(17) is ytos in component Ito (dimensionless).
C STATES(18) is xtof in component Ito (dimensionless).
C STATES(19) is ytof in component Ito (dimensionless).
C CONSTS(37) is gtos in component Ito (uS_per_nF).
C CONSTS(38) is gtof in component Ito (uS_per_nF).
C ALGBRC(7) is rt1 in component Ito (dimensionless).
C ALGBRC(54) is rt2 in component Ito (dimensionless).
C ALGBRC(56) is rt3 in component Ito (dimensionless).
C ALGBRC(14) is rt4 in component Ito (dimensionless).
C ALGBRC(58) is rt5 in component Ito (dimensionless).
C ALGBRC(18) is xtos_inf in component Ito (dimensionless).
C ALGBRC(60) is ytos_inf in component Ito (dimensionless).
C ALGBRC(22) is xtof_inf in component Ito (dimensionless).
C ALGBRC(62) is ytof_inf in component Ito (dimensionless).
C ALGBRC(55) is rs_inf in component Ito (dimensionless).
C ALGBRC(23) is txs in component Ito (ms).
C ALGBRC(63) is tys in component Ito (ms).
C ALGBRC(25) is txf in component Ito (ms).
C ALGBRC(65) is tyf in component Ito (ms).
C STATES(20) is Na_i in component Na (mM).
C CONSTS(39) is gNaK in component INaK (nA_per_nF).
C CONSTS(40) is xkmko in component INaK (mM).
C CONSTS(41) is xkmnai in component INaK (mM).
C CONSTS(78) is sigma in component INaK (dimensionless).
C ALGBRC(64) is fNaK in component INaK (dimensionless).
C STATES(21) is Ca_submem in component Ca (uM).
C CONSTS(42) is gNaCa in component INaCa (uM_per_ms).
C ALGBRC(68) is aloss in component INaCa (dimensionless).
C ALGBRC(98) is yz1 in component INaCa (mM4).
C ALGBRC(99) is yz2 in component INaCa (mM4).
C ALGBRC(69) is yz3 in component INaCa (mM4).
C ALGBRC(100) is yz4 in component INaCa (mM4).
C ALGBRC(96) is zw3 in component INaCa (mM4).
C ALGBRC(67) is zw4 in component INaCa (dimensionless).
C ALGBRC(101) is zw8 in component INaCa (mM4).
C ALGBRC(102) is jNaCa in component INaCa (uM_per_ms).
C CONSTS(43) is xkdna in component INaCa (uM).
C CONSTS(44) is xmcao in component INaCa (mM).
C CONSTS(45) is xmnao in component INaCa (mM).
C CONSTS(46) is xmnai in component INaCa (mM).
C CONSTS(47) is xmcai in component INaCa (mM).
C STATES(22) is Ca_NSR in component Ca (uM).
C ALGBRC(88) is dCa_JSR in component Ca (uM_per_ms).
C CONSTS(48) is cstar in component Irel (uM).
C STATES(23) is Ca_JSR in component Irel (uM).
C CONSTS(49) is gryr in component Irel (per_ms).
C CONSTS(50) is gbarsr in component Irel (dimensionless).
C CONSTS(51) is gdyad in component Irel (mmole_per_coulomb_cm).
C CONSTS(52) is ax in component Irel (per_mV).
C CONSTS(53) is ay in component Irel (per_mV).
C CONSTS(54) is av in component Irel (per_ms).
C CONSTS(79) is bv in component Irel (uM_per_ms).
C ALGBRC(70) is Qr0 in component Irel (uM_per_ms).
C ALGBRC(71) is Qr in component Irel (uM_per_ms).
C ALGBRC(72) is sparkV in component Irel (dimensionless).
C ALGBRC(92) is spark_rate in component Irel (per_ms).
C CONSTS(55) is taua in component Irel (ms).
C CONSTS(56) is taur in component Irel (ms).
C ALGBRC(93) is xirp in component Irel (uM_per_ms).
C ALGBRC(95) is xicap in component Irel (uM_per_ms).
C ALGBRC(97) is xiryr in component Irel (uM_per_ms).
C STATES(24) is xir in component Irel (uM_per_ms).
C ALGBRC(73) is jup in component Ileak_Iup_Ixfer (uM_per_ms).
C ALGBRC(74) is jleak in component Ileak_Iup_Ixfer (uM_per_ms).
C CONSTS(57) is cup in component Ileak_Iup_Ixfer (uM).
C CONSTS(58) is kj in component Ileak_Iup_Ixfer (uM).
C CONSTS(59) is vup in component Ileak_Iup_Ixfer (uM_per_ms).
C CONSTS(60) is gleak in component Ileak_Iup_Ixfer (per_ms).
C CONSTS(61) is bcal in component Ca (uM).
C CONSTS(62) is xkcal in component Ca (uM).
C CONSTS(63) is srmax in component Ca (uM).
C CONSTS(64) is srkd in component Ca (uM).
C CONSTS(65) is bmem in component Ca (uM).
C CONSTS(66) is kmem in component Ca (uM).
C CONSTS(67) is bsar in component Ca (uM).
C CONSTS(68) is ksar in component Ca (uM).
C ALGBRC(75) is bpxs in component Ca (dimensionless).
C ALGBRC(76) is spxs in component Ca (dimensionless).
C ALGBRC(77) is mempxs in component Ca (dimensionless).
C ALGBRC(78) is sarpxs in component Ca (dimensionless).
C ALGBRC(79) is dcsib in component Ca (dimensionless).
C ALGBRC(80) is bpxi in component Ca (dimensionless).
C ALGBRC(81) is spxi in component Ca (dimensionless).
C ALGBRC(82) is mempxi in component Ca (dimensionless).
C ALGBRC(83) is sarpxi in component Ca (dimensionless).
C ALGBRC(84) is dciib in component Ca (dimensionless).
C CONSTS(69) is xkon in component Ca (per_uM_per_ms).
C CONSTS(70) is xkoff in component Ca (per_ms).
C CONSTS(71) is btrop in component Ca (uM).
C ALGBRC(87) is xbi in component Ca (uM_per_ms).
C ALGBRC(86) is xbs in component Ca (uM_per_ms).
C STATES(25) is tropi in component Ca (uM).
C STATES(26) is trops in component Ca (uM).
C CONSTS(72) is taud in component Ca (ms).
C CONSTS(73) is taups in component Ca (ms).
C ALGBRC(85) is jd in component Ca (uM_per_ms).
C CONSTS(74) is K_i in component reversal_potentials (mM).
C CONSTS(75) is prNaK in component reversal_potentials (dimensionless).
C RATES(1) is d/dt V in component cell (mV).
C RATES(3) is d/dt xh in component INa (dimensionless).
C RATES(4) is d/dt xj in component INa (dimensionless).
C RATES(2) is d/dt xm in component INa (dimensionless).
C RATES(6) is d/dt c1 in component ICaL (dimensionless).
C RATES(7) is d/dt c2 in component ICaL (dimensionless).
C RATES(8) is d/dt xi1ca in component ICaL (dimensionless).
C RATES(9) is d/dt xi1ba in component ICaL (dimensionless).
C RATES(10) is d/dt xi2ca in component ICaL (dimensionless).
C RATES(11) is d/dt xi2ba in component ICaL (dimensionless).
C RATES(12) is d/dt xr in component IKr (dimensionless).
C RATES(14) is d/dt xs1 in component IKs (dimensionless).
C RATES(15) is d/dt xs2 in component IKs (dimensionless).
C RATES(16) is d/dt xtos in component Ito (dimensionless).
C RATES(17) is d/dt ytos in component Ito (dimensionless).
C RATES(18) is d/dt xtof in component Ito (dimensionless).
C RATES(19) is d/dt ytof in component Ito (dimensionless).
C RATES(23) is d/dt Ca_JSR in component Irel (uM).
C RATES(24) is d/dt xir in component Irel (uM_per_ms).
C RATES(20) is d/dt Na_i in component Na (mM).
C RATES(5) is d/dt Ca_dyad in component Ca (uM).
C RATES(21) is d/dt Ca_submem in component Ca (uM).
C RATES(13) is d/dt Ca_i in component Ca (uM).
C RATES(22) is d/dt Ca_NSR in component Ca (uM).
C RATES(25) is d/dt tropi in component Ca (uM).
C RATES(26) is d/dt trops in component Ca (uM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 8.314472
      CONSTS(2) = 308
      CONSTS(3) = 96.4853415
      CONSTS(4) = 5.4
      CONSTS(5) = 1.8
      CONSTS(6) = 136
      STATES(1) = -87.169816169406
      CONSTS(7) = 8
      CONSTS(8) = 0
      CONSTS(9) = 400
      CONSTS(10) = 3
      CONSTS(11) = -15
      STATES(2) = 0.001075453357
      STATES(3) = 0.990691306716
      STATES(4) = 0.993888937283
      CONSTS(12) = 12
      STATES(5) = 1.716573130685
      STATES(6) = 0.000018211252
      STATES(7) = 0.979322592773
      STATES(8) = 0.001208153482
      STATES(9) = 0.000033616596
      STATES(10) = 0.004173008466
      STATES(11) = 0.015242594688
      CONSTS(13) = 182
      CONSTS(14) = 0.00054
      CONSTS(15) = 0
      CONSTS(16) = 8
      CONSTS(17) = -40
      CONSTS(18) = 3
      CONSTS(19) = -40
      CONSTS(20) = 4
      CONSTS(21) = -40
      CONSTS(22) = 11.32
      CONSTS(23) = 3
      CONSTS(24) = 6.09365
      CONSTS(25) = 1.03615e-4
      CONSTS(26) = 0.00413
      CONSTS(27) = 0.00224
      CONSTS(28) = 0.3
      CONSTS(29) = 3
      CONSTS(30) = 0.00195
      CONSTS(31) = 78.0329
      CONSTS(32) = 1
      CONSTS(33) = 3
      CONSTS(34) = 0.3
      STATES(12) = 0.007074239331
      CONSTS(35) = 0.0125
      STATES(13) = 0.256752008084
      CONSTS(36) = 0.1386
      STATES(14) = 0.048267587131
      STATES(15) = 0.105468807033
      STATES(16) = 0.00364776906
      STATES(17) = 0.174403618112
      STATES(18) = 0.003643592594
      STATES(19) = 0.993331326442
      CONSTS(37) = 0.04
      CONSTS(38) = 0.11
      STATES(20) = 11.441712311614
      CONSTS(39) = 1.5
      CONSTS(40) = 1.5
      CONSTS(41) = 12
      STATES(21) = 0.226941113355
      CONSTS(42) = 0.84
      CONSTS(43) = 0.3
      CONSTS(44) = 1.3
      CONSTS(45) = 87.5
      CONSTS(46) = 12.3
      CONSTS(47) = 0.0036
      STATES(22) = 104.450004990523
      CONSTS(48) = 90
      STATES(23) = 97.505463697266
      CONSTS(49) = 2.58079
      CONSTS(50) = 26841.8
      CONSTS(51) = 9000
      CONSTS(52) = 0.3576
      CONSTS(53) = 0.05
      CONSTS(54) = 11.3
      CONSTS(55) = 100
      CONSTS(56) = 30
      STATES(24) = 0.006679257264
      CONSTS(57) = 0.5
      CONSTS(58) = 50
      CONSTS(59) = 0.4
      CONSTS(60) = 0.00002069
      CONSTS(61) = 24
      CONSTS(62) = 7
      CONSTS(63) = 47
      CONSTS(64) = 0.6
      CONSTS(65) = 15
      CONSTS(66) = 0.3
      CONSTS(67) = 42
      CONSTS(68) = 13
      CONSTS(69) = 0.0327
      CONSTS(70) = 0.0196
      CONSTS(71) = 70
      STATES(25) = 22.171689894953
      STATES(26) = 19.864701949854
      CONSTS(72) = 4
      CONSTS(73) = 0.5
      CONSTS(74) = 140
      CONSTS(75) = 0.01833
      CONSTS(76) = CONSTS(3)/( CONSTS(1)*CONSTS(2))
      CONSTS(77) = ( (( CONSTS(30)*CONSTS(28))/CONSTS(29))*CONSTS(27))/CONSTS(26)
      CONSTS(78) = (EXP(CONSTS(6)/67.3000) - 1.00000)/7.00000
      CONSTS(79) =  (1.00000 - CONSTS(54))*CONSTS(48) - 50.0000
      CONSTS(80) =  (1.00000/CONSTS(76))*log(CONSTS(4)/CONSTS(74))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(23) = (STATES(22) - STATES(23))/CONSTS(55)
      ALGBRC(2) = TERNRY(STATES(1).LT.- 40.0000,  0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
      ALGBRC(9) = TERNRY(STATES(1).LT.- 40.0000,  3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)), 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))))
      RATES(3) =  ALGBRC(2)*(1.00000 - STATES(3)) -  ALGBRC(9)*STATES(3)
      ALGBRC(3) = TERNRY(STATES(1).LT.- 40.0000, ( ( - 127140.*EXP( 0.244400*STATES(1)) -  3.47400e-05*EXP( - 0.0439100*STATES(1)))*1.00000*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))), 0.00000)
      ALGBRC(10) = TERNRY(STATES(1).LT.- 40.0000, ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))), ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))
      RATES(4) =  ALGBRC(3)*(1.00000 - STATES(4)) -  ALGBRC(10)*STATES(4)
      ALGBRC(1) = TERNRY(ABS(STATES(1)+47.1300).GT.0.00100000, ( 0.320000*1.00000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300))), 3.20000)
      ALGBRC(8) =  0.0800000*EXP(- STATES(1)/11.0000)
      RATES(2) =  ALGBRC(1)*(1.00000 - STATES(2)) -  ALGBRC(8)*STATES(2)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
      ALGBRC(17) = TERNRY(ABS(STATES(1)+30.0000).LT.0.00100000/0.0687000, 1.00000/(7.19000e-05/0.148000+0.000131000/0.0687000), 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(14) = (ALGBRC(6) - STATES(14))/ALGBRC(17)
      ALGBRC(5) = 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(12) = TERNRY(ABS(STATES(1)+10.0000).GT.0.00100000, ( 0.000610000*1.00000*(STATES(1)+10.0000))/(EXP( 0.145000*(STATES(1)+10.0000)) - 1.00000), 0.000610000/0.145000)
      ALGBRC(16) = 1.00000/(ALGBRC(5)+ALGBRC(12))
      ALGBRC(20) = 1.00000/(1.00000+EXP(- (STATES(1)+50.0000)/7.50000))
      RATES(12) = (ALGBRC(20) - STATES(12))/ALGBRC(16)
      ALGBRC(13) = ALGBRC(6)
      ALGBRC(21) =  4.00000*ALGBRC(17)
      RATES(15) = (ALGBRC(13) - STATES(15))/ALGBRC(21)
      ALGBRC(7) = - (STATES(1)+3.00000)/15.0000
      ALGBRC(18) = 1.00000/(1.00000+EXP(ALGBRC(7)))
      ALGBRC(23) = 9.00000/(1.00000+EXP(- ALGBRC(7)))+0.500000
      RATES(16) = (ALGBRC(18) - STATES(16))/ALGBRC(23)
      ALGBRC(22) = ALGBRC(18)
      ALGBRC(14) = ( (- STATES(1)/30.0000)*STATES(1))/30.0000
      ALGBRC(25) =  3.50000*EXP(ALGBRC(14))+1.50000
      RATES(18) = (ALGBRC(22) - STATES(18))/ALGBRC(25)
      ALGBRC(19) = 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(15))/CONSTS(16)))
      ALGBRC(24) = ALGBRC(19)/CONSTS(32)
      ALGBRC(26) = (1.00000 - ALGBRC(19))/CONSTS(32)
      ALGBRC(27) = 1.00000/(1.00000+CONSTS(23)/STATES(5) ** 3.00000)
      ALGBRC(35) = 10.0000+ 4954.00*EXP(STATES(1)/15.6000)
      ALGBRC(36) = CONSTS(31)/(1.00000+STATES(5)/CONSTS(24) ** 4.00000)+0.100000
      ALGBRC(34) = 1.00000 - 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(19))/CONSTS(20)))
      ALGBRC(37) =  (ALGBRC(35) - ALGBRC(36))*ALGBRC(34)+ALGBRC(36)
      ALGBRC(39) = 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(21))/CONSTS(22)))
      ALGBRC(40) = ( ALGBRC(27)*ALGBRC(39))/ALGBRC(37)
      ALGBRC(41) = (1.00000 - ALGBRC(39))/ALGBRC(37)
      ALGBRC(38) =  (ALGBRC(35) - 450.000)*ALGBRC(34)+450.000
      ALGBRC(42) = ALGBRC(39)/ALGBRC(38)
      ALGBRC(43) = (1.00000 - ALGBRC(39))/ALGBRC(38)
      RATES(7) = ( ALGBRC(26)*STATES(6)+ ALGBRC(41)*STATES(10)+ ALGBRC(43)*STATES(11)) -  (ALGBRC(40)+ALGBRC(42)+ALGBRC(24))*STATES(7)
      ALGBRC(31) = 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(17))/CONSTS(18)))
      ALGBRC(32) = (1.00000 - ALGBRC(31))/CONSTS(33)
      ALGBRC(29) =  0.0241680*ALGBRC(27)
      ALGBRC(44) = ( (( (( ALGBRC(32)*ALGBRC(24))/ALGBRC(26))*ALGBRC(29))/CONSTS(25))*ALGBRC(41))/ALGBRC(40)
      RATES(10) = ( ALGBRC(32)*STATES(8)+ ALGBRC(40)*STATES(7)) -  (ALGBRC(41)+ALGBRC(44))*STATES(10)
      ALGBRC(33) = ALGBRC(32)
      ALGBRC(45) = ( (( (( ALGBRC(33)*ALGBRC(24))/ALGBRC(26))*CONSTS(26))/CONSTS(27))*ALGBRC(43))/ALGBRC(42)
      RATES(11) = ( ALGBRC(33)*STATES(9)+ ALGBRC(42)*STATES(7)) -  (ALGBRC(43)+ALGBRC(45))*STATES(11)
      ALGBRC(46) = (((((1.00000 - STATES(8)) - STATES(10)) - STATES(9)) - STATES(11)) - STATES(6)) - STATES(7)
      RATES(6) = ( ALGBRC(24)*STATES(7)+ CONSTS(25)*STATES(8)+ CONSTS(27)*STATES(9)+ CONSTS(29)*ALGBRC(46)) -  (ALGBRC(26)+CONSTS(28)+CONSTS(26)+ALGBRC(29))*STATES(6)
      ALGBRC(28) =  0.0182688*ALGBRC(27)
      ALGBRC(30) = ( (( ALGBRC(28)*CONSTS(28))/CONSTS(29))*CONSTS(25))/ALGBRC(29)
      RATES(8) = ( ALGBRC(29)*STATES(6)+ ALGBRC(44)*STATES(10)+ ALGBRC(28)*ALGBRC(46)) -  (ALGBRC(32)+CONSTS(25)+ALGBRC(30))*STATES(8)
      RATES(9) = ( CONSTS(26)*STATES(6)+ ALGBRC(45)*STATES(11)+ CONSTS(30)*ALGBRC(46)) -  (ALGBRC(33)+CONSTS(27)+CONSTS(77))*STATES(9)
      ALGBRC(54) = (STATES(1)+33.5000)/10.0000
      ALGBRC(60) = 1.00000/(1.00000+EXP(ALGBRC(54)))
      ALGBRC(56) = (STATES(1)+60.0000)/10.0000
      ALGBRC(63) = 3000.00/(1.00000+EXP(ALGBRC(56)))+30.0000
      RATES(17) = (ALGBRC(60) - STATES(17))/ALGBRC(63)
      ALGBRC(62) = ALGBRC(60)
      ALGBRC(58) = (STATES(1)+33.5000)/10.0000
      ALGBRC(65) = 20.0000/(1.00000+EXP(ALGBRC(58)))+20.0000
      RATES(19) = (ALGBRC(62) - STATES(19))/ALGBRC(65)
      ALGBRC(73) = ( CONSTS(59)*STATES(13)*STATES(13))/( STATES(13)*STATES(13)+ CONSTS(57)*CONSTS(57))
      ALGBRC(74) =  (( CONSTS(60)*STATES(22)*STATES(22))/( STATES(22)*STATES(22)+ CONSTS(58)*CONSTS(58)))*( STATES(22)*16.6670 - STATES(13))
      ALGBRC(80) = ( CONSTS(61)*CONSTS(62))/( (CONSTS(62)+STATES(13))*(CONSTS(62)+STATES(13)))
      ALGBRC(81) = ( CONSTS(63)*CONSTS(64))/( (CONSTS(64)+STATES(13))*(CONSTS(64)+STATES(13)))
      ALGBRC(82) = ( CONSTS(65)*CONSTS(66))/( (CONSTS(66)+STATES(13))*(CONSTS(66)+STATES(13)))
      ALGBRC(83) = ( CONSTS(67)*CONSTS(68))/( (CONSTS(68)+STATES(13))*(CONSTS(68)+STATES(13)))
      ALGBRC(84) = 1.00000/(1.00000+ALGBRC(80)+ALGBRC(81)+ALGBRC(82)+ALGBRC(83))
      ALGBRC(87) =  CONSTS(69)*STATES(13)*(CONSTS(71) - STATES(25)) -  CONSTS(70)*STATES(25)
      ALGBRC(85) = (STATES(21) - STATES(13))/CONSTS(72)
      RATES(13) =  ALGBRC(84)*(((ALGBRC(85) - ALGBRC(73))+ALGBRC(74)) - ALGBRC(87))
      RATES(25) = ALGBRC(87)
      ALGBRC(86) =  CONSTS(69)*STATES(21)*(CONSTS(71) - STATES(26)) -  CONSTS(70)*STATES(26)
      RATES(26) = ALGBRC(86)
      ALGBRC(88) = (- STATES(24)+ALGBRC(73)) - ALGBRC(74)
      RATES(22) = ALGBRC(88)
      ALGBRC(70) = TERNRY(STATES(23).GT.50.0000.AND.STATES(23).LT.CONSTS(48), (STATES(23) - 50.0000)/1.00000, TERNRY(STATES(23).GE.CONSTS(48),  CONSTS(54)*STATES(23)+CONSTS(79), 0.00000)
      ALGBRC(71) = ( STATES(22)*ALGBRC(70))/CONSTS(48)
      ALGBRC(89) = STATES(21)/1000.00
      ALGBRC(15) =  STATES(1)*2.00000*CONSTS(76)
      ALGBRC(90) = TERNRY(ABS(ALGBRC(15)).LT.0.00100000, ( 4.00000*CONSTS(14)*CONSTS(3)*CONSTS(76)*( ALGBRC(89)*EXP(ALGBRC(15)) -  0.341000*CONSTS(5)))/( 2.00000*CONSTS(76)), ( 4.00000*CONSTS(14)*STATES(1)*CONSTS(3)*CONSTS(76)*( ALGBRC(89)*EXP(ALGBRC(15)) -  0.341000*CONSTS(5)))/(EXP(ALGBRC(15)) - 1.00000))
      ALGBRC(72) = EXP( - CONSTS(53)*(STATES(1)+30.0000))/(1.00000+EXP( - CONSTS(53)*(STATES(1)+30.0000)))
      ALGBRC(92) =  (CONSTS(49)/1.00000)*ALGBRC(46)*ABS(ALGBRC(90))*ALGBRC(72)
      RATES(24) =  ALGBRC(92)*ALGBRC(71) - ( STATES(24)*(1.00000 - ( CONSTS(56)*ALGBRC(88))/STATES(22)))/CONSTS(56)
      ALGBRC(93) = ( (( ALGBRC(46)*ALGBRC(71)*ABS(ALGBRC(90))*CONSTS(50))/1.00000)*EXP( - CONSTS(52)*(STATES(1)+30.0000)))/(1.00000+EXP( - CONSTS(52)*(STATES(1)+30.0000)))
      ALGBRC(95) =  ALGBRC(46)*CONSTS(51)*ABS(ALGBRC(90))
      ALGBRC(97) = ALGBRC(93)+ALGBRC(95)
      RATES(5) = ALGBRC(97) - (STATES(5) - STATES(21))/CONSTS(73)
      ALGBRC(91) =  CONSTS(13)*ALGBRC(46)*ALGBRC(90)
      ALGBRC(68) = 1.00000/(1.00000+CONSTS(43)/STATES(21) ** 3.00000)
      ALGBRC(96) =  STATES(20) ** 3.00000*CONSTS(5)*EXP( STATES(1)*0.350000*CONSTS(76)) -  CONSTS(6) ** 3.00000*ALGBRC(89)*EXP( STATES(1)*(0.350000 - 1.00000)*CONSTS(76))
      ALGBRC(67) = 1.00000+ 0.200000*EXP( STATES(1)*(0.350000 - 1.00000)*CONSTS(76))
      ALGBRC(98) =  CONSTS(44)*STATES(20) ** 3.00000+ CONSTS(45) ** 3.00000*ALGBRC(89)
      ALGBRC(99) =  CONSTS(46) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(89)/CONSTS(47))
      ALGBRC(69) =  CONSTS(47)*CONSTS(6) ** 3.00000*(1.00000+STATES(20)/CONSTS(46) ** 3.00000)
      ALGBRC(100) =  STATES(20) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(89)
      ALGBRC(101) = ALGBRC(98)+ALGBRC(99)+ALGBRC(69)+ALGBRC(100)
      ALGBRC(102) = ( CONSTS(42)*ALGBRC(68)*ALGBRC(96))/( ALGBRC(67)*ALGBRC(101))
      ALGBRC(75) = ( CONSTS(61)*CONSTS(62))/( (CONSTS(62)+STATES(21))*(CONSTS(62)+STATES(21)))
      ALGBRC(76) = ( CONSTS(63)*CONSTS(64))/( (CONSTS(64)+STATES(21))*(CONSTS(64)+STATES(21)))
      ALGBRC(77) = ( CONSTS(65)*CONSTS(66))/( (CONSTS(66)+STATES(21))*(CONSTS(66)+STATES(21)))
      ALGBRC(78) = ( CONSTS(67)*CONSTS(68))/( (CONSTS(68)+STATES(21))*(CONSTS(68)+STATES(21)))
      ALGBRC(79) = 1.00000/(1.00000+ALGBRC(75)+ALGBRC(76)+ALGBRC(77)+ALGBRC(78))
      RATES(21) =  ALGBRC(79)*( 50.0000*(((STATES(24) - ALGBRC(85)) - ALGBRC(91))+ALGBRC(102)) - ALGBRC(86))
      ALGBRC(64) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(1)*CONSTS(76))+ 0.0365000*CONSTS(78)*EXP( - STATES(1)*CONSTS(76)))
      ALGBRC(66) = ( (( CONSTS(39)*ALGBRC(64)*STATES(20))/(STATES(20)+CONSTS(41)))*CONSTS(4))/(CONSTS(4)+CONSTS(40))
      ALGBRC(103) =  CONSTS(7)*ALGBRC(102)
      ALGBRC(106) =  (1.00000/CONSTS(76))*log(CONSTS(6)/STATES(20))
      ALGBRC(107) =  CONSTS(12)*STATES(3)*STATES(4)*STATES(2)*STATES(2)*STATES(2)*(STATES(1) - ALGBRC(106))
      RATES(20) = - (ALGBRC(107)+ 3.00000*ALGBRC(66)+ 3.00000*ALGBRC(103))/( CONSTS(7)*1000.00)
      ALGBRC(47) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - CONSTS(80)) - 59.2150)))
      ALGBRC(48) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(80))+5.47600))+ 1.00000*EXP( 0.0617500*((STATES(1) - CONSTS(80)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(80))+4.75300)))
      ALGBRC(49) = ALGBRC(47)/(ALGBRC(47)+ALGBRC(48))
      ALGBRC(50) =  CONSTS(34)* (CONSTS(4)/5.40000) ** (1.0 / 2)*ALGBRC(49)*(STATES(1) - CONSTS(80))
      ALGBRC(55) = 1.00000/(1.00000+EXP(ALGBRC(54)))
      ALGBRC(57) =  CONSTS(37)*STATES(16)*(STATES(17)+ 0.500000*ALGBRC(55))*(STATES(1) - CONSTS(80))
      ALGBRC(59) =  CONSTS(38)*STATES(18)*STATES(19)*(STATES(1) - CONSTS(80))
      ALGBRC(61) = ALGBRC(57)+ALGBRC(59)
      ALGBRC(94) =  2.00000*CONSTS(7)*ALGBRC(91)
      ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(1)+33.0000)/22.4000))
      ALGBRC(52) =  CONSTS(35)* (CONSTS(4)/5.40000) ** (1.0 / 2)*STATES(12)*ALGBRC(51)*(STATES(1) - CONSTS(80))
      ALGBRC(104) =  (1.00000/CONSTS(76))*log((CONSTS(4)+ CONSTS(75)*CONSTS(6))/(CONSTS(74)+ CONSTS(75)*STATES(20)))
      ALGBRC(53) = 1.00000+0.800000/(1.00000+0.500000/STATES(13) ** 3.00000)
      ALGBRC(105) =  CONSTS(36)*ALGBRC(53)*STATES(14)*STATES(15)*(STATES(1) - ALGBRC(104))
      ALGBRC(4) =  INT(VOI/CONSTS(9))*CONSTS(9)
      ALGBRC(11) = TERNRY(VOI - ALGBRC(4).GE.CONSTS(8).AND.VOI - ALGBRC(4).LE.CONSTS(8)+CONSTS(10), CONSTS(11), 0.00000)
      ALGBRC(108) = - (ALGBRC(107)+ALGBRC(50)+ALGBRC(52)+ALGBRC(105)+ALGBRC(61)+ALGBRC(103)+ALGBRC(94)+ALGBRC(66)+ALGBRC(11))
      RATES(1) = ALGBRC(108)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = TERNRY(STATES(1).LT.- 40.0000,  0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
      ALGBRC(9) = TERNRY(STATES(1).LT.- 40.0000,  3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)), 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))))
      ALGBRC(3) = TERNRY(STATES(1).LT.- 40.0000, ( ( - 127140.*EXP( 0.244400*STATES(1)) -  3.47400e-05*EXP( - 0.0439100*STATES(1)))*1.00000*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))), 0.00000)
      ALGBRC(10) = TERNRY(STATES(1).LT.- 40.0000, ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))), ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))
      ALGBRC(1) = TERNRY(ABS(STATES(1)+47.1300).GT.0.00100000, ( 0.320000*1.00000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300))), 3.20000)
      ALGBRC(8) =  0.0800000*EXP(- STATES(1)/11.0000)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
      ALGBRC(17) = TERNRY(ABS(STATES(1)+30.0000).LT.0.00100000/0.0687000, 1.00000/(7.19000e-05/0.148000+0.000131000/0.0687000), 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(5) = 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(12) = TERNRY(ABS(STATES(1)+10.0000).GT.0.00100000, ( 0.000610000*1.00000*(STATES(1)+10.0000))/(EXP( 0.145000*(STATES(1)+10.0000)) - 1.00000), 0.000610000/0.145000)
      ALGBRC(16) = 1.00000/(ALGBRC(5)+ALGBRC(12))
      ALGBRC(20) = 1.00000/(1.00000+EXP(- (STATES(1)+50.0000)/7.50000))
      ALGBRC(13) = ALGBRC(6)
      ALGBRC(21) =  4.00000*ALGBRC(17)
      ALGBRC(7) = - (STATES(1)+3.00000)/15.0000
      ALGBRC(18) = 1.00000/(1.00000+EXP(ALGBRC(7)))
      ALGBRC(23) = 9.00000/(1.00000+EXP(- ALGBRC(7)))+0.500000
      ALGBRC(22) = ALGBRC(18)
      ALGBRC(14) = ( (- STATES(1)/30.0000)*STATES(1))/30.0000
      ALGBRC(25) =  3.50000*EXP(ALGBRC(14))+1.50000
      ALGBRC(19) = 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(15))/CONSTS(16)))
      ALGBRC(24) = ALGBRC(19)/CONSTS(32)
      ALGBRC(26) = (1.00000 - ALGBRC(19))/CONSTS(32)
      ALGBRC(27) = 1.00000/(1.00000+CONSTS(23)/STATES(5) ** 3.00000)
      ALGBRC(35) = 10.0000+ 4954.00*EXP(STATES(1)/15.6000)
      ALGBRC(36) = CONSTS(31)/(1.00000+STATES(5)/CONSTS(24) ** 4.00000)+0.100000
      ALGBRC(34) = 1.00000 - 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(19))/CONSTS(20)))
      ALGBRC(37) =  (ALGBRC(35) - ALGBRC(36))*ALGBRC(34)+ALGBRC(36)
      ALGBRC(39) = 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(21))/CONSTS(22)))
      ALGBRC(40) = ( ALGBRC(27)*ALGBRC(39))/ALGBRC(37)
      ALGBRC(41) = (1.00000 - ALGBRC(39))/ALGBRC(37)
      ALGBRC(38) =  (ALGBRC(35) - 450.000)*ALGBRC(34)+450.000
      ALGBRC(42) = ALGBRC(39)/ALGBRC(38)
      ALGBRC(43) = (1.00000 - ALGBRC(39))/ALGBRC(38)
      ALGBRC(31) = 1.00000/(1.00000+EXP(- (STATES(1) - CONSTS(17))/CONSTS(18)))
      ALGBRC(32) = (1.00000 - ALGBRC(31))/CONSTS(33)
      ALGBRC(29) =  0.0241680*ALGBRC(27)
      ALGBRC(44) = ( (( (( ALGBRC(32)*ALGBRC(24))/ALGBRC(26))*ALGBRC(29))/CONSTS(25))*ALGBRC(41))/ALGBRC(40)
      ALGBRC(33) = ALGBRC(32)
      ALGBRC(45) = ( (( (( ALGBRC(33)*ALGBRC(24))/ALGBRC(26))*CONSTS(26))/CONSTS(27))*ALGBRC(43))/ALGBRC(42)
      ALGBRC(46) = (((((1.00000 - STATES(8)) - STATES(10)) - STATES(9)) - STATES(11)) - STATES(6)) - STATES(7)
      ALGBRC(28) =  0.0182688*ALGBRC(27)
      ALGBRC(30) = ( (( ALGBRC(28)*CONSTS(28))/CONSTS(29))*CONSTS(25))/ALGBRC(29)
      ALGBRC(54) = (STATES(1)+33.5000)/10.0000
      ALGBRC(60) = 1.00000/(1.00000+EXP(ALGBRC(54)))
      ALGBRC(56) = (STATES(1)+60.0000)/10.0000
      ALGBRC(63) = 3000.00/(1.00000+EXP(ALGBRC(56)))+30.0000
      ALGBRC(62) = ALGBRC(60)
      ALGBRC(58) = (STATES(1)+33.5000)/10.0000
      ALGBRC(65) = 20.0000/(1.00000+EXP(ALGBRC(58)))+20.0000
      ALGBRC(73) = ( CONSTS(59)*STATES(13)*STATES(13))/( STATES(13)*STATES(13)+ CONSTS(57)*CONSTS(57))
      ALGBRC(74) =  (( CONSTS(60)*STATES(22)*STATES(22))/( STATES(22)*STATES(22)+ CONSTS(58)*CONSTS(58)))*( STATES(22)*16.6670 - STATES(13))
      ALGBRC(80) = ( CONSTS(61)*CONSTS(62))/( (CONSTS(62)+STATES(13))*(CONSTS(62)+STATES(13)))
      ALGBRC(81) = ( CONSTS(63)*CONSTS(64))/( (CONSTS(64)+STATES(13))*(CONSTS(64)+STATES(13)))
      ALGBRC(82) = ( CONSTS(65)*CONSTS(66))/( (CONSTS(66)+STATES(13))*(CONSTS(66)+STATES(13)))
      ALGBRC(83) = ( CONSTS(67)*CONSTS(68))/( (CONSTS(68)+STATES(13))*(CONSTS(68)+STATES(13)))
      ALGBRC(84) = 1.00000/(1.00000+ALGBRC(80)+ALGBRC(81)+ALGBRC(82)+ALGBRC(83))
      ALGBRC(87) =  CONSTS(69)*STATES(13)*(CONSTS(71) - STATES(25)) -  CONSTS(70)*STATES(25)
      ALGBRC(85) = (STATES(21) - STATES(13))/CONSTS(72)
      ALGBRC(86) =  CONSTS(69)*STATES(21)*(CONSTS(71) - STATES(26)) -  CONSTS(70)*STATES(26)
      ALGBRC(88) = (- STATES(24)+ALGBRC(73)) - ALGBRC(74)
      ALGBRC(70) = TERNRY(STATES(23).GT.50.0000.AND.STATES(23).LT.CONSTS(48), (STATES(23) - 50.0000)/1.00000, TERNRY(STATES(23).GE.CONSTS(48),  CONSTS(54)*STATES(23)+CONSTS(79), 0.00000)
      ALGBRC(71) = ( STATES(22)*ALGBRC(70))/CONSTS(48)
      ALGBRC(89) = STATES(21)/1000.00
      ALGBRC(15) =  STATES(1)*2.00000*CONSTS(76)
      ALGBRC(90) = TERNRY(ABS(ALGBRC(15)).LT.0.00100000, ( 4.00000*CONSTS(14)*CONSTS(3)*CONSTS(76)*( ALGBRC(89)*EXP(ALGBRC(15)) -  0.341000*CONSTS(5)))/( 2.00000*CONSTS(76)), ( 4.00000*CONSTS(14)*STATES(1)*CONSTS(3)*CONSTS(76)*( ALGBRC(89)*EXP(ALGBRC(15)) -  0.341000*CONSTS(5)))/(EXP(ALGBRC(15)) - 1.00000))
      ALGBRC(72) = EXP( - CONSTS(53)*(STATES(1)+30.0000))/(1.00000+EXP( - CONSTS(53)*(STATES(1)+30.0000)))
      ALGBRC(92) =  (CONSTS(49)/1.00000)*ALGBRC(46)*ABS(ALGBRC(90))*ALGBRC(72)
      ALGBRC(93) = ( (( ALGBRC(46)*ALGBRC(71)*ABS(ALGBRC(90))*CONSTS(50))/1.00000)*EXP( - CONSTS(52)*(STATES(1)+30.0000)))/(1.00000+EXP( - CONSTS(52)*(STATES(1)+30.0000)))
      ALGBRC(95) =  ALGBRC(46)*CONSTS(51)*ABS(ALGBRC(90))
      ALGBRC(97) = ALGBRC(93)+ALGBRC(95)
      ALGBRC(91) =  CONSTS(13)*ALGBRC(46)*ALGBRC(90)
      ALGBRC(68) = 1.00000/(1.00000+CONSTS(43)/STATES(21) ** 3.00000)
      ALGBRC(96) =  STATES(20) ** 3.00000*CONSTS(5)*EXP( STATES(1)*0.350000*CONSTS(76)) -  CONSTS(6) ** 3.00000*ALGBRC(89)*EXP( STATES(1)*(0.350000 - 1.00000)*CONSTS(76))
      ALGBRC(67) = 1.00000+ 0.200000*EXP( STATES(1)*(0.350000 - 1.00000)*CONSTS(76))
      ALGBRC(98) =  CONSTS(44)*STATES(20) ** 3.00000+ CONSTS(45) ** 3.00000*ALGBRC(89)
      ALGBRC(99) =  CONSTS(46) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(89)/CONSTS(47))
      ALGBRC(69) =  CONSTS(47)*CONSTS(6) ** 3.00000*(1.00000+STATES(20)/CONSTS(46) ** 3.00000)
      ALGBRC(100) =  STATES(20) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(89)
      ALGBRC(101) = ALGBRC(98)+ALGBRC(99)+ALGBRC(69)+ALGBRC(100)
      ALGBRC(102) = ( CONSTS(42)*ALGBRC(68)*ALGBRC(96))/( ALGBRC(67)*ALGBRC(101))
      ALGBRC(75) = ( CONSTS(61)*CONSTS(62))/( (CONSTS(62)+STATES(21))*(CONSTS(62)+STATES(21)))
      ALGBRC(76) = ( CONSTS(63)*CONSTS(64))/( (CONSTS(64)+STATES(21))*(CONSTS(64)+STATES(21)))
      ALGBRC(77) = ( CONSTS(65)*CONSTS(66))/( (CONSTS(66)+STATES(21))*(CONSTS(66)+STATES(21)))
      ALGBRC(78) = ( CONSTS(67)*CONSTS(68))/( (CONSTS(68)+STATES(21))*(CONSTS(68)+STATES(21)))
      ALGBRC(79) = 1.00000/(1.00000+ALGBRC(75)+ALGBRC(76)+ALGBRC(77)+ALGBRC(78))
      ALGBRC(64) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(1)*CONSTS(76))+ 0.0365000*CONSTS(78)*EXP( - STATES(1)*CONSTS(76)))
      ALGBRC(66) = ( (( CONSTS(39)*ALGBRC(64)*STATES(20))/(STATES(20)+CONSTS(41)))*CONSTS(4))/(CONSTS(4)+CONSTS(40))
      ALGBRC(103) =  CONSTS(7)*ALGBRC(102)
      ALGBRC(106) =  (1.00000/CONSTS(76))*log(CONSTS(6)/STATES(20))
      ALGBRC(107) =  CONSTS(12)*STATES(3)*STATES(4)*STATES(2)*STATES(2)*STATES(2)*(STATES(1) - ALGBRC(106))
      ALGBRC(47) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - CONSTS(80)) - 59.2150)))
      ALGBRC(48) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(80))+5.47600))+ 1.00000*EXP( 0.0617500*((STATES(1) - CONSTS(80)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(80))+4.75300)))
      ALGBRC(49) = ALGBRC(47)/(ALGBRC(47)+ALGBRC(48))
      ALGBRC(50) =  CONSTS(34)* (CONSTS(4)/5.40000) ** (1.0 / 2)*ALGBRC(49)*(STATES(1) - CONSTS(80))
      ALGBRC(55) = 1.00000/(1.00000+EXP(ALGBRC(54)))
      ALGBRC(57) =  CONSTS(37)*STATES(16)*(STATES(17)+ 0.500000*ALGBRC(55))*(STATES(1) - CONSTS(80))
      ALGBRC(59) =  CONSTS(38)*STATES(18)*STATES(19)*(STATES(1) - CONSTS(80))
      ALGBRC(61) = ALGBRC(57)+ALGBRC(59)
      ALGBRC(94) =  2.00000*CONSTS(7)*ALGBRC(91)
      ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(1)+33.0000)/22.4000))
      ALGBRC(52) =  CONSTS(35)* (CONSTS(4)/5.40000) ** (1.0 / 2)*STATES(12)*ALGBRC(51)*(STATES(1) - CONSTS(80))
      ALGBRC(104) =  (1.00000/CONSTS(76))*log((CONSTS(4)+ CONSTS(75)*CONSTS(6))/(CONSTS(74)+ CONSTS(75)*STATES(20)))
      ALGBRC(53) = 1.00000+0.800000/(1.00000+0.500000/STATES(13) ** 3.00000)
      ALGBRC(105) =  CONSTS(36)*ALGBRC(53)*STATES(14)*STATES(15)*(STATES(1) - ALGBRC(104))
      ALGBRC(4) =  INT(VOI/CONSTS(9))*CONSTS(9)
      ALGBRC(11) = TERNRY(VOI - ALGBRC(4).GE.CONSTS(8).AND.VOI - ALGBRC(4).LE.CONSTS(8)+CONSTS(10), CONSTS(11), 0.00000)
      ALGBRC(108) = - (ALGBRC(107)+ALGBRC(50)+ALGBRC(52)+ALGBRC(105)+ALGBRC(61)+ALGBRC(103)+ALGBRC(94)+ALGBRC(66)+ALGBRC(11))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END