trapezoid.m

用于求积分

function I=trapezoid(fun,a,b,npanel)
n=npanel+1;
h=(b-a)/(n-1);
x=a:h:b;
f=feval(fun,x);
I=h*(0.5*f(1)+sum(f(2:n-1))+0.5*f(n));