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 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