>>> # compute the convex union of two polyhedrons
>>> Polyhedron(square1 | sqaure2)
And(Ge(x, 0), Ge(y, 0), Ge(-y + 3, 0), Ge(-x + 3, 0), Ge(x - y + 2, 0), Ge(-x + y + 2, 0))
>>> # compute the convex union of two polyhedrons
>>> Polyhedron(square1 | sqaure2)
And(Ge(x, 0), Ge(y, 0), Ge(-y + 3, 0), Ge(-x + 3, 0), Ge(x - y + 2, 0), Ge(-x + y + 2, 0))