simpson4.asv

来自「梯形与simpson求积公式 分别是梯形,simpson,和3/8simpso」· ASV 代码 · 共 9 行

function I = simpson4(fun,a,b,npanel)

n = 4*npanel + 1;
h = (b-a)/(n-1);
x = a:h:b;
f = feval(fun,x);

I = (3*h/8) * ( f(1) + 4*sum(f(2:2:n-1)) + 2*sum(f(3:2:n-2)) + f(n) );
%              f(a)    f_even             f_odd              f(b)