用MATLAB编写一个函数，功能为采用辛普森积分公式对已知被积函数求数值定积分

可以编写一个 MATLAB 函数，实现对已知被积函数的数值定积分，具体实现方式如下：

function I = Simpson(f, a, b, n)
% 输入参数：
% f - 被积函数的句柄
% a - 积分下限
% b - 积分上限
% n - 区间分割数（必须为偶数）
% 输出参数：
% I - 数值定积分结果

if mod(n, 2) ~= 0 % 区间分割数必须为偶数
    error('区间分割数必须为偶数');
end

h = (b - a) / n; % 子区间宽度
x = a:h:b; % 子区间节点
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

使用时，可以先定义被积函数的句柄，再调用 Simpson 函数进行数值计算，例如：

f = @(x) sin(x); % 定义被积函数
a = 0; % 积分下限
b = pi / 2; % 积分上限
n = 100; % 区间分割数
I = Simpson(f, a, b, n); % 计算数值定积分
disp(I); % 显示积分结果

需要注意的是，由于 Simpson 公式要求区间分割数为偶数，因此输入参数 n 必须为偶数，否则会报错。