|
#This example takes quite some time!
n:=50;; Betti2:=[];; Betti3:=[];;
for p in [10..133]*(1/200) do
Print("p= ",p,"\n");
Y:=RandomSimplicialTwoComplex(n,p);;
K:=CliqueComplex(Y,4);;
W:=RegularCWComplex(K);;
W:=ContractedComplex(W);;
b2:=BettiNumber(W,2);;
Add(Betti2,[p,b2]);;
b3:=BettiNumber(W,3);;
Add(Betti3,[p,b3]);;
od;
ScatterPlot(Betti2);
ScatterPlot(Betti3);
[ Dauer der Verarbeitung: 0.14 Sekunden
(vorverarbeitet)
]
|
