A=[2 1 1 0; 1 3 1 1; 1 1 4 1; 0 1 1 5];
y0=[1 1 1 1]';

mocninna_metoda_norm(A, y0, 0.001, 50);