function raiz=Newton_composite(fun,x,tol)
% Newton composite: misma derivada cada dos iteraciones
 in=1; x1=0; x2=2;
 if nargin<3, tol=eps^0.5; end
 if nargin==2, x1=0; x2=x; end
 while abs(x2-x1)>tol
    x1=x2;
    if in
       [f df] = fun(x1); in = 0;
    else
       [f] = fun(x1); in = 1;
    end
    x2=x1-f/df;
    fprintf('%18.15f\n', x2);
 end
 raiz=x2;
end
