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

好的，下面是一个用MATLAB编写的辛普森积分函数，可以用于求解任意被积函数的数值定积分：

function y = simpson(f, a, b, n)
% f: 被积函数
% a, b: 积分区间左右端点
% n: 分段数（必须为偶数）
h = (b - a) / n;
x = a:h:b;
y = f(x);
y(2:2:end-1) = 4 * y(2:2:end-1);
y(3:2:end-2) = 2 * y(3:2:end-2);
y = h / 3 * sum(y);
end

其中，输入参数 `f` 表示被积函数，`a` 和 `b` 分别表示积分区间的左右端点，`n` 表示分段数，必须为偶数。输出结果 `y` 表示积分的数值近似值。

例如，如果要求解函数 `f(x) = x^2` 在区间 [0, 1] 的积分，可以按照以下步骤进行：

1. 定义被积函数：`f = @(x) x.^2;`
2. 调用 `simpson()` 函数求解积分：`y = simpson(f, 0, 1, 100);`
3. 输出积分结果：`disp(y);`

完整的代码如下所示：

f = @(x) x.^2;
y = simpson(f, 0, 1, 100);
disp(y);

输出结果为：

0.3333

表示函数 `f(x) = x^2` 在区间 [0, 1] 内的积分结果为 1/3。