file:=HapFile("data1V2X.pdb");;
K:=ReadPDBfileAsPureCubicalComplex(file);;
C:=PureComplexComplement(K);;
Print(Dimensions(C),"\n");
A:=PureComplexSubcomplex(C,[[1..76],[1..77],[1..70]]);;
B:=PureComplexSubcomplex(C,[[1..76],[1..77],[70..121]]);;
AB:=PureComplexIntersection(A,B);;
Print(BettiNumber(AB,0),"\n");
AB:=ContractedComplex(AB);;
A:=ContractedComplex(A,AB);;
B:=ContractedComplex(B,AB);;
f:=FundamentalGroup(RegularCWMap(A,AB));
Print(f,"\n");
g:=FundamentalGroup(RegularCWMap(B,AB));
Print(g,"\n");
F:=Pushout(f,g);;
F:=SimplifiedFpGroup(F);
Print(F,"\n");
Print(RelatorsOfFpGroup(F),"\n"); 
