read "itall.txt":
T:=gasket(8):

spike1:=1:
spike2:=16:

readFunc(T,12,cat("../data/8phat4spike.",spike1)):
readFunc(T,13,cat("../data/8phat4spike.",spike2)):
funcDifference(T,6,12,13,14):
funcQuotient(T,6,12,13,15):
ourPlot2(T,6,12,cat("spike ",spike1));
ourPlot2(T,6,13,cat("spike ",spike2));
ourPlot2(T,6,14,cat("spike ",spike1," - spike ",spike2));
ourPlot2(T,6,15,cat("spike ",spike1," / spike ",spike2));