复合辛普森求积matlab

以下是复合辛普森求积的Matlab实现代码：

function I = simpson(f, a, b, n)
% f: 被积函数
% a: 积分下限
% b: 积分上限
% n: 区间数，必须为偶数
h = (b - a) / n;
x = a:h:b;
y = f(x);
I = h / 3 * (y(1) + 4 * sum(y(2:2:end-1)) + 2 * sum(y(3:2:end-2)) + y(end));
end

其中，`f`为被积函数，`a`和`b`为积分下限和上限，`n`为区间数，必须为偶数。函数返回积分结果`I`。

使用方法如下：

1.定义被积函数，例如：

f = @(x) exp(-x.^2);

2.调用`simpson`函数进行积分计算，例如：

I = simpson(f, 0, 1, 100);

这将计算函数`f`在区间`[0,1]`上的积分，使用100个区间进行复合辛普森求积。