diff --git a/triangle_analysis.m b/triangle_analysis.m
index 5cd064dd2563dbfddce90d6a3d9d2d2b5c66bdf6..6fc0b5eadaadb18af35e1603f44ed44e66ac3f85 100644
--- a/triangle_analysis.m
+++ b/triangle_analysis.m
@@ -4,6 +4,7 @@ Quit[]
 
 
 Needs["ErrorBarPlots`"]
+Needs["MaTeX`"]
 
 
 (* ::Section:: *)
@@ -179,8 +180,9 @@ plot4=ListPlot[coarseData[[16;;20]],Joined->True,PlotRange->{0*minCount,maxCount
 
 
 (* For export *)
-numPlots=20;
-selectedData=gdata[[2,1]][[-numPlots;;-1]];
+numPlots=25;
+selectedPlots={6,7,8,11,12,13,16,17,18,21,22,23};
+selectedData=gdata[[2,1]][[selectedPlots]];
 measureSkip=1;
 minCount=Min[Map[Min[#[[2]]]&,selectedData]];
 maxCount=Max[Map[Max[#[[2]]]&,selectedData]];
@@ -189,7 +191,10 @@ maxTime=Max[Map[Length[#[[2]]]&,selectedData]];
 skipPts=Max[1,Round[maxTime/5000]]; (* Plotting every point is slow. Plot only once per `skipPts` timesteps *)
 coarseData=Map[#[[2,1;;maxTime;;skipPts]]&,selectedData];
 labels=Map["{n,tau} = "<>ToString[#[[1]]]&,selectedData];
-plotTimeEvol=ListPlot[coarseData,Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime},Frame->True,FrameLabel->{"timesteps","number of triangles"},ImageSize->300]
+plotTimeEvol=ListPlot[coarseData,Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime},
+Frame->True,FrameLabel->{MaTeX["\\text{timesteps}"],MaTeX["\\text{number of triangles}"]},
+PlotLabel->MaTeX["n=1000,\\; \\tau = 2.2"],
+ImageSize->250]
 (* Map[ListPlot[#,Joined->True,PlotRange\[Rule]{minCount,maxCount},DataRange\[Rule]{0,maxTime}]&,coarseData] *)