
load fellowship_record
fr = fellowship_record; % abbreviate

n = size(fellowship_record,2); % number of students in sample

% see nutshell "Intro to Proba..." Section 4

eventW1W2 = fr(1,:) & fr(2,:) ; 
% eventW1W2 will be assigned a 1 in index i only if both fr(1,i) == 1 and fr(2,i) == 1 
% (i.e., if a student wins both fellowships); otherwise fr(1,i) will be assigned a 0 in index i

% recall ":" expands to "[1:end]" so fr(1,:) is a single-index row array
% which is just the first row of fr --- hence, the results for fellowship 1 for the n
% students in the sample

phiW1W2 = sum(eventW1W2)/n; % the "sum" counts ones in the logical array eventW1W2
% to calculate the number function; the 1/n then normalizes to calculate frequency

% now proceed to the other events requested and create logical arrays and then frequencies


