Monte Carlo Method and pi

Karl Pech KarlPech at
Thu Jul 8 23:30:01 CEST 2004


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

print 4*sy/n                             # compute pi

Unfortunately, even for n = 2000000 the result is ~2,13... .
It converges very slow to pi and I don't know why.

Please help me!


More information about the Python-list mailing list