a = sprandsym(n,p);
a(find(a)) = 1;
a(1:(n+1):end) = 0;
degrees = sum(a);
x = min(degrees):1:max(degrees);
hist(degrees, x);
xlabel('Links');
ylabel('Nodes');
title('Food Web Simulation');