使用蒙特卡罗技术进行数值积分的简单工具。

mcint的Python项目详细描述


一个可以用来计算多维积分的蒙特卡罗积分器。结果是数值逼近,即BR/>依赖于使用随机数产生。BR/>实例1<BR/>(X** 2)
BR/> DEF样本():描述Monte-Carlo样本是如何取而代之的:
产量。随机():BR/> BR/>结果,误差=McCIT.积分(整合式,SAMPPLER(),测度=1,n=100)
< BR/>打印“x** 2在0和1之间的积分近似为”,result

integrate()函数的第二个参数应该是iterable
表达式,在本例中,它是一个生成器。我们可以使用以下方法消除这个采样器:


result,error=mcint.integrate(integrand,iter(random.random,-1),measure=1.0,n=100)


这将从random.random()函数创建一个iterable对象,该函数将连续调用random.random(),直到返回-1(它永远不会返回介于0.0和1.0之间的值)。


dy`::

import random
import math


def integrand(x):
返回(x[0]**2+x[1]**2)


def sampler():
为真时:
y=random.random()
x=random.random()
如果x**2+y**2<;=1:
屈服(x,y)

result,error=mcint.integrate(integrand,sampler(),measure=math.pi/4)

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java ParsePushReceiver参数   java如何从设备读取完整数据?   java将java_设置为home,但忽略错误   java如何从歌曲中对专辑进行排序?   java libnaude+windows 10 x64+Eclipse   java如何将maven目标更改为如图所示的想法中的quickicon?   java swing布局中心面板,可滚动显示多个窗口窗格   使用MOSQUITO代理的mqtt中的java SSL   java如何通过属性值获取XML字符串   java在服务器每次启动时停止GWT编译   java如何让javac搜索类路径的子目录?   可以比较java中的两个不同类吗?   JAVAAndroid活动内部类中的lang.NoClassDefFoundError   java HttpServletRequest获取请求头参数块   C++socket与java客户端的连接   java如何在Apache commons http客户端上使用SSL客户端证书   使用预编译正则表达式模式提高java速度   JavaRhino将两个已编译脚本合并为一个脚本