Don't panic!* Morten Dam Jørgensen - Computerfysik.dk

n = 100; % grid size (n x n)
X = randi(2,n,n)-1; % MAGIC GOES HERE!!!!
Y = X; colormap bone
for t=1:1000000 % number of steps
image(Y) % Tegn gitter
for x=1:size(X,1)
for y=1:size(X,2)
if x == 1
if y == 1
s = sum(sum(X(x:x+1, y:y+1))) - X(x,y);
B
elseif y == size(X,2)
s = sum(sum(X(x:x+1, y-1:y))) - X(x,y);
else
E
s = sum(sum(X(x:x+1, y-1:y+1))) - X(x,y);
end
elseif x == size(X,1)
K
if y == 1
s = sum(sum(X(x-1:x, y:y+1))) - X(x,y);
elseif y == size(X,2)
L
s = sum(sum(X(x-1:x, y-1:y))) - X(x,y);
else
s = sum(sum(X(x-1:x, y-1:y+1))) - X(x,y);
A
end
else
if y == 1
s = sum(sum(X(x-1:x+1, y:y+1))) - X(x,y);
G
elseif y == size(X,2)
s = sum(sum(X(x-1:x+1, y-1:y))) - X(x,y);
else
E
s = sum(sum(X(x-1:x+1, y-1:y+1))) - X(x,y);
end
end
R
if X(x,y) == 1 && s < 2 % Too few neighbours, underpopulation
Y(x,y) = 0;
elseif X(x,y) == 1 && (s == 2 || s == 3) % balance
Y(x,y) = X(x,y);
elseif X(x,y) == 1 && s > 3 % over-population
Y(x,y) = 0;
elseif X(x,y) == 0 && s == 3 % Spawn new cell
Y(x,y) = 1;
end
end
end
pause(0.01)
X = Y;
end