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.009940445423126221
      CONSTS(1) = 25
      CONSTS(2) = 1
      STATES(2) = 0.01715379953384399
      STATES(3) = 0.2954076826572418
      STATES(4) = 0.1668413728475571
      CONSTS(3) = 0
      CONSTS(4) = 0
      STATES(5) = 0.009814870543777943
      STATES(6) = 0
      STATES(7) = 0.07760512828826904
      CONSTS(5) = 1
      CONSTS(6) = 50
      STATES(8) = 0.3117263317108154
      CONSTS(7) = 0.05
      STATES(9) = 1.174216866493225
      STATES(10) = 0.6605867147445679
      CONSTS(8) = 0.5
      CONSTS(9) = 0.25
      CONSTS(10) = 0.15
      STATES(11) = 0.6716265678405762
      STATES(12) = 0.01855352707207203
      CONSTS(11) = 0.005
      CONSTS(12) = 0.005
      CONSTS(13) = 0
      CONSTS(14) = 0.15
      CONSTS(15) = 1
      CONSTS(16) = 1
      CONSTS(17) = 0.005
      CONSTS(18) = 0.005
      CONSTS(19) = 1
      CONSTS(20) = 0.07
      CONSTS(21) = 0.18
      CONSTS(22) = 0.01
      CONSTS(23) = 0.01
      STATES(13) = 0
      CONSTS(24) = 1
      STATES(14) = 0.9992357492446899
      CONSTS(25) = 0.01
      CONSTS(26) = 0.01
      CONSTS(27) = 0.004951
      CONSTS(28) = 10000
      CONSTS(29) = 0
      CONSTS(30) = 0.2
      CONSTS(31) = 0.5
      CONSTS(32) = 3
      CONSTS(33) = 0.25
      CONSTS(34) = 0.1
      CONSTS(35) = 0.1
      CONSTS(36) = 0.01
      CONSTS(37) = 0.01
      CONSTS(38) = 0.01
      CONSTS(39) = 0.1
      CONSTS(40) = 0.5
      CONSTS(41) = 0.5
      CONSTS(42) = 0.02
      CONSTS(43) = 2
      CONSTS(44) = 0
      CONSTS(45) = 88
      CONSTS(46) = 88
      CONSTS(47) = 88
      CONSTS(48) = 88
      CONSTS(49) = 88
      CONSTS(50) = 0.18
      CONSTS(51) = 3.5
      CONSTS(52) = 0
      CONSTS(53) = 0.2
      CONSTS(54) = 1
      CONSTS(55) = 0.1
      CONSTS(56) = 0
      CONSTS(57) = 1.8
      CONSTS(58) = 0
      CONSTS(59) = 0.8
      CONSTS(60) = 5
      CONSTS(61) = 5
      CONSTS(62) = 0
      CONSTS(63) = 5
      CONSTS(64) = 0
      CONSTS(65) = 1
      CONSTS(66) = 0.1
      CONSTS(67) = 0.1
      CONSTS(68) = 0.1
      CONSTS(69) = 0.01
      CONSTS(70) = 0.03
      CONSTS(71) = 0
      CONSTS(72) = 0.025
      CONSTS(73) = 0.01
      CONSTS(74) = 0.18
      CONSTS(75) = 0.005
      CONSTS(76) = 2
      CONSTS(77) = 0.1
      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