gap> W:=PureComplexComplement(FiltrationTerm(T,25));;
gap> L:=[];;
gap> for i in [1..PathComponentOfPureComplex(W,0)] do
gap> P:=PathComponentOfPureComplex(W,i);;
gap> Q:=ThickeningFiltration(P,50,M);;
gap> Add(L,Q);;
gap> od;;

gap> B:=Basins(L);
gap> Display(B);
