f.m文件：

function f = f(x)f = (x^2)*sin(x);end

复化梯形公式

Tn.m文件：

function Tn=Tn(n) % n代表区间数a = -2; % 区间下界b = 2; % 区间下界h=(b-a)/n;sum=0;for k=1:n-1sum=sum+f(a+k.*h);endTn=(f(a)+2*sum+f(b))*h/2;end

复化Simpson公式

Sn.m文件：

function Sn = Sn(n)a = -2;b = 2;h = (b-a)/n;sum1 = 0;sum2 = 0;for i = 0:n-1sum1 = sum1 + f(a+(i+1/2).*h);endfor j = 1:n-1sum2 = sum2 + f(a+j.*h);endSn = h/6*(f(a)+4*sum1+2*sum2+f(b));

如果觉得《MATLAB 复化梯形公式 复化Simpson公式》对你有帮助，请点赞、收藏，并留下你的观点哦！