file:=HapFile("data253a.txt");;
Read(file);
G:=SymmetricMatrixToFilteredGraph(D,20);
K:=FilteredRegularCWComplex(CliqueComplex(G,2));
P:=PersistentBettiNumbers(K,0);;
BarCodeCompactDisplay(P);
P:=PersistentBettiNumbers(K,1);;
BarCodeCompactDisplay(P);

Y:=FiltrationTerm(K,12);;
Display(HomotopyGraph(Y));;
