Hi,

I am trying to construct a semialgebraic set using package SumOfSquares, which involves a mixture of equality and inequality that would be easier to describe by a for loop. Is there an easy way for me to take the intersection of semialgebraic sets and algebraic sets?

For example, I’d like to achieve the following

```
using JuMP
using DynamicPolynomials
using SumOfSquares
@polyvar x y
S1 = @set x^2 - x == 0
S2 = @set x + y <= 1
intersect!(S1, S2)
```

Many thanks in advance!