写一个关于复合梯形求积公式和复合辛普生求积公式的matlab代码

在MATLAB中，我们可以编写函数来计算复合梯形法和复合辛普森法则的积分。这里是一个简单的例子：

function result = composite_trapezoidal(f, a, b, n)
% 复合梯形法求积
h = (b - a) / n;
area_sum = f(a) + f(b);
for i = 1:n-1
    area_sum += 2 * f(a + i * h);
end
result = h * area_sum;

function result = composite_simpson(f, a, b, n)
% 复合辛普森法则求积
if rem(n, 2) == 0 % 如果n是偶数
    h = (b - a) / n;
    sub_areas = 4 * [f(a), cumsum(f(a + 2:2:(b - 2)))]';
    result = h / 3 * sum(sub_areas);
else % 如果n是奇数，需要处理左端点的一个额外项
    h = (b - a) / n;
    sub_areas = 4 * [f(a), cumsum(f(a + 1:2:(b - 1)))]';
    extra_area = f(a) + f(b);
    result = h / 3 * (sub_areas(end - 1) + 2*extra_area);
end

在这个例子中，`f` 是你要积分的函数，`a` 和 `b` 是积分区间的边界，`n` 是分割的子区间数量。注意，实际应用中可能需要对函数连续性和奇偶性做一些检查。