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