%
% Example 5 from Sec 5.3
%
A = [5 -8 1; 0 0 7; 0 0 -2]
[P,D] = eig(A)
det(P)
P * D * inv(P)
