Is there a way to identify the eigenvectors of a non-normal matrix using QR?

Question

After applying the Hessenburg reduction to A. I\'m trying to use the Wilkinson shift to calculate the eigenvalues/eigenvectors of a matrix A.

I found a unitary matri