function x = GS(A,b,x0,TOL,Nmax)
    % Your function here
end
