clear all
M=zeros(52);
for run=1:7000
n=ceil(2500*rand(1));
V=(rand(1,8)>0.5);
blackNumEdge=0;
whiteNumEdge=0;

for j=1:2:7
    if V(j)==1
        blackNumEdge=blackNumEdge+1;
    end
end


for k=2:8
    if V(k)+V(k-1)==2
        blackNumEdge=blackNumEdge+1;
    elseif V(k)+V(k-1)==0
        whiteNumEdge=whiteNumEdge+1;
    end
end


if V(1)+V(8)==2
    blackNumEdge=blackNumEdge+1;
elseif V(1)+V(8)==0
    whiteNumEdge=whiteNumEdge+1;
end


if blackNumEdge < whiteNumEdge
    V=(V==0);
elseif blackNumEdge == whiteNumEdge
    V=0.*V;
else
    M(ceil(n/50)+1,mod(n,50)+2)=100;
end

M(ceil(n/50),mod(n,50)+2)=V(1);
M(ceil(n/50),mod(n,50)+3)=V(2);
M(ceil(n/50)+1,mod(n,50)+3)=V(3);
M(ceil(n/50)+2,mod(n,50)+3)=V(4);
M(ceil(n/50)+2,mod(n,50)+2)=V(5);
M(ceil(n/50)+2,mod(n,50)+1)=V(6);
M(ceil(n/50)+1,mod(n,50)+1)=V(7);
M(ceil(n/50),mod(n,50)+1)=V(8);

image(100.*M)
title(run)
drawnow
end