RosettaCodeData/Task/Monte-Carlo-methods/Sidef/monte-carlo-methods.sidef
2017-09-25 22:28:19 +02:00

9 lines
214 B
Text

func monteCarloPi(nthrows) {
4 * (^nthrows -> count_by {
hypot(1.rand(2) - 1, 1.rand(2) - 1) < 1
}) / nthrows
}
for n in [1e2, 1e3, 1e4, 1e5, 1e6] {
printf("%9d: %07f\n", n, monteCarloPi(n))
}