import numpy as np
def pi():
n=100
t=0
for i in range(1,n+1):
x=np.random.random()
y=np.random.random()
if x**2+y**2<=1:
t=t+1
return 4*t/n
res = []
for _ in range(100):
res.append(pi())
print(f'平均:{np.mean(res)}')
print(f'標準偏差:{np.std(res)}')
aW1wb3J0IG51bXB5IGFzIG5wCgoKZGVmIHBpKCk6CiAgICBuPTEwMAogICAgdD0wCiAgICBmb3IgaSBpbiByYW5nZSgxLG4rMSk6CiAgICAgICAgeD1ucC5yYW5kb20ucmFuZG9tKCkKICAgICAgICB5PW5wLnJhbmRvbS5yYW5kb20oKQogICAgICAgIGlmIHgqKjIreSoqMjw9MTogCiAgICAgICAgICAgIHQ9dCsxCiAgICByZXR1cm4gNCp0L24KCgpyZXMgPSBbXQpmb3IgXyBpbiByYW5nZSgxMDApOgogICAgcmVzLmFwcGVuZChwaSgpKQoKCnByaW50KGYn5bmz5Z2H77yae25wLm1lYW4ocmVzKX0nKQpwcmludChmJ+aomea6luWBj+W3ru+8mntucC5zdGQocmVzKX0nKQo=