next up previous index
Next: Set Facts, Venn Diagrams Up: Sample Space 09/24 Previous: Problem Session,09/25/98: De Méré's

Another Matlab Example for dice game

Suppose we play a dice game where at each throw of a dice I win the number of points if the throw is even and lose the number of points if the throw is odd.

How much will I win or lose over 1000 games?

How much does one win or lose on average?

In order to solve this by simulation we write a separate matlab function by editing a file called dice.m accessible from where you are working from.

%%%%%%%%%%%dice.m file%%%%%

function [win,aver]=dice(B)
%Play the dice game B times
gains=[-1,2,-3,4,-5,6];
plays=unidrnd(6,B,1);
win=sum(gains(plays));
aver=win
%Commands effect:
gains =    -1     2    -3     4    -5     6;
gains([1,6,3,3,2,1]) 
ans =-1     6    -3    -3     2    -1

%Call from matlab:
[w,a]=dice(1000)
w =   563
a =    0.5630
[w,a]=dice(1000)
w =    506
a =    0.5060
[w,a]=dice(10000)
w =5469
a =0.5469
[w,a]=dice(1000)
w =   484
a =    0.4840
[w,a]=dice(10000)
w =        4971
a =    0.4971



Susan Holmes
1998-12-07