simpson.asv

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

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

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

I = (h/3) * ( f(1) + sum(f(2:n-1)) + 0.5*f(n) );