trapezoid.m

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

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));