function [val,der,H] = maxrouting(x,m,nf)eps = 1e-6;f = x(1:nf);val = -sum(log(eps+f));if nargout > 1	der = [-1./(eps+f);zeros(m,1)];	if nargout > 2		H = diag(1./((eps+f).^2));		H = [H,zeros(nf,m);zeros(m,nf),zeros(m,m)];	endend