read "itall.txt":

maxlvl:=4:
path:="../data/":
lambdatol:=.0000000001:
bmtol:=.0000000001:


lambdalist:=[];

for lvl from 1 to maxlvl do
	lpath:=cat(path,"level",lvl,"/"):
	numvertices:=(3^(lvl+1)-3)/2:
	fd:=fopen(cat(lpath,lvl,"lambdalist.txt"),READ):
	for i from 1 to numvertices do
		lambdalist:=[op(lambdalist),fscanf(fd,"%g\n")[1]];
	end do;
	fclose(fd);
end do;

l:=nops(lambdalist);
i:=1;

while (i < l) do
	badlist:=[];
	for j from (i+1) to nops(lambdalist) do
		if (abs(lambdalist[j]-lambdalist[i]) < lambdatol) then
			badlist:=[op(badlist),j=NULL];
		end if;
	end do;
	lambdalist:=subsop(op(badlist),lambdalist);
	l:=nops(lambdalist);
	i:=i+1;
end do;

fd:=fopen(cat(path,lvl,"bm.txt"),WRITE);

fprintf(fd,"lambda,\t\t\t\tbm\t(level %d)\n",lvl);

for i from 1 to nops(lambdalist) do
	fprintf(fd,"%20.20g\t\t%20.20g\n",lambdalist[i],calcbm(lambdalist[i],[],bmtol));
end do;

fclose(fd);