蒙特卡洛（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)