matlab复化simpson积分

以下是matlab实现复合辛普森积分的代码：

clc; clear;
a = 0; b = 1; % 积分范围
num = 1000; % 积分准确度
h = (b-a)/(2*num);
f = @(x) exp(-x); % 积分表达式
I = 0; % 积分结果
I = f(a) - f(b);
for i = 1:num
    I = I + (2*f(a+2*i*h) + 4*f(a+(2*i-1)*h));
end
I = I * h / 3;
disp(I);

其中，`a`和`b`分别为积分区间的上下限，`num`为积分准确度，`f`为被积函数，`I`为积分结果。在循环中，根据复合辛普森公式计算积分结果。最后将积分结果乘以`h/3`即可得到最终结果。