Digits:=20:
read "everything.txt":
read "func_file.txt":
read "funcarith.txt":

lvl:=4;
path:="../data/";
numEigens:=120;
startEigen:=12;
for i from 1 to numEigens do
	readFunc(T,i+startEigen-1,cat(path,"n4eigen.",i));
end do;
if 0=1
then
for i from 1 to numEigens do
	for j from i to numEigens do
		m:=discInnerProduct(T,lvl,i+startEigen-1,j+startEigen-1);
                #fd:=fopen(cat(path,"neigensvalues.txt"),APPEND);
                printf("(%d,%d)=%20.20g\n",i,j,m);
                #fclose(fd);
	end do;
end do;
end if;
for i from 1 to numEigens do
	discLaplac(T,lvl,i+startEigen-1,startEigen+numEigens);
	funcQuotient(T,lvl,startEigen+numEigens,i+startEigen-1,startEigen+numEigens);
	m:=findFuncMaxDifferenceWithoutZero(T,lvl,startEigen+numEigens,.01);
        #fd:=fopen(cat(path,"neigensvalues.txt"),APPEND);
        printf("max variance of (laplacian(%d)/%d) =%20.20g\n",i,i,m);
        #fclose(fd);
end do;