3d Plot working
[linpy.git] / examples / diamond.py
1 #!/usr/bin/env python3
2
3 from pypol import *
4
5 x, y, z = symbols('x y z')
6 diam = Ge(y, x - 1) & Le(y, x + 1) & Ge(y, -x - 1) & Le(y, -x + 1)
7 print('diamond:', diam)
8 print()
9 rhom1 = Le(0, x) & Le(x, 3) & Le(0, y) & Le(y, 3) & Le(0, z) & Le(z, 3) \
10 & Le(z - 2, x) & Ge(z + 2, x) & Ge(z - 1, -x) & Le(z - 5, -x) \
11 & Le(z - 2, y) & Ge(z + 2, y) & Ge(z - 1, -y) & Le(z - 5, -y) \
12 & Le(y - 2, x) & Ge(y + 2, x) & Ge(y - 1, -x) & Le(y - 5, -x)
13 rhom1.plot()
14 rhom2 = rhom1 & Le(x + y + z, 7) & Ge(-2, -x - y - z ) \
15 & Le(x + y - z, 4) & Ge(x + y - z, -1) \
16 & Le(x - y + z, 4) & Ge(x - y + z, -1) \
17 & Le(-x + y + z, 4) & Ge(-x + y + z, -1)
18 rhom2.plot()