read "itall.txt":

lvl:=1:
phatlvl:=8:
numpsis:=(3^(lvl+1)-3)/2:
path:=cat("../data/","level",lvl,"/"):
T:=gasket(phatlvl):

vl:=ourVertices(lvl):
vl:=subsop(1=NULL,2=NULL,3=NULL,vl):



for j from 1 to nops(vl) do
ser:=0;
fd:=fopen(cat(path,lvl,"lambdalist.txt"), READ);
for i from 1 to numpsis do
	ulambda:=fscanf(fd,"%g\n")[1];
	print(i);print(vl[j]);
	print("ulambda=");
	print(ulambda);
	print("bm=");
	print(calcbm(ulambda,[],.000000000001));
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"eigen.",i));
	print("val=");
	print(T[vl[j],12]^2);
	ser:=ser+T[vl[j],12]^2*(3^(lvl+1))/(2*calcbm(ulambda,[],.000000000001));	
	print("sum");print(i);print(ser);
end do;
print("sum=");
fclose(fd);
print(ser);
end do;