计算一维积分的蒙特卡罗算法

用蒙特卡罗法计算定积分：

我们可以用c语言来实现，这里我们主要介绍如何用matlab来实现算法：

具体实现代码如下：

function s=monte_carlo(a,b,n)
t=rand(1,n);
x=a+(b-a)*t;
s=sum(monte_carlo_f(x));
s=s*(b-a)/n;
end

function f=monte_carlo_f(x)
f=x.^2+sin(x);
end

分别在matlab中运行如下两个代码，则在命令行窗口中会出现如下图所示：这时会出现输入参数不足的问题，我们只需要在命令行窗口中输入变量并赋值便可得到结果：

好了，今天的分享就到这里了，希望对大家有帮助。