from numpy import random
import math

iterations = 200
samples = 1000

def M(x,y):
    z = x+y*1j
    k = 0
    while abs(z)<2 and k<iterations:
      z = z*z+ x+y*1j
      k += 1
    return math.floor(k/iterations)

sum = 0
n = 0
while n < samples:
   n += 1
   sum += M(-2 + 3*random.rand(), -1.5 +3*random.rand())

final = 9*sum/samples
print(final);