matlab辛普森法求积分

使用 MATLAB 中的 Simpson 函数来计算函数的定积分。

假设要计算函数 f 在区间 [a, b] 上的定积分，可以使用以下代码：

a = 0; % 区间左端点
b = pi/2; % 区间右端点
n = 100; % 将区间分成 n 个子区间
h = (b-a)/n; % 子区间宽度

% 定义函数 f，例如 f = @(x) sin(x)
% 计算 Simpson 积分
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));

其中，h 是每个子区间的宽度，x 是子区间的节点，y 是函数在节点处的取值，I 是 Simpson 积分的结果。

需要注意的是，Simpson 积分要求区间分成偶数个子区间，因此在计算时需要将 n 转化为偶数。