Read(“AS.txt”);
a0 := IdentityMat(Rationals, 756);;
allones := 0*a0+1;;
a0+a1+a2+a3+a4 = allones;
pairs := Combinations([a1,a2,a3,a4],2);;
ForAll(pairs, t -> t[1]*t[2] = t[2]*t[1]);
vec := VectorSpace(Rationals,[a0,a1,a2,a3,a4]);
ForAll(pairs, t -> t[1]*t[2] in vec);