read "itall.txt":

lvl:=3:
phatlvl:=8:
numpsi:=(3^(lvl+1)-3)/2:
startFunc:=12:
path:="../data/":

for i from 1 to numpsi do
	readFunc(T,i+startFunc-1,cat(path,phatlvl,"phat",lvl,"psi.",i));
end do;
for i from 1 to numpsi do
	for j from i to numpsi do
		funcProduct(T,phatlvl,i+startFunc-1,j+startFunc-1,startFunc+numpsi);
		err:=simp_rule(T,phatlvl,startFunc+numpsi);
		fd:=fopen(cat(path,phatlvl,"phat",lvl,"psiinner.txt"),APPEND);
		fprintf(fd,"%d %d %20.20g\n",i,j,err);
		fclose(fd);
	end do;
end do;