复化辛普森公式matlab

在MATLAB中，可以使用以下代码实现复化辛普森公式：

假设要求函数f(x)在区间[a,b]上的积分，将区间[a,b]分成n个小区间，则每个小区间的长度为h=(b-a)/n。

代码如下：

function res = simpson(f, a, b, n)
% 复化辛普森公式求函数f在区间[a,b]上的积分
% 输入参数：
% f：被积函数句柄
% a：积分下限
% b：积分上限
% n：小区间数
% 输出参数：
% res：积分结果

h = (b-a)/n;
x = a:h:b;
y = f(x);
res = h/3*(y(1) + y(end) + 4*sum(y(2:2:end-1)) + 2*sum(y(3:2:end-2)));
end

其中，h表示每个小区间的长度，x表示小区间的端点，y表示小区间端点对应的函数值。

函数f需要通过句柄传入，例如：

f = @(x) exp(-x.^2);
a = 0;
b = 1;
n = 100;
res = simpson(f, a, b, n);

以上代码计算了函数exp(-x^2)在区间[0,1]上的积分，n取100。