srule.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 10 行

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;