trap.m

来自「基于matlab的数值计算试验程序。equ_direct」· M 代码 · 共 8 行

function I=Trap(f,a,b,n)
%复化梯形公式
h=(b-a)/n;S=feval(f,a);
for i=1:n-1
    x(i)=a+h*i;
    S=S+2*feval(f,x(i));
end
S=S+feval(f,b);I=h*S/2;