% input della matrice
a=[2 -5;1 -4]
% calcolo degli autovettori normalizzati (colonne della matrice x)
% e degli autovalori (diagonale della matrice lambda)
[x,lambda]=eig(a)
% estraiamo dalla matrice x il primo autovettore v1
v1=x(:,1)
% estraiamo dalla matrice x il primo autovettore v1
v1=x(:,2)
% estraiamo dalla matrice lambda la diagonale
diag(lambda)
% lambda e' la diagonalizzazione della matrice a
% vale  a = x*lambda*inv(x)  e  lambda = inv(x)*a*x:
x*lambda*inv(x)
a
inv(x)*a*x
% (che e' uguale a lambda salvo piccoli errori di calcolo)
