MATLAB: Monte carlo simulation in matlab two dices roll

homeworkmonte carlo simulation in matlab two dices roll

if 2 dices were thrown & there top value were added ,what is the probability of getting a sum of 7 ?

Best Answer

n        = 10;               % Number or trials
throw    = randi(6, n, 2);   % Value of throws, 2 dice
SumThrow = sum(throw, 2);    % Sum of both dice
Match    = (SumThrow == 7);  % Logical array: 1 if sum of values is 7
count    = sum(Match);       % Number of matching throws

Or with modifications of your code:

count = 0;
for i=1:10
  throws1  = randi(6, 1);
  throws2  = randi(6, 1);
  sumThrow = throws1 + throws2;   % Do not use "sum" as name of a variable
  if sumThrow == 7
    fprintf('sum is 7\n')
    count = count + 1;
  else
    fprintf('false\n')
  end
end
fprintf('Number of throws with sum=7: %d\n', count);

Related Solutions

MATLAB: Write a program to calculate the probability of two fair dice having sum equal to 7. There are two fair dice and we are interested in seeing whether the sum of the two fair dice having sum equal to 7. Use Monte Carlo Simulation.

Everything is fine with your code (although the calculation of n7 and nAll could be done with no loop at all), except for your calculation of prob where for some reason you multiply two probabilities together, the theoretical probabilty (n7/nAll) and your monte-carlo probability (count/ntrials), which of course doesn't make much sense.

MATLAB: Two Dice Monte Carlo Simulation

You don't have your for-loops set up properly. You are trying to find the probability that five consecutive 7's occur in 100 rolls of two dice. So, one trial consists of 100 rolls of two dice. Your Monte Carlo simulation should have an outer loop running the number of trials, and inside this loop is one trial consisting of 100 rolls of the dice where you look for five consecutive 7's. An outline something like this based on your nomenclature:

p=1/6;  % probability of getting a 7 with two dice
ntrials = 100000;  % pick some large number of trials
N=100;  % the number of dice rolls in one trial
sumseven=zeros(ntrials,1);  % changed N to ntrials


for m=1:ntrials  % changed N to ntrials
    % inside this loop is the code for one trial, which is 100 rolls of the dice
    l=0;
    for k=1:N  % changed 5 to N
        % your code here for rolling two dice and seeing if you have 5 consecutive 7's
    end
end
p_sumseven=sum(sumseven(:))/ntrials;  % changed N to ntrials
% etc.

Best Answer

Related Solutions

MATLAB: Write a program to calculate the probability of two fair dice having sum equal to 7. There are two fair dice and we are interested in seeing whether the sum of the two fair dice having sum equal to 7. Use Monte Carlo Simulation.

MATLAB: Two Dice Monte Carlo Simulation

Related Question