状态：=选择（proc（x）（x[1]=x[2]或x[5]=1）和（x[2]=x[3]或x[6]=1）和（x[3]=x[4]或x[7]=1）结束proc，[seq（seq（seq（seq（seq（seq（seq（[a，b，c，d，e，f，g]，g=0..1），f=0..1），e=0..1），d=0..1），c=0..1），b=0..1），a=0..1）]）：
T： =矩阵（54，54，proc（i，j）局部k；
如果add（states[j，k]-states[i，k]，k=1..4）>0，则返回0fi；
如果状态[j，5]>状态[i，5]或状态[j、6]>状态[i，6]或状态[j，7]>状态[1，7]，则返回0 fi；
如果状态[i，1]>=状态[i，2]和状态[j，5]<>状态[i，5]，则返回0fi；
如果状态[i，2]>=状态[i、3]和状态[j，6]<>状态[i和6]，则返回0 fi；
如果状态[i，3]>=状态[i、4]和状态[j，7]<>状态[i和7]，则返回0 fi；
1
结束进程）：
U： =矢量（54，1）：
E[0]：=向量（54）：E[0][1]：=1：
对于从1到25的k，做E[k]：=T。E[k-1]od（E[k-1]od）：
seq（U^%T.E[j]，j=1..25）； #罗伯特·伊斯雷尔2019年9月8日
