function [x,dnor,iter]=Newt_Rap(funcion,x0,tol,maxiter)
% Newton-Raphson para cualquier sistema no lineal 
%
 if nargin<4, maxiter=99999; end
 if nargin<3, tol=sqrt(eps); end
 dnor=1.0; iter=0; A=[]; Dnor=[];
 while dnor>tol && iter<maxiter
    iter=iter+1;
    [f,J]=funcion(x0);
    p=J\f';
    x=x0-p; [f]=funcion(x);
    dnor=norm(f,inf); A=[A;x']; x0=x; Dnor=[Dnor;dnor];
 end
 s=' %15.10e'; for i=1:size(x0)-1 s=[s ' %15.10e']; end
 for i=1:iter
    fprintf([s ' %15.10e\n'], A(i,:),Dnor(i));
 end
end
