%-- 2011/06/20 13:41 --% A = [4, -3; -3, 3] eig(A) (7+sqrt(37))/2 (7-sqrt(37))/2 A = [-7, -9; -9, -9] eig(A) A = [-7, -18; 0, -9] eig(A) tmp = rand(5,5) A = (tmp+tmp') [V,D] = eig(A) tmp = rand(3,3) [V,D] = eig(A) A = [1, 2, 3; 2, 4, 5; 3, 5, 9] [V,D] = eig(A) V\D*V V*D/V V*V' V'*V