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