% ################################################################
% Supporting file for top level program percolation.m.
% ###############################################################
function ret = f(x, b, n)
ret = sum(b*x.^b.*log(x))./sum(x.^b) - sum(b*log(x))/n -1;
return