方法一：
pi=0
N=100
for i in range(N):
pi+=1/pow(16,i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6))
print(pi)
方法二：
import random
N=1000*1000
k=0
for i in range(N):
x,y=(),()
dist=pow(x**2+y**2,0.5)
if dist<=1.0:
k+=1
pi=4*(k/N)
print(pi)