|
%CA setupn=200;%initialize the arraysz=zeros(n,n);cells=z;sum=z;%set a few cells to onecells(100,100)=1;%index definition for cell updatex=2:1:199;y=2:1:199;((x.^2)+(y.^2))<=100^2;%nearest neighbor sumsum(x,y)=cells(x,y+1)+cells(x,y-1)+cells(x+1,y)+cells(x-1,y);% The CA rulecells=(sum==1);plot(x,y |
