用matlab编写复合梯形公式Tn=h/2*[f(a)+f(b)+2∑f(xk)]的程序

可以按照以下方式编写复合梯形公式Tn的Matlab程序：

function T = compositeTrapezoidal(f, a, b, n)
% f: 被积函数
% a, b: 积分区间
% n: 子区间个数

h = (b - a) / n; % 子区间宽度
x = linspace(a, b, n+1); % 子区间节点
y = f(x); % 子区间节点处的函数值
T = h/2 * (y(1) + 2*sum(y(2:end-1)) + y(end)); % 复合梯形公式
end

其中，`f`是被积函数，`a`和`b`是积分区间的端点，`n`是子区间的个数。函数中使用`linspace`函数生成了节点向量`x`，并且通过`f`函数计算出节点处的函数值`y`。最后，根据复合梯形公式计算出积分值`T`。