count_gp.asv

for double color balls, it can help you make the list with high p

function COUNT=count_gp(k,X,flag)
% flag = 0, count and sort in descending order
% flag = 1, count and sort in ascending order
% k --- the count of items
% X --- Sources
COUNT=1:1:k;
for m=1:k
    COUNT(m)=length(find(m*7<X<(m+1)*7));
end
[COUNT,IX]=sort(COUNT,2);
COUNT=[IX;COUNT];
if(flag==1)
    COUNT=fliplr(COUNT);
end