До версии 1.3 этот код работал хорошо:
from sympy import Intersection;
from sympy import solveset;
from sympy import S; from sympy.abc import x;
from sympy.functions.elementary.miscellaneous import Min, Max;
print Intersection([solveset(p, x, S.Reals) for p in [((((x + 2.0000) * 3.0000)+18.000000) > 0.000), ((((x + 2.0000) * 3.0000)+18.000000) < 1.000)]])
Я не знаю, как использовать новый Intersection
, ср https://github.com/sympy/sympy/pull/16344