Approximate solution of normal equations by eigenvalue decomposition