Python-蒙特卡罗应用

第1关：蒙特卡洛方法求定积分一

任务描述

本关任务：用蒙特卡罗方法求函数f(x)=（x/25+1/5），在区间[a,b]中定积分。假设a=0,b=1，即求

 

要求将函数f(x)定义为匿名函数，求出的积分保留5位小数，a,b的值可以自定义。

import numpy as np
#将积分函数f定义成匿名函数
# # # # # # # # # # #begin # # # # # # # # # # # #
f = lambda x: x / 25 + 1 / 5

# # # # # # # # # # #end # # # # # # # # # # # # # #
#在一行输入积分区间【 a, b】 和实验点数n， 用空格分隔
# # # # # # # # # # # # # # # # # #begin # # # # # # # # # # # # # # #
m = input()
m = list(map(int, m.split()))

a = int(m[0])

b = int(m[1])
n = int(m[2])
# # # # # # # # # # # # # # # # # # #end # # # # # # # # # # # # # # # #
#矩形区域为：[a, b, fmin, fmax], 矩形面积 = (b - a) * (fmax - fim)
fmin = 0
fmax = f(b)
np.random.seed(0)
#利用均匀分布产生n个在【