function [M,alpha] = economics_polyfit(PDF)

number_rounds = size(PDF,2);

for i = 1:number_rounds
polyfit_x(i,:) = log(PDF{i}(1,end-3:end));
polyfit_y(i,:) = log(PDF{i}(2,end-3:end));

polyfit_M(i,:) = polyfit(polyfit_x(i,:),polyfit_y(i,:),1);
end

M = mean(polyfit_M(:,1));
%alpha = -M - 2; THIS IS THE DEFINITION USED IN EARLIER PAPERS; WE WILL
%AVOID THAT HENCEFORTH
alpha=-M;

end