add ability to specify which numerical method to use in solve_kepler_E

[mdbernard]

In kepler.py, the function solve_kepler_E uses Newton's Method. Two changes should be made:

• Change the function to use numerical.newton
• Add the ability to specify which numerical method to use