matlab复化辛普森法求定积分

以下是使用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); % 输出积分结果