read "itall.txt":

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

ver:=vertices(phatlvl):
ver:=subsop(1=NULL,2=NULL,3=NULL,ver):

readFunc(T,12,cat(path,phatlvl,"phatsumabspsi.txt")):
ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"sumabspsi"));
ma:=LocFindFuncMax2(T,phatlvl,ver,12):
mi:=LocFindFuncMin2(T,phatlvl,ver,12):
printf("min=%20.20g at %a\n",mi[1],mi[2]);
printf("max=%20.20g at %a\n",ma[1],ma[2]);