怎么用MATLAB代码实现，取n=2,3,4,分别用复化梯形公式和复化Simpson公式计算积分

1. 复合梯形公式：

function I = composite_trapezoidal(f, a, b, n)
% f: 被积函数
% a, b: 积分区间 [a,b]
% n: 子区间数
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

2. 复合Simpson公式：

function I = composite_simpson(f, a, b, n)
% f: 被积函数
% a, b: 积分区间 [a,b]
% n: 子区间数（必须为偶数）
h = (b-a)/n; % 步长
x = linspace(a, b, n+1); % 子区间节点
y = f(x); % 子区间节点上的函数值
I = h/3 * (y(1) + 4*sum(y(2:2:end-1)) + 2*sum(y(3:2:end-2)) + y(end)); % Simpson公式
end

使用示例：

f = @(x) sin(x);
a = 0;
b = pi;
n = 2;
I1 = composite_trapezoidal(f, a, b, n);
I2 = composite_simpson(f, a, b, n);