A = [4 2 4; 2 0 5; 5 1 -6];
lambda = eig(A)'
for k = 1:150
    a(k) = norm((A^k)^(1/k),1);
    b(k) = norm((A^k)^(1/k),2);
    c(k) = norm((A^k)^(1/k),inf);
end
k = 1:150; plot(k,a,k,b,'--r',k,c,'-.k')
legend('1-norm','2-norm','\infty-norm')
