trapr1.m

来自「MATLAB程序设计教程的源码」· M 代码 · 共 9 行

function s=trapr1(f,a,b,n)
h=(b-a)/n;
s=0;
for k=1:n-1
    x=a+h*k;
    s=s+feval(f,x);
end
s=h*(feval(f,a)+feval(f,b))/2+h*s;