traprl.m

来自「用于求积分」· M 代码 · 共 10 行

function s=traprl(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;