On the solution of normal equations