function X=gausseidel(A,b,x0,tol)
[m,n]=size(A);
xn0=x0;
xk=x0;
paro=Inf;
X=zeros(n,1);
itergaus=0;
while paro>tol
    xn0=xk;
    for i=1:n
        suma1=A(i,1:i-1)*xn0(1:i-1);
        suma2=A(i,i+1:n)*xk(i+1:n);
        xk(i)=(b(i)-suma1-suma2)/A(i,i);
        paro=norm(xk-xn0);
    end
    itergaus=itergaus+1
end
X=xk;
end