simpson.m

来自「简单实现了复化梯形公式的求解」· M 代码 · 共 21 行

%复化Simpson公式
function s=fsimpson(f,a,b,n)
%f表示被积函数句柄
%a,b表示被积区间[a,b]的端点
%n表示区间个数
%s是用复化simpson公式求得的积分值

h=(b-a)/n;
fa=feval(f,a);
fb=feval(f,b);
s=fb+fa;
x=a;
for i=1:n
x=x+h/2;
fx=feval(f,x);
s=s+4*fx;
fx=feval(f,x);
s=s+2*fx;
end
s=h*s/6;