A=[1 1 0 0; 1 2 0 1; 0 0 3 3; 0 1 2 3]; X=[1;1;0;1]; n=4; tol=0.001; k_old=9999999; k_new=0; while abs(k_old-k_new)>tol y=A*X; k_old=k_new; k_new=abs(y(1)); for i=1:n if k_new<abs(y(i)) k_new=abs(y(i)); end end X=y/k_new; end k_new X