 Unshifted QR algorithm

The unshifted QR algorithm factors $A_k=Q_k{R}_k$ and forms $A_{k+1}=R_k{Q}_k=Q_k^T{A}_k{Q}_k$, preserving eigenvalues while driving the matrix toward block diagonal form.