Monte Carlo Method and pi

[Karl Pech]

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

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!

Regards
Karl