function p = BR_MC(num_sims,num_trials)
%

p = zeros(num_sims,1);
for s = 1:num_sims
    s
    T_A = rand(num_trials,1)<=0.95;
    T_B = rand(num_trials,1)<=0.75;
    T_C = rand(num_trials,1)<=0.7512;
    indexes = find(T_A==0|T_B==0|T_C==1);
    p(s) = length(indexes)/num_trials;
end
