The Newton-Raphson second order method for a multivariables function can be written with the equation (I) form.
The solution is obtained with an efficient linear solving method like the partial pivoting gauss algorithm.
The second order Taylor's development gives the exact approximation of the function.
We have found with several examples that the above method converges in two iterations and the solution is obtained in the first iteration.
The conclusion is that this solving method corresponds to the exact one.
Another way is to obtain the solution with the inverting process as it is shown in equation (II) : this method can not always converge easily and it depends on the Hessian matrix conditioning because the inverse is a source of roundoff errors.