read "itall.txt":

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

setoptions3d(shading=ZHUE);

fd:=fopen(cat(path,lvl,"limlambdalist.txt"),READ);

#l:=.4;
for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"eigen.",i));
	plotsetup(gif,plotoutput=cat("/tmp/richard/epics/eigen",i,".gif"));
	#ourPlotFancy(T,phatlvl,12,cat(phatlvl,"phat",lvl,"psi.",i));
	l:=fscanf(fd,"%g\n")[1];
	ourPlotFancy(T,phatlvl,12,sprintf("p=%g",l));
end do;
fclose(fd);