- /
-
cellular automaton
on 4 Oct 2021
- 1
- 28
- 0
- 0
- 171
rng(19);
N=60;
M=rand(N)>0.5;A=[1,1,1;1,0,1;1,1,1;];
for t=1:100
B=conv2(M,A);
S=B(2:N+1,2:N+1);
M = (S==3)+M.*(S==2);
end
colormap gray
imagesc(M);
pbaspect([1,1,1]);
ax=gca;
ax.Visible='off';