蒙特卡洛法计算π

# 设一个面积为4的正方形，和一个直径为1的圆形，其中心位于xy坐标轴原点
# 圆形面积/正方形面积=π/4
# 因此π=圆形面积/正方形面积*4
# 而圆形面积/正方形面积可以通过投点概率表示

import random
import math

m =0
	# 用于计数：在圆内的点数
i =0
	# 用于计数：在正方形内的点数
while i < 1000000:
		#投点1000000次
    a =random.uniform (-1,1)
    b =random.uniform (-1,1)
    	#投点范围在xy坐标轴上的限定
    c = math.sqrt(a*a+b*b)
    	#计算是否在圆内
    if c < 1:
        m =m+1
        	#计算在圆内的点数
    i =i +1 
    		#计算投点次数
print (m/i*4)
		#计算圆的面积