python实验——第二次

使用蒙特卡洛方法计算圆周率近似值

蒙特·卡罗方法是一种通过概率统计来得到问题近似解的方法,在很多领域都有重要的应用,其中就包括圆周率近似值的计算问题。假设有一块边长为2的正方木板，上面画一个单位圆,然后随意往木板上掷飞镖,落点坐标必然在木板上(更多的时候是落在单位圆内）， 如果掷的次数足够多，那么落在单位圆内的次圆周率近似值的方法,如图3.1所示。编写程序,模拟蒙特·卡罗计算圆周率近似值的方法,输入掷飞镖次数,然后输出圆周率 。观察实验结果，理解实验结果随着模拟次数增多越来越接近圆周率的原因。

import random;
import time;

start=time.perf_counter()
n = 0;
m = 100000000;
i = 0;
while i <= m:
    x = random.uniform(-1.000001,1.000001);
    y = random.uniform(-1.000001,1.000001);
    if x*x+y*y<=1:
        n += 1;
    i += 1;

print(4*n/m);
print("程序运行时间为{}s".format(time.perf_counter()-start))