function ExpectedValueDie(it)
%function ExpectedValueDie(it)
%
% it = number of iterations (default = 1000)

if nargin < 1
    it = 1000;
end

value = zeros(1,it);
EX = zeros(1,it);

for i=1:it
    
    r=randperm(6);
    
    value(i) = r(1);
    
    EX(i) = sum(value(1:i))/i;
    
    figure(1);clf;
    colordef black;
    set(gcf,'doublebuffer','on');
    plot(1:i,EX(1:i),'r-','linewidth',2);
    drawnow;
    
end