% Setting up data;

Y = hist(degMat,4);
X = 1:length(Y);

% Plotting the data
plot(X,Y,'+r'), hold on
% Polynomial curve fitting of log values so that we have a linear equation
% Simplifying, log(Y) = log(a*X.^k) = log(X)*k + log(a)
p = polyfit(log(X),log(Y),1); 
% Evaluating Coefficients
k = p(1);
% Accounting for the log transformation
a = exp(p(2));
% Final plot
fplot(@(X) a*X.^k,[X(1) X(end)])