C
C There are a total of 23 entries in the algebraic variable array.
C There are a total of 14 entries in each of the rate and state variable arrays.
C There are a total of 90 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is GLCo in component GLCo (millimolar).
C STATES(1) is GLCi in component GLCi (millimolar).
C ALGBRC(10) is v_GLT in component v_GLT (flux).
C ALGBRC(9) is v_GLK in component v_GLK (flux).
C STATES(2) is G6P in component G6P (millimolar).
C ALGBRC(6) is v_PGI in component v_PGI (flux).
C CONSTS(2) is v_trehalose in component v_trehalose (flux).
C CONSTS(3) is v_glycogen in component v_glycogen (flux).
C STATES(3) is F6P in component F6P (millimolar).
C ALGBRC(21) is v_PFK in component v_PFK (flux).
C STATES(4) is F16P in component F16P (millimolar).
C ALGBRC(14) is v_ALD in component v_ALD (flux).
C CONSTS(4) is F26P in component F26P (millimolar).
C STATES(5) is TRIO in component TRIO (millimolar).
C ALGBRC(11) is v_GAPDH in component v_GAPDH (flux).
C ALGBRC(22) is v_G3PDH in component v_G3PDH (flux).
C STATES(6) is BPG in component BPG (millimolar).
C ALGBRC(12) is v_PGK in component v_PGK (flux).
C STATES(7) is P3G in component P3G (millimolar).
C ALGBRC(7) is v_PGM in component v_PGM (flux).
C STATES(8) is P2G in component P2G (millimolar).
C ALGBRC(8) is v_ENO in component v_ENO (flux).
C STATES(9) is PEP in component PEP (millimolar).
C ALGBRC(13) is v_PYK in component v_PYK (flux).
C STATES(10) is PYR in component PYR (millimolar).
C ALGBRC(15) is v_PDC in component v_PDC (flux).
C STATES(11) is ACE in component ACE (millimolar).
C ALGBRC(16) is v_ADH in component v_ADH (flux).
C ALGBRC(23) is v_succinate in component v_succinate (flux).
C STATES(12) is P in component P (millimolar).
C ALGBRC(17) is v_ATP in component v_ATP (flux).
C STATES(13) is NADH in component NADH (millimolar).
C STATES(14) is NAD in component NAD (millimolar).
C ALGBRC(1) is DHAP in component DHAP (millimolar).
C CONSTS(5) is Keq_TPI in component DHAP (dimensionless).
C ALGBRC(2) is GAP in component GAP (millimolar).
C ALGBRC(3) is ADP in component ADP (millimolar).
C CONSTS(6) is SUM_P in component ADP (millimolar).
C CONSTS(7) is Keq_AK in component ADP (dimensionless).
C ALGBRC(4) is ATP in component ATP (millimolar).
C ALGBRC(5) is AMP in component AMP (millimolar).
C CONSTS(8) is ETOH in component ETOH (millimolar).
C CONSTS(9) is GLY in component GLY (millimolar).
C CONSTS(10) is Km_G6P in component v_PGI (millimolar).
C CONSTS(11) is Km_F6P in component v_PGI (millimolar).
C CONSTS(12) is Keq_PGI in component v_PGI (dimensionless).
C CONSTS(13) is v_max in component v_PGI (flux).
C CONSTS(14) is Km_P3G in component v_PGM (millimolar).
C CONSTS(15) is Km_P2G in component v_PGM (millimolar).
C CONSTS(16) is Keq_PGM in component v_PGM (dimensionless).
C CONSTS(17) is v_max in component v_PGM (flux).
C CONSTS(18) is Km_P2G in component v_ENO (millimolar).
C CONSTS(19) is Km_PEP in component v_ENO (millimolar).
C CONSTS(20) is Keq_ENO in component v_ENO (dimensionless).
C CONSTS(21) is v_max in component v_ENO (flux).
C CONSTS(22) is Km_GLCi in component v_GLK (millimolar).
C CONSTS(23) is Km_G6P in component v_GLK (millimolar).
C CONSTS(24) is Km_ATP in component v_GLK (millimolar).
C CONSTS(25) is Km_ADP in component v_GLK (millimolar).
C CONSTS(26) is Keq_GLK in component v_GLK (dimensionless).
C CONSTS(27) is v_max in component v_GLK (flux).
C CONSTS(28) is Km_GAP in component v_GAPDH (millimolar).
C CONSTS(29) is Km_BPG in component v_GAPDH (millimolar).
C CONSTS(30) is Km_NAD in component v_GAPDH (millimolar).
C CONSTS(31) is Km_NADH in component v_GAPDH (millimolar).
C CONSTS(32) is v_max_f in component v_GAPDH (flux).
C CONSTS(33) is v_max_r in component v_GAPDH (flux).
C CONSTS(34) is Km_BPG in component v_PGK (millimolar).
C CONSTS(35) is Km_P3G in component v_PGK (millimolar).
C CONSTS(36) is Km_ADP in component v_PGK (millimolar).
C CONSTS(37) is Km_ATP in component v_PGK (millimolar).
C CONSTS(38) is Keq_PGK in component v_PGK (dimensionless).
C CONSTS(39) is v_max in component v_PGK (flux).
C CONSTS(40) is Km_PEP in component v_PYK (millimolar).
C CONSTS(41) is Km_PYR in component v_PYK (millimolar).
C CONSTS(42) is Km_ADP in component v_PYK (millimolar).
C CONSTS(43) is Km_ATP in component v_PYK (millimolar).
C CONSTS(44) is Keq_PYK in component v_PYK (dimensionless).
C CONSTS(45) is v_max in component v_PYK (flux).
C CONSTS(46) is Km_Glco in component v_GLT (millimolar).
C CONSTS(47) is Km_Glci in component v_GLT (millimolar).
C CONSTS(48) is Keq_Glc in component v_GLT (dimensionless).
C CONSTS(49) is v_max in component v_GLT (flux).
C CONSTS(50) is Km_F16P in component v_ALD (millimolar).
C CONSTS(51) is Km_GAP in component v_ALD (millimolar).
C CONSTS(52) is Km_GAPi in component v_ALD (millimolar).
C CONSTS(53) is Km_DHAP in component v_ALD (millimolar).
C CONSTS(54) is Keq_ALD in component v_ALD (millimolar).
C CONSTS(55) is v_max in component v_ALD (flux).
C CONSTS(56) is nPDC in component v_PDC (dimensionless).
C CONSTS(57) is Km_PYR in component v_PDC (millimolar).
C CONSTS(58) is v_max in component v_PDC (flux).
C CONSTS(59) is Km_NAD in component v_ADH (millimolar).
C CONSTS(60) is Km_NADH in component v_ADH (millimolar).
C CONSTS(61) is Km_ACE in component v_ADH (millimolar).
C CONSTS(62) is Km_ETOH in component v_ADH (millimolar).
C CONSTS(63) is Ki_NADH in component v_ADH (millimolar).
C CONSTS(64) is Ki_NAD in component v_ADH (millimolar).
C CONSTS(65) is Ki_ACE in component v_ADH (millimolar).
C CONSTS(66) is Ki_ETOH in component v_ADH (millimolar).
C CONSTS(67) is v_max in component v_ADH (flux).
C CONSTS(68) is Keq_ADH in component v_ADH (dimensionless).
C CONSTS(69) is Km_ATPase in component v_ATP (first_order_rate_constant).
C CONSTS(70) is v_max in component v_PFK (flux).
C CONSTS(71) is gR in component v_PFK (dimensionless).
C CONSTS(72) is L0 in component v_PFK (per_millimolar2).
C CONSTS(73) is Km_F6P in component v_PFK (millimolar).
C CONSTS(74) is Km_ATP in component v_PFK (millimolar).
C CONSTS(75) is CATP in component v_PFK (millimolar).
C CONSTS(76) is Km_AMP in component v_PFK (millimolar).
C CONSTS(77) is CAMP in component v_PFK (dimensionless).
C CONSTS(78) is Ki_ATP in component v_PFK (millimolar).
C CONSTS(79) is Ci_ATP in component v_PFK (dimensionless).
C CONSTS(80) is Km_F26BP in component v_PFK (millimolar).
C CONSTS(81) is CF26BP in component v_PFK (dimensionless).
C CONSTS(82) is Km_F16BP in component v_PFK (millimolar).
C CONSTS(83) is CF16BP in component v_PFK (dimensionless).
C ALGBRC(19) is L in component v_PFK (per_millimolar2).
C ALGBRC(18) is R in component v_PFK (dimensionless).
C ALGBRC(20) is T in component v_PFK (millimolar).
C CONSTS(84) is Km_DHAP in component v_G3PDH (millimolar).
C CONSTS(85) is Km_GLY in component v_G3PDH (millimolar).
C CONSTS(86) is Km_NADH in component v_G3PDH (millimolar).
C CONSTS(87) is Km_NAD in component v_G3PDH (millimolar).
C CONSTS(88) is Keq_G3PDH in component v_G3PDH (dimensionless).
C CONSTS(89) is v_max in component v_G3PDH (flux).
C CONSTS(90) is Km_succinate in component v_succinate (first_order_rate_constant).
C RATES(1) is d/dt GLCi in component GLCi (millimolar).
C RATES(2) is d/dt G6P in component G6P (millimolar).
C RATES(3) is d/dt F6P in component F6P (millimolar).
C RATES(4) is d/dt F16P in component F16P (millimolar).
C RATES(5) is d/dt TRIO in component TRIO (millimolar).
C RATES(6) is d/dt BPG in component BPG (millimolar).
C RATES(7) is d/dt P3G in component P3G (millimolar).
C RATES(8) is d/dt P2G in component P2G (millimolar).
C RATES(9) is d/dt PEP in component PEP (millimolar).
C RATES(10) is d/dt PYR in component PYR (millimolar).
C RATES(11) is d/dt ACE in component ACE (millimolar).
C RATES(12) is d/dt P in component P (millimolar).
C RATES(13) is d/dt NADH in component NADH (millimolar).
C RATES(14) is d/dt NAD in component NAD (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 50.0
      STATES(1) = 0.087
      STATES(2) = 2.45
      CONSTS(2) = 2.4
      CONSTS(3) = 6.0
      STATES(3) = 0.62
      STATES(4) = 5.51
      CONSTS(4) = 0.02
      STATES(5) = 0.96
      STATES(6) = 0.0
      STATES(7) = 0.9
      STATES(8) = 0.12
      STATES(9) = 0.07
      STATES(10) = 1.85
      STATES(11) = 0.17
      STATES(12) = 6.31
      STATES(13) = 0.39
      STATES(14) = 1.20
      CONSTS(5) = 0.045
      CONSTS(6) = 4.1
      CONSTS(7) = 0.450
      CONSTS(8) = 50.0
      CONSTS(9) = 0.15
      CONSTS(10) = 1.4
      CONSTS(11) = 0.3
      CONSTS(12) = 0.314
      CONSTS(13) = 339.677
      CONSTS(14) = 1.2
      CONSTS(15) = 0.08
      CONSTS(16) = 0.19
      CONSTS(17) = 2525.81
      CONSTS(18) = 0.04
      CONSTS(19) = 0.5
      CONSTS(20) = 6.7
      CONSTS(21) = 365.806
      CONSTS(22) = 0.08
      CONSTS(23) = 30.0
      CONSTS(24) = 0.15
      CONSTS(25) = 0.23
      CONSTS(26) = 3800
      CONSTS(27) = 226.452
      CONSTS(28) = 0.21
      CONSTS(29) = 0.0098
      CONSTS(30) = 0.09
      CONSTS(31) = 0.06
      CONSTS(32) = 1184.52
      CONSTS(33) = 6549.68
      CONSTS(34) = 0.003
      CONSTS(35) = 0.53
      CONSTS(36) = 0.2
      CONSTS(37) = 0.3
      CONSTS(38) = 3200
      CONSTS(39) = 1306.45
      CONSTS(40) = 0.14
      CONSTS(41) = 21.0
      CONSTS(42) = 0.53
      CONSTS(43) = 1.5
      CONSTS(44) = 6500
      CONSTS(45) = 1088.71
      CONSTS(46) = 1.1918
      CONSTS(47) = 1.1918
      CONSTS(48) = 1
      CONSTS(49) = 97.264
      CONSTS(50) = 0.3
      CONSTS(51) = 2.0
      CONSTS(52) = 10.0
      CONSTS(53) = 2.4
      CONSTS(54) = 0.069
      CONSTS(55) = 322.258
      CONSTS(56) = 1.9
      CONSTS(57) = 4.33
      CONSTS(58) = 174.194
      CONSTS(59) = 0.17
      CONSTS(60) = 0.11
      CONSTS(61) = 1.11
      CONSTS(62) = 17.0
      CONSTS(63) = 0.031
      CONSTS(64) = 0.92
      CONSTS(65) = 1.1
      CONSTS(66) = 90.0
      CONSTS(67) = 810.0
      CONSTS(68) = 0.000069
      CONSTS(69) = 33.7
      CONSTS(70) = 182.903
      CONSTS(71) = 5.12
      CONSTS(72) = 0.66
      CONSTS(73) = 0.1
      CONSTS(74) = 0.71
      CONSTS(75) = 3.0
      CONSTS(76) = 0.0995
      CONSTS(77) = 0.0845
      CONSTS(78) = 0.65
      CONSTS(79) = 100.0
      CONSTS(80) = 0.000682
      CONSTS(81) = 0.0174
      CONSTS(82) = 0.111
      CONSTS(83) = 0.397
      CONSTS(84) = 0.4
      CONSTS(85) = 1
      CONSTS(86) = 0.023
      CONSTS(87) = 0.93
      CONSTS(88) = 4300
      CONSTS(89) = 70.15
      CONSTS(90) = 21.4
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(7) = ( (CONSTS(17)/CONSTS(14))*(STATES(7) - STATES(8)/CONSTS(16)))/(1.00000+STATES(7)/CONSTS(14)+STATES(8)/CONSTS(15))
      ALGBRC(8) = ( (CONSTS(21)/CONSTS(18))*(STATES(8) - STATES(9)/CONSTS(20)))/(1.00000+STATES(8)/CONSTS(18)+STATES(9)/CONSTS(19))
      RATES(8) = ALGBRC(7) - ALGBRC(8)
      ALGBRC(3) = (CONSTS(6) -  ( STATES(12) ** 2.00000*(1.00000 -  4.00000*CONSTS(7))+ 2.00000*CONSTS(6)*STATES(12)*( 4.00000*CONSTS(7) - 1.00000)+CONSTS(6) ** 2.00000) ** (1.0 / 2))/(1.00000 -  4.00000*CONSTS(7))
      ALGBRC(4) = (STATES(12) - ALGBRC(3))/2.00000
      ALGBRC(9) = ( (CONSTS(27)/( CONSTS(22)*CONSTS(24)))*( STATES(1)*ALGBRC(4) - ( STATES(2)*ALGBRC(3))/CONSTS(26)))/( (1.00000+STATES(1)/CONSTS(22)+STATES(2)/CONSTS(23))*(1.00000+ALGBRC(4)/CONSTS(24)+ALGBRC(3)/CONSTS(25)))
      ALGBRC(6) = ( (CONSTS(13)/CONSTS(10))*(STATES(2) - STATES(3)/CONSTS(12)))/(1.00000+STATES(2)/CONSTS(10)+STATES(3)/CONSTS(11))
      RATES(2) = ALGBRC(9) - (ALGBRC(6)+ 2.00000*CONSTS(2)+CONSTS(3))
      ALGBRC(10) = ( CONSTS(49)*(CONSTS(1) - STATES(1)/CONSTS(48)))/( CONSTS(46)*(1.00000+CONSTS(1)/CONSTS(46)+STATES(1)/CONSTS(47)+( 0.910000*CONSTS(1)*STATES(1))/( CONSTS(47)*CONSTS(46))))
      RATES(1) = ALGBRC(10) - ALGBRC(9)
      ALGBRC(11) = (( CONSTS(32)*(CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5)*STATES(14))/( CONSTS(28)*CONSTS(30)) - ( CONSTS(33)*STATES(6)*STATES(13))/( CONSTS(29)*CONSTS(31)))/( (1.00000+( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(28)+STATES(6)/CONSTS(29))*(1.00000+STATES(14)/CONSTS(30)+STATES(13)/CONSTS(31)))
      ALGBRC(12) = ( (CONSTS(39)/( CONSTS(35)*CONSTS(37)))*( CONSTS(38)*STATES(6)*ALGBRC(3) -  STATES(7)*ALGBRC(4)))/( (1.00000+STATES(6)/CONSTS(34)+STATES(7)/CONSTS(35))*(1.00000+ALGBRC(4)/CONSTS(37)+ALGBRC(3)/CONSTS(36)))
      RATES(6) = ALGBRC(11) - ALGBRC(12)
      RATES(7) = ALGBRC(12) - ALGBRC(7)
      ALGBRC(13) = ( (CONSTS(45)/( CONSTS(40)*CONSTS(42)))*( STATES(9)*ALGBRC(3) - ( STATES(10)*ALGBRC(4))/CONSTS(44)))/( (1.00000+STATES(9)/CONSTS(40)+STATES(10)/CONSTS(41))*(1.00000+ALGBRC(4)/CONSTS(43)+ALGBRC(3)/CONSTS(42)))
      RATES(9) = ALGBRC(8) - ALGBRC(13)
      ALGBRC(15) = ( CONSTS(58)*(STATES(10) ** CONSTS(56)/CONSTS(57) ** CONSTS(56)))/(1.00000+STATES(10) ** CONSTS(56)/CONSTS(57) ** CONSTS(56))
      RATES(10) = ALGBRC(13) - ALGBRC(15)
      ALGBRC(5) = CONSTS(6) - (ALGBRC(4)+ALGBRC(3))
      ALGBRC(19) =  CONSTS(72)*(1.00000+ CONSTS(79)*(ALGBRC(4)/CONSTS(78)))/(1.00000+ALGBRC(4)/CONSTS(78)) ** 2.00000*(1.00000+ CONSTS(77)*(ALGBRC(5)/CONSTS(76)))/(1.00000+ALGBRC(5)/CONSTS(76)) ** 2.00000*(1.00000+ CONSTS(81)*(CONSTS(4)/CONSTS(80))+ CONSTS(83)*(STATES(4)/CONSTS(82)))/(1.00000+CONSTS(4)/CONSTS(80)+STATES(4)/CONSTS(82)) ** 2.00000
      ALGBRC(18) = 1.00000+STATES(3)/CONSTS(73)+ALGBRC(4)/CONSTS(74)+ CONSTS(71)*(STATES(3)/CONSTS(73))*(ALGBRC(4)/CONSTS(74))
      ALGBRC(20) = 1.00000+ CONSTS(75)*(ALGBRC(4)/CONSTS(74))
      ALGBRC(21) =  CONSTS(70)*CONSTS(71)*(STATES(3)/CONSTS(73))*(ALGBRC(4)/CONSTS(74))*(ALGBRC(18)/(ALGBRC(18) ** 2.00000+ ALGBRC(19)*ALGBRC(20) ** 2.00000))
      RATES(3) = ALGBRC(6) - ALGBRC(21)
      ALGBRC(14) = ( (CONSTS(55)/CONSTS(50))*(STATES(4) - ( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5)*(1.00000/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(54)))/(1.00000+STATES(4)/CONSTS(50)+( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(51)+( (1.00000/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(53)+( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5)*(1.00000/(1.00000+CONSTS(5)))*STATES(5))/( CONSTS(51)*CONSTS(53))+( STATES(4)*(CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5))/( CONSTS(52)*CONSTS(50)))
      RATES(4) = ALGBRC(21) - ALGBRC(14)
      ALGBRC(22) = ( (CONSTS(89)/( CONSTS(84)*CONSTS(86)))*( (1.00000/(1.00000+CONSTS(5)))*STATES(5)*STATES(13) - ( CONSTS(9)*STATES(14))/CONSTS(88)))/( (1.00000+( (1.00000/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(84)+CONSTS(9)/CONSTS(85))*(1.00000+STATES(13)/CONSTS(86)+STATES(14)/CONSTS(87)))
      RATES(5) =  2.00000*ALGBRC(14) - (ALGBRC(11)+ALGBRC(22))
      ALGBRC(16) = - ( (CONSTS(67)/( CONSTS(64)*CONSTS(62)))*(( STATES(14)*CONSTS(8) - ( STATES(13)*STATES(11))/CONSTS(68))/(1.00000+STATES(14)/CONSTS(64)+( CONSTS(59)*CONSTS(8))/( CONSTS(64)*CONSTS(62))+( CONSTS(60)*STATES(11))/( CONSTS(63)*CONSTS(61))+STATES(13)/CONSTS(63)+( STATES(14)*CONSTS(8))/( CONSTS(64)*CONSTS(62))+( CONSTS(60)*STATES(14)*STATES(11))/( CONSTS(64)*CONSTS(63)*CONSTS(61))+( CONSTS(59)*CONSTS(8)*STATES(13))/( CONSTS(64)*CONSTS(62)*CONSTS(63))+( STATES(13)*STATES(11))/( CONSTS(63)*CONSTS(61))+( STATES(14)*CONSTS(8)*STATES(11))/( CONSTS(64)*CONSTS(62)*CONSTS(65))+( CONSTS(8)*STATES(13)*STATES(11))/( CONSTS(66)*CONSTS(63)*CONSTS(61)))))
      ALGBRC(23) =  CONSTS(90)*STATES(11)
      RATES(11) = ALGBRC(15) - (ALGBRC(16)+ 2.00000*ALGBRC(23))
      ALGBRC(17) =  CONSTS(69)*ALGBRC(4)
      RATES(12) = (ALGBRC(12)+ALGBRC(13)) - (ALGBRC(9)+ALGBRC(21)+ALGBRC(17)+CONSTS(2)+CONSTS(3)+ 4.00000*ALGBRC(23))
      RATES(13) = (ALGBRC(11)+ 3.00000*ALGBRC(23)) - (ALGBRC(16)+ALGBRC(22))
      RATES(14) = (ALGBRC(16)+ALGBRC(22)) - (ALGBRC(11)+ 3.00000*ALGBRC(23))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(7) = ( (CONSTS(17)/CONSTS(14))*(STATES(7) - STATES(8)/CONSTS(16)))/(1.00000+STATES(7)/CONSTS(14)+STATES(8)/CONSTS(15))
      ALGBRC(8) = ( (CONSTS(21)/CONSTS(18))*(STATES(8) - STATES(9)/CONSTS(20)))/(1.00000+STATES(8)/CONSTS(18)+STATES(9)/CONSTS(19))
      ALGBRC(3) = (CONSTS(6) -  ( STATES(12) ** 2.00000*(1.00000 -  4.00000*CONSTS(7))+ 2.00000*CONSTS(6)*STATES(12)*( 4.00000*CONSTS(7) - 1.00000)+CONSTS(6) ** 2.00000) ** (1.0 / 2))/(1.00000 -  4.00000*CONSTS(7))
      ALGBRC(4) = (STATES(12) - ALGBRC(3))/2.00000
      ALGBRC(9) = ( (CONSTS(27)/( CONSTS(22)*CONSTS(24)))*( STATES(1)*ALGBRC(4) - ( STATES(2)*ALGBRC(3))/CONSTS(26)))/( (1.00000+STATES(1)/CONSTS(22)+STATES(2)/CONSTS(23))*(1.00000+ALGBRC(4)/CONSTS(24)+ALGBRC(3)/CONSTS(25)))
      ALGBRC(6) = ( (CONSTS(13)/CONSTS(10))*(STATES(2) - STATES(3)/CONSTS(12)))/(1.00000+STATES(2)/CONSTS(10)+STATES(3)/CONSTS(11))
      ALGBRC(10) = ( CONSTS(49)*(CONSTS(1) - STATES(1)/CONSTS(48)))/( CONSTS(46)*(1.00000+CONSTS(1)/CONSTS(46)+STATES(1)/CONSTS(47)+( 0.910000*CONSTS(1)*STATES(1))/( CONSTS(47)*CONSTS(46))))
      ALGBRC(11) = (( CONSTS(32)*(CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5)*STATES(14))/( CONSTS(28)*CONSTS(30)) - ( CONSTS(33)*STATES(6)*STATES(13))/( CONSTS(29)*CONSTS(31)))/( (1.00000+( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(28)+STATES(6)/CONSTS(29))*(1.00000+STATES(14)/CONSTS(30)+STATES(13)/CONSTS(31)))
      ALGBRC(12) = ( (CONSTS(39)/( CONSTS(35)*CONSTS(37)))*( CONSTS(38)*STATES(6)*ALGBRC(3) -  STATES(7)*ALGBRC(4)))/( (1.00000+STATES(6)/CONSTS(34)+STATES(7)/CONSTS(35))*(1.00000+ALGBRC(4)/CONSTS(37)+ALGBRC(3)/CONSTS(36)))
      ALGBRC(13) = ( (CONSTS(45)/( CONSTS(40)*CONSTS(42)))*( STATES(9)*ALGBRC(3) - ( STATES(10)*ALGBRC(4))/CONSTS(44)))/( (1.00000+STATES(9)/CONSTS(40)+STATES(10)/CONSTS(41))*(1.00000+ALGBRC(4)/CONSTS(43)+ALGBRC(3)/CONSTS(42)))
      ALGBRC(15) = ( CONSTS(58)*(STATES(10) ** CONSTS(56)/CONSTS(57) ** CONSTS(56)))/(1.00000+STATES(10) ** CONSTS(56)/CONSTS(57) ** CONSTS(56))
      ALGBRC(5) = CONSTS(6) - (ALGBRC(4)+ALGBRC(3))
      ALGBRC(19) =  CONSTS(72)*(1.00000+ CONSTS(79)*(ALGBRC(4)/CONSTS(78)))/(1.00000+ALGBRC(4)/CONSTS(78)) ** 2.00000*(1.00000+ CONSTS(77)*(ALGBRC(5)/CONSTS(76)))/(1.00000+ALGBRC(5)/CONSTS(76)) ** 2.00000*(1.00000+ CONSTS(81)*(CONSTS(4)/CONSTS(80))+ CONSTS(83)*(STATES(4)/CONSTS(82)))/(1.00000+CONSTS(4)/CONSTS(80)+STATES(4)/CONSTS(82)) ** 2.00000
      ALGBRC(18) = 1.00000+STATES(3)/CONSTS(73)+ALGBRC(4)/CONSTS(74)+ CONSTS(71)*(STATES(3)/CONSTS(73))*(ALGBRC(4)/CONSTS(74))
      ALGBRC(20) = 1.00000+ CONSTS(75)*(ALGBRC(4)/CONSTS(74))
      ALGBRC(21) =  CONSTS(70)*CONSTS(71)*(STATES(3)/CONSTS(73))*(ALGBRC(4)/CONSTS(74))*(ALGBRC(18)/(ALGBRC(18) ** 2.00000+ ALGBRC(19)*ALGBRC(20) ** 2.00000))
      ALGBRC(14) = ( (CONSTS(55)/CONSTS(50))*(STATES(4) - ( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5)*(1.00000/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(54)))/(1.00000+STATES(4)/CONSTS(50)+( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(51)+( (1.00000/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(53)+( (CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5)*(1.00000/(1.00000+CONSTS(5)))*STATES(5))/( CONSTS(51)*CONSTS(53))+( STATES(4)*(CONSTS(5)/(1.00000+CONSTS(5)))*STATES(5))/( CONSTS(52)*CONSTS(50)))
      ALGBRC(22) = ( (CONSTS(89)/( CONSTS(84)*CONSTS(86)))*( (1.00000/(1.00000+CONSTS(5)))*STATES(5)*STATES(13) - ( CONSTS(9)*STATES(14))/CONSTS(88)))/( (1.00000+( (1.00000/(1.00000+CONSTS(5)))*STATES(5))/CONSTS(84)+CONSTS(9)/CONSTS(85))*(1.00000+STATES(13)/CONSTS(86)+STATES(14)/CONSTS(87)))
      ALGBRC(16) = - ( (CONSTS(67)/( CONSTS(64)*CONSTS(62)))*(( STATES(14)*CONSTS(8) - ( STATES(13)*STATES(11))/CONSTS(68))/(1.00000+STATES(14)/CONSTS(64)+( CONSTS(59)*CONSTS(8))/( CONSTS(64)*CONSTS(62))+( CONSTS(60)*STATES(11))/( CONSTS(63)*CONSTS(61))+STATES(13)/CONSTS(63)+( STATES(14)*CONSTS(8))/( CONSTS(64)*CONSTS(62))+( CONSTS(60)*STATES(14)*STATES(11))/( CONSTS(64)*CONSTS(63)*CONSTS(61))+( CONSTS(59)*CONSTS(8)*STATES(13))/( CONSTS(64)*CONSTS(62)*CONSTS(63))+( STATES(13)*STATES(11))/( CONSTS(63)*CONSTS(61))+( STATES(14)*CONSTS(8)*STATES(11))/( CONSTS(64)*CONSTS(62)*CONSTS(65))+( CONSTS(8)*STATES(13)*STATES(11))/( CONSTS(66)*CONSTS(63)*CONSTS(61)))))
      ALGBRC(23) =  CONSTS(90)*STATES(11)
      ALGBRC(17) =  CONSTS(69)*ALGBRC(4)
      ALGBRC(1) = STATES(5)/(1.00000+CONSTS(5))
      ALGBRC(2) = ( STATES(5)*CONSTS(5))/(1.00000+CONSTS(5))
      RETURN
      END