function [manHist maxHist] = getGamma(Number,gam_mode,gam_mu)


if nargin < 3
    gam_mu=3.6;
    gam_mode=1.6;
    Number = 120;
elseif nargin < 2
    gam_mode=1.6;
    Number = 120;
elseif nargin < 1
    Number = 120;
end


k=1/(gam_mu-gam_mode); % k = .5
m=gam_mu/(gam_mu-gam_mode); % m = 1.8

x=gamm_rnd(Number,m,k);

% for i=1:4
%     b=size(find(ceil(x)==i),1);
%     if b>=27
%         b
%     end
% end


%hist(ceil(x),25)

manHist = 0;
x = ceil(x);

for i = 1:max(x)
    
   manHist(i) = length(find(x==i));
    
end

maxHist = max(manHist);