Percolation/Mean run density: Difference between revisions
Rename Perl 6 -> Raku, alphabetize, minor clean-up
Thundergnat (talk | contribs) (Rename Perl 6 -> Raku, alphabetize, minor clean-up) |
|||
Line 95:
0.9 100000 0.0900 0.0900 +0.0000 (+0.03%)
</pre>
=={{header|C++}}==
<lang cpp>#include <algorithm>
Line 956 ⟶ 957:
R(n, p)= 0.096200
R(n, p)= 0.091730</pre>
=={{header|Perl 6}}==▼
<lang perl6>sub R($n, $p) { [+] ((rand < $p) xx $n).squish }▼
say 't= ', constant t = 100;▼
for .1, .3 ... .9 -> $p {▼
say "p= $p, K(p)= {$p*(1-$p)}";▼
for 10, 100, 1000 -> $n {▼
printf " R(%6d, p)= %f\n", $n, t R/ [+] R($n, $p)/$n xx t▼
}▼
}</lang>▼
{{out}}▼
<pre>t= 100▼
p= 0.1, K(p)= 0.09▼
R( 10, p)= 0.088000▼
R( 100, p)= 0.085600▼
R( 1000, p)= 0.089150▼
p= 0.3, K(p)= 0.21▼
R( 10, p)= 0.211000▼
R( 100, p)= 0.214600▼
R( 1000, p)= 0.211160▼
p= 0.5, K(p)= 0.25▼
R( 10, p)= 0.279000▼
R( 100, p)= 0.249200▼
R( 1000, p)= 0.250870▼
p= 0.7, K(p)= 0.21▼
R( 10, p)= 0.258000▼
R( 100, p)= 0.215400▼
R( 1000, p)= 0.209560▼
p= 0.9, K(p)= 0.09▼
R( 10, p)= 0.181000▼
R( 100, p)= 0.094500▼
R( 1000, p)= 0.091330</pre>▼
=={{header|Phix}}==
Line 1,053 ⟶ 1,020:
0.9 10000 0.0901 0.0900 +0.0001 (+0.11%)
0.9 100000 0.0900 0.0900 +0.0000 (+0.03%)
</pre>
=={{header|Python}}==
Line 1,101 ⟶ 1,068:
t=500 p=0.90 n= 4096 p(1-p)=0.090 sim=0.090 delta=0.2%
t=500 p=0.90 n=16384 p(1-p)=0.090 sim=0.090 delta=0.0%</pre>
=={{header|Racket}}==
Line 1,237 ⟶ 1,203:
10000 0.0906 0.000568
100000 0.0908 0.000775</pre>
(formerly Perl 6)
▲<lang perl6>sub R($n, $p) { [+] ((rand < $p) xx $n).squish }
▲say 't= ', constant t = 100;
▲for .1, .3 ... .9 -> $p {
▲ say "p= $p, K(p)= {$p*(1-$p)}";
▲ for 10, 100, 1000 -> $n {
▲ printf " R(%6d, p)= %f\n", $n, t R/ [+] R($n, $p)/$n xx t
▲ }
▲}</lang>
▲{{out}}
▲<pre>t= 100
▲p= 0.1, K(p)= 0.09
▲ R( 10, p)= 0.088000
▲ R( 100, p)= 0.085600
▲ R( 1000, p)= 0.089150
▲p= 0.3, K(p)= 0.21
▲ R( 10, p)= 0.211000
▲ R( 100, p)= 0.214600
▲ R( 1000, p)= 0.211160
▲p= 0.5, K(p)= 0.25
▲ R( 10, p)= 0.279000
▲ R( 100, p)= 0.249200
▲ R( 1000, p)= 0.250870
▲p= 0.7, K(p)= 0.21
▲ R( 10, p)= 0.258000
▲ R( 100, p)= 0.215400
▲ R( 1000, p)= 0.209560
▲p= 0.9, K(p)= 0.09
▲ R( 10, p)= 0.181000
▲ R( 100, p)= 0.094500
▲ R( 1000, p)= 0.091330</pre>
=={{header|Sidef}}==
| |||