MontyHall problemi

Herkese iyi günler dilerim
Yapmak istediğim deneyin açıklama videosunu aşağıya ekledim : 
https://www.youtube.com/watch?v=vZzFjyjJUt8
Bu deneyi oluştururken kullandığım kodları sorunun en altına ekliyorum.
Benim kodumda her tekrarda kazanma olasılğını hesaplamak yerine bana çıktı olarak genel olasılığı veriyor. Bunun yerine 100000 tekrar yaptırmak ve her tekrarda bulduğum olasılıkları bir vektöre aktarmak istiyorum. Ve bunun yanında olasılık tekrar grafiğini çizmek istiyorum.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function montyHall(kapiSayisi,tekrarSayisi)

assert(kapiSayisi > 2);

function num = randInt(n)
num = floor( n*rand()+1 );
end

%The first column will tallie wins, the second losses
kapiyiDegis = [0 0];
ayniKapi = [0 0];

for i = (1:tekrarSayisi)

secilebilirKapi = (1:kapiSayisi); %Preallocate the available doors
arabaKapi = randInt(kapiSayisi); %Define the winning door
ilkSecim = randInt(kapiSayisi); %The player picks his initial choice

secilebilirKapi(ilkSecim) = []; %Remove the players choice from the available doors

%Pick the door to open from the available doors
acKapi = secilebilirKapi(randperm(numel(secilebilirKapi))); %Sort the available doors randomly
acKapi(acKapi == arabaKapi) = []; %Make sure Monty doesn't open the winning door
acKapi = acKapi(randInt(numel(acKapi))); %Choose a random door to open

secilebilirKapi(secilebilirKapi==acKapi) = []; %Remove the open door from the available doors
secilebilirKapi(end+1) = ilkSecim; %Put the player's original choice back into the pool of available doors
secilebilirKapi = sort(secilebilirKapi);

yeniSecim = secilebilirKapi(randInt(numel(secilebilirKapi))); %Pick one of the available doors

if yeniSecim == ilkSecim
switch yeniSecim == arabaKapi
case true
ayniKapi(1) = ayniKapi(1) + 1;
case false
ayniKapi(2) = ayniKapi(2) + 1;
otherwise
error 'ERROR'
end
else
switch yeniSecim == arabaKapi
case true
kapiyiDegis(1) = kapiyiDegis(1) + 1;
case false
kapiyiDegis(2) = kapiyiDegis(2) + 1;
otherwise
error 'ERROR'
end
end
end

disp(sprintf('İlk seçim değiştirilirse kazanma yüzdedi: %f%%\nİlk seçimde kalınırsa kazanma yüzdesi: %f%%\n', [kapiyiDegis(1)/sum(kapiyiDegis),ayniKapi(1)/sum(ayniKapi)] * 100));

end