"Source: Text/Geometry/Polyhedra.text";
"Line: 1953";
"Date: Thu Sep 25 14:47:45 2025";
"Main: Fri Sep 26 14:29:19 2025";
// original file: Text/Geometry/Polyhedra.text, line: 1953
// Example: H156E15 ()
print "Example: H156E15";
ei := GetEchoInput();
SetEchoInput(true);
P:=Polytope([
    [ -3, -2, -1, 3 ],
    [ 2, 2, 0, -3 ],
    [ 3, 0, -2, 0 ],
    [ 4, 1, 61, -121 ]
]);
Dimension(P);
assert $1 eq 3;
bool,u,F,phi,v:=IsPyramid(P);
bool;
assert $1;
F;
u;
Image(phi,P) + v;
F:=PolyhedronInSublattice(F);
F;
NumberOfInteriorPoints(F);
assert $1 eq 0;
S:=Polytope([
    [ 0, 0 ],
    [ 1, 0 ],
    [ 0, 1 ]
]);
bool:=IsEquivalent(F,S);
bool;
assert $1;
EhrhartDeltaVector(P);
SetEchoInput(ei);
