Revision history [back]
 | 1 | initial version |  |
Like this, maybe:
sage: x, y = polygens(QQ, 'x,y')
sage: R = x.parent()
sage: f = x^3*y+2*x^2*y^2+x*y^3+y^7+x^5
sage: P = f.newton_polytope()
sage: [R({exp:cf for exp,cf in f.dict().items() if exp in facet}) for facet in P.facets()]
[x^5 + x^3*y, y^7 + x^5, y^7 + x*y^3, x^3*y + 2*x^2*y^2 + x*y^3]