function prd = prdistribution(x)

N = size(x,1);
prd = zeros(ceil(max(x))*(10^5),1);
for i = 1:N
  prd(ceil((x(i))*(10^5)),1) = prd(ceil((x(i))*(10^5)),1) + 1;
end;

loglog(prd,'r.-');
%loglog([1:size(prd,1)],prd,'r.-',[1:100:1000],([1:100:1000].^(-2.1))*10000,'k.-');
grid on;
title('Stanford-Berkeley Web');
xlabel('PageRank values (x10^-^5)');
ylabel('Number of pages');
