C
C There are a total of 24 entries in the algebraic variable array.
C There are a total of 24 entries in each of the rate and state variable arrays.
C There are a total of 58 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is Glc_x_0 in component Glc_x_0 (millimolar).
C ALGBRC(1) is inGlc in component inGlc (flux).
C STATES(2) is Glc_x in component Glc_x (millimolar).
C ALGBRC(2) is GlcTrans in component GlcTrans (flux).
C STATES(3) is Glc in component Glc (millimolar).
C ALGBRC(3) is HK in component HK (flux).
C STATES(4) is G6P in component G6P (millimolar).
C ALGBRC(4) is PGI in component PGI (flux).
C ALGBRC(11) is storage in component storage (flux).
C STATES(5) is F6P in component F6P (millimolar).
C ALGBRC(5) is PFK in component PFK (flux).
C STATES(6) is FBP in component FBP (millimolar).
C ALGBRC(6) is ALD in component ALD (flux).
C STATES(7) is GAP in component GAP (millimolar).
C ALGBRC(8) is GAPDH in component GAPDH (flux).
C ALGBRC(7) is TIM in component TIM (flux).
C STATES(8) is DHAP in component DHAP (millimolar).
C ALGBRC(15) is lpGlyc in component lpGlyc (flux).
C STATES(9) is BPG in component BPG (millimolar).
C ALGBRC(9) is lpPEP in component lpPEP (flux).
C STATES(10) is PEP in component PEP (millimolar).
C ALGBRC(10) is PK in component PK (flux).
C STATES(11) is Pyr in component Pyr (millimolar).
C ALGBRC(12) is PDC in component PDC (flux).
C STATES(12) is ATP in component ATP (millimolar).
C ALGBRC(14) is consum in component consum (flux).
C ALGBRC(16) is AK in component AK (flux).
C STATES(13) is ADP in component ADP (millimolar).
C STATES(14) is AMP in component AMP (millimolar).
C STATES(15) is CN_x_0 in component CN_x_0 (millimolar).
C ALGBRC(24) is inCN in component inCN (flux).
C STATES(16) is CN_x in component CN_x (millimolar).
C ALGBRC(22) is lacto in component lacto (flux).
C STATES(17) is ACA in component ACA (millimolar).
C ALGBRC(19) is difACA in component difACA (flux).
C ALGBRC(13) is ADH in component ADH (flux).
C STATES(18) is ACA_x in component ACA_x (millimolar).
C ALGBRC(17) is outACA in component outACA (flux).
C STATES(19) is EtOH in component EtOH (millimolar).
C ALGBRC(20) is difEtOH in component difEtOH (flux).
C STATES(20) is EtOH_x in component EtOH_x (millimolar).
C ALGBRC(18) is outEtOH in component outEtOH (flux).
C STATES(21) is Glyc in component Glyc (millimolar).
C ALGBRC(23) is difGlyc in component difGlyc (flux).
C STATES(22) is Glyc_x in component Glyc_x (millimolar).
C ALGBRC(21) is outGlyc in component outGlyc (flux).
C STATES(23) is NADH in component NADH (millimolar).
C STATES(24) is NAD in component NAD (millimolar).
C CONSTS(1) is k0 in component model_parameters (first_order_rate_constant).
C CONSTS(2) is y_vol in component model_parameters (dimensionless).
C CONSTS(3) is K2_Glc in component GlcTrans (millimolar).
C CONSTS(4) is K2_IG6P in component GlcTrans (millimolar).
C CONSTS(5) is K2_IIG6P in component GlcTrans (millimolar).
C CONSTS(6) is V_2m in component GlcTrans (flux).
C CONSTS(7) is P2 in component GlcTrans (dimensionless).
C CONSTS(8) is K3_DGlc in component HK (millimolar).
C CONSTS(9) is K3_Glc in component HK (millimolar).
C CONSTS(10) is K3_ATP in component HK (millimolar).
C CONSTS(11) is V_3m in component HK (flux).
C CONSTS(12) is V_4m in component PGI (flux).
C CONSTS(13) is K4_G6P in component PGI (millimolar).
C CONSTS(14) is K4_F6P in component PGI (millimolar).
C CONSTS(15) is K4_eq in component PGI (dimensionless).
C CONSTS(16) is K5 in component PFK (millimolar2).
C CONSTS(17) is kappa5 in component PFK (dimensionless).
C CONSTS(18) is V_5m in component PFK (flux).
C CONSTS(19) is K6_eq in component ALD (millimolar).
C CONSTS(20) is K6_FBP in component ALD (millimolar).
C CONSTS(21) is K6_DHAP in component ALD (millimolar).
C CONSTS(22) is K6_GAP in component ALD (millimolar).
C CONSTS(23) is K6_IGAP in component ALD (millimolar).
C CONSTS(24) is V_6r in component ALD (flux).
C CONSTS(25) is V_6f in component ALD (flux).
C CONSTS(26) is K7_eq in component TIM (dimensionless).
C CONSTS(27) is K7_DHAP in component TIM (millimolar).
C CONSTS(28) is K7_GAP in component TIM (millimolar).
C CONSTS(29) is V_7m in component TIM (flux).
C CONSTS(30) is K8_NAD in component GAPDH (millimolar).
C CONSTS(31) is K8_NADH in component GAPDH (millimolar).
C CONSTS(32) is K8_GAP in component GAPDH (millimolar).
C CONSTS(33) is K8_BPG in component GAPDH (millimolar).
C CONSTS(34) is K8_eq in component GAPDH (dimensionless).
C CONSTS(35) is V_8m in component GAPDH (flux).
C CONSTS(36) is k9f in component lpPEP (second_order_rate_constant).
C CONSTS(37) is k9r in component lpPEP (second_order_rate_constant).
C CONSTS(38) is K10_PEP in component PK (millimolar).
C CONSTS(39) is K10_ADP in component PK (millimolar).
C CONSTS(40) is V_10m in component PK (flux).
C CONSTS(41) is K11 in component PDC (millimolar).
C CONSTS(42) is V_11m in component PDC (flux).
C CONSTS(43) is K12_NADH in component ADH (millimolar).
C CONSTS(44) is K12_ACA in component ADH (millimolar).
C CONSTS(45) is V_12m in component ADH (flux).
C CONSTS(46) is k13 in component difEtOH (first_order_rate_constant).
C CONSTS(47) is K15_NADH in component lpGlyc (millimolar).
C CONSTS(48) is K15_INADH in component lpGlyc (millimolar).
C CONSTS(49) is K15_INAD in component lpGlyc (millimolar).
C CONSTS(50) is K15_DHAP in component lpGlyc (millimolar).
C CONSTS(51) is V_15m in component lpGlyc (flux).
C CONSTS(52) is k16 in component difGlyc (first_order_rate_constant).
C CONSTS(53) is k18 in component difACA (first_order_rate_constant).
C CONSTS(54) is k20 in component lacto (second_order_rate_constant).
C CONSTS(55) is k22 in component storage (second_order_rate_constant).
C CONSTS(56) is k23 in component consum (first_order_rate_constant).
C CONSTS(57) is k24f in component AK (second_order_rate_constant).
C CONSTS(58) is k24r in component AK (second_order_rate_constant).
C RATES(1) is d/dt Glc_x_0 in component Glc_x_0 (millimolar).
C RATES(2) is d/dt Glc_x in component Glc_x (millimolar).
C RATES(3) is d/dt Glc in component Glc (millimolar).
C RATES(4) is d/dt G6P in component G6P (millimolar).
C RATES(5) is d/dt F6P in component F6P (millimolar).
C RATES(6) is d/dt FBP in component FBP (millimolar).
C RATES(7) is d/dt GAP in component GAP (millimolar).
C RATES(8) is d/dt DHAP in component DHAP (millimolar).
C RATES(9) is d/dt BPG in component BPG (millimolar).
C RATES(10) is d/dt PEP in component PEP (millimolar).
C RATES(11) is d/dt Pyr in component Pyr (millimolar).
C RATES(12) is d/dt ATP in component ATP (millimolar).
C RATES(13) is d/dt ADP in component ADP (millimolar).
C RATES(14) is d/dt AMP in component AMP (millimolar).
C RATES(15) is d/dt CN_x_0 in component CN_x_0 (millimolar).
C RATES(16) is d/dt CN_x in component CN_x (millimolar).
C RATES(17) is d/dt ACA in component ACA (millimolar).
C RATES(18) is d/dt ACA_x in component ACA_x (millimolar).
C RATES(19) is d/dt EtOH in component EtOH (millimolar).
C RATES(20) is d/dt EtOH_x in component EtOH_x (millimolar).
C RATES(21) is d/dt Glyc in component Glyc (millimolar).
C RATES(22) is d/dt Glyc_x in component Glyc_x (millimolar).
C RATES(23) is d/dt NADH in component NADH (millimolar).
C RATES(24) is d/dt NAD in component NAD (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 24.0
      STATES(2) = 6.7
      STATES(3) = 0.573074
      STATES(4) = 4.2
      STATES(5) = 0.49
      STATES(6) = 4.64
      STATES(7) = 0.115
      STATES(8) = 2.95
      STATES(9) = 0.00027
      STATES(10) = 0.04
      STATES(11) = 8.7
      STATES(12) = 2.1
      STATES(13) = 1.5
      STATES(14) = 0.33
      STATES(15) = 5.60
      STATES(16) = 5.20358
      STATES(17) = 1.48153
      STATES(18) = 1.28836
      STATES(19) = 19.2379
      STATES(20) = 16.4514
      STATES(21) = 4.196
      STATES(22) = 1.68478
      STATES(23) = 0.33
      STATES(24) = 0.65
      CONSTS(1) = 0.048
      CONSTS(2) = 59.0
      CONSTS(3) = 1.7
      CONSTS(4) = 1.2
      CONSTS(5) = 7.2
      CONSTS(6) = 1014.96
      CONSTS(7) = 1.0
      CONSTS(8) = 0.37
      CONSTS(9) = 0.0
      CONSTS(10) = 0.1
      CONSTS(11) = 5.7547
      CONSTS(12) = 496.042
      CONSTS(13) = 0.8
      CONSTS(14) = 0.15
      CONSTS(15) = 0.13
      CONSTS(16) = 0.021
      CONSTS(17) = 0.15
      CONSTS(18) = 45.4327
      CONSTS(19) = 0.081
      CONSTS(20) = 0.3
      CONSTS(21) = 2.0
      CONSTS(22) = 4.0
      CONSTS(23) = 10.0
      CONSTS(24) = 1.10391E4
      CONSTS(25) = 2.20782E3
      CONSTS(26) = 0.055
      CONSTS(27) = 1.23
      CONSTS(28) = 1.27
      CONSTS(29) = 1.16365E2
      CONSTS(30) = 0.1
      CONSTS(31) = 0.06
      CONSTS(32) = 0.6
      CONSTS(33) = 0.01
      CONSTS(34) = 0.0055
      CONSTS(35) = 8.33858E2
      CONSTS(36) = 4.43866E5
      CONSTS(37) = 1.52862E3
      CONSTS(38) = 0.2
      CONSTS(39) = 0.17
      CONSTS(40) = 3.43096E2
      CONSTS(41) = 0.3
      CONSTS(42) = 5.31328E1
      CONSTS(43) = 0.1
      CONSTS(44) = 0.71
      CONSTS(45) = 8.98023E1
      CONSTS(46) = 16.72
      CONSTS(47) = 0.13
      CONSTS(48) = 0.034
      CONSTS(49) = 0.13
      CONSTS(50) = 25.0
      CONSTS(51) = 8.14797E1
      CONSTS(52) = 1.9
      CONSTS(53) = 24.7
      CONSTS(54) = 2.83828E-3
      CONSTS(55) = 2.25932
      CONSTS(56) = 3.20760
      CONSTS(57) = 4.32900E2
      CONSTS(58) = 1.33333E2
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*(STATES(1) - STATES(2))
      RATES(1) = ALGBRC(1)
      ALGBRC(2) =  (CONSTS(6)/CONSTS(2))*((STATES(2)/CONSTS(3))/(1.00000+STATES(2)/CONSTS(3)+ (( CONSTS(7)*(STATES(2)/CONSTS(3))+1.00000)/( CONSTS(7)*(STATES(3)/CONSTS(3))+1.00000))*(1.00000+STATES(3)/CONSTS(3))+STATES(4)/CONSTS(4)+( STATES(3)*STATES(4))/( CONSTS(3)*CONSTS(5)))) -  (CONSTS(6)/CONSTS(2))*((STATES(3)/CONSTS(3))/(1.00000+STATES(3)/CONSTS(3)+ (( CONSTS(7)*(STATES(3)/CONSTS(3))+1.00000)/( CONSTS(7)*(STATES(2)/CONSTS(3))+1.00000))*(1.00000+STATES(2)/CONSTS(3))+STATES(4)/CONSTS(4)+( STATES(3)*STATES(4))/( CONSTS(3)*CONSTS(5))))
      RATES(2) = ALGBRC(1)+ALGBRC(2)
      ALGBRC(3) = ( CONSTS(11)*STATES(12)*STATES(3))/( CONSTS(8)*CONSTS(10)+ CONSTS(9)*STATES(12)+ CONSTS(10)*STATES(3)+ STATES(3)*STATES(12))
      RATES(3) = ALGBRC(2)+ALGBRC(3)
      ALGBRC(4) = ( CONSTS(12)*STATES(4))/(CONSTS(13)+STATES(4)+ (CONSTS(13)/CONSTS(14))*STATES(5)) - ( CONSTS(12)*(STATES(5)/CONSTS(15)))/(CONSTS(13)+STATES(4)+ (CONSTS(13)/CONSTS(14))*STATES(5))
      ALGBRC(5) = ( CONSTS(18)*STATES(5) ** 2.00000)/( CONSTS(16)*(1.00000+ CONSTS(17)*(STATES(12)/STATES(14)))+STATES(5) ** 2.00000)
      RATES(5) = ALGBRC(4)+ALGBRC(5)
      ALGBRC(6) = ( CONSTS(25)*STATES(6))/(CONSTS(20)+STATES(6)+( STATES(7)*CONSTS(21)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(8)*CONSTS(22)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(6)*STATES(7))/CONSTS(23)+( STATES(7)*STATES(8)*CONSTS(25))/( CONSTS(19)*CONSTS(24))) - ( CONSTS(25)*(( STATES(7)*STATES(8))/CONSTS(19)))/(CONSTS(20)+STATES(6)+( STATES(7)*CONSTS(21)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(8)*CONSTS(22)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(6)*STATES(7))/CONSTS(23)+( STATES(7)*STATES(8)*CONSTS(25))/( CONSTS(19)*CONSTS(24)))
      RATES(6) = ALGBRC(5)+ALGBRC(6)
      ALGBRC(8) = ( CONSTS(35)*STATES(7)*STATES(24))/( CONSTS(32)*CONSTS(30)*(1.00000+STATES(7)/CONSTS(32)+STATES(9)/CONSTS(33))*(1.00000+STATES(24)/CONSTS(30)+STATES(23)/CONSTS(31))) - ( CONSTS(35)*(( STATES(9)*STATES(23))/CONSTS(34)))/( CONSTS(32)*CONSTS(30)*(1.00000+STATES(7)/CONSTS(32)+STATES(9)/CONSTS(33))*(1.00000+STATES(24)/CONSTS(30)+STATES(23)/CONSTS(31)))
      ALGBRC(7) = ( CONSTS(29)*STATES(8))/(CONSTS(27)+STATES(8)+ (CONSTS(27)/CONSTS(28))*STATES(7)) - ( CONSTS(29)*(STATES(7)/CONSTS(26)))/(CONSTS(27)+STATES(8)+ (CONSTS(27)/CONSTS(28))*STATES(7))
      RATES(7) = ALGBRC(6)+ALGBRC(7)+ALGBRC(8)
      ALGBRC(9) =  CONSTS(36)*STATES(9)*STATES(13) -  CONSTS(37)*STATES(10)*STATES(12)
      RATES(9) = ALGBRC(8)+ALGBRC(9)
      ALGBRC(10) = ( CONSTS(40)*STATES(13)*STATES(10))/( (CONSTS(38)+STATES(10))*(CONSTS(39)+STATES(13)))
      RATES(10) = ALGBRC(9)+ALGBRC(10)
      ALGBRC(11) =  CONSTS(55)*STATES(12)*STATES(4)
      RATES(4) = ALGBRC(3)+ALGBRC(4)+ALGBRC(11)
      ALGBRC(12) = ( CONSTS(42)*STATES(11))/(CONSTS(41)+STATES(11))
      RATES(11) = ALGBRC(10)+ALGBRC(12)
      ALGBRC(15) = ( CONSTS(51)*STATES(8))/( CONSTS(50)*(1.00000+ (CONSTS(48)/STATES(23))*(1.00000+STATES(24)/CONSTS(49)))+ STATES(8)*(1.00000+ (CONSTS(47)/STATES(23))*(1.00000+STATES(24)/CONSTS(49))))
      RATES(8) = ALGBRC(6)+ALGBRC(7)+ALGBRC(15)
      ALGBRC(14) =  CONSTS(56)*STATES(12)
      ALGBRC(16) =  CONSTS(57)*STATES(14)*STATES(12) -  CONSTS(58)*STATES(13) ** 2.00000
      RATES(12) = ALGBRC(5)+ALGBRC(10)+ALGBRC(9)+ALGBRC(11)+ALGBRC(3)+ALGBRC(14)+ALGBRC(16)
      RATES(13) = ALGBRC(5)+ALGBRC(10)+ALGBRC(9)+ALGBRC(11)+ALGBRC(3)+ALGBRC(14)+ALGBRC(16)
      RATES(14) = ALGBRC(16)
      ALGBRC(13) = ( CONSTS(45)*STATES(17)*STATES(23))/( (CONSTS(43)+STATES(23))*(CONSTS(44)+STATES(17)))
      RATES(23) = ALGBRC(8)+ALGBRC(13)+ALGBRC(15)
      RATES(24) = ALGBRC(8)+ALGBRC(13)+ALGBRC(15)
      ALGBRC(19) =  (CONSTS(53)/CONSTS(2))*(STATES(17) - STATES(18))
      RATES(17) = ALGBRC(19)+ALGBRC(12)+ALGBRC(13)
      ALGBRC(20) =  (CONSTS(46)/CONSTS(2))*(STATES(19) - STATES(20))
      RATES(19) = ALGBRC(20)+ALGBRC(13)
      ALGBRC(18) =  CONSTS(1)*STATES(20)
      RATES(20) = ALGBRC(20)+ALGBRC(18)
      ALGBRC(22) =  CONSTS(54)*STATES(18)*STATES(16)
      ALGBRC(17) =  CONSTS(1)*STATES(18)
      RATES(18) = ALGBRC(22)+ALGBRC(19)+ALGBRC(17)
      ALGBRC(23) =  (CONSTS(52)/CONSTS(2))*(STATES(21) - STATES(22))
      RATES(21) = ALGBRC(23)+ALGBRC(15)
      ALGBRC(21) =  CONSTS(1)*STATES(22)
      RATES(22) = ALGBRC(23)+ALGBRC(21)
      ALGBRC(24) =  CONSTS(1)*(STATES(15) - STATES(16))
      RATES(15) = ALGBRC(24)
      RATES(16) = ALGBRC(22)+ALGBRC(24)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*(STATES(1) - STATES(2))
      ALGBRC(2) =  (CONSTS(6)/CONSTS(2))*((STATES(2)/CONSTS(3))/(1.00000+STATES(2)/CONSTS(3)+ (( CONSTS(7)*(STATES(2)/CONSTS(3))+1.00000)/( CONSTS(7)*(STATES(3)/CONSTS(3))+1.00000))*(1.00000+STATES(3)/CONSTS(3))+STATES(4)/CONSTS(4)+( STATES(3)*STATES(4))/( CONSTS(3)*CONSTS(5)))) -  (CONSTS(6)/CONSTS(2))*((STATES(3)/CONSTS(3))/(1.00000+STATES(3)/CONSTS(3)+ (( CONSTS(7)*(STATES(3)/CONSTS(3))+1.00000)/( CONSTS(7)*(STATES(2)/CONSTS(3))+1.00000))*(1.00000+STATES(2)/CONSTS(3))+STATES(4)/CONSTS(4)+( STATES(3)*STATES(4))/( CONSTS(3)*CONSTS(5))))
      ALGBRC(3) = ( CONSTS(11)*STATES(12)*STATES(3))/( CONSTS(8)*CONSTS(10)+ CONSTS(9)*STATES(12)+ CONSTS(10)*STATES(3)+ STATES(3)*STATES(12))
      ALGBRC(4) = ( CONSTS(12)*STATES(4))/(CONSTS(13)+STATES(4)+ (CONSTS(13)/CONSTS(14))*STATES(5)) - ( CONSTS(12)*(STATES(5)/CONSTS(15)))/(CONSTS(13)+STATES(4)+ (CONSTS(13)/CONSTS(14))*STATES(5))
      ALGBRC(5) = ( CONSTS(18)*STATES(5) ** 2.00000)/( CONSTS(16)*(1.00000+ CONSTS(17)*(STATES(12)/STATES(14)))+STATES(5) ** 2.00000)
      ALGBRC(6) = ( CONSTS(25)*STATES(6))/(CONSTS(20)+STATES(6)+( STATES(7)*CONSTS(21)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(8)*CONSTS(22)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(6)*STATES(7))/CONSTS(23)+( STATES(7)*STATES(8)*CONSTS(25))/( CONSTS(19)*CONSTS(24))) - ( CONSTS(25)*(( STATES(7)*STATES(8))/CONSTS(19)))/(CONSTS(20)+STATES(6)+( STATES(7)*CONSTS(21)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(8)*CONSTS(22)*CONSTS(25))/( CONSTS(19)*CONSTS(24))+( STATES(6)*STATES(7))/CONSTS(23)+( STATES(7)*STATES(8)*CONSTS(25))/( CONSTS(19)*CONSTS(24)))
      ALGBRC(8) = ( CONSTS(35)*STATES(7)*STATES(24))/( CONSTS(32)*CONSTS(30)*(1.00000+STATES(7)/CONSTS(32)+STATES(9)/CONSTS(33))*(1.00000+STATES(24)/CONSTS(30)+STATES(23)/CONSTS(31))) - ( CONSTS(35)*(( STATES(9)*STATES(23))/CONSTS(34)))/( CONSTS(32)*CONSTS(30)*(1.00000+STATES(7)/CONSTS(32)+STATES(9)/CONSTS(33))*(1.00000+STATES(24)/CONSTS(30)+STATES(23)/CONSTS(31)))
      ALGBRC(7) = ( CONSTS(29)*STATES(8))/(CONSTS(27)+STATES(8)+ (CONSTS(27)/CONSTS(28))*STATES(7)) - ( CONSTS(29)*(STATES(7)/CONSTS(26)))/(CONSTS(27)+STATES(8)+ (CONSTS(27)/CONSTS(28))*STATES(7))
      ALGBRC(9) =  CONSTS(36)*STATES(9)*STATES(13) -  CONSTS(37)*STATES(10)*STATES(12)
      ALGBRC(10) = ( CONSTS(40)*STATES(13)*STATES(10))/( (CONSTS(38)+STATES(10))*(CONSTS(39)+STATES(13)))
      ALGBRC(11) =  CONSTS(55)*STATES(12)*STATES(4)
      ALGBRC(12) = ( CONSTS(42)*STATES(11))/(CONSTS(41)+STATES(11))
      ALGBRC(15) = ( CONSTS(51)*STATES(8))/( CONSTS(50)*(1.00000+ (CONSTS(48)/STATES(23))*(1.00000+STATES(24)/CONSTS(49)))+ STATES(8)*(1.00000+ (CONSTS(47)/STATES(23))*(1.00000+STATES(24)/CONSTS(49))))
      ALGBRC(14) =  CONSTS(56)*STATES(12)
      ALGBRC(16) =  CONSTS(57)*STATES(14)*STATES(12) -  CONSTS(58)*STATES(13) ** 2.00000
      ALGBRC(13) = ( CONSTS(45)*STATES(17)*STATES(23))/( (CONSTS(43)+STATES(23))*(CONSTS(44)+STATES(17)))
      ALGBRC(19) =  (CONSTS(53)/CONSTS(2))*(STATES(17) - STATES(18))
      ALGBRC(20) =  (CONSTS(46)/CONSTS(2))*(STATES(19) - STATES(20))
      ALGBRC(18) =  CONSTS(1)*STATES(20)
      ALGBRC(22) =  CONSTS(54)*STATES(18)*STATES(16)
      ALGBRC(17) =  CONSTS(1)*STATES(18)
      ALGBRC(23) =  (CONSTS(52)/CONSTS(2))*(STATES(21) - STATES(22))
      ALGBRC(21) =  CONSTS(1)*STATES(22)
      ALGBRC(24) =  CONSTS(1)*(STATES(15) - STATES(16))
      RETURN
      END