给出复合辛普森公式求积分的matlab实现代码

复合辛普森公式是一种数值积分方法，可以用来求解函数的定积分。下面是一个基于 Matlab 的实现代码：

function result = composite_simpson(f, a, b, n)
% f: 被积函数
% a, b: 积分区间
% n: 划分数
% result: 积分结果

% 计算步长
h = (b - a) / n;

% 计算x和f(x)的值
x = a:h:b;
fx = f(x);

% 使用复合辛普森公式计算积分
sum1 = sum(fx(2:2:end-2));
sum2 = sum(fx(3:2:end-1));
result = h / 3 * (fx(1) + 4 * sum1 + 2 * sum2 + fx(end));
end

函数 `composite_simpson` 接收 4 个参数，分别是被积函数 `f`、积分区间 `[a, b]`、划分数 `n`，以及返回的积分结果 `result`。其中，输入的被积函数 `f` 必须是一个函数句柄或者函数字符串，可以是任何定义在 Matlab 中的函数；输入的积分区间 `[a, b]` 必须是一个 1x2 的向量，分别表示积分区间的起点和终点；输入的划分数 `n` 必须是一个正整数，表示将积分区间分成多少个小区间。

在函数体内，首先计算步长 `h`，然后根据步长计算积分区间中的离散点 `x` 和对应的函数值 `fx`。接下来使用复合辛普森公式计算积分结果，并将结果存储在 `result` 中，最后返回给调用者。