I am trying to implement the newton method for maximization in higher dimensions and I was wondering if there exists any solvers for this in Python? In Scipy there is a solver for the 1-dimensional case, but I do not see one for the multi-dimensional case. I suppose that it is possible to implement it using the Hessian and Gradient solvers in Numdifftools
EDIT: It looks like scipy.optimize.minimize does this. I was looking under the multi-dimensional heading and it wasn't there, that's why I missed it. It was under the general-purpose heading
