func5ok.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 8 行

n=500
a=linspace(0,1,n)
s=0;
for i=2:n-1       %i can't start from 1,else the denominator will be 0
s=s+func5(a(i))+4*func5(a(i)+1/n/2)+func5(a(i+1))
disp(s)
end
s=s*1/n/6