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-1
sum=sum+f(a+k.*h);
end
Tn=(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-1
sum1 = sum1 + f(a+(i+1/2).*h);
end
for j = 1:n-1
sum2 = sum2 + f(a+j.*h);
end
Sn = h/6*(f(a)+4*sum1+2*sum2+f(b));