牛顿-柯特斯公式matlab

牛顿-柯特斯公式是一种数值积分方法，可以用于计算函数在一定区间上的定积分。在MATLAB中，可以使用以下代码实现牛顿-柯特斯公式：

function I = ncint(f, a, b, m, n)
% f为被积函数，a、b为积分区间，m为最高阶数，n为划分数
h = (b-a)/n; % 计算划分宽度
x = a:h:b; % 划分区间
y = f(x); % 计算函数值
% 计算系数c
c = zeros(m+1,1);
c(1) = 1/2;
for k = 1:m
    c(k+1) = 1/(4^k-1)*(4^k*c(k)-sum(c(1:k)));
end
% 计算积分值
I = h*sum(c.*y(1:n+1)');

其中，f为被积函数，a、b为积分区间，m为最高阶数，n为划分数。函数先计算出划分区间和函数值，然后根据公式计算系数c，最后根据公式计算积分值。