gap> K:=PureCubicalKnot(3,1);
prime knot 1 with 3 crossings

gap> L:=ReflectedCubicalKnot(K);
Reflected( prime knot 1 with 3 crossings )

gap> M:=KnotSum(K,L);
prime knot 1 with 3 crossings + Reflected( prime knot 1 with 3 crossings )

gap> Display(M);
