使用蒙托卡罗方法求解图形面积

求解该图形面积主要关键点在于求解左下角的直线下面部分和圆外的小面积s,阴影部分可以看做一个10*10的正方形和圆围成的圆外面积减去开始提到的小面积s,求解代码如下:

#使用蒙托卡罗方法求解图形面积
import math
import random
import time

def TEST_DO_NOT_CHANGE():
    result = -1
    ########## 开始作答
    s = 10000000
    fit = 0
    git = 0
    #start = time.perf_counter()
    for i in range(s):
        x = random.uniform(0,10)
        y = random.uniform(0,10)
        z = math.sqrt((x-5)**2+(y-5)**2)
        if z > 5:
            fit = fit + 1
        if z >5 and y < 0.5 *x:
            git = git +1
    git = git - (fit/4)
    #end = time.perf_counter()
    #print("样本点有",s)
    #print("持续时间",end-start)
    result = ((fit-git) /s)*100     
        

    ########## 结束作答  
    return result        
print (TEST_DO_NOT_CHANGE())

 

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值