Anonymous user
Percolation/Mean run density: Difference between revisions
no edit summary
(Go solution) |
No edit summary |
||
Line 545:
</pre>
=={{header|Mathematica}}==
<lang Mathematica>meanRunDensity[p_, len_, trials_] :=
Mean[Length[Cases[Split@#, {1, ___}]] & /@
Unitize[Chop[RandomReal[1, {trials, len}], 1 - p]]]/len
Column@Table[
Grid[Join[{{p, n, K, diff}},
Table[{q, n, x = meanRunDensity[q, n, 100] // N,
q (1 - q) - x}, {n, {100, 1000, 10000, 100000}}], {}],
Alignment -> Left], {q, {.1, .3, .5, .7, .9}}]</lang>
{{out}}
<pre>
p n K diff
0.1 100 0.0905 -0.0005
0.1 1000 0.0900 -0.00001
0.1 10000 0.0902 -0.00015
0.1 100000 0.0901 -0.0001265
p n K diff
0.3 100 0.2088 0.0012
0.3 1000 0.2101 -0.00011
0.3 10000 0.2099 0.000049
0.3 100000 0.2100 -0.0000352
p n K diff
0.5 100 0.2533 -0.0033
0.5 1000 0.2515 -0.00146
0.5 10000 0.2501 -0.000131
0.5 100000 0.2500 -0.0000425
p n K diff
0.7 100 0.2172 -0.0072
0.7 1000 0.2106 -0.0006
0.7 10000 0.2098 0.000194
0.7 100000 0.2102 -0.0002176
p n K diff
0.9 100 0.0924 -0.0024
0.9 1000 0.0895 0.00049
0.9 10000 0.0899 0.00013
0.9 100000 0.0900 -0.0000144
</pre>
=={{header|Perl}}==
{{trans|Perl 6}}
| |||