ftrapz.m

来自「Matlab求解方程」· M 代码 · 共 9 行

function I=ftrapz(fh,a,b,n)
h=(b-a)/n;
x=a:h:b;
np1=n+1;
c=ones(np1,1);
c(1,1)=0.5;
c(np1,1)=0.5;
f=feval(fh,x);
I=h*f*c;