rei.py

a=[0]
c=0
n=200000
import numpy
for i in range(2*n):
    b=numpy.random.rand(1,2)
    a.append(b)
for j in range(i,n):
    d=a[j]**2+a[2*j]**2
    if d <= 1:
        c+=1
pai=4*c/n
print(a[3])
print(a[3][0,1])