蒙特卡洛法求pi

swimminginsky · 发表于 2020-3-10 22:09:54

马上注册，结交更多好友，享用更多功能^_^

您需要登录才可以下载或查看，没有账号？立即注册

x

import math
import random
import time

countIn = 0
count = 0
i = 0
rounds = 20000
startTime = time.perf_counter()
while i <= rounds:
x = random.randint(0, 10000)
#print(x)
y = random.randint(0, 10000)
#print(y)
count += 1
if math.sqrt(x * x + y * y) <= 10000:
      countIn += 1
      #print('pi  = ', float(countIn/count) * 4)
endTime = time.perf_counter()
i += 1

print('pi  = ', float(countIn/count) * 4, '  the escaped time is %f seconds'%(endTime - startTime))

请大家多多指教

wp231957 · 发表于 2020-3-10 22:15:31

要是没看到最后一行字，我还以为是秀代码的

账号		自动登录	找回密码
密码			立即注册

[学习笔记] 蒙特卡洛法求pi

马上注册，结交更多好友，享用更多功能^_^