read "itall.txt":

phatlvl:=[5,8,8,8,8];

vertex:=[0,2];
startlvl:=1;
endlvl:=5;

T:=gasket(8);

#readFunc(T,12,"../data/level1/5phat1psi.2");
#copyRotateScaleChunk(T,1,sphatlvl,12,13,vertex);
#ourPlot(T,6,13,"blah");

spikelvl:=[];

for j from startlvl to endlvl do
	readFunc(T,12 -1 +j,cat("../data/level",j,"/",phatlvl[j],"phat",j,"spike.2"));	
	spikelvl:=[op(spikelvl),(phatlvl[j]+2-j)];
	ourPlot2(T,spikelvl[j],12 -1 +j,cat(phatlvl[j],"phat",j,"spike.2"));
end do;

for j from startlvl to endlvl do
	funcDifference(T,spikelvl[endlvl],12 -1 +endlvl,12 -1 +j,11);
	ourPlot2(T,spikelvl[endlvl],11,cat("boundary function ",j));
	print(cat("Func Max is ",findFuncMax(T,spikelvl[endlvl],11)));
	print(cat("Func Min is ",findFuncMin(T,spikelvl[endlvl],11)));

end do;