matlab使用复化simpson公式计算的代码案例

以下是使用复化Simpson公式计算函数f(x)在区间[a,b]上的定积分的MATLAB代码案例：

function I = simpson(f,a,b,n)
% f: 被积函数
% a: 积分下限
% b: 积分上限
% n: 区间数（必须为偶数）
% I: 积分结果

% 计算区间宽度
h = (b-a)/n;

% 计算节点
x = linspace(a,b,n+1);

% 计算f(x)的值
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

使用方法：将要求的被积函数f(x)定义为一个函数句柄，然后调用simpson函数即可。例如，计算函数f(x)=x^2在区间[0,1]上的定积分，可以这样做：

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

其中，"^"表示乘方运算。