Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

For an answer of the interpretation of eqns see your first post.

in order to calculate the intersection of your two planes and the cube, you can simply put all conditions in a new Polyhedron

cube = polytopes.n_cube(3)
plane1 = Polyhedron(eqns=[(0,1,0,0)])
plane2 = Polyhedron(eqns=[(1,0,0,-1)])
intersec=Polyhedron(eqns=plane1.equations()+plane2.equations(), ieqs=cube.inequalities())
print intersec.Hrepresentation()