C
C There are a total of 79 entries in the algebraic variable array.
C There are a total of 11 entries in each of the rate and state variable arrays.
C There are a total of 75 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is Vm in component membrane (millivolt).
C CONSTS(1) is cm in component membrane (femtofarad).
C ALGBRC(4) is Ica in component Ica (femtoampere).
C ALGBRC(2) is Ik in component Ik (femtoampere).
C ALGBRC(5) is Ikca in component Ikca (femtoampere).
C ALGBRC(44) is Ikatp in component Ikatp (femtoampere).
C CONSTS(2) is gK in component Ik (picosiemens).
C CONSTS(3) is VK in component model_parameters (millivolt).
C STATES(2) is n in component n (dimensionless).
C ALGBRC(1) is n_infinity in component n (dimensionless).
C CONSTS(4) is tau_n in component n (millisecond).
C CONSTS(5) is vn in component n (millivolt).
C CONSTS(6) is sn in component n (millivolt).
C CONSTS(7) is gCa in component Ica (picosiemens).
C CONSTS(8) is VCa in component model_parameters (millivolt).
C ALGBRC(3) is m_infinity in component m (dimensionless).
C CONSTS(9) is v in component m (millivolt).
C CONSTS(10) is sm in component m (millivolt).
C CONSTS(11) is gkCa in component Ikca (picosiemens).
C CONSTS(12) is kd in component Ikca (micromolar).
C STATES(3) is c in component c (micromolar).
C CONSTS(13) is gkATP_ in component Ikatp (picosiemens).
C ALGBRC(42) is katpo in component Ikatp (dimensionless).
C ALGBRC(7) is topo in component Ikatp (dimensionless).
C ALGBRC(40) is bottomo in component Ikatp (dimensionless).
C ALGBRC(6) is mgadp in component Ikatp (micromolar).
C ALGBRC(8) is adp3m in component Ikatp (micromolar).
C ALGBRC(37) is atp4m in component Ikatp (micromolar).
C ALGBRC(36) is atp in component atp (micromolar).
C STATES(4) is adp in component adp (micromolar).
C ALGBRC(9) is JGPDH in component JGPDH (micromolar_millisecond).
C CONSTS(14) is kGPDH in component JGPDH (micromolar_millisecond).
C STATES(5) is FBP in component FBP (micromolar).
C STATES(6) is G6P in component G6P (micromolar).
C CONSTS(75) is JGK_ms in component JGK (micromolar_millisecond).
C ALGBRC(79) is JPFK_ms in component JPFK (micromolar_millisecond).
C ALGBRC(10) is F6P in component F6P (micromolar).
C ALGBRC(78) is JPFK in component JPFK (micromolar_second).
C CONSTS(15) is bottom1 in component JPFK (dimensionless).
C CONSTS(16) is topa1 in component JPFK (dimensionless).
C CONSTS(17) is k1 in component JPFK (micromolar).
C CONSTS(18) is k2 in component JPFK (micromolar).
C CONSTS(19) is k3 in component JPFK (micromolar).
C CONSTS(20) is k4 in component JPFK (micromolar).
C CONSTS(21) is VmaxPFK in component JPFK (micromolar_millisecond).
C ALGBRC(38) is weight2 in component JPFK (dimensionless).
C CONSTS(72) is topa2 in component JPFK (dimensionless).
C ALGBRC(41) is bottom2 in component JPFK (dimensionless).
C ALGBRC(14) is topa3 in component JPFK (dimensionless).
C ALGBRC(12) is weight3 in component JPFK (dimensionless).
C ALGBRC(43) is bottom3 in component JPFK (dimensionless).
C CONSTS(22) is f13 in component JPFK (dimensionless).
C CONSTS(23) is f43 in component JPFK (dimensionless).
C CONSTS(24) is f23 in component JPFK (dimensionless).
C CONSTS(25) is f42 in component JPFK (dimensionless).
C CONSTS(26) is f41 in component JPFK (dimensionless).
C ALGBRC(45) is weight4 in component JPFK (dimensionless).
C ALGBRC(46) is topa4 in component JPFK (dimensionless).
C ALGBRC(47) is bottom4 in component JPFK (dimensionless).
C ALGBRC(16) is weight5 in component JPFK (dimensionless).
C ALGBRC(48) is topa5 in component JPFK (dimensionless).
C ALGBRC(49) is bottom5 in component JPFK (dimensionless).
C ALGBRC(50) is weight6 in component JPFK (dimensionless).
C ALGBRC(51) is topa6 in component JPFK (dimensionless).
C ALGBRC(52) is bottom6 in component JPFK (dimensionless).
C ALGBRC(17) is weight7 in component JPFK (dimensionless).
C ALGBRC(53) is topa7 in component JPFK (dimensionless).
C ALGBRC(54) is bottom7 in component JPFK (dimensionless).
C ALGBRC(55) is weight8 in component JPFK (dimensionless).
C ALGBRC(56) is topa8 in component JPFK (dimensionless).
C ALGBRC(57) is bottom8 in component JPFK (dimensionless).
C CONSTS(74) is weight9 in component JPFK (dimensionless).
C ALGBRC(58) is topa9 in component JPFK (dimensionless).
C ALGBRC(59) is bottom9 in component JPFK (dimensionless).
C ALGBRC(60) is weight10 in component JPFK (dimensionless).
C ALGBRC(61) is topa10 in component JPFK (dimensionless).
C ALGBRC(62) is bottom10 in component JPFK (dimensionless).
C ALGBRC(18) is weight11 in component JPFK (dimensionless).
C ALGBRC(63) is topa11 in component JPFK (dimensionless).
C ALGBRC(64) is bottom11 in component JPFK (dimensionless).
C ALGBRC(65) is weight12 in component JPFK (dimensionless).
C ALGBRC(66) is topa12 in component JPFK (dimensionless).
C ALGBRC(67) is bottom12 in component JPFK (dimensionless).
C ALGBRC(19) is weight13 in component JPFK (dimensionless).
C ALGBRC(68) is topa13 in component JPFK (dimensionless).
C ALGBRC(69) is bottom13 in component JPFK (dimensionless).
C ALGBRC(70) is weight14 in component JPFK (dimensionless).
C ALGBRC(71) is topa14 in component JPFK (dimensionless).
C ALGBRC(72) is bottom14 in component JPFK (dimensionless).
C ALGBRC(20) is weight15 in component JPFK (dimensionless).
C ALGBRC(73) is topa15 in component JPFK (dimensionless).
C ALGBRC(74) is bottom15 in component JPFK (dimensionless).
C ALGBRC(75) is weight16 in component JPFK (dimensionless).
C ALGBRC(76) is topa16 in component JPFK (dimensionless).
C ALGBRC(77) is bottom16 in component JPFK (dimensionless).
C ALGBRC(21) is topb in component JPFK (dimensionless).
C CONSTS(27) is AMP in component model_parameters (micromolar).
C CONSTS(28) is lambda in component JPFK (dimensionless).
C CONSTS(29) is kappa in component model_parameters (dimensionless).
C ALGBRC(15) is JPDH in component JPDH (micromolar_millisecond).
C CONSTS(30) is p1 in component JPDH (dimensionless).
C CONSTS(31) is p2 in component JPDH (dimensionless).
C CONSTS(32) is p3 in component JPDH (micromolar).
C CONSTS(33) is JGPDHbas in component JPDH (micromolar_millisecond).
C ALGBRC(13) is NADm in component NADm (millimolar).
C STATES(7) is Cam in component Cam (micromolar).
C STATES(8) is NADHm in component NADHm (millimolar).
C CONSTS(34) is gamma in component model_parameters (dimensionless).
C ALGBRC(11) is JO in component JO (micromolar_millisecond).
C CONSTS(35) is p4 in component JO (micromolar_millisecond).
C CONSTS(36) is p5 in component JO (millimolar).
C CONSTS(37) is p6 in component JO (millivolt).
C CONSTS(38) is p7 in component JO (millivolt).
C STATES(9) is delta_psi in component delta_psi (millivolt).
C CONSTS(39) is NADmtot in component NADm (millimolar).
C CONSTS(40) is Cmito in component delta_psi (micromolar_millivolt).
C ALGBRC(22) is JHres in component JHres (micromolar_millisecond).
C ALGBRC(30) is JHatp in component JHatp (micromolar_millisecond).
C ALGBRC(34) is JANT in component JANT (micromolar_millisecond).
C ALGBRC(23) is JHleak in component JHleak (micromolar_millisecond).
C ALGBRC(25) is JNaCa in component JNaCa (micromolar_millisecond).
C ALGBRC(24) is Juni in component Juni (micromolar_millisecond).
C CONSTS(41) is p8 in component JHres (micromolar_millisecond).
C CONSTS(42) is p9 in component JHres (millimolar).
C CONSTS(43) is p10 in component JHres (millivolt).
C CONSTS(44) is p11 in component JHres (millivolt).
C ALGBRC(28) is JF1F0 in component JF1F0 (micromolar_millisecond).
C CONSTS(45) is p13 in component JF1F0 (millimolar).
C CONSTS(46) is p14 in component JF1F0 (millivolt).
C CONSTS(47) is p15 in component JF1F0 (millivolt).
C CONSTS(48) is p16 in component JF1F0 (micromolar_millisecond).
C ALGBRC(26) is ATPm in component ATPm (millimolar).
C CONSTS(49) is JGK in component JGK (micromolar_second).
C CONSTS(50) is p17 in component JHleak (micromolar_millisecond_millivolt).
C CONSTS(51) is p18 in component JHleak (micromolar_millisecond).
C CONSTS(52) is p19 in component JANT (micromolar_millisecond).
C CONSTS(53) is p20 in component JANT (dimensionless).
C CONSTS(54) is FRT in component JANT (per_millivolt).
C ALGBRC(32) is RATm in component RATm (dimensionless).
C STATES(10) is ADPm in component ADPm (millimolar).
C CONSTS(55) is p21 in component Juni (per_micromolar_millisecond_millivolt).
C CONSTS(56) is p22 in component Juni (second_order_rate_constant).
C CONSTS(57) is p23 in component JNaCa (micromolar_millisecond).
C CONSTS(58) is p24 in component JNaCa (per_millivolt).
C CONSTS(59) is fmito in component Cam (dimensionless).
C ALGBRC(27) is Jmito in component Jmito (micromolar_millisecond).
C CONSTS(60) is Amtot in component ATPm (millimolar).
C CONSTS(73) is delta in component model_parameters (dimensionless).
C ALGBRC(39) is Jhyd in component Jhyd (micromolar_millisecond).
C CONSTS(61) is khyd in component Jhyd (second_order_rate_constant).
C CONSTS(62) is khydbas in component Jhyd (first_order_rate_constant).
C CONSTS(63) is atot in component atp (micromolar).
C CONSTS(64) is fcyt in component c (dimensionless).
C ALGBRC(35) is Jer in component Jer (micromolar_millisecond).
C ALGBRC(29) is Jmem in component Jmem (micromolar_millisecond).
C CONSTS(65) is kPMCA in component Jmem (first_order_rate_constant).
C CONSTS(66) is alpha in component Jmem (micromolar_millisecond_femtoampere).
C CONSTS(67) is Cbas in component Jmem (micromolar).
C ALGBRC(31) is Jleak in component Jleak (micromolar_millisecond).
C CONSTS(68) is pleak in component Jleak (first_order_rate_constant).
C STATES(11) is Caer in component Caer (micromolar).
C ALGBRC(33) is JSERCA in component JSERCA (micromolar_millisecond).
C CONSTS(69) is kSERCA in component JSERCA (first_order_rate_constant).
C CONSTS(70) is fer in component Caer (dimensionless).
C CONSTS(71) is Vc_Ver in component Caer (dimensionless).
C RATES(1) is d/dt Vm in component membrane (millivolt).
C RATES(2) is d/dt n in component n (dimensionless).
C RATES(6) is d/dt G6P in component G6P (micromolar).
C RATES(5) is d/dt FBP in component FBP (micromolar).
C RATES(8) is d/dt NADHm in component NADHm (millimolar).
C RATES(9) is d/dt delta_psi in component delta_psi (millivolt).
C RATES(7) is d/dt Cam in component Cam (micromolar).
C RATES(10) is d/dt ADPm in component ADPm (millimolar).
C RATES(4) is d/dt adp in component adp (micromolar).
C RATES(3) is d/dt c in component c (micromolar).
C RATES(11) is d/dt Caer in component Caer (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -60.0
      CONSTS(1) = 5300.0
      CONSTS(2) = 2700.0
      CONSTS(3) = -75.0
      STATES(2) = 0.0
      CONSTS(4) = 20.0
      CONSTS(5) = -16.0
      CONSTS(6) = 5.0
      CONSTS(7) = 1000.0
      CONSTS(8) = 25.0
      CONSTS(9) = -20.0
      CONSTS(10) = 12.0
      CONSTS(11) = 300.0
      CONSTS(12) = 0.5
      STATES(3) = 0.17
      CONSTS(13) = 16000.0
      STATES(4) = 1137.0
      CONSTS(14) = 0.0005
      STATES(5) = 2.16
      STATES(6) = 301.0
      CONSTS(15) = 1
      CONSTS(16) = 0
      CONSTS(17) = 30
      CONSTS(18) = 1
      CONSTS(19) = 50000
      CONSTS(20) = 1000
      CONSTS(21) = 5.0
      CONSTS(22) = 0.02
      CONSTS(23) = 20
      CONSTS(24) = 0.2
      CONSTS(25) = 20
      CONSTS(26) = 20
      CONSTS(27) = 500.0
      CONSTS(28) = 0.06
      CONSTS(29) = 0.001
      CONSTS(30) = 400.0
      CONSTS(31) = 1.0
      CONSTS(32) = 0.01
      CONSTS(33) = 0.0005
      STATES(7) = 0.2
      STATES(8) = 0.4
      CONSTS(34) = 0.001
      CONSTS(35) = 0.6
      CONSTS(36) = 0.1
      CONSTS(37) = 177.0
      CONSTS(38) = 5.0
      STATES(9) = 164.0
      CONSTS(39) = 10.0
      CONSTS(40) = 1.8
      CONSTS(41) = 7.0
      CONSTS(42) = 0.1
      CONSTS(43) = 177.0
      CONSTS(44) = 5.0
      CONSTS(45) = 10.0
      CONSTS(46) = 190.0
      CONSTS(47) = 8.5
      CONSTS(48) = 35.0
      CONSTS(49) = 0.4
      CONSTS(50) = 0.002
      CONSTS(51) = -0.03
      CONSTS(52) = 0.35
      CONSTS(53) = 2.0
      CONSTS(54) = 0.037410133
      STATES(10) = 11.1
      CONSTS(55) = 0.04
      CONSTS(56) = 1.1
      CONSTS(57) = 0.01
      CONSTS(58) = 0.016
      CONSTS(59) = 0.01
      CONSTS(60) = 15.0
      CONSTS(61) = 0.00005
      CONSTS(62) = 0.00005
      CONSTS(63) = 2500.0
      CONSTS(64) = 0.01
      CONSTS(65) = 0.1
      CONSTS(66) = 4.5E-6
      CONSTS(67) = 0.05
      CONSTS(68) = 0.0002
      STATES(11) = 345.0
      CONSTS(69) = 0.4
      CONSTS(70) = 0.01
      CONSTS(71) = 31.0
      CONSTS(72) = CONSTS(16)
      CONSTS(73) = 3.90000/53.2000
      CONSTS(74) = CONSTS(27)/CONSTS(17)
      CONSTS(75) =  CONSTS(29)*CONSTS(49)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = 1.00000/(1.00000+EXP((CONSTS(5) - STATES(1))/CONSTS(6)))
      RATES(2) = (ALGBRC(1) - STATES(2))/CONSTS(4)
      ALGBRC(9) =  CONSTS(14)* ABS(STATES(5)/1.00000) ** (1.0 / 2)
      ALGBRC(13) = CONSTS(39) - STATES(8)
      ALGBRC(15) =  (CONSTS(30)/(CONSTS(31)+STATES(8)/ALGBRC(13)))*(STATES(7)/(CONSTS(32)+STATES(7)))*(ALGBRC(9)+CONSTS(33))
      ALGBRC(11) = ( CONSTS(35)*(STATES(8)/(CONSTS(36)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(37))/CONSTS(38)))
      RATES(8) =  CONSTS(34)*(ALGBRC(15) - ALGBRC(11))
      ALGBRC(25) =  (CONSTS(57)/STATES(3))*STATES(7)*EXP( CONSTS(58)*STATES(9))
      ALGBRC(24) =  ( CONSTS(55)*STATES(9) - CONSTS(56))*STATES(3) ** 2.00000
      ALGBRC(27) = ALGBRC(25) - ALGBRC(24)
      RATES(7) =  - CONSTS(59)*ALGBRC(27)
      ALGBRC(22) = ( CONSTS(41)*(STATES(8)/(CONSTS(42)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(43))/CONSTS(44)))
      ALGBRC(26) = CONSTS(60) - STATES(10)
      ALGBRC(28) = (( CONSTS(48)*CONSTS(45))/(CONSTS(45)+ALGBRC(26)))/(1.00000+EXP((CONSTS(46) - STATES(9))/CONSTS(47)))
      ALGBRC(30) =  3.00000*ALGBRC(28)
      ALGBRC(32) = ALGBRC(26)/STATES(10)
      ALGBRC(34) =  CONSTS(52)*((ALGBRC(32)/(ALGBRC(32)+CONSTS(53)))/EXP( -0.500000*CONSTS(54)*STATES(9)))
      ALGBRC(23) =  CONSTS(50)*STATES(9)+CONSTS(51)
      RATES(9) = (ALGBRC(22) - (ALGBRC(30)+ALGBRC(34)+ALGBRC(23)+ALGBRC(25)+ 2.00000*ALGBRC(24)))/CONSTS(40)
      RATES(10) =  CONSTS(34)*(ALGBRC(34) - ALGBRC(28))
      ALGBRC(31) =  CONSTS(68)*(STATES(11) - STATES(3))
      ALGBRC(33) =  CONSTS(69)*STATES(3)
      ALGBRC(35) = ALGBRC(31) - ALGBRC(33)
      ALGBRC(3) = 1.00000/(1.00000+EXP((CONSTS(9) - STATES(1))/CONSTS(10)))
      ALGBRC(4) =  CONSTS(7)*ALGBRC(3)*(STATES(1) - CONSTS(8))
      ALGBRC(29) = - ( CONSTS(66)*ALGBRC(4)+ CONSTS(65)*(STATES(3) - CONSTS(67)))
      RATES(3) =  CONSTS(64)*(ALGBRC(29)+ALGBRC(35)+ CONSTS(73)*ALGBRC(27))
      RATES(11) =  - CONSTS(70)*CONSTS(71)*ALGBRC(35)
      ALGBRC(36) = CONSTS(63) - STATES(4)
      ALGBRC(39) =  ( CONSTS(61)*STATES(3)+CONSTS(62))*ALGBRC(36)
      RATES(4) =  - CONSTS(73)*ALGBRC(34)+ALGBRC(39)
      ALGBRC(2) =  CONSTS(2)*STATES(2)*(STATES(1) - CONSTS(3))
      ALGBRC(5) =  (CONSTS(11)/(1.00000+CONSTS(12)/STATES(3) ** 2.00000))*(STATES(1) - CONSTS(3))
      ALGBRC(6) =  0.165000*STATES(4)
      ALGBRC(7) =  0.0800000*(1.00000+( 2.00000*ALGBRC(6))/17.0000)+ 0.890000*ALGBRC(6)/17.0000 ** 2.00000
      ALGBRC(8) =  0.135000*STATES(4)
      ALGBRC(37) =  0.0500000*ALGBRC(36)
      ALGBRC(40) =  1.00000+ALGBRC(6)/17.0000 ** 2.00000*(1.00000+ALGBRC(8)/26.0000+ALGBRC(37)/1.00000)
      ALGBRC(42) = ALGBRC(7)/ALGBRC(40)
      ALGBRC(44) =  CONSTS(13)*ALGBRC(42)*(STATES(1) - CONSTS(3))
      RATES(1) = - (ALGBRC(2)+ALGBRC(4)+ALGBRC(5)+ALGBRC(44))/CONSTS(1)
      ALGBRC(10) =  0.300000*STATES(6)
      ALGBRC(12) = ALGBRC(10) ** 2.00000/( CONSTS(19)*1.00000)
      ALGBRC(14) = CONSTS(72)+ALGBRC(12)
      ALGBRC(45) =  ALGBRC(10)*ALGBRC(36) ** 2.00000/( CONSTS(23)*CONSTS(19)*CONSTS(20)*1.00000 ** 2.00000)
      ALGBRC(46) = ALGBRC(14)+ALGBRC(45)
      ALGBRC(48) = ALGBRC(46)
      ALGBRC(51) = ALGBRC(48)
      ALGBRC(17) = ( STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(24)*1.00000)
      ALGBRC(53) = ALGBRC(51)+ALGBRC(17)
      ALGBRC(55) = ( STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(25)*CONSTS(23)*1.00000 ** 2.00000)
      ALGBRC(56) = ALGBRC(53)+ALGBRC(55)
      ALGBRC(58) = ALGBRC(56)
      ALGBRC(61) = ALGBRC(58)
      ALGBRC(18) = ( CONSTS(27)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(22)*1.00000)
      ALGBRC(63) = ALGBRC(61)+ALGBRC(18)
      ALGBRC(65) = ( CONSTS(27)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(20)*CONSTS(22)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
      ALGBRC(66) = ALGBRC(63)+ALGBRC(65)
      ALGBRC(68) = ALGBRC(66)
      ALGBRC(71) = ALGBRC(68)
      ALGBRC(73) = ALGBRC(71)
      ALGBRC(75) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(22)*CONSTS(25)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
      ALGBRC(76) = ALGBRC(73)+ALGBRC(75)
      ALGBRC(38) = ALGBRC(36) ** 2.00000/( CONSTS(20)*1.00000)
      ALGBRC(41) = CONSTS(15)+ALGBRC(38)
      ALGBRC(43) = ALGBRC(41)+ALGBRC(12)
      ALGBRC(47) = ALGBRC(43)+ALGBRC(45)
      ALGBRC(16) = STATES(5)/CONSTS(18)
      ALGBRC(49) = ALGBRC(47)+ALGBRC(16)
      ALGBRC(50) = ( STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(20)*CONSTS(25)*1.00000)
      ALGBRC(52) = ALGBRC(49)+ALGBRC(50)
      ALGBRC(54) = ALGBRC(52)+ALGBRC(17)
      ALGBRC(57) = ALGBRC(54)+ALGBRC(55)
      ALGBRC(59) = ALGBRC(57)+CONSTS(74)
      ALGBRC(60) = ( CONSTS(27)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(20)*CONSTS(26)*1.00000)
      ALGBRC(62) = ALGBRC(59)+ALGBRC(60)
      ALGBRC(64) = ALGBRC(62)+ALGBRC(18)
      ALGBRC(67) = ALGBRC(64)+ALGBRC(65)
      ALGBRC(19) = ( CONSTS(27)*STATES(5))/( CONSTS(17)*CONSTS(18))
      ALGBRC(69) = ALGBRC(67)+ALGBRC(19)
      ALGBRC(70) = ( CONSTS(27)*STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(20)*CONSTS(25)*CONSTS(26)*1.00000)
      ALGBRC(72) = ALGBRC(69)+ALGBRC(70)
      ALGBRC(20) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(24)*CONSTS(22)*1.00000)
      ALGBRC(74) = ALGBRC(72)+ALGBRC(20)
      ALGBRC(77) = ALGBRC(74)+ALGBRC(75)
      ALGBRC(21) = ALGBRC(20)
      ALGBRC(78) = ( CONSTS(28)*CONSTS(21)*ALGBRC(76)+ CONSTS(21)*ALGBRC(21))/ALGBRC(77)
      ALGBRC(79) =  CONSTS(29)*ALGBRC(78)
      RATES(6) = CONSTS(75) - ALGBRC(79)
      RATES(5) = ALGBRC(79) -  0.500000*ALGBRC(9)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = 1.00000/(1.00000+EXP((CONSTS(5) - STATES(1))/CONSTS(6)))
      ALGBRC(9) =  CONSTS(14)* ABS(STATES(5)/1.00000) ** (1.0 / 2)
      ALGBRC(13) = CONSTS(39) - STATES(8)
      ALGBRC(15) =  (CONSTS(30)/(CONSTS(31)+STATES(8)/ALGBRC(13)))*(STATES(7)/(CONSTS(32)+STATES(7)))*(ALGBRC(9)+CONSTS(33))
      ALGBRC(11) = ( CONSTS(35)*(STATES(8)/(CONSTS(36)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(37))/CONSTS(38)))
      ALGBRC(25) =  (CONSTS(57)/STATES(3))*STATES(7)*EXP( CONSTS(58)*STATES(9))
      ALGBRC(24) =  ( CONSTS(55)*STATES(9) - CONSTS(56))*STATES(3) ** 2.00000
      ALGBRC(27) = ALGBRC(25) - ALGBRC(24)
      ALGBRC(22) = ( CONSTS(41)*(STATES(8)/(CONSTS(42)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(43))/CONSTS(44)))
      ALGBRC(26) = CONSTS(60) - STATES(10)
      ALGBRC(28) = (( CONSTS(48)*CONSTS(45))/(CONSTS(45)+ALGBRC(26)))/(1.00000+EXP((CONSTS(46) - STATES(9))/CONSTS(47)))
      ALGBRC(30) =  3.00000*ALGBRC(28)
      ALGBRC(32) = ALGBRC(26)/STATES(10)
      ALGBRC(34) =  CONSTS(52)*((ALGBRC(32)/(ALGBRC(32)+CONSTS(53)))/EXP( -0.500000*CONSTS(54)*STATES(9)))
      ALGBRC(23) =  CONSTS(50)*STATES(9)+CONSTS(51)
      ALGBRC(31) =  CONSTS(68)*(STATES(11) - STATES(3))
      ALGBRC(33) =  CONSTS(69)*STATES(3)
      ALGBRC(35) = ALGBRC(31) - ALGBRC(33)
      ALGBRC(3) = 1.00000/(1.00000+EXP((CONSTS(9) - STATES(1))/CONSTS(10)))
      ALGBRC(4) =  CONSTS(7)*ALGBRC(3)*(STATES(1) - CONSTS(8))
      ALGBRC(29) = - ( CONSTS(66)*ALGBRC(4)+ CONSTS(65)*(STATES(3) - CONSTS(67)))
      ALGBRC(36) = CONSTS(63) - STATES(4)
      ALGBRC(39) =  ( CONSTS(61)*STATES(3)+CONSTS(62))*ALGBRC(36)
      ALGBRC(2) =  CONSTS(2)*STATES(2)*(STATES(1) - CONSTS(3))
      ALGBRC(5) =  (CONSTS(11)/(1.00000+CONSTS(12)/STATES(3) ** 2.00000))*(STATES(1) - CONSTS(3))
      ALGBRC(6) =  0.165000*STATES(4)
      ALGBRC(7) =  0.0800000*(1.00000+( 2.00000*ALGBRC(6))/17.0000)+ 0.890000*ALGBRC(6)/17.0000 ** 2.00000
      ALGBRC(8) =  0.135000*STATES(4)
      ALGBRC(37) =  0.0500000*ALGBRC(36)
      ALGBRC(40) =  1.00000+ALGBRC(6)/17.0000 ** 2.00000*(1.00000+ALGBRC(8)/26.0000+ALGBRC(37)/1.00000)
      ALGBRC(42) = ALGBRC(7)/ALGBRC(40)
      ALGBRC(44) =  CONSTS(13)*ALGBRC(42)*(STATES(1) - CONSTS(3))
      ALGBRC(10) =  0.300000*STATES(6)
      ALGBRC(12) = ALGBRC(10) ** 2.00000/( CONSTS(19)*1.00000)
      ALGBRC(14) = CONSTS(72)+ALGBRC(12)
      ALGBRC(45) =  ALGBRC(10)*ALGBRC(36) ** 2.00000/( CONSTS(23)*CONSTS(19)*CONSTS(20)*1.00000 ** 2.00000)
      ALGBRC(46) = ALGBRC(14)+ALGBRC(45)
      ALGBRC(48) = ALGBRC(46)
      ALGBRC(51) = ALGBRC(48)
      ALGBRC(17) = ( STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(24)*1.00000)
      ALGBRC(53) = ALGBRC(51)+ALGBRC(17)
      ALGBRC(55) = ( STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(25)*CONSTS(23)*1.00000 ** 2.00000)
      ALGBRC(56) = ALGBRC(53)+ALGBRC(55)
      ALGBRC(58) = ALGBRC(56)
      ALGBRC(61) = ALGBRC(58)
      ALGBRC(18) = ( CONSTS(27)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(22)*1.00000)
      ALGBRC(63) = ALGBRC(61)+ALGBRC(18)
      ALGBRC(65) = ( CONSTS(27)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(20)*CONSTS(22)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
      ALGBRC(66) = ALGBRC(63)+ALGBRC(65)
      ALGBRC(68) = ALGBRC(66)
      ALGBRC(71) = ALGBRC(68)
      ALGBRC(73) = ALGBRC(71)
      ALGBRC(75) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(22)*CONSTS(25)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
      ALGBRC(76) = ALGBRC(73)+ALGBRC(75)
      ALGBRC(38) = ALGBRC(36) ** 2.00000/( CONSTS(20)*1.00000)
      ALGBRC(41) = CONSTS(15)+ALGBRC(38)
      ALGBRC(43) = ALGBRC(41)+ALGBRC(12)
      ALGBRC(47) = ALGBRC(43)+ALGBRC(45)
      ALGBRC(16) = STATES(5)/CONSTS(18)
      ALGBRC(49) = ALGBRC(47)+ALGBRC(16)
      ALGBRC(50) = ( STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(20)*CONSTS(25)*1.00000)
      ALGBRC(52) = ALGBRC(49)+ALGBRC(50)
      ALGBRC(54) = ALGBRC(52)+ALGBRC(17)
      ALGBRC(57) = ALGBRC(54)+ALGBRC(55)
      ALGBRC(59) = ALGBRC(57)+CONSTS(74)
      ALGBRC(60) = ( CONSTS(27)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(20)*CONSTS(26)*1.00000)
      ALGBRC(62) = ALGBRC(59)+ALGBRC(60)
      ALGBRC(64) = ALGBRC(62)+ALGBRC(18)
      ALGBRC(67) = ALGBRC(64)+ALGBRC(65)
      ALGBRC(19) = ( CONSTS(27)*STATES(5))/( CONSTS(17)*CONSTS(18))
      ALGBRC(69) = ALGBRC(67)+ALGBRC(19)
      ALGBRC(70) = ( CONSTS(27)*STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(20)*CONSTS(25)*CONSTS(26)*1.00000)
      ALGBRC(72) = ALGBRC(69)+ALGBRC(70)
      ALGBRC(20) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(24)*CONSTS(22)*1.00000)
      ALGBRC(74) = ALGBRC(72)+ALGBRC(20)
      ALGBRC(77) = ALGBRC(74)+ALGBRC(75)
      ALGBRC(21) = ALGBRC(20)
      ALGBRC(78) = ( CONSTS(28)*CONSTS(21)*ALGBRC(76)+ CONSTS(21)*ALGBRC(21))/ALGBRC(77)
      ALGBRC(79) =  CONSTS(29)*ALGBRC(78)
      RETURN
      END