fsimpson.asv

来自「数值分析各种公式的解法」· ASV 代码 · 共 12 行

function f=fx
f=exp(sin(x));
function fsimpson(a,b,n)
h=(b-a)/n;
t=h*(f(a)+f(b))/6;
for i=0:(n-1)
    t=t+h*(4*f(a+(i+1/2)*h)+2*f(a+i*h))/6;
end
I=t
function f=f(x)
f=exp(sin(x));