simpson积分matlab程序

以下是一个使用Simpson积分法计算函数在给定区间上的积分的MATLAB程序：

function I = simpson(f, a, b, n)
% f: 被积函数
% a: 积分下限
% b: 积分上限
% n: 区间数（必须为偶数）
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));
end

该程序定义了一个名为simpson的函数，该函数接受四个参数：被积函数f、积分下限a、积分上限b和区间数n。它使用Simpson积分法计算函数f在区间[a, b]上的积分，并返回结果I。

使用该程序的示例：

>> f = @(x) exp(-x.^2); % 定义被积函数
>> a = 0; b = 1; % 定义积分区间
>> n = 100; % 定义区间数
>> I = simpson(f, a, b, n) % 计算积分
I =
    0.7468