Modified random distribution: Difference between revisions
m
→{{header|Julia}}
mNo edit summary |
|||
Line 29:
modifier(x) = (y = 2x - 1; y < 0 ? -y : y)
modrands(rands1, rands2) = [x for (i, x) in enumerate(rands1) if rands2[i] < modifier(x)]
histogram(modrands(rand(10000), rand(10000)), nbins = 15
</lang>{{out}}
<pre>
[0.0 , 0.
[0.05, 0.1 ) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 442
[0.1 , 0.15) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 389
[0.
[0.
[0.
[0.
[0.
[0.
[0.45, 0.5 ) ┤▇ 20
[0.5 , 0.55) ┤▇▇ 24
[0.55, 0.6 ) ┤▇▇▇▇▇ 66
[0.6 , 0.65) ┤▇▇▇▇▇▇▇▇▇ 119
[0.65, 0.7 ) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇ 174
[0.7 , 0.75) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 232
[0.75, 0.8 ) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 288
[0.8 , 0.85) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 319
[0.85, 0.9 ) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 374
[0.9 , 0.95) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 459
[0.95, 1.0 ) ┤▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 433
└ ┘
Frequency
</pre>
|