function [xn,err,k] =gaussseidel(A, b, x0, tol, maxiter) 

% x0 e um vector coluna

    % n e a dimensao da matriz

    n=length(A);

    % definir um matriz coluna Xn de zeros 

    xn=zeros(n,1);

    %inicialmente o erro tem que ser maior que tol

    erro(1)=tol+1;

    % K e igual a 2 pois o erro so e calculado a partir de xn(2)
    k=2;


    while erro(k-1)>tol & k<maxiter

        for i=1:n

            s=0;

            % soma da primeira parcela da matriz dividida por xn

            % e mais preciso do que o jacobi, pois utiliza o xn ja calculado

            for j=1:i-1

                s=s+A(i,j)*xn(j);


            end
            % soma da segunda parcela da matriz dividida por x0

            for j=i+1:n

                s=s+A(i,j)*x0(j);

            end

            % achar o vector xn

            xn(i)=(b(i)-s)/A(i,i);

        end

        % norma infinito

        erro(k)=norm(xn-x0,inf)/norm(xn,inf);

        % o x0 passa a ser o xn

        x0=xn;

        % acrescentar 1 unidade a k

        k=k+1;

    end
return