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