function [res] = ranking_degree( input_matrix )
% Plot node degree ranking
%   input_matrix: symmetric contact matrix

s = size(input_matrix);

res=zeros(1,s);

for ii= 1:s
    for jj=1:s
        if(ii<jj)
            if(input_matrix(ii,jj)>0)
                res(ii)=res(ii)+1;
            end
        end
    end
end

res=sort(res,'descend');
plot(res,'ko','MarkerSize',10);

xlabel('Ranking');
ylabel('Degree')

return 
end