Separate out the matrices from the model code.
Squash a few bugs...
Scripts to run the Huisman model, including paper figures