Least squares solution of large Toeplitz linear systems