results = zeros(2,10);  % first value = first draw, 10th = 10th draw
for i = 1:100000 % 100000 trials
   balls = randperm(16,16);  % if the value is less than <= 6 it is white
                        % if the value is > 6 it is red
   if(balls(1) > 6)    % check for value of the first ball
    results(1,i) = 1;  % store a 1 if it is red
   end 
   if(balls(10) > 6)  % check for value of the tenth ball
     results(2,i) = 1;  % store a 1 if it is red
   end
end


probs = sum(results,2)/i; % divide the sum of the total number of red balls drawn for each ball position, and divide by total trials
disp('The suggested probability of the first ball being red (found through sampling) is')
disp(probs(1))
disp('The suggested probability of the tenth ball being red (found through sampling) is')
disp(probs(2))