Made original script into Ipython Notebook, analysis now includes rho