\[GothicCapitalR] = {{{1, 0, 0}, {0, 1, 0}, {0, 0, 1}}, {{0, 1,
0}, {0, 0, 1}, {1, 0, 0}}, {{0, 0, 1}, {1, 0, 0}, {0, 1,
0}}, {{1, 0, 0}, {0, 0, 1}, {0, 1, 0}}, {{0, 0, 1}, {0, 1,
0}, {1, 0, 0}}, {{0, 1, 0}, {1, 0, 0}, {0, 0, 1}}};
i = 1;
j = 1;
det = 1;
a = Subsets[Range[6], {3}];
v = {x, y, z};
k = \[GothicCapitalR].v;
r = k[[a[[i]]]];
Det[r];
While[i < 21,
r = k[[a[[i]]]];
det = Factor[Det[r]];
Print[det]; i++]
-(x+y+z) (x^2-x y+y^2-x z-y z+z^2)
(x-y) (y-z) (x+y+z)
-(x-z) (y-z) (x+y+z)
-(x-y) (x-z) (x+y+z)
(x-z) (y-z) (x+y+z)
(x-y) (x-z) (x+y+z)
-(x-y) (y-z) (x+y+z)
-(x-z) (y-z) (x+y+z)
-(x-y) (y-z) (x+y+z)
-(x-y) (x-z) (x+y+z)
-(x-y) (x-z) (x+y+z)
(x-y) (y-z) (x+y+z)
-(x-z) (y-z) (x+y+z)
(x-y) (y-z) (x+y+z)
(x-y) (x-z) (x+y+z)
-(x-z) (y-z) (x+y+z)
-(x-y) (x-z) (x+y+z)
(x-z) (y-z) (x+y+z)
(x-y) (y-z) (x+y+z)
-(x+y+z) (x^2-x y+y^2-x z-y z+z^2)
Anyone have any suggestions on how to get the output of this loop to be a list of polynomials so that I can use the PolynomialLCM
function on it? I will be adapting this loop to larger dimensions and will need to be able to have Mathematica perform this calculation for me.