Computational Physics, Linear Algebra, Math, Mathematics, Science.
A array([[6, 2], [2, 3]]) eigVals, eigVecs = np.linalg.eig(A) eigVecs array([[ 0.89442719, -0.4472136 ], [ 0.4472136 , 0.89442719]]) eigVals = np.diag(eigVals) eigVals array([[7., 0.], [0., 2.]]) eigVecs.dot(eigVals).dot(eigVecs.T) array([[6., 2.], [2., 3.]]) We can see that the result corresponds