srule.m

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

function Z=srule(f1421,a0,b0)
%  f是被积函数。
%  a0和b0是被积区间的上下限。
%  输出的值是所求的积分值。
h=(b0-a0)/2;
C=zeros(1,3);
C=[feval('f1421',a0),feval('f1421',(a0+b0)/2),feval('f1421',b0)];
S=h*(C(1)+4*C(2)+C(3))/3;
Z=S;