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