用matlab写一个复合牛顿-柯特斯公式的代码

以下是用matlab写的复合牛顿-柯特斯公式的代码：

function [y] = newton_cotes(f, a, b, n)
% f: 被积函数, a, b: 积分区间, n: 区间数
x = linspace(a, b, n+1); % 将积分区间分成n个小区间
h = (b - a) / n; % 小区间长度
y = f(x(1)) + f(x(end)); % 初始化积分结果
for i = 1:n-1 % 遍历每个小区间
    % 使用3次牛顿-柯特斯公式进行积分
    y = y + 2 * f(x(i+1)) + 3 * f((2*x(i+1) + x(i))/3) + 3 * f((x(i+1) + 2*x(i))/3);
end
y = y * 3 * h / 8; % 将积分结果乘以系数
end

使用方法：

定义被积函数f(x)，调用函数newton_cotes(f, a, b, n)即可求出在[a, b]上的积分值。其中，n为区间数，可以根据需要进行调整。