function x=JacobiIt2(A,b,omega)
D=diag(diag(A));
U=triu(A,1);
L=tril(A,-1);
if det(D)==0
    error('Van a főátlóban 0')
end
B=eye(size(A))*(1-omega)-omega*inv(D)*(L+U);
f=omega*inv(D)*b;
xpontos=linsolve(A,b);
x=zeros(size(b));
for k=1:10
    x=B*x+f;
    maradek(k)=norm(x-xpontos);
end
rho=max(abs(eig(B)));
fprintf('A spektrálsugár logaritmusa: %4.2f\n',log(rho));
figure(2)
plot(1:10, log(maradek),'r+')