Quick Start
-----------
1. Change Matlab's current folder to this one.
2. Type 'simulateTrajectories' in Matlab's command window to generate Fig. 2.
3. Type 'simulateTrajectories('pk', '3dbox', 90, 3000)' for another
   simulation.

About
-----
These are accompanying files to "Benchmark Problem: A PK/PD Model and Safety 
Constraints for Anesthesia Delivery" by Victor Gan, Guy A. Dumont and Ian M.
Mitchell for the Workshop on Applied Verification for Continuous and Hybrid
Systems.

The simulation allows for generating the trajectories for a set of initial
points obeying a predefined DE (equation 1 or equation 5). A safe region of the
state space is also defined. If any of the trajectories leave the safe region,
the trajectory is coloured dashed red. Alternatively, the trajectory is coloured
solid green if it stays exclusively within the safe region.

The objective of this benchmark is to find the maximal set of states in the safe
region where there exists an input signal u that leads to a trajectory
completely enclosed in the safe region for a given time horizon (a "solid green"
trajectory).

Requirements
------------
Tested with MATLAB R2012a / version 7.14

File Descriptions
-----------------
simulateTrajectories.m  The main function
getTrajectories.m       Finds the trajectories from a set of initial points
getViablePoints.m       Indicates whether trajectories stay within a safe region
pkModel.m               The DE for the PK model (equation 1)
pkpdModel.m             The DE for the simplified PKPD model (equation 5)
plotConcentrationTime.m Plots concentration-time graphs for each state
plotTrajectories.m      Plots the trajectories in a 3-dimensional state space
README.txt              This file
uPaper.m                An example input function used in Fig. 2.

Contact
-------
Let me know if you run into any problems. 
Victor Gan vhg@cs.ubc.ca
