function x=seidel(A,b,x0,n)
  h=length(x0);
  for i=1:h
    B(i,:)=A(i,:)/(-A(i,i));
    B(i,i)=0;
    c(i)=b(i)/A(i,i);
  endfor
  c=c';
  for i=1:n
    x=zeros(h,1)
    x(1)=B(1,:)*x0+c(1)
    for j=2:h
      x(j)=B(j,1:j-1)*x(1:j-1)+B(j,j:h)*x0(j:h)+c(j)
    endfor
    x0=x
  endfor