Iterative solutions of normal equations