read "SG_FEM/everything.txt":
read "func_file.txt":
read "funcarith.txt":
read "algo24.txt":

calclevel:=6;
T:=gasket(calclevel);
for i from 12 to 23 do
	readFunc(T,i,cat("output/psi8.",i));
end do;
for i from 12 to 23 do
	discLaplac(T,calclevel,i,i+12);
end do;
for i from 12 to 23 do
	funcQuotient(T,calclevel,i+12,i,i+24);
end do;
for i from 12 to 23 do
	ourPlot(T,calclevel,i+24,cat("eigenvalue",i));
end do;