|
MATHEMATICA
|
lon[n_] := Length@matrices[n]; matrices[n_] := Flatten[Table[{{a, b}, {c, d}}, {a, n}, {b, n}, {c, n}, {d, n}], 3]; com[n_] := Sum[If[Mod[matrices[n][[i]].matrices[n][[j]] - matrices[n][[j]].matrices[n][[i]], n] == {{0, 0}, {0, 0}}, 1, 0], {i, lon[n]}, {j, i, lon[n]}]; Table[2*com[i] - i^4, {i, 4}]
|