% f=inline('(x-1).^3'); df=inline('3*(x-1).^2'); x0=1.5;
f=inline('x.^2-2'); df=inline('2*x'); x0=1;
tolstep=10^(-15); maxit=10000;
[x,ko]=newton2011(x0,tolstep,maxit,f,df);
format long e; x