everything up to and including CS281A final report simulations and plotting; port grid_model to PyMC3 with HMC and MH sampling methods