matlab复合梯形公式代码

以下是使用复合梯形公式计算定积分的 MATLAB 代码：

function I = comp_trapezoidal(f, a, b, n)
% 使用复合梯形公式计算定积分
% f: 被积函数
% a, b: 积分区间
% n: 区间数
% I: 积分值

h = (b - a) / n; % 子区间宽度
x = linspace(a, b, n+1); % 子区间节点
y = f(x); % 子区间端点函数值
I = h/2 * (y(1) + 2*sum(y(2:end-1)) + y(end)); % 计算积分值
end

其中，`f` 是被积函数，`a` 和 `b` 是积分区间的端点，`n` 是子区间数。函数首先计算子区间宽度 `h` 和节点 `x`，然后计算每个子区间的端点函数值 `y`，最后使用复合梯形公式计算积分值 `I`。