C++ code converted to python

Basic tests with planet at y=0 work - now to write
more stringent tests