python 蒙特卡洛算法求圆周率

import random
k=0
sum=0
n=int(input('请输入要做的实验次数：'))
for i in range(n):
    x,y=random.uniform(-1,1),random.uniform(-1,1)
    if x**2+y**2<=1:
        k+=1
        sum+=1
    else:
        sum+=1
pi=(k/sum)*4
print('经过%d次实验之后所得圆周率的值为%f'%(n,pi))