电脑
python
整体代码，请下载于：pan.baidu.com/s/1bpbL0nZDQBuICD8PWo26oQ，提取码: 4e43。
元胞自动机需要一个初始状态，这里给出一个初始状态：a=np.zeros((100,136))u,v=a.shapefor i in range(u): for j in range(v): if (i+j)%10<2 or (i-j)%10<2: a[i,j]=1只有0和1的简单情形。本文只考虑这种简单情形。
它的前36步变化的动画如下。
前100步变化。
换一种初始状态：a=np.zeros((100,136))u,v=a.shapefor i in range(u): for j in range(v): if (i+j)%10<2 or (i-j)%10<3: a[i,j]=1前100步动画如下。
a=np.zeros((100,136))u,v=a.shapefor i in range(u): for j in range(v): if (i**2+j**2)%10<2: a[i,j]=1
a=np.zeros((100,136))u,v=a.shapefor i in range(u): for j in range(v): if 20<(i**2+j**2)%36<30: a[i,j]=1
a=np.zeros((100,136))u,v=a.shapefor i in range(u): for j in range(v): if i==68 or j==68: a[i,j]=1