On Tue, 11 Dec 2012 19:14:38 +0100, Jaime Villate <villate at fe.up.pt> wrote: > Is there a way to restrict solve or to_poly_solve to only real solutions? realp(x) := if imagpart(x) = 0 then true else false$ solution: solve(...)$ sublist(solution, realp) Regards Enrique