Unbias a random generator: Difference between revisions
m
→{{header|Mathematica}}
(Added Fōrmulæ solution) |
|||
Line 1,309:
</pre>
=={{header|Mathematica}}/{{header|Wolfram Language}}==
<lang Mathematica>rand[bias_, n_] := 1 - Unitize@RandomInteger[bias - 1, n]
▲ DeleteCases[rand[bias, {n, 2}], {a_, a_}][[All, 1]]</lang>
▲<pre>count = 1000000;
TableForm[
Table[{n, Total[rand[n, count]]/count // N,
Total[#]/Length[#] &@unbiased[n, count] // N}, {n, 3, 6}],
TableHeadings -> {None, {n, "biased", "unbiased"}}]</lang>
{{out}}
<pre>n biased unbiased
3 0.33312 0.500074
4 0.24932 0.499883
5 0.1998 0.498421
6 0.16620 0.49805</pre>
=={{header|NetRexx}}==
|