Z := IntegerRing();
n := 11;
P<x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11> := PolynomialRing(Z,n);
X := [x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11];
A := Matrix(P,n,n,[0 : i in [1..n^2]]);
A;
for i in [1..n] do
    for j in [1..n] do A[i,j] := X[i]^(j-1); end for;
end for;
A;
time D := Determinant(A);
Length(D);
time F := Factorization(D);
F;