蒙特卡洛模拟MATLAB代码

以下是一个简单的蒙特卡洛模拟的MATLAB代码示例：

% 定义模拟次数
n = 10000;

% 定义问题的参数
a = 2;
b = 3;

% 定义函数
f = @(x) x.^2 + a*x + b;

% 定义积分区间
xmin = -5;
xmax = 5;

% 计算积分
sum = 0;
for i = 1:n
    % 随机生成x
    x = xmin + (xmax - xmin) * rand();
    
    % 计算函数值
    y = f(x);
    
    % 累加
    sum = sum + y;
end

% 计算积分值
integral = (xmax - xmin) * sum / n;
disp(integral);

这个示例计算了函数 $f(x) = x^2 + ax + b$ 在区间 $[-5, 5]$ 上的积分，其中 $a=2$，$b=3$。它使用了 $n=10000$ 次的随机模拟来估算积分值。