C = [1 0 0; 0 -2 1; 1 2 1]
eig(C)
[P, L] = eig(C)