Advertisement
Guest User

Untitled

a guest
Sep 9th, 2014
245
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.30 KB | None | 0 0
  1. n = 60000; q = 0.495; NSim = 10^5; sigma = Sqrt[4 q (1 - q)];
  2. sims = Table[Total[RandomChoice[{q, 1 - q} -> {1, -1}, n]], {NSim}];
  3. distr = HistogramDistribution[sims];
  4. Show[
  5. DiscretePlot[PDF[distr, x], {x, -1200, 600}],
  6. Plot[PDF[NormalDistribution[-600, Sqrt[n]*sigma], x], {x, -1200,
  7. 600}]
  8. ]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement