function ExpectedValueTwoCoins(it);

%function ExpectedValueTwoCoins(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(2);
    r1 = r(1);
    r=randperm(2);
    r2=r(1);
    
    if r1 == 1 & r2 == 1
        value(i) = 0;
    elseif (r1 == 1 & r2 == 2) | (r1 == 2 & r2 == 1)
        value(i) = 1;
    elseif r1 == 2 & r2 == 2
        value(i) = 2;
    end
        
    
    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