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

func:=12;
 
lvl:=2;


for k from 1 to (3^(lvl-1)+3)/2 do 		
		print(cat("../data/ortheigen5level",lvl,".",k));
		T[[3],func]:=readFunc(T,func, cat("../data/ortheigen5level",lvl,".",k));
		func:=func+1;
end do;

#for j from 2 to 5 do 
#        for k from 1 to (3^(j)-3)/2 do 
#                T[[3],func]:=readFunc(T,func, cat("../data/ortheigen6level",j,".",k));
#                func:=func+1;
#        end do;
#end do;

for i from 12 to (func -1) do

        discLaplac(T,T[[3],i],i,func+2);
        funcQuotient(T,T[[3],i],func+2,i,func+2);
        m:=findFuncMaxDifferenceWithoutZero(T,T[[3],i],func+2,.00000001);
        #fd:=fopen(cat(path,"neigensvalues.txt"),APPEND);
        printf("max variance of (laplacian(%d)/%d) =%20.20g\n",i,i,m);
        #fclose(fd);


end do;