C
C There are a total of 2 entries in the algebraic variable array.
C There are a total of 5 entries in each of the rate and state variable arrays.
C There are a total of 20 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is a_D in component parameters (first_order_rate_constant).
C CONSTS(2) is a_E in component parameters (first_order_rate_constant).
C CONSTS(3) is a_X in component parameters (first_order_rate_constant).
C CONSTS(4) is k in component parameters (dimensionless).
C CONSTS(5) is q_D in component parameters (dimensionless).
C CONSTS(6) is q_E in component parameters (dimensionless).
C CONSTS(7) is q_X in component parameters (dimensionless).
C CONSTS(8) is f in component parameters (first_order_rate_constant).
C CONSTS(9) is g in component parameters (first_order_rate_constant).
C CONSTS(10) is p_S in component parameters (first_order_rate_constant).
C CONSTS(11) is p_D in component parameters (first_order_rate_constant).
C CONSTS(12) is p_E in component parameters (first_order_rate_constant).
C CONSTS(13) is p_X in component parameters (first_order_rate_constant).
C CONSTS(14) is d_D in component parameters (first_order_rate_constant).
C CONSTS(15) is d_E in component parameters (first_order_rate_constant).
C CONSTS(16) is d_X in component parameters (first_order_rate_constant).
C CONSTS(17) is a_f in component parameters (first_order_rate_constant).
C CONSTS(18) is R_T in component parameters (dimensionless).
C CONSTS(19) is theta in component parameters (dimensionless).
C CONSTS(20) is GF in component parameters (dimensionless).
C STATES(1) is D in component D (dimensionless).
C STATES(2) is E in component E (dimensionless).
C STATES(3) is R_S in component R_S (dimensionless).
C STATES(4) is X in component X (dimensionless).
C STATES(5) is R in component R (dimensionless).
C ALGBRC(1) is unpho_RB in component unpho_RB (dimensionless).
C ALGBRC(2) is free_E2F in component free_E2F (dimensionless).
C RATES(1) is d/dt D in component D (dimensionless).
C RATES(2) is d/dt E in component E (dimensionless).
C RATES(5) is d/dt R in component R (dimensionless).
C RATES(3) is d/dt R_S in component R_S (dimensionless).
C RATES(4) is d/dt X in component X (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.4
      CONSTS(2) = 0.16
      CONSTS(3) = 0.08
      CONSTS(4) = 0.05
      CONSTS(5) = 0.6
      CONSTS(6) = 0.6
      CONSTS(7) = 0.8
      CONSTS(8) = 0.2
      CONSTS(9) = 0.528
      CONSTS(10) = 0.6
      CONSTS(11) = 0.48
      CONSTS(12) = 0.096
      CONSTS(13) = 0.48
      CONSTS(14) = 0.4
      CONSTS(15) = 0.2
      CONSTS(16) = 1.04
      CONSTS(17) = 0.9
      CONSTS(18) = 2.5
      CONSTS(19) = 1.5
      CONSTS(20) = 6.3
      STATES(1) = 0.1
      STATES(2) = 0.6
      STATES(3) = 1
      STATES(4) = 0.7
      STATES(5) = 0.5
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ( CONSTS(1)*CONSTS(4)*CONSTS(20))/(1.00000+ CONSTS(4)*CONSTS(20)) -  CONSTS(14)*STATES(2)*STATES(1)
      RATES(2) =  CONSTS(2)*(1.00000+ CONSTS(17)*(CONSTS(19) - STATES(3)))*1.00000 -  CONSTS(15)*STATES(4)*STATES(2)
      RATES(5) = ( CONSTS(13)*STATES(4)*((CONSTS(18) - STATES(3)) - STATES(5)))/(CONSTS(7)+((CONSTS(18) - STATES(3)) - STATES(5))+STATES(4)) -  CONSTS(10)*(CONSTS(19) - STATES(3))*STATES(5)
      RATES(3) = ( CONSTS(10)*(CONSTS(19) - STATES(3))*STATES(5) - ( CONSTS(11)*STATES(3)*STATES(1))/(CONSTS(5)+STATES(3)+STATES(1))) - ( CONSTS(12)*STATES(3)*STATES(2))/(CONSTS(6)+STATES(3)+STATES(2))
      RATES(4) = ( CONSTS(3)*STATES(2)+ CONSTS(8)*(CONSTS(19) - STATES(3))+ CONSTS(9)*STATES(4) ** 2.00000*STATES(2)) -  CONSTS(16)*STATES(4)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = STATES(5)+STATES(3)
      ALGBRC(2) = CONSTS(19) - STATES(3)
      RETURN
      END