changed to:
solve(sin(x)<1, assume=Q.real(x)) ==  Ne(x,pi/2)