function x=seidel(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(1)=B(1,:)*x0+c(1)
    for j=2:n
      x(j)=B(j,1:j-1)*x(1:j-1)'+B(j,j:n)*x0(j:n)+c(j)
    endfor
    x=x'
    while q/(1-q)*norm(x-x0,inf)>hiba
      i=i+1
      x0=x;
      x(1)=B(1,:)*x0+c(1);
      for j=2:n
        x(j)=B(j,1:j-1)*x(1:j-1)+B(j,j:n)*x0(j:n)+c(j);
      endfor
      x
    endwhile
    else disp('az együtthatómátrix nem négyzetes vagy diagolásian domináns')
  endif
endfunction