Add demonstration file
[linpy.git] / demo.py
1 #!/usr/bin/env python3
2
3 from pypol import *
4
5 x, y = symbols('x y')
6
7 sq1 = Le(0, x) & Le(x, 2) & Le(0, y) & Le(y, 2)
8 sq2 = Le(2, x) & Le(x, 4) & Le(2, y) & Le(y, 4)
9
10 print('sq1 =', sq1)
11 print('sq2 =', sq2)
12 print()
13 print('¬sq1 =', ~sq1)
14 print()
15 print('sq1 - sq2 =', sq1 - sq2)
16 print('sq1 - sq2 =', Polyhedron(sq1 - sq2))
17 print()
18 print('sq1 ∩ sq2 =', sq1 & sq2)
19 print('sq1 ∪ sq2 =', sq1 | sq2)
20 print()
21 print('sq1 ⊔ sq2 =', Polyhedron(sq1 | sq2))