x = bisect(@ff,1,2);
disp(x);

[xstar, gx] = fixpoint(@gg,4);
disp(xstar);
disp(gx);

x = newton(@ff2,1);
disp(x);

x = secante(@ff,1);
disp(x);

[x,fval,exitflag] = fzero(@ff,1)

%opts = optimset('Display','iter');
%x = fzero(@ff,1,opts);

x0 = [0,0];
[x, f, ef] = fsolve(@sistemanl,x0)

sistemanl(x)

% opts = optimset('Display','iter');
% x0 = [0,0];
% [x, f, flag]= fsolve(@sistemanl,x0,opts)