trapr1.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 13 行

function s=trapr1(f1431,a,b,n)
%  f是被积函数。
%  a,b分别为积分的上下限。
%  n是子区间的个数。
%  s是梯形总面积。
h=(b-a)/n;
s=0;
for k=1:(n-1)
    x=a+h*k;
    s=s+feval('f1431',x);
end
s=h*(feval('f1431',a)+feval('f1431',b))/2+h*s;