function [] = chap_3_prob_13()
% 
% Written by:
% -- 
% John L. Weatherwax                2006-08-28
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----
  
%format rat; 

p1 = ( 4 * nchoosek(48,12) )/ nchoosek(52,13)
p2 = ( 3 * nchoosek(36,12) )/ nchoosek(39,13)
p3 = ( 2 * nchoosek(24,12) )/ nchoosek(26,13)
p4 = 1;

p1*p2*p3*p4