added the instructional pdf that accompanies this project
implemented a better equilibrium test. It works well and outputs a nice results summary. Need to consider how size of prevE arrays and max_pc_diff var affect the point at which equilibrium is reached
updated the alias file for project-b
added functions to output mesh at intervals for visualising
removed old equilibrium code that was commented out. Began refactoring to run multiple simulations for a range of beta vars. Need to get hot and cold running simulataneously
added formulae to calculate more system values. Need to devise a clear plan for rest of project