function optomega=optomegasz(A)
legjobb=1;
D=diag(diag(A));
U=triu(A,1);
L=tril(A,-1);
hold on
for omega=0.01:0.01:2
B=eye(size(A))*(1-omega)-omega*inv(D)*(L+U);
plot(omega,max(abs(eig(B))),'r+')
if max(abs(eig(B)))<legjobb
        optomega=omega;
        legjobb=max(abs(eig(B)));
    end
end