Computational Thinking, Linear Algebra, Python, Jupyter Notebook, Math.
np.linalg.eigvals(A) array([ 2.78023+0.j , -0.11832+0.j , 0.26911+0.44246j, 0.26911–0.44246j, 0.07454+0.49287j, 0.07454–0.49287j]) Check that Q is orthogonal: np.allclose(np.eye(n), Q @ Q.T), np.allclose(np.eye(n), Q.T @ Q) (True, True) This is really really slow. Practical QR (QR with shifts) The