(%i1) load(fourier_elim); (%i2) fourier_elim([x^2 < y , x^2 > y],[x,y]); 2 2 (%o2) [- (y - x ) > 0, y - x > 0] The empty set should be displayed as emptset not [- (y - x ) > 0, y - x > 0]. Is this a bug in fourier_elim or a limitation in the underlying algorithm?