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.
monte carlo simulatio
Related Solutions
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') endendfprintf('Number of throws with sum=7: %d\n', count);Nice try, but why not just use the sum() function - that's what it's there for:
% Roll the dice "numberOfRolls" times
numberOfRolls = 200; % Number of times you roll all 6 dice.
n = 10; % Number of dice.
maxFaceValue = 6;rolls = randi(maxFaceValue, n, numberOfRolls)% Sum up dice values for each roll.
columnSums = sum(rolls, 1)% Find out how many times each sum occurred.
edges = min(columnSums):max(columnSums)counts = histc(columnSums, edges)% Normalize
grandTotalSum = sum(counts(:))normalizedCountSums = counts / grandTotalSumbar(edges, normalizedCountSums, 'BarWidth', 1);grid on;% Enlarge figure to full screen.
set(gcf, 'units','normalized','outerposition',[0 0 1 1]);% Give a name to the title bar.
set(gcf,'name','Demo by ImageAnalyst','numbertitle','off') title('Frequency of Roll Sums', 'FontSize', 40);
Best Answer