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