Mathematica script
Wed Sep 18 11:11:27 BST 2013