function [ Pi ] = Pi( idx, Clusters )
    seq=size(idx, 2);
    Pi=ones(1, max(Clusters));
    for i=1:seq
        Pi(Clusters(idx(i)))=Pi(Clusters(idx(i)))+1;
    end
    
    disp(['     Valeur de Pi avant normatisation:  ' num2str(Pi)]);
    % Normalise Pi
    Pi=Pi/sum(Pi);
    disp(['     Probabilité de Pi:  ' num2str(Pi)]);
end

