# 2*x + y >= 0 ; -x - y + 2 - z >= 0 ; z >= 0 ; -x + 1 >= 0 ; -y + 2 >= 0
ineqs = [[0,2,1,0], [2,-1,-1,-1], [0,0,0,1], [1,-1,0,0], [2,0,-1,0]]

P = Polyhedron(ieqs = ineqs)

to_plot = P.projection().render_solid_3d(opacity=.7)
to_plot.show()

#to plot the outline of a 3D polyhedron
to_plot = P.projection().render_wireframe_3d(color = 'blue', alpha = 0.8)
to_plot.show()