蒙特卡洛（Monte Carlo）方法求面積

如圖，刷微博時，看到一個問題，第一個想到的就是用蒙特卡洛方法求解，當時正在練python，於是嘗試用python編寫程序。

 

import random
# 先求s1
k=0
n=100000000
for i in range(n):
    x=random.uniform(0,10)
    y=random.uniform(0,10)
    if ((x-5)**2+(y-5)**2>25) and (y<-2*x+20):
        k=k+1
    else: k=k

s1=(k/n)*100

#求s2
import math
s2=(10*10-math.pi*25)/4

s=s1+s2
print("k=",k)
print("s2=",s2)
print("s1=",s1)
print("s=",s)