%%
% The |eig| command will give just the eigenvalues if one output is given,
% or both V and D if two are given.
A = ones(2,2);
lambda = eig(A)

%%
[V,D] = eig(A)

%%
% We can check the fact that this is an EVD.
norm( A - V*D*V' )