RosettaCodeData/Task/Monte-Carlo-methods/Julia/monte-carlo-methods.julia
2020-02-17 23:21:07 -08:00

11 lines
260 B
Text

using Printf
function monteπ(n)
s = count(rand() ^ 2 + rand() ^ 2 < 1 for _ in 1:n)
return 4s / n
end
for n in 10 .^ (3:8)
p = monteπ(n)
println("$(lpad(n, 9)): π ≈ $(lpad(p, 10)), pct.err = ", @sprintf("%2.5f%%", abs(p - π) / π))
end