A modified Newton method with cubic convergence: the multivariate case