Root finding with scipy.optimize: Halley's comet from Kepler's equation

Afterwards you can solve an equation in one variable with scipy.optimize, choose a bracket for brentq or a start for newton, and solve a whole array at once.