SAGE/Python code written by NoviceMS, Copyright December 2015