Digits:=20:
read "everything.txt":
read "func_file.txt":
read "funcarith.txt":
read "genlevel5.txt":
read "autoset.txt":
read "algo24.txt":


lvl:=4;
path:="../data/";
startEigen:=12;

lasteigen:=makeAllRings(T,lvl,startEigen);
lasteigen:=lasteigen[-1];
makeEdge5(T,lvl,lasteigen,[1,2]);
lasteigen:=lasteigen+1;
makeEdge5(T,lvl,lasteigen,[0,1]);
print("lasteigen= "); print(lasteigen);


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 12 to lasteigen do
	discLaplac(T,lvl,i,lasteigen+1);
	funcQuotient(T,lvl,lasteigen+1,i,lasteigen+1);
	ourPlot(T,lvl,lasteigen+1,"blah");
	m:=findFuncMaxDifferenceWithoutZero(T,lvl,lasteigen+1,.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;