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)) }