Multi-dimensional Newton-Raphson extension with global convergence

No implementation of the paper was online so far, so here we go