Anonymous user
Monte Carlo methods: Difference between revisions
added Ursala
(added Ursala) |
|||
Line 767:
1e6 => 3.141924
1e8 => 3.14167724</pre>
=={{header|Ursala}}==
<lang Ursala>#import std
#import flo
mcp "n" = times/4. div\float"n" (rep"n" (fleq/.5+ sqrt+ plus+ ~~ sqr+ minus/.5+ rand)?/~& plus/1.) 0.</lang>
test program:
<lang Ursala>#cast %eL
pis = mcp* <10,100,1000,10000,100000,1000000></lang>
output:
<pre><
2.800000e+00,
3.600000e+00,
3.164000e+00,
3.118800e+00,
3.144480e+00,
3.141668e+00></pre>
| |||