function sample2dice(n)
%function sample2dice(n)
%n = number of samples.

X=zeros(1,n);
Y=zeros(1,n);
EX=zeros(1,n);
EY=zeros(1,n);
EXY=zeros(1,n);

SX=zeros(1,n);
SY=zeros(1,n);
SXY=zeros(1,n);

for i=1:n
    
    r1 = randperm(6);
    r2 = randperm(6);
    X(i) = r1(1);
    Y(i) = r2(1);
    
    EX(i) = sum(X(1:i))/i;
    EY(i) = sum(Y(1:i))/i;
    SX(i) = sqrt(sum((X(1:i)-EX(i)).^2)/i);
    SY(i) = sqrt(sum((Y(1:i)-EY(i)).^2)/i);
    XY(i) = X(i)+Y(i);
    EXY(i) = sum(XY(1:i))/i;
    SXY(i) = sqrt(sum((XY(1:i)-EXY(i)).^2)/i);
    
    figure(1);clf;
     set(gcf,'doublebuffer','on');
    subplot(1,3,1);
    hist(X(1:i));
    subplot(1,3,2);
    hist(Y(1:i));
    subplot(1,3,3);
    hist(XY(1:i));
    drawnow;
    
    figure(2);clf;
     set(gcf,'doublebuffer','on');
    colordef black;
    plot(1:i,EX(1:i),'r-','linewidth',2);hold on;
    plot(1:i,EY(1:i),'g-','linewidth',2);hold on;
    plot(1:i,EXY(1:i),'y-','linewidth',2);hold on;
    plot(1:i,SX(1:i),'r--','linewidth',2);hold on;
    plot(1:i,SY(1:i),'g--','linewidth',2);hold on;
    plot(1:i,SXY(1:i),'y--','linewidth',2);hold on;
    drawnow;
    
end