read "itall.txt":

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

ver:=subsop(1=NULL,2=NULL,3=NULL,vertices(lvl)):
linever:=our_line_vertices(phatlvl,[2],2,4):
for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"dirichlet.",i)):
	#ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"psi.",i));
	#final_points:=SG_plot_line(T,8,12,linever,10):
	print(i);
	print(ver[i]);	
	pointplot(SG_plot_line(T,8,12,linever,10),color=black,style=LINE);
end do;
print("done printing psis");