diff --git a/triangle_exponent_plots.m b/triangle_exponent_plots.m
index 6c3db7cc1c48c7859cf557ae0ec1e19386e076fc..33c4d4f69d668b1422514564df1dbd7b89a92ceb 100644
--- a/triangle_exponent_plots.m
+++ b/triangle_exponent_plots.m
@@ -20,6 +20,7 @@ Needs["ErrorBarPlots`"]
 (* graphdata_exponent_highN.m *)
 (* graphdata_properties2.m *)
 (* graphdata_canonical_properties.m *)
+(* graphdata_exponent_accurate1.m *)
 gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent_hightau.m"];
 gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *)
 averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
@@ -69,6 +70,7 @@ mediansFitsExtra=Map[LinearModelFit[#,logn,logn]&,mediansLoglogdata];
 
 canonicalloglogdata=Log[canonicalDatapoints[[All,nRange]]];
 canonicalFitsExtra=Map[LinearModelFit[#,logn,logn]&,canonicalloglogdata];
+canonicalNLMs=Map[NonlinearModelFit[#,factor*n^exponent,{factor,exponent},n]&,canonicalDatapoints];
 
 
 averagesFitsExtra[[1]]["ParameterConfidenceIntervalTable"]
@@ -94,21 +96,28 @@ mediansExponents=Transpose[{tauValues,mediansFits[[All,2,1]]}];
 Show[ListPlot[{averagesExponents,mediansExponents},Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","exponent"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]]
 
 
+canonicalNLMs[[1]]["ParameterConfidenceIntervals"]
+
+
 (* ::Subsection:: *)
 (*T(\[Tau]) including error bars*)
 
 
 (* For visual, shift the tau values slightly left or right to distinguish the two datasets *)
 tauValues=averagesGrouped[[All,1,1,1,2]];
+canonicalTauValues=averagesGrouped2[[All,1,1,1,2]];
 averagesExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
 Transpose[{tauValues-0.003,averagesFitsExtra}]];
 mediansExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
 Transpose[{tauValues+0.003,mediansFitsExtra}]];
 canonicalExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
-Transpose[{tauValues+0.000,canonicalFitsExtra}]];
+Transpose[{canonicalTauValues+0.000,canonicalFitsExtra}]];
+canonicalExponentsErrorBars2=Map[{{#[[1]],exponent/.#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-exponent/.#[[2]]["BestFitParameters"]]}&,
+Transpose[{canonicalTauValues+0.000,canonicalNLMs}]];
 plot2=Show[
 ErrorListPlot[{averagesExponentsErrorBars,mediansExponentsErrorBars,canonicalExponentsErrorBars},
 Joined->True,PlotMarkers->Automatic,
+PlotLegends->Placed[{"average","median","canonical"},{Left,Bottom}],
 Frame->True,FrameLabel->{"tau","triangle exponent"},
 PlotRange->{{2,3},{0,1.6}},
 ImageSize->300],