fact: Definition root_seq : {poly F} -> seq F := let fix loop (p : {poly F}) (n : nat) := if n is n.+1 then if (size p != 1%N) =P true is ReflectT p_neq0 then let x := projT1 (sigW (closed_rootP p p_neq0)) in x :: loop (p %/ ('X - x%:P)) n else [::] else [::] in fun p => loop p (size p).
imports: mathcomp(all_ssreflect) mathcomp(all_algebra) mathcomp(all_real_closed) CoqEAL(ssrcomplements)