复合梯形公式代码matlab

以下是复合梯形公式的 MATLAB 代码示例：

function [I] = composite_trapezoidal(f, a, b, n)
% f: 被积函数
% a，b：积分区间
% n：区间数

h = (b - a) / n; % 子区间长度
x = a:h:b; % 子区间节点
y = f(x); % 子区间的函数值
y(1) = f(a); % 强制更新左端点的函数值
y(n+1) = f(b); % 强制更新右端点的函数值
I = h * sum(y) - h / 2 * (y(1) + y(n+1)); % 计算积分值
end

该函数接受四个输入参数，分别是被积函数 f、积分区间的左右端点 a 和 b，以及区间数 n。函数返回复合梯形公式的积分值。