function x=jacobi(A,b,x0,hiba)
  [n,m]=size(A)
  for i=1:n
    B(i,:)=A(i,:)/(-A(i,i))
    B(i,i)=0
    c(i)=b(i)/A(i,i)
  endfor
  c=c'
  q=norm(B,inf)
  if q<1 & n==m
    i=1
    x=B*x0+c
    while q/(1-q)*norm(x-x0,inf)>hiba
      i=i+1
      x0=x;
      x=B*x0+c
    endwhile
    else disp('az együtthatómátrix nem négyzetes vagy diagolásian domináns')
  endif
endfunction