 LIB "polymake.lib";
// the lattice points of the unit square in the plane
list points=intvec(0,0),intvec(0,1),intvec(1,0),intvec(1,1);
// the secondary polytope of this lattice point configuration is computed
intmat secpoly=secondaryPolytope(points)[1];
list np=polymakePolytope(secpoly);
// the vertices of the secondary polytope are:
np[1];
// its dimension is
np[2];
// np[3] contains information how the vertices are connected to each other,
// e.g. the first vertex (number 0) is connected to the second one
np[3][1];
// the affine hull has the equation
ring r=0,x(1..4),dp;
matrix M[5][1]=1,x(1),x(2),x(3),x(4);
intmat(np[4])*M;
