A=[1 1 0;1 2 1;0 1 -1];
%x=[1 1 1]';
evector=[1 1 1]' %initial value
tolerance=1e-6;
err=1;

while (max(err)>tolerance)
dummy=evector;
product=A*evector;
evalue=product(1);
evector=product./evalue;
err=abs((evector-dummy)./evector);
end

evalue
evector