This is a pep-compliant version of Dan Lenski's
proposal. See for example

https://stackoverflow.com/a/33181083
https://mail.python.org/pipermail/stdlib-sig/2015-July/000990.html

for the implementation in python 3

This is taken from

 https://mathematica.stackexchange.com/a/104696