function [P,S] = BR_rand_MP(num_samps)
%

samps = zeros(num_samps,2);

for s = 1:num_samps
    r = rand;
    if r<0.15
        samps(s,:) = [0,0];
    elseif r<0.3
        samps(s,:) = [0,1];
    elseif r<0.6
        samps(s,:) = [1,0];
    else
        samps(s,:) = [1,1];
    end
end

num_P = length(find(samps(:,1)==1));
num_nPvQ = length(find(samps(:,1)==0|samps(:,2)==1));
P(1) = length(find(samps(:,1)==1));
P(2) = length(find(samps(:,2)==1));
P = P/num_samps;
S = [num_P;num_nPvQ]/num_samps;
