natrapz.m

来自「Matlab数值解算法实现代码」· M 代码 · 共 8 行

function  t=natrapz(fname,a,b,n)
%用途：定步长梯形法求函数的积分
%格式：t=natrapz(fname,a,b,n) fname是被积函数，a,b分别为下上限，n为等份数
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));