improved_regions.sagews

1+1

2

P=Polyhedron(vertices=((129/10, 3/20, -3),(123/10, 3/20, -3), (12, 3/16, -57/20),(117/10, 3/20, -27/10)),base_ring=QQ)

Q=Polyhedron(vertices=[[1033/170, -9/170, 19/17], [61/10, -1/20, 11/10], [91/20, -1/20, 11/10], [771/170, -9/170, 19/17]],rays=[[-3, 0, 1]],base_ring=QQ)

P.intersection(Q)

The empty polyhedron in QQ^3

P.plot()+Q.plot()

3D rendering not yet implemented