Digits:=20;
read "SG_FEM/everything.txt":
read "func_file.txt":
read "funcarith.txt":

for i from 24 to 35 do
	readFunc(T,i,cat("output/psi8.",i));
end do;
for i from 24 to 35 do
	for j from i to 35 do
		funcProduct(T,8,i,j,36);
		err:=simp_rule(T,8,36);
		fd:=fopen("output/psiinnerprod.txt",APPEND);
		fprintf(fd,"%d %d %20.20g\n",i,j,err);
		fclose(fd);
	end do;
end do;