A robust and efficient iterative method for the numerical solution of