Using IPython notebook of course:

http://nbviewer.ipython.org/3843014/