function gausseseidel(A,B,x0)
 x=x0;
 len=length(A);
for i=1:10
    for j=1:len
        sum=0;
        for k=1:len
            if(k==j)
            else
                sum=sum+A(j,k)*x(k);
            end
        end
        x(j)=(1/A(j,j))*(B(j)-sum);
       fprintf('i:%d ',i);
       for k=1:len
           fprintf('x%d:%d ',k,x(k));
       end
       fprintf('\n');
    end
end
end