function y=gseidel(A, b,x0,tol)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Aproxima una raiz del sistema Ax=b por el metodo de Gauss-seidel.
%% Las variables de entrada son:la matriz de coeficientes A, 
%el término independiente b, una estimación inicial x0 y una 
%cierta tolerancia tol para el resultado de sustituir la iteración  
%en el sistema.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

U=triu(A,1);
L=tril(A,-1);
d=diag(diag(A));
x=x0;
i=1;
while max(abs(A*x-b))>tol
   x=(L+d)\(b-U*x);
   y=[i,x'];
   i=i+1;
end