Use ZZ/(2)[x[1..3]];
M:=[[x[3]^3], [x[1], x[2]], [x[2]^2], [x[1]^2], [x[2], x[3], x[1]],
[x[1], x[3], x[1]], [x[1], x[3]^2], [x[2], x[3], x[2], x[3]], [x[1],
x[3], x[2], x[3]]];
NC.AdMatrix(M);
[[[x[3], x[2], x[1]], [x[3]^2, x[1]], [x[1], x[3], x[2]], [x[2], x[3], x[2]],
[x[3]^2, x[2]], [x[2], x[1], x[3]], [x[3], x[1], x[3]], [x[3], x[2], x[3]],
[x[2], x[3]^2]], Mat([
[0, 0, 0, 0, 0, 1, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 1, 0, 0],
[1, 0, 0, 0, 0, 0, 0, 0, 0],
[1, 0, 0, 0, 0, 0, 0, 0, 0],
[1, 0, 0, 0, 0, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 0, 0, 0, 0],
[0, 0, 1, 0, 0, 0, 0, 0, 0],
[0, 0, 0, 1, 0, 0, 0, 0, 1],
[0, 1, 0, 0, 1, 0, 0, 0, 0]
])]
-------------------------------
|