count_gp.m

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-1)*7<X & X<m*7));
end
%[COUNT,IX]=sort(COUNT,2);
%COUNT=[IX;COUNT];
%if(flag==1)
 %   COUNT=fliplr(COUNT);
 %end