diff --git a/triangle_spectrum_plots.m b/triangle_spectrum_plots.m
new file mode 100644
index 0000000000000000000000000000000000000000..2fbbd2afd660a1d7fb4b82a8b5c7240057f93cf4
--- /dev/null
+++ b/triangle_spectrum_plots.m
@@ -0,0 +1,27 @@
+(* ::Package:: *)
+
+gsraw=Import[NotebookDirectory[]<>"cpp/graphdata_spectrum.m"];
+(* gsraw=SortBy[gsraw,{#[[1,1]]&,#[[1,2]]&}]; (* Sort by n and then by tau. The {} forces a *stable* sort because otherwise Mathematica sorts also on triangle count and other things. *) *)
+
+
+gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
+(* Data format: *)
+(* gdata[[ tau index, n index, run index , datatype index ]] *)
+(* datatype index:
+1: {n,tau}
+2: etmt
+*)
+tauvalues=gdata[[All,1,1,1,2]];
+nlabels=Map["n = "<>ToString[#]&,gdata[[1,All,1,1,1]]];
+taulabels=Map["\[Tau] = "<>ToString[#]&,gdata[[All,1,1,1,2]]];
+
+
+ListPlot[gdata[[1,1,1,{4,6}]]]
+ListPlot[gdata[[1,1,1,{5,7}]]]
+
+
+Histogram[gdata[[1,1,1,{4,6}]],50]
+Histogram[gdata[[1,1,1,{5,7}]],50]
+
+
+