file:=HapFile("25spheres.txt");;
Read(file);
K:=PureCubicalComplexToCubicalComplex(M);;
Print(Size(K),"\n");
L:=DVFReducedCubicalComplex(K);;
D:=ChainComplex(L);;
Print(List([0..3],D!.dimension),"\n");
