% Written by:
% -- 
% John L. Weatherwax                2005-08-14
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

clc; 

pE7 = nchoosek(13,7)*nchoosek(39,6)/nchoosek(52,13)

% compute P(F): 
% 
t_sum = 0.0;
for ii=1:12, 
  t_sum = t_sum + (ii/13) * nchoosek(13,ii) * nchoosek(52-13,13-ii) / nchoosek(52,13); 
end
t_sum = t_sum + 1 / nchoosek(52,13)

% compute P(E_7|F)
%
( 7 / 13 ) * pE7 / t_sum