clear
A = [2 -1 0;-1 2 -1;0 -1 2]
pause
[u  d] = eig(A)
pause
disp('column one of u times k = 1/u(2,1)')
u(:,1)/u(2,1)
pause
disp('column two of u times k = 1/u(1,2)')
u(:,2)/u(1,2)
pause
disp('column three of u times k = 1/u(2,3)')
u(:,3)/u(2,3)
pause
A = [ -1 2; -1 -3]
pause
[u d] = eig(A)
pause
disp('column one of u times k = 1/u(1,1)')
u(:,1)/u(1,1)
pause
disp('column two of u times k = 1/u(2,2)')
u(:,2)/u(1,2)
pause
disp('eigenvalues and eigenvectors of the three-tank mixing problem')
A=[-10  4 0;10 -12 2;0 8 -8]
pause
[u d ] = eig(A)