SetDirectory["/Users/takasu/home/Granada/数値計算一般/乱数発生/Random_Mersenne-Twister/build/Development/"]

results = ReadList["data-binomial",Real];
Length[results]

Take[results,10]

gSimHisto=Histogram[results,10,PlotRange->{{0,10},All}]

freq=Reverse[Sort[Tally[results]],2]

(* P = 1/2 の二項分布の定義*)
prob[i_, n_]:=n!/i!/(n-i)!/2^n

prob[2,10]

numHead = Table[i,{i,0,10}]

probSeq = Map[prob[#,10]&,numHead]

points = Transpose[{numHead+0.5,1000probSeq}]

gSol=ListPlot[points,Joined->True]

Show[gSimHisto,gSol,PlotRange->All]