ftrapz.m

Matlab求解方程

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;