MC算法的MATLAB程序

以下是一个简单的MC算法的MATLAB程序，用于求解定积分：

function [I, E] = monte_carlo(f, a, b, N)
% f: 被积函数
% a, b: 积分区间
% N: 抽样点数

% 计算积分区间的长度
L = b - a;

% 生成均匀分布的随机数
x = a + L * rand(1, N);

% 计算被积函数值
fx = feval(f, x);

% 计算积分值
I = L * sum(fx) / N;

% 计算误差
E = L * std(fx) / sqrt(N);

使用方法：

1. 定义被积函数`f`，例如`f = @(x)exp(-x.^2)`表示积分区间内的高斯函数。
2. 定义积分区间`a`和`b`，例如`a = 0`，`b = 1`表示对`[0,1]`区间内的函数进行积分。
3. 定义抽样点数`N`，例如`N = 1000`表示使用1000个随机点进行积分。
4. 调用`monte_carlo`函数，例如`[I, E] = monte_carlo(f, a, b, N)`表示计算积分值和误差。