C
C There are a total of 33 entries in the algebraic variable array.
C There are a total of 10 entries in each of the rate and state variable arrays.
C There are a total of 61 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is Vc_Vg in component volume_ratio (dimensionless).
C CONSTS(2) is Glc_o in component Glc_o (millimolar).
C STATES(1) is Glc_i in component Glc_i (millimolar).
C ALGBRC(23) is V_HK in component V_HK (flux).
C ALGBRC(20) is V_glucose_transport in component V_glucose_transport (flux).
C ALGBRC(1) is Glc_6_P_g in component Glc_6_P_g (millimolar).
C CONSTS(3) is Keq_PGI in component Glc_6_P_g (dimensionless).
C ALGBRC(2) is Fru_6_P_g in component Fru_6_P_g (millimolar).
C STATES(2) is hexose_P_g in component hexose_P_g (millimolar).
C ALGBRC(30) is V_PFK in component V_PFK (flux).
C STATES(3) is Fru_1_6_BP_g in component Fru_1_6_BP_g (millimolar).
C ALGBRC(33) is V_ALD in component V_ALD (flux).
C ALGBRC(9) is GA_3_P_g in component GA_3_P_g (millimolar).
C CONSTS(4) is Keq_TIM in component DHAP (dimensionless).
C ALGBRC(7) is DHAP_g in component DHAP_g (millimolar).
C STATES(4) is triose_P in component triose_P (millimolar).
C ALGBRC(24) is V_GAPDH in component V_GAPDH (flux).
C ALGBRC(28) is V_GDH in component V_GDH (flux).
C ALGBRC(22) is V_GPO in component V_GPO (flux).
C STATES(5) is one_three_BPGA_g in component one_three_BPGA_g (millimolar).
C ALGBRC(25) is V_PGK in component V_PGK (flux).
C ALGBRC(3) is three_PGA in component three_PGA (millimolar).
C STATES(6) is N in component N (millimolar).
C CONSTS(5) is Keq_PGM in component three_PGA (dimensionless).
C CONSTS(6) is Keq_ENO in component three_PGA (dimensionless).
C ALGBRC(4) is two_PGA_c in component two_PGA_c (millimolar).
C ALGBRC(29) is V_PYK in component V_PYK (flux).
C ALGBRC(5) is PEP_c in component PEP_c (millimolar).
C STATES(7) is PYR_c in component PYR_c (millimolar).
C ALGBRC(21) is V_pyruvate_transport in component V_pyruvate_transport (flux).
C CONSTS(7) is glycerol_g in component glycerol_g (millimolar).
C ALGBRC(6) is DHAP in component DHAP (millimolar).
C ALGBRC(8) is DHAP_c in component DHAP_c (millimolar).
C ALGBRC(14) is Gly_3_P in component Gly_3_P (millimolar).
C ALGBRC(11) is ATP_g in component ATP_g (millimolar).
C ALGBRC(13) is ADP_g in component ADP_g (millimolar).
C CONSTS(8) is C4 in component C4 (millimolar).
C ALGBRC(15) is Gly_3_P_c in component Gly_3_P_c (millimolar).
C ALGBRC(16) is Gly_3_P_g in component Gly_3_P_g (millimolar).
C ALGBRC(10) is NAD_g in component NAD_g (millimolar).
C STATES(8) is NADH_g in component NADH_g (millimolar).
C CONSTS(9) is C3 in component C3 (millimolar).
C STATES(9) is P_g in component P_g (millimolar).
C ALGBRC(26) is V_GK in component V_GK (flux).
C STATES(10) is P_c in component P_c (millimolar).
C ALGBRC(31) is V_ATP_utilisation in component V_ATP_utilisation (flux).
C CONSTS(10) is Keq_AK in component ATP_g (dimensionless).
C CONSTS(11) is C1 in component C1 (millimolar).
C ALGBRC(12) is ATP_c in component ATP_c (millimolar).
C CONSTS(12) is Keq_AK in component ATP_c (dimensionless).
C CONSTS(13) is C2 in component C2 (millimolar).
C ALGBRC(17) is ADP_c in component ADP_c (millimolar).
C ALGBRC(18) is AMP_g in component AMP_g (millimolar).
C ALGBRC(19) is AMP_c in component AMP_c (millimolar).
C CONSTS(14) is K_Glc in component V_glucose_transport (millimolar).
C CONSTS(15) is alpha in component V_glucose_transport (dimensionless).
C CONSTS(16) is V_glucose_transport_max in component V_glucose_transport (flux).
C CONSTS(17) is K_pyruvate in component V_pyruvate_transport (millimolar).
C CONSTS(18) is V_pyruvate_transport_max in component V_pyruvate_transport (flux).
C CONSTS(19) is K_Gly_3_P in component V_GPO (millimolar).
C CONSTS(20) is V_GPO_max in component V_GPO (flux).
C CONSTS(21) is K_Glc_i in component V_HK (millimolar).
C CONSTS(22) is K_ATP in component V_HK (millimolar).
C CONSTS(23) is K_ADP in component V_HK (millimolar).
C CONSTS(24) is V_HK_max in component V_HK (flux).
C CONSTS(25) is K_NAD in component V_GAPDH (millimolar).
C CONSTS(26) is K_GA_3_P in component V_GAPDH (millimolar).
C CONSTS(27) is K_1_3_BPGA in component V_GAPDH (millimolar).
C CONSTS(28) is K_NADH in component V_GAPDH (millimolar).
C CONSTS(29) is V_GAPDH_max_plus in component V_GAPDH (flux).
C CONSTS(30) is V_GAPDH_max_ratio in component V_GAPDH (dimensionless).
C CONSTS(31) is K_ADP in component V_PGK (millimolar).
C CONSTS(32) is K_1_3_BPGA in component V_PGK (millimolar).
C CONSTS(33) is K_3_PGA in component V_PGK (millimolar).
C CONSTS(34) is K_ATP in component V_PGK (millimolar).
C CONSTS(35) is V_PGK_max_plus in component V_PGK (flux).
C CONSTS(36) is V_PGK_max_ratio in component V_PGK (dimensionless).
C CONSTS(37) is K_ADP in component V_GK (millimolar).
C CONSTS(38) is K_Gly_3_P in component V_GK (millimolar).
C CONSTS(39) is K_glycerol in component V_GK (millimolar).
C CONSTS(40) is K_ATP in component V_GK (millimolar).
C CONSTS(41) is V_GK_max_plus in component V_GK (flux).
C CONSTS(42) is V_GK_max_ratio in component V_GK (dimensionless).
C CONSTS(43) is K_NADH in component V_GDH (millimolar).
C CONSTS(44) is K_Gly_3_P in component V_GDH (millimolar).
C CONSTS(45) is K_DHAP in component V_GDH (millimolar).
C CONSTS(46) is K_NAD in component V_GDH (millimolar).
C CONSTS(47) is V_GDH_max_plus in component V_GDH (flux).
C CONSTS(48) is V_GDH_max_ratio in component V_GDH (dimensionless).
C CONSTS(49) is n in component V_PFK (dimensionless).
C CONSTS(50) is Km_Fru_6_P in component V_PFK (millimolar).
C CONSTS(51) is Km_ATP in component V_PFK (millimolar).
C CONSTS(52) is V_PFK_max in component V_PFK (flux).
C ALGBRC(27) is Km_PEP in component V_PYK (millimolar).
C CONSTS(53) is Km_ADP in component V_PYK (millimolar).
C CONSTS(54) is n in component V_PYK (dimensionless).
C CONSTS(55) is V_PYK_max in component V_PYK (flux).
C ALGBRC(32) is Km_Fru_1_6_BP in component V_ALD (millimolar).
C CONSTS(56) is Km_GA_3_P in component V_ALD (millimolar).
C CONSTS(57) is Ki_GA_3_P in component V_ALD (millimolar).
C CONSTS(58) is Km_DHAP in component V_ALD (millimolar).
C CONSTS(59) is V_ALD_max_plus in component V_ALD (flux).
C CONSTS(60) is V_ALD_max_ratio in component V_ALD (dimensionless).
C CONSTS(61) is k in component V_ATP_utilisation (flux).
C RATES(1) is d/dt Glc_i in component Glc_i (millimolar).
C RATES(2) is d/dt hexose_P_g in component hexose_P_g (millimolar).
C RATES(3) is d/dt Fru_1_6_BP_g in component Fru_1_6_BP_g (millimolar).
C RATES(4) is d/dt triose_P in component triose_P (millimolar).
C RATES(5) is d/dt one_three_BPGA_g in component one_three_BPGA_g (millimolar).
C RATES(6) is d/dt N in component N (millimolar).
C RATES(7) is d/dt PYR_c in component PYR_c (millimolar).
C RATES(8) is d/dt NADH_g in component NADH_g (millimolar).
C RATES(9) is d/dt P_g in component P_g (millimolar).
C RATES(10) is d/dt P_c in component P_c (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 22.3
      CONSTS(2) = 5.0
      STATES(1) = 0.01
      CONSTS(3) = 0.29
      STATES(2) = 0.01
      STATES(3) = 0.01
      CONSTS(4) = 0.045
      STATES(4) = 0.01
      STATES(5) = 0.01
      STATES(6) = 0.01
      CONSTS(5) = 0.187
      CONSTS(6) = 6.7
      STATES(7) = 0.01
      CONSTS(7) = 0.00
      CONSTS(8) = 120.0
      STATES(8) = 0.01
      CONSTS(9) = 4.0
      STATES(9) = 0.01
      STATES(10) = 0.01
      CONSTS(10) = 0.442
      CONSTS(11) = 3.9
      CONSTS(12) = 0.442
      CONSTS(13) = 3.9
      CONSTS(14) = 2.0
      CONSTS(15) = 0.75
      CONSTS(16) = 106.2
      CONSTS(17) = 1.96
      CONSTS(18) = 160.0
      CONSTS(19) = 1.7
      CONSTS(20) = 200.0
      CONSTS(21) = 0.1
      CONSTS(22) = 0.116
      CONSTS(23) = 0.126
      CONSTS(24) = 625.0
      CONSTS(25) = 0.45
      CONSTS(26) = 0.15
      CONSTS(27) = 0.1
      CONSTS(28) = 0.02
      CONSTS(29) = 1470.0
      CONSTS(30) = 0.67
      CONSTS(31) = 0.1
      CONSTS(32) = 0.05
      CONSTS(33) = 1.62
      CONSTS(34) = 0.29
      CONSTS(35) = 640.0
      CONSTS(36) = 0.029
      CONSTS(37) = 0.12
      CONSTS(38) = 5.1
      CONSTS(39) = 0.12
      CONSTS(40) = 0.19
      CONSTS(41) = 0.0
      CONSTS(42) = 167.0
      CONSTS(43) = 0.01
      CONSTS(44) = 2.0
      CONSTS(45) = 0.1
      CONSTS(46) = 0.4
      CONSTS(47) = 533.0
      CONSTS(48) = 0.28
      CONSTS(49) = 1.2
      CONSTS(50) = 0.82
      CONSTS(51) = 2.6E-2
      CONSTS(52) = 780.0
      CONSTS(53) = 0.114
      CONSTS(54) = 2.5
      CONSTS(55) = 2.6E3
      CONSTS(56) = 6.7E-2
      CONSTS(57) = 9.8E-2
      CONSTS(58) = 1.5E-2
      CONSTS(59) = 780.0
      CONSTS(60) = 1.19
      CONSTS(61) = 50
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(11) = ((CONSTS(11) - - ( STATES(9)*(1.00000 -  4.00000*CONSTS(10))))+CONSTS(11) - - ( STATES(9)*(1.00000 -  4.00000*CONSTS(10))) ** 2.00000 -  4.00000*(1.00000 -  4.00000*CONSTS(10))*( - CONSTS(10)*STATES(9) ** 2.00000) ** 0.500000)/( 2.00000*(1.00000 -  4.00000*CONSTS(10)))
      ALGBRC(13) = STATES(9) -  2.00000*ALGBRC(11)
      ALGBRC(23) =  CONSTS(24)*(( (ALGBRC(11)/CONSTS(22))*(STATES(1)/CONSTS(21)))/( (1.00000+ALGBRC(11)/CONSTS(22)+ALGBRC(13)/CONSTS(23))*(1.00000+STATES(1)/CONSTS(21))))
      ALGBRC(20) =  CONSTS(16)*((CONSTS(2) - STATES(1))/(CONSTS(14)+CONSTS(2)+STATES(1)+ CONSTS(15)*CONSTS(2)*(STATES(1)/CONSTS(14))))
      RATES(1) = ALGBRC(20) - ALGBRC(23)
      ALGBRC(6) = ( STATES(4)*(1.00000+CONSTS(1)))/(1.00000+CONSTS(1)+CONSTS(4))
      ALGBRC(7) = ALGBRC(6)
      ALGBRC(9) =  CONSTS(4)*ALGBRC(7)
      ALGBRC(10) = CONSTS(9) - STATES(8)
      ALGBRC(24) =  CONSTS(29)*(( (ALGBRC(9)/CONSTS(26))*(ALGBRC(10)/CONSTS(25) - CONSTS(30))*(STATES(5)/CONSTS(27))*(STATES(8)/CONSTS(28)))/( (1.00000+ALGBRC(9)/CONSTS(26)+STATES(5)/CONSTS(27))*(1.00000+ALGBRC(10)/CONSTS(25)+STATES(8)/CONSTS(28))))
      ALGBRC(3) = ( STATES(6)*(1.00000+CONSTS(1)))/(1.00000+ (1.00000+CONSTS(5)+ CONSTS(5)*CONSTS(6))*CONSTS(1))
      ALGBRC(25) =  CONSTS(35)*(( (STATES(5)/CONSTS(32))*(ALGBRC(13)/CONSTS(31) - CONSTS(36))*(ALGBRC(3)/CONSTS(33))*(ALGBRC(11)/CONSTS(34)))/( (1.00000+STATES(5)/CONSTS(32)+ALGBRC(3)/CONSTS(33))*(1.00000+ALGBRC(13)/CONSTS(31)+ALGBRC(11)/CONSTS(34))))
      RATES(5) = ALGBRC(24) - ALGBRC(25)
      ALGBRC(4) =  CONSTS(5)*ALGBRC(3)
      ALGBRC(5) =  CONSTS(6)*ALGBRC(4)
      ALGBRC(12) = ((CONSTS(13) - - ( STATES(10)*(1.00000 -  4.00000*CONSTS(12))))+CONSTS(13) - - ( STATES(10)*(1.00000 -  4.00000*CONSTS(12))) ** 2.00000 -  4.00000*(1.00000 -  4.00000*CONSTS(12))*( - CONSTS(12)*STATES(10) ** 2.00000) ** 0.500000)/( 2.00000*(1.00000 -  4.00000*CONSTS(12)))
      ALGBRC(17) = STATES(10) -  2.00000*ALGBRC(12)
      ALGBRC(27) =  0.340000*(1.00000+ALGBRC(12)/0.570000+ALGBRC(17)/0.640000)
      ALGBRC(29) =  CONSTS(55)*(( ALGBRC(5)/ALGBRC(27) ** CONSTS(54)*(ALGBRC(17)/CONSTS(53)))/( (1.00000+ALGBRC(5)/ALGBRC(27) ** CONSTS(54))*(1.00000+ALGBRC(17)/CONSTS(53))))
      RATES(6) = ALGBRC(25) - ALGBRC(29)
      ALGBRC(21) =  CONSTS(18)*(STATES(7)/CONSTS(17))*(1.00000+STATES(7)/CONSTS(17))
      RATES(7) = ALGBRC(29) - ALGBRC(21)
      ALGBRC(1) = STATES(2)/CONSTS(3)
      ALGBRC(2) = STATES(2) - ALGBRC(1)
      ALGBRC(14) = (CONSTS(8) - (ALGBRC(1)+ALGBRC(2)+ 2.00000*STATES(3)+ALGBRC(9)+STATES(5)+ 2.00000*ALGBRC(11)+ALGBRC(13)))/(1.00000+CONSTS(1)) - ALGBRC(6)
      ALGBRC(16) = ALGBRC(14)
      ALGBRC(28) =  CONSTS(47)*(( (ALGBRC(7)/CONSTS(45))*(STATES(8)/CONSTS(43) - CONSTS(48))*(ALGBRC(16)/CONSTS(44))*(ALGBRC(10)/CONSTS(46)))/( (1.00000+ALGBRC(7)/CONSTS(45)+ALGBRC(16)/CONSTS(44))*(1.00000+STATES(8)/CONSTS(43)+ALGBRC(10)/CONSTS(46))))
      RATES(8) = ALGBRC(24) - ALGBRC(28)
      ALGBRC(30) =  CONSTS(52)*(( ALGBRC(2)/CONSTS(50) ** CONSTS(49)*(ALGBRC(11)/CONSTS(51)))/( (1.00000+ALGBRC(2)/CONSTS(50) ** CONSTS(49))*(1.00000+ALGBRC(11)/CONSTS(51))))
      RATES(2) = ALGBRC(23) - ALGBRC(30)
      ALGBRC(26) =  CONSTS(41)*(( (ALGBRC(16)/CONSTS(38))*(ALGBRC(13)/CONSTS(37) - CONSTS(42))*(CONSTS(7)/CONSTS(39))*(ALGBRC(11)/CONSTS(40)))/( (1.00000+ALGBRC(16)/CONSTS(38)+CONSTS(7)/CONSTS(39))*(1.00000+ALGBRC(13)/CONSTS(37)+ALGBRC(11)/CONSTS(40))))
      RATES(9) = (ALGBRC(25)+ALGBRC(26)) - (ALGBRC(23)+ALGBRC(30))
      ALGBRC(31) =  CONSTS(61)*(ALGBRC(12)/ALGBRC(17))
      RATES(10) = ALGBRC(29) - ALGBRC(31)
      ALGBRC(18) = CONSTS(11) - (ALGBRC(11)+ALGBRC(13))
      ALGBRC(32) =  0.00900000*(1.00000+ALGBRC(11)/0.680000+ALGBRC(13)/1.51000+ALGBRC(18)/3.65000)
      ALGBRC(33) =  CONSTS(59)*((STATES(3)/ALGBRC(32) -  CONSTS(60)*(( ALGBRC(9)*ALGBRC(7))/( CONSTS(56)*CONSTS(58))))/(1.00000+STATES(3)/ALGBRC(32)+ALGBRC(9)/CONSTS(56)+ALGBRC(7)/CONSTS(58)+( STATES(3)*ALGBRC(9))/( ALGBRC(32)*CONSTS(57))+( ALGBRC(7)*ALGBRC(9))/( CONSTS(58)*CONSTS(56))))
      RATES(3) = ALGBRC(30) - ALGBRC(33)
      ALGBRC(15) = ALGBRC(14)
      ALGBRC(22) =  CONSTS(20)*(ALGBRC(15)/CONSTS(19))*(1.00000+ALGBRC(15)/CONSTS(19))
      RATES(4) = ( 2.00000*ALGBRC(33)+ALGBRC(22)) - (ALGBRC(24)+ALGBRC(28))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(11) = ((CONSTS(11) - - ( STATES(9)*(1.00000 -  4.00000*CONSTS(10))))+CONSTS(11) - - ( STATES(9)*(1.00000 -  4.00000*CONSTS(10))) ** 2.00000 -  4.00000*(1.00000 -  4.00000*CONSTS(10))*( - CONSTS(10)*STATES(9) ** 2.00000) ** 0.500000)/( 2.00000*(1.00000 -  4.00000*CONSTS(10)))
      ALGBRC(13) = STATES(9) -  2.00000*ALGBRC(11)
      ALGBRC(23) =  CONSTS(24)*(( (ALGBRC(11)/CONSTS(22))*(STATES(1)/CONSTS(21)))/( (1.00000+ALGBRC(11)/CONSTS(22)+ALGBRC(13)/CONSTS(23))*(1.00000+STATES(1)/CONSTS(21))))
      ALGBRC(20) =  CONSTS(16)*((CONSTS(2) - STATES(1))/(CONSTS(14)+CONSTS(2)+STATES(1)+ CONSTS(15)*CONSTS(2)*(STATES(1)/CONSTS(14))))
      ALGBRC(6) = ( STATES(4)*(1.00000+CONSTS(1)))/(1.00000+CONSTS(1)+CONSTS(4))
      ALGBRC(7) = ALGBRC(6)
      ALGBRC(9) =  CONSTS(4)*ALGBRC(7)
      ALGBRC(10) = CONSTS(9) - STATES(8)
      ALGBRC(24) =  CONSTS(29)*(( (ALGBRC(9)/CONSTS(26))*(ALGBRC(10)/CONSTS(25) - CONSTS(30))*(STATES(5)/CONSTS(27))*(STATES(8)/CONSTS(28)))/( (1.00000+ALGBRC(9)/CONSTS(26)+STATES(5)/CONSTS(27))*(1.00000+ALGBRC(10)/CONSTS(25)+STATES(8)/CONSTS(28))))
      ALGBRC(3) = ( STATES(6)*(1.00000+CONSTS(1)))/(1.00000+ (1.00000+CONSTS(5)+ CONSTS(5)*CONSTS(6))*CONSTS(1))
      ALGBRC(25) =  CONSTS(35)*(( (STATES(5)/CONSTS(32))*(ALGBRC(13)/CONSTS(31) - CONSTS(36))*(ALGBRC(3)/CONSTS(33))*(ALGBRC(11)/CONSTS(34)))/( (1.00000+STATES(5)/CONSTS(32)+ALGBRC(3)/CONSTS(33))*(1.00000+ALGBRC(13)/CONSTS(31)+ALGBRC(11)/CONSTS(34))))
      ALGBRC(4) =  CONSTS(5)*ALGBRC(3)
      ALGBRC(5) =  CONSTS(6)*ALGBRC(4)
      ALGBRC(12) = ((CONSTS(13) - - ( STATES(10)*(1.00000 -  4.00000*CONSTS(12))))+CONSTS(13) - - ( STATES(10)*(1.00000 -  4.00000*CONSTS(12))) ** 2.00000 -  4.00000*(1.00000 -  4.00000*CONSTS(12))*( - CONSTS(12)*STATES(10) ** 2.00000) ** 0.500000)/( 2.00000*(1.00000 -  4.00000*CONSTS(12)))
      ALGBRC(17) = STATES(10) -  2.00000*ALGBRC(12)
      ALGBRC(27) =  0.340000*(1.00000+ALGBRC(12)/0.570000+ALGBRC(17)/0.640000)
      ALGBRC(29) =  CONSTS(55)*(( ALGBRC(5)/ALGBRC(27) ** CONSTS(54)*(ALGBRC(17)/CONSTS(53)))/( (1.00000+ALGBRC(5)/ALGBRC(27) ** CONSTS(54))*(1.00000+ALGBRC(17)/CONSTS(53))))
      ALGBRC(21) =  CONSTS(18)*(STATES(7)/CONSTS(17))*(1.00000+STATES(7)/CONSTS(17))
      ALGBRC(1) = STATES(2)/CONSTS(3)
      ALGBRC(2) = STATES(2) - ALGBRC(1)
      ALGBRC(14) = (CONSTS(8) - (ALGBRC(1)+ALGBRC(2)+ 2.00000*STATES(3)+ALGBRC(9)+STATES(5)+ 2.00000*ALGBRC(11)+ALGBRC(13)))/(1.00000+CONSTS(1)) - ALGBRC(6)
      ALGBRC(16) = ALGBRC(14)
      ALGBRC(28) =  CONSTS(47)*(( (ALGBRC(7)/CONSTS(45))*(STATES(8)/CONSTS(43) - CONSTS(48))*(ALGBRC(16)/CONSTS(44))*(ALGBRC(10)/CONSTS(46)))/( (1.00000+ALGBRC(7)/CONSTS(45)+ALGBRC(16)/CONSTS(44))*(1.00000+STATES(8)/CONSTS(43)+ALGBRC(10)/CONSTS(46))))
      ALGBRC(30) =  CONSTS(52)*(( ALGBRC(2)/CONSTS(50) ** CONSTS(49)*(ALGBRC(11)/CONSTS(51)))/( (1.00000+ALGBRC(2)/CONSTS(50) ** CONSTS(49))*(1.00000+ALGBRC(11)/CONSTS(51))))
      ALGBRC(26) =  CONSTS(41)*(( (ALGBRC(16)/CONSTS(38))*(ALGBRC(13)/CONSTS(37) - CONSTS(42))*(CONSTS(7)/CONSTS(39))*(ALGBRC(11)/CONSTS(40)))/( (1.00000+ALGBRC(16)/CONSTS(38)+CONSTS(7)/CONSTS(39))*(1.00000+ALGBRC(13)/CONSTS(37)+ALGBRC(11)/CONSTS(40))))
      ALGBRC(31) =  CONSTS(61)*(ALGBRC(12)/ALGBRC(17))
      ALGBRC(18) = CONSTS(11) - (ALGBRC(11)+ALGBRC(13))
      ALGBRC(32) =  0.00900000*(1.00000+ALGBRC(11)/0.680000+ALGBRC(13)/1.51000+ALGBRC(18)/3.65000)
      ALGBRC(33) =  CONSTS(59)*((STATES(3)/ALGBRC(32) -  CONSTS(60)*(( ALGBRC(9)*ALGBRC(7))/( CONSTS(56)*CONSTS(58))))/(1.00000+STATES(3)/ALGBRC(32)+ALGBRC(9)/CONSTS(56)+ALGBRC(7)/CONSTS(58)+( STATES(3)*ALGBRC(9))/( ALGBRC(32)*CONSTS(57))+( ALGBRC(7)*ALGBRC(9))/( CONSTS(58)*CONSTS(56))))
      ALGBRC(15) = ALGBRC(14)
      ALGBRC(22) =  CONSTS(20)*(ALGBRC(15)/CONSTS(19))*(1.00000+ALGBRC(15)/CONSTS(19))
      ALGBRC(8) = ALGBRC(6)
      ALGBRC(19) = CONSTS(13) - (ALGBRC(12)+ALGBRC(17))
      RETURN
      END