Base Polynomial found using census program: ( also in the Salem list as #23): Clear[h, k, a, b, c, d, e, f, g, h, j, k, l, m, p, x]; h[x_] = a + b*x + c*x^2 + d*x^3 + e*x^4 + f*x^5 + g*x^6 + j*x^7 + x^8; n0 = 8; Union[Flatten[Table[If [Abs[x /. NSolve[h[x] == 0, x][[n0]]] - Max[Abs[Table[ x /. NSolve[h[x] == 0, x][[n]], {n,1, n0}]]] == 0 && Im[x /. NSolve[h[x] == 0, x][[n0]]] ==0 && Sort[Abs[Table[x /.NSolve[h[x] == 0, x][[n]], {n, 1, n0}]]][[n0 - 1]] <= 1 &&Sort[Abs[Table[x /. NSolve[h[x] == 0,x][[n]], {n, 1, n0}]]][[1]] > 0, {Sort[Abs[Table[x /.NSolve[h[x] == 0, x][[n]], {n, 1, n0}]]], h[x]}, {}], {a, -1,1}, {b, -1, 1}, {c, -1, 1}, {d, -1, 1}, {e, -1, 1}, {f, -1,1}, {g, -1, 1}, {j, -1, 1}], 7]]