matlab三门问题

%%2019年8月5号
%%概率论的三羊问题
%%n代表随机次数
n = 100000;
%%不改变注意
nochange = 0;
%%改变注意
change = 0;
for i= 1 : n
%%随机生成一个正数x表示车可能出现在第x扇门后
    x = randi([1,3],1);
%%随机生成一个正数y表示自己选的门
	y = randi([1,3],1);
%%如果刚好选到是车 而且不改变注意 此时能选到车 还是能中奖的
if x == y
      nochange = nochange + 1;
end
%%如果不能选到车  已经知道有一扇门后面是羊 只有改变注意才能中奖 
if x ~= y 
    change = change + 1;
end
end
%%输出不改变注意时的获奖概率
disp(nochange/n);
%%输出已经打开一扇门是山羊 
disp(change/n);