C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 6 entries in each of the rate and state variable arrays.
C There are a total of 9 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is x_1 in component x_1 (picomolar).
C CONSTS(1) is k_t in component model_parameters (per_minute).
C CONSTS(2) is B_max in component model_parameters (picomolar).
C CONSTS(3) is k_on in component model_parameters (per_picomolar_per_minute).
C STATES(2) is x_2 in component x_2 (picomolar).
C CONSTS(4) is k_off in component model_parameters (per_minute).
C STATES(3) is x_3 in component x_3 (picomolar).
C CONSTS(5) is k_ex in component model_parameters (per_minute).
C STATES(4) is x_4 in component x_4 (picomolar).
C CONSTS(6) is k_mob in component model_parameters (per_minute).
C CONSTS(7) is k_e in component model_parameters (per_minute).
C CONSTS(8) is k_di in component model_parameters (per_minute).
C CONSTS(9) is k_de in component model_parameters (per_minute).
C STATES(5) is x_5 in component x_5 (picomolar).
C STATES(6) is x_6 in component x_6 (picomolar).
C RATES(1) is d/dt x_1 in component x_1 (picomolar).
C RATES(2) is d/dt x_2 in component x_2 (picomolar).
C RATES(3) is d/dt x_3 in component x_3 (picomolar).
C RATES(4) is d/dt x_4 in component x_4 (picomolar).
C RATES(5) is d/dt x_5 in component x_5 (picomolar).
C RATES(6) is d/dt x_6 in component x_6 (picomolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 516
      CONSTS(1) = 0.03294
      CONSTS(2) = 129
      CONSTS(3) = 0.10496e-3
      STATES(2) = 2010.19
      CONSTS(4) = 0.01721
      STATES(3) = 0
      CONSTS(5) = 0.00994
      STATES(4) = 0
      CONSTS(6) = 0.4178
      CONSTS(7) = 0.07483
      CONSTS(8) = 0.003179
      CONSTS(9) = 0.0164
      STATES(5) = 0
      STATES(6) = 0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ((( CONSTS(1)*CONSTS(2)+ CONSTS(6)*STATES(3)) -  CONSTS(1)*STATES(1)) -  CONSTS(3)*STATES(1)*STATES(2))+ CONSTS(4)*STATES(3)+ CONSTS(5)*STATES(4)
      RATES(2) =  - CONSTS(3)*STATES(1)*STATES(2)+ CONSTS(4)*STATES(3)+ CONSTS(5)*STATES(4)
      RATES(3) = ( CONSTS(3)*STATES(1)*STATES(2) -  CONSTS(4)*STATES(3)) -  CONSTS(7)*STATES(3)
      RATES(4) = (( CONSTS(7)*STATES(3) -  CONSTS(5)*STATES(4)) -  CONSTS(8)*STATES(4)) -  CONSTS(9)*STATES(4)
      RATES(5) =  CONSTS(8)*STATES(4)
      RATES(6) =  CONSTS(9)*STATES(4)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END