trap.m

基于matlab的数值计算试验程序。equ_direct

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;