diff --git a/triangle_spectrum_plots.m b/triangle_spectrum_plots.m
index e525f1aa7a49ff7c4cda6e15d7c26d53a00a96b5..7b0ed5b396fe37e67b9b6249c597017a008cbb91 100644
--- a/triangle_spectrum_plots.m
+++ b/triangle_spectrum_plots.m
@@ -21,8 +21,28 @@ nlabels=Map["n = "<>ToString[#]&,gdata[[1,All,1,1,1]]];
 taulabels=Map["\[Tau] = "<>ToString[#]&,gdata[[All,1,1,1,2]]];
 
 
-Histogram[gdata[[1,1,All,5,3,1]]]
-Histogram[gdata[[1,1,All,5,3,2]]]
+histograms=Map[Histogram[#[[All,5,2,1]],{0,#[[1,1,1]],Automatic}]&,gdata,{2}];
+
+
+TableForm[histograms,TableHeadings->{taulabels,nlabels}]
+
+
+initialLaplacianData=gdata[[All,All,All,5,2]];
+mixedLaplacianData=gdata[[All,All,All,5,4]];
+lmaxDifference=Map[#[[5,2,1]]-#[[5,4,1]]&,gdata,{3}];
+lmaxRatio=Map[#[[5,2,1]]/#[[5,4,1]]&,gdata,{3}];
+lgapDifference=Map[#[[5,2,2]]-#[[5,4,2]]&,gdata,{3}];
+lgapRatio=Map[#[[5,2,2]]/#[[5,4,2]]&,gdata,{3}];
+lratioDifference=Map[#[[5,2,3]]-#[[5,4,3]]&,gdata,{3}];
+lratioRatio=Map[#[[5,2,3]]/#[[5,4,3]]&,gdata,{3}];
+
+
+Histogram[lmaxDifference[[1,1,All]]]
+Histogram[lmaxRatio[[1,1,All]]]
+Histogram[lgapDifference[[1,1,All]]]
+Histogram[lgapRatio[[1,1,All]]]
+Histogram[lratioDifference[[1,1,All]]]
+Histogram[lratioRatio[[1,1,All]]]