sites=Range[-100,100];spread[p_]:=Sqrt[p.sites^2-(p.sites)^2];evolution=quantumWalk[100];ListLinePlot[{Table[{t,spread[evolution〚t+1〛]},{t,0,100}],Table[{t,Sqrt[t]},{t,0,100}]},PlotLegends{"quantum (~ t)","classical (~ Sqrt[t])"},AxesLabel{"steps","std. dev."}]