trapr1.m

余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程的数值解法；3插值法；4.常微分方程的解法

function s=trapr1(f1431,a,b,n)
%  f是被积函数。
%  a,b分别为积分的上下限。
%  n是子区间的个数。
%  s是梯形总面积。
h=(b-a)/n;
s=0;
for k=1:(n-1)
    x=a+h*k;
    s=s+feval('f1431',x);
end
s=h*(feval('f1431',a)+feval('f1431',b))/2+h*s;