N = 100;
x = ones(N,1);
tol = [1e-6,1e-6];
[sol, it_hist, ierr] = nsolgm(x,@f,tol);
semilogy(it_hist(:,2),it_hist(:,1)/it_hist(1,1),'o-.');
hold on;
%grid on;