function raiz=Newton_sec(fun,x,tol)
%Newton secante para una variable
 x1=x; x2=x-0.1; if nargin<3, tol=eps^0.5; end
 while abs(x2-x1)>tol
    x0=x1; x1=x2;
    [f df] = fun(x0,x1);
    x2=x1-f/df;
    fprintf(' %18.15f\n', x2);
 end
 raiz=x2;
end