Scattered corrections and improvements, and tests for python and julia versions

I have added a number of scattered corrections and improvements in both versions, including updates in the Readme.md files. I have also included tests in julia and made as compatible as possible with those in python. Still, there are some things to check.
Once you have restricted `a` to the correct domain, you should just be able to call the power function on the new interval, no?