(* ::Package:: *)

(* ::Section:: *)

(*P(Z^n | start 011..11) on the chain {1,...,n}*)

rawdata=Import[NotebookDirectory[]<>"data/data_reach_chain_end.m"];

data=GatherBy[rawdata,{#[[2]]}&];

(* data[[n index, prob index, p / n / prob / runs within timelimit / total runs]] *)

nlabels=Map["n = "<>ToString[#[[1,2]]]&,data];

plotZoomed=ListPlot[data[[All,All,{1,3}]],

Joined->True,

PlotMarkers->Automatic,

PlotRange->{{0.6,0.7},{0,0.5}},

Frame->True,

ImageSize->250,

Epilog->{Dashed,Line[{{1-1/E,0},{1-1/E,1}}]}

];

plotBoth=ListPlot[data[[All,All,{1,3}]],

Joined->True,

FrameLabel->{"p (probability of BAD)","P(\!\(\*SuperscriptBox[\(Z\), \(n\)]\) | start 011..11)"},

PlotMarkers->Automatic,

PlotLegends->nlabels,

PlotRange->{{0,1},{0,1}},

Frame->True,

ImageSize->500,

Epilog->Inset[plotZoomed,Scaled[{0.32,0.6}]]]