function LSQ_New_1(fx,x)         
% Mínimos cuadrados Newton
%
 tol=sqrt(eps); p=ones(2,1); [f grad Hess] = fx(x);
 while abs(grad'*p)>tol
    p=-Hess\grad;
    x=x+p;
    [f grad Hess] = fx(x);
    fprintf(' %15.10e %15.10e %15.10e\n',x,f);
 end
end
