function oProb = oHyper(N, K, P)
%OHYPER

oProb = zeros(K,1);
for k = 1:K
oProb(k) = nchoosek(K,k)*P^k * nchoosek(N-K,K-k)*(1-P)^(K-k);
end
oProb = oProb / sum(oProb);

end

clear;clc
format LONG

%%
oProb = oHyper(30, 15, 0.5);
hProb = hygepdf(1:15,30,15,15)';
disp([oProb hProb]);

%%
oProb = oHyper(30, 15, 2/3);
fprintf('\\nOracle distribution for P = 2/3\\n')
disp(oProb);
fprintf('\\nOracle prob for 8 or less successes %5f\\n', sum(oProb(1:8)));

Oracle distribution for P = 2/3
0.000000010080107
0.000000987850518
0.000037099275003
0.000667786950053
0.006464177676510
0.035912098202832
0.118729794058344
0.237459588116687
0.287296785622658
0.206853685648314
0.085476729606741
0.018994828801498
0.002023117860515
0.000082576239205
0.000000734011015

Oracle prob for 8 or less successes 0.399272