function [x,pr]=aggregate(x,T)
% produces bins and probabilities for convoluation of simple random sample (with replacement) from vector x
% taken T at a time
s=[]
for i=1:10000;
   s=[s sum(sample(x,T))];
end
[pr,x]=hist(s,50);
pr=pr/sum(pr);