gap> p:=2;;k:=2;;
gap> P:=PersistentHomologyOfFilteredChainComplex(C,k,p);;
gap> BarCodeDisplay(P);
