%%
clear
clc

%%
% 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')