TITLE Two Compartment Model
:  National Biomedical Simulation Resource, 8-4-89
:  Reference - Randall, Microcomputers and Physiological Simulation,
:    Second Edition, Raven Press, 1987, Chapter 11.
:  Default output - ConcPlasma decreases from about 17 mg/l at time=0 to
:    9 mg/l at time=30 in a double exponential curve.  ConcTissue rises
:    from zero at time=0 to a broad maximum of about 27 mg/l at time=5.4
:    and falls to 13 mg/l at time=30.
: $Header: /usr/local/src/oldscop/scop/examples/RCS/twocomp.mod,v 13.3 1997/02/17 22:31:46 kootsey Exp $

PARAMETER
{
  VolPlasma = 3         (L)
  VolTissue = 9         (L)
  RateP = .5            (L/hr)
  RatePT = 1.0          (L/hr)
  VolInject = 5         (mL)
  ConcInject = 100      (mg/mL)
  dose_rate = 0.1       (mg/hr)
  dose_dur = 0.5        (hr)
  dose_per = 4          (hr)
}

INDEPENDENT
{
  time FROM 0 TO 30 WITH 100  (hr)
}

PLOT AmtTissue VS time

ASSIGNED
{
  ConcPlasma FROM 0 TO 200   (mg/L)
  ConcTissue FROM 0 TO 200   (mg/L)
}

STATE
{
  AmtPlasma  FROM 0 TO 500   (mg)
  AmtTissue  FROM 0 TO 500   (mg)
}

INITIAL
{
  AmtPlasma = VolInject * ConcInject
  ConcPlasma = AmtPlasma / VolPlasma
  ConcTissue = AmtTissue / VolTissue
}

DERIVATIVE flows
{
  AmtPlasma' = - RateP * AmtPlasma/VolPlasma + RatePT *
		(AmtTissue/VolTissue - AmtPlasma/VolPlasma)
                + perpulse(time, 0, dose_rate, dose_dur, dose_per)
  AmtTissue' = - RatePT * (AmtTissue/VolTissue - AmtPlasma/VolPlasma)
}

BREAKPOINT
{
  SOLVE flows
  ConcPlasma = AmtPlasma / VolPlasma
  ConcTissue = AmtTissue / VolTissue
}