
(* Final all numbers with small Mahler measure, not Salem, *)
(* but still with a root on the real axis *)

<<salem.m

ps=<<small_mahler.m
ps=Select[ps,!Salem[#] && !Salem[#/.x->-x]&]

HasRealRoot[p_]:=Min[Abs[Im[Rts[p]]]] < 0.00001

qs=Select[ps,HasRealRoot]
