gap> p:=Random([2,3,5,7,11,13,17]);;
gap> q:=Random([1..p-1]);;
gap> L:=ThreeManifoldViaDehnSurgery([[1,2],[1,2]],p,q);
Regular CW-complex of dimension 3
