hypergeometric

x = 0:4; % number of parts from local supplier
m = 300; % Total number of parts
k = 100; % Total parts from local supplier
n = 4; % number of parts selected
local = hygepdf(4,m,k,n); % all parts from local supplier
r=hygecdf(4,m,k,n)-hygecdf(1,m,k,n); %two or more parts from local supplier
Px = hygepdf(x,m,k,n); %P(parts from local supplier)
plot(x,Px,'r');
title("Hypergeometric Probability Distribution");