function[V,L]=myeig(A)
L=eig(A);
n=size(L,1);
V=zeros(size(A));
for i=1:n
    V(:,i)=null(L(i)*eye(size(A))-A,'r');
end
end