polymake.lib
LIB "polymake.lib"; ring r=0,(x,y,z),dp; poly f=y3+x2+xy+2xz+yz+z2+1; // the lattice points of the Newton polytope of f are newtonPolytopeLP(f);