application 'polytope';
$Q =cube(2,1);
$Q2=cube(2,2);
$QxQ2=product($Q,$Q2);
$Q2xQ=product($Q2,$Q);
map { print $_->VERTICES->[0], "\n" } ($QxQ2, $Q2xQ);
$P=conv($QxQ2,$Q2xQ); 
print $P->F_VECTOR;
$P->SCHLEGEL(FACET=>3);
$C5=cube(5);
print isomorphic($P->GRAPH->ADJACENCY,$C5->GRAPH->ADJACENCY);