gap> L:=[[1,2,3,4,5,6],[8,9,10,11,12,7]];;
gap> M:=[[1,7,8,2],[11,10,4,5]];;
gap> N:=[[2,8,9,3],[12,11,5,6]];;
gap> P:=[[3,9,10,4],[7,12,6,1]];;
gap> Y:=PoincarePrismCWComplex(L,M,N,P);;
gap> IsClosedManifold(Y);
true
gap> StructureDescription(FundamentalGroup(Y));
"C3 : C4"
