复化辛普森公式MATLAB代码

复化辛普森公式是一种数值积分方法，用于计算函数在给定区间上的定积分。下面是使用MATLAB实现复化辛普森公式的代码：

function result = composite_simpson(f, a, b, n)
    % f: 被积函数
    % a: 积分下限
    % b: 积分上限
    % n: 划分的子区间数
    
    h = (b - a) / n; % 子区间宽度
    x = a:h:b; % 子区间节点
    
    % 计算每个子区间的积分值
    y = f(x);
    
    % 计算复化辛普森公式的近似积分值
    result = (h/3) * (y(1) + 4*sum(y(2:2:end-1)) + 2*sum(y(3:2:end-2)) + y(end));
end

使用方法：
1. 定义被积函数f(x)，例如：`f = @(x) x.^2;`
2. 调用`composite_simpson`函数计算积分值，例如：`result = composite_simpson(f, a, b, n);`

其中，`f`为被积函数，`a`和`b`为积分区间的下限和上限，`n`为划分的子区间数。返回值`result`为复化辛普森公式的近似积分值。