 [x, f, ok, count] = newtonND('fun', 'grad', 'hessian', [4; 5], 1e-5, 20)