蒙特卡罗的方法估计自然常数 matlab代码

下面是一个使用蒙特卡罗方法估计自然常数e的Matlab代码示例：

N = 1000000; % 模拟次数
count = 0; % 计数器
for i = 1:N
    x = rand; % 生成0-1之间的随机数
    if x <= exp(1) % 如果在曲线下面就计数器加1
        count = count + 1;
    end
end
e_est = count/N*(1/mean(diff(linspace(0,1)))) % 输出估计值

这段代码中，我们生成了N个随机点，然后检查每个点是否在曲线y=e的下面，如果是，计数器就加1。最后，我们将计数器除以总模拟次数，再乘以区间长度，得到e的估计值。