A = [1 2 0; 2 5 -1; 4 10 -1] B = A' C = A * B X = inv(A) I = inv(A) * A eig(A) svd(A) whos