read "itall.txt":

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

fd:=fopen(cat(path,phatlvl,"phat",lvl,"dirintegrals.txt"),WRITE);
fprintf(fd,"l1, l2, and integral data for the %d phat %d dirichelt\n",phatlvl,lvl);
fclose(fd);

for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"dirichlet.",i));
	funcProduct(T,phatlvl,12,12,13);
	l2:=simp_rule(T,phatlvl,13);
	funcAbs(T,phatlvl,12,13);
	l1:=simp_rule(T,phatlvl,13);
	integ:=simp_rule(T,phatlvl,12);
	fd:=fopen(cat(path,phatlvl,"phat",lvl,"dirintegrals.txt"),APPEND);
	fprintf(fd,"%d:L1=%20.20g\tL2=%20.20g\tintegral=%20.20g\n",i,l1,l2,integ);
	fclose(fd);
end do;