Revision as of 10:06, 6 March 2010 (view source) Christian (talk \| contribs) (PureBasic) ← Older edit		Revision as of 14:49, 1 April 2010 (view source) 79.37.200.72 (talk) (→‎{{header\|D}}) Newer edit →
Line 243: =={{header\|D}}== D V.2. <lang d>▼ <lang d>import std.stdio:, ~~writefln~~std.random, std.math; ~~import std.random: rand;~~ double pi(int nthrows) { int inside; ~~for~~foreach (~~int~~ i; i0 <.. nthrows~~; i++~~) { ~~double~~if r1(uniform(0,1.0) =^^ ~~rand~~2 + uniform(0,1.0) /^^ ~~cast(double~~2 <= 1)~~uint.max;~~ ~~double r2 = rand() / cast(double)uint.max;~~ ~~if (r1r1 + r2r2 <= 1.0)~~ inside++; } return 4.0 * inside / nthrows; } void main() { foreach (p; 1 .. 9) ~~for (int n = 10_000; n <= 100_000_000; n *= 10)~~ writefln("%9d11s: %07f", n10 ^^ p, pi(n10 ^^ p)); ▲}</lang d> } ~~</lang>~~ Sample output: <pre> ~~100000~~ 10: 32.~~137760~~400000 ~~1000000~~ 100: 3.~~143508~~200000 ~~10000000~~ 1000: 3.~~140836~~136000 ~~100000000~~ 10000: 3.~~141666~~170000 100000: 3.141680 ~~</pre>~~ 1000000: 3.141884 ~~For much faster results you can use the pseudorandom generator of Tango, or another one like R250-521.~~ 10000000: 3.141182 100000000: 3.141838</pre> =={{header\|E}}==

Monte Carlo methods (view source)