gap> ap:=ArcPresentation(PureCubicalKnot(3,1));;
gap> p:=17;;q:=16;;
gap> W:=ThreeManifoldViaDehnSurgery(ap,p,q);
Regular CW-complex of dimension 3
