Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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