fork download
  1. import numpy as np
  2.  
  3.  
  4. def pi():
  5. n=100
  6. t=0
  7. for i in range(1,n+1):
  8. x=np.random.random()
  9. y=np.random.random()
  10. if x**2+y**2<=1:
  11. t=t+1
  12. return 4*t/n
  13.  
  14.  
  15. res = []
  16. for _ in range(100):
  17. res.append(pi())
  18.  
  19.  
  20. print(f'平均:{np.mean(res)}')
  21. print(f'標準偏差:{np.std(res)}')
  22.  
Success #stdin #stdout 0.18s 28816KB
stdin
Standard input is empty
stdout
平均:3.1463999999999994
標準偏差:0.16230539116123038