simprl.m

来自「用于求积分」· M 代码 · 共 15 行

function s=simprl(f,a,b,n)
h=(b-a)/(2*n);
s1=0;
s2=0;
for k=1:n
    x=a+h*(2*k-1);
    q=feval(f,x);
    s1=s1+feval(f,x);
end
for k=1:(n-1)
    x=a+h*2*k;
    w=feval(f,x);
    s2=s2+feval(f,x);
end
s=h*(feval(f,a)+feval(f,b)+4*s1+2*s2)/3;