Monte Carlo Method and pi

Nick Smallbone nick at nick8325.freeserve.co.uk
Thu Jul 8 23:53:12 CEST 2004


"Karl Pech" <KarlPech at users.sf.net> wrote in message
news:ccke8p$c0h$06$1 at news.t-online.com...
> Hi,
>
>
> I got a task there I have to compute pi using the
> Method above.
>
> So I wrote the following program:
> ---
> import random
> import math
>
> n = long(raw_input("Please enter the number of iterations: "))
> sy = 0                                   # sum of the function-values
>
> for i in range(0, n):
>   x = random.random()                    # coordinates
>   sy += math.sqrt(1-math.sqrt(x))        # calculate y and add to sy

This line is wrong. I think it should be sy += math.sqrt(1 - x ** 2). Doing
that gives me ~3.1 for 1000 iterations.

Nick





More information about the Python-list mailing list