C
C There are a total of 42 entries in the algebraic variable array.
C There are a total of 31 entries in each of the rate and state variable arrays.
C There are a total of 57 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is Casp8 in component Casp8 (micromolar).
C STATES(2) is Casp8Bid in component Casp8 (micromolar).
C STATES(3) is Bid in component Bid (micromolar).
C ALGBRC(3) is J_0 in component Casp8 (flux).
C ALGBRC(1) is J_f0 in component Casp8 (flux).
C ALGBRC(2) is J_Casp8 in component Casp8 (flux).
C CONSTS(1) is k_f0 in component model_constant (first_order_rate_constant).
C CONSTS(2) is k_10 in component model_constant (second_order_rate_constant).
C CONSTS(3) is k_00 in component model_constant (first_order_rate_constant).
C CONSTS(4) is u in component model_constant (first_order_rate_constant).
C STATES(4) is Apaf_1 in component Apaf_1 (micromolar).
C STATES(5) is CytcApaf_1 in component CytcApaf_1 (micromolar).
C STATES(6) is Cytc in component Cytc (micromolar).
C ALGBRC(4) is J_Apaf_1 in component Apaf_1 (flux).
C ALGBRC(5) is J_1 in component Apaf_1 (flux).
C CONSTS(5) is P_Apaf_1 in component model_constant (flux).
C CONSTS(6) is k_11 in component model_constant (second_order_rate_constant).
C CONSTS(7) is k_01 in component model_constant (first_order_rate_constant).
C STATES(7) is Apop in component Apop (micromolar).
C ALGBRC(7) is J_1b in component CytcApaf_1 (flux).
C CONSTS(8) is k_11b in component model_constant (rate).
C CONSTS(9) is k_01b in component model_constant (first_order_rate_constant).
C CONSTS(10) is p in component model_constant (dimensionless).
C ALGBRC(10) is J_14 in component Cytc_mito (flux).
C ALGBRC(6) is J_Cytc in component Cytc (flux).
C STATES(8) is Cytc_mito in component Cytc_mito (micromolar).
C STATES(9) is Bax_2 in component Bax_2 (micromolar).
C ALGBRC(8) is J_Cytc_mito in component Cytc_mito (flux).
C CONSTS(11) is P_Cytc_mito in component model_constant (flux).
C CONSTS(12) is k14 in component model_constant (second_order_rate_constant).
C ALGBRC(16) is J_12b in component tBidBax (flux).
C ALGBRC(9) is J_Bax_2 in component Bax_2 (flux).
C STATES(10) is tBid_mito in component tBid_mito (micromolar).
C STATES(11) is Bax in component Bax (micromolar).
C ALGBRC(13) is J_11 in component tBid (flux).
C ALGBRC(12) is J_12a in component tBid_mito (flux).
C ALGBRC(11) is J_tBid_mito in component tBid_mito (flux).
C CONSTS(13) is k12a in component model_constant (second_order_rate_constant).
C STATES(12) is tBid in component tBid (micromolar).
C ALGBRC(24) is J_f8 in component Casp3Bid (flux).
C ALGBRC(14) is J_tBid in component tBid (flux).
C CONSTS(14) is k11 in component model_constant (first_order_rate_constant).
C STATES(13) is tBidBax in component tBidBax (micromolar).
C ALGBRC(15) is J_tBidBax in component tBidBax (flux).
C CONSTS(15) is k12b in component model_constant (second_order_rate_constant).
C STATES(14) is Bcl_2 in component Bcl_2 (micromolar).
C ALGBRC(18) is J_13 in component Bax (flux).
C ALGBRC(17) is J_Bax in component Bax (flux).
C CONSTS(56) is P_Bax in component Bax (flux).
C CONSTS(16) is k13 in component model_constant (second_order_rate_constant).
C CONSTS(17) is P_oBax in component model_constant (flux).
C CONSTS(18) is p53 in component model_constant (micromolar).
C CONSTS(19) is p53_thresh in component model_constant (micromolar).
C CONSTS(20) is u_Bax in component model_constant (first_order_rate_constant).
C STATES(15) is Casp3 in component Casp3 (micromolar).
C STATES(16) is Casp3Bcl_2 in component Casp3Bcl_2 (micromolar).
C ALGBRC(19) is J_9 in component Bcl_2 (flux).
C ALGBRC(21) is J_Bcl_2 in component Bcl_2 (flux).
C CONSTS(21) is k_19 in component model_constant (second_order_rate_constant).
C CONSTS(22) is k_09 in component model_constant (first_order_rate_constant).
C CONSTS(57) is P_Bcl_2 in component Bcl_2 (flux).
C CONSTS(23) is P_oBcl_2 in component model_constant (flux).
C CONSTS(24) is u_Bcl_2 in component model_constant (first_order_rate_constant).
C ALGBRC(20) is J_f9 in component Casp3Bcl_2 (flux).
C CONSTS(25) is k_f9 in component model_constant (first_order_rate_constant).
C STATES(17) is Casp3Bid in component Casp3Bid (micromolar).
C ALGBRC(22) is J_8 in component Casp3Bid (flux).
C CONSTS(26) is k_18 in component model_constant (second_order_rate_constant).
C CONSTS(27) is k_08 in component model_constant (first_order_rate_constant).
C CONSTS(28) is k_f8 in component model_constant (first_order_rate_constant).
C ALGBRC(23) is J_Bid in component Bid (flux).
C CONSTS(29) is P_Bid in component model_constant (flux).
C STATES(18) is IAP in component IAP (micromolar).
C STATES(19) is Casp3IAP in component Casp3 (micromolar).
C ALGBRC(33) is J_f6 in component Casp9Pro3 (flux).
C ALGBRC(32) is J_f6b in component ApopCasp9_2Pro3 (flux).
C ALGBRC(26) is J_7 in component Casp3 (flux).
C ALGBRC(25) is J_Casp3 in component Casp3 (flux).
C CONSTS(30) is k_17 in component model_constant (second_order_rate_constant).
C CONSTS(31) is k_07 in component model_constant (first_order_rate_constant).
C ALGBRC(28) is J_2 in component Pro9 (flux).
C ALGBRC(41) is J_4b in component ApopCasp9 (flux).
C STATES(20) is Pro9 in component Pro9 (micromolar).
C STATES(21) is ApopPro9 in component Pro9 (micromolar).
C STATES(22) is ApopPro9_2 in component ApopPro9_2 (micromolar).
C ALGBRC(29) is J_3 in component Pro9 (flux).
C ALGBRC(27) is J_Pro9 in component Pro9 (flux).
C CONSTS(32) is k_12 in component model_constant (second_order_rate_constant).
C CONSTS(33) is k_02 in component model_constant (first_order_rate_constant).
C CONSTS(34) is k_13 in component model_constant (second_order_rate_constant).
C CONSTS(35) is k_03 in component model_constant (first_order_rate_constant).
C CONSTS(36) is P_Pro9 in component model_constant (flux).
C ALGBRC(30) is J_f3 in component ApopPro9_2 (flux).
C CONSTS(37) is k_f3 in component model_constant (first_order_rate_constant).
C STATES(23) is ApopCasp9_2Pro3 in component ApopCasp9_2Pro3 (micromolar).
C STATES(24) is ApopCasp9_2 in component ApopCasp9_2 (micromolar).
C STATES(25) is Pro3 in component Pro3 (micromolar).
C ALGBRC(31) is J_6b in component ApopCasp9_2Pro3 (flux).
C CONSTS(38) is k_16b in component model_constant (second_order_rate_constant).
C CONSTS(39) is k_06b in component model_constant (first_order_rate_constant).
C CONSTS(40) is k_f6b in component model_constant (first_order_rate_constant).
C STATES(26) is Casp9Pro3 in component Casp9Pro3 (micromolar).
C ALGBRC(34) is J_6 in component Pro3 (flux).
C CONSTS(41) is k_f6 in component model_constant (first_order_rate_constant).
C STATES(27) is Casp9 in component Casp9 (micromolar).
C ALGBRC(35) is J_Pro3 in component Pro3 (flux).
C CONSTS(42) is k_16 in component model_constant (second_order_rate_constant).
C CONSTS(43) is k_06 in component model_constant (first_order_rate_constant).
C CONSTS(44) is P_Pro3 in component model_constant (flux).
C STATES(28) is Casp9IAP in component IAP (micromolar).
C STATES(29) is ApopCasp9IAP in component IAP (micromolar).
C STATES(30) is ApopCasp9_2IAP in component IAP (micromolar).
C STATES(31) is ApopCasp9 in component ApopCasp9 (micromolar).
C ALGBRC(36) is J_5 in component IAP (flux).
C ALGBRC(37) is J_5b in component IAP (flux).
C ALGBRC(38) is J_5c in component IAP (flux).
C ALGBRC(40) is J_IAP in component IAP (flux).
C CONSTS(45) is P_IAP in component model_constant (flux).
C CONSTS(46) is k_15 in component model_constant (second_order_rate_constant).
C CONSTS(47) is k_05 in component model_constant (first_order_rate_constant).
C CONSTS(48) is k_15b in component model_constant (second_order_rate_constant).
C CONSTS(49) is k_05b in component model_constant (first_order_rate_constant).
C CONSTS(50) is k_15c in component model_constant (second_order_rate_constant).
C CONSTS(51) is k_05c in component model_constant (first_order_rate_constant).
C ALGBRC(39) is J_4 in component ApopCasp9 (flux).
C CONSTS(52) is k_14 in component model_constant (first_order_rate_constant).
C CONSTS(53) is k_04 in component model_constant (second_order_rate_constant).
C CONSTS(54) is k_14b in component model_constant (first_order_rate_constant).
C CONSTS(55) is k_04b in component model_constant (second_order_rate_constant).
C ALGBRC(42) is J_Casp9 in component Casp9 (flux).
C RATES(1) is d/dt Casp8 in component Casp8 (micromolar).
C RATES(2) is d/dt Casp8Bid in component Casp8 (micromolar).
C RATES(4) is d/dt Apaf_1 in component Apaf_1 (micromolar).
C RATES(5) is d/dt CytcApaf_1 in component CytcApaf_1 (micromolar).
C RATES(6) is d/dt Cytc in component Cytc (micromolar).
C RATES(8) is d/dt Cytc_mito in component Cytc_mito (micromolar).
C RATES(9) is d/dt Bax_2 in component Bax_2 (micromolar).
C RATES(10) is d/dt tBid_mito in component tBid_mito (micromolar).
C RATES(12) is d/dt tBid in component tBid (micromolar).
C RATES(13) is d/dt tBidBax in component tBidBax (micromolar).
C RATES(11) is d/dt Bax in component Bax (micromolar).
C RATES(14) is d/dt Bcl_2 in component Bcl_2 (micromolar).
C RATES(16) is d/dt Casp3Bcl_2 in component Casp3Bcl_2 (micromolar).
C RATES(17) is d/dt Casp3Bid in component Casp3Bid (micromolar).
C RATES(3) is d/dt Bid in component Bid (micromolar).
C RATES(19) is d/dt Casp3IAP in component Casp3 (micromolar).
C RATES(15) is d/dt Casp3 in component Casp3 (micromolar).
C RATES(7) is d/dt Apop in component Apop (micromolar).
C RATES(20) is d/dt Pro9 in component Pro9 (micromolar).
C RATES(21) is d/dt ApopPro9 in component Pro9 (micromolar).
C RATES(22) is d/dt ApopPro9_2 in component ApopPro9_2 (micromolar).
C RATES(23) is d/dt ApopCasp9_2Pro3 in component ApopCasp9_2Pro3 (micromolar).
C RATES(26) is d/dt Casp9Pro3 in component Casp9Pro3 (micromolar).
C RATES(25) is d/dt Pro3 in component Pro3 (micromolar).
C RATES(28) is d/dt Casp9IAP in component IAP (micromolar).
C RATES(29) is d/dt ApopCasp9IAP in component IAP (micromolar).
C RATES(30) is d/dt ApopCasp9_2IAP in component IAP (micromolar).
C RATES(18) is d/dt IAP in component IAP (micromolar).
C RATES(31) is d/dt ApopCasp9 in component ApopCasp9 (micromolar).
C RATES(27) is d/dt Casp9 in component Casp9 (micromolar).
C RATES(24) is d/dt ApopCasp9_2 in component ApopCasp9_2 (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0
      STATES(2) = 0
      STATES(3) = 0.004
      CONSTS(1) = 0.1
      CONSTS(2) = 10
      CONSTS(3) = 0.5
      CONSTS(4) = 0.006
      STATES(4) = 0.004
      STATES(5) = 0
      STATES(6) = 0
      CONSTS(5) = 0.0003
      CONSTS(6) = 5
      CONSTS(7) = 0.5
      STATES(7) = 0
      CONSTS(8) = 50000
      CONSTS(9) = 0.5
      CONSTS(10) = 4
      STATES(8) = 0.004
      STATES(9) = 0
      CONSTS(11) = 0.0003
      CONSTS(12) = 10
      STATES(10) = 0
      STATES(11) = 0.004
      CONSTS(13) = 10
      STATES(12) = 0
      CONSTS(14) = 10
      STATES(13) = 0
      CONSTS(15) = 10
      STATES(14) = 0.004
      CONSTS(16) = 10
      CONSTS(17) = 0.00003
      CONSTS(18) = 0.0066
      CONSTS(19) = 0.004
      CONSTS(20) = 0.006
      STATES(15) = 0
      STATES(16) = 0
      CONSTS(21) = 10
      CONSTS(22) = 0.5
      CONSTS(23) = 0.00008
      CONSTS(24) = 0.006
      CONSTS(25) = 0.1
      STATES(17) = 0
      CONSTS(26) = 10
      CONSTS(27) = 0.5
      CONSTS(28) = 0.1
      CONSTS(29) = 0.00003
      STATES(18) = 0.004
      STATES(19) = 0
      CONSTS(30) = 5
      CONSTS(31) = 0.0035
      STATES(20) = 0.004
      STATES(21) = 0
      STATES(22) = 0
      CONSTS(32) = 10
      CONSTS(33) = 0.5
      CONSTS(34) = 10
      CONSTS(35) = 0.5
      CONSTS(36) = 0.0003
      CONSTS(37) = 0.1
      STATES(23) = 0
      STATES(24) = 0
      STATES(25) = 0.004
      CONSTS(38) = 10
      CONSTS(39) = 0.5
      CONSTS(40) = 0.1
      STATES(26) = 0
      CONSTS(41) = 0.001
      STATES(27) = 0
      CONSTS(42) = 10
      CONSTS(43) = 0.5
      CONSTS(44) = 0.0003
      STATES(28) = 0
      STATES(29) = 0
      STATES(30) = 0
      STATES(31) = 0
      CONSTS(45) = 0.00003
      CONSTS(46) = 5
      CONSTS(47) = 0.0035
      CONSTS(48) = 5
      CONSTS(49) = 0.0035
      CONSTS(50) = 5
      CONSTS(51) = 0.0035
      CONSTS(52) = 5
      CONSTS(53) = 0.5
      CONSTS(54) = 5
      CONSTS(55) = 0.5
      CONSTS(56) =  CONSTS(17)*(1.00000+CONSTS(18) ** 4.00000/(CONSTS(18) ** 4.00000+CONSTS(19) ** 4.00000))
      CONSTS(57) = ( CONSTS(23)*CONSTS(19) ** 4.00000)/(CONSTS(18) ** 4.00000+CONSTS(19) ** 4.00000)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) =  CONSTS(2)*STATES(1)*STATES(3) -  CONSTS(3)*STATES(2)
      ALGBRC(1) =  CONSTS(1)*STATES(2)
      ALGBRC(2) =  - CONSTS(4)*STATES(1)
      RATES(1) = - ALGBRC(3)+ALGBRC(1)+ALGBRC(2)
      RATES(2) = ALGBRC(3) - ALGBRC(1)
      ALGBRC(4) = CONSTS(5) -  CONSTS(4)*STATES(4)
      ALGBRC(5) =  CONSTS(6)*STATES(6)*STATES(4) -  CONSTS(7)*STATES(5)
      RATES(4) = - ALGBRC(5)+ALGBRC(4)
      ALGBRC(7) =  CONSTS(8)*STATES(5) ** CONSTS(10) -  CONSTS(9)*STATES(7)
      RATES(5) = ALGBRC(5) -  7.00000*ALGBRC(7)
      ALGBRC(10) =  CONSTS(12)*STATES(9)*STATES(8)
      ALGBRC(6) =  - CONSTS(4)*STATES(6)
      RATES(6) = (ALGBRC(10) - ALGBRC(5))+ALGBRC(6)
      ALGBRC(8) = CONSTS(11) -  CONSTS(4)*STATES(8)
      RATES(8) = - ALGBRC(10)+ALGBRC(8)
      ALGBRC(13) =  CONSTS(14)*STATES(12)
      ALGBRC(12) =  CONSTS(13)*STATES(10)*STATES(11)
      ALGBRC(11) =  - CONSTS(4)*STATES(10)
      RATES(10) = (ALGBRC(13) - ALGBRC(12))+ALGBRC(11)
      ALGBRC(16) =  CONSTS(15)*STATES(13)*STATES(11)
      ALGBRC(9) =  - CONSTS(4)*STATES(9)
      RATES(9) = ALGBRC(16)+ALGBRC(9)
      ALGBRC(15) =  - CONSTS(4)*STATES(13)
      RATES(13) = (ALGBRC(12) - ALGBRC(16))+ALGBRC(15)
      ALGBRC(18) =  CONSTS(16)*STATES(14)*STATES(11)
      ALGBRC(17) = CONSTS(56) -  CONSTS(20)*STATES(11)
      RATES(11) = ((- ALGBRC(12) - ALGBRC(16)) - ALGBRC(18))+ALGBRC(17)
      ALGBRC(19) =  CONSTS(21)*STATES(15)*STATES(14) -  CONSTS(22)*STATES(16)
      ALGBRC(21) = CONSTS(57) -  CONSTS(24)*STATES(14)
      RATES(14) = (- ALGBRC(19) - ALGBRC(18))+ALGBRC(21)
      ALGBRC(20) =  CONSTS(25)*STATES(16)
      RATES(16) = ALGBRC(19) - ALGBRC(20)
      ALGBRC(24) =  CONSTS(28)*STATES(17)
      ALGBRC(14) =  - CONSTS(4)*STATES(12)
      RATES(12) = ((ALGBRC(1)+ALGBRC(24)) - ALGBRC(13))+ALGBRC(16)+ALGBRC(14)
      ALGBRC(22) =  CONSTS(26)*STATES(15)*STATES(3) -  CONSTS(27)*STATES(17)
      RATES(17) = ALGBRC(22) - ALGBRC(24)
      ALGBRC(23) = CONSTS(29) -  CONSTS(4)*STATES(3)
      RATES(3) = (- ALGBRC(3) - ALGBRC(22))+ALGBRC(23)
      ALGBRC(26) =  CONSTS(30)*STATES(15)*STATES(18) -  CONSTS(31)*STATES(19)
      RATES(19) = ALGBRC(26)
      ALGBRC(28) =  CONSTS(32)*STATES(7)*STATES(20) -  CONSTS(33)*STATES(21)
      ALGBRC(29) =  CONSTS(34)*STATES(21)*STATES(20) -  CONSTS(35)*STATES(22)
      ALGBRC(27) = CONSTS(36) -  CONSTS(4)*STATES(20)
      RATES(20) = (- ALGBRC(28) - ALGBRC(29))+ALGBRC(27)
      RATES(21) = ALGBRC(28) - ALGBRC(29)
      ALGBRC(30) =  CONSTS(37)*STATES(22)
      RATES(22) = ALGBRC(29) - ALGBRC(30)
      ALGBRC(32) =  CONSTS(40)*STATES(23)
      ALGBRC(31) =  CONSTS(38)*STATES(24)*STATES(25) -  CONSTS(39)*STATES(23)
      RATES(23) = ALGBRC(31) - ALGBRC(32)
      ALGBRC(33) =  CONSTS(41)*STATES(26)
      ALGBRC(25) =  - CONSTS(4)*STATES(15)
      RATES(15) = (((((ALGBRC(33)+ALGBRC(32)) - ALGBRC(26)) - ALGBRC(22))+ALGBRC(24)) - ALGBRC(19))+ALGBRC(20)+ALGBRC(25)
      ALGBRC(34) =  CONSTS(42)*STATES(27)*STATES(25) -  CONSTS(43)*STATES(26)
      RATES(26) = ALGBRC(34) - ALGBRC(33)
      ALGBRC(35) = CONSTS(44) -  CONSTS(4)*STATES(25)
      RATES(25) = (- ALGBRC(34) - ALGBRC(31))+ALGBRC(35)
      ALGBRC(36) =  CONSTS(46)*STATES(27)*STATES(18) -  CONSTS(47)*STATES(28)
      RATES(28) = ALGBRC(36)
      ALGBRC(37) =  CONSTS(48)*STATES(31)*STATES(18) -  CONSTS(49)*STATES(29)
      RATES(29) = ALGBRC(37)
      ALGBRC(38) =  CONSTS(50)*STATES(24)*STATES(18) -  CONSTS(51)*STATES(30)
      RATES(30) = ALGBRC(38)
      ALGBRC(40) = CONSTS(45) -  CONSTS(4)*STATES(18)
      RATES(18) = (((- ALGBRC(36) - ALGBRC(37)) - ALGBRC(38)) - ALGBRC(26))+ALGBRC(40)
      ALGBRC(39) =  CONSTS(52)*STATES(24) -  CONSTS(53)*STATES(31)*STATES(27)
      RATES(24) = (((ALGBRC(30) - ALGBRC(39)) - ALGBRC(38)) - ALGBRC(31))+ALGBRC(32)
      ALGBRC(41) =  CONSTS(54)*STATES(31) -  CONSTS(55)*STATES(7)*STATES(27)
      RATES(7) = (ALGBRC(7) - ALGBRC(28))+ALGBRC(41)
      RATES(31) = (ALGBRC(39) - ALGBRC(41)) - ALGBRC(37)
      ALGBRC(42) =  - CONSTS(4)*STATES(27)
      RATES(27) = (((ALGBRC(39)+ALGBRC(41)) - ALGBRC(36)) - ALGBRC(34))+ALGBRC(33)+ALGBRC(42)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) =  CONSTS(2)*STATES(1)*STATES(3) -  CONSTS(3)*STATES(2)
      ALGBRC(1) =  CONSTS(1)*STATES(2)
      ALGBRC(2) =  - CONSTS(4)*STATES(1)
      ALGBRC(4) = CONSTS(5) -  CONSTS(4)*STATES(4)
      ALGBRC(5) =  CONSTS(6)*STATES(6)*STATES(4) -  CONSTS(7)*STATES(5)
      ALGBRC(7) =  CONSTS(8)*STATES(5) ** CONSTS(10) -  CONSTS(9)*STATES(7)
      ALGBRC(10) =  CONSTS(12)*STATES(9)*STATES(8)
      ALGBRC(6) =  - CONSTS(4)*STATES(6)
      ALGBRC(8) = CONSTS(11) -  CONSTS(4)*STATES(8)
      ALGBRC(13) =  CONSTS(14)*STATES(12)
      ALGBRC(12) =  CONSTS(13)*STATES(10)*STATES(11)
      ALGBRC(11) =  - CONSTS(4)*STATES(10)
      ALGBRC(16) =  CONSTS(15)*STATES(13)*STATES(11)
      ALGBRC(9) =  - CONSTS(4)*STATES(9)
      ALGBRC(15) =  - CONSTS(4)*STATES(13)
      ALGBRC(18) =  CONSTS(16)*STATES(14)*STATES(11)
      ALGBRC(17) = CONSTS(56) -  CONSTS(20)*STATES(11)
      ALGBRC(19) =  CONSTS(21)*STATES(15)*STATES(14) -  CONSTS(22)*STATES(16)
      ALGBRC(21) = CONSTS(57) -  CONSTS(24)*STATES(14)
      ALGBRC(20) =  CONSTS(25)*STATES(16)
      ALGBRC(24) =  CONSTS(28)*STATES(17)
      ALGBRC(14) =  - CONSTS(4)*STATES(12)
      ALGBRC(22) =  CONSTS(26)*STATES(15)*STATES(3) -  CONSTS(27)*STATES(17)
      ALGBRC(23) = CONSTS(29) -  CONSTS(4)*STATES(3)
      ALGBRC(26) =  CONSTS(30)*STATES(15)*STATES(18) -  CONSTS(31)*STATES(19)
      ALGBRC(28) =  CONSTS(32)*STATES(7)*STATES(20) -  CONSTS(33)*STATES(21)
      ALGBRC(29) =  CONSTS(34)*STATES(21)*STATES(20) -  CONSTS(35)*STATES(22)
      ALGBRC(27) = CONSTS(36) -  CONSTS(4)*STATES(20)
      ALGBRC(30) =  CONSTS(37)*STATES(22)
      ALGBRC(32) =  CONSTS(40)*STATES(23)
      ALGBRC(31) =  CONSTS(38)*STATES(24)*STATES(25) -  CONSTS(39)*STATES(23)
      ALGBRC(33) =  CONSTS(41)*STATES(26)
      ALGBRC(25) =  - CONSTS(4)*STATES(15)
      ALGBRC(34) =  CONSTS(42)*STATES(27)*STATES(25) -  CONSTS(43)*STATES(26)
      ALGBRC(35) = CONSTS(44) -  CONSTS(4)*STATES(25)
      ALGBRC(36) =  CONSTS(46)*STATES(27)*STATES(18) -  CONSTS(47)*STATES(28)
      ALGBRC(37) =  CONSTS(48)*STATES(31)*STATES(18) -  CONSTS(49)*STATES(29)
      ALGBRC(38) =  CONSTS(50)*STATES(24)*STATES(18) -  CONSTS(51)*STATES(30)
      ALGBRC(40) = CONSTS(45) -  CONSTS(4)*STATES(18)
      ALGBRC(39) =  CONSTS(52)*STATES(24) -  CONSTS(53)*STATES(31)*STATES(27)
      ALGBRC(41) =  CONSTS(54)*STATES(31) -  CONSTS(55)*STATES(7)*STATES(27)
      ALGBRC(42) =  - CONSTS(4)*STATES(27)
      RETURN
      END