Monte Carlo methods: Difference between revisions
Content added Content deleted
No edit summary |
m (correct percent versus ratio) |
||
Line 1,865: | Line 1,865: | ||
for n in 10 .^ (3:8) |
for n in 10 .^ (3:8) |
||
p = monteπ(n) |
p = monteπ(n) |
||
println("$(lpad(n, 9)): π ≈ $(lpad(p, 10)), pct.err = ", @sprintf("%2.5f%%", abs(p - π) / π)) |
println("$(lpad(n, 9)): π ≈ $(lpad(p, 10)), pct.err = ", @sprintf("%2.5f%%", 100 * abs(p - π) / π)) |
||
end</syntaxhighlight> |
end</syntaxhighlight> |
||
{{out}} |
{{out}} |
||
<pre> 1000: π ≈ 3.224, pct.err = 0.02623% |
<pre> 1000: π ≈ 3.224, pct.err = 0.02623% |
||
10000: π ≈ 3.1336, pct.err = 0. |
10000: π ≈ 3.1336, pct.err = 0.254% |
||
100000: π ≈ 3.13468, pct.err = 0. |
100000: π ≈ 3.13468, pct.err = 0.220% |
||
1000000: π ≈ 3.14156, pct.err = 0. |
1000000: π ≈ 3.14156, pct.err = 0.001% |
||
10000000: π ≈ 3.1412348, pct.err = 0. |
10000000: π ≈ 3.1412348, pct.err = 0.011% |
||
100000000: π ≈ 3.14123216, pct.err = 0. |
100000000: π ≈ 3.14123216, pct.err = 0.011%</pre> |
||
=={{header|K}}== |
=={{header|K}}== |