Monte Carlo methods: Difference between revisions

</pre>

<lang lisp>(def pi-seq

(pmap first

(iterate

(fn [[pi within total]]

[ (double (* (/ within total) 4))

(if (<= (+ (Math/pow (rand) 2) (Math/pow (rand) 2)) 1) (inc within) within)

(inc total) ] )

[0.0 0N 1N])))

(nth pi-seq 1000000)

</lang>

{{out}}

<pre>

3.140492