function Newtrp_df_1(fx,x)
% Newton-Raphson por diferencias finitas para sistemas
%
 global h
 tol=sqrt(eps); dnor=1.0; h=tol;
 while dnor >tol
    [f J]=fx(x);
    p=J\f';
    x=x-p;
    dnor=norm(fx(x));
    fprintf('%15.10e %15.10e %15.10e %15.10e\n', x,dnor);
 end
end
