new method: Region.diff
rm Polytope methods __add__, __sub__, __and__

because they are not closed: they return a Region instead of Polytope in certain cases.
Reflect changes where polytope.polytope is used.