function  t=natrapz(fname,a,b,n)
%用途：定步長梯形法求函數(shù)的積分
%格式：t=natrapz(fname,a,b,n) fname是被積函數(shù)，a,b分別為下上限，n為等份數(shù)
h=(b-a)/n;
fa=feval(fname,a);fb=feval(fname,b);f=feval(fname,a+h:h:b-h+eps);
t=h*(0.5*(fa+fb)+sum(f));