Gauss-Seidel Iterative

function apX=GaussSeidelIterative(n,A,b,XO,TOL,N)
k=1;
x=zeros(n,1);

while k<N
    for i=1:n
        x(i)=((-1)*(A(i,1:i-1)*x(1:i-1)+A(i,i+1:n)*XO(i+1:n))+b(i))/A(i,i);
    end
    k=k+1;
    for i=1:n
        XO(i)=x(i);
    end
   k
   x
end