A=[2 8 10; 8 3 4 ;10 4 7]
xk=ones(3,1);
err=1;
tol=1e-5;

while(err>tol)
yk=A*xk;
lambda=yk(1);
yk=yk./yk(1);
err=abs((yk-xk)./yk)*100;
err=max(err);
xk=yk;
end
lambda
evector=xk