+

来自「此文挡是用Matlab开发的几种计算数值微积分的方法」· 代码 · 共 31 行

part1
function s=simpr1(f,a,b,N)
h=(b-a)/(2*N);
s1=0;
s2=0;
for k=1:N
    x=a+h*(2*k-1);
    s1=s1+feval(f,x);
end
for k=1:(N-1)
    x=a+h*2*k;
    s2=s2+feval(f,x);
end
s=h*(feval(f,a)+feval(f,b)+4*s1+2*s2)/3;
************************************
part2
function y1=f(x)
y1=exp(x);
************************************
test
>> s=simpr1('f',0,1,1)

s =

    1.7189

>> vpa(s)
 
ans =
 
1.71886