A=[2 0 0; 0 1 0; 0 -2 3];  
y0=[1 1 1]'; 

mocninna_metoda(A, y0, 0.001, 50);

[Vlastni_vektory,Spektralni_matice] = eig(A)
