closes #24. Brings script for fig. 3 into the code

[eric-pedersen]

This code addresses issue #24 by bringing the script for this figure into the main code.