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