蒙特卡羅方法python簡單實現

蒙特卡羅方法是一種計算方法。原理是通過大量隨機樣本，去了解一個系統，進而得到所要計算的值。

1、估算pi值

import random
n = 10000
k = 0
for i in range(n):
    x = random.uniform(-1,1)
    y = random.uniform(-1,1)
    if x**2 + y**2 <1:
        k +=1
print(4*float(k)/float(n))

2、計算y=x^2在[0,1]的積分

import random
n = 1000
k = 0
for i in range(n):
    x = random.uniform(0,1)
    y = random.uniform(0,1)
    if x**2 > y:
        k+=1
print(float(k)/float(n))