function Broyden_3(fx,x)
% Método cuasi Newton con la fórmula de Broyden
%
 tol=sqrt(eps); n=length(x); dnor=1.0; J=zeros(n,n); f=fx(x); 
 J(1,1)=3; J(2,2)=-178.2; J(3,3)=20;   % Trampa para empezar mejor
 while dnor >tol
    p=J\f';
    x1=x-p;
    f=fx(x1);
    dnor=norm(f);
    fprintf('  %15.10e %15.10e %15.10e %15.10e\n',x,dnor);
    J=broy(J,f,p);
    x=x1;
 end
end

function J=broy(J,f,p)
 prod=p'*p;
 J=J-(1/prod)*f'*p';
end
