用mathematica写一个元胞自动机的程序

下面是一个使用 Mathematica 编写的元胞自动机程序的示例：

(* 定义规则集 *)
rule = {1, 0, 1, 1, 0, 1, 1, 0};

(* 初始化网格 *)
grid = Table[RandomInteger[], {i, 1, 100}, {j, 1, 100}];

(* 进行演化 *)
Do[
  grid = Map[
    (1 &) /@ (rule[[# + 1]] &) /@ (grid[[#]] &) /@ (Range[-1, 1] + # &) /@
        Flatten[Table[{i, j}, {i, 1, 100}, {j, 1, 100}], 1],
    {2}
  ],
  {n, 1, 100}
];

(* 显示网格 *)
ArrayPlot[grid]

上面的程序实现了一个简单的元胞自动机，其中 `rule` 定义了规则集，`grid` 初始化了网格，`Do` 进行了演化，最后使用 `ArrayPlot` 显示网格。