要求：

（1）使用蒙特卡羅方法

（2）打印進度條

（3）計算用時

蒙特卡羅方法（向一個單位圓及其單位正方形拋點,取其四分之一來算，結果再乘上4即可）

---

話不多說，直接上代碼

import time import random dian=pow(10,6) #拋一百萬個點,別拋太多，會用時過多，當拋到一定數量的點后，其精確值都差不多了 youxiaodian=0 #落在圓內的點 start=time.perf_counter() #計時 for i in range(1,dian+1):a='*'*int((i/dian)*100) #a,b,c都是打印進度條用b='.'*int(((dian-i)/dian)*100)c=int((i/dian)*100)print("\r當前進度：{:^3.0f}% [{}->{}]".format(c, a, b),end="") #打印進度條x,y=random.random(),random.random() #x，y取0-1的隨機數（因為是單位圓）if pow((x*x+y*y),1/2)<1: #r=1,點(x,y)到圓心的距離=x平方+y平方開根號youxiaodian+=1 #小于半徑的話就表示落在1/4圓內，園內點+1 pi=4*(youxiaodian/dian) #剛剛計算的只是1/4圓周率，所以×4 print("\n") print("圓周率pi接近于{}, 計算用時{}秒".format(pi,time.perf_counter()-start).center(100,"="))

---

運行結果：

當前進度：100% [****************************************************************************************************->]==================================圓周率pi接近于3.13702, 計算用時11.517367秒==================================

注：每次運行的pi值并不是一樣的，因為引用了random函數，結果只能接近于真實的pi值，要想運行的結果每次都一樣，可以采用random.seed（n）獲取隨機種子。

總結

以上是生活随笔為你收集整理的圆周率计算——蒙特卡罗方法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。