C
C There are a total of 21 entries in the algebraic variable array.
C There are a total of 21 entries in each of the rate and state variable arrays.
C There are a total of 116 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is GLY in component GLY (millimolar).
C ALGBRC(4) is flux_GP in component flux_GP (flux).
C STATES(2) is G1P in component G1P (millimolar).
C ALGBRC(5) is V_PGLM in component V_PGLM (flux).
C STATES(3) is G6P in component G6P (millimolar).
C ALGBRC(6) is V_PGI in component V_PGI (flux).
C STATES(4) is F6P in component F6P (millimolar).
C ALGBRC(10) is V_PFK in component V_PFK (flux).
C STATES(5) is FBP in component FBP (millimolar).
C ALGBRC(11) is V_ALD in component V_ALD (flux).
C STATES(6) is DHAP in component DHAP (millimolar).
C ALGBRC(12) is V_TPI in component V_TPI (flux).
C STATES(7) is GAP in component GAP (millimolar).
C ALGBRC(14) is V_GAPDH in component V_GAPDH (flux).
C STATES(8) is Thirteen_BPG in component Thirteen_BPG (millimolar).
C ALGBRC(15) is V_PGK in component V_PGK (flux).
C STATES(9) is three_PG in component three_PG (millimolar).
C ALGBRC(16) is V_PGM in component V_PGM (flux).
C STATES(10) is two_PG in component two_PG (millimolar).
C ALGBRC(17) is V_ENOL in component V_ENOL (flux).
C STATES(11) is PEP in component PEP (millimolar).
C ALGBRC(18) is V_PK in component V_PK (flux).
C STATES(12) is PYR in component PYR (millimolar).
C ALGBRC(19) is V_LDH in component V_LDH (flux).
C STATES(13) is LAC in component LAC (millimolar).
C ALGBRC(1) is output in component LAC (flux).
C STATES(14) is Pi in component Pi (millimolar).
C CONSTS(1) is V_ATPase in component V_ATPase (flux).
C STATES(15) is ADP in component ADP (millimolar).
C ALGBRC(21) is V_ADK in component V_ADK (flux).
C ALGBRC(20) is V_CK in component V_CK (flux).
C STATES(16) is ATP in component ATP (millimolar).
C STATES(17) is AMP in component AMP (millimolar).
C STATES(18) is PCr in component PCr (millimolar).
C STATES(19) is Cr in component Cr (millimolar).
C STATES(20) is NADH in component NADH (millimolar).
C STATES(21) is NAD in component NAD (millimolar).
C CONSTS(2) is frac_a in component flux_GP (dimensionless).
C CONSTS(3) is frac_b in component flux_GP (dimensionless).
C ALGBRC(2) is V_GPa in component V_GPa (flux).
C ALGBRC(3) is V_GPb in component V_GPb (flux).
C CONSTS(4) is Ki_GLY_f in component V_GPa (millimolar).
C CONSTS(5) is Ki_GLY_b in component V_GPa (millimolar).
C CONSTS(6) is K_GLY_f in component V_GPa (millimolar).
C CONSTS(7) is K_GLY_b in component V_GPa (millimolar).
C CONSTS(8) is Keq_GP in component V_GPa (dimensionless).
C CONSTS(9) is Ki_Pi in component V_GPa (millimolar).
C CONSTS(10) is K_Pi in component V_GPa (millimolar).
C CONSTS(11) is Ki_G1P in component V_GPa (millimolar).
C CONSTS(101) is V_max_r in component V_GPa (flux).
C CONSTS(12) is V_max_f in component V_GPa (flux).
C CONSTS(13) is Ki_GLY_f in component V_GPb (millimolar).
C CONSTS(14) is Ki_GLY_b in component V_GPb (millimolar).
C CONSTS(15) is K_GLY_f in component V_GPb (millimolar).
C CONSTS(16) is Keq_GP in component V_GPb (dimensionless).
C CONSTS(17) is nH in component V_GPb (dimensionless).
C CONSTS(18) is K_AMP in component V_GPb (millimolar_1_75).
C CONSTS(19) is K_GLY_b in component V_GPb (millimolar).
C CONSTS(20) is Ki_Pi in component V_GPb (millimolar).
C CONSTS(21) is K_Pi in component V_GPb (millimolar).
C CONSTS(22) is Ki_G1P in component V_GPb (millimolar).
C CONSTS(23) is K_G1P in component V_GPb (millimolar).
C CONSTS(102) is V_max_r in component V_GPb (flux).
C CONSTS(24) is V_max_f in component V_GPb (flux).
C CONSTS(25) is K_G1P in component V_PGLM (millimolar).
C CONSTS(26) is K_G6P in component V_PGLM (millimolar).
C CONSTS(27) is Keq_PGLM in component V_PGLM (dimensionless).
C CONSTS(103) is V_max_r in component V_PGLM (flux).
C CONSTS(28) is V_max_f in component V_PGLM (flux).
C CONSTS(29) is K_G6P in component V_PGI (millimolar).
C CONSTS(30) is K_F6P in component V_PGI (millimolar).
C CONSTS(31) is Keq_PGI in component V_PGI (dimensionless).
C CONSTS(32) is V_max_r in component V_PGI (flux).
C CONSTS(104) is V_max_f in component V_PGI (flux).
C CONSTS(33) is K_FBP in component V_PFK (millimolar).
C CONSTS(34) is K_FBP_ in component V_PFK (millimolar).
C CONSTS(35) is K_F6P in component V_PFK (millimolar).
C CONSTS(36) is K_ATP in component V_PFK (millimolar).
C CONSTS(37) is K_F6P_ in component V_PFK (millimolar).
C CONSTS(38) is K_ATP_ in component V_PFK (millimolar).
C CONSTS(39) is K_ADP in component V_PFK (millimolar).
C CONSTS(40) is K_ADP_ in component V_PFK (millimolar).
C CONSTS(41) is Ki_ATP in component V_PFK (millimolar).
C CONSTS(42) is Ka_AMP in component V_PFK (millimolar).
C CONSTS(105) is V_max_r in component V_PFK (flux).
C CONSTS(43) is V_max_f in component V_PFK (flux).
C CONSTS(106) is alpha in component V_PFK (dimensionless).
C ALGBRC(9) is L in component V_PFK (dimensionless).
C CONSTS(44) is Lo in component V_PFK (dimensionless).
C ALGBRC(7) is delta in component V_PFK (dimensionless).
C ALGBRC(8) is delta_ in component V_PFK (dimensionless).
C CONSTS(45) is d in component V_PFK (dimensionless).
C CONSTS(46) is var_e in component V_PFK (dimensionless).
C CONSTS(47) is K_GAP in component V_ALD (millimolar).
C CONSTS(48) is K_FBP in component V_ALD (millimolar).
C CONSTS(49) is K_DHAP in component V_ALD (millimolar).
C CONSTS(50) is Keq_ALD in component V_ALD (millimolar).
C CONSTS(107) is V_max_r in component V_ALD (flux).
C CONSTS(51) is V_max_f in component V_ALD (flux).
C CONSTS(52) is K_GAP in component V_TPI (millimolar).
C CONSTS(53) is K_DHAP in component V_TPI (millimolar).
C CONSTS(54) is Keq_TPI in component V_TPI (dimensionless).
C CONSTS(108) is V_max_r in component V_TPI (flux).
C CONSTS(55) is V_max_f in component V_TPI (flux).
C ALGBRC(13) is D_GAPDH in component V_GAPDH (dimensionless).
C CONSTS(56) is K_GAP in component V_GAPDH (millimolar).
C CONSTS(57) is K_Pi in component V_GAPDH (millimolar).
C CONSTS(58) is K_NAD in component V_GAPDH (millimolar).
C CONSTS(59) is K_NADH in component V_GAPDH (millimolar).
C CONSTS(60) is K_Thirteen_BPG in component V_GAPDH (millimolar).
C CONSTS(61) is Keq_GAPDH in component V_GAPDH (dimensionless).
C CONSTS(109) is V_max_r in component V_GAPDH (flux).
C CONSTS(62) is V_max_f in component V_GAPDH (flux).
C CONSTS(63) is K_three_PG in component V_PGK (millimolar).
C CONSTS(64) is K_ADP in component V_PGK (millimolar).
C CONSTS(65) is K_ATP in component V_PGK (millimolar).
C CONSTS(66) is K_Thirteen_BPG in component V_PGK (millimolar).
C CONSTS(67) is Keq_PGK in component V_PGK (dimensionless).
C CONSTS(68) is V_max_r in component V_PGK (flux).
C CONSTS(110) is V_max_f in component V_PGK (flux).
C CONSTS(69) is K_two_PG in component V_PGM (millimolar).
C CONSTS(70) is K_three_PG in component V_PGM (millimolar).
C CONSTS(71) is Keq_PGM in component V_PGM (dimensionless).
C CONSTS(111) is V_max_r in component V_PGM (flux).
C CONSTS(72) is V_max_f in component V_PGM (flux).
C CONSTS(73) is K_PEP in component V_ENOL (millimolar).
C CONSTS(74) is K_two_PG in component V_ENOL (millimolar).
C CONSTS(75) is Keq_ENOL in component V_ENOL (dimensionless).
C CONSTS(112) is V_max_r in component V_ENOL (flux).
C CONSTS(76) is V_max_f in component V_ENOL (flux).
C CONSTS(77) is K_PYR in component V_PK (millimolar).
C CONSTS(78) is K_ADP in component V_PK (millimolar).
C CONSTS(79) is K_ATP in component V_PK (millimolar).
C CONSTS(80) is K_PEP in component V_PK (millimolar).
C CONSTS(81) is Keq_PK in component V_PK (dimensionless).
C CONSTS(113) is V_max_r in component V_PK (flux).
C CONSTS(82) is V_max_f in component V_PK (flux).
C CONSTS(83) is K_LAC in component V_LDH (millimolar).
C CONSTS(84) is K_NADH in component V_LDH (millimolar).
C CONSTS(85) is K_NAD in component V_LDH (millimolar).
C CONSTS(86) is K_PYR in component V_LDH (millimolar).
C CONSTS(87) is Keq_LDH in component V_LDH (dimensionless).
C CONSTS(114) is V_max_r in component V_LDH (flux).
C CONSTS(88) is V_max_f in component V_LDH (flux).
C CONSTS(89) is Ki_ADP in component V_CK (millimolar).
C CONSTS(90) is K_Cr in component V_CK (millimolar).
C CONSTS(91) is K_PCr in component V_CK (millimolar).
C CONSTS(92) is Ki_PCr in component V_CK (millimolar).
C CONSTS(93) is Ki_ATP in component V_CK (millimolar).
C CONSTS(94) is Keq_CK in component V_CK (dimensionless).
C CONSTS(95) is V_max_r in component V_CK (flux).
C CONSTS(115) is V_max_f in component V_CK (flux).
C CONSTS(96) is K_ADP in component V_ADK (millimolar).
C CONSTS(97) is K_AMP in component V_ADK (millimolar).
C CONSTS(98) is K_ATP in component V_ADK (millimolar).
C CONSTS(99) is Keq_ADK in component V_ADK (dimensionless).
C CONSTS(116) is V_max_r in component V_ADK (flux).
C CONSTS(100) is V_max_f in component V_ADK (flux).
C RATES(1) is d/dt GLY in component GLY (millimolar).
C RATES(2) is d/dt G1P in component G1P (millimolar).
C RATES(3) is d/dt G6P in component G6P (millimolar).
C RATES(4) is d/dt F6P in component F6P (millimolar).
C RATES(5) is d/dt FBP in component FBP (millimolar).
C RATES(6) is d/dt DHAP in component DHAP (millimolar).
C RATES(7) is d/dt GAP in component GAP (millimolar).
C RATES(8) is d/dt Thirteen_BPG in component Thirteen_BPG (millimolar).
C RATES(9) is d/dt three_PG in component three_PG (millimolar).
C RATES(10) is d/dt two_PG in component two_PG (millimolar).
C RATES(11) is d/dt PEP in component PEP (millimolar).
C RATES(12) is d/dt PYR in component PYR (millimolar).
C RATES(13) is d/dt LAC in component LAC (millimolar).
C RATES(14) is d/dt Pi in component Pi (millimolar).
C RATES(15) is d/dt ADP in component ADP (millimolar).
C RATES(16) is d/dt ATP in component ATP (millimolar).
C RATES(17) is d/dt AMP in component AMP (millimolar).
C RATES(18) is d/dt PCr in component PCr (millimolar).
C RATES(19) is d/dt Cr in component Cr (millimolar).
C RATES(20) is d/dt NADH in component NADH (millimolar).
C RATES(21) is d/dt NAD in component NAD (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 112.0
      STATES(2) = 0.0589
      STATES(3) = 0.75
      STATES(4) = 0.228
      STATES(5) = 0.0723
      STATES(6) = 0.0764
      STATES(7) = 0.0355
      STATES(8) = 0.065
      STATES(9) = 0.052
      STATES(10) = 0.005
      STATES(11) = 0.0194
      STATES(12) = 0.0994
      STATES(13) = 1.3
      STATES(14) = 4.1
      CONSTS(1) = 600.0
      STATES(15) = 0.013
      STATES(16) = 8.2
      STATES(17) = 0.99
      STATES(18) = 34.67
      STATES(19) = 35.0
      STATES(20) = 0.001
      STATES(21) = 0.5
      CONSTS(2) = 0.4
      CONSTS(3) = 0.6
      CONSTS(4) = 2.0
      CONSTS(5) = 2.0
      CONSTS(6) = 1.7
      CONSTS(7) = 0.15
      CONSTS(8) = 0.42
      CONSTS(9) = 4.7
      CONSTS(10) = 4.0
      CONSTS(11) = 10.1
      CONSTS(12) = 20.0
      CONSTS(13) = 15.0
      CONSTS(14) = 4.4
      CONSTS(15) = 1.7
      CONSTS(16) = 0.42
      CONSTS(17) = 1.75
      CONSTS(18) = 1.9E-6
      CONSTS(19) = 0.15
      CONSTS(20) = 4.6
      CONSTS(21) = 0.2
      CONSTS(22) = 7.4
      CONSTS(23) = 1.5
      CONSTS(24) = 30.0
      CONSTS(25) = 0.063
      CONSTS(26) = 0.03
      CONSTS(27) = 16.62
      CONSTS(28) = 480.0
      CONSTS(29) = 0.18
      CONSTS(30) = 0.119
      CONSTS(31) = 0.45
      CONSTS(32) = 880.0
      CONSTS(33) = 4.02
      CONSTS(34) = 4.02
      CONSTS(35) = 0.18
      CONSTS(36) = 0.08
      CONSTS(37) = 20.0
      CONSTS(38) = 0.25
      CONSTS(39) = 2.7
      CONSTS(40) = 2.7
      CONSTS(41) = 0.87
      CONSTS(42) = 0.06
      CONSTS(43) = 56.0
      CONSTS(44) = 13.0
      CONSTS(45) = 0.01
      CONSTS(46) = 0.01
      CONSTS(47) = 1.0
      CONSTS(48) = 0.05
      CONSTS(49) = 2.0
      CONSTS(50) = 9.5E-5
      CONSTS(51) = 104.0
      CONSTS(52) = 0.32
      CONSTS(53) = 0.61
      CONSTS(54) = 0.052
      CONSTS(55) = 12000.0
      CONSTS(56) = 0.0025
      CONSTS(57) = 0.29
      CONSTS(58) = 0.09
      CONSTS(59) = 0.0033
      CONSTS(60) = 0.0008
      CONSTS(61) = 0.089
      CONSTS(62) = 1265.0
      CONSTS(63) = 1.2
      CONSTS(64) = 0.008
      CONSTS(65) = 0.35
      CONSTS(66) = 0.002
      CONSTS(67) = 57109.0
      CONSTS(68) = 1120.0
      CONSTS(69) = 0.014
      CONSTS(70) = 0.2
      CONSTS(71) = 0.18
      CONSTS(72) = 1120.0
      CONSTS(73) = 0.37
      CONSTS(74) = 0.1
      CONSTS(75) = 0.49
      CONSTS(76) = 192.0
      CONSTS(77) = 7.05
      CONSTS(78) = 0.3
      CONSTS(79) = 1.13
      CONSTS(80) = 0.08
      CONSTS(81) = 10304.0
      CONSTS(82) = 1440.0
      CONSTS(83) = 17.0
      CONSTS(84) = 0.002
      CONSTS(85) = 0.849
      CONSTS(86) = 0.335
      CONSTS(87) = 16198.0
      CONSTS(88) = 1920.0
      CONSTS(89) = 0.135
      CONSTS(90) = 3.8
      CONSTS(91) = 1.11
      CONSTS(92) = 3.9
      CONSTS(93) = 3.5
      CONSTS(94) = 233.0
      CONSTS(95) = 500.0
      CONSTS(96) = 0.35
      CONSTS(97) = 0.32
      CONSTS(98) = 0.27
      CONSTS(99) = 2.21
      CONSTS(100) = 880.0
      CONSTS(101) = ( CONSTS(12)*CONSTS(7)*CONSTS(11))/( CONSTS(4)*CONSTS(10)*CONSTS(8))
      CONSTS(102) = ( CONSTS(24)*CONSTS(14)*CONSTS(23))/( CONSTS(15)*CONSTS(20)*CONSTS(16))
      CONSTS(103) = ( CONSTS(28)*CONSTS(26))/( CONSTS(25)*CONSTS(27))
      CONSTS(104) = ( CONSTS(32)*CONSTS(29)*CONSTS(31))/CONSTS(30)
      CONSTS(105) = ( CONSTS(43)*CONSTS(39)*CONSTS(33))/( CONSTS(36)*CONSTS(35))
      CONSTS(106) = ( CONSTS(35)*CONSTS(36))/( CONSTS(37)*CONSTS(38))
      CONSTS(107) = ( CONSTS(51)*CONSTS(49)*CONSTS(47))/( CONSTS(48)*CONSTS(50))
      CONSTS(108) = ( CONSTS(55)*CONSTS(53))/( CONSTS(52)*CONSTS(54))
      CONSTS(109) = ( CONSTS(62)*1.00000*CONSTS(60)*CONSTS(59))/( CONSTS(56)*CONSTS(58)*CONSTS(57)*CONSTS(61))
      CONSTS(110) = ( CONSTS(68)*CONSTS(66)*CONSTS(64)*CONSTS(67))/( CONSTS(63)*CONSTS(65))
      CONSTS(111) = ( CONSTS(72)*CONSTS(69))/( CONSTS(70)*CONSTS(71))
      CONSTS(112) = ( CONSTS(76)*CONSTS(73))/( CONSTS(74)*CONSTS(75))
      CONSTS(113) = ( CONSTS(82)*CONSTS(79)*CONSTS(77))/( CONSTS(80)*CONSTS(78)*CONSTS(81))
      CONSTS(114) = ( CONSTS(88)*CONSTS(83)*CONSTS(85))/( CONSTS(86)*CONSTS(84)*CONSTS(87))
      CONSTS(115) = ( CONSTS(95)*CONSTS(93)*CONSTS(90)*CONSTS(94))/( CONSTS(89)*CONSTS(91))
      CONSTS(116) = ( CONSTS(100)*CONSTS(96) ** 2.00000)/( CONSTS(98)*CONSTS(97)*CONSTS(99))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = ( CONSTS(12)*(( STATES(1)*STATES(14))/( CONSTS(4)*CONSTS(10))) -  CONSTS(101)*(( STATES(1)*STATES(2))/( CONSTS(5)*CONSTS(11))))/(1.00000+STATES(1)/CONSTS(4)+STATES(14)/CONSTS(9)+STATES(1)/CONSTS(5)+STATES(2)/CONSTS(11)+( STATES(1)*STATES(14))/( CONSTS(6)*CONSTS(9))+( STATES(1)*STATES(2))/( CONSTS(7)*CONSTS(11)))
      ALGBRC(3) =  (( CONSTS(24)*(( STATES(1)*STATES(14))/( CONSTS(13)*CONSTS(21))) -  CONSTS(102)*(( STATES(1)*STATES(2))/( CONSTS(14)*CONSTS(22))))/(1.00000+STATES(1)/CONSTS(13)+STATES(14)/CONSTS(20)+STATES(1)/CONSTS(14)+STATES(2)/CONSTS(22)+( STATES(1)*STATES(14))/( CONSTS(15)*CONSTS(20))+( STATES(1)*STATES(2))/( CONSTS(19)*CONSTS(22))))*((STATES(17) ** CONSTS(17)/CONSTS(18))/(1.00000+STATES(17) ** CONSTS(17)/CONSTS(18)))
      ALGBRC(4) =  CONSTS(2)*ALGBRC(2)+ CONSTS(3)*ALGBRC(3)
      RATES(1) = - ALGBRC(4)
      ALGBRC(5) = ( CONSTS(28)*(STATES(2)/CONSTS(25)) -  CONSTS(103)*(STATES(3)/CONSTS(26)))/(1.00000+STATES(2)/CONSTS(25)+STATES(3)/CONSTS(26))
      RATES(2) = ALGBRC(4) - ALGBRC(5)
      ALGBRC(6) = ( CONSTS(104)*(STATES(3)/CONSTS(29)) -  CONSTS(32)*(STATES(4)/CONSTS(30)))/(1.00000+STATES(3)/CONSTS(29)+STATES(4)/CONSTS(30))
      RATES(3) = ALGBRC(5) - ALGBRC(6)
      ALGBRC(9) =  CONSTS(44)* ((1.00000+STATES(16)/CONSTS(41))/(1.00000+ CONSTS(45)*(STATES(16)/CONSTS(41))))*((1.00000+ CONSTS(46)*(STATES(17)/CONSTS(42)))/(1.00000+STATES(17)/CONSTS(42))) ** 4.00000
      ALGBRC(7) =  (1.00000+STATES(4)/CONSTS(35))*(1.00000+STATES(16)/CONSTS(36))+STATES(15)/CONSTS(39)+ (STATES(5)/CONSTS(33))*(1.00000+STATES(15)/CONSTS(39))
      ALGBRC(8) =  (1.00000+STATES(4)/CONSTS(37))*(1.00000+STATES(16)/CONSTS(38))+STATES(15)/CONSTS(40)+ (STATES(5)/CONSTS(34))*(1.00000+STATES(15)/CONSTS(40))
      ALGBRC(10) =  (( CONSTS(43)*(( STATES(16)*STATES(4))/( CONSTS(36)*CONSTS(35))) -  CONSTS(105)*(( STATES(15)*STATES(5))/( CONSTS(39)*CONSTS(33))))/ALGBRC(7))*((1.00000+ CONSTS(106)*ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 3.00000)/(1.00000+ ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 4.00000))
      RATES(4) = ALGBRC(6) - ALGBRC(10)
      ALGBRC(11) = ( CONSTS(51)*(STATES(5)/CONSTS(48)) -  CONSTS(107)*(( STATES(6)*STATES(7))/( CONSTS(49)*CONSTS(47))))/(1.00000+STATES(5)/CONSTS(48)+STATES(6)/CONSTS(49)+STATES(7)/CONSTS(47))
      RATES(5) = ALGBRC(10) - ALGBRC(11)
      ALGBRC(12) = ( CONSTS(55)*(STATES(7)/CONSTS(52)) -  CONSTS(108)*(STATES(6)/CONSTS(53)))/(1.00000+STATES(6)/CONSTS(53)+STATES(7)/CONSTS(52))
      RATES(6) = ALGBRC(11)+ALGBRC(12)
      ALGBRC(13) = 1.00000+STATES(7)/CONSTS(56)+STATES(21)/CONSTS(58)+STATES(14)/CONSTS(57)+( STATES(7)*STATES(21))/( CONSTS(56)*CONSTS(58))+( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))+STATES(8)/CONSTS(60)+STATES(20)/CONSTS(59)+( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))
      ALGBRC(14) = ( CONSTS(62)*(( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))) -  CONSTS(109)*(( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))))/ALGBRC(13)
      RATES(7) = ALGBRC(11) - (ALGBRC(12)+ALGBRC(14))
      RATES(14) = CONSTS(1) - (ALGBRC(4)+ALGBRC(14))
      ALGBRC(15) = ( CONSTS(110)*(( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))) -  CONSTS(68)*(( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65))))/(1.00000+STATES(8)/CONSTS(66)+STATES(15)/CONSTS(64)+( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))+STATES(9)/CONSTS(63)+STATES(16)/CONSTS(65)+( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65)))
      RATES(8) = ALGBRC(14) - ALGBRC(15)
      ALGBRC(16) = ( CONSTS(72)*(STATES(9)/CONSTS(70)) -  CONSTS(111)*(STATES(10)/CONSTS(69)))/(1.00000+STATES(9)/CONSTS(70)+STATES(10)/CONSTS(69))
      RATES(9) = ALGBRC(15) - ALGBRC(16)
      ALGBRC(17) = ( CONSTS(76)*(STATES(10)/CONSTS(74)) -  CONSTS(112)*(STATES(11)/CONSTS(73)))/(1.00000+STATES(10)/CONSTS(74)+STATES(11)/CONSTS(73))
      RATES(10) = ALGBRC(16) - ALGBRC(17)
      ALGBRC(18) = ( CONSTS(82)*(( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))) -  CONSTS(113)*(( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79))))/(1.00000+STATES(11)/CONSTS(80)+STATES(15)/CONSTS(78)+( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))+STATES(12)/CONSTS(77)+STATES(16)/CONSTS(79)+( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79)))
      RATES(11) = ALGBRC(17) - ALGBRC(18)
      ALGBRC(19) = ( CONSTS(88)*(( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))) -  CONSTS(114)*(( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85))))/(1.00000+STATES(12)/CONSTS(86)+STATES(20)/CONSTS(84)+( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))+STATES(13)/CONSTS(83)+STATES(21)/CONSTS(85)+( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85)))
      RATES(12) = ALGBRC(18) - ALGBRC(19)
      ALGBRC(1) =  0.200000*STATES(13)
      RATES(13) = ALGBRC(19) - ALGBRC(1)
      ALGBRC(20) = ( CONSTS(115)*(( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90))) -  CONSTS(95)*(( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))))/(1.00000+STATES(15)/CONSTS(89)+STATES(18)/CONSTS(92)+( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))+STATES(16)/CONSTS(93)+( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90)))
      RATES(18) = ALGBRC(20)
      RATES(19) = - ALGBRC(20)
      RATES(20) = ALGBRC(14) - ALGBRC(19)
      RATES(21) = ALGBRC(19) - ALGBRC(14)
      ALGBRC(21) = ( CONSTS(100)*(( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))) -  CONSTS(116)*(STATES(15) ** 2.00000/CONSTS(96) ** 2.00000))/(1.00000+STATES(16)/CONSTS(98)+STATES(17)/CONSTS(97)+( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))+( 2.00000*STATES(15))/CONSTS(96)+STATES(15) ** 2.00000/CONSTS(96) ** 2.00000)
      RATES(15) = ALGBRC(10)+ 2.00000*ALGBRC(21)+ALGBRC(20)+CONSTS(1)+- (ALGBRC(15)+ALGBRC(18))
      RATES(16) = ALGBRC(15)+ALGBRC(18)+- (ALGBRC(10)+ALGBRC(21)+ALGBRC(20)+CONSTS(1))
      RATES(17) = - ALGBRC(21)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = ( CONSTS(12)*(( STATES(1)*STATES(14))/( CONSTS(4)*CONSTS(10))) -  CONSTS(101)*(( STATES(1)*STATES(2))/( CONSTS(5)*CONSTS(11))))/(1.00000+STATES(1)/CONSTS(4)+STATES(14)/CONSTS(9)+STATES(1)/CONSTS(5)+STATES(2)/CONSTS(11)+( STATES(1)*STATES(14))/( CONSTS(6)*CONSTS(9))+( STATES(1)*STATES(2))/( CONSTS(7)*CONSTS(11)))
      ALGBRC(3) =  (( CONSTS(24)*(( STATES(1)*STATES(14))/( CONSTS(13)*CONSTS(21))) -  CONSTS(102)*(( STATES(1)*STATES(2))/( CONSTS(14)*CONSTS(22))))/(1.00000+STATES(1)/CONSTS(13)+STATES(14)/CONSTS(20)+STATES(1)/CONSTS(14)+STATES(2)/CONSTS(22)+( STATES(1)*STATES(14))/( CONSTS(15)*CONSTS(20))+( STATES(1)*STATES(2))/( CONSTS(19)*CONSTS(22))))*((STATES(17) ** CONSTS(17)/CONSTS(18))/(1.00000+STATES(17) ** CONSTS(17)/CONSTS(18)))
      ALGBRC(4) =  CONSTS(2)*ALGBRC(2)+ CONSTS(3)*ALGBRC(3)
      ALGBRC(5) = ( CONSTS(28)*(STATES(2)/CONSTS(25)) -  CONSTS(103)*(STATES(3)/CONSTS(26)))/(1.00000+STATES(2)/CONSTS(25)+STATES(3)/CONSTS(26))
      ALGBRC(6) = ( CONSTS(104)*(STATES(3)/CONSTS(29)) -  CONSTS(32)*(STATES(4)/CONSTS(30)))/(1.00000+STATES(3)/CONSTS(29)+STATES(4)/CONSTS(30))
      ALGBRC(9) =  CONSTS(44)* ((1.00000+STATES(16)/CONSTS(41))/(1.00000+ CONSTS(45)*(STATES(16)/CONSTS(41))))*((1.00000+ CONSTS(46)*(STATES(17)/CONSTS(42)))/(1.00000+STATES(17)/CONSTS(42))) ** 4.00000
      ALGBRC(7) =  (1.00000+STATES(4)/CONSTS(35))*(1.00000+STATES(16)/CONSTS(36))+STATES(15)/CONSTS(39)+ (STATES(5)/CONSTS(33))*(1.00000+STATES(15)/CONSTS(39))
      ALGBRC(8) =  (1.00000+STATES(4)/CONSTS(37))*(1.00000+STATES(16)/CONSTS(38))+STATES(15)/CONSTS(40)+ (STATES(5)/CONSTS(34))*(1.00000+STATES(15)/CONSTS(40))
      ALGBRC(10) =  (( CONSTS(43)*(( STATES(16)*STATES(4))/( CONSTS(36)*CONSTS(35))) -  CONSTS(105)*(( STATES(15)*STATES(5))/( CONSTS(39)*CONSTS(33))))/ALGBRC(7))*((1.00000+ CONSTS(106)*ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 3.00000)/(1.00000+ ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 4.00000))
      ALGBRC(11) = ( CONSTS(51)*(STATES(5)/CONSTS(48)) -  CONSTS(107)*(( STATES(6)*STATES(7))/( CONSTS(49)*CONSTS(47))))/(1.00000+STATES(5)/CONSTS(48)+STATES(6)/CONSTS(49)+STATES(7)/CONSTS(47))
      ALGBRC(12) = ( CONSTS(55)*(STATES(7)/CONSTS(52)) -  CONSTS(108)*(STATES(6)/CONSTS(53)))/(1.00000+STATES(6)/CONSTS(53)+STATES(7)/CONSTS(52))
      ALGBRC(13) = 1.00000+STATES(7)/CONSTS(56)+STATES(21)/CONSTS(58)+STATES(14)/CONSTS(57)+( STATES(7)*STATES(21))/( CONSTS(56)*CONSTS(58))+( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))+STATES(8)/CONSTS(60)+STATES(20)/CONSTS(59)+( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))
      ALGBRC(14) = ( CONSTS(62)*(( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))) -  CONSTS(109)*(( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))))/ALGBRC(13)
      ALGBRC(15) = ( CONSTS(110)*(( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))) -  CONSTS(68)*(( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65))))/(1.00000+STATES(8)/CONSTS(66)+STATES(15)/CONSTS(64)+( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))+STATES(9)/CONSTS(63)+STATES(16)/CONSTS(65)+( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65)))
      ALGBRC(16) = ( CONSTS(72)*(STATES(9)/CONSTS(70)) -  CONSTS(111)*(STATES(10)/CONSTS(69)))/(1.00000+STATES(9)/CONSTS(70)+STATES(10)/CONSTS(69))
      ALGBRC(17) = ( CONSTS(76)*(STATES(10)/CONSTS(74)) -  CONSTS(112)*(STATES(11)/CONSTS(73)))/(1.00000+STATES(10)/CONSTS(74)+STATES(11)/CONSTS(73))
      ALGBRC(18) = ( CONSTS(82)*(( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))) -  CONSTS(113)*(( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79))))/(1.00000+STATES(11)/CONSTS(80)+STATES(15)/CONSTS(78)+( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))+STATES(12)/CONSTS(77)+STATES(16)/CONSTS(79)+( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79)))
      ALGBRC(19) = ( CONSTS(88)*(( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))) -  CONSTS(114)*(( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85))))/(1.00000+STATES(12)/CONSTS(86)+STATES(20)/CONSTS(84)+( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))+STATES(13)/CONSTS(83)+STATES(21)/CONSTS(85)+( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85)))
      ALGBRC(1) =  0.200000*STATES(13)
      ALGBRC(20) = ( CONSTS(115)*(( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90))) -  CONSTS(95)*(( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))))/(1.00000+STATES(15)/CONSTS(89)+STATES(18)/CONSTS(92)+( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))+STATES(16)/CONSTS(93)+( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90)))
      ALGBRC(21) = ( CONSTS(100)*(( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))) -  CONSTS(116)*(STATES(15) ** 2.00000/CONSTS(96) ** 2.00000))/(1.00000+STATES(16)/CONSTS(98)+STATES(17)/CONSTS(97)+( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))+( 2.00000*STATES(15))/CONSTS(96)+STATES(15) ** 2.00000/CONSTS(96) ** 2.00000)
      RETURN
      END