gap> F:=ThickeningFiltration(K,15);;
gap> F:=FilteredPureCubicalComplexToCubicalComplex(F);;
gap> F:=FilteredCubicalComplexToFilteredRegularCWComplex(F);;
gap> P:=PersistentBettiNumbersAlt(F,1);;
gap> BarCodeCompactDisplay(P);
