C
C There are a total of 66 entries in the algebraic variable array.
C There are a total of 42 entries in each of the rate and state variable arrays.
C There are a total of 94 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 FeLnNO in component FeLn (micromolar).
C STATES(4) 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 ALGBRC(4) is r_20NO in component Casp8 (flux).
C ALGBRC(62) is r_19NO in component PTPC (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 CONSTS(5) is k_20NO in component model_constant (second_order_rate_constant).
C STATES(5) is Apaf_1 in component Apaf_1 (micromolar).
C STATES(6) is CytcApaf_1 in component CytcApaf_1 (micromolar).
C STATES(7) is Cytc in component Cytc (micromolar).
C ALGBRC(5) is J_Apaf_1 in component Apaf_1 (flux).
C ALGBRC(6) is J_1 in component Apaf_1 (flux).
C CONSTS(6) is P_Apaf_1 in component model_constant (flux).
C CONSTS(7) is k_11 in component model_constant (second_order_rate_constant).
C CONSTS(8) is k_01 in component model_constant (first_order_rate_constant).
C STATES(8) is Apop in component Apop (micromolar).
C ALGBRC(8) is J_1b in component CytcApaf_1 (flux).
C CONSTS(9) is k_11b in component model_constant (rate).
C CONSTS(10) is k_01b in component model_constant (first_order_rate_constant).
C CONSTS(11) is p in component model_constant (dimensionless).
C STATES(9) is Cytc_mito in component Cytc_mito (micromolar).
C STATES(10) is PTPC in component PTPC (micromolar).
C ALGBRC(9) is PTPC_act in component Cytc (micromolar).
C CONSTS(12) is PTPC_0 in component model_constant (micromolar).
C ALGBRC(13) is J_14 in component Cytc_mito (flux).
C ALGBRC(7) is J_Cytc in component Cytc (flux).
C CONSTS(13) is k_1 in component model_constant (second_order_rate_constant).
C STATES(11) is Bax_2 in component Bax_2 (micromolar).
C ALGBRC(11) is J_Cytc_mito in component Cytc_mito (flux).
C CONSTS(14) is P_Cytc_mito in component model_constant (flux).
C CONSTS(15) is k14 in component model_constant (second_order_rate_constant).
C ALGBRC(18) is J_12b in component tBidBax (flux).
C ALGBRC(10) is J_Bax_2 in component Bax_2 (flux).
C STATES(12) is tBid_mito in component tBid_mito (micromolar).
C STATES(13) is Bax in component Bax (micromolar).
C ALGBRC(15) is J_11 in component tBid (flux).
C ALGBRC(14) is J_12a in component tBid_mito (flux).
C ALGBRC(12) is J_tBid_mito in component tBid_mito (flux).
C CONSTS(16) is k12a in component model_constant (second_order_rate_constant).
C STATES(14) is tBid in component tBid (micromolar).
C ALGBRC(26) is J_f8 in component Casp3Bid (flux).
C ALGBRC(16) is J_tBid in component tBid (flux).
C CONSTS(17) is k11 in component model_constant (first_order_rate_constant).
C STATES(15) is tBidBax in component tBidBax (micromolar).
C ALGBRC(17) is J_tBidBax in component tBidBax (flux).
C CONSTS(18) is k12b in component model_constant (second_order_rate_constant).
C STATES(16) is Bcl_2 in component Bcl_2 (micromolar).
C ALGBRC(20) is J_13 in component Bax (flux).
C ALGBRC(19) is J_Bax in component Bax (flux).
C CONSTS(92) is P_Bax in component Bax (flux).
C CONSTS(19) is k13 in component model_constant (second_order_rate_constant).
C CONSTS(20) is P_oBax in component model_constant (flux).
C CONSTS(21) is p53 in component model_constant (micromolar).
C CONSTS(22) is p53_thresh in component model_constant (micromolar).
C CONSTS(23) is u_Bax in component model_constant (first_order_rate_constant).
C STATES(17) is Casp3 in component Casp3 (micromolar).
C STATES(18) is Casp3Bcl_2 in component Casp3Bcl_2 (micromolar).
C ALGBRC(21) is J_9 in component Bcl_2 (flux).
C ALGBRC(23) is J_Bcl_2 in component Bcl_2 (flux).
C CONSTS(24) is k_19 in component model_constant (second_order_rate_constant).
C CONSTS(25) is k_09 in component model_constant (first_order_rate_constant).
C CONSTS(93) is P_Bcl_2 in component Bcl_2 (flux).
C CONSTS(26) is P_oBcl_2 in component model_constant (flux).
C CONSTS(27) is u_Bcl_2 in component model_constant (first_order_rate_constant).
C ALGBRC(22) is J_f9 in component Casp3Bcl_2 (flux).
C CONSTS(28) is k_f9 in component model_constant (first_order_rate_constant).
C STATES(19) is Casp3Bid in component Casp3Bid (micromolar).
C ALGBRC(24) is J_8 in component Casp3Bid (flux).
C CONSTS(29) is k_18 in component model_constant (second_order_rate_constant).
C CONSTS(30) is k_08 in component model_constant (first_order_rate_constant).
C CONSTS(31) is k_f8 in component model_constant (first_order_rate_constant).
C ALGBRC(25) is J_Bid in component Bid (flux).
C CONSTS(32) is P_Bid in component model_constant (flux).
C STATES(20) is IAP in component IAP (micromolar).
C STATES(21) is Casp3IAP in component Casp3 (micromolar).
C ALGBRC(35) is J_f6 in component Casp9Pro3 (flux).
C ALGBRC(34) is J_f6b in component ApopCasp9_2Pro3 (flux).
C ALGBRC(28) is J_7 in component Casp3 (flux).
C ALGBRC(27) is J_Casp3 in component Casp3 (flux).
C ALGBRC(66) is r_22NO in component FeLn (flux).
C CONSTS(33) is k_17 in component model_constant (second_order_rate_constant).
C CONSTS(34) is k_07 in component model_constant (first_order_rate_constant).
C ALGBRC(30) is J_2 in component Pro9 (flux).
C ALGBRC(43) is J_4b in component ApopCasp9 (flux).
C STATES(22) is Pro9 in component Pro9 (micromolar).
C STATES(23) is ApopPro9 in component Pro9 (micromolar).
C STATES(24) is ApopPro9_2 in component ApopPro9_2 (micromolar).
C ALGBRC(31) is J_3 in component Pro9 (flux).
C ALGBRC(29) is J_Pro9 in component Pro9 (flux).
C CONSTS(35) is k_12 in component model_constant (second_order_rate_constant).
C CONSTS(36) is k_02 in component model_constant (first_order_rate_constant).
C CONSTS(37) is k_13 in component model_constant (second_order_rate_constant).
C CONSTS(38) is k_03 in component model_constant (first_order_rate_constant).
C CONSTS(39) is P_Pro9 in component model_constant (flux).
C ALGBRC(32) is J_f3 in component ApopPro9_2 (flux).
C CONSTS(40) is k_f3 in component model_constant (first_order_rate_constant).
C STATES(25) is ApopCasp9_2Pro3 in component ApopCasp9_2Pro3 (micromolar).
C STATES(26) is ApopCasp9_2 in component ApopCasp9_2 (micromolar).
C STATES(27) is Pro3 in component Pro3 (micromolar).
C ALGBRC(33) is J_6b in component ApopCasp9_2Pro3 (flux).
C CONSTS(41) is k_16b in component model_constant (second_order_rate_constant).
C CONSTS(42) is k_06b in component model_constant (first_order_rate_constant).
C CONSTS(43) is k_f6b in component model_constant (first_order_rate_constant).
C STATES(28) is Casp9Pro3 in component Casp9Pro3 (micromolar).
C ALGBRC(36) is J_6 in component Pro3 (flux).
C CONSTS(44) is k_f6 in component model_constant (first_order_rate_constant).
C STATES(29) is Casp9 in component Casp9 (micromolar).
C ALGBRC(37) is J_Pro3 in component Pro3 (flux).
C CONSTS(45) is k_16 in component model_constant (second_order_rate_constant).
C CONSTS(46) is k_06 in component model_constant (first_order_rate_constant).
C CONSTS(47) is P_Pro3 in component model_constant (flux).
C STATES(30) is Casp9IAP in component IAP (micromolar).
C STATES(31) is ApopCasp9IAP in component IAP (micromolar).
C STATES(32) is ApopCasp9_2IAP in component IAP (micromolar).
C STATES(33) is ApopCasp9 in component ApopCasp9 (micromolar).
C ALGBRC(38) is J_5 in component IAP (flux).
C ALGBRC(39) is J_5b in component IAP (flux).
C ALGBRC(40) is J_5c in component IAP (flux).
C ALGBRC(42) is J_IAP in component IAP (flux).
C CONSTS(48) is P_IAP in component model_constant (flux).
C CONSTS(49) is k_15 in component model_constant (second_order_rate_constant).
C CONSTS(50) is k_05 in component model_constant (first_order_rate_constant).
C CONSTS(51) is k_15b in component model_constant (second_order_rate_constant).
C CONSTS(52) is k_05b in component model_constant (first_order_rate_constant).
C CONSTS(53) is k_15c in component model_constant (second_order_rate_constant).
C CONSTS(54) is k_05c in component model_constant (first_order_rate_constant).
C ALGBRC(41) is J_4 in component ApopCasp9 (flux).
C CONSTS(55) is k_14 in component model_constant (first_order_rate_constant).
C CONSTS(56) is k_04 in component model_constant (second_order_rate_constant).
C CONSTS(57) is k_14b in component model_constant (first_order_rate_constant).
C CONSTS(58) is k_04b in component model_constant (second_order_rate_constant).
C ALGBRC(44) is J_Casp9 in component Casp9 (flux).
C ALGBRC(65) is r_21NO in component FeLn (flux).
C STATES(34) is NO in component NO (micromolar).
C STATES(35) is O_2m in component O_2m (micromolar).
C CONSTS(59) is O_2 in component NO (micromolar).
C STATES(36) is NO_2 in component NO (micromolar).
C STATES(37) is N2O3 in component N2O3 (micromolar).
C STATES(38) is GSNO in component GSNO (micromolar).
C STATES(39) is CcOX in component NO (micromolar).
C STATES(40) is FeLn in component FeLn (micromolar).
C CONSTS(91) is r_1NO in component NO (flux).
C ALGBRC(45) is r_4NO in component NO (flux).
C ALGBRC(46) is r_12aNO in component NO (flux).
C ALGBRC(47) is r_12bNOp in component NO (flux).
C ALGBRC(48) is r_12bNOm in component NO (flux).
C ALGBRC(49) is r_14NO in component NO (flux).
C ALGBRC(50) is r_15NO in component NO (flux).
C ALGBRC(51) is r_16NO in component NO (flux).
C CONSTS(60) is k_1NO in component model_constant (flux).
C CONSTS(61) is k_4NO in component model_constant (second_order_rate_constant).
C CONSTS(62) is k_12aNO in component model_constant (rate2).
C CONSTS(63) is k_12bNOp in component model_constant (second_order_rate_constant).
C CONSTS(64) is k_12bNOm in component model_constant (first_order_rate_constant).
C CONSTS(65) is k_14NO in component model_constant (first_order_rate_constant).
C CONSTS(66) is k_15NO in component model_constant (second_order_rate_constant).
C CONSTS(67) is k_16NO in component model_constant (second_order_rate_constant).
C CONSTS(68) is SOD in component O_2m (micromolar).
C CONSTS(94) is r_2NO in component O_2m (flux).
C ALGBRC(52) is r_5NO in component O_2m (flux).
C ALGBRC(53) is r_10NO in component O_2m (flux).
C CONSTS(69) is k_2NO in component model_constant (flux).
C CONSTS(70) is k_5NO in component model_constant (second_order_rate_constant).
C CONSTS(71) is k_10NO in component model_constant (rate2).
C STATES(41) is ONOO_m in component ONOO_m (micromolar).
C STATES(42) is GSH in component GSH (micromolar).
C CONSTS(72) is GPX in component ONOO_m (micromolar).
C CONSTS(73) is CO_2 in component ONOO_m (micromolar).
C CONSTS(74) is Cyt_c in component ONOO_m (micromolar).
C ALGBRC(54) is r_6NO in component ONOO_m (flux).
C ALGBRC(55) is r_7NO in component ONOO_m (flux).
C ALGBRC(58) is r_8NO in component ONOO_m (flux).
C ALGBRC(61) is r_9NO in component ONOO_m (flux).
C ALGBRC(64) is r_18NO in component ONOO_m (flux).
C CONSTS(75) is k_6NO in component model_constant (second_order_rate_constant).
C CONSTS(76) is k_7NO in component model_constant (second_order_rate_constant).
C CONSTS(77) is k_8NO in component model_constant (second_order_rate_constant).
C CONSTS(78) is k_9NO in component model_constant (second_order_rate_constant).
C CONSTS(79) is k_18NO in component model_constant (second_order_rate_constant).
C ALGBRC(56) is GSSG in component GSH (micromolar).
C CONSTS(80) is FeLn_0 in component model_constant (micromolar).
C CONSTS(81) is GSH_0 in component model_constant (micromolar).
C ALGBRC(57) is r_11NO in component GSH (flux).
C ALGBRC(60) is r_m in component GSH (flux).
C ALGBRC(63) is r_17NO in component GSH (flux).
C CONSTS(82) is k_11NO in component model_constant (second_order_rate_constant).
C CONSTS(83) is v_m in component model_constant (flux).
C CONSTS(84) is k_m in component model_constant (micromolar).
C CONSTS(85) is k_17NO in component model_constant (second_order_rate_constant).
C ALGBRC(59) is r_13NO in component N2O3 (flux).
C CONSTS(86) is k_13NO in component model_constant (first_order_rate_constant).
C CONSTS(87) is k_21NO in component model_constant (second_order_rate_constant).
C CONSTS(88) is k_22NO in component model_constant (second_order_rate_constant).
C CONSTS(89) is k_19NO in component model_constant (second_order_rate_constant).
C CONSTS(90) is k_17bNO in component model_constant (second_order_rate_constant).
C RATES(1) is d/dt Casp8 in component Casp8 (micromolar).
C RATES(2) is d/dt Casp8Bid in component Casp8 (micromolar).
C RATES(5) is d/dt Apaf_1 in component Apaf_1 (micromolar).
C RATES(6) is d/dt CytcApaf_1 in component CytcApaf_1 (micromolar).
C RATES(7) is d/dt Cytc in component Cytc (micromolar).
C RATES(9) is d/dt Cytc_mito in component Cytc_mito (micromolar).
C RATES(11) is d/dt Bax_2 in component Bax_2 (micromolar).
C RATES(12) is d/dt tBid_mito in component tBid_mito (micromolar).
C RATES(14) is d/dt tBid in component tBid (micromolar).
C RATES(15) is d/dt tBidBax in component tBidBax (micromolar).
C RATES(13) is d/dt Bax in component Bax (micromolar).
C RATES(16) is d/dt Bcl_2 in component Bcl_2 (micromolar).
C RATES(18) is d/dt Casp3Bcl_2 in component Casp3Bcl_2 (micromolar).
C RATES(19) is d/dt Casp3Bid in component Casp3Bid (micromolar).
C RATES(4) is d/dt Bid in component Bid (micromolar).
C RATES(21) is d/dt Casp3IAP in component Casp3 (micromolar).
C RATES(17) is d/dt Casp3 in component Casp3 (micromolar).
C RATES(8) is d/dt Apop in component Apop (micromolar).
C RATES(22) is d/dt Pro9 in component Pro9 (micromolar).
C RATES(23) is d/dt ApopPro9 in component Pro9 (micromolar).
C RATES(24) is d/dt ApopPro9_2 in component ApopPro9_2 (micromolar).
C RATES(25) is d/dt ApopCasp9_2Pro3 in component ApopCasp9_2Pro3 (micromolar).
C RATES(28) is d/dt Casp9Pro3 in component Casp9Pro3 (micromolar).
C RATES(27) is d/dt Pro3 in component Pro3 (micromolar).
C RATES(30) is d/dt Casp9IAP in component IAP (micromolar).
C RATES(31) is d/dt ApopCasp9IAP in component IAP (micromolar).
C RATES(32) is d/dt ApopCasp9_2IAP in component IAP (micromolar).
C RATES(20) is d/dt IAP in component IAP (micromolar).
C RATES(33) is d/dt ApopCasp9 in component ApopCasp9 (micromolar).
C RATES(29) is d/dt Casp9 in component Casp9 (micromolar).
C RATES(26) is d/dt ApopCasp9_2 in component ApopCasp9_2 (micromolar).
C RATES(34) is d/dt NO in component NO (micromolar).
C RATES(39) is d/dt CcOX in component NO (micromolar).
C RATES(36) is d/dt NO_2 in component NO (micromolar).
C RATES(35) is d/dt O_2m in component O_2m (micromolar).
C RATES(41) is d/dt ONOO_m in component ONOO_m (micromolar).
C RATES(42) is d/dt GSH in component GSH (micromolar).
C RATES(38) is d/dt GSNO in component GSNO (micromolar).
C RATES(37) is d/dt N2O3 in component N2O3 (micromolar).
C RATES(40) is d/dt FeLn in component FeLn (micromolar).
C RATES(3) is d/dt FeLnNO in component FeLn (micromolar).
C RATES(10) is d/dt PTPC in component PTPC (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0
      STATES(2) = 0
      STATES(3) = 0
      STATES(4) = 0.004
      CONSTS(1) = 0.1
      CONSTS(2) = 10
      CONSTS(3) = 0.5
      CONSTS(4) = 0.006
      CONSTS(5) = 66
      STATES(5) = 0.004
      STATES(6) = 0
      STATES(7) = 0
      CONSTS(6) = 0.0003
      CONSTS(7) = 5
      CONSTS(8) = 0.5
      STATES(8) = 0
      CONSTS(9) = 50000
      CONSTS(10) = 0.5
      CONSTS(11) = 4
      STATES(9) = 0.004
      STATES(10) = 0
      CONSTS(12) = 0
      CONSTS(13) = 1
      STATES(11) = 0
      CONSTS(14) = 0.0003
      CONSTS(15) = 10
      STATES(12) = 0
      STATES(13) = 0.004
      CONSTS(16) = 10
      STATES(14) = 0
      CONSTS(17) = 10
      STATES(15) = 0
      CONSTS(18) = 10
      STATES(16) = 0.004
      CONSTS(19) = 10
      CONSTS(20) = 0.00003
      CONSTS(21) = 0.0066
      CONSTS(22) = 0.004
      CONSTS(23) = 0.006
      STATES(17) = 0
      STATES(18) = 0
      CONSTS(24) = 10
      CONSTS(25) = 0.5
      CONSTS(26) = 0.00008
      CONSTS(27) = 0.006
      CONSTS(28) = 0.1
      STATES(19) = 0
      CONSTS(29) = 10
      CONSTS(30) = 0.5
      CONSTS(31) = 0.1
      CONSTS(32) = 0.00003
      STATES(20) = 0.004
      STATES(21) = 0
      CONSTS(33) = 5
      CONSTS(34) = 0.0035
      STATES(22) = 0.004
      STATES(23) = 0
      STATES(24) = 0
      CONSTS(35) = 10
      CONSTS(36) = 0.5
      CONSTS(37) = 10
      CONSTS(38) = 0.5
      CONSTS(39) = 0.0003
      CONSTS(40) = 0.1
      STATES(25) = 0
      STATES(26) = 0
      STATES(27) = 0.004
      CONSTS(41) = 10
      CONSTS(42) = 0.5
      CONSTS(43) = 0.1
      STATES(28) = 0
      CONSTS(44) = 0.001
      STATES(29) = 0
      CONSTS(45) = 10
      CONSTS(46) = 0.5
      CONSTS(47) = 0.0003
      STATES(30) = 0
      STATES(31) = 0
      STATES(32) = 0
      STATES(33) = 0
      CONSTS(48) = 0.00003
      CONSTS(49) = 5
      CONSTS(50) = 0.0035
      CONSTS(51) = 5
      CONSTS(52) = 0.0035
      CONSTS(53) = 5
      CONSTS(54) = 0.0035
      CONSTS(55) = 5
      CONSTS(56) = 0.5
      CONSTS(57) = 5
      CONSTS(58) = 0.5
      STATES(34) = 0
      STATES(35) = 35
      CONSTS(59) = 35
      STATES(36) = 0
      STATES(37) = 0
      STATES(38) = 0
      STATES(39) = 0.1
      STATES(40) = 0.05
      CONSTS(60) = 1
      CONSTS(61) = 6700
      CONSTS(62) = 0.000006
      CONSTS(63) = 1100
      CONSTS(64) = 81000
      CONSTS(65) = 0.0002
      CONSTS(66) = 100
      CONSTS(67) = 1.21
      CONSTS(68) = 10
      CONSTS(69) = 0.1
      CONSTS(70) = 2400
      CONSTS(71) = 0.0006
      STATES(41) = 0
      STATES(42) = 10000
      CONSTS(72) = 5.8
      CONSTS(73) = 1000
      CONSTS(74) = 400
      CONSTS(75) = 0.00135
      CONSTS(76) = 2
      CONSTS(77) = 0.058
      CONSTS(78) = 0.025
      CONSTS(79) = 1
      CONSTS(80) = 0.05
      CONSTS(81) = 10000
      CONSTS(82) = 66
      CONSTS(83) = 320
      CONSTS(84) = 50
      CONSTS(85) = 66
      CONSTS(86) = 1600
      CONSTS(87) = 66
      CONSTS(88) = 66
      CONSTS(89) = 10
      CONSTS(90) = 0.0002
      CONSTS(91) = CONSTS(60)
      CONSTS(92) =  CONSTS(20)*(1.00000+CONSTS(21) ** 4.00000/(CONSTS(21) ** 4.00000+CONSTS(22) ** 4.00000))
      CONSTS(93) = ( CONSTS(26)*CONSTS(22) ** 4.00000)/(CONSTS(21) ** 4.00000+CONSTS(22) ** 4.00000)
      CONSTS(94) = CONSTS(69)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) =  CONSTS(2)*STATES(1)*STATES(4) -  CONSTS(3)*STATES(2)
      ALGBRC(1) =  CONSTS(1)*STATES(2)
      RATES(2) = ALGBRC(3) - ALGBRC(1)
      ALGBRC(5) = CONSTS(6) -  CONSTS(4)*STATES(5)
      ALGBRC(6) =  CONSTS(7)*STATES(7)*STATES(5) -  CONSTS(8)*STATES(6)
      RATES(5) = - ALGBRC(6)+ALGBRC(5)
      ALGBRC(8) =  CONSTS(9)*STATES(6) ** CONSTS(11) -  CONSTS(10)*STATES(8)
      RATES(6) = ALGBRC(6) -  7.00000*ALGBRC(8)
      ALGBRC(9) = CONSTS(12) - STATES(10)
      ALGBRC(13) =  CONSTS(15)*STATES(11)*STATES(9)
      ALGBRC(7) =  - CONSTS(4)*STATES(7)
      RATES(7) = (ALGBRC(13) - ALGBRC(6))+ALGBRC(7)+ CONSTS(13)*ALGBRC(9)*STATES(9)
      ALGBRC(11) = CONSTS(14) -  CONSTS(4)*STATES(9)
      RATES(9) = - ALGBRC(13)+ALGBRC(11)
      ALGBRC(15) =  CONSTS(17)*STATES(14)
      ALGBRC(14) =  CONSTS(16)*STATES(12)*STATES(13)
      ALGBRC(12) =  - CONSTS(4)*STATES(12)
      RATES(12) = (ALGBRC(15) - ALGBRC(14))+ALGBRC(12)
      ALGBRC(18) =  CONSTS(18)*STATES(15)*STATES(13)
      ALGBRC(10) =  - CONSTS(4)*STATES(11)
      RATES(11) = ALGBRC(18)+ALGBRC(10)
      ALGBRC(17) =  - CONSTS(4)*STATES(15)
      RATES(15) = (ALGBRC(14) - ALGBRC(18))+ALGBRC(17)
      ALGBRC(20) =  CONSTS(19)*STATES(16)*STATES(13)
      ALGBRC(19) = CONSTS(92) -  CONSTS(23)*STATES(13)
      RATES(13) = ((- ALGBRC(14) - ALGBRC(18)) - ALGBRC(20))+ALGBRC(19)
      ALGBRC(21) =  CONSTS(24)*STATES(17)*STATES(16) -  CONSTS(25)*STATES(18)
      ALGBRC(23) = CONSTS(93) -  CONSTS(27)*STATES(16)
      RATES(16) = (- ALGBRC(21) - ALGBRC(20))+ALGBRC(23)
      ALGBRC(22) =  CONSTS(28)*STATES(18)
      RATES(18) = ALGBRC(21) - ALGBRC(22)
      ALGBRC(26) =  CONSTS(31)*STATES(19)
      ALGBRC(16) =  - CONSTS(4)*STATES(14)
      RATES(14) = ((ALGBRC(1)+ALGBRC(26)) - ALGBRC(15))+ALGBRC(18)+ALGBRC(16)
      ALGBRC(24) =  CONSTS(29)*STATES(17)*STATES(4) -  CONSTS(30)*STATES(19)
      RATES(19) = ALGBRC(24) - ALGBRC(26)
      ALGBRC(25) = CONSTS(32) -  CONSTS(4)*STATES(4)
      RATES(4) = (- ALGBRC(3) - ALGBRC(24))+ALGBRC(25)
      ALGBRC(28) =  CONSTS(33)*STATES(17)*STATES(20) -  CONSTS(34)*STATES(21)
      RATES(21) = ALGBRC(28)
      ALGBRC(30) =  CONSTS(35)*STATES(8)*STATES(22) -  CONSTS(36)*STATES(23)
      ALGBRC(31) =  CONSTS(37)*STATES(23)*STATES(22) -  CONSTS(38)*STATES(24)
      ALGBRC(29) = CONSTS(39) -  CONSTS(4)*STATES(22)
      RATES(22) = (- ALGBRC(30) - ALGBRC(31))+ALGBRC(29)
      RATES(23) = ALGBRC(30) - ALGBRC(31)
      ALGBRC(32) =  CONSTS(40)*STATES(24)
      RATES(24) = ALGBRC(31) - ALGBRC(32)
      ALGBRC(34) =  CONSTS(43)*STATES(25)
      ALGBRC(33) =  CONSTS(41)*STATES(26)*STATES(27) -  CONSTS(42)*STATES(25)
      RATES(25) = ALGBRC(33) - ALGBRC(34)
      ALGBRC(35) =  CONSTS(44)*STATES(28)
      ALGBRC(36) =  CONSTS(45)*STATES(29)*STATES(27) -  CONSTS(46)*STATES(28)
      RATES(28) = ALGBRC(36) - ALGBRC(35)
      ALGBRC(37) = CONSTS(47) -  CONSTS(4)*STATES(27)
      RATES(27) = (- ALGBRC(36) - ALGBRC(33))+ALGBRC(37)
      ALGBRC(38) =  CONSTS(49)*STATES(29)*STATES(20) -  CONSTS(50)*STATES(30)
      RATES(30) = ALGBRC(38)
      ALGBRC(39) =  CONSTS(51)*STATES(33)*STATES(20) -  CONSTS(52)*STATES(31)
      RATES(31) = ALGBRC(39)
      ALGBRC(40) =  CONSTS(53)*STATES(26)*STATES(20) -  CONSTS(54)*STATES(32)
      RATES(32) = ALGBRC(40)
      ALGBRC(42) = CONSTS(48) -  CONSTS(4)*STATES(20)
      RATES(20) = (((- ALGBRC(38) - ALGBRC(39)) - ALGBRC(40)) - ALGBRC(28))+ALGBRC(42)
      ALGBRC(41) =  CONSTS(55)*STATES(26) -  CONSTS(56)*STATES(33)*STATES(29)
      RATES(26) = (((ALGBRC(32) - ALGBRC(41)) - ALGBRC(40)) - ALGBRC(33))+ALGBRC(34)
      ALGBRC(43) =  CONSTS(57)*STATES(33) -  CONSTS(58)*STATES(8)*STATES(29)
      RATES(8) = (ALGBRC(8) - ALGBRC(30))+ALGBRC(43)
      RATES(33) = (ALGBRC(41) - ALGBRC(43)) - ALGBRC(39)
      ALGBRC(46) =  CONSTS(62)*STATES(34) ** 2.00000*CONSTS(59)
      ALGBRC(47) =  CONSTS(63)*STATES(36)*STATES(34)
      ALGBRC(48) =  CONSTS(64)*STATES(37)
      RATES(36) = ( 2.00000*ALGBRC(46) - ALGBRC(47))+ALGBRC(48)
      ALGBRC(50) =  CONSTS(66)*STATES(39)*STATES(34)
      RATES(39) = - ALGBRC(50)
      ALGBRC(45) =  CONSTS(61)*STATES(34)*STATES(35)
      ALGBRC(49) =  CONSTS(65)*STATES(38)
      ALGBRC(51) =  CONSTS(67)*STATES(40)*STATES(34)
      RATES(34) = (((((CONSTS(91) - ALGBRC(45)) -  2.00000*ALGBRC(46)) - ALGBRC(47))+ALGBRC(48)+ALGBRC(49)) - ALGBRC(50)) - ALGBRC(51)
      ALGBRC(52) =  CONSTS(70)*CONSTS(68)*STATES(35)
      ALGBRC(53) =  CONSTS(71)*STATES(38) ** 2.00000*STATES(35)
      RATES(35) = ((CONSTS(94) - ALGBRC(45)) - ALGBRC(52)) - ALGBRC(53)
      ALGBRC(2) =  - CONSTS(4)*STATES(1)
      ALGBRC(4) =  CONSTS(5)*STATES(3)*STATES(1)
      ALGBRC(62) =  CONSTS(89)*STATES(37)*STATES(1)
      RATES(1) = ((- ALGBRC(3)+ALGBRC(1)+ALGBRC(2)) - ALGBRC(62)) - ALGBRC(4)
      ALGBRC(54) =  CONSTS(75)*STATES(41)*STATES(42)
      ALGBRC(55) =  CONSTS(76)*STATES(41)*CONSTS(72)
      ALGBRC(58) =  CONSTS(77)*STATES(41)*CONSTS(73)
      ALGBRC(61) =  CONSTS(78)*STATES(41)*CONSTS(74)
      ALGBRC(64) =  CONSTS(79)*STATES(41)*STATES(10)
      RATES(41) = ((((ALGBRC(45) - ALGBRC(54)) - ALGBRC(55)) - ALGBRC(58)) - ALGBRC(61)) - ALGBRC(64)
      ALGBRC(57) =  CONSTS(82)*STATES(37)*STATES(42)
      ALGBRC(56) = ((CONSTS(81) - STATES(42)) - STATES(38))/2.00000
      ALGBRC(60) = ( CONSTS(83)*ALGBRC(56))/(CONSTS(84)+ALGBRC(56))
      ALGBRC(63) =  CONSTS(85)*STATES(3)*STATES(42)
      RATES(42) = ((- ALGBRC(54) - ALGBRC(57))+ 2.00000*ALGBRC(60)) - ALGBRC(63)
      RATES(38) = (((ALGBRC(54) -  2.00000*ALGBRC(53))+ALGBRC(57)) - ALGBRC(49))+ALGBRC(63)
      ALGBRC(59) =  CONSTS(86)*STATES(37)
      RATES(37) = (((- ALGBRC(57)+ALGBRC(47)) - ALGBRC(48)) - ALGBRC(59)) - ALGBRC(62)
      RATES(10) = - ALGBRC(62)
      ALGBRC(44) =  - CONSTS(4)*STATES(29)
      ALGBRC(65) =  CONSTS(87)*STATES(3)*STATES(29)
      RATES(29) = ((((ALGBRC(41)+ALGBRC(43)) - ALGBRC(38)) - ALGBRC(36))+ALGBRC(35)+ALGBRC(44)) - ALGBRC(65)
      ALGBRC(27) =  - CONSTS(4)*STATES(17)
      ALGBRC(66) =  CONSTS(88)*STATES(3)*STATES(17)
      RATES(17) = ((((((ALGBRC(35)+ALGBRC(34)) - ALGBRC(28)) - ALGBRC(24))+ALGBRC(26)) - ALGBRC(21))+ALGBRC(22)+ALGBRC(27)) - ALGBRC(66)
      RATES(40) = - ALGBRC(51)+ALGBRC(63)+ALGBRC(4)+ALGBRC(65)+ALGBRC(66)
      RATES(3) = (((ALGBRC(51) - ALGBRC(63)) - ALGBRC(4)) - ALGBRC(65)) - ALGBRC(66)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) =  CONSTS(2)*STATES(1)*STATES(4) -  CONSTS(3)*STATES(2)
      ALGBRC(1) =  CONSTS(1)*STATES(2)
      ALGBRC(5) = CONSTS(6) -  CONSTS(4)*STATES(5)
      ALGBRC(6) =  CONSTS(7)*STATES(7)*STATES(5) -  CONSTS(8)*STATES(6)
      ALGBRC(8) =  CONSTS(9)*STATES(6) ** CONSTS(11) -  CONSTS(10)*STATES(8)
      ALGBRC(9) = CONSTS(12) - STATES(10)
      ALGBRC(13) =  CONSTS(15)*STATES(11)*STATES(9)
      ALGBRC(7) =  - CONSTS(4)*STATES(7)
      ALGBRC(11) = CONSTS(14) -  CONSTS(4)*STATES(9)
      ALGBRC(15) =  CONSTS(17)*STATES(14)
      ALGBRC(14) =  CONSTS(16)*STATES(12)*STATES(13)
      ALGBRC(12) =  - CONSTS(4)*STATES(12)
      ALGBRC(18) =  CONSTS(18)*STATES(15)*STATES(13)
      ALGBRC(10) =  - CONSTS(4)*STATES(11)
      ALGBRC(17) =  - CONSTS(4)*STATES(15)
      ALGBRC(20) =  CONSTS(19)*STATES(16)*STATES(13)
      ALGBRC(19) = CONSTS(92) -  CONSTS(23)*STATES(13)
      ALGBRC(21) =  CONSTS(24)*STATES(17)*STATES(16) -  CONSTS(25)*STATES(18)
      ALGBRC(23) = CONSTS(93) -  CONSTS(27)*STATES(16)
      ALGBRC(22) =  CONSTS(28)*STATES(18)
      ALGBRC(26) =  CONSTS(31)*STATES(19)
      ALGBRC(16) =  - CONSTS(4)*STATES(14)
      ALGBRC(24) =  CONSTS(29)*STATES(17)*STATES(4) -  CONSTS(30)*STATES(19)
      ALGBRC(25) = CONSTS(32) -  CONSTS(4)*STATES(4)
      ALGBRC(28) =  CONSTS(33)*STATES(17)*STATES(20) -  CONSTS(34)*STATES(21)
      ALGBRC(30) =  CONSTS(35)*STATES(8)*STATES(22) -  CONSTS(36)*STATES(23)
      ALGBRC(31) =  CONSTS(37)*STATES(23)*STATES(22) -  CONSTS(38)*STATES(24)
      ALGBRC(29) = CONSTS(39) -  CONSTS(4)*STATES(22)
      ALGBRC(32) =  CONSTS(40)*STATES(24)
      ALGBRC(34) =  CONSTS(43)*STATES(25)
      ALGBRC(33) =  CONSTS(41)*STATES(26)*STATES(27) -  CONSTS(42)*STATES(25)
      ALGBRC(35) =  CONSTS(44)*STATES(28)
      ALGBRC(36) =  CONSTS(45)*STATES(29)*STATES(27) -  CONSTS(46)*STATES(28)
      ALGBRC(37) = CONSTS(47) -  CONSTS(4)*STATES(27)
      ALGBRC(38) =  CONSTS(49)*STATES(29)*STATES(20) -  CONSTS(50)*STATES(30)
      ALGBRC(39) =  CONSTS(51)*STATES(33)*STATES(20) -  CONSTS(52)*STATES(31)
      ALGBRC(40) =  CONSTS(53)*STATES(26)*STATES(20) -  CONSTS(54)*STATES(32)
      ALGBRC(42) = CONSTS(48) -  CONSTS(4)*STATES(20)
      ALGBRC(41) =  CONSTS(55)*STATES(26) -  CONSTS(56)*STATES(33)*STATES(29)
      ALGBRC(43) =  CONSTS(57)*STATES(33) -  CONSTS(58)*STATES(8)*STATES(29)
      ALGBRC(46) =  CONSTS(62)*STATES(34) ** 2.00000*CONSTS(59)
      ALGBRC(47) =  CONSTS(63)*STATES(36)*STATES(34)
      ALGBRC(48) =  CONSTS(64)*STATES(37)
      ALGBRC(50) =  CONSTS(66)*STATES(39)*STATES(34)
      ALGBRC(45) =  CONSTS(61)*STATES(34)*STATES(35)
      ALGBRC(49) =  CONSTS(65)*STATES(38)
      ALGBRC(51) =  CONSTS(67)*STATES(40)*STATES(34)
      ALGBRC(52) =  CONSTS(70)*CONSTS(68)*STATES(35)
      ALGBRC(53) =  CONSTS(71)*STATES(38) ** 2.00000*STATES(35)
      ALGBRC(2) =  - CONSTS(4)*STATES(1)
      ALGBRC(4) =  CONSTS(5)*STATES(3)*STATES(1)
      ALGBRC(62) =  CONSTS(89)*STATES(37)*STATES(1)
      ALGBRC(54) =  CONSTS(75)*STATES(41)*STATES(42)
      ALGBRC(55) =  CONSTS(76)*STATES(41)*CONSTS(72)
      ALGBRC(58) =  CONSTS(77)*STATES(41)*CONSTS(73)
      ALGBRC(61) =  CONSTS(78)*STATES(41)*CONSTS(74)
      ALGBRC(64) =  CONSTS(79)*STATES(41)*STATES(10)
      ALGBRC(57) =  CONSTS(82)*STATES(37)*STATES(42)
      ALGBRC(56) = ((CONSTS(81) - STATES(42)) - STATES(38))/2.00000
      ALGBRC(60) = ( CONSTS(83)*ALGBRC(56))/(CONSTS(84)+ALGBRC(56))
      ALGBRC(63) =  CONSTS(85)*STATES(3)*STATES(42)
      ALGBRC(59) =  CONSTS(86)*STATES(37)
      ALGBRC(44) =  - CONSTS(4)*STATES(29)
      ALGBRC(65) =  CONSTS(87)*STATES(3)*STATES(29)
      ALGBRC(27) =  - CONSTS(4)*STATES(17)
      ALGBRC(66) =  CONSTS(88)*STATES(3)*STATES(17)
      RETURN
      END