I want to see a 1000 results in a Histogram

I have created a script. I wish to loop this a 1000 times and place the result in a Histogram. Below is the script.

k= min(10,randi([1,13]));
a=0;
while k < 17
    r = min(10,randi([1,13]));
if r==1
   a= a+1;
    r=11;
    k=k+r;
    else
    k=k+r;
      end
if k>21 && a > 0
    k = k - 10;
      a=a-1;
end
end
k

I can run the script as a loop, but that does not help me at all with creating a Histogram.

[Adam]

Why don't you just use

doc histogram

?

Answer 1

[Image Analyst]

Try this. Put all the code below (both functions) into one m-file, say, test1.m

function test1
% Initialization / clean-up code.
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
workspace;  % Make sure the workspace panel is showing.
fontSize = 20;
% Call Thomas's code 1000 times, then histogram it.
for kk = 1 : 1000
  k(kk) = ComputeK(kk);
end
histogram(k);
grid on;
title('Histogram of k', 'FontSize', fontSize);
xlabel('k', 'FontSize', fontSize);
ylabel('Counts', 'FontSize', fontSize);

% Some function to do stuff but
% appears to generate k as the output.
function k = ComputeK(k)
k = min(10,randi([1,13]));
a = 0;
while k < 17
  r = min(10, randi([1, 13]));
  if r == 1
    a = a+1;
    r=11;
    k=k+r;
  else
    k=k+r;
  end
  if k > 21 && a > 0
    k = k - 10;
    a = a-1;
  end
end
% k

%

[Thomas Verhaegen]

Hi,

Thanks for the help. I am so much further now. I almost have all things figured out. Now based on your file.

function test1
% Initialization / clean-up code.
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
workspace;  % Make sure the workspace panel is showing.
fontSize = 20;
% Call Thomas's code 1000 times, then histogram it.
for kk = 1 : 1000
  k(kk) = ComputeK(kk);
end
histogram(k);
grid on;
title('Histogram of k', 'FontSize', fontSize);
xlabel('k', 'FontSize', fontSize);
ylabel('Counts', 'FontSize', fontSize);
% Some function to do stuff but
% appears to generate k as the output.
function k = ComputeK(k)
k = I want k to go from 1 to 10;
a = 0;
while k < 17
  r = min(10, randi([1, 13]));
  if r == 1
    a = a+1;
    r=11;
    k=k+r;
  else
    k=k+r;
  end
  if k > 21 && a > 0
    k = k - 10;
    a = a-1;
  end
end
% k

Now I want to run every value from 1 to 10 as k instead of the original random. Now the assigment is to gather how many times I go above 21 for a certain value. by example if I enter k=1 I will get 230/1000 above 21.

Now I need to make a histogram. With all the results with k going from 1 to 10. There lies my problem. I can only make a safe or bust file. as you can see.

function Opdrachtc
clear;
clc;   
fontSize = 15;
for kk = 1 : 1000
  h(kk) = MDealer(kk);
    H=categorical(h,[0,1],{'Safe','Bust'});
end
histogram(H,'FaceColor','k','Edgecolor','r','LineWidth',2);
title('BlackJack', 'FontSize', fontSize);
fontSize=12;
xlabel('Mogelijkheden', 'FontSize', fontSize);
ylabel('Kans op 1000', 'FontSize', fontSize);

      function k= MDealer(~)

%   Berekenen van mogelijke kaarten van de dealer
%   Een kaart met waarde tussen 1 en 10 wordt toegevoegd, zolang de waarde
%   van de kaarten onder de 17 blijft. Zolang men niet boven de 21 komt
%   zal een 1 tellen als een 11. Indien men boven de 21 komt zakt de waarde
%   terug van 11 naar 1.

k=8 ;
a = 0;
if k == 1
    k=11;
    a=a+1;
 end
while k < 17
  r = min(10, randi([1, 13]));
  if r == 1
    a = a+1;
    r=11;
    k=k+r;
  else
    k=k+r;
  end
  if k > 21 && a > 0
    k = k - 10;
    a = a-1;
    end

end
if k>21
    k=1;
else
    k=0;
end

wich is nice but not what I want. I want the x-axes to display 1 - 10 and if able leave out the safe.

[Image Analyst]

If you want k to just do from 1 to 10, use a for loop instead of a while loop:

for k = 1 : 10

[Thomas Verhaegen]

function Opdrachtc
clear;
clc;   
fontSize = 15;
for kk = 1 : 10
  h(kk) = MDealer(kk)
    H=categorical(h,[1,2,3,4,5,6,7,8,9,10]);
end
histogram(H,'FaceColor','k','Edgecolor','r','LineWidth',2);
title('BlackJack', 'FontSize', fontSize);
fontSize=12;
xlabel('Mogelijkheden', 'FontSize', fontSize);
ylabel('Kans op 1000', 'FontSize', fontSize);

      function k= MDealer(~)

%   Berekenen van mogelijke kaarten van de dealer
%   Een kaart met waarde tussen 1 en 10 wordt toegevoegd, zolang de waarde
%   van de kaarten onder de 17 blijft. Zolang men niet boven de 21 komt
%   zal een 1 tellen als een 11. Indien men boven de 21 komt zakt de waarde
%   terug van 11 naar 1.

for l=1:10
k=l ;
a = 0;
if k == 1
    k=11;
    a=a+1;
 end
while k < 17
  r = min(10, randi([1, 13]));
  if r == 1
    a = a+1;
    r=11;
    k=k+r;
  else
    k=k+r;
  end
  if k > 21 && a > 0
    k = k - 10;
    a = a-1;
    end
end
if k>21 && l==1
    k=1;
elseif k>21 &&l==2
    k=2;
elseif k>21 && l==3
    k=3;
elseif k>21 && l==4
    k=4;
    elseif k>21 &&l==5
    k=5;
elseif k>21 && l==6
    k=6;
elseif k>21 && l==7
    k=7;
elseif k>21 &&l==8
    k=8;
elseif k>21 && l==9
    k=9;
elseif k>21 && l==10
    k=10;
else k=0;
end

x(l)=k;
end
k=x;
k