simpson4.m

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

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

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

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