> For an example, construct_domain([Integral(x/2, (x, 0, 2))],composite=None) gives (EX, [EX(Integral(x/2, (x, 0, 2)))]) 

I don't know sympy.polys well, can you explain why this is "wrong"?