m=50;
X=rand(m,m)>0.5;
ip=[m 1:m-1];
im=[2:m 1];
while any(any(X))
  % somme des huit voisins de X dans le sens trigo
  N=X+X(ip,:)+X(ip,ip)+X(:,ip)+X(im,ip)+X(im,:)+X(im,im)+X(:,im)+X(ip,im);

  % majorité biaisée
  % X vaut 4 ou est strictement plus grand que 5
  X=(N>5)|(N==4);
  %affichage
  pcolor(X)
  axis off image
  shading flat %remove the lattice
  drawnow
end

