I had implemented it once in 2 lines of Python. some 25 years ago. If I can shake the rust maybe I'll be able to rederive it. 

Please do! I would love to see a concise and readable implementation in Python.