- /
-
Cellular automaton, rule 150
on 29 Oct 2021
- 1
- 23
- 0
- 0
- 142
N=bitget(150,1:8);
D=1000;
A=zeros(D,2*D);
A(1,D)=1;
for n=1:D-1
for m=2:(2*D-1)
A(n+1,m)=N( A(n,[m-1 m m+1])*[4 2 1]'+1 );
end
end
imagesc(A);
axis off equal