
% f=inline('(x-1).^3'); a=0; b=1.5;
f=inline('x.^2-2'); a=1; b=2;
tolres=10^(-15);tolintv=10^(-15);maxit=10000;
[a,b]=bisezione2011(a,b,tolintv,tolres,maxit,f);
format long e; [a b]