% Program to Forecast Denver Broncos Number of Wins 

clc;
p=zeros(16,2);
alpha = 0.5;
n = 16;
for x  = (1:1:n)
    p(x,1) = x;
    p(x,2)=(factorial(n)/(factorial(x)*factorial(n-x)))*(alpha.^x)*((1-alpha).^(n-x));
    %p(x,2)=nchoosek(n,x)*(alpha.^x)*((1-alpha).^(n-x));
end

disp(['    Wins    Probability']);
disp([p]);