数值分析复化梯形公式matlab,数值分析复化梯形公式,复化Simpson公式MATLAB程序

以下是MATLAB代码实现复化梯形公式和复化Simpson公式：

复化梯形公式：

function [T, h] = composite_trapezoidal(f, a, b, n)
% f: 被积函数
% a, b: 积分区间端点
% n: 区间数
% T: 积分近似值
% h: 区间长度

h = (b - a) / n;
x = a:h:b;
y = f(x);
T = h/2 * (y(1) + 2*sum(y(2:n)) + y(n+1));
end

复化Simpson公式：

function [S, h] = composite_simpson(f, a, b, n)
% f: 被积函数
% a, b: 积分区间端点
% n: 区间数
% S: 积分近似值
% h: 区间长度

h = (b - a) / n;
x = a:h:b;
y = f(x);
S = h/3 * (y(1) + 4*sum(y(2:2:n)) + 2*sum(y(3:2:n-1)) + y(n+1));
end

其中，`f`为被积函数，`a`和`b`为积分区间的端点，`n`为区间数。输出的`T`和`S`为积分的近似值，`h`为区间长度。