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