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 k from 1 to (3^(lvl)-3)/2 do 
	print(cat("../data/ortheigen6level",lvl,".",k)):
	T[[3],func]:=readFunc(T,func, cat("../data/ortheigen6level",lvl,".",k)):
	func:=func+1:
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):
#ver:=vertices(lvl):
#for vertex in ver do
#	print(vertex):
#        print(T[vertex,func+2]):
#end do:
        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: