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 30 entries in the algebraic variable array.
C There are a total of 14 entries in each of the rate and state variable arrays.
C There are a total of 93 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is CycA in component CycA (dimensionless).
C CONSTS(1) is k_assa in component kinetic_parameters (first_order_rate_constant).
C CONSTS(2) is k_dissa in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(21) is V_di in component V_di (first_order_rate_constant).
C ALGBRC(22) is V_sa in component V_sa (first_order_rate_constant).
C ALGBRC(15) is V_da in component V_da (first_order_rate_constant).
C STATES(2) is Tri_A in component Tri_A (dimensionless).
C STATES(3) is CKI in component CKI (dimensionless).
C STATES(4) is CycB in component CycB (dimensionless).
C CONSTS(3) is k_assb in component kinetic_parameters (first_order_rate_constant).
C CONSTS(4) is k_dissb in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(25) is V_sb in component V_sb (first_order_rate_constant).
C ALGBRC(26) is V_db in component V_db (first_order_rate_constant).
C ALGBRC(30) is V_25 in component V_25 (first_order_rate_constant).
C ALGBRC(28) is V_wee in component V_wee (first_order_rate_constant).
C STATES(5) is pB in component pB (dimensionless).
C STATES(6) is BCKI in component BCKI (dimensionless).
C STATES(7) is CycE in component CycE (dimensionless).
C CONSTS(5) is k_disse in component kinetic_parameters (first_order_rate_constant).
C CONSTS(6) is k_asse in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(24) is V_se in component V_se (first_order_rate_constant).
C ALGBRC(14) is V_de in component V_de (first_order_rate_constant).
C STATES(8) is Tri_E in component Tri_E (dimensionless).
C ALGBRC(2) is CycD in component CycD (dimensionless).
C CONSTS(7) is CycD_0 in component CycD (dimensionless).
C STATES(9) is mass in component mass (dimensionless).
C STATES(10) is Cdc20_A in component Cdc20_A (dimensionless).
C CONSTS(8) is k_a20 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(9) is k_i20 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(10) is k_d20 in component kinetic_parameters (first_order_rate_constant).
C STATES(11) is APCP in component APCP (dimensionless).
C STATES(12) is Cdc20_i in component Cdc20_i (dimensionless).
C CONSTS(11) is J_a20 in component Cdc20_A (dimensionless).
C CONSTS(12) is J_i20 in component Cdc20_A (dimensionless).
C CONSTS(13) is k_s20p in component kinetic_parameters (first_order_rate_constant).
C CONSTS(14) is k_s20pp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(15) is n20 in component Cdc20_i (dimensionless).
C CONSTS(16) is J_20 in component Cdc20_i (dimensionless).
C CONSTS(17) is J_a20 in component Cdc20_i (dimensionless).
C CONSTS(18) is J_i20 in component Cdc20_i (dimensionless).
C ALGBRC(1) is APC in component APC (dimensionless).
C CONSTS(19) is APC_T in component APC (dimensionless).
C CONSTS(20) is k_aie in component kinetic_parameters (first_order_rate_constant).
C CONSTS(21) is k_iie in component kinetic_parameters (first_order_rate_constant).
C CONSTS(22) is J_aie in component APCP (dimensionless).
C CONSTS(23) is J_iie in component APCP (dimensionless).
C STATES(13) is pBCKI in component pBCKI (dimensionless).
C ALGBRC(19) is V_si in component V_si (first_order_rate_constant).
C ALGBRC(10) is Cdh1_i in component Cdh1_i (dimensionless).
C CONSTS(24) is Cdh1_T in component Cdh1_i (dimensionless).
C STATES(14) is Cdh1 in component Cdh1 (dimensionless).
C ALGBRC(18) is V_ah1 in component V_ah1 (first_order_rate_constant).
C ALGBRC(20) is V_ih1 in component V_ih1 (first_order_rate_constant).
C CONSTS(25) is J_ah1 in component Cdh1 (dimensionless).
C CONSTS(26) is J_ih1 in component Cdh1 (dimensionless).
C CONSTS(27) is mu in component mass (first_order_rate_constant).
C CONSTS(28) is maxmass in component mass (dimensionless).
C ALGBRC(11) is V_atf in component V_atf (first_order_rate_constant).
C CONSTS(29) is k_atfp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(30) is k_atfapp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(31) is k_atfepp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(32) is k_atfdpp in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(12) is V_itf in component V_itf (first_order_rate_constant).
C CONSTS(33) is k_itfp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(34) is k_itfapp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(35) is k_itfbpp in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(13) is TF_E in component TF_E (dimensionless).
C CONSTS(36) is J_itf in component TF_E (dimensionless).
C CONSTS(37) is J_atf in component TF_E (dimensionless).
C CONSTS(38) is k_dep in component kinetic_parameters (first_order_rate_constant).
C CONSTS(39) is k_deepp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(40) is k_deapp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(41) is k_debpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(42) is k_dap in component kinetic_parameters (first_order_rate_constant).
C CONSTS(43) is k_dapp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(44) is k_dappp in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(17) is TF_I in component TF_I (dimensionless).
C CONSTS(45) is k_afi in component kinetic_parameters (first_order_rate_constant).
C CONSTS(46) is k_ifip in component kinetic_parameters (first_order_rate_constant).
C CONSTS(47) is k_ifibpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(48) is J_ifi in component TF_I (dimensionless).
C CONSTS(49) is J_afi in component TF_I (dimensionless).
C ALGBRC(16) is Cdc14 in component TF_I (dimensionless).
C CONSTS(50) is k_ah1p in component kinetic_parameters (first_order_rate_constant).
C CONSTS(51) is k_ah1pp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(52) is k_ih1p in component kinetic_parameters (first_order_rate_constant).
C CONSTS(53) is k_ih1app in component kinetic_parameters (first_order_rate_constant).
C CONSTS(54) is k_ih1bpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(55) is k_ih1epp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(56) is k_ih1dpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(57) is k_sip in component kinetic_parameters (first_order_rate_constant).
C CONSTS(58) is k_sipp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(59) is k_dip in component kinetic_parameters (first_order_rate_constant).
C CONSTS(60) is k_diapp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(61) is k_diepp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(62) is k_didpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(63) is k_dibpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(64) is k_14di in component kinetic_parameters (dimensionless).
C ALGBRC(3) is preMPF in component preMPF (dimensionless).
C ALGBRC(4) is CycBT in component CycBT (dimensionless).
C ALGBRC(5) is CycAT in component CycAT (dimensionless).
C ALGBRC(6) is CycET in component CycET (dimensionless).
C ALGBRC(7) is Tri_B in component Tri_B (dimensionless).
C ALGBRC(8) is CKIT in component CKIT (dimensionless).
C ALGBRC(9) is Cdc20_T in component Cdc20_T (dimensionless).
C ALGBRC(23) is TF_B in component TF_B (dimensionless).
C CONSTS(65) is k_afb in component kinetic_parameters (first_order_rate_constant).
C CONSTS(66) is k_ifb in component kinetic_parameters (first_order_rate_constant).
C CONSTS(67) is J_ifb in component TF_B (dimensionless).
C CONSTS(68) is J_afb in component TF_B (dimensionless).
C CONSTS(69) is k_sbp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(70) is k_sbpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(71) is k_sap in component kinetic_parameters (first_order_rate_constant).
C CONSTS(72) is k_sapp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(73) is k_sep in component kinetic_parameters (first_order_rate_constant).
C CONSTS(74) is k_sepp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(75) is k_dbp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(76) is k_dbhpp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(77) is k_dbcpp in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(27) is Wee1 in component Wee1 (dimensionless).
C CONSTS(78) is k_aweep in component kinetic_parameters (first_order_rate_constant).
C CONSTS(79) is k_aweepp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(80) is k_iweep in component kinetic_parameters (first_order_rate_constant).
C CONSTS(81) is k_iweepp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(82) is J_iwee in component Wee1 (dimensionless).
C CONSTS(83) is J_awee in component Wee1 (dimensionless).
C CONSTS(84) is k_weep in component kinetic_parameters (first_order_rate_constant).
C CONSTS(85) is k_weepp in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(29) is Cdc25 in component Cdc25 (dimensionless).
C CONSTS(86) is k_a25p in component kinetic_parameters (first_order_rate_constant).
C CONSTS(87) is k_a25pp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(88) is k_i25p in component kinetic_parameters (first_order_rate_constant).
C CONSTS(89) is k_i25pp in component kinetic_parameters (first_order_rate_constant).
C CONSTS(90) is J_i25 in component Cdc25 (dimensionless).
C CONSTS(91) is J_a25 in component Cdc25 (dimensionless).
C CONSTS(92) is k_25p in component kinetic_parameters (first_order_rate_constant).
C CONSTS(93) is k_25pp in component kinetic_parameters (first_order_rate_constant).
C RATES(1) is d/dt CycA in component CycA (dimensionless).
C RATES(4) is d/dt CycB in component CycB (dimensionless).
C RATES(7) is d/dt CycE in component CycE (dimensionless).
C RATES(10) is d/dt Cdc20_A in component Cdc20_A (dimensionless).
C RATES(12) is d/dt Cdc20_i in component Cdc20_i (dimensionless).
C RATES(11) is d/dt APCP in component APCP (dimensionless).
C RATES(5) is d/dt pB in component pB (dimensionless).
C RATES(6) is d/dt BCKI in component BCKI (dimensionless).
C RATES(13) is d/dt pBCKI in component pBCKI (dimensionless).
C RATES(2) is d/dt Tri_A in component Tri_A (dimensionless).
C RATES(8) is d/dt Tri_E in component Tri_E (dimensionless).
C RATES(3) is d/dt CKI in component CKI (dimensionless).
C RATES(14) is d/dt Cdh1 in component Cdh1 (dimensionless).
C RATES(9) is d/dt mass in component mass (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0.09450308233499527
CONSTS(1) = 50
CONSTS(2) = 0.06
STATES(2) = 0.3492224216461182
STATES(3) = 0.02882070094347
STATES(4) = 0.1903585940599442
CONSTS(3) = 60
CONSTS(4) = 0.05
STATES(5) = 0.01234426256269217
STATES(6) = 0.679449200630188
STATES(7) = 0.2092024385929108
CONSTS(5) = 0
CONSTS(6) = 0
STATES(8) = 0
CONSTS(7) = 0.108
STATES(9) = 1.338267803192139
STATES(10) = 0.3572723865509033
CONSTS(8) = 1
CONSTS(9) = 0.16
CONSTS(10) = 0.05
STATES(11) = 0.07591402530670166
STATES(12) = 0.7702387571334839
CONSTS(11) = 1
CONSTS(12) = 1
CONSTS(13) = 0.001
CONSTS(14) = 1
CONSTS(15) = 1
CONSTS(16) = 10
CONSTS(17) = 1
CONSTS(18) = 1
CONSTS(19) = 1
CONSTS(20) = 0.1
CONSTS(21) = 0.15
CONSTS(22) = 0.1
CONSTS(23) = 0.1
STATES(13) = 0.04795938357710838
CONSTS(24) = 1
STATES(14) = 0.7189393639564514
CONSTS(25) = 0.03
CONSTS(26) = 0.03
CONSTS(27) = 0.005776
CONSTS(28) = 10000
CONSTS(29) = 0
CONSTS(30) = 1.5
CONSTS(31) = 0.38
CONSTS(32) = 3
CONSTS(33) = 0.75
CONSTS(34) = 0
CONSTS(35) = 8
CONSTS(36) = 0.01
CONSTS(37) = 0.01
CONSTS(38) = 0.12
CONSTS(39) = 0
CONSTS(40) = 0
CONSTS(41) = 0
CONSTS(42) = 0.01
CONSTS(43) = 0.16
CONSTS(44) = 0
CONSTS(45) = 6
CONSTS(46) = 0.008
CONSTS(47) = 0.05
CONSTS(48) = 2
CONSTS(49) = 1
CONSTS(50) = 0.02
CONSTS(51) = 0.8
CONSTS(52) = 0.001
CONSTS(53) = 0.35
CONSTS(54) = 0.1
CONSTS(55) = 0.06
CONSTS(56) = 0.005
CONSTS(57) = 0.018
CONSTS(58) = 0.18
CONSTS(59) = 0.02
CONSTS(60) = 0.1
CONSTS(61) = 0.12
CONSTS(62) = 0.1
CONSTS(63) = 0.8
CONSTS(64) = 12
CONSTS(65) = 1
CONSTS(66) = 0.15
CONSTS(67) = 0.1
CONSTS(68) = 0.1
CONSTS(69) = 0.004
CONSTS(70) = 0.04
CONSTS(71) = 0.0015
CONSTS(72) = 0.01
CONSTS(73) = 0
CONSTS(74) = 0.15
CONSTS(75) = 0.003
CONSTS(76) = 0.4
CONSTS(77) = 0.15
CONSTS(78) = 0.3
CONSTS(79) = 0
CONSTS(80) = 0
CONSTS(81) = 1
CONSTS(82) = 0.05
CONSTS(83) = 0.05
CONSTS(84) = 0.02
CONSTS(85) = 0.2
CONSTS(86) = 0
CONSTS(87) = 1
CONSTS(88) = 0.3
CONSTS(89) = 0
CONSTS(90) = 0.1
CONSTS(91) = 0.1
CONSTS(92) = 0.01
CONSTS(93) = 5
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(10) = (( CONSTS(8)*STATES(11)*STATES(12))/(CONSTS(11)+STATES(12)) - ( STATES(10)*CONSTS(9))/(CONSTS(12)+STATES(10))) - STATES(10)*CONSTS(10)
RATES(12) = (((CONSTS(13)+ CONSTS(14)*STATES(4) ** CONSTS(15))/(CONSTS(16) ** CONSTS(15)+STATES(4) ** CONSTS(15))+( STATES(10)*CONSTS(9))/(CONSTS(18)+STATES(10))) - CONSTS(10)*STATES(12)) - ( CONSTS(8)*STATES(11)*STATES(12))/(CONSTS(17)+STATES(12))
RATES(9) = CONSTS(27)*STATES(9)*(1.00000 - STATES(9)/CONSTS(28))
ALGBRC(1) = (CONSTS(19) - STATES(11))/1.00000
RATES(11) = ( CONSTS(20)*STATES(4)*ALGBRC(1))/(CONSTS(22)+ALGBRC(1)) - ( CONSTS(21)*STATES(11))/(CONSTS(23)+STATES(11))
ALGBRC(10) = (CONSTS(24) - STATES(14))/1.00000
ALGBRC(16) = STATES(10)
ALGBRC(18) = CONSTS(50)+ CONSTS(51)*ALGBRC(16)
ALGBRC(2) = CONSTS(7)*STATES(9)
ALGBRC(20) = CONSTS(52)+ CONSTS(53)*STATES(1)+ CONSTS(54)*STATES(4)+ CONSTS(55)*STATES(7)+ CONSTS(56)*ALGBRC(2)
RATES(14) = ( ALGBRC(10)*ALGBRC(18))/(CONSTS(25)+ALGBRC(10)) - ( STATES(14)*ALGBRC(20))/(CONSTS(26)+STATES(14))
ALGBRC(21) = (CONSTS(59)+ CONSTS(60)*STATES(1)+ CONSTS(63)*STATES(4)+ CONSTS(61)*STATES(7)+ CONSTS(62)*ALGBRC(2))/(1.00000+ CONSTS(64)*ALGBRC(16))
ALGBRC(15) = CONSTS(42)+ (CONSTS(43)+CONSTS(44))*STATES(10)+ CONSTS(44)*STATES(12)
RATES(2) = (( CONSTS(1)*STATES(3)*STATES(1) - CONSTS(2)*STATES(2)) - ALGBRC(21)*STATES(2)) - ALGBRC(15)*STATES(2)
ALGBRC(14) = CONSTS(38)+ CONSTS(39)*STATES(7)+ CONSTS(40)*STATES(1)+ CONSTS(41)*STATES(4)
RATES(8) = (( CONSTS(6)*STATES(3)*STATES(7) - CONSTS(5)*STATES(8)) - ALGBRC(21)*STATES(8)) - ALGBRC(14)*STATES(8)
ALGBRC(11) = CONSTS(29)+ CONSTS(30)*STATES(1)+ CONSTS(31)*STATES(7)+ CONSTS(32)*ALGBRC(2)
ALGBRC(12) = CONSTS(33)+ CONSTS(34)*STATES(1)+ CONSTS(35)*STATES(4)
ALGBRC(13) = ( 2.00000*ALGBRC(11)*CONSTS(36))/((ALGBRC(12) - ALGBRC(11))+ CONSTS(37)*ALGBRC(12)+ CONSTS(36)*ALGBRC(11)+ ((ALGBRC(12) - ALGBRC(11))+ CONSTS(37)*ALGBRC(12)+ CONSTS(36)*ALGBRC(11) ** 2.00000 - 4.00000*(ALGBRC(12) - ALGBRC(11))*CONSTS(36)*ALGBRC(11)) ** (1.0 / 2))
ALGBRC(22) = STATES(9)*(CONSTS(71)+ CONSTS(72)*ALGBRC(13))
RATES(1) = (( CONSTS(2)*STATES(2)+ ALGBRC(21)*STATES(2)+ALGBRC(22)) - ALGBRC(15)*STATES(1)) - CONSTS(1)*STATES(3)*STATES(1)
ALGBRC(24) = STATES(9)*(CONSTS(73)+ CONSTS(74)*ALGBRC(13))
RATES(7) = (( CONSTS(5)*STATES(8)+ ALGBRC(21)*STATES(8)+ALGBRC(24)) - ALGBRC(14)*STATES(7)) - CONSTS(6)*STATES(3)*STATES(7)
ALGBRC(26) = CONSTS(75)+ CONSTS(76)*STATES(14)+ CONSTS(77)*STATES(10)
ALGBRC(17) = ( 2.00000*CONSTS(45)*ALGBRC(16)*CONSTS(48))/(((CONSTS(46)+ CONSTS(47)*STATES(4)) - CONSTS(45)*ALGBRC(16))+ CONSTS(49)*(CONSTS(46)+ CONSTS(47)*STATES(4))+ CONSTS(48)*CONSTS(45)*ALGBRC(16)+ (((CONSTS(46)+ CONSTS(47)*STATES(4)) - CONSTS(45)*ALGBRC(16))+ CONSTS(49)*(CONSTS(46)+ CONSTS(47)*STATES(4))+ CONSTS(48)*CONSTS(45)*ALGBRC(16) ** 2.00000 - 4.00000*((CONSTS(46)+ CONSTS(47)*STATES(4)) - CONSTS(45)*ALGBRC(16))*CONSTS(48)*CONSTS(45)*ALGBRC(16)) ** (1.0 / 2))
ALGBRC(19) = CONSTS(57)+ CONSTS(58)*ALGBRC(17)
RATES(3) = ((((((( - CONSTS(3)*STATES(4)*STATES(3)+ CONSTS(4)*STATES(6)) - CONSTS(3)*STATES(5)*STATES(3))+ CONSTS(4)*STATES(13)+ ALGBRC(26)*STATES(6)+ ALGBRC(26)*STATES(13)+ALGBRC(19)) - ALGBRC(21)*STATES(3)) - CONSTS(1)*STATES(3)*STATES(1))+ CONSTS(2)*STATES(2)+ ALGBRC(15)*STATES(2)) - CONSTS(6)*STATES(3)*STATES(7))+ CONSTS(5)*STATES(8)+ ALGBRC(14)*STATES(8)
ALGBRC(23) = ( 2.00000*CONSTS(65)*STATES(4)*CONSTS(67))/((CONSTS(66) - CONSTS(65)*STATES(4))+ CONSTS(68)*CONSTS(66)+ CONSTS(67)*CONSTS(65)*STATES(4)+ ((CONSTS(66) - CONSTS(65)*STATES(4))+ CONSTS(68)*CONSTS(66)+ CONSTS(67)*CONSTS(65)*STATES(4) ** 2.00000 - 4.00000*(CONSTS(66) - CONSTS(65)*STATES(4))*CONSTS(65)*STATES(4)*CONSTS(67)) ** (1.0 / 2))
ALGBRC(25) = STATES(9)*(CONSTS(69)+ CONSTS(70)*ALGBRC(23))
ALGBRC(29) = ( 2.00000*(CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90))/(((CONSTS(88)+ CONSTS(89)*ALGBRC(16)) - (CONSTS(86)+ CONSTS(87)*STATES(4)))+ CONSTS(91)*(CONSTS(88)+ CONSTS(89)*ALGBRC(16))+ (CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90)+ (((CONSTS(88)+ CONSTS(89)*ALGBRC(16)) - (CONSTS(86)+ CONSTS(87)*STATES(4)))+ CONSTS(91)*(CONSTS(88)+ CONSTS(89)*ALGBRC(16))+ (CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90) ** 2.00000 - 4.00000*((CONSTS(88)+ CONSTS(89)*ALGBRC(16)) - (CONSTS(86)+ CONSTS(87)*STATES(4)))*(CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90)) ** (1.0 / 2))
ALGBRC(30) = CONSTS(92)+ CONSTS(93)*ALGBRC(29)
ALGBRC(27) = ( 2.00000*(CONSTS(78)+ CONSTS(79)*ALGBRC(16))*CONSTS(82))/(((CONSTS(80)+ CONSTS(81)*STATES(4)) - (CONSTS(78)+ CONSTS(79)*ALGBRC(16)))+ CONSTS(83)*(CONSTS(80)+ CONSTS(81)*STATES(4))+ CONSTS(82)*(CONSTS(78)+ CONSTS(79)*ALGBRC(16))+ (((CONSTS(80)+ CONSTS(81)*STATES(4)) - (CONSTS(78)+ CONSTS(79)*ALGBRC(16)))+ CONSTS(83)*(CONSTS(80)+ CONSTS(81)*STATES(4))+ CONSTS(82)*(CONSTS(78)+ CONSTS(79)*ALGBRC(16)) ** 2.00000 - 4.00000*((CONSTS(80)+ CONSTS(81)*STATES(4)) - (CONSTS(78)+ CONSTS(79)*ALGBRC(16)))*(CONSTS(78)+ CONSTS(79)*ALGBRC(16))*CONSTS(82)) ** (1.0 / 2))
ALGBRC(28) = CONSTS(84)+ CONSTS(85)*ALGBRC(27)
RATES(4) = ((((ALGBRC(25) - ALGBRC(26)*STATES(4))+ ALGBRC(30)*STATES(5)) - ALGBRC(28)*STATES(4)) - CONSTS(3)*STATES(4)*STATES(3))+ CONSTS(4)*STATES(6)+ ALGBRC(21)*STATES(6)
RATES(5) = ((( ALGBRC(28)*STATES(4) - ALGBRC(26)*STATES(5)) - CONSTS(3)*STATES(5)*STATES(3))+ CONSTS(4)*STATES(13)+ ALGBRC(21)*STATES(13)) - ALGBRC(30)*STATES(5)
RATES(6) = (((( CONSTS(3)*STATES(4)*STATES(3) - CONSTS(4)*STATES(6))+ ALGBRC(30)*STATES(13)) - ALGBRC(28)*STATES(6)) - ALGBRC(26)*STATES(6)) - ALGBRC(21)*STATES(6)
RATES(13) = (((( CONSTS(3)*STATES(5)*STATES(3) - CONSTS(4)*STATES(13)) - ALGBRC(30)*STATES(13))+ ALGBRC(28)*STATES(6)) - ALGBRC(26)*STATES(13)) - ALGBRC(21)*STATES(13)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = (CONSTS(19) - STATES(11))/1.00000
ALGBRC(10) = (CONSTS(24) - STATES(14))/1.00000
ALGBRC(16) = STATES(10)
ALGBRC(18) = CONSTS(50)+ CONSTS(51)*ALGBRC(16)
ALGBRC(2) = CONSTS(7)*STATES(9)
ALGBRC(20) = CONSTS(52)+ CONSTS(53)*STATES(1)+ CONSTS(54)*STATES(4)+ CONSTS(55)*STATES(7)+ CONSTS(56)*ALGBRC(2)
ALGBRC(21) = (CONSTS(59)+ CONSTS(60)*STATES(1)+ CONSTS(63)*STATES(4)+ CONSTS(61)*STATES(7)+ CONSTS(62)*ALGBRC(2))/(1.00000+ CONSTS(64)*ALGBRC(16))
ALGBRC(15) = CONSTS(42)+ (CONSTS(43)+CONSTS(44))*STATES(10)+ CONSTS(44)*STATES(12)
ALGBRC(14) = CONSTS(38)+ CONSTS(39)*STATES(7)+ CONSTS(40)*STATES(1)+ CONSTS(41)*STATES(4)
ALGBRC(11) = CONSTS(29)+ CONSTS(30)*STATES(1)+ CONSTS(31)*STATES(7)+ CONSTS(32)*ALGBRC(2)
ALGBRC(12) = CONSTS(33)+ CONSTS(34)*STATES(1)+ CONSTS(35)*STATES(4)
ALGBRC(13) = ( 2.00000*ALGBRC(11)*CONSTS(36))/((ALGBRC(12) - ALGBRC(11))+ CONSTS(37)*ALGBRC(12)+ CONSTS(36)*ALGBRC(11)+ ((ALGBRC(12) - ALGBRC(11))+ CONSTS(37)*ALGBRC(12)+ CONSTS(36)*ALGBRC(11) ** 2.00000 - 4.00000*(ALGBRC(12) - ALGBRC(11))*CONSTS(36)*ALGBRC(11)) ** (1.0 / 2))
ALGBRC(22) = STATES(9)*(CONSTS(71)+ CONSTS(72)*ALGBRC(13))
ALGBRC(24) = STATES(9)*(CONSTS(73)+ CONSTS(74)*ALGBRC(13))
ALGBRC(26) = CONSTS(75)+ CONSTS(76)*STATES(14)+ CONSTS(77)*STATES(10)
ALGBRC(17) = ( 2.00000*CONSTS(45)*ALGBRC(16)*CONSTS(48))/(((CONSTS(46)+ CONSTS(47)*STATES(4)) - CONSTS(45)*ALGBRC(16))+ CONSTS(49)*(CONSTS(46)+ CONSTS(47)*STATES(4))+ CONSTS(48)*CONSTS(45)*ALGBRC(16)+ (((CONSTS(46)+ CONSTS(47)*STATES(4)) - CONSTS(45)*ALGBRC(16))+ CONSTS(49)*(CONSTS(46)+ CONSTS(47)*STATES(4))+ CONSTS(48)*CONSTS(45)*ALGBRC(16) ** 2.00000 - 4.00000*((CONSTS(46)+ CONSTS(47)*STATES(4)) - CONSTS(45)*ALGBRC(16))*CONSTS(48)*CONSTS(45)*ALGBRC(16)) ** (1.0 / 2))
ALGBRC(19) = CONSTS(57)+ CONSTS(58)*ALGBRC(17)
ALGBRC(23) = ( 2.00000*CONSTS(65)*STATES(4)*CONSTS(67))/((CONSTS(66) - CONSTS(65)*STATES(4))+ CONSTS(68)*CONSTS(66)+ CONSTS(67)*CONSTS(65)*STATES(4)+ ((CONSTS(66) - CONSTS(65)*STATES(4))+ CONSTS(68)*CONSTS(66)+ CONSTS(67)*CONSTS(65)*STATES(4) ** 2.00000 - 4.00000*(CONSTS(66) - CONSTS(65)*STATES(4))*CONSTS(65)*STATES(4)*CONSTS(67)) ** (1.0 / 2))
ALGBRC(25) = STATES(9)*(CONSTS(69)+ CONSTS(70)*ALGBRC(23))
ALGBRC(29) = ( 2.00000*(CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90))/(((CONSTS(88)+ CONSTS(89)*ALGBRC(16)) - (CONSTS(86)+ CONSTS(87)*STATES(4)))+ CONSTS(91)*(CONSTS(88)+ CONSTS(89)*ALGBRC(16))+ (CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90)+ (((CONSTS(88)+ CONSTS(89)*ALGBRC(16)) - (CONSTS(86)+ CONSTS(87)*STATES(4)))+ CONSTS(91)*(CONSTS(88)+ CONSTS(89)*ALGBRC(16))+ (CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90) ** 2.00000 - 4.00000*((CONSTS(88)+ CONSTS(89)*ALGBRC(16)) - (CONSTS(86)+ CONSTS(87)*STATES(4)))*(CONSTS(86)+ CONSTS(87)*STATES(4))*CONSTS(90)) ** (1.0 / 2))
ALGBRC(30) = CONSTS(92)+ CONSTS(93)*ALGBRC(29)
ALGBRC(27) = ( 2.00000*(CONSTS(78)+ CONSTS(79)*ALGBRC(16))*CONSTS(82))/(((CONSTS(80)+ CONSTS(81)*STATES(4)) - (CONSTS(78)+ CONSTS(79)*ALGBRC(16)))+ CONSTS(83)*(CONSTS(80)+ CONSTS(81)*STATES(4))+ CONSTS(82)*(CONSTS(78)+ CONSTS(79)*ALGBRC(16))+ (((CONSTS(80)+ CONSTS(81)*STATES(4)) - (CONSTS(78)+ CONSTS(79)*ALGBRC(16)))+ CONSTS(83)*(CONSTS(80)+ CONSTS(81)*STATES(4))+ CONSTS(82)*(CONSTS(78)+ CONSTS(79)*ALGBRC(16)) ** 2.00000 - 4.00000*((CONSTS(80)+ CONSTS(81)*STATES(4)) - (CONSTS(78)+ CONSTS(79)*ALGBRC(16)))*(CONSTS(78)+ CONSTS(79)*ALGBRC(16))*CONSTS(82)) ** (1.0 / 2))
ALGBRC(28) = CONSTS(84)+ CONSTS(85)*ALGBRC(27)
ALGBRC(3) = STATES(5)+STATES(6)
ALGBRC(4) = STATES(4)+STATES(5)+STATES(6)+STATES(13)
ALGBRC(5) = STATES(1)+STATES(2)
ALGBRC(6) = STATES(7)+STATES(8)
ALGBRC(7) = STATES(6)+STATES(13)
ALGBRC(8) = STATES(3)+STATES(6)+STATES(13)+STATES(2)+STATES(8)
ALGBRC(9) = STATES(12)+STATES(10)
RETURN
END