read "itall.txt":

PRINTEIGENS:=0:
PRINTPSIS:=0:

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

linever:=our_line_vertices(phatlvl,[2],2,4):

#if (PRINTEIGENS=1) then 
	#for i from 1 to numpsis do
	for i from 8 to 8 do
	        readFunc(T,12,cat(path,phatlvl,"phat",lvl,"eigen.",i)):
		#pointplot(SG_plot_line(T,8,12,linever,Pi),color=blue,style=LINE,title=cat(phatlvl,"phat",lvl,"eigen.",i));
		#display([pointplot(SG_plot_line(T,8,12,linever,Pi),color=blue,style=LINE),plot(sin(x)*findFuncMax(T,8,12),x=0..Pi,color=red)],title=cat(phatlvl,"phat",lvl,"eigen.",i));
		points:=SG_plot_line(T,8,12,linever,Pi);
		ma:=findFuncMax(T,8,12);
		val:=[seq([points[i][1],points[i][2]-sin(points[i][1])*ma],i=1..nops(points))];
		pointplot(val,color=blue,style=LINE);
	 #ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"eigen.",i));
	end do;
	print("done printing eigens");
#end if;



#if(PRINTPSIS=1) then 
#	for i from 1 to numpsis do
#		readFunc(T,12,cat(path,phatlvl,"phat",lvl,"psi.",i));
#		ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"psi.",i));
#	end do;
#	print("done printing psis");
#end if;