diff --git a/triangle_canonical_timeevol.m b/triangle_canonical_timeevol.m
new file mode 100644
index 0000000000000000000000000000000000000000..42bc8457f9a4c824dd0a8ff169c63b7a9b750c5d
--- /dev/null
+++ b/triangle_canonical_timeevol.m
@@ -0,0 +1,22 @@
+(* ::Package:: *)
+
+gsraw=Import[NotebookDirectory[]<>"data/graphdata_canonical_timeevol.m"];
+gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
+(* Data format: *)
+(* gdata[[ tau index, n index, run index, datatype index ]] *)
+(* datatype index:
+1: {n,tau}
+2: { triangles }
+*)
+
+
+tauChoice=3;
+runs=1;;5;
+Map[
+ListPlot[#[[1,runs,2,1;;3000]],
+Joined->True,
+PlotRange->Automatic,
+Frame->True,
+FrameLabel->{"timesteps","number of triangles"},
+ImageSize->300
+]&,gdata]