simpr1.m

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

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