C
C There are a total of 2 entries in the algebraic variable array.
C There are a total of 3 entries in each of the rate and state variable arrays.
C There are a total of 13 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is C in component C (micromolar).
C CONSTS(1) is v_i in component C (micromolar_per_minute).
C CONSTS(2) is v_d in component C (micromolar_per_minute).
C STATES(2) is X in component X (micromolar).
C CONSTS(3) is K_d in component C (micromolar).
C CONSTS(4) is k_d in component C (per_minute).
C STATES(3) is M in component M (micromolar).
C ALGBRC(1) is V_1 in component V_1 (micromolar_per_minute).
C CONSTS(5) is K_1 in component M (per_minute).
C CONSTS(6) is V_2 in component M (per_minute).
C CONSTS(7) is K_2 in component M (per_minute).
C ALGBRC(2) is V_3 in component V_3 (micromolar_per_minute).
C CONSTS(8) is K_3 in component X (per_minute).
C CONSTS(9) is V_4 in component X (per_minute).
C CONSTS(10) is K_4 in component X (per_minute).
C CONSTS(11) is K_c in component V_1 (micromolar).
C CONSTS(12) is V_M1 in component V_1 (per_minute).
C CONSTS(13) is V_M3 in component V_3 (per_minute).
C RATES(1) is d/dt C in component C (micromolar).
C RATES(3) is d/dt M in component M (micromolar).
C RATES(2) is d/dt X in component X (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.1
      CONSTS(1) = 0.05
      CONSTS(2) = 0.25
      STATES(2) = 0.1
      CONSTS(3) = 0.02
      CONSTS(4) = 0.01
      STATES(3) = 0.1
      CONSTS(5) = 0.01
      CONSTS(6) = 1.5
      CONSTS(7) = 0.01
      CONSTS(8) = 0.01
      CONSTS(9) = 0.5
      CONSTS(10) = 0.01
      CONSTS(11) = 0.5
      CONSTS(12) = 3
      CONSTS(13) = 1
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = (CONSTS(1) -  CONSTS(2)*STATES(2)*(STATES(1)/(CONSTS(3)+STATES(1)))) -  CONSTS(4)*STATES(1)
      ALGBRC(1) =  (STATES(1)/(CONSTS(11)+STATES(1)))*CONSTS(12)
      RATES(3) =  ALGBRC(1)*((1.00000 - STATES(3))/(CONSTS(5)+(1.00000 - STATES(3)))) -  CONSTS(6)*(STATES(3)/(CONSTS(7)+STATES(3)))
      ALGBRC(2) =  STATES(3)*CONSTS(13)
      RATES(2) =  ALGBRC(2)*((1.00000 - STATES(2))/(CONSTS(8)+(1.00000 - STATES(2)))) -  CONSTS(9)*(STATES(2)/(CONSTS(10)+STATES(2)))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  (STATES(1)/(CONSTS(11)+STATES(1)))*CONSTS(12)
      ALGBRC(2) =  STATES(3)*CONSTS(13)
      RETURN
      END