C
C There are a total of 53 entries in the algebraic variable array.
C There are a total of 22 entries in each of the rate and state variable arrays.
C There are a total of 83 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C CONSTS(1) is eps in component E (dimensionless).
C CONSTS(2) is K1 in component A (per_hour).
C CONSTS(3) is K1a in component A (per_hour).
C CONSTS(4) is K2 in component A (per_hour).
C CONSTS(5) is K2a in component A (per_hour).
C CONSTS(6) is K2aa in component A (per_hour).
C CONSTS(7) is K3 in component A (per_hour).
C CONSTS(8) is K3a in component A (per_hour).
C CONSTS(9) is K4 in component A (per_hour).
C CONSTS(10) is J1 in component A (dimensionless).
C CONSTS(11) is J3 in component A (dimensionless).
C CONSTS(12) is J4 in component A (dimensionless).
C ALGBRC(2) is v_20 in component A (per_hour).
C ALGBRC(22) is v_21 in component A (per_hour).
C ALGBRC(1) is v_42 in component A (per_hour).
C ALGBRC(23) is v_19 in component A (per_hour).
C CONSTS(13) is GA in component A (dimensionless).
C CONSTS(14) is GB in component A (dimensionless).
C CONSTS(15) is GE in component A (dimensionless).
C ALGBRC(13) is V_2 in component A (per_hour).
C ALGBRC(14) is V_4 in component A (per_hour).
C STATES(1) is Cdh1 in component A (dimensionless).
C CONSTS(16) is Cdc20 in component D (dimensionless).
C STATES(2) is CYCA in component D (dimensionless).
C STATES(3) is CYCB in component A (dimensionless).
C STATES(4) is CYCE in component D (dimensionless).
C CONSTS(17) is k15 in component B (per_hour).
C CONSTS(18) is k16 in component B (per_hour).
C CONSTS(19) is k17 in component B (per_hour).
C CONSTS(20) is k17a in component B (per_hour).
C CONSTS(21) is k18 in component B (per_hour).
C CONSTS(22) is K9 in component B (per_hour).
C CONSTS(23) is J17 in component B (dimensionless).
C CONSTS(24) is J15 in component B (dimensionless).
C ALGBRC(15) is v_1 in component B (per_hour).
C ALGBRC(16) is v_2 in component B (per_hour).
C ALGBRC(3) is v_34 in component B (per_hour).
C ALGBRC(6) is v_39 in component B (per_hour).
C ALGBRC(4) is v_41 in component B (per_hour).
C STATES(5) is DRG in component B (dimensionless).
C STATES(6) is ERG in component B (dimensionless).
C CONSTS(25) is K19 in component C (per_hour).
C CONSTS(26) is K19a in component C (per_hour).
C CONSTS(27) is K20 in component C (per_hour).
C CONSTS(28) is K21 in component C (dimensionless).
C CONSTS(29) is K22 in component C (per_hour).
C CONSTS(30) is K23a in component C (per_hour).
C CONSTS(31) is K23 in component C (per_hour).
C CONSTS(32) is K26 in component C (per_hour).
C CONSTS(33) is K26R in component C (per_hour).
C ALGBRC(17) is v_29 in component C (per_hour).
C ALGBRC(24) is v_30 in component C (per_hour).
C ALGBRC(26) is v_43 in component C (per_hour).
C ALGBRC(28) is v_44 in component C (per_hour).
C ALGBRC(30) is v_45 in component C (per_hour).
C ALGBRC(32) is v_46 in component C (per_hour).
C ALGBRC(34) is v_47 in component C (per_hour).
C ALGBRC(36) is v_48 in component C (per_hour).
C ALGBRC(38) is v_49 in component C (per_hour).
C ALGBRC(41) is v_50 in component C (per_hour).
C ALGBRC(43) is v_51 in component C (per_hour).
C ALGBRC(45) is v_52 in component C (per_hour).
C CONSTS(34) is LA in component C (dimensionless).
C CONSTS(35) is LB in component C (dimensionless).
C CONSTS(36) is LD in component C (dimensionless).
C CONSTS(37) is LE in component C (dimensionless).
C CONSTS(38) is FE in component C (dimensionless).
C CONSTS(39) is FB in component C (dimensionless).
C STATES(7) is Rb in component C (dimensionless).
C STATES(8) is PPRb in component C (dimensionless).
C STATES(9) is E2F in component C (dimensionless).
C STATES(10) is PE2F in component C (dimensionless).
C STATES(11) is E2FRb in component C (dimensionless).
C STATES(12) is PE2FRb in component C (dimensionless).
C ALGBRC(5) is PP1A in component C (dimensionless).
C CONSTS(40) is PP1T in component C (dimensionless).
C CONSTS(41) is CYCDT in component C (dimensionless).
C CONSTS(42) is K5 in component D (per_hour).
C CONSTS(43) is K6a in component D (per_hour).
C CONSTS(44) is K6 in component D (per_hour).
C CONSTS(45) is K7a in component D (per_hour).
C CONSTS(46) is K7 in component D (per_hour).
C CONSTS(47) is K8a in component D (per_hour).
C CONSTS(48) is K8 in component D (per_hour).
C CONSTS(49) is K10 in component D (per_hour).
C CONSTS(50) is k24 in component D (per_hour).
C CONSTS(51) is k24r in component D (per_hour).
C CONSTS(52) is K25 in component D (per_hour).
C CONSTS(53) is K25R in component D (per_hour).
C CONSTS(54) is K29 in component D (per_hour).
C CONSTS(55) is K30 in component D (per_hour).
C CONSTS(56) is J8 in component D (dimensionless).
C ALGBRC(27) is v_3 in component D (per_hour).
C ALGBRC(29) is v_4 in component D (per_hour).
C ALGBRC(31) is v_5 in component D (per_hour).
C ALGBRC(33) is v_6 in component D (per_hour).
C ALGBRC(35) is v_7 in component D (per_hour).
C ALGBRC(37) is v_8 in component D (per_hour).
C ALGBRC(39) is v_9 in component D (per_hour).
C ALGBRC(42) is v_10 in component D (per_hour).
C ALGBRC(44) is v_11 in component D (per_hour).
C ALGBRC(46) is v_12 in component D (per_hour).
C ALGBRC(47) is v_13 in component D (per_hour).
C ALGBRC(50) is v_14 in component D (per_hour).
C ALGBRC(49) is v_15 in component D (per_hour).
C ALGBRC(52) is v_16 in component D (per_hour).
C ALGBRC(40) is v_17 in component D (per_hour).
C ALGBRC(48) is v_18 in component D (per_hour).
C ALGBRC(51) is v_36 in component D (per_hour).
C ALGBRC(53) is v_38 in component D (per_hour).
C CONSTS(80) is v_40 in component D (per_hour).
C CONSTS(57) is HA in component D (dimensionless).
C CONSTS(58) is HB in component D (dimensionless).
C CONSTS(59) is HE in component D (dimensionless).
C CONSTS(60) is YE in component D (dimensionless).
C CONSTS(61) is YB in component D (dimensionless).
C STATES(13) is p27 in component D (dimensionless).
C ALGBRC(18) is V_6 in component D (per_hour).
C ALGBRC(25) is V_8 in component D (per_hour).
C STATES(14) is MASS in component E (dimensionless).
C STATES(15) is CA in component D (dimensionless).
C STATES(16) is CD in component D (dimensionless).
C STATES(17) is CE in component D (dimensionless).
C STATES(18) is CYCD in component D (dimensionless).
C CONSTS(62) is CYCET in component D (dimensionless).
C CONSTS(63) is K27 in component E (per_hour).
C CONSTS(64) is K28 in component E (per_hour).
C ALGBRC(7) is v_31 in component E (per_hour).
C ALGBRC(19) is v_32 in component E (per_hour).
C ALGBRC(8) is v_33 in component E (per_hour).
C CONSTS(65) is r31switch in component E (dimensionless).
C CONSTS(66) is MU in component E (per_hour).
C STATES(19) is GM in component E (dimensionless).
C CONSTS(67) is K11a in component F (per_hour).
C CONSTS(68) is K11 in component F (per_hour).
C CONSTS(69) is K12 in component F (per_hour).
C CONSTS(70) is K13 in component F (per_hour).
C CONSTS(71) is K14 in component F (per_hour).
C CONSTS(72) is K31 in component F (per_hour).
C CONSTS(73) is K32 in component F (per_hour).
C CONSTS(74) is K33 in component F (per_hour).
C CONSTS(75) is K34 in component F (per_hour).
C CONSTS(76) is J13 in component F (dimensionless).
C CONSTS(77) is J14 in component F (dimensionless).
C CONSTS(78) is J31 in component F (dimensionless).
C CONSTS(79) is J32 in component F (dimensionless).
C ALGBRC(9) is v_22 in component F (per_hour).
C ALGBRC(10) is v_23 in component F (per_hour).
C ALGBRC(20) is v_24 in component F (per_hour).
C ALGBRC(12) is v_25 in component F (per_hour).
C ALGBRC(11) is v_26 in component F (per_hour).
C CONSTS(81) is v_27 in component F (per_hour).
C CONSTS(82) is v_28 in component F (per_hour).
C ALGBRC(21) is v_35 in component F (per_hour).
C CONSTS(83) is v_37 in component F (per_hour).
C STATES(20) is IEP in component F (dimensionless).
C STATES(21) is PPX in component F (dimensionless).
C STATES(22) is Cdc20T in component F (dimensionless).
C RATES(1) is d/dt Cdh1 in component A (dimensionless).
C RATES(3) is d/dt CYCB in component A (dimensionless).
C RATES(5) is d/dt DRG in component B (dimensionless).
C RATES(6) is d/dt ERG in component B (dimensionless).
C RATES(8) is d/dt PPRb in component C (dimensionless).
C RATES(9) is d/dt E2F in component C (dimensionless).
C RATES(10) is d/dt PE2F in component C (dimensionless).
C RATES(7) is d/dt Rb in component C (dimensionless).
C RATES(11) is d/dt E2FRb in component C (dimensionless).
C RATES(12) is d/dt PE2FRb in component C (dimensionless).
C RATES(15) is d/dt CA in component D (dimensionless).
C RATES(16) is d/dt CD in component D (dimensionless).
C RATES(17) is d/dt CE in component D (dimensionless).
C RATES(2) is d/dt CYCA in component D (dimensionless).
C RATES(18) is d/dt CYCD in component D (dimensionless).
C RATES(4) is d/dt CYCE in component D (dimensionless).
C RATES(13) is d/dt p27 in component D (dimensionless).
C RATES(19) is d/dt GM in component E (dimensionless).
C RATES(14) is d/dt MASS in component E (dimensionless).
C RATES(22) is d/dt Cdc20T in component F (dimensionless).
C RATES(20) is d/dt IEP in component F (dimensionless).
C RATES(21) is d/dt PPX in component F (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      CONSTS(2) = 0.6
      CONSTS(3) = 0.1
      CONSTS(4) = 20
      CONSTS(5) = 0.05
      CONSTS(6) = 1
      CONSTS(7) = 140
      CONSTS(8) = 7.5
      CONSTS(9) = 40
      CONSTS(10) = 0.1
      CONSTS(11) = 0.01
      CONSTS(12) = 0.04
      CONSTS(13) = 0.3
      CONSTS(14) = 1
      CONSTS(15) = 0
      STATES(1) = 0.000653278
      CONSTS(16) = 0.00220177
      STATES(2) = 1.4094
      STATES(3) = 2.72898
      STATES(4) = 0.0229112
      CONSTS(17) = 0.25
      CONSTS(18) = 0.25
      CONSTS(19) = 10
      CONSTS(20) = 0.35
      CONSTS(21) = 10
      CONSTS(22) = 2.5
      CONSTS(23) = 0.3
      CONSTS(24) = 0.1
      STATES(5) = 0.900533
      STATES(6) = 0.0121809
      CONSTS(25) = 20
      CONSTS(26) = 0
      CONSTS(27) = 10
      CONSTS(28) = 1
      CONSTS(29) = 1
      CONSTS(30) = 0.005
      CONSTS(31) = 1
      CONSTS(32) = 10000
      CONSTS(33) = 200
      CONSTS(34) = 3
      CONSTS(35) = 5
      CONSTS(36) = 3.3
      CONSTS(37) = 5
      CONSTS(38) = 25
      CONSTS(39) = 2
      STATES(7) = 0.000190871
      STATES(8) = 9.97574
      STATES(9) = 0.989986
      STATES(10) = 3.98594
      STATES(11) = 0.00478911
      STATES(12) = 0.0192822
      CONSTS(40) = 1
      CONSTS(41) = 0.010976
      CONSTS(42) = 20
      CONSTS(43) = 10
      CONSTS(44) = 100
      CONSTS(45) = 0
      CONSTS(46) = 0.6
      CONSTS(47) = 0.1
      CONSTS(48) = 2
      CONSTS(49) = 5
      CONSTS(50) = 1000
      CONSTS(51) = 10
      CONSTS(52) = 1000
      CONSTS(53) = 10
      CONSTS(54) = 0.05
      CONSTS(55) = 20
      CONSTS(56) = 0.1
      CONSTS(57) = 0.5
      CONSTS(58) = 1
      CONSTS(59) = 0.5
      CONSTS(60) = 1
      CONSTS(61) = 0.05
      STATES(13) = 0.00922806
      STATES(14) = 1.68776
      STATES(15) = 0.0356927
      STATES(16) = 0.010976
      STATES(17) = 0.000542587
      STATES(18) = 0.43929
      CONSTS(62) = 0.000542587
      CONSTS(63) = 0.2
      CONSTS(64) = 0.2
      CONSTS(65) = 1
      CONSTS(66) = 0.061
      STATES(19) = 1.35565
      CONSTS(67) = 0
      CONSTS(68) = 1.5
      CONSTS(69) = 1.5
      CONSTS(70) = 5
      CONSTS(71) = 2.5
      CONSTS(72) = 0.7
      CONSTS(73) = 1.8
      CONSTS(74) = 0.05
      CONSTS(75) = 0.05
      CONSTS(76) = 0.005
      CONSTS(77) = 0.005
      CONSTS(78) = 0.01
      CONSTS(79) = 0.01
      STATES(20) = 0.154655
      STATES(21) = 1
      STATES(22) = 2.36733
      CONSTS(80) =  CONSTS(1)*CONSTS(42)
      CONSTS(81) = ( CONSTS(71)*CONSTS(16))/(CONSTS(77)+CONSTS(16))
      CONSTS(82) =  CONSTS(69)*CONSTS(16)
      CONSTS(83) =  CONSTS(1)*CONSTS(74)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(8) =  CONSTS(1)*CONSTS(66)*STATES(19)
      RATES(14) = ALGBRC(8)
      ALGBRC(9) =  CONSTS(75)*STATES(21)
      RATES(21) = CONSTS(83) - ALGBRC(9)
      ALGBRC(16) =  CONSTS(21)*STATES(5)
      ALGBRC(4) =  CONSTS(1)*(( CONSTS(19)*STATES(5)/CONSTS(23) ** 2.00000)/(STATES(5)/CONSTS(23) ** 2.00000+1.00000)+ CONSTS(20)*STATES(6))
      RATES(5) = ALGBRC(4) - ALGBRC(16)
      ALGBRC(15) =  CONSTS(18)*STATES(6)
      ALGBRC(3) = ( CONSTS(1)*CONSTS(17))/(STATES(5)/CONSTS(24) ** 2.00000+1.00000)
      RATES(6) = ALGBRC(3) - ALGBRC(15)
      ALGBRC(7) =  CONSTS(63)*STATES(14)*CONSTS(65)
      ALGBRC(19) =  CONSTS(64)*STATES(19)
      RATES(19) = ALGBRC(7) - ALGBRC(19)
      ALGBRC(12) =  CONSTS(69)*STATES(22)
      ALGBRC(21) =  CONSTS(1)*CONSTS(67)+ CONSTS(68)*STATES(3)
      RATES(22) = ALGBRC(21) - ALGBRC(12)
      ALGBRC(10) = ( CONSTS(72)*STATES(3)*(1.00000 - STATES(20)))/((CONSTS(78) - STATES(20))+1.00000)
      ALGBRC(20) = ( CONSTS(73)*STATES(21)*STATES(20))/(CONSTS(79)+STATES(20))
      RATES(20) = ALGBRC(10) - ALGBRC(20)
      ALGBRC(2) = ( (CONSTS(8)+ CONSTS(7)*CONSTS(16))*(1.00000 - STATES(1)))/((CONSTS(11) - STATES(1))+1.00000)
      ALGBRC(14) =  CONSTS(9)*( CONSTS(13)*STATES(2)+ CONSTS(14)*STATES(3)+ CONSTS(15)*STATES(4))
      ALGBRC(22) = ( ALGBRC(14)*STATES(1))/(CONSTS(12)+STATES(1))
      RATES(1) = ALGBRC(2) - ALGBRC(22)
      ALGBRC(1) =  CONSTS(1)*(( CONSTS(2)*STATES(3)/CONSTS(10) ** 2.00000)/(STATES(3)/CONSTS(10) ** 2.00000+1.00000)+CONSTS(3))
      ALGBRC(13) =  CONSTS(6)*CONSTS(16)+ CONSTS(5)*(1.00000 - STATES(1))+ CONSTS(4)*STATES(1)
      ALGBRC(23) =  ALGBRC(13)*STATES(3)
      RATES(3) = ALGBRC(1) - ALGBRC(23)
      ALGBRC(17) =  STATES(11)*CONSTS(27)*( CONSTS(41)*CONSTS(36)+ CONSTS(34)*STATES(2)+ CONSTS(35)*STATES(3)+ CONSTS(37)*STATES(4))
      ALGBRC(24) =  STATES(12)*CONSTS(27)*( CONSTS(41)*CONSTS(36)+ CONSTS(34)*STATES(2)+ CONSTS(35)*STATES(3)+ CONSTS(37)*STATES(4))
      ALGBRC(26) =  STATES(7)*CONSTS(27)*( CONSTS(41)*CONSTS(36)+ CONSTS(34)*STATES(2)+ CONSTS(35)*STATES(3)+ CONSTS(37)*STATES(4))
      ALGBRC(5) = CONSTS(40)/( CONSTS(28)*CONSTS(38)*(STATES(2)+STATES(4)+ CONSTS(39)*STATES(3))+1.00000)
      ALGBRC(28) =  STATES(8)*( CONSTS(26)*(CONSTS(40) - ALGBRC(5))+ CONSTS(25)*ALGBRC(5))
      RATES(8) = (ALGBRC(17)+ALGBRC(24)+ALGBRC(26)) - ALGBRC(28)
      ALGBRC(30) =  STATES(11)*CONSTS(33)
      ALGBRC(32) =  STATES(9)*(CONSTS(30)+ CONSTS(31)*(STATES(2)+STATES(3)))
      ALGBRC(34) =  STATES(10)*CONSTS(29)
      ALGBRC(36) =  STATES(9)*STATES(7)*CONSTS(32)
      RATES(9) = ((ALGBRC(17)+ALGBRC(30)+ALGBRC(34)) - ALGBRC(32)) - ALGBRC(36)
      ALGBRC(27) =  CONSTS(49)*STATES(16)
      ALGBRC(35) =  CONSTS(50)*STATES(13)*STATES(18)
      ALGBRC(37) =  CONSTS(51)*STATES(16)
      ALGBRC(18) = CONSTS(43)+ CONSTS(44)*( CONSTS(57)*STATES(2)+ CONSTS(58)*STATES(3)+ CONSTS(59)*STATES(4))
      ALGBRC(40) =  ALGBRC(18)*STATES(16)
      RATES(16) = ((ALGBRC(35) - ALGBRC(37)) - ALGBRC(40)) - ALGBRC(27)
      ALGBRC(6) =  CONSTS(1)*CONSTS(22)*STATES(5)
      ALGBRC(29) =  CONSTS(49)*STATES(18)
      RATES(18) = ((ALGBRC(6)+ALGBRC(40)+ALGBRC(37)) - ALGBRC(35)) - ALGBRC(29)
      ALGBRC(38) =  STATES(12)*CONSTS(33)
      ALGBRC(41) =  STATES(7)*STATES(10)*CONSTS(32)
      RATES(10) = ((ALGBRC(24)+ALGBRC(38)+ALGBRC(32)) - ALGBRC(34)) - ALGBRC(41)
      RATES(7) = (((ALGBRC(28)+ALGBRC(30)+ALGBRC(38)) - ALGBRC(36)) - ALGBRC(41)) - ALGBRC(26)
      ALGBRC(43) =  STATES(12)*CONSTS(29)
      ALGBRC(45) =  STATES(11)*(CONSTS(30)+ CONSTS(31)*(STATES(2)+STATES(3)))
      RATES(11) = (((ALGBRC(43)+ALGBRC(36)) - ALGBRC(45)) - ALGBRC(17)) - ALGBRC(30)
      RATES(12) = (((ALGBRC(45)+ALGBRC(41)) - ALGBRC(43)) - ALGBRC(24)) - ALGBRC(38)
      ALGBRC(33) =  CONSTS(52)*STATES(13)*STATES(2)
      ALGBRC(42) =  CONSTS(55)*CONSTS(16)*STATES(15)
      ALGBRC(46) =  CONSTS(53)*STATES(15)
      ALGBRC(48) =  ALGBRC(18)*STATES(15)
      RATES(15) = ((ALGBRC(33) - ALGBRC(46)) - ALGBRC(48)) - ALGBRC(42)
      ALGBRC(39) =  CONSTS(55)*CONSTS(16)*STATES(2)
      ALGBRC(51) =  CONSTS(1)*CONSTS(54)*STATES(9)*STATES(14)
      RATES(2) = ((ALGBRC(51) - ALGBRC(39)) - ALGBRC(33))+ALGBRC(46)+ALGBRC(48)
      ALGBRC(31) =  CONSTS(52)*STATES(13)*STATES(4)
      ALGBRC(44) =  CONSTS(53)*STATES(17)
      ALGBRC(25) = ( CONSTS(48)*( CONSTS(60)*(STATES(2)+STATES(4))+ CONSTS(61)*STATES(3)))/(CONSTS(62)+CONSTS(56))+CONSTS(47)
      ALGBRC(47) =  ALGBRC(25)*STATES(17)
      ALGBRC(49) =  ALGBRC(18)*STATES(13)
      RATES(13) = (((((CONSTS(80)+ALGBRC(27)+ALGBRC(37)) - ALGBRC(49)) - ALGBRC(31)) - ALGBRC(33)) - ALGBRC(35))+ALGBRC(44)+ALGBRC(46)+ALGBRC(47)+ALGBRC(42)
      ALGBRC(52) =  ALGBRC(18)*STATES(17)
      RATES(17) = ((ALGBRC(31) - ALGBRC(44)) - ALGBRC(47)) - ALGBRC(52)
      ALGBRC(50) =  ALGBRC(25)*STATES(4)
      ALGBRC(53) =  CONSTS(1)*(CONSTS(45)+ CONSTS(46)*STATES(9))
      RATES(4) = ((ALGBRC(53) - ALGBRC(50)) - ALGBRC(31))+ALGBRC(44)+ALGBRC(52)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(8) =  CONSTS(1)*CONSTS(66)*STATES(19)
      ALGBRC(9) =  CONSTS(75)*STATES(21)
      ALGBRC(16) =  CONSTS(21)*STATES(5)
      ALGBRC(4) =  CONSTS(1)*(( CONSTS(19)*STATES(5)/CONSTS(23) ** 2.00000)/(STATES(5)/CONSTS(23) ** 2.00000+1.00000)+ CONSTS(20)*STATES(6))
      ALGBRC(15) =  CONSTS(18)*STATES(6)
      ALGBRC(3) = ( CONSTS(1)*CONSTS(17))/(STATES(5)/CONSTS(24) ** 2.00000+1.00000)
      ALGBRC(7) =  CONSTS(63)*STATES(14)*CONSTS(65)
      ALGBRC(19) =  CONSTS(64)*STATES(19)
      ALGBRC(12) =  CONSTS(69)*STATES(22)
      ALGBRC(21) =  CONSTS(1)*CONSTS(67)+ CONSTS(68)*STATES(3)
      ALGBRC(10) = ( CONSTS(72)*STATES(3)*(1.00000 - STATES(20)))/((CONSTS(78) - STATES(20))+1.00000)
      ALGBRC(20) = ( CONSTS(73)*STATES(21)*STATES(20))/(CONSTS(79)+STATES(20))
      ALGBRC(2) = ( (CONSTS(8)+ CONSTS(7)*CONSTS(16))*(1.00000 - STATES(1)))/((CONSTS(11) - STATES(1))+1.00000)
      ALGBRC(14) =  CONSTS(9)*( CONSTS(13)*STATES(2)+ CONSTS(14)*STATES(3)+ CONSTS(15)*STATES(4))
      ALGBRC(22) = ( ALGBRC(14)*STATES(1))/(CONSTS(12)+STATES(1))
      ALGBRC(1) =  CONSTS(1)*(( CONSTS(2)*STATES(3)/CONSTS(10) ** 2.00000)/(STATES(3)/CONSTS(10) ** 2.00000+1.00000)+CONSTS(3))
      ALGBRC(13) =  CONSTS(6)*CONSTS(16)+ CONSTS(5)*(1.00000 - STATES(1))+ CONSTS(4)*STATES(1)
      ALGBRC(23) =  ALGBRC(13)*STATES(3)
      ALGBRC(17) =  STATES(11)*CONSTS(27)*( CONSTS(41)*CONSTS(36)+ CONSTS(34)*STATES(2)+ CONSTS(35)*STATES(3)+ CONSTS(37)*STATES(4))
      ALGBRC(24) =  STATES(12)*CONSTS(27)*( CONSTS(41)*CONSTS(36)+ CONSTS(34)*STATES(2)+ CONSTS(35)*STATES(3)+ CONSTS(37)*STATES(4))
      ALGBRC(26) =  STATES(7)*CONSTS(27)*( CONSTS(41)*CONSTS(36)+ CONSTS(34)*STATES(2)+ CONSTS(35)*STATES(3)+ CONSTS(37)*STATES(4))
      ALGBRC(5) = CONSTS(40)/( CONSTS(28)*CONSTS(38)*(STATES(2)+STATES(4)+ CONSTS(39)*STATES(3))+1.00000)
      ALGBRC(28) =  STATES(8)*( CONSTS(26)*(CONSTS(40) - ALGBRC(5))+ CONSTS(25)*ALGBRC(5))
      ALGBRC(30) =  STATES(11)*CONSTS(33)
      ALGBRC(32) =  STATES(9)*(CONSTS(30)+ CONSTS(31)*(STATES(2)+STATES(3)))
      ALGBRC(34) =  STATES(10)*CONSTS(29)
      ALGBRC(36) =  STATES(9)*STATES(7)*CONSTS(32)
      ALGBRC(27) =  CONSTS(49)*STATES(16)
      ALGBRC(35) =  CONSTS(50)*STATES(13)*STATES(18)
      ALGBRC(37) =  CONSTS(51)*STATES(16)
      ALGBRC(18) = CONSTS(43)+ CONSTS(44)*( CONSTS(57)*STATES(2)+ CONSTS(58)*STATES(3)+ CONSTS(59)*STATES(4))
      ALGBRC(40) =  ALGBRC(18)*STATES(16)
      ALGBRC(6) =  CONSTS(1)*CONSTS(22)*STATES(5)
      ALGBRC(29) =  CONSTS(49)*STATES(18)
      ALGBRC(38) =  STATES(12)*CONSTS(33)
      ALGBRC(41) =  STATES(7)*STATES(10)*CONSTS(32)
      ALGBRC(43) =  STATES(12)*CONSTS(29)
      ALGBRC(45) =  STATES(11)*(CONSTS(30)+ CONSTS(31)*(STATES(2)+STATES(3)))
      ALGBRC(33) =  CONSTS(52)*STATES(13)*STATES(2)
      ALGBRC(42) =  CONSTS(55)*CONSTS(16)*STATES(15)
      ALGBRC(46) =  CONSTS(53)*STATES(15)
      ALGBRC(48) =  ALGBRC(18)*STATES(15)
      ALGBRC(39) =  CONSTS(55)*CONSTS(16)*STATES(2)
      ALGBRC(51) =  CONSTS(1)*CONSTS(54)*STATES(9)*STATES(14)
      ALGBRC(31) =  CONSTS(52)*STATES(13)*STATES(4)
      ALGBRC(44) =  CONSTS(53)*STATES(17)
      ALGBRC(25) = ( CONSTS(48)*( CONSTS(60)*(STATES(2)+STATES(4))+ CONSTS(61)*STATES(3)))/(CONSTS(62)+CONSTS(56))+CONSTS(47)
      ALGBRC(47) =  ALGBRC(25)*STATES(17)
      ALGBRC(49) =  ALGBRC(18)*STATES(13)
      ALGBRC(52) =  ALGBRC(18)*STATES(17)
      ALGBRC(50) =  ALGBRC(25)*STATES(4)
      ALGBRC(53) =  CONSTS(1)*(CONSTS(45)+ CONSTS(46)*STATES(9))
      ALGBRC(11) = ( CONSTS(70)*STATES(20)*(STATES(22) - CONSTS(16)))/((CONSTS(76) - CONSTS(16))+STATES(22))
      RETURN
      END