C
C There are a total of 0 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 36 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is x1 in component x1 (microg_l).
C CONSTS(1) is lambda_1 in component model_parameters (first_order_rate_constant).
C STATES(2) is x3 in component x3 (microg_l).
C CONSTS(2) is a1 in component model_parameters (flux).
C CONSTS(3) is a2 in component model_parameters (flux).
C CONSTS(4) is a3 in component model_parameters (first_order_rate_constant).
C CONSTS(5) is a4 in component model_parameters (second_order_rate_constant).
C CONSTS(6) is a5 in component model_parameters (per_microg_l).
C CONSTS(7) is a6 in component model_parameters (per_microg_l2).
C CONSTS(8) is a7 in component model_parameters (per_microg_l).
C CONSTS(9) is a8 in component model_parameters (per_microg_l2).
C STATES(3) is x2 in component x2 (microg_l).
C CONSTS(10) is lambda_2 in component model_parameters (first_order_rate_constant).
C CONSTS(11) is a9 in component model_parameters (flux).
C CONSTS(12) is a10 in component model_parameters (first_order_rate_constant).
C CONSTS(13) is a11 in component model_parameters (second_order_rate_constant).
C CONSTS(14) is a12 in component model_parameters (per_microg_l).
C CONSTS(15) is a13 in component model_parameters (per_microg_l2).
C CONSTS(16) is a14 in component model_parameters (per_microg_l).
C CONSTS(17) is a15 in component model_parameters (per_microg_l2).
C CONSTS(34) is lambda_3_ in component model_parameters (first_order_rate_constant).
C STATES(4) is x4 in component x4 (microg_l).
C STATES(5) is x5 in component x5 (microg_l).
C CONSTS(18) is a16 in component model_parameters (flux).
C CONSTS(19) is a17 in component model_parameters (first_order_rate_constant).
C CONSTS(20) is a18 in component model_parameters (second_order_rate_constant).
C CONSTS(21) is a19 in component model_parameters (first_order_rate_constant).
C CONSTS(22) is a20 in component model_parameters (second_order_rate_constant).
C CONSTS(23) is a21 in component model_parameters (per_microg_l).
C CONSTS(24) is a22 in component model_parameters (per_microg_l2).
C CONSTS(25) is a23 in component model_parameters (per_microg_l).
C CONSTS(26) is a24 in component model_parameters (per_microg_l2).
C CONSTS(27) is a25 in component model_parameters (first_order_rate_constant).
C CONSTS(28) is a26 in component model_parameters (first_order_rate_constant).
C CONSTS(35) is lambda_4_ in component model_parameters (first_order_rate_constant).
C CONSTS(29) is a27 in component model_parameters (first_order_rate_constant).
C CONSTS(36) is lambda_5_ in component model_parameters (first_order_rate_constant).
C CONSTS(30) is a28 in component model_parameters (first_order_rate_constant).
C CONSTS(31) is lambda_3 in component model_parameters (first_order_rate_constant).
C CONSTS(32) is lambda_4 in component model_parameters (first_order_rate_constant).
C CONSTS(33) is lambda_5 in component model_parameters (first_order_rate_constant).
C RATES(1) is d/dt x1 in component x1 (microg_l).
C RATES(3) is d/dt x2 in component x2 (microg_l).
C RATES(2) is d/dt x3 in component x3 (microg_l).
C RATES(4) is d/dt x4 in component x4 (microg_l).
C RATES(5) is d/dt x5 in component x5 (microg_l).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.01067
      CONSTS(1) = 0.059
      STATES(2) = 6.51
      CONSTS(2) = 0.000017
      CONSTS(3) = 0.0023
      CONSTS(4) = 0.6
      CONSTS(5) = 45
      CONSTS(6) = 36
      CONSTS(7) = 216
      CONSTS(8) = 0.28
      CONSTS(9) = 0.36
      STATES(3) = 0.04665
      CONSTS(10) = 0.028
      CONSTS(11) = 0.0003
      CONSTS(12) = 0.18
      CONSTS(13) = 150
      CONSTS(14) = 18
      CONSTS(15) = 460
      CONSTS(16) = 0.46
      CONSTS(17) = 0.1
      STATES(4) = 60.61
      STATES(5) = 12.61
      CONSTS(18) = 0.04
      CONSTS(19) = 150
      CONSTS(20) = 3800
      CONSTS(21) = 57
      CONSTS(22) = 2600
      CONSTS(23) = 200
      CONSTS(24) = 9400
      CONSTS(25) = 10
      CONSTS(26) = 320
      CONSTS(27) = 0.04
      CONSTS(28) = 0.00097
      CONSTS(29) = 0.57
      CONSTS(30) = 0.0017
      CONSTS(31) = 0.0986
      CONSTS(32) = 0.024
      CONSTS(33) = 3e-5
      CONSTS(34) = CONSTS(31)+CONSTS(29)+CONSTS(30)
      CONSTS(35) = CONSTS(32)+CONSTS(27)
      CONSTS(36) = CONSTS(33)+CONSTS(28)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = (CONSTS(2)+(CONSTS(3)+ CONSTS(4)*STATES(1)+ CONSTS(5)*STATES(1) ** 2.00000)/(1.00000+ CONSTS(6)*STATES(1)+ CONSTS(7)*STATES(1) ** 2.00000+ CONSTS(8)*STATES(2)+ CONSTS(9)*STATES(2) ** 2.00000)) -  CONSTS(1)*STATES(1)
      RATES(3) = (CONSTS(11)+ CONSTS(12)*STATES(1)+ CONSTS(13)*STATES(1) ** 2.00000)/(1.00000+ CONSTS(14)*STATES(1)+ CONSTS(15)*STATES(1) ** 2.00000+ CONSTS(16)*STATES(2)+ CONSTS(17)*STATES(2) ** 2.00000) -  CONSTS(10)*STATES(3)
      RATES(2) = (CONSTS(18)+( CONSTS(19)*STATES(1)+ CONSTS(20)*STATES(1) ** 2.00000+ CONSTS(21)*STATES(3)+ CONSTS(22)*STATES(3) ** 2.00000)/(1.00000+ CONSTS(23)*STATES(1)+ CONSTS(24)*STATES(1) ** 2.00000+ CONSTS(25)*STATES(3)+ CONSTS(26)*STATES(3) ** 2.00000)+ CONSTS(27)*STATES(4)+ CONSTS(28)*STATES(5)) -  CONSTS(34)*STATES(2)
      RATES(4) =  CONSTS(29)*STATES(2) -  CONSTS(35)*STATES(4)
      RATES(5) =  CONSTS(30)*STATES(2) -  CONSTS(36)*STATES(5)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END