
%% Il problema si risolve 
% 1) calcolando Prob(zero componenti difettosi) su quattro estrazioni
% 2) calcolando Prob(4 componenti funzionanti) su quattro estrazioni

N=100; %100 componenti
K=4;    % componenti estratti

%% Soluzione con ipergeometrica

% la soluzione del problema e': Prob(zero componenti difettosi)
M=5;   % 5 difettosi
p_comp_difettoso = hygepdf(0, N, M, K)

%oppure: Prob(4 componenti funzionanti)
M=95;   % 95 funzionanti
p_comp_funzionante = hygepdf(4, N, M, K)

%% Soluzione senza ipergeometrica con regola del prodotto
% calcolando Prob(4 componenti funzionanti) su quattro estrazioni
% 
M=95; % 95 funzionanti

p_comp_funzionante = (M/N) * (M-1)/(N-1) * (M-2)/(N-2) * (M-3)/(N-3)

% oppure in forma iterativa

p_comp_funzionante = 1;
for k=0:3
    p_comp_funzionante = p_comp_funzionante * (M-k)/(N-k);
end
p_comp_funzionante


