C C There are a total of 0 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 11 entries in the constant variable array. C C C VOI is time in component environment (minute). C STATES(1) is X1 in component X1 (ng_ml). C CONSTS(1) is a0 in component X1 (ng_ml_min). C CONSTS(2) is a1 in component X1 (ng_ml_min). C CONSTS(3) is a2 in component X1 (dl_microg). C CONSTS(4) is b1 in component X1 (first_order_rate_constant). C STATES(2) is X3 in component X3 (microg_dl). C STATES(3) is X2 in component X2 (pg_ml). C CONSTS(5) is a3 in component X2 (pg_ml_min). C CONSTS(6) is a4 in component X2 (first_order_rate_constant). C CONSTS(7) is a5 in component X2 (dl_microg). C CONSTS(8) is b2 in component X2 (first_order_rate_constant). C CONSTS(9) is a6 in component X3 (microg_dl_min). C CONSTS(10) is a7 in component X3 (first_order_rate_constant). C CONSTS(11) is b3 in component X3 (first_order_rate_constant). C RATES(1) is d/dt X1 in component X1 (ng_ml). C RATES(3) is d/dt X2 in component X2 (pg_ml). C RATES(2) is d/dt X3 in component X3 (microg_dl). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 10.0 CONSTS(1) = 0.0014 CONSTS(2) = 0.000517 CONSTS(3) = 0.0164 CONSTS(4) = 0.0598 STATES(2) = 0.0 STATES(3) = 0.0 CONSTS(5) = 1.38 CONSTS(6) = 0.60 CONSTS(7) = 0.00498 CONSTS(8) = 0.053 CONSTS(9) = 0.0084 CONSTS(10) = 0.0081 CONSTS(11) = 0.0138 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(1) = (CONSTS(1)+CONSTS(2)/(1.00000+ CONSTS(3)*STATES(2))) - CONSTS(4)*STATES(1) RATES(3) = (CONSTS(5)+ CONSTS(6)*STATES(1))/(1.00000+ CONSTS(7)*STATES(2)) - CONSTS(8)*STATES(3) RATES(2) = (CONSTS(9)+ CONSTS(10)*STATES(3)) - CONSTS(11)*STATES(2) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RETURN END