function x = Jacobi(A,b,x0,TOL,Nmax)
    % your function here
end
