gap> L13_1:=ThreeManifoldViaDehnSurgery([[1,2],[1,2]],13,1);;
gap> DijkgraafWittenInvariant(L13_1,CyclicGroup(13));
[ g1^12, g1^12, g1^10, g1^10, g1^9, g1^9, g1^4, g1^4, g1^3, g1^3, g1, g1, id ]
gap> L13_2:=ThreeManifoldViaDehnSurgery([[1,2],[1,2]],13,2);;
gap> DijkgraafWittenInvariant(L13_2,CyclicGroup(13));
[ g1^11, g1^11, g1^8, g1^8, g1^7, g1^7, g1^6, g1^6, g1^5, g1^5, g1^2, g1^2, 
  id ]
