func4.m

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

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